FairRoot/PandaRoot
Functions
sim_complete_pbarA.C File Reference

Go to the source code of this file.

Functions

int sim_complete_pbarA (Int_t nEvents=100, UInt_t seed=-1, TString SimEngine="TGeant3", Float_t mom=8.0, Float_t atarg=20.0)
 

Function Documentation

int sim_complete_pbarA ( Int_t  nEvents = 100,
UInt_t  seed = -1,
TString  SimEngine = "TGeant3",
Float_t  mom = 8.0,
Float_t  atarg = 20.0 
)

Definition at line 8 of file sim_complete_pbarA.C.

References allDigiFile, BeamMomentum, Bool_t, boxGen, Cave, ctime, digiFile, Dipole, Double_t, Drc, Dsk, Emc, emcHitProd, fField, fIni, fRun, FTof, Fts, gDebug, Gem, kParameterMerged, mom, Muo, Mvd, nEvents, output, parIo1, Pipe, primGen, printf(), rtdb, rtime, seed, PndMdt::SetBarrel(), PndMdt::SetEndcap(), PndMdt::SetForward(), PndEmc::SetGeometryVersion(), PndMdt::SetMdtCoil(), PndMdt::SetMdtMagnet(), PndMdt::SetMdtMFIron(), PndMdt::SetMuonFilter(), PndDrc::SetRunCherenkov(), PndEmc::SetStorageOfData(), PndDsk::SetStoreCerenkovs(), PndDsk::SetStoreTrackPoints(), PndEvtGenDirect::SetStoreTree(), sqrt(), Stt, timer, and TString.

9 {
10  //-----User Settings:-----------------------------------------------
11  TString OutputFile ="sim_complete.root";
12  TString ParOutputfile ="simparams.root";
13  TString MediaFile ="media_pnd.geo";
14  gDebug = 0;
15  TString digiFile = "all.par"; //The emc run the hit producer directly
16  // choose your event generator
17  Bool_t UseEvtGenDirect =kTRUE;
18  Bool_t UseDpm =kFALSE;
19  Bool_t UseFtf =kFALSE;
20  Bool_t UseBoxGenerator =kFALSE;
21 
22  Double_t BeamMomentum = 0.; // beam momentum ONLY for the scaling of the dipole field.
23  if (UseBoxGenerator)
24  {
25  BeamMomentum =15.0; // ** change HERE if you run Box generator
26  }
27  else
28  {
29  BeamMomentum = mom; // for DPM/EvtGen BeamMomentum is always = mom
30  }
31  //------------------------------------------------------------------
32  //TLorentzVector fIni(0, 0, mom, sqrt(mom*mom+9.3827203e-01*9.3827203e-01)+9.3827203e-01);
33  TLorentzVector fIni(0, 0, mom, sqrt(mom*mom+9.3827203e-01*9.3827203e-01)+atarg*0.931494);
34  TDatabasePDG::Instance()->AddParticle("pbarASystem","pbarASystem",fIni.M(),kFALSE,0.1,0, "",88890);
35  TDatabasePDG::Instance()->AddParticle("pbarpSystem","pbarpSystem",5.5,kFALSE,0.1,0, "",88888);
36  TDatabasePDG::Instance()->AddParticle("pbarnSystem","pbarnSystem",5.5,kFALSE,0.1,-3, "",88887);
37  TDatabasePDG::Instance()->AddParticle("A-1System","pbarASystem",1.8,kFALSE,0.1,0, "",88891);
38  //------------------------------------------------------------------
39  TStopwatch timer;
40  timer.Start();
41  //gRandom->SetSeed();
42  gRandom->SetSeed(seed);
43 
44  // Create the Simulation run manager--------------------------------
45  FairRunSim *fRun = new FairRunSim();
46  fRun->SetName(SimEngine.Data() );
47  fRun->SetOutputFile(OutputFile.Data());
48  fRun->SetGenerateRunInfo(kFALSE);
49  fRun->SetBeamMom(BeamMomentum);
50  fRun->SetMaterials(MediaFile.Data());
51  fRun->SetUseFairLinks(kTRUE);
52  FairRuntimeDb *rtdb=fRun->GetRuntimeDb();
53 
54  // Set the parameters
55  //-------------------------------
56  TString allDigiFile = gSystem->Getenv("VMCWORKDIR");
57  allDigiFile += "/macro/params/";
58  allDigiFile += digiFile;
59 
60 
61  //-------Set the parameter output --------------------
62  FairParAsciiFileIo* parIo1 = new FairParAsciiFileIo();
63  parIo1->open(allDigiFile.Data(),"in");
64  rtdb->setFirstInput(parIo1);
65 
66  //---------------------Set Parameter output ----------
67  Bool_t kParameterMerged=kTRUE;
68  FairParRootFileIo* output=new FairParRootFileIo(kParameterMerged);
69  output->open(ParOutputfile.Data());
70  rtdb->setOutput(output);
71 
72  // Create and add detectors
73 
74  //------------------------- CAVE -----------------
75 
76  FairModule *Cave= new PndCave("CAVE");
77  Cave->SetGeometryFileName("pndcave.geo");
78  fRun->AddModule(Cave);
79  //------------------------- Magnet -----------------
80  //FairModule *Magnet= new PndMagnet("MAGNET");
81  //Magnet->SetGeometryFileName("FullSolenoid_V842.root");
82  //Magnet->SetGeometryFileName("FullSuperconductingSolenoid_v831.root");
83  //fRun->AddModule(Magnet);
84  FairModule *Dipole= new PndMagnet("MAGNET");
85  Dipole->SetGeometryFileName("dipole.geo");
86  fRun->AddModule(Dipole);
87  //------------------------- Pipe -----------------
88  FairModule *Pipe= new PndPipe("PIPE");
89  Pipe->SetGeometryFileName("beampipe_201309.root");
90  fRun->AddModule(Pipe);
91  //------------------------- STT -----------------
92  FairDetector *Stt= new PndStt("STT", kTRUE);
93  Stt->SetGeometryFileName("straws_skewed_blocks_35cm_pipe.geo");
94  fRun->AddModule(Stt);
95  //------------------------- MVD -----------------
96  FairDetector *Mvd = new PndMvdDetector("MVD", kTRUE);
97  Mvd->SetGeometryFileName("Mvd-2.1_FullVersion.root");
98  fRun->AddModule(Mvd);
99  //------------------------- GEM -----------------
100  FairDetector *Gem = new PndGemDetector("GEM", kTRUE);
101  Gem->SetGeometryFileName("gem_3Stations_Tube.root");
102  fRun->AddModule(Gem);
103  //------------------------- EMC -----------------
104  PndEmc *Emc = new PndEmc("EMC",kTRUE);
105  Emc->SetGeometryVersion(1);
106  Emc->SetStorageOfData(kFALSE);
107  fRun->AddModule(Emc);
108  //------------------------- SCITIL -----------------
109  FairDetector *SciT = new PndSciT("SCIT",kTRUE);
110  SciT->SetGeometryFileName("SciTil_201601.root");
111  fRun->AddModule(SciT);
112  //------------------------- DRC -----------------
113  PndDrc *Drc = new PndDrc("DIRC", kTRUE);
114  Drc->SetGeometryFileName("dirc_l0_p0_updated.root");
115  Drc->SetRunCherenkov(kFALSE);
116  fRun->AddModule(Drc);
117  //------------------------- DISC -----------------
118  PndDsk* Dsk = new PndDsk("DSK", kTRUE);
119  Dsk->SetStoreCerenkovs(kFALSE);
120  Dsk->SetStoreTrackPoints(kFALSE);
121  fRun->AddModule(Dsk);
122  //------------------------- MDT -----------------
123  PndMdt *Muo = new PndMdt("MDT",kTRUE);
124  Muo->SetBarrel("fast");
125  Muo->SetEndcap("fast");
126  Muo->SetMuonFilter("fast");
127  Muo->SetForward("fast");
128  Muo->SetMdtMagnet(kTRUE);
129  Muo->SetMdtCoil(kTRUE);
130  Muo->SetMdtMFIron(kTRUE);
131  fRun->AddModule(Muo);
132  //------------------------- FTS -----------------
133  FairDetector *Fts= new PndFts("FTS", kTRUE);
134  Fts->SetGeometryFileName("fts.geo");
135  fRun->AddModule(Fts);
136  //------------------------- FTOF -----------------
137  FairDetector *FTof = new PndFtof("FTOF",kTRUE);
138  FTof->SetGeometryFileName("ftofwall.root");
139  fRun->AddModule(FTof);
140  //------------------------- RICH ----------------
141  FairDetector *Rich= new PndRich("RICH",kFALSE);
142  Rich->SetGeometryFileName("rich_v2_shift.geo");
143  fRun->AddModule(Rich);
144 
145  // Create and Set Event Generator
146  //-------------------------------
147  FairPrimaryGenerator* primGen = new FairPrimaryGenerator();
148  fRun->SetGenerator(primGen);
149 
150  if(UseBoxGenerator){ // Box Generator
151  FairBoxGenerator* boxGen = new FairBoxGenerator(22, 5); // 13 = muon; 1 = multipl.
152  boxGen->SetPRange(mom,mom); // GeV/c
153  boxGen->SetPhiRange(0., 360.); // Azimuth angle range [degree]
154  boxGen->SetThetaRange(0., 90.); // Polar angle in lab system range [degree]
155  boxGen->SetXYZ(0., 0., 0.); // cm
156  primGen->AddGenerator(boxGen);
157  }
158  if(UseDpm){
159  PndDpmDirect *Dpm= new PndDpmDirect(mom,1);
160  primGen->AddGenerator(Dpm);
161  }
162  if(UseFtf){
163  // TString macfile = gSystem->Getenv("VMCWORKDIR");
164  // macfile += "/pgenerators/FtfEvtGen/PbarP.mac";
165  // PndFtfDirect *Ftf = new PndFtfDirect(macfile.Data());
166  PndFtfDirect *Ftf = new PndFtfDirect("anti_proton", "G4_H", 1, "ftfp", mom, 123456);
167  primGen->AddGenerator(Ftf);
168  }
169  if(UseEvtGenDirect){
170  TString EvtInput =gSystem->Getenv("VMCWORKDIR");
171  //EvtInput+="/macro/run/psi2s_Jpsi2pi_Jpsi_mumu.dec";
172  EvtInput+="/mysim/pbarNe20_0/pbarNe20_pbarpX.dec";
173  TString defDECAY =gSystem->Getenv("VMCWORKDIR");
174  defDECAY+="/pgenerators/EvtGen/EvtGen/Private/DECAY.DEC";
175  TString defPDL =gSystem->Getenv("VMCWORKDIR");
176  defPDL+="/pgenerators/EvtGen/EvtGen/Private/evt.pdl";
177  //PndEvtGenDirect *EvtGen = new PndEvtGenDirect("pbarpSystem", EvtInput.Data(), mom);
178  PndEvtGenDirect *EvtGen = new PndEvtGenDirect("pbarASystem",EvtInput.Data(),mom,seed,defDECAY.Data(),defPDL.Data(),atarg);
179  EvtGen->SetStoreTree(kTRUE);
180  primGen->AddGenerator(EvtGen);
181  }
182 
183  //---------------------Create and Set the Field(s)----------
184  PndMultiField *fField= new PndMultiField("AUTO");
185  fRun->SetField(fField);
186 
187  // EMC Hit producer
188  //-------------------------------
189  PndEmcHitProducer* emcHitProd = new PndEmcHitProducer();
190  fRun->AddTask(emcHitProd);
191 
192  //------------------------- Initialize the RUN -----------------
193  fRun->Init();
194  //------------------------- Run the Simulation -----------------
195  fRun->Run(nEvents);
196  //------------------------- Save the parameters -----------------
197  rtdb->saveOutput();
198  //------------------------Print some info and exit----------------
199  timer.Stop();
200  Double_t rtime = timer.RealTime();
201  Double_t ctime = timer.CpuTime();
202  printf("RealTime=%f seconds, CpuTime=%f seconds\n",rtime,ctime);
203 
204  cout << " Test passed" << endl;
205  cout << " All ok " << endl;
206 
207  //exit(0);
208  return 0;
209 
210 };
Drc
PndDrc * Drc
Definition: sim_emc_apd.C:75
fField
PndMultiField * fField
Definition: sim_emc_apd.C:97
PndMdt::SetForward
void SetForward(TString name)
Definition: PndMdt.h:34
printf
printf("RealTime=%f seconds, CpuTime=%f seconds\n", rtime, ctime)
FTof
FairDetector * FTof
Definition: sim_ftof.C:49
kParameterMerged
Bool_t kParameterMerged
Definition: sim_emc_apd.C:113
sqrt
friend F32vec4 sqrt(const F32vec4 &a)
Definition: P4_F32vec4.h:29
PndEmcHitProducer
creates PndEmcHits from PndEmcPoints
Definition: PndEmcHitProducer.h:37
Emc
PndEmc * Emc
Definition: sim_emc_apd.C:55
PndMdt::SetMdtCoil
void SetMdtCoil(bool opt=false)
Definition: PndMdt.h:28
digiFile
TString digiFile
Definition: bump_emc.C:20
Mvd
FairDetector * Mvd
Definition: sim_emc_apd.C:51
PndMdt::SetMdtMFIron
void SetMdtMFIron(bool opt=false)
Definition: PndMdt.h:29
mom
Double_t mom
Definition: plot_dirc.C:14
PndEmc::SetStorageOfData
void SetStorageOfData(Bool_t val)
Definition: PndEmc.cxx:941
BeamMomentum
double BeamMomentum
Definition: sim_ftof_stof.C:17
allDigiFile
TString allDigiFile
Definition: hit_muo.C:36
output
FairParRootFileIo * output
Definition: sim_emc_apd.C:120
Muo
PndMdt * Muo
Definition: sim_emc_apd.C:67
PndEmc
Simulation of EMC.
Definition: PndEmc.h:26
primGen
FairPrimaryGenerator * primGen
Definition: sim_emc_apd.C:81
Dsk
FairDetector * Dsk
Definition: run_DpmSim.C:66
fRun
FairRunAna * fRun
Definition: hit_dirc.C:58
Gem
FairDetector * Gem
Definition: runJohan.C:71
PndMdt::SetMdtMagnet
void SetMdtMagnet(bool opt=false)
Definition: PndMdt.h:27
PndDsk::SetStoreTrackPoints
void SetStoreTrackPoints(Bool_t storeTrackPoints)
Definition: PndDsk.h:148
PndEvtGenDirect::SetStoreTree
void SetStoreTree(Bool_t store=true)
Definition: PndEvtGenDirect.h:47
Stt
FairDetector * Stt
Definition: sim_emc_apd.C:47
PndDsk::SetStoreCerenkovs
void SetStoreCerenkovs(Bool_t storeCerenkovs)
Definition: PndDsk.h:146
Double_t
Double_t
Definition: SimCompleteLinkDef.h:6
Dipole
FairModule * Dipole
Definition: sim_emc_apd.C:40
Cave
FairModule * Cave
Definition: sim_emc_apd.C:32
PndDrc
Definition: PndDrc.h:31
nEvents
Int_t nEvents
Definition: hit_dirc.C:11
seed
unsigned int seed
Definition: createPndLmdMisalignmentMatrices.C:5
timer
TStopwatch timer
Definition: hit_dirc.C:51
PndEmc::SetGeometryVersion
virtual void SetGeometryVersion(const Int_t GeoNumber)
Definition: PndEmc.cxx:966
PndMdt::SetBarrel
void SetBarrel(TString name)
Definition: PndMdt.h:31
rtdb
FairRuntimeDb * rtdb
Definition: hit_dirc.C:66
gDebug
gDebug
Definition: sim_emc_apd.C:6
emcHitProd
PndEmcHitProducer * emcHitProd
Definition: full_tmpWaveform_emc.C:67
ctime
Double_t ctime
Definition: hit_dirc.C:114
parIo1
FairParAsciiFileIo * parIo1
Definition: bump_emc.C:53
boxGen
FairBoxGenerator * boxGen
Definition: sim_emc_apd.C:85
PndMdt::SetEndcap
void SetEndcap(TString name)
Definition: PndMdt.h:32
PndMdt::SetMuonFilter
void SetMuonFilter(TString name)
Definition: PndMdt.h:33
PndDrc::SetRunCherenkov
void SetRunCherenkov(Bool_t ch)
Definition: PndDrc.h:222
PndStt
Definition: PndStt.h:34
PndMdt
Definition: PndMdt.h:20
fIni
TLorentzVector fIni
Definition: full_core_ntp.C:29
Pipe
FairModule * Pipe
Definition: sim_emc_apd.C:44
rtime
Double_t rtime
Definition: hit_dirc.C:113
Fts
FairDetector * Fts
Definition: sim_ftof_stof.C:58
PndDsk
Definition: PndDsk.h:23
PndFts
Definition: PndFts.h:25
PndCave
Definition: PndCave.h:8
Generated on Wed Apr 3 2019 08:44:24 for FairRoot/PandaRoot by   doxygen 1.8.5