Compilation warnings
[u/mrichter/AliRoot.git] / FMD / Config.C
1 //____________________________________________________________________
2 //
3 // $Id$
4 //
5 // One can use the configuration macro in compiled mode by
6 // root [0] gSystem->Load("libgeant321");
7 // root [0] gSystem->SetIncludePath("-I$ROOTSYS/include -I$ALICE_ROOT/include\
8 //                   -I$ALICE_ROOT -I$ALICE/geant3/TGeant3");
9 // root [0] .x grun.C(1,"ConfigPPR.C++")
10 //
11 /** @file    Config.C
12     @author  Christian Holm Christensen <cholm@nbi.dk>
13     @date    Mon Mar 27 12:50:29 2006
14     @brief   Simulation configuration script
15 */
16
17 //____________________________________________________________________
18 // 
19 // Generator types 
20 //
21 enum EG_t {
22   test50,
23   kParam_8000,                  //
24   kParam_4000,                  //
25   kParam_2000,                  //
26   kParam_fmd,                   //
27   kHijing_cent1,                //
28   kHijing_cent2,                //
29   kHijing_per1,                 //
30   kHijing_per2,                 //
31   kHijing_per3,                 //
32   kHijing_per4,                 //
33   kHijing_per5,                 //
34   kHijing_jj25,                 //
35   kHijing_jj50,                 //
36   kHijing_jj75,                 //
37   kHijing_jj100,                //
38   kHijing_jj200,                //
39   kHijing_gj25,                 //
40   kHijing_gj50,                 //
41   kHijing_gj75,                 //
42   kHijing_gj100,                //
43   kHijing_gj200,                //
44   kHijing_pA,                   //
45   kPythia6,                     //
46   kPythia6Jets20_24,            //
47   kPythia6Jets24_29,            //
48   kPythia6Jets29_35,            //
49   kPythia6Jets35_42,            //
50   kPythia6Jets42_50,            //
51   kPythia6Jets50_60,            //
52   kPythia6Jets60_72,            //
53   kPythia6Jets72_86,            //
54   kPythia6Jets86_104,           //
55   kPythia6Jets104_125,          //
56   kPythia6Jets125_150,          //
57   kPythia6Jets150_180,          //
58   kD0PbPb5500,                  //
59   kCharmSemiElPbPb5500,         //
60   kBeautySemiElPbPb5500,        //
61   kCocktailTRD,                 //
62   kPyJJ,                        //
63   kPyGJ,                        //
64   kMuonCocktailCent1,           //
65   kMuonCocktailPer1,            //
66   kMuonCocktailPer4,            //
67   kMuonCocktailCent1HighPt,     //
68   kMuonCocktailPer1HighPt,      //
69   kMuonCocktailPer4HighPt,      //
70   kMuonCocktailCent1Single,     //
71   kMuonCocktailPer1Single,      //
72   kMuonCocktailPer4Single,
73   kFMD1Flat, 
74   kFMD2Flat, 
75   kFMD3Flat,
76   kFMDFlat,
77   kEgMax
78 };
79
80 //____________________________________________________________________
81 // 
82 // Generator types names
83 //
84 const char* egName[kEgMax] = {
85   "test50",
86   "kParam_8000",                //
87   "kParam_4000",                //
88   "kParam_2000",                //
89   "kParam_fmd",                 //
90   "kHijing_cent1",              //
91   "kHijing_cent2",              //
92   "kHijing_per1",               //
93   "kHijing_per2",               //
94   "kHijing_per3",               //
95   "kHijing_per4",               //
96   "kHijing_per5",               //
97   "kHijing_jj25",               //
98   "kHijing_jj50",               //
99   "kHijing_jj75",               //
100   "kHijing_jj100",              //
101   "kHijing_jj200",              //
102   "kHijing_gj25",               //
103   "kHijing_gj50",               //
104   "kHijing_gj75",               //
105   "kHijing_gj100",              //
106   "kHijing_gj200",              //
107   "kHijing_pA",                 //
108   "kPythia6",                   //
109   "kPythia6Jets20_24",          //
110   "kPythia6Jets24_29",          //
111   "kPythia6Jets29_35",          //
112   "kPythia6Jets35_42",          //
113   "kPythia6Jets42_50",          //
114   "kPythia6Jets50_60",          //
115   "kPythia6Jets60_72",          //
116   "kPythia6Jets72_86",          //
117   "kPythia6Jets86_104",         //
118   "kPythia6Jets104_125",        //
119   "kPythia6Jets125_150",        //
120   "kPythia6Jets150_180",        //
121   "kD0PbPb5500",                //
122   "kCharmSemiElPbPb5500",       //
123   "kBeautySemiElPbPb5500",      //
124   "kCocktailTRD",               //
125   "kPyJJ",                      //
126   "kPyGJ",                      //
127   "kMuonCocktailCent1",         //
128   "kMuonCocktailPer1",          //
129   "kMuonCocktailPer4",          //
130   "kMuonCocktailCent1HighPt",   //
131   "kMuonCocktailPer1HighPt",    //
132   "kMuonCocktailPer4HighPt",    //
133   "kMuonCocktailCent1Single",   //
134   "kMuonCocktailPer1Single",    //
135   "kMuonCocktailPer4Single",
136   "kFMD1Flat",
137   "kFMD2Flat",
138   "kFMD3Flat",
139   "kFMDFlat"
140 };
141
142 //____________________________________________________________________
143 enum Geo_t {
144   kHoles,                       //
145   kNoHoles                      //
146 };
147
148 //____________________________________________________________________
149 enum Rad_t {
150   kGluonRadiation,              //
151   kNoGluonRadiation             //
152 };
153
154 //____________________________________________________________________
155 enum MC_t {
156   kFLUKA, 
157   kGEANT3, 
158   kGEANT4, 
159   kGEANT3TGEO,
160 };
161
162 //____________________________________________________________________
163 // Functions
164 Float_t       EtaToTheta(Float_t eta);
165 EG_t          LookupEG(const Char_t* name);
166 AliGenerator* GeneratorFactory(EG_t eg, Rad_t rad);
167 AliGenHijing* HijingStandard();
168 void          ProcessEnvironmentVars(EG_t& eg, Int_t& seed);
169
170 //____________________________________________________________________
171 void 
172 Config()
173 {
174   //____________________________________________________________________
175   // This part for configuration    
176   EG_t  eg   = kPythia6;
177   Geo_t geo  = kNoHoles;
178   Rad_t rad  = kGluonRadiation;
179   AliMagF::BMap_t mag  = AliMagF::k5kG;
180   Int_t seed = 12345; //Set 0 to use the current time
181   MC_t  mc   = kGEANT3TGEO;
182   
183   //____________________________________________________________________
184   // Get settings from environment variables
185   ProcessEnvironmentVars(eg, seed);
186
187   //____________________________________________________________________
188   // Set Random Number seed
189   gRandom->SetSeed(seed);
190   cout<<"Seed for random number generation= "<<gRandom->GetSeed()<<endl; 
191
192
193   //__________________________________________________________________
194   switch (mc) {
195   case kFLUKA: 
196     // 
197     // libraries required by fluka21
198     // 
199     gSystem->Load("libGeom");
200     cout << "\t* Loading TFluka..." << endl;  
201     gSystem->Load("libTFluka");    
202     gSystem->MakeDirectory("peg");
203     // 
204     // FLUKA MC
205     //
206     cout << "\t* Instantiating TFluka..." << endl;
207     new TFluka("C++ Interface to Fluka", 0/*verbosity*/);
208     break;
209   case kGEANT3: 
210     {
211       //
212       // Libraries needed by GEANT 3.21 
213       //
214       gSystem->Load("$ALICE_ROOT/lib/tgt_$ALICE_TARGET/libpythia6.so");
215       gSystem->Load("libgeant321");
216       
217       // 
218       // GEANT 3.21 MC 
219       // 
220       TGeant3* gmc = new TGeant3("C++ Interface to Geant3");
221       gmc->SetSWIT(4, 1000);
222     }
223     break;
224   case kGEANT3TGEO:
225     {
226       //
227       // Libraries needed by GEANT 3.21 
228       //
229       gSystem->Load("$ALICE_ROOT/lib/tgt_$ALICE_TARGET/liblhapdf.so");
230       gSystem->Load("$ALICE_ROOT/lib/tgt_$ALICE_TARGET/libpythia6.so");
231       gSystem->Load("libEGPythia6.so"); //<- For non-debian (sigh!)
232       // gSystem->Load("EGPythia6.so");
233       gSystem->Load("libgeant321");
234     
235       // 
236       // GEANT 3.21 MC 
237       // 
238       TGeant3TGeo* gmc  = new TGeant3TGeo("C++ Interface to Geant3");
239       gmc->SetSWIT(4, 1000);
240       Printf("Making a TGeant3TGeo objet");
241     }
242     break;
243   default:
244     gAlice->Fatal("Config.C", "No MC type chosen");
245     return;
246   }
247
248   //__________________________________________________________________
249   AliRunLoader* rl = 0;
250
251   cout<<"Config.C: Creating Run Loader ..."<<endl;
252   rl = AliRunLoader::Open("galice.root",
253                           AliConfig::GetDefaultEventFolderName(),
254                           "recreate");
255   if (!rl) {
256     gAlice->Fatal("Config.C","Can not instatiate the Run Loader");
257     return;
258   }
259   rl->SetCompressionLevel(2);
260   rl->SetNumberOfEventsPerFile(3);
261   gAlice->SetRunLoader(rl);
262
263   //__________________________________________________________________
264   // For FLUKA 
265   switch (mc) {
266   case kFLUKA: 
267     {
268       //
269       // Use kTRUE as argument to generate alice.pemf first
270       //
271       TString alice_pemf(gSystem->Which(".", "peg/mat17.pemf"));
272       if (!alice_pemf.IsNull()) 
273         ((TFluka*)gMC)->SetGeneratePemf(kFALSE);
274       else
275         ((TFluka*)gMC)->SetGeneratePemf(kTRUE);
276       TString flupro(gSystem->Getenv("FLUPRO"));
277       if (flupro.IsNull()) 
278         Fatal("Config.C", "Environment variable FLUPRO not set");
279 #if 0
280       char* files[] = { "brems_fin.bin", "cohff.bin", "elasct.bin", 
281                         "gxsect.bin", "nuclear.bin", "sigmapi.bin", 
282                         0 };
283       char* file = files[0];
284       while (file) {
285         TString which(gSystem->Which(".", file));
286         if (which.IsNull()) {
287           if (gSystem->Symlink(Form("%s/%s", flupro.Data(), file), file)!=0) 
288             Fatal("Config.C", "Couldn't link $(FLUPRO)/%s -> .", file);
289         }
290         file++;
291       }
292 #endif
293       TString neuxsc(gSystem->Which(".", "neuxsc.bin"));
294       if (neuxsc.IsNull()) 
295         gSystem->Symlink(Form("%s/neuxsc_72.bin", flupro.Data()), 
296                          "neuxsc.bin"); 
297       gSystem->CopyFile("$(FLUPRO)/random.dat", "old.seed", kTRUE);
298     }
299     break;
300   }
301
302   //__________________________________________________________________
303   //
304   // Set External decayer
305 #if 0
306   AliDecayer *decayer = new AliDecayerPythia();
307   switch (eg) {
308   case kD0PbPb5500:           decayer->SetForceDecay(kHadronicD);      break;
309   case kCharmSemiElPbPb5500:  decayer->SetForceDecay(kSemiElectronic); break;
310   case kBeautySemiElPbPb5500: decayer->SetForceDecay(kSemiElectronic); break;
311   default:                    decayer->SetForceDecay(kAll);            break;
312   }
313   decayer->Init();
314   gMC->SetExternalDecayer(decayer);
315 #endif
316
317   //__________________________________________________________________
318   // *********** STEERING parameters FOR ALICE SIMULATION ************
319   // - Specify event type to be tracked through the ALICE setup
320   // - All positions are in cm, angles in degrees, and P and E in GeV 
321   gMC->SetProcess("DCAY",1);
322   gMC->SetProcess("PAIR",1);
323   gMC->SetProcess("COMP",1);
324   gMC->SetProcess("PHOT",1);
325   gMC->SetProcess("PFIS",0);
326   gMC->SetProcess("DRAY",0);
327   gMC->SetProcess("ANNI",1);
328   gMC->SetProcess("BREM",1);
329   gMC->SetProcess("MUNU",1);
330   gMC->SetProcess("CKOV",1);
331   gMC->SetProcess("HADR",1);
332   gMC->SetProcess("LOSS",1); // 0:none 1,3:dray 2:nodray 4:nofluct (def:2)
333   gMC->SetProcess("MULS",1);
334   gMC->SetProcess("RAYL",1);
335
336   Float_t cut = 1.e-3;        // 1MeV cut by default
337   Float_t tofmax = 1.e10;
338
339   gMC->SetCut("CUTGAM", cut);
340   gMC->SetCut("CUTELE", cut);
341   gMC->SetCut("CUTNEU", cut);
342   gMC->SetCut("CUTHAD", cut);
343   gMC->SetCut("CUTMUO", cut);
344   gMC->SetCut("BCUTE",  cut); 
345   gMC->SetCut("BCUTM",  cut); 
346   gMC->SetCut("DCUTE",  cut); 
347   gMC->SetCut("DCUTM",  cut); 
348   gMC->SetCut("PPCUTM", cut);
349   gMC->SetCut("TOFMAX", tofmax); 
350
351   
352   //__________________________________________________________________
353   // Generator Configuration
354   AliGenerator* gener = GeneratorFactory(eg, rad);
355   gener->SetOrigin(0, 0, 0);    // vertex position
356   gener->SetSigma(0, 0, 5.3);   // Sigma in (X,Y,Z) (cm) on IP position
357   gener->SetCutVertexZ(1.);     // Truncate at 1 sigma
358   gener->SetVertexSmear(kPerEvent); 
359   gener->SetTrackingFlag(1);
360   gener->Init();
361     
362   //__________________________________________________________________
363   // Field (L3 0.4 T)
364   AliMagF* field = new AliMagF("Maps","Maps",2,1.,1., 10.,AliMagF::k5kG);
365   // field->SetL3ConstField(0); //Using const. field in the barrel
366   TGeoGlobalMagField::Instance()->SetField(field);
367
368   rl->CdGAFile();
369
370   TFile* magF = TFile::Open("mag.root", "RECREATE");
371   field->Write("mag");
372   magF->Close();
373
374   //__________________________________________________________________
375   // 
376   // Used detectors 
377   // 
378   Bool_t useABSO  = kFALSE; 
379   Bool_t useACORDE= kFALSE; 
380   Bool_t useDIPO  = kFALSE; 
381   Bool_t useFMD   = kTRUE; 
382   Bool_t useFRAME = kFALSE; 
383   Bool_t useHALL  = kFALSE; 
384   Bool_t useITS   = kTRUE;
385   Bool_t useMAG   = kFALSE; 
386   Bool_t useMUON  = kFALSE; 
387   Bool_t usePHOS  = kFALSE; 
388   Bool_t usePIPE  = kTRUE; 
389   Bool_t usePMD   = kFALSE; 
390   Bool_t useHMPID = kFALSE; 
391   Bool_t useSHIL  = kFALSE; 
392   Bool_t useT0    = kTRUE; 
393   Bool_t useTOF   = kFALSE; 
394   Bool_t useTPC   = kFALSE;
395   Bool_t useTRD   = kFALSE; 
396   Bool_t useZDC   = kFALSE; 
397   Bool_t useEMCAL = kFALSE; 
398   Bool_t useVZERO = kTRUE;
399
400   cout << "\t* Creating the detectors ..." << endl;
401   // ================= Alice BODY parameters =========================
402   AliBODY *BODY = new AliBODY("BODY", "Alice envelop");
403   
404   
405   if (useMAG)   AliMAG  *MAG  = new AliMAG("MAG", "Magnet");
406   if (useABSO)  AliABSO *ABSO = new AliABSOv3("ABSO", "Muon Absorber");
407   if (useDIPO)  AliDIPO *DIPO = new AliDIPOv3("DIPO", "Dipole version 3");
408   if (useHALL)  AliHALL *HALL = new AliHALLv3("HALL", "Alice Hall");
409   if (useFRAME) {
410     AliFRAMEv2 *FRAME = new AliFRAMEv2("FRAME", "Space Frame");
411     switch (geo) {
412     case kHoles: FRAME->SetHoles(1); break;
413     default:     FRAME->SetHoles(0); break;
414     }
415   }
416   if (useSHIL)   AliSHIL   *SHIL   = new AliSHILv3("SHIL", "Shielding v3");
417   if (usePIPE)   AliPIPE   *PIPE   = new AliPIPEv3("PIPE", "Beam Pipe");
418   if (useITS)    AliITS    *ITS    = new AliITSv11Hybrid("ITS","ITS v11Hybrid");
419   if (useTPC)    AliTPC    *TPC    = new AliTPCv2("TPC", "Default");
420   if (useTOF)    AliTOF    *TOF    = new AliTOFv6T0("TOF", "normal TOF");
421   if (useHMPID)  AliHMPID  *HMPID  = new AliHMPIDv1("HMPID", "normal HMPID");
422   if (useZDC)    AliZDC    *ZDC    = new AliZDCv3("ZDC", "normal ZDC");
423   if (useTRD)    AliTRD    *TRD    = new AliTRDv1("TRD", "TRD slow simulator");
424   if (useFMD)    AliFMD    *FMD    = new AliFMDv1("FMD", "normal FMD");
425   if (useMUON)   AliMUON   *MUON   = new AliMUONv1("MUON", "default");
426   if (usePHOS)   AliPHOS   *PHOS   = new AliPHOSv1("PHOS", "IHEP");
427   if (usePMD)    AliPMD    *PMD    = new AliPMDv1("PMD", "normal PMD");
428   if (useT0)     AliT0     *T0     = new AliT0v1("T0", "T0 Detector");
429   if (useEMCAL)  AliEMCAL  *EMCAL  = new AliEMCALv2("EMCAL", "EMCAL_COMPLETE");
430   if (useACORDE) AliACORDE *ACORDE = new AliACORDEv1("ACORDE", "normal ACORDE");
431   if (useVZERO)  AliVZERO  *VZERO  = new AliVZEROv7("VZERO", "normal VZERO");
432 }
433
434 //____________________________________________________________________
435 Float_t EtaToTheta(Float_t arg)
436 {
437   return (180./TMath::Pi())*2.*TMath::ATan(TMath::Exp(-arg));
438 }
439
440 //____________________________________________________________________
441 Int_t 
442 LookupEG(const Char_t* name) 
443 {
444   TString n(name);
445   for (Int_t i = 0; i < kEgMax; i++) {
446     if (n == egName[i]) return i;
447   }
448   return -1;
449 }
450
451 //____________________________________________________________________  
452 AliGenerator* 
453 GeneratorFactory(EG_t eg, Rad_t rad)  
454 {
455   Int_t isw = 3;
456   if (rad == kNoGluonRadiation) isw = 0;
457   
458   // Possibly load libAliPythia6
459   switch (eg) {
460   case kPythia6:
461   case kPythia6Jets20_24:
462   case kPythia6Jets24_29:
463   case kPythia6Jets29_35:
464   case kPythia6Jets35_42:
465   case kPythia6Jets42_50:
466   case kPythia6Jets50_60:
467   case kPythia6Jets60_72:
468   case kPythia6Jets72_86:
469   case kPythia6Jets86_104:
470   case kPythia6Jets104_125:
471   case kPythia6Jets125_150:
472   case kPythia6Jets150_180:
473     gSystem->Load("liblhapdf.so");
474     // gSystem->Load("/usr/lib/libpythia.so");
475     // gSystem->ListLibraries();
476     gSystem->Load("EGPythia6.so");
477     gSystem->Load("libAliPythia6");
478     break;
479   }
480   
481   AliGenerator * gGener = 0;
482   switch (eg) {
483   case test50:
484     {
485       AliGenHIJINGpara *gener = new AliGenHIJINGpara(50);
486       gener->SetMomentumRange(0, 999999.);
487       gener->SetPhiRange(0., 360.);
488       // Set pseudorapidity range from -8 to 8.
489       Float_t thmin = EtaToTheta(8);   // theta min. <---> eta max
490       Float_t thmax = EtaToTheta(-8);  // theta max. <---> eta min 
491       gener->SetThetaRange(thmin,thmax);
492       gGener=gener;
493     }
494     break;
495   case kParam_8000:
496     {
497       AliGenHIJINGpara *gener = new AliGenHIJINGpara(86030);
498       gener->SetMomentumRange(0, 999999.);
499       gener->SetPhiRange(0., 360.);
500       // Set pseudorapidity range from -8 to 8.
501       Float_t thmin = EtaToTheta(8);   // theta min. <---> eta max
502       Float_t thmax = EtaToTheta(-8);  // theta max. <---> eta min 
503       gener->SetThetaRange(thmin,thmax);
504       gGener=gener;
505     }
506     break;
507   case kParam_4000:
508     {
509       AliGenHIJINGpara *gener = new AliGenHIJINGpara(43015);
510       gener->SetMomentumRange(0, 999999.);
511       gener->SetPhiRange(0., 360.);
512       // Set pseudorapidity range from -8 to 8.
513       Float_t thmin = EtaToTheta(8);   // theta min. <---> eta max
514       Float_t thmax = EtaToTheta(-8);  // theta max. <---> eta min 
515       gener->SetThetaRange(thmin,thmax);
516       gGener=gener;
517     }
518     break;
519   case kParam_2000:
520     {
521       AliGenHIJINGpara *gener = new AliGenHIJINGpara(21507);
522       gener->SetMomentumRange(0, 999999.);
523       gener->SetPhiRange(0., 360.);
524       // Set pseudorapidity range from -8 to 8.
525       Float_t thmin = EtaToTheta(8);   // theta min. <---> eta max
526       Float_t thmax = EtaToTheta(-8);  // theta max. <---> eta min 
527       gener->SetThetaRange(thmin,thmax);
528       gGener=gener;
529     }
530     break;
531   case kParam_fmd:
532     {
533       AliGenHIJINGpara *gener = new AliGenHIJINGpara(500);
534       gener->SetMomentumRange(0, 999999.);
535       gener->SetPhiRange(0., 360.);
536       // Set pseudorapidity range from -8 to 8.
537       Float_t thmin = EtaToTheta(6);   // theta min. <---> eta max
538       Float_t thmax = EtaToTheta(2);  // theta max. <---> eta min 
539       gener->SetThetaRange(thmin,thmax);
540       gGener=gener;
541     }
542     break;
543     //
544     //  Hijing Central
545     //
546   case kHijing_cent1:
547     {
548       AliGenHijing *gener = HijingStandard();
549       // impact parameter range
550       gener->SetImpactParameterRange(0., 5.);
551       gGener=gener;
552     }
553     break;
554   case kHijing_cent2:
555     {
556       AliGenHijing *gener = HijingStandard();
557       // impact parameter range
558       gener->SetImpactParameterRange(0., 2.);
559       gGener=gener;
560     }
561     break;
562     //
563     // Hijing Peripheral 
564     //
565   case kHijing_per1:
566     {
567       AliGenHijing *gener = HijingStandard();
568       // impact parameter range
569       gener->SetImpactParameterRange(5., 8.6);
570       gGener=gener;
571     }
572     break;
573   case kHijing_per2:
574     {
575       AliGenHijing *gener = HijingStandard();
576       // impact parameter range
577       gener->SetImpactParameterRange(8.6, 11.2);
578       gGener=gener;
579     }
580     break;
581   case kHijing_per3:
582     {
583       AliGenHijing *gener = HijingStandard();
584       // impact parameter range
585       gener->SetImpactParameterRange(11.2, 13.2);
586       gGener=gener;
587     }
588     break;
589   case kHijing_per4:
590     {
591       AliGenHijing *gener = HijingStandard();
592       // impact parameter range
593       gener->SetImpactParameterRange(13.2, 15.);
594       gGener=gener;
595     }
596     break;
597   case kHijing_per5:
598     {
599       AliGenHijing *gener = HijingStandard();
600       // impact parameter range
601       gener->SetImpactParameterRange(15., 100.);
602       gGener=gener;
603     }
604     break;
605     //
606     //  Jet-Jet
607     //
608   case kHijing_jj25:
609     {
610       AliGenHijing *gener = HijingStandard();
611       // impact parameter range
612       gener->SetImpactParameterRange(0., 5.);
613       // trigger
614       gener->SetTrigger(1);
615       gener->SetPtJet(25.);
616       gener->SetRadiation(isw);
617       gener->SetSimpleJets(!isw);
618       gener->SetJetEtaRange(-0.3,0.3);
619       gener->SetJetPhiRange(75., 165.);   
620       gGener=gener;
621     }
622     break;
623
624   case kHijing_jj50:
625     {
626       AliGenHijing *gener = HijingStandard();
627       // impact parameter range
628       gener->SetImpactParameterRange(0., 5.);
629       // trigger
630       gener->SetTrigger(1);
631       gener->SetPtJet(50.);
632       gener->SetRadiation(isw);
633       gener->SetSimpleJets(!isw);
634       gener->SetJetEtaRange(-0.3,0.3);
635       gener->SetJetPhiRange(75., 165.);   
636       gGener=gener;
637     }
638     break;
639
640   case kHijing_jj75:
641     {
642       AliGenHijing *gener = HijingStandard();
643       // impact parameter range
644       gener->SetImpactParameterRange(0., 5.);
645       // trigger
646       gener->SetTrigger(1);
647       gener->SetPtJet(75.);
648       gener->SetRadiation(isw);
649       gener->SetSimpleJets(!isw);
650       gener->SetJetEtaRange(-0.3,0.3);
651       gener->SetJetPhiRange(75., 165.);   
652       gGener=gener;
653     }
654     break;
655
656   case kHijing_jj100:
657     {
658       AliGenHijing *gener = HijingStandard();
659       // impact parameter range
660       gener->SetImpactParameterRange(0., 5.);
661       // trigger
662       gener->SetTrigger(1);
663       gener->SetPtJet(100.);
664       gener->SetRadiation(isw);
665       gener->SetSimpleJets(!isw);
666       gener->SetJetEtaRange(-0.3,0.3);
667       gener->SetJetPhiRange(75., 165.);   
668       gGener=gener;
669     }
670     break;
671
672   case kHijing_jj200:
673     {
674       AliGenHijing *gener = HijingStandard();
675       // impact parameter range
676       gener->SetImpactParameterRange(0., 5.);
677       // trigger
678       gener->SetTrigger(1);
679       gener->SetPtJet(200.);
680       gener->SetRadiation(isw);
681       gener->SetSimpleJets(!isw);
682       gener->SetJetEtaRange(-0.3,0.3);
683       gener->SetJetPhiRange(75., 165.);   
684       gGener=gener;
685     }
686     break;
687     //
688     // Gamma-Jet
689     //
690   case kHijing_gj25:
691     {
692       AliGenHijing *gener = HijingStandard();
693       // impact parameter range
694       gener->SetImpactParameterRange(0., 5.);
695       // trigger
696       gener->SetTrigger(2);
697       gener->SetPtJet(25.);
698       gener->SetRadiation(isw);
699       gener->SetSimpleJets(!isw);
700       gener->SetJetEtaRange(-0.12, 0.12);
701       gener->SetJetPhiRange(220., 320.);
702       gGener=gener;
703     }
704     break;
705
706   case kHijing_gj50:
707     {
708       AliGenHijing *gener = HijingStandard();
709       // impact parameter range
710       gener->SetImpactParameterRange(0., 5.);
711       // trigger
712       gener->SetTrigger(2);
713       gener->SetPtJet(50.);
714       gener->SetRadiation(isw);
715       gener->SetSimpleJets(!isw);
716       gener->SetJetEtaRange(-0.12, 0.12);
717       gener->SetJetPhiRange(220., 320.);
718       gGener=gener;
719     }
720     break;
721
722   case kHijing_gj75:
723     {
724       AliGenHijing *gener = HijingStandard();
725       // impact parameter range
726       gener->SetImpactParameterRange(0., 5.);
727       // trigger
728       gener->SetTrigger(2);
729       gener->SetPtJet(75.);
730       gener->SetRadiation(isw);
731       gener->SetSimpleJets(!isw);
732       gener->SetJetEtaRange(-0.12, 0.12);
733       gener->SetJetPhiRange(220., 320.);
734       gGener=gener;
735     }
736     break;
737
738   case kHijing_gj100:
739     {
740       AliGenHijing *gener = HijingStandard();
741       // impact parameter range
742       gener->SetImpactParameterRange(0., 5.);
743       // trigger
744       gener->SetTrigger(2);
745       gener->SetPtJet(100.);
746       gener->SetRadiation(isw);
747       gener->SetSimpleJets(!isw);
748       gener->SetJetEtaRange(-0.12, 0.12);
749       gener->SetJetPhiRange(220., 320.);
750       gGener=gener;
751     }
752     break;
753
754   case kHijing_gj200:
755     {
756       AliGenHijing *gener = HijingStandard();
757       // impact parameter range
758       gener->SetImpactParameterRange(0., 5.);
759       // trigger
760       gener->SetTrigger(2);
761       gener->SetPtJet(200.);
762       gener->SetRadiation(isw);
763       gener->SetSimpleJets(!isw);
764       gener->SetJetEtaRange(-0.12, 0.12);
765       gener->SetJetPhiRange(220., 320.);
766       gGener=gener;
767     }
768     break;
769   case kHijing_pA:
770     {
771
772       AliGenCocktail *gener  = new AliGenCocktail();
773
774       AliGenHijing   *hijing = new AliGenHijing(-1);
775       // centre of mass energy 
776       hijing->SetEnergyCMS(TMath::Sqrt(82./208.) * 14000.);
777       // impact parameter range
778       hijing->SetImpactParameterRange(0., 15.);
779       // reference frame
780       hijing->SetReferenceFrame("CMS");
781       hijing->SetBoostLHC(1);
782       // projectile
783       hijing->SetProjectile("P", 1, 1);
784       hijing->SetTarget    ("A", 208, 82);
785       // tell hijing to keep the full parent child chain
786       hijing->KeepFullEvent();
787       // enable jet quenching
788       hijing->SetJetQuenching(0);
789       // enable shadowing
790       hijing->SetShadowing(1);
791       // Don't track spectators
792       hijing->SetSpectators(0);
793       // kinematic selection
794       hijing->SetSelectAll(0);
795       //
796       AliGenSlowNucleons*  gray    = new AliGenSlowNucleons(1);
797       AliSlowNucleonModel* model   = new AliSlowNucleonModelExp();
798       gray->SetSlowNucleonModel(model);
799       gray->SetDebug(1);
800       gener->AddGenerator(hijing,"Hijing pPb", 1);
801       gener->AddGenerator(gray,  "Gray Particles",1);
802       gGener=gener;
803     }
804     break;
805   case kPythia6:
806     {
807       AliGenPythia *gener = new AliGenPythia(-1); 
808       gener->SetMomentumRange(0,999999);
809       gener->SetThetaRange(0., 180.);
810       gener->SetYRange(-12,12);
811       gener->SetPtRange(0,1000);
812       gener->SetProcess(kPyMb);
813       gener->SetEnergyCMS(14000.);
814       gGener=gener;
815     }
816     break;
817   case kPythia6Jets20_24:
818     {
819       AliGenPythia * gener = new AliGenPythia(-1);
820       gener->SetEnergyCMS(5500.);//        Centre of mass energy
821       gener->SetProcess(kPyJets);//        Process type
822       gener->SetJetEtaRange(-0.5, 0.5);//  Final state kinematic cuts
823       gener->SetJetPhiRange(0., 360.);
824       gener->SetJetEtRange(10., 1000.);
825       gener->SetGluonRadiation(1,1);
826       //    gener->SetPtKick(0.);
827       //   Structure function
828       gener->SetStrucFunc(kCTEQ4L);
829       gener->SetPtHard(20., 24.);// Pt transfer of the hard scattering
830       gener->SetPycellParameters(2., 274, 432, 0., 4., 5., 1.0);
831       gener->SetForceDecay(kAll);//  Decay type (semielectronic, etc.)
832       gGener=gener;
833     }
834     break;
835   case kPythia6Jets24_29:
836     {
837       AliGenPythia * gener = new AliGenPythia(-1);
838       gener->SetEnergyCMS(5500.);//        Centre of mass energy
839       gener->SetProcess(kPyJets);//        Process type
840       gener->SetJetEtaRange(-0.5, 0.5);//  Final state kinematic cuts
841       gener->SetJetPhiRange(0., 360.);
842       gener->SetJetEtRange(10., 1000.);
843       gener->SetGluonRadiation(1,1);
844       //    gener->SetPtKick(0.);
845       //   Structure function
846       gener->SetStrucFunc(kCTEQ4L);
847       gener->SetPtHard(24., 29.);// Pt transfer of the hard scattering
848       gener->SetPycellParameters(2., 274, 432, 0., 4., 5., 1.0);
849       gener->SetForceDecay(kAll);//  Decay type (semielectronic, etc.)
850       gGener=gener;
851     }
852     break;
853   case kPythia6Jets29_35:
854     {
855       AliGenPythia * gener = new AliGenPythia(-1);
856       gener->SetEnergyCMS(5500.);//        Centre of mass energy
857       gener->SetProcess(kPyJets);//        Process type
858       gener->SetJetEtaRange(-0.5, 0.5);//  Final state kinematic cuts
859       gener->SetJetPhiRange(0., 360.);
860       gener->SetJetEtRange(10., 1000.);
861       gener->SetGluonRadiation(1,1);
862       //    gener->SetPtKick(0.);
863       //   Structure function
864       gener->SetStrucFunc(kCTEQ4L);
865       gener->SetPtHard(29., 35.);// Pt transfer of the hard scattering
866       gener->SetPycellParameters(2., 274, 432, 0., 4., 5., 1.0);
867       gener->SetForceDecay(kAll);//  Decay type (semielectronic, etc.)
868       gGener=gener;
869     }
870     break;
871   case kPythia6Jets35_42:
872     {
873       AliGenPythia * gener = new AliGenPythia(-1);
874       gener->SetEnergyCMS(5500.);//        Centre of mass energy
875       gener->SetProcess(kPyJets);//        Process type
876       gener->SetJetEtaRange(-0.5, 0.5);//  Final state kinematic cuts
877       gener->SetJetPhiRange(0., 360.);
878       gener->SetJetEtRange(10., 1000.);
879       gener->SetGluonRadiation(1,1);
880       //    gener->SetPtKick(0.);
881       //   Structure function
882       gener->SetStrucFunc(kCTEQ4L);
883       gener->SetPtHard(35., 42.);// Pt transfer of the hard scattering
884       gener->SetPycellParameters(2., 274, 432, 0., 4., 5., 1.0);
885       gener->SetForceDecay(kAll);//  Decay type (semielectronic, etc.)
886       gGener=gener;
887     }
888     break;
889   case kPythia6Jets42_50:
890     {
891       AliGenPythia * gener = new AliGenPythia(-1);
892       gener->SetEnergyCMS(5500.);//        Centre of mass energy
893       gener->SetProcess(kPyJets);//        Process type
894       gener->SetJetEtaRange(-0.5, 0.5);//  Final state kinematic cuts
895       gener->SetJetPhiRange(0., 360.);
896       gener->SetJetEtRange(10., 1000.);
897       gener->SetGluonRadiation(1,1);
898       //    gener->SetPtKick(0.);
899       //   Structure function
900       gener->SetStrucFunc(kCTEQ4L);
901       gener->SetPtHard(42., 50.);// Pt transfer of the hard scattering
902       gener->SetPycellParameters(2., 274, 432, 0., 4., 5., 1.0);
903       gener->SetForceDecay(kAll);//  Decay type (semielectronic, etc.)
904       gGener=gener;
905     }
906     break;
907   case kPythia6Jets50_60:
908     {
909       AliGenPythia * gener = new AliGenPythia(-1);
910       gener->SetEnergyCMS(5500.);//        Centre of mass energy
911       gener->SetProcess(kPyJets);//        Process type
912       gener->SetJetEtaRange(-0.5, 0.5);//  Final state kinematic cuts
913       gener->SetJetPhiRange(0., 360.);
914       gener->SetJetEtRange(10., 1000.);
915       gener->SetGluonRadiation(1,1);
916       //    gener->SetPtKick(0.);
917       //   Structure function
918       gener->SetStrucFunc(kCTEQ4L);
919       gener->SetPtHard(50., 60.);// Pt transfer of the hard scattering
920       gener->SetPycellParameters(2., 274, 432, 0., 4., 5., 1.0);
921       gener->SetForceDecay(kAll);//  Decay type (semielectronic, etc.)
922       gGener=gener;
923     }
924     break;
925   case kPythia6Jets60_72:
926     {
927       AliGenPythia * gener = new AliGenPythia(-1);
928       gener->SetEnergyCMS(5500.);//        Centre of mass energy
929       gener->SetProcess(kPyJets);//        Process type
930       gener->SetJetEtaRange(-0.5, 0.5);//  Final state kinematic cuts
931       gener->SetJetPhiRange(0., 360.);
932       gener->SetJetEtRange(10., 1000.);
933       gener->SetGluonRadiation(1,1);
934       //    gener->SetPtKick(0.);
935       //   Structure function
936       gener->SetStrucFunc(kCTEQ4L);
937       gener->SetPtHard(60., 72.);// Pt transfer of the hard scattering
938       gener->SetPycellParameters(2., 274, 432, 0., 4., 5., 1.0);
939       gener->SetForceDecay(kAll);//  Decay type (semielectronic, etc.)
940       gGener=gener;
941     }
942     break;
943   case kPythia6Jets72_86:
944     {
945       AliGenPythia * gener = new AliGenPythia(-1);
946       gener->SetEnergyCMS(5500.);//        Centre of mass energy
947       gener->SetProcess(kPyJets);//        Process type
948       gener->SetJetEtaRange(-0.5, 0.5);//  Final state kinematic cuts
949       gener->SetJetPhiRange(0., 360.);
950       gener->SetJetEtRange(10., 1000.);
951       gener->SetGluonRadiation(1,1);
952       //    gener->SetPtKick(0.);
953       //   Structure function
954       gener->SetStrucFunc(kCTEQ4L);
955       gener->SetPtHard(72., 86.);// Pt transfer of the hard scattering
956       gener->SetPycellParameters(2., 274, 432, 0., 4., 5., 1.0);
957       gener->SetForceDecay(kAll);//  Decay type (semielectronic, etc.)
958       gGener=gener;
959     }
960     break;
961   case kPythia6Jets86_104:
962     {
963       AliGenPythia * gener = new AliGenPythia(-1);
964       gener->SetEnergyCMS(5500.);//        Centre of mass energy
965       gener->SetProcess(kPyJets);//        Process type
966       gener->SetJetEtaRange(-0.5, 0.5);//  Final state kinematic cuts
967       gener->SetJetPhiRange(0., 360.);
968       gener->SetJetEtRange(10., 1000.);
969       gener->SetGluonRadiation(1,1);
970       //    gener->SetPtKick(0.);
971       //   Structure function
972       gener->SetStrucFunc(kCTEQ4L);
973       gener->SetPtHard(86., 104.);// Pt transfer of the hard scattering
974       gener->SetPycellParameters(2., 274, 432, 0., 4., 5., 1.0);
975       gener->SetForceDecay(kAll);//  Decay type (semielectronic, etc.)
976       gGener=gener;
977     }
978     break;
979   case kPythia6Jets104_125:
980     {
981       AliGenPythia * gener = new AliGenPythia(-1);
982       gener->SetEnergyCMS(5500.);//        Centre of mass energy
983       gener->SetProcess(kPyJets);//        Process type
984       gener->SetJetEtaRange(-0.5, 0.5);//  Final state kinematic cuts
985       gener->SetJetPhiRange(0., 360.);
986       gener->SetJetEtRange(10., 1000.);
987       gener->SetGluonRadiation(1,1);
988       //    gener->SetPtKick(0.);
989       //   Structure function
990       gener->SetStrucFunc(kCTEQ4L);
991       gener->SetPtHard(104., 125.);// Pt transfer of the hard scattering
992       gener->SetPycellParameters(2., 274, 432, 0., 4., 5., 1.0);
993       gener->SetForceDecay(kAll);//  Decay type (semielectronic, etc.)
994       gGener=gener;
995     }
996     break;
997   case kPythia6Jets125_150:
998     {
999       AliGenPythia * gener = new AliGenPythia(-1);
1000       gener->SetEnergyCMS(5500.);//        Centre of mass energy
1001       gener->SetProcess(kPyJets);//        Process type
1002       gener->SetJetEtaRange(-0.5, 0.5);//  Final state kinematic cuts
1003       gener->SetJetPhiRange(0., 360.);
1004       gener->SetJetEtRange(10., 1000.);
1005       gener->SetGluonRadiation(1,1);
1006       //    gener->SetPtKick(0.);
1007       //   Structure function
1008       gener->SetStrucFunc(kCTEQ4L);
1009       gener->SetPtHard(125., 150.);// Pt transfer of the hard scattering
1010       gener->SetPycellParameters(2., 274, 432, 0., 4., 5., 1.0);
1011       gener->SetForceDecay(kAll);//  Decay type (semielectronic, etc.)
1012       gGener=gener;
1013     }
1014     break;
1015   case kPythia6Jets150_180:
1016     {
1017       AliGenPythia * gener = new AliGenPythia(-1);
1018       gener->SetEnergyCMS(5500.);//        Centre of mass energy
1019       gener->SetProcess(kPyJets);//        Process type
1020       gener->SetJetEtaRange(-0.5, 0.5);//  Final state kinematic cuts
1021       gener->SetJetPhiRange(0., 360.);
1022       gener->SetJetEtRange(10., 1000.);
1023       gener->SetGluonRadiation(1,1);
1024       //    gener->SetPtKick(0.);
1025       //   Structure function
1026       gener->SetStrucFunc(kCTEQ4L);
1027       gener->SetPtHard(150., 180.);// Pt transfer of the hard scattering
1028       gener->SetPycellParameters(2., 274, 432, 0., 4., 5., 1.0);
1029       gener->SetForceDecay(kAll);//  Decay type (semielectronic, etc.)
1030       gGener=gener;
1031     }
1032     break;
1033   case kD0PbPb5500:
1034     {
1035       AliGenPythia * gener = new AliGenPythia(10);
1036       gener->SetProcess(kPyD0PbPbMNR);
1037       gener->SetStrucFunc(kCTEQ4L);
1038       gener->SetPtHard(2.1,-1.0);
1039       gener->SetEnergyCMS(5500.);
1040       gener->SetNuclei(208,208);
1041       gener->SetForceDecay(kHadronicD);
1042       gener->SetYRange(-2,2);
1043       gener->SetFeedDownHigherFamily(kFALSE);
1044       gener->SetStackFillOpt(AliGenPythia::kParentSelection);
1045       gener->SetCountMode(AliGenPythia::kCountParents);
1046       gGener=gener;
1047     }
1048     break;
1049   case kCharmSemiElPbPb5500:
1050     {
1051       AliGenPythia * gener = new AliGenPythia(10);
1052       gener->SetProcess(kPyCharmPbPbMNR);
1053       gener->SetStrucFunc(kCTEQ4L);
1054       gener->SetPtHard(2.1,-1.0);
1055       gener->SetEnergyCMS(5500.);
1056       gener->SetNuclei(208,208);
1057       gener->SetForceDecay(kSemiElectronic);
1058       gener->SetYRange(-2,2);
1059       gener->SetFeedDownHigherFamily(kFALSE);
1060       gener->SetCountMode(AliGenPythia::kCountParents);
1061       gGener=gener;
1062     }
1063     break;
1064   case kBeautySemiElPbPb5500:
1065     {
1066       AliGenPythia *gener = new AliGenPythia(10);
1067       gener->SetProcess(kPyBeautyPbPbMNR);
1068       gener->SetStrucFunc(kCTEQ4L);
1069       gener->SetPtHard(2.75,-1.0);
1070       gener->SetEnergyCMS(5500.);
1071       gener->SetNuclei(208,208);
1072       gener->SetForceDecay(kSemiElectronic);
1073       gener->SetYRange(-2,2);
1074       gener->SetFeedDownHigherFamily(kFALSE);
1075       gener->SetCountMode(AliGenPythia::kCountParents);
1076       gGener=gener;
1077     }
1078     break;
1079   case kCocktailTRD:
1080     {
1081       AliGenCocktail *gener  = new AliGenCocktail();
1082
1083       AliGenParam *jpsi = new AliGenParam(10,
1084                                           new AliGenMUONlib(),
1085                                           AliGenMUONlib::kJpsiFamily,
1086                                           "Vogt PbPb");
1087
1088       jpsi->SetPtRange(0, 100);
1089       jpsi->SetYRange(-1., +1.);
1090       jpsi->SetForceDecay(kDiElectron);
1091
1092       AliGenParam *ups = new AliGenParam(10,
1093                                          new AliGenMUONlib(),
1094                                          AliGenMUONlib::kUpsilonFamily,
1095                                          "Vogt PbPb");
1096       ups->SetPtRange(0, 100);
1097       ups->SetYRange(-1., +1.);
1098       ups->SetForceDecay(kDiElectron);
1099         
1100       AliGenParam *charm = new AliGenParam(10,
1101                                            new AliGenMUONlib(), 
1102                                            AliGenMUONlib::kCharm,
1103                                            "central");
1104       charm->SetPtRange(0, 100);
1105       charm->SetYRange(-1.5, +1.5);
1106       charm->SetForceDecay(kSemiElectronic);
1107         
1108         
1109       AliGenParam *beauty = new AliGenParam(10,
1110                                             new AliGenMUONlib(), 
1111                                             AliGenMUONlib::kBeauty,
1112                                             "central");
1113       beauty->SetPtRange(0, 100);
1114       beauty->SetYRange(-1.5, +1.5);
1115       beauty->SetForceDecay(kSemiElectronic);
1116
1117       gener->AddGenerator(jpsi,"J/psi",1);
1118       gener->AddGenerator(ups,"Upsilon",1);
1119       gener->AddGenerator(charm,"Charm",1);
1120       gener->AddGenerator(beauty,"Beauty",1);
1121       gGener=gener;
1122     }
1123     break;
1124   case kPyJJ:
1125     {
1126       AliGenPythia *gener = new AliGenPythia(-1);
1127       gener->SetEnergyCMS(5500.);
1128       gener->SetProcess(kPyJets);
1129       Double_t ptHardMin=10.0, ptHardMax=-1.0;
1130       gener->SetPtHard(ptHardMin,ptHardMax);
1131       gener->SetYHard(-0.7,0.7);
1132       gener->SetJetEtaRange(-0.2,0.2);
1133       gener->SetEventListRange(0,1);
1134       gGener=gener;
1135     }
1136     break;
1137   case kPyGJ:
1138     {
1139       AliGenPythia *gener = new AliGenPythia(-1);
1140       gener->SetEnergyCMS(5500.);
1141       gener->SetProcess(kPyDirectGamma);
1142       Double_t ptHardMin=10.0, ptHardMax=-1.0;
1143       gener->SetPtHard(ptHardMin,ptHardMax);
1144       gener->SetYHard(-1.0,1.0);
1145       gener->SetGammaEtaRange(-0.13,0.13);
1146       gener->SetGammaPhiRange(210.,330.);
1147       gener->SetEventListRange(0,1);
1148       gGener=gener;
1149     }
1150     break;
1151   case kMuonCocktailCent1:
1152     {
1153       AliGenMUONCocktail * gener = new AliGenMUONCocktail();
1154       gener->SetPtRange(1.0,100.);       // Transverse momentum range   
1155       gener->SetPhiRange(0.,360.);    // Azimuthal angle range  
1156       gener->SetYRange(-4.0,-2.4);
1157       gener->SetMuonPtCut(0.8);
1158       gener->SetMuonThetaCut(171.,178.);
1159       gener->SetMuonMultiplicity(2);
1160       gener->SetNumberOfCollisions(1626.);  //Centrality class Cent1 for PDC04
1161       gener->SetNumberOfParticipants(359.4);//Centrality class Cent1 for PDC04
1162       gGener=gener;
1163     }
1164     break;
1165   case kMuonCocktailPer1:
1166     {
1167       AliGenMUONCocktail * gener = new AliGenMUONCocktail();
1168       gener->SetPtRange(1.0,100.);       // Transverse momentum range   
1169       gener->SetPhiRange(0.,360.);    // Azimuthal angle range  
1170       gener->SetYRange(-4.0,-2.4);
1171       gener->SetMuonPtCut(0.8);
1172       gener->SetMuonThetaCut(171.,178.);
1173       gener->SetMuonMultiplicity(2);
1174       gener->SetNumberOfCollisions(820.0);//Centrality class Per1 for PDC04
1175       gener->SetNumberOfParticipants(229.3);//Centrality class Per1 for PDC04
1176       gGener=gener;
1177     }
1178     break;
1179   case kMuonCocktailPer4:
1180     {
1181       AliGenMUONCocktail * gener = new AliGenMUONCocktail();
1182       gener->SetPtRange(1.0,100.);       // Transverse momentum range   
1183       gener->SetPhiRange(0.,360.);    // Azimuthal angle range  
1184       gener->SetYRange(-4.0,-2.4);
1185       gener->SetMuonPtCut(0.8);
1186       gener->SetMuonThetaCut(171.,178.);
1187       gener->SetMuonMultiplicity(2);
1188       gener->SetNumberOfCollisions(13.6);//Centrality class Per4 for PDC04
1189       gener->SetNumberOfParticipants(13.3);//Centrality class Per4 for PDC04
1190       gGener=gener;
1191     }
1192     break;
1193   case kMuonCocktailCent1HighPt:
1194     {
1195       AliGenMUONCocktail * gener = new AliGenMUONCocktail();
1196       gener->SetPtRange(1.0,100.);       // Transverse momentum range   
1197       gener->SetPhiRange(0.,360.);    // Azimuthal angle range  
1198       gener->SetYRange(-4.0,-2.4);
1199       gener->SetMuonPtCut(2.5);
1200       gener->SetMuonThetaCut(171.,178.);
1201       gener->SetMuonMultiplicity(2);
1202       gener->SetNumberOfCollisions(1626.);  //Centrality class Cent1 for PDC04
1203       gener->SetNumberOfParticipants(359.4);//Centrality class Cent1 for PDC04
1204       gGener=gener;
1205     }
1206     break;
1207   case kMuonCocktailPer1HighPt :
1208     {
1209       AliGenMUONCocktail * gener = new AliGenMUONCocktail();
1210       gener->SetPtRange(1.0,100.);       // Transverse momentum range   
1211       gener->SetPhiRange(0.,360.);    // Azimuthal angle range  
1212       gener->SetYRange(-4.0,-2.4);
1213       gener->SetMuonPtCut(2.5);
1214       gener->SetMuonThetaCut(171.,178.);
1215       gener->SetMuonMultiplicity(2);
1216       gener->SetNumberOfCollisions(820.0);//Centrality class Per1 for PDC04
1217       gener->SetNumberOfParticipants(229.3);//Centrality class Per1 for PDC04
1218       gGener=gener;
1219     }
1220     break;
1221   case kMuonCocktailPer4HighPt:
1222     {
1223       AliGenMUONCocktail * gener = new AliGenMUONCocktail();
1224       gener->SetPtRange(1.0,100.);       // Transverse momentum range   
1225       gener->SetPhiRange(0.,360.);    // Azimuthal angle range  
1226       gener->SetYRange(-4.0,-2.4);
1227       gener->SetMuonPtCut(2.5);
1228       gener->SetMuonThetaCut(171.,178.);
1229       gener->SetMuonMultiplicity(2);
1230       gener->SetNumberOfCollisions(13.6);//Centrality class Per4 for PDC04
1231       gener->SetNumberOfParticipants(13.3);//Centrality class Per4 for PDC04
1232       gGener=gener;
1233     }
1234     break;
1235   case kMuonCocktailCent1Single:
1236     {
1237       AliGenMUONCocktail * gener = new AliGenMUONCocktail();
1238       gener->SetPtRange(1.0,100.);       // Transverse momentum range   
1239       gener->SetPhiRange(0.,360.);    // Azimuthal angle range  
1240       gener->SetYRange(-4.0,-2.4);
1241       gener->SetMuonPtCut(0.8);
1242       gener->SetMuonThetaCut(171.,178.);
1243       gener->SetMuonMultiplicity(1);
1244       gener->SetNumberOfCollisions(1626.);  //Centrality class Cent1 for PDC04
1245       gener->SetNumberOfParticipants(359.4);//Centrality class Cent1 for PDC04
1246       gGener=gener;
1247     }
1248     break;
1249   case kMuonCocktailPer1Single :
1250     {
1251       AliGenMUONCocktail * gener = new AliGenMUONCocktail();
1252       gener->SetPtRange(1.0,100.);       // Transverse momentum range   
1253       gener->SetPhiRange(0.,360.);    // Azimuthal angle range  
1254       gener->SetYRange(-4.0,-2.4);
1255       gener->SetMuonPtCut(0.8);
1256       gener->SetMuonThetaCut(171.,178.);
1257       gener->SetMuonMultiplicity(1);
1258       gener->SetNumberOfCollisions(820.0);//Centrality class Per1 for PDC04
1259       gener->SetNumberOfParticipants(229.3);//Centrality class Per1 for PDC04
1260       gGener=gener;
1261     }
1262     break;
1263   case kMuonCocktailPer4Single:
1264     {
1265       AliGenMUONCocktail * gener = new AliGenMUONCocktail();
1266       gener->SetPtRange(1.0,100.);       // Transverse momentum range   
1267       gener->SetPhiRange(0.,360.);    // Azimuthal angle range  
1268       gener->SetYRange(-4.0,-2.4);
1269       gener->SetMuonPtCut(0.8);
1270       gener->SetMuonThetaCut(171.,178.);
1271       gener->SetMuonMultiplicity(1);
1272       gener->SetNumberOfCollisions(13.6);//Centrality class Per4 for PDC04
1273       gener->SetNumberOfParticipants(13.3);//Centrality class Per4 for PDC04
1274       gGener=gener;
1275     }
1276     break;
1277   case kFMD1Flat: 
1278     {
1279       AliGenBox* gener = new AliGenBox(2000);
1280       gener->SetPart(kPiPlus);
1281       gener->SetMomentumRange(3,4);
1282       gener->SetPhiRange(0, 360);
1283       gener->SetThetaRange(0.77, 3.08);
1284       gGener = gener;
1285     }
1286     break;
1287   case kFMD2Flat: 
1288     {
1289       AliGenBox* gener = new AliGenBox(100);
1290       gener->SetPart(kPiPlus);
1291       gener->SetMomentumRange(3,4);
1292       gener->SetPhiRange(0, 360);
1293       gener->SetThetaRange(2.95, 20.42);
1294       gGener = gener;
1295     }
1296     break;
1297   case kFMD3Flat: 
1298     {
1299       AliGenBox* gener = new AliGenBox(2000);
1300       gener->SetPart(kPiPlus);
1301       gener->SetMomentumRange(3,4);
1302       gener->SetPhiRange(0, 360);
1303       gener->SetThetaRange(155.97, 176.73);
1304       gGener = gener;
1305     }
1306     break;
1307   case kFMDFlat:
1308     {
1309       AliGenCocktail* gener = new AliGenCocktail();
1310       gener->SetMomentumRange(3,4);
1311       gener->SetPhiRange(0, 360);
1312       AliGenBox* gener3 = new AliGenBox(2000);
1313       gener3->SetThetaRange(155.97, 176.73);
1314       gener3->SetPart(kPiPlus);
1315       gener->AddGenerator(gener3, "FMD3", .33);
1316       AliGenBox* gener2 = new AliGenBox(2000);
1317       gener2->SetThetaRange(2.95, 20.42);
1318       gener2->SetPart(kPiPlus);
1319       gener->AddGenerator(gener2, "FMD2", .33);
1320       AliGenBox* gener1 = new AliGenBox(2000);
1321       gener1->SetThetaRange(0.77, 3.08);
1322       gener1->SetPart(kPiPlus);
1323       gener->AddGenerator(gener1, "FMD1", .34);
1324       gGener = gener;
1325     }
1326     break;
1327     
1328   default: break;
1329   }
1330   return gGener;
1331 }
1332
1333 //____________________________________________________________________
1334 AliGenHijing* 
1335 HijingStandard()
1336 {
1337   AliGenHijing *gener = new AliGenHijing(-1);
1338   // centre of mass energy 
1339   gener->SetEnergyCMS(5500.);
1340   // reference frame
1341   gener->SetReferenceFrame("CMS");
1342   // projectile
1343   gener->SetProjectile("A", 208, 82);
1344   gener->SetTarget    ("A", 208, 82);
1345   // tell hijing to keep the full parent child chain
1346   gener->KeepFullEvent();
1347   // enable jet quenching
1348   gener->SetJetQuenching(1);
1349   // enable shadowing
1350   gener->SetShadowing(1);
1351   // neutral pion and heavy particle decays switched off
1352   gener->SetDecaysOff(1);
1353   // Don't track spectators
1354   gener->SetSpectators(0);
1355   // kinematic selection
1356   gener->SetSelectAll(0);
1357   return gener;
1358 }
1359
1360
1361 //____________________________________________________________________
1362 void 
1363 ProcessEnvironmentVars(EG_t& eg, Int_t& seed)
1364 {
1365   // Run type
1366   if (gSystem->Getenv("CONFIG_RUN_TYPE")) {
1367     Int_t eg1 = LookupEG(gSystem->Getenv("CONFIG_RUN_TYPE"));
1368     if  (eg1 >= 0) eg = EG_t(eg1);
1369   }
1370   // Random Number seed
1371   if (gSystem->Getenv("CONFIG_SEED")) {
1372     seed = atoi(gSystem->Getenv("CONFIG_SEED"));
1373   }
1374 }
1375
1376 //____________________________________________________________________
1377 //
1378 // EOF
1379 //