{ // ======================================================================== // Verbosity level (0=quiet, 1=event level, 2=track level, 3=debug) Int_t iVerbose = 3; Int_t nEvents = 1000; // ---- Load libraries ------------------------------------------------- gROOT->Macro("../Libs.C"); gSystem->Load("libGeane"); // ------------------------------------------------------------------------ // Output file TString parFile = "../data/MvdTrackingParams.root"; PndMvdFileNameCreator namecreator("../data/MvdTracking.root"); TString MCFile = namecreator.GetSimFileName(false).c_str(); TString DigiFile = namecreator.GetDigiFileName(false).c_str(); TString RecoFile = namecreator.GetRecoFileName(false).c_str(); TString TrackFile = namecreator.GetTrackFindingFileName(false).c_str(); TString outFile = namecreator.GetKalmanFileName(false).c_str(); std::cout << "DigiFile: " << DigiFile.Data()<< std::endl; std::cout << "RecoFile: " << RecoFile.Data()<< std::endl; std::cout << "TrackFinderFile: " << TrackFile.Data()<< std::endl; std::cout << "KalmanFile: " << outFile.Data() << std::endl; // --- Now choose concrete engines for the different tasks ------------- // ------------------------------------------------------------------------ // In general, the following parts need not be touched // ======================================================================== // ----- Timer -------------------------------------------------------- TStopwatch timer; timer.Start(); // ------------------------------------------------------------------------ // ----- Reconstruction run ------------------------------------------- CbmRunAna *fRun= new CbmRunAna(); fRun->SetInputFile(MCFile); // fRun->AddFriend(DigiFile); fRun->AddFriend(RecoFile); fRun->AddFriend(TrackFile); fRun->SetOutputFile(outFile); // ------------------------------------------------------------------------ // THIS IS STRONGLY NEEDED CbmGeane *Geane = new CbmGeane(MCFile); // PndMultiField *fField= new PndMultiField(); // PndTransMap *map= new PndTransMap("TransMap", "R"); // PndDipoleMap *map1= new PndDipoleMap("DipoleMap", "R"); // PndSolenoidMap *map2= new PndSolenoidMap("SolenoidMap", "R"); // fField->AddField(map); // fField->AddField(map1); // fField->AddField(map2); // Geane->SetField(fField); PndConstField *fMagField=new PndConstField(); fMagField->SetField(0, 0 ,20. ); // values are in kG // MinX=-75, MinY=-40,MinZ=-12 ,MaxX=75, MaxY=40 ,MaxZ=124 ); // values are in cm fMagField->SetFieldRegion(-500, 500,-500, 500, -200, 200); Geane->SetField(fMagField); // ----- Parameter database -------------------------------------------- CbmRuntimeDb* rtdb = fRun->GetRuntimeDb(); // CbmParRootFileIo* parInput1 = new CbmParRootFileIo(); CbmParAsciiFileIo* parInput1 = new CbmParAsciiFileIo(); parInput1->open(parFile.Data(),"in"); rtdb->setFirstInput(parInput1); /*Bool_t kParameterMerged=kTRUE; CbmParRootFileIo* output=new CbmParRootFileIo(kParameterMerged); output->open(parOutFile); rtdb->setOutput(output); */ fRun->LoadGeometry(); // ------------------------------------------------------------------------ // ========================================================================= // ====== Hit Producers ====== // ========================================================================= // ----- MVD hit producer -------------------------------------------- PndMvdKalmanTask* mvdKalman = new PndMvdKalmanTask(); mvdKalman->SetVerbose(iVerbose); fRun->AddTask(mvdKalman); // CbmParRootFileIo* output=new CbmParRootFileIo(kTRUE); // output->open(parOutFile.Data()); // rtdb->setOutput(output); rtdb->print(); // ===== End of HitProducers ===== // ========================================================================= PndMvdGeoPar* geoPar = (PndMvdGeoPar*)(rtdb->getContainer("PndMvdGeoPar")); // ----- Intialise and run -------------------------------------------- fRun->Init(); fRun->Run(0,nEvents); // ------------------------------------------------------------------------ mvdKalman->WriteHistograms("MvdKalmanHistos.root"); TFile histos("MvdKalmanHistos.root","READ"); // rtdb->saveOutput(); rtdb->print(); // ----- Finish ------------------------------------------------------- timer.Stop(); Double_t rtime = timer.RealTime(); Double_t ctime = timer.CpuTime(); cout << endl << endl; cout << "Macro finished succesfully." << endl; cout << "Output file is " << outFile << endl; cout << "Parameter file is " << parFile << endl; cout << "Real time " << rtime << " s, CPU time " << ctime << " s" << endl; cout << endl; // ------------------------------------------------------------------------ }