// Macro for running Cbm with Geant3 or Geant4 (M. Al-Turany , D. Bertini) // Modified 22/06/2005 D.Bertini { TStopwatch timer; timer.Start(); gDebug=0; int verboseLevel = 0; Int_t nEvents = 1000; //FileNames TString simOutput="MvdTpc_D+D-_10G_1000.root"; // TString simOutput="Mvd_DPMfixed_4GeV_10000.root"; TString parOutput="MvdParams.root"; // Load basic libraries // gROOT->Macro("$VMCWORKDIR/gconfig/rootlogon.C"); gROOT->Macro("$VMCWORKDIR/gconfig/rootlogon.C"); FairRunSim *fRun = new FairRunSim(); // set the MC version used // ------------------------ fRun->SetName("TGeant4"); // Choose the Geant Navigation System fRun->SetOutputFile(simOutput); // Set Material file Name //----------------------- fRun->SetMaterials("media_pnd.geo"); // Create and add detectors //------------------------- FairModule *Cave= new PndCave("CAVE"); Cave->SetGeometryFileName("pndcave.geo"); fRun->AddModule(Cave); FairModule *Magnet= new PndMagnet("MAGNET"); Magnet->SetGeometryFileName("FullSolenoid.root"); fRun->AddModule(Magnet); FairModule *Pipe= new PndPipe("PIPE"); //Pipe->SetGeometryFileName("pipebeamtarget.geo"); fRun->AddModule(Pipe); FairDetector *Mvd = new PndMvdDetector("MVD", kTRUE); Mvd->SetGeometryFileName("MVD_v1.0_woPassiveTraps.root"); // Mvd->SetGeometryFileName("MVD14.root"); Mvd->SetVerboseLevel(verboseLevel); fRun->AddModule(Mvd); FairDetector *Tpc = new PndTpcDetector("TPC", kTRUE); Tpc->SetGeometryFileName("tpc.geo"); fRun->AddModule(Tpc); // FairDetector *Stt= new PndStt("STT", kTRUE); // Stt->SetGeometryFileName("straws_skewed_blocks_35cm_pipe.geo"); // fRun->AddModule(Stt); PndEmc *Emc = new PndEmc("EMC",kTRUE); Emc->SetGeometryFileNameDouble("emc_module1245.dat","emc_module3new.root"); fRun->AddModule(Emc); /* FairDetector *Tof = new PndTof("TOF",kTRUE); Tof->SetGeometryFileName("tofbarrel.geo"); fRun->AddModule(Tof); PndDrc *Drc = new PndDrc("DIRC", kTRUE); Drc->SetRunCherenkov(kFALSE); // for fast sim Cherenkov -> kFALSE fRun->AddModule(Drc); */ // FairDetector *Emc = new CbmEmc("EMC",kTRUE); // Emc->SetGeometryFileName("emc_module12345.dat"); // fRun->AddModule(Emc); // FairDetector *Drc = new CbmDrc("DIRC", kTRUE); // Drc->SetGeometryFileName("dirc.geo"); // fRun->AddModule(Drc); FairPrimaryGenerator* primGen = new FairPrimaryGenerator(); fRun->SetGenerator(primGen); // Particle Generator (pdgid,mult, px,py,pz, vx,vy,vz) // single pions for testing // FairParticleGenerator* partGenX = new FairParticleGenerator(211,1, 1.,0.,0., 0.,0.,0.); // FairParticleGenerator* partGenY = new FairParticleGenerator(211,1, 0.,1.,0., 0.,0.,0.); // FairParticleGenerator* partGenZ = new FairParticleGenerator(211,1, 0.,0.1,1., 0.,0.,0.); // FairParticleGenerator* partGenXYZ = new FairParticleGenerator(211,1, 1.,1.,5., 0.,0.,0.); // primGen->AddGenerator(partGenX); // primGen->AddGenerator(partGenY); // primGen->AddGenerator(partGenZ); // primGen->AddGenerator(partGenXYZ); // Ion Generator // FairIonGenerator *fIongen= new FairIonGenerator(79, 197,79,1, 0.,0., 25, 0.,0.,-1.); // primGen->AddGenerator(fIongen); // Box Generator //Pions /* FairBoxGenerator *fBox = new FairBoxGenerator(211, 3); fBox->SetPRange(0.1,2.5); fBox->SetThetaRange(5,150); fBox->SetPhiRange(0.,360); fBox->SetCosTheta(); primGen->AddGenerator(fBox); */ //EvtGen Generator FairEvtGenGenerator* evtGen = new FairEvtGenGenerator("D+D-_10G_1000.evt"); primGen->AddGenerator(evtGen); // Urqmd Generator // FairUrqmdGenerator* urqmdGen = new FairUrqmdGenerator("../../input/00-03fm.100ev.f14"); // primGen->AddGenerator(urqmdGen); // DPM Generator //PndDpmGenerator* dpmGen = new PndDpmGenerator("/home/kph/stockman/fairroot/cbmroot/pandaroot/pgenerators/DpmEvtGen/Background-micro.root"); //primGen->AddGenerator(dpmGen); // Field Map Definition // -------------------- // 1- Reading the new field map in the old format fRun->SetBeamMom(15); PndMultiField *fField= new PndMultiField("FULL"); fRun->SetField(fField); fRun->SetStoreTraj(kTRUE); // toggle this for use with EVE fRun->Init(); //fRun->SetStoreTraj(kFALSE); // -Trajectories Visualization (TGeoManager Only ) // ----------------------------------------------- // Set cuts for storing the trajectpries // FairTrajFilter* trajFilter = FairTrajFilter::Instance(); // trajFilter->SetStepSizeCut(0.01); // 1 cm // trajFilter->SetVertexCut(-200., -200., -200, 200., 200., 200.); // trajFilter->SetMomentumCutP(10e-3); // p_lab > 10 MeV // trajFilter->SetEnergyCut(0., 1.02); // 0 < Etot < 1.04 GeV // trajFilter->SetStorePrimaries(kTRUE); // trajFilter->SetStoreSecondaries(kTRUE); // Fill the Parameter containers for this run //------------------------------------------- FairRuntimeDb *rtdb=fRun->GetRuntimeDb(); Bool_t kParameterMerged=kTRUE; FairParRootFileIo* output=new FairParRootFileIo(kParameterMerged); output->open(parOutput.Data(),"RECREATE"); rtdb->setOutput(output); // Transport nEvents // ----------------- fRun->Run(nEvents); // rtdb->saveOutput(); rtdb->print(); timer.Stop(); Double_t rtime = timer.RealTime(); Double_t ctime = timer.CpuTime(); printf("RealTime=%f seconds, CpuTime=%f seconds\n",rtime,ctime); }