]> git.uio.no Git - u/mrichter/AliRoot.git/blame - EMCAL/Config.C
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Config.C for EMCAL production.
[u/mrichter/AliRoot.git] / EMCAL / Config.C
CommitLineData
b7d89bab 1enum PprRun_t
2{
3 kTest50,
4 kParam_8000, kParam_4000, kParam_2000,
5 kHijing_cent1, kHijing_cent2,
6 kHijing_per1, kHijing_per2, kHijing_per3, kHijing_per4, kHijing_per5,
7 kHijing_jj25, kHijing_jj50, kHijing_jj75, kHijing_jj100, kHijing_jj125,
8 kHijing_gj25, kHijing_gj50, kHijing_gj75, kHijing_gj100, kHijing_gj125,
9 kJetPlusBg, kGammaPlusBg, kParam_8000_Ecal,
10 kJets_100,
11 kGammaGun, kGammaBox
12};
13
14enum PprGeo_t
15{
16 kHoles, kNoHoles
17};
18
19enum PprRad_t
20{
21 kGluonRadiation, kNoGluonRadiation
22};
23
24// This part for configuration
25//static PprRun_t run = kParam_8000_Ecal;
26static PprRun_t run = kJets_100;
27static PprGeo_t geo = kHoles;
28static PprRad_t rad = kNoGluonRadiation;
29// Comment line
30static TString comment;
31
32void Config()
33{
34
35 // 7-DEC-2000 09:00
36 // Switch on Transition adiation simulation. 6/12/00 18:00
37 // iZDC=1 7/12/00 09:00
38 // ThetaRange is (0., 180.). It was (0.28,179.72) 7/12/00 09:00
39 // Theta range given through pseudorapidity limits 22/6/2001
40
41 // Set Random Number seed
42 // gRandom->SetSeed(12345);
43
44 new AliGeant3("C++ Interface to Geant3");
45
46 if (!gSystem->Getenv("CONFIG_FILE"))
47 {
48 TFile *rootfile = new TFile("galice.root", "recreate");
49
50 rootfile->SetCompressionLevel(2);
51 }
52
53 TGeant3 *geant3 = (TGeant3 *) gMC;
54
55 //
56 // Set External decayer
57 AliDecayer *decayer = new AliDecayerPythia();
58
59 decayer->SetForceDecay(kAll);
60 decayer->Init();
61 gMC->SetExternalDecayer(decayer);
62 //
63 //
64 //=======================================================================
65 // ******* GEANT STEERING parameters FOR ALICE SIMULATION *******
66 geant3->SetTRIG(1); //Number of events to be processed
67 geant3->SetSWIT(4, 10);
68 geant3->SetDEBU(0, 0, 1);
69 geant3->SetDCAY(1);
70 geant3->SetPAIR(1);
71 geant3->SetCOMP(1);
72 geant3->SetPHOT(1);
73 geant3->SetPFIS(0);
74 geant3->SetDRAY(0);
75 geant3->SetANNI(1);
76 geant3->SetBREM(1);
77 geant3->SetMUNU(1);
78 geant3->SetCKOV(1);
79 geant3->SetHADR(1); //Select pure GEANH (HADR 1) or GEANH/NUCRIN (HADR 3)
80 geant3->SetLOSS(2);
81 geant3->SetMULS(1);
82 geant3->SetRAYL(1);
83 geant3->SetAUTO(1); //Select automatic STMIN etc... calc. (AUTO 1) or manual (AUTO 0)
84 geant3->SetABAN(0); //Restore 3.16 behaviour for abandoned tracks
85 geant3->SetOPTI(2); //Select optimisation level for GEANT geometry searches (0,1,2)
86 geant3->SetERAN(5.e-7);
87
88 Float_t cut = 1.e-3; // 1MeV cut by default
89 Float_t tofmax = 1.e10;
90
91 // GAM ELEC NHAD CHAD MUON EBREM MUHAB EDEL MUDEL MUPA TOFMAX
92 geant3->SetCUTS(cut, cut, cut, cut, cut, cut, cut, cut, cut, cut, tofmax);
93 gAlice->TrackingLimits(900, 900);
94
95 //
96 //=======================================================================
97 // ************* STEERING parameters FOR ALICE SIMULATION **************
98 // --- Specify event type to be tracked through the ALICE setup
99 // --- All positions are in cm, angles in degrees, and P and E in GeV
100
101
102// Generator Configuration
103 gAlice->SetDebug(1);
104
105 AliGenerator* gener = GeneratorFactory(run);
106 gener->Init();
107 gener->SetTrackingFlag(1);
108
109
110
111 if (rad == kGluonRadiation)
112 {
113 //coment= comment.Append(" | Gluon Radiation On");
114
115 } else {
116 //coment= comment.Append(" | Gluon Radiation Off");
117 }
118
119 if (geo == kHoles)
120 {
121 //coment= comment.Append(" | Holes for PHOS/RICH");
122
123 } else {
124 //coment= comment.Append(" | No holes for PHOS/RICH");
125 }
126
127 printf("\n \n Comment: %s \n \n", (char*) comment);
128
129
130// Field (L3 0.4 T)
131
132 AliMagFCM* field = new AliMagFCM(
133 "Map2","$(ALICE_ROOT)/data/field01.dat", 2, 1., 10.);
134 field->SetSolenoidField(4.);
135 gAlice->SetField(field);
136
137 Int_t iMAG = 0;
138 Int_t iITS = 0;
139 Int_t iTPC = 0;
140 Int_t iTOF = 0;
141 Int_t iRICH = 0;
142 Int_t iZDC = 0;
143 Int_t iCASTOR = 0;
144 Int_t iTRD = 0;
145 Int_t iABSO = 0;
146 Int_t iDIPO = 0;
147 Int_t iHALL = 0;
148 Int_t iFRAME = 0;
149 Int_t iSHIL = 0;
150 Int_t iPIPE = 0;
151 Int_t iFMD = 0;
152 Int_t iMUON = 0;
153 Int_t iPHOS = 0;
154 Int_t iPMD = 0;
155 Int_t iSTART = 0;
156 Int_t iVZERO = 0;
157 Int_t iEMCAL = 1;
158//
159// enable/disable StepManager()
160//
161
162 Int_t enableABSO = 0;
163 Int_t enableCASTOR = 0;
164 Int_t enableDIPO = 0;
165 Int_t enableFMD = 0;
166 Int_t enableFRAME = 0;
167 Int_t enableHALL = 0;
168 Int_t enableITS = 0;
169 Int_t enableMAG = 0;
170 Int_t enableMUON = 0;
171 Int_t enablePHOS = 1;
172 Int_t enablePIPE = 0;
173 Int_t enablePMD = 0;
174 Int_t enableRICH = 0;
175 Int_t enableSHIL = 0;
176 Int_t enableSTART = 0;
177 Int_t enableTOF = 0;
178 Int_t enableTPC = 0;
179 Int_t enableTRD = 0;
180 Int_t enableZDC = 0;
181 Int_t enableEMCAL = 1;
182
183 //=================== Alice BODY parameters =============================
184 AliBODY *BODY = new AliBODY("BODY", "Alice envelop");
185
186
187 if (iMAG)
188 {
189 //=================== MAG parameters ============================
190 // --- Start with Magnet since detector layouts may be depending ---
191 // --- on the selected Magnet dimensions ---
192 AliMAG *MAG = new AliMAG("MAG", "Magnet");
193 }
194
195
196 if (iABSO)
197 {
198 //=================== ABSO parameters ============================
199 AliABSO *ABSO = new AliABSOv0("ABSO", "Muon Absorber");
200 }
201
202 if (iDIPO)
203 {
204 //=================== DIPO parameters ============================
205
206 AliDIPO *DIPO = new AliDIPOv2("DIPO", "Dipole version 2");
207 }
208
209 if (iHALL)
210 {
211 //=================== HALL parameters ============================
212
213 AliHALL *HALL = new AliHALL("HALL", "Alice Hall");
214 }
215
216
217 if (iFRAME)
218 {
219 //=================== FRAME parameters ============================
220
221 AliFRAMEv2 *FRAME = new AliFRAMEv2("FRAME", "Space Frame");
222 if (geo == kHoles) {
223 FRAME->SetHoles(1);
224 } else {
225 FRAME->SetHoles(0);
226 }
227 }
228
229 if (iSHIL)
230 {
231 //=================== SHIL parameters ============================
232
233 AliSHIL *SHIL = new AliSHILv0("SHIL", "Shielding");
234 }
235
236
237 if (iPIPE)
238 {
239 //=================== PIPE parameters ============================
240
241 AliPIPE *PIPE = new AliPIPEv0("PIPE", "Beam Pipe");
242 }
243
244 if(iITS) {
245
246//=================== ITS parameters ============================
247 //
248 // As the innermost detector in ALICE, the Inner Tracking System "impacts" on
249 // almost all other detectors. This involves the fact that the ITS geometry
250 // still has several options to be followed in parallel in order to determine
251 // the best set-up which minimizes the induced background. All the geometries
252 // available to date are described in the following. Read carefully the comments
253 // and use the default version (the only one uncommented) unless you are making
254 // comparisons and you know what you are doing. In this case just uncomment the
255 // ITS geometry you want to use and run Aliroot.
256 //
257 // Detailed geometries:
258 //
259 //
260 //AliITS *ITS = new AliITSv5symm("ITS","Updated ITS TDR detailed version with symmetric services");
261 //
262 //AliITS *ITS = new AliITSv5asymm("ITS","Updates ITS TDR detailed version with asymmetric services");
263 //
264 AliITSvPPRasymm *ITS = new AliITSvPPRasymm("ITS","New ITS PPR detailed version with asymmetric services");
265 ITS->SetMinorVersion(2); // don't touch this parameter if you're not an ITS developer
266 ITS->SetReadDet(kFALSE); // don't touch this parameter if you're not an ITS developer
267 // ITS->SetWriteDet("$ALICE_ROOT/ITS/ITSgeometry_vPPRasymm2.det"); // don't touch this parameter if you're not an ITS developer
268 ITS->SetThicknessDet1(200.); // detector thickness on layer 1 must be in the range [100,300]
269 ITS->SetThicknessDet2(200.); // detector thickness on layer 2 must be in the range [100,300]
270 ITS->SetThicknessChip1(200.); // chip thickness on layer 1 must be in the range [150,300]
271 ITS->SetThicknessChip2(200.); // chip thickness on layer 2 must be in the range [150,300]
272 ITS->SetRails(1); // 1 --> rails in ; 0 --> rails out
273 ITS->SetCoolingFluid(1); // 1 --> water ; 0 --> freon
274 //
275 //AliITSvPPRsymm *ITS = new AliITSvPPRsymm("ITS","New ITS PPR detailed version with symmetric services");
276 //ITS->SetMinorVersion(2); // don't touch this parameter if you're not an ITS developer
277 //ITS->SetReadDet(kFALSE); // don't touch this parameter if you're not an ITS developer
278 //ITS->SetWriteDet("$ALICE_ROOT/ITS/ITSgeometry_vPPRsymm2.det"); // don't touch this parameter if you're not an ITS developer
279 //ITS->SetThicknessDet1(200.); // detector thickness on layer 1 must be in the range [100,300]
280 //ITS->SetThicknessDet2(200.); // detector thickness on layer 2 must be in the range [100,300]
281 //ITS->SetThicknessChip1(200.); // chip thickness on layer 1 must be in the range [150,300]
282 //ITS->SetThicknessChip2(200.); // chip thickness on layer 2 must be in the range [150,300]
283 //ITS->SetRails(1); // 1 --> rails in ; 0 --> rails out
284 //ITS->SetCoolingFluid(1); // 1 --> water ; 0 --> freon
285 //
286 //
287 // Coarse geometries (warning: no hits are produced with these coarse geometries and they unuseful
288 // for reconstruction !):
289 //
290 //
291 //AliITSvPPRcoarseasymm *ITS = new AliITSvPPRcoarseasymm("ITS","New ITS PPR coarse version with asymmetric services");
292 //ITS->SetRails(1); // 1 --> rails in ; 0 --> rails out
293 //ITS->SetSupportMaterial(0); // 0 --> Copper ; 1 --> Aluminum ; 2 --> Carbon
294 //
295 //AliITS *ITS = new AliITSvPPRcoarsesymm("ITS","New ITS PPR coarse version with symmetric services");
296 //ITS->SetRails(1); // 1 --> rails in ; 0 --> rails out
297 //ITS->SetSupportMaterial(0); // 0 --> Copper ; 1 --> Aluminum ; 2 --> Carbon
298 //
299 //
300 //
301 // Geant3 <-> EUCLID conversion
302 // ============================
303 //
304 // SetEUCLID is a flag to output (=1) or not to output (=0) both geometry and
305 // media to two ASCII files (called by default ITSgeometry.euc and
306 // ITSgeometry.tme) in a format understandable to the CAD system EUCLID.
307 // The default (=0) means that you dont want to use this facility.
308 //
309 ITS->SetEUCLID(0);
310 if (!enableITS) ITS->DisableStepManager();
311 }
312
313
314 if (iTPC)
315 {
316 //============================ TPC parameters ================================
317 // --- This allows the user to specify sectors for the SLOW (TPC geometry 2)
318 // --- Simulator. SecAL (SecAU) <0 means that ALL lower (upper)
319 // --- sectors are specified, any value other than that requires at least one
320 // --- sector (lower or upper)to be specified!
321 // --- Reminder: sectors 1-24 are lower sectors (1-12 -> z>0, 13-24 -> z<0)
322 // --- sectors 25-72 are the upper ones (25-48 -> z>0, 49-72 -> z<0)
323 // --- SecLows - number of lower sectors specified (up to 6)
324 // --- SecUps - number of upper sectors specified (up to 12)
325 // --- Sens - sensitive strips for the Slow Simulator !!!
326 // --- This does NOT work if all S or L-sectors are specified, i.e.
327 // --- if SecAL or SecAU < 0
328 //
329 //
330 //-----------------------------------------------------------------------------
331
332 // gROOT->LoadMacro("SetTPCParam.C");
333 // AliTPCParam *param = SetTPCParam();
334 AliTPC *TPC = new AliTPCv2("TPC", "Default");
335
336 // All sectors included
337 TPC->SetSecAL(-1);
338 TPC->SetSecAU(-1);
339 if (!enableTPC) TPC->DisableStepManager();
340 }
341
342
343 if (iRICH)
344 {
345 //=================== RICH parameters ===========================
346 AliRICH *RICH = new AliRICHv1("RICH", "normal RICH");
347 if (!enableRICH) RICH->DisableStepManager();
348 }
349
350
351 if (iZDC)
352 {
353 //=================== ZDC parameters ============================
354
355 AliZDC *ZDC = new AliZDCv2("ZDC", "normal ZDC");
356 if (!enableZDC) ZDC->DisableStepManager();
357 }
358
359 if (iCASTOR)
360 {
361 //=================== CASTOR parameters ============================
362
363 AliCASTOR *CASTOR = new AliCASTORv1("CASTOR", "normal CASTOR");
364 if (!enableCASTOR) CASTOR->DisableStepManager();
365 }
366
367 if (iTRD)
368 {
369 //=================== TRD parameters ============================
370
371 AliTRD *TRD = new AliTRDv1("TRD", "TRD slow simulator");
372
373 // Select the gas mixture (0: 97% Xe + 3% isobutane, 1: 90% Xe + 10% CO2)
374 TRD->SetGasMix(1);
375 if (geo == kHoles) {
376 // With hole in front of PHOS
377 TRD->SetPHOShole();
378 // With hole in front of RICH
379 TRD->SetRICHhole();
380 }
381 // Switch on TR
382 AliTRDsim *TRDsim = TRD->CreateTR();
383 if (!enableTRD) TRD->DisableStepManager();
384 }
385
386 if (iFMD)
387 {
388 //=================== FMD parameters ============================
389
390 AliFMD *FMD = new AliFMDv1("FMD", "normal FMD");
391 FMD->SetRingsSi1(256);
392 FMD->SetRingsSi2(64);
393 FMD->SetSectorsSi1(20);
394 FMD->SetSectorsSi2(24);
395 if (!enableFMD) FMD->DisableStepManager();
396 }
397
398 if (iMUON)
399 {
400 //=================== MUON parameters ===========================
401
402 AliMUON *MUON = new AliMUONv1("MUON", "default");
403 if (!enableMUON) MUON->DisableStepManager();
404 }
405 //=================== PHOS parameters ===========================
406
407 if (iPHOS)
408 {
409 AliPHOS *PHOS = new AliPHOSv1("PHOS", "IHEP");
410 if (!enablePHOS) PHOS->DisableStepManager();
411 }
412
413
414 if (iPMD)
415 {
416 //=================== PMD parameters ============================
417 AliPMD *PMD = new AliPMDv1("PMD", "normal PMD");
418 PMD->SetPAR(1., 1., 0.8, 0.02);
419 PMD->SetIN(6., 18., -580., 27., 27.);
420 PMD->SetGEO(0.0, 0.2, 4.);
421 PMD->SetPadSize(0.8, 1.0, 1.0, 1.5);
422 if (!enablePMD) PMD->DisableStepManager();
423 }
424
425 if (iTOF) {
426 if (geo == kHoles) {
427 //=================== TOF parameters ============================
428 AliTOF *TOF = new AliTOFv2FHoles("TOF", "TOF with Holes");
429 } else {
430 AliTOF *TOF = new AliTOFv4("TOF", "normal TOF");
431 }
432 if (!enableTOF) TOF->DisableStepManager();
433 }
434
435 if (iSTART)
436 {
437 //=================== START parameters ============================
438 AliSTART *START = new AliSTARTv1("START", "START Detector");
439 if (!enableSTART) START->DisableStepManager();
440 }
441
442 if (iEMCAL && !iRICH)
443 {
444 //=================== EMCAL parameters ============================
445 AliEMCAL *EMCAL = new AliEMCALv1("EMCAL", "EMCALArch1b");
446 if (!enableEMCAL) EMCAL->DisableStepManager();
447 }
448}
449
450Float_t EtaToTheta(Float_t arg){
451 return (180./TMath::Pi())*2.*atan(exp(-arg));
452}
453
454
455
456AliGenerator* GeneratorFactory(PprRun_t run) {
457 Int_t isw = 3;
458 if (rad == kNoGluonRadiation) isw = 0;
459
460
461 switch (run) {
462 case kGammaGun:
463 gener = new AliGenFixed(1);
464 gener->SetMomentum(100.);
465 gener->SetPhi(240.);
466 gener->SetTheta(91.);
467 gener->SetPart(kElectron);
468 break;
469
470 case kGammaBox:
471 gener = new AliGenBox(100);
472 gener->SetMomentumRange(100., 100.1);
473 gener->SetPhiRange(0,360);
474 gener->SetThetaRange(45., 135.);
475 gener->SetPart(kElectron);
476 break;
477
478 case kTest50:
479 //coment= comment.Append(":HIJINGparam test 50 particles");
480 gener = new AliGenHIJINGparaBa(50);
481 gener->SetMomentumRange(0, 999999.);
482 gener->SetPhiRange(-180., 180.);
483 // Set pseudorapidity range from -8 to 8.
484 Float_t thmin = EtaToTheta(8); // theta min. <---> eta max
485 Float_t thmax = EtaToTheta(-8); // theta max. <---> eta min
486 gener->SetThetaRange(thmin,thmax);
487 break;
488
489 case kParam_8000:
490 //coment= comment.Append(":HIJINGparam N=8000");
491 AliGenHIJINGpara *gener = new AliGenHIJINGpara(86030);
492 gener->SetMomentumRange(0, 999999.);
493 gener->SetPhiRange(-180., 180.);
494 // Set pseudorapidity range from -8 to 8.
495 Float_t thmin = EtaToTheta(8); // theta min. <---> eta max
496 Float_t thmax = EtaToTheta(-8); // theta max. <---> eta min
497 gener->SetThetaRange(thmin,thmax);
498 break;
499 case kParam_4000:
500 //coment= comment.Append("HIJINGparam N=4000");
501 AliGenHIJINGpara *gener = new AliGenHIJINGpara(43015);
502 gener->SetMomentumRange(0, 999999.);
503 gener->SetPhiRange(-180., 180.);
504 // Set pseudorapidity range from -8 to 8.
505 Float_t thmin = EtaToTheta(8); // theta min. <---> eta max
506 Float_t thmax = EtaToTheta(-8); // theta max. <---> eta min
507 gener->SetThetaRange(thmin,thmax);
508 break;
509 case kParam_2000:
510 //coment= comment.Append("HIJINGparam N=2000");
511 AliGenHIJINGpara *gener = new AliGenHIJINGpara(21507);
512 gener->SetMomentumRange(0, 999999.);
513 gener->SetPhiRange(-180., 180.);
514 // Set pseudorapidity range from -8 to 8.
515 Float_t thmin = EtaToTheta(8); // theta min. <---> eta max
516 Float_t thmax = EtaToTheta(-8); // theta max. <---> eta min
517 gener->SetThetaRange(thmin,thmax);
518 break;
519 case kParam_8000_Ecal:
520 //coment= comment.Append(":HIJINGparam N=8000 baryons");
521 AliGenHIJINGparaBa *gener = new AliGenHIJINGparaBa(82534);
522 gener->SetMomentumRange(0, 999999.);
523 gener->SetPhiRange(-180., 180.);
524 // Set pseudorapidity range from -8 to 8.
525 Float_t thmin = EtaToTheta( 5); // theta min. <---> eta max
526 Float_t thmax = EtaToTheta(-5); // theta max. <---> eta min
527 gener->SetThetaRange(thmin,thmax);
528 break;
529//
530// Hijing Central
531//
532 case kHijing_cent1:
533 //coment= comment.Append("HIJING cent1");
534 AliGenHijing *gener = HijingStandard();
535// impact parameter range
536 gener->SetImpactParameterRange(0., 5.);
537 break;
538 case kHijing_cent2:
539 //coment= comment.Append("HIJING cent2");
540 AliGenHijing *gener = HijingStandard();
541// impact parameter range
542 gener->SetImpactParameterRange(0., 2.);
543 break;
544//
545// Hijing Peripheral
546//
547 case kHijing_per1:
548 //coment= comment.Append("HIJING per1");
549 AliGenHijing *gener = HijingStandard();
550// impact parameter range
551 gener->SetImpactParameterRange(5., 8.6);
552 break;
553 case kHijing_per2:
554 //coment= comment.Append("HIJING per2");
555 AliGenHijing *gener = HijingStandard();
556// impact parameter range
557 gener->SetImpactParameterRange(8.6, 11.2);
558 break;
559 case kHijing_per3:
560 //coment= comment.Append("HIJING per3");
561 AliGenHijing *gener = HijingStandard();
562// impact parameter range
563 gener->SetImpactParameterRange(11.2, 13.2);
564 break;
565 case kHijing_per4:
566 //coment= comment.Append("HIJING per4");
567 AliGenHijing *gener = HijingStandard();
568// impact parameter range
569 gener->SetImpactParameterRange(13.2, 15.);
570 break;
571 case kHijing_per5:
572 //coment= comment.Append("HIJING per5");
573 AliGenHijing *gener = HijingStandard();
574// impact parameter range
575 gener->SetImpactParameterRange(15., 100.);
576 break;
577//
578// Jet-Jet
579//
580 case kHijing_jj25:
581 //coment= comment.Append("HIJING Jet 25 GeV");
582 AliGenHijing *gener = HijingStandard();
583// impact parameter range
584 gener->SetImpactParameterRange(0., 5.);
585 // trigger
586 gener->SetTrigger(1);
587 gener->SetPtMinJet(25.);
588 gener->SetRadiation(isw);
589 gener->SetJetEtaRange(-0.3,0.3);
590 gener->SetJetPhiRange(15.,105.);
591 break;
592
593 case kHijing_jj50:
594 //coment= comment.Append("HIJING Jet 50 GeV");
595 AliGenHijing *gener = HijingStandard();
596// impact parameter range
597 gener->SetImpactParameterRange(0., 5.);
598 // trigger
599 gener->SetTrigger(1);
600 gener->SetPtMinJet(50.);
601 gener->SetRadiation(isw);
602 gener->SetJetEtaRange(-0.3,0.3);
603 gener->SetJetPhiRange(15.,105.);
604 break;
605
606 case kHijing_jj75:
607 //coment= comment.Append("HIJING Jet 75 GeV");
608 AliGenHijing *gener = HijingStandard();
609// impact parameter range
610 gener->SetImpactParameterRange(0., 5.);
611 // trigger
612 gener->SetTrigger(1);
613 gener->SetPtMinJet(75.);
614 gener->SetRadiation(isw);
615 gener->SetJetEtaRange(-0.3,0.3);
616 gener->SetJetPhiRange(15.,105.);
617 break;
618
619 case kHijing_jj100:
620 //coment= comment.Append("HIJING Jet 100 GeV");
621 AliGenHijing *gener = HijingStandard();
622// impact parameter range
623 gener->SetImpactParameterRange(10., 10.1);
624 // trigger
625 gener->SetTrigger(1);
626 gener->SetPtMinJet(100.);
627 gener->SetRadiation(isw);
628 gener->SetJetEtaRange(-0.3,0.3);
629 gener->SetJetPhiRange(15.,105.);
630 break;
631
632 case kHijing_jj125:
633 //coment= comment.Append("HIJING Jet 125 GeV");
634 AliGenHijing *gener = HijingStandard();
635// impact parameter range
636 gener->SetImpactParameterRange(0., 5.);
637 // trigger
638 gener->SetTrigger(1);
639 gener->SetPtMinJet(125.);
640 gener->SetRadiation(isw);
641 gener->SetJetEtaRange(-0.3,0.3);
642 gener->SetJetPhiRange(15.,105.);
643 break;
644//
645// Gamma-Jet
646//
647 case kHijing_gj25:
648 //coment= comment.Append("HIJING Gamma 25 GeV");
649 AliGenHijing *gener = HijingStandard();
650// impact parameter range
651 gener->SetImpactParameterRange(0., 5.);
652 // trigger
653 gener->SetTrigger(2);
654 gener->SetPtMinJet(25.);
655 gener->SetRadiation(isw);
656 gener->SetJetEtaRange(-0.3,0.3);
657 gener->SetJetPhiRange(15.,105.);
658 break;
659
660 case kHijing_gj50:
661 //coment= comment.Append("HIJING Gamma 50 GeV");
662 AliGenHijing *gener = HijingStandard();
663// impact parameter range
664 gener->SetImpactParameterRange(0., 5.);
665 // trigger
666 gener->SetTrigger(2);
667 gener->SetPtMinJet(50.);
668 gener->SetRadiation(isw);
669 gener->SetJetEtaRange(-0.3,0.3);
670 gener->SetJetPhiRange(15.,105.);
671 break;
672
673 case kHijing_gj75:
674 //coment= comment.Append("HIJING Gamma 75 GeV");
675 AliGenHijing *gener = HijingStandard();
676// impact parameter range
677 gener->SetImpactParameterRange(0., 5.);
678 // trigger
679 gener->SetTrigger(2);
680 gener->SetPtMinJet(75.);
681 gener->SetRadiation(isw);
682 gener->SetJetEtaRange(-0.3,0.3);
683 gener->SetJetPhiRange(15.,105.);
684 break;
685
686 case kHijing_gj100:
687 //coment= comment.Append("HIJING Gamma 100 GeV");
688 AliGenHijing *gener = HijingStandard();
689// impact parameter range
690 gener->SetImpactParameterRange(0., 5.);
691 // trigger
692 gener->SetTrigger(2);
693 gener->SetPtMinJet(100.);
694 gener->SetRadiation(isw);
695 gener->SetJetEtaRange(-0.3,0.3);
696 gener->SetJetPhiRange(15.,105.);
697 break;
698
699 case kHijing_gj125:
700 //coment= comment.Append("HIJING Gamma 125 GeV");
701 AliGenHijing *gener = HijingStandard();
702// impact parameter range
703 gener->SetImpactParameterRange(0., 5.);
704 // trigger
705 gener->SetTrigger(2);
706 gener->SetPtMinJet(125.);
707 gener->SetRadiation(isw);
708 gener->SetJetEtaRange(-0.3,0.3);
709 gener->SetJetPhiRange(15.,105.);
710 break;
711 case kJetPlusBg:
712 AliGenCocktail *gener = new AliGenCocktail();
713 gener->SetMomentumRange(0, 999999.);
714 gener->SetPhiRange(-180., 180.);
715 // Set pseudorapidity range from -8 to 8.
716 Float_t thmin = EtaToTheta( 5.); // theta min. <---> eta max
717 Float_t thmax = EtaToTheta(-5.); // theta max. <---> eta min
718 gener->SetThetaRange(thmin,thmax);
719
720//
721// Underlying Event
722//
723// AliGenHIJINGparaBa *bg = new AliGenHIJINGparaBa(82534);
724 AliGenHIJINGparaBa *bg = new AliGenHIJINGparaBa(10);
725//
726// Jets from Pythia
727//
728 AliGenPythia *jets = new AliGenPythia(-1);
729// Centre of mass energy
730 jets->SetEnergyCMS(5500.);
731// Process type
732 jets->SetProcess(kPyJets);
733// final state kinematic cuts
734 jets->SetJetEtaRange(-0.3, 0.3);
735 jets->SetJetPhiRange(15., 105.);
736// Structure function
737 jets->SetStrucFunc(kGRV_LO_98);
738//
739// Pt transfer of the hard scattering
740 jets->SetPtHard(100.,100.1);
741// Decay type (semielectronic, semimuonic, nodecay)
742 jets->SetForceDecay(kAll);
743//
744// Add all to cockail ...
745//
746 gener->AddGenerator(bg,"Underlying Event", 1);
747// gener->AddGenerator(jets,"Jets",1);
748 break;
749 case kGammaPlusBg:
750 AliGenCocktail *gener = new AliGenCocktail();
751 gener->SetMomentumRange(0, 999999.);
752 gener->SetPhiRange(-180., 180.);
753 // Set pseudorapidity range from -8 to 8.
754 Float_t thmin = EtaToTheta( 5.); // theta min. <---> eta max
755 Float_t thmax = EtaToTheta(-5.); // theta max. <---> eta min
756 gener->SetThetaRange(thmin,thmax);
757
758//
759// Underlying Event
760//
761 AliGenHIJINGparaBa *bg = new AliGenHIJINGparaBa(82534);
762//
763// Jets from Pythia
764//
765 AliGenPythia *jets = new AliGenPythia(-1);
766// Centre of mass energy
767 jets->SetEnergyCMS(5500.);
768// Process type
769 jets->SetProcess(kPyDirectGamma);
770// final state kinematic cuts
771 jets->SetJetEtaRange(-0.3, 0.3);
772 jets->SetJetPhiRange(15., 105.);
773 jets->SetGammaEtaRange(-0.12, 0.12);
774 jets->SetGammaPhiRange(220., 320.);
775// Structure function
776 jets->SetStrucFunc(kGRV_LO_98);
777//
778// Pt transfer of the hard scattering
779 jets->SetPtHard(100.,100.1);
780// Decay type (semielectronic, semimuonic, nodecay)
781 jets->SetForceDecay(kAll);
782//
783// Add all to cockail ...
784//
785 gener->AddGenerator(bg,"Underlying Event", 1);
786 gener->AddGenerator(jets,"Jets",1);
787 break;
788 case kJets_100:
789// 100 GeV Jets
790 AliGenPythia *gener = PythiaJets(100.);
791 break;
792 }
793
794 return gener;
795}
796
797AliGenHijing* HijingStandard()
798{
799 AliGenHijing *gener = new AliGenHijing(-1);
800// centre of mass energy
801 gener->SetEnergyCMS(5500.);
802// reference frame
803 gener->SetReferenceFrame("CMS");
804// projectile
805 gener->SetProjectile("A", 208, 82);
806 gener->SetTarget ("A", 208, 82);
807// tell hijing to keep the full parent child chain
808 gener->KeepFullEvent();
809// enable jet quenching
810 gener->SetJetQuenching(1);
811// enable shadowing
812 gener->SetShadowing(1);
813// neutral pion and heavy particle decays switched off
814 gener->SetDecaysOff(1);
815// Don't track spectators
816 gener->SetSpectators(0);
817// kinematic selection
818 gener->SetSelectAll(0);
819 return gener;
820}
821
822AliGenPythia* PythiaJets(Float_t energy)
823{
824 AliGenPythia *gener = new AliGenPythia(-1);
825// Centre of mass energy
826 gener->SetEnergyCMS(5500.);
827// Process type
828 gener->SetProcess(kPyJets);
829// final state kinematic cuts
830 gener->SetJetEtaRange(-0.3, 0.3);
831 gener->SetJetPhiRange(15., 105.);
832// Structure function
833 gener->SetStrucFunc(kGRV_LO_98);
834//
835// Pt transfer of the hard scattering
836 gener->SetPtHard(energy, energy+0.1);
837// Decay type (semielectronic, semimuonic, nodecay)
838 gener->SetForceDecay(kAll);
839//
840 return gener;
841}
842
843
844AliGenPythia* PythiaGamma(Float_t energy)
845{
846 AliGenPythia *gener = new AliGenPythia(-1);
847// Centre of mass energy
848 gener->SetEnergyCMS(5500.);
849// Process type
850 gener->SetProcess(kPyDirectGamma);
851// final state kinematic cuts
852 gener->SetJetEtaRange(-0.3, 0.3);
853 gener->SetJetPhiRange(15., 105.);
854 gener->SetGammaEtaRange(-0.12, 0.12);
855 gener->SetGammaPhiRange(220., 320.);
856// Structure function
857 gener->SetStrucFunc(kGRV_LO_98);
858//
859// Pt transfer of the hard scattering
860 gener->SetPtHard(energy, energy+0.1);
861// Decay type (semielectronic, semimuonic, nodecay)
862 gener->SetForceDecay(kAll);
863//
864 return gener;
865}
866
867
868
869
870
871
872
873
874
875
876