1 #if !defined(__CINT__) || defined(__MAKECINT__)
6 #include <TVirtualMC.h>
8 #include "STEER/AliRunLoader.h"
9 #include "STEER/AliRun.h"
10 #include "STEER/AliConfig.h"
11 #include "PYTHIA6/AliDecayerPythia.h"
12 #include "EVGEN/AliGenCocktail.h"
13 #include "EVGEN/AliGenFixed.h"
14 #include "EVGEN/AliGenBox.h"
15 #include "EVGEN/AliGenScan.h"
16 #include "EVGEN/AliGenHIJINGpara.h"
17 #include "THijing/AliGenHijing.h"
18 #include "PYTHIA6/AliGenPythia.h"
19 #include "THerwig/AliGenHerwig.h"
20 #include "TIsajet/AliGenIsajet.h"
21 #include "TDPMjet/AliGenDPMjet.h"
22 #include "EVGEN/AliGenParam.h"
23 #include "EVGEN/AliGenMUONlib.h"
24 #include "EVGEN/AliGenPHOSlib.h"
25 #include "EVGEN/AliGenGSIlib.h"
26 #include "EVGEN/AliGenFLUKAsource.h"
27 #include "EVGEN/AliGenExtFile.h"
28 #include "EVGEN/AliGenHalo.h"
29 #include "EVGEN/AliGenReaderTreeK.h"
30 #include "STRUCT/AliBODY.h"
31 #include "STRUCT/AliMAG.h"
35 enum gentype_t {hijing, hijingParam, gun, box, pythia, herwig,
36 param1, param2, param3, param4,
37 cocktail, fluka, halo, ntuple, scan, isajet, dpmjet};
39 gentype_t gentype = dpmjet;
46 // Set Random Number seed
47 gRandom->SetSeed(12345); //Set 0 to use the current time
48 cout<<"Seed for random number generation= "<<gRandom->GetSeed()<<endl;
51 // libraries required by geant321
53 gSystem->Load("libgeant321");
55 gSystem->Load("libherwig.so");
56 gSystem->Load("libTHerwig.so");
57 gSystem->Load("libisajet.so");
58 gSystem->Load("libTIsajet.so");
59 gSystem->Load("libdpmjet.so");
60 gSystem->Load("libTDPMjet.so");
62 new TGeant3TGeo("C++ Interface to Geant3");
64 //=======================================================================
65 // Create the output file
69 cout<<"Config.C: Creating Run Loader ..."<<endl;
70 rl = AliRunLoader::Open("galice.root",
71 AliConfig::GetDefaultEventFolderName(),
75 gAlice->Fatal("Config.C","Can not instatiate the Run Loader");
78 rl->SetCompressionLevel(2);
79 rl->SetNumberOfEventsPerFile(100);
80 gAlice->SetRunLoader(rl);
83 // Set External decayer
84 AliDecayer* decayer = new AliDecayerPythia();
85 decayer->SetForceDecay(kAll);
87 gMC->SetExternalDecayer(decayer);
91 //=======================================================================
92 // ************* STEERING parameters FOR ALICE SIMULATION **************
93 // --- Specify event type to be tracked through the ALICE setup
94 // --- All positions are in cm, angles in degrees, and P and E in GeV
97 gMC->SetProcess("DCAY",1);
98 gMC->SetProcess("PAIR",1);
99 gMC->SetProcess("COMP",1);
100 gMC->SetProcess("PHOT",1);
101 gMC->SetProcess("PFIS",0);
102 gMC->SetProcess("DRAY",0);
103 gMC->SetProcess("ANNI",1);
104 gMC->SetProcess("BREM",1);
105 gMC->SetProcess("MUNU",1);
106 gMC->SetProcess("CKOV",1);
107 gMC->SetProcess("HADR",1);
108 gMC->SetProcess("LOSS",2);
109 gMC->SetProcess("MULS",1);
110 gMC->SetProcess("RAYL",1);
112 Float_t cut = 1.e-3; // 1MeV cut by default
113 Float_t tofmax = 1.e10;
115 gMC->SetCut("CUTGAM", cut);
116 gMC->SetCut("CUTELE", cut);
117 gMC->SetCut("CUTNEU", cut);
118 gMC->SetCut("CUTHAD", cut);
119 gMC->SetCut("CUTMUO", cut);
120 gMC->SetCut("BCUTE", cut);
121 gMC->SetCut("BCUTM", cut);
122 gMC->SetCut("DCUTE", cut);
123 gMC->SetCut("DCUTM", cut);
124 gMC->SetCut("PPCUTM", cut);
125 gMC->SetCut("TOFMAX", tofmax);
128 AliGenerator * gGener = 0x0;
132 //*********************************************
133 // Example for Fixed Particle Gun
134 //*********************************************
136 AliGenFixed *gener = new AliGenFixed(ntracks);
137 gener->SetMomentum(50);
140 gener->SetOrigin(0,0,0); //vertex position
141 gener->SetPart(kNeutron);
146 //*********************************************
147 // Example for Moving Particle Gun *
148 //*********************************************
150 AliGenBox *gener = new AliGenBox(ntracks);
151 gener->SetMomentumRange(3,4);
152 gener->SetPhiRange(0,360);
153 gener->SetThetaRange(90, 180. );
154 gener->SetOrigin(0,0,0);
156 gener->SetSigma(0,0,0); //Sigma in (X,Y,Z) (cm) on IP position
157 gener->SetPart(kMuonPlus);
162 //*********************************************
163 // Scanning on a grid *
164 //*********************************************
166 AliGenScan *gener = new AliGenScan(-1);
167 gener->SetMomentumRange(4,4);
168 gener->SetPhiRange(0,360);
169 gener->SetThetaRange(9,9);
171 gener->SetSigma(6,6,0); //Sigma in (X,Y,Z) (cm) on IP position
172 gener->SetPart(kMuonPlus);
173 gener->SetRange(20, -100, 100, 20, -100, 100, 1, 500, 500);
180 AliGenHIJINGpara *gener = new AliGenHIJINGpara(ntracks);
181 gener->SetMomentumRange(0,999);
182 gener->SetPhiRange(0,360);
183 gener->SetThetaRange(2,10);
184 gener->SetOrigin(0,0,0); //vertex position
185 gener->SetSigma(0,0,0); //Sigma in (X,Y,Z) (cm) on IP position
191 AliGenHijing *gener = new AliGenHijing(-1);
192 // centre of mass energy
193 gener->SetEnergyCMS(5500);
195 gener->SetReferenceFrame("CMS ");
197 gener->SetProjectile("A ", 208, 82);
198 gener->SetTarget ("A ", 208, 82);
199 // impact parameter range
200 gener->SetImpactParameterRange(0, 3.);
201 // evaluate cross section before run
202 gener->SetEvaluate(0);
203 // tell hijing to keep the full parent child chain
204 gener->KeepFullEvent();
205 // enable jet quenching
206 gener->SetJetQuenching(1);
208 gener->SetShadowing(1);
209 // neutral pion and heavy particle decays switched off
210 gener->SetDecaysOff(1);
212 gener->SetTrigger(0);
213 // kinematic selection
214 gener->SetSelectAll(0);
216 gener->SetMomentumRange(0,999);
218 gener->SetPhiRange(0.,360.);
220 gener->SetThetaRange(0,180.);
221 // select flavor (0: no, 4: charm+beauty, 5:beauty)
224 gener->SetOrigin(0., 0.0 ,0);
225 gener->SetSigma(0,0,5.3);
226 gener->SetVertexSmear(kPerEvent);
228 gener->SetTrackingFlag(0);
234 //********************************************
235 // Example for Charm Production with Pythia *
236 //********************************************
238 AliGenPythia *gener = new AliGenPythia(-1);
239 // final state kinematic cuts
240 gener->SetMomentumRange(0,999);
241 gener->SetPhiRange(0. ,360.);
242 gener->SetThetaRange(0., 180.);
243 gener->SetYRange(-10,10);
244 gener->SetPtRange(0,100);
245 // vertex position and smearing
246 gener->SetOrigin(0,0,0); // vertex position
247 gener->SetVertexSmear(kPerEvent);
248 gener->SetSigma(0,0,5.6); // Sigma in (X,Y,Z) (cm) on IP position
249 // Structure function. See the list in EVGEN/AliStructFuncType.h
250 gener->SetStrucFunc(kGRVHO);
251 // Select corection for nuclear structure functions
252 // gener->SetNuclei(208,208);
254 // Process type. See the list in PYTHIA6/AliPythia.h
255 gener->SetProcess(kPyBeauty);
257 // Pt transfer of the hard scattering
258 gener->SetPtHard(0.,5.);
259 // Decay type (semielectronic, semimuonic, nodecay)
260 gener->SetForceDecay(kSemiElectronic);
261 // Centre of mass energy
262 gener->SetEnergyCMS(5500.);
264 gener->SetTrackingFlag(0);
270 //********************************************
271 // Example for Charm Production with Pythia *
272 //********************************************
274 AliGenHerwig *gener = new AliGenHerwig(-1);
275 // final state kinematic cuts
276 gener->SetMomentumRange(0,7000);
277 gener->SetPhiRange(0. ,360.);
278 gener->SetThetaRange(0., 180.);
279 gener->SetYRange(-10,10);
280 gener->SetPtRange(0,7000);
281 // vertex position and smearing
282 gener->SetOrigin(0,0,0); // vertex position
283 gener->SetVertexSmear(kPerEvent);
284 gener->SetSigma(0,0,5.6); // Sigma in (X,Y,Z) (cm) on IP position
286 gener->SetBeamMomenta(7000,7000);
288 gener->SetProjectile("P");
289 gener->SetTarget("P");
290 // Structure function
291 gener->SetStrucFunc(kGRVHO);
293 gener->SetPtHardMin(200);
296 gener->SetProcess(8000);
298 gener->SetTrackingFlag(0);
303 AliGenIsajet *gener = new AliGenIsajet(-1);
307 AliGenDPMjet *gener = new AliGenDPMjet(-1);
308 //kDpmSingleDiffr, kDpmDoubleDiffr, kDpmDiffr, kDpmNonDiffr
309 gener->SetProcess(kDpmMb);
316 //*******************************************************
317 // Example for J/psi Production from Parameterisation
318 // using default library (AliMUONlib)
319 //*******************************************************
322 new AliGenParam(ntracks, AliGenMUONlib::kUpsilon);
323 gener->SetMomentumRange(0,999);
324 gener->SetPtRange(0,999);
325 gener->SetPhiRange(0. , 360.);
326 gener->SetYRange(2.5,4);
327 gener->SetCutOnChild(1);
328 gener->SetChildThetaRange(2,9);
329 gener->SetOrigin(0,0,0); //vertex position
330 gener->SetSigma(0,0,5.3); //Sigma in (X,Y,Z) (cm) on IP position
331 gener->SetForceDecay(kDiMuon);
332 gener->SetTrackingFlag(0);
338 //*******************************************************
339 // Example for Omega Production from Parameterisation
340 // specifying library.
341 //*******************************************************
343 AliGenParam *gener = new AliGenParam(1000,new AliGenPHOSlib(),
344 AliGenPHOSlib::kOmega);
345 gener->SetWeighting(kNonAnalog);
346 gener->SetForceDecay(kNoDecay);
347 gener->SetPtRange(0,100);
348 gener->SetThetaRange(45,135);
349 gener->SetTrackingFlag(0);
355 //*******************************************************
356 // Example for Upsilon Production from Parameterisation
357 // specifying library.
359 //*******************************************************
361 AliGenParam *gener = new AliGenParam(1000,new AliGenGSIlib(),
362 AliGenGSIlib::kUpsilon, "MUON");
363 gener->SetMomentumRange(0,999);
364 gener->SetPtRange(0,999);
365 gener->SetPhiRange(0., 360.);
366 gener->SetYRange(2.5,4);
367 gener->SetCutOnChild(1);
368 gener->SetChildThetaRange(2,9);
369 gener->SetOrigin(0,0,0); //vertex position
370 gener->SetSigma(0,0,5.3); //Sigma in (X,Y,Z) (cm) on IP position
371 gener->SetForceDecay(kDiMuon);
372 gener->SetTrackingFlag(0);
378 //*******************************************************
379 // Example for Omega Production from Parameterisation
380 // specifying library.
381 // The alternative way.
382 //*******************************************************
384 AliGenLib* Lib=new AliGenPHOSlib();
385 Int_t iOmega = AliGenPHOSlib::kOmega;
386 AliGenParam *gener = new AliGenParam(50, iOmega,
387 Lib->GetPt(iOmega, ""),
388 Lib->GetY (iOmega, ""),
389 Lib->GetIp(iOmega, ""));
390 gener->SetPtRange(0,999);
391 gener->SetWeighting(kNonAnalog);
392 gener->SetForceDecay(kNoDecay);
393 gener->SetTrackingFlag(0);
399 //*******************************************************
400 // Example for a FLUKA Boundary Source *
401 //*******************************************************
403 AliGenFLUKAsource *gener = new AliGenFLUKAsource(-1);
404 gener->SetFileName("$(ALICE_ROOT)/data/all32.root");
405 gener->SetPartFlag(9);
406 gener->SetAgeMax(1.e-5);
408 gener->SetFraction(0.0315);
409 // gener->SetFraction(0.75*0.0315);
411 // gener->SetPartFlag(10);
412 gener->SetMomentumRange(0,999);
413 gener->SetPhiRange(0.,360.);
414 gener->SetThetaRange(0., 180.);
415 gener->SetAgeMax(1.e-5);
418 // gener->SetFraction(0.0315);
424 //*******************************************************
425 // Example for reading from a external file *
426 //*******************************************************
428 AliGenExtFile *gener = new AliGenExtFile(-1);
429 gener->SetVertexSmear(kPerEvent);
430 gener->SetTrackingFlag(1);
432 AliGenReaderTreeK * reader = new AliGenReaderTreeK();
433 reader->SetFileName("$(ALICE_ROOT)/data/dtujet93.root");
434 gener->SetReader(reader);
440 //*******************************************************
441 // Example for Tunnel Halo Source *
442 //*******************************************************
444 AliGenHalo *gener = new AliGenHalo(ntracks);
445 gener->SetFileName("/h1/morsch/marsip/marsip5.mu");
451 //*******************************************************
452 // Example for a Cocktail *
453 //*******************************************************
455 AliGenCocktail *gener = new AliGenCocktail();
457 gener->SetPhiRange(0,360);
458 gener->SetYRange(2.5,4);
459 gener->SetThetaRange(2,9);
460 gener->SetPtRange(0,10);
461 gener->SetOrigin(0,0,0); //vertex position
462 gener->SetSigma(0,0,0); //Sigma in (X,Y,Z) (cm) on IP position
463 gener->SetMomentumRange(0,999);
465 AliGenParam *jpsi = new AliGenParam(1,AliGenMUONlib::kJpsi);
466 jpsi->SetForceDecay(kDiMuon);
467 jpsi->SetCutOnChild(1);
470 AliGenFLUKAsource *bg = new AliGenFLUKAsource(-1);
471 bg->AddFile("$(ALICE_ROOT)/data/all32.root");
474 bg->SetAgeMax(1.e-5);
476 // gener->SetFraction(0.0315);
477 bg->SetFraction(0.01*0.0315);
479 gener->AddGenerator(jpsi,"J/Psi", 1);
480 gener->AddGenerator(bg,"Background",1);
487 // Activate this line if you want the vertex smearing to happen
490 // gener->SetVertexSmear(kPerTrack);
494 gAlice->SetField(-999,2); //Specify maximum magnetic field in Tesla (neg. ==> default field)
499 //=================== Alice BODY parameters =============================
500 AliBODY *BODY = new AliBODY("BODY","Alice envelop");
504 //=================== MAG parameters ============================
505 // --- Start with Magnet since detector layouts may be depending ---
506 // --- on the selected Magnet dimensions ---
507 AliMAG *MAG = new AliMAG("MAG","Magnet");