In vmctest/production:
[u/mrichter/AliRoot.git] / test / vmctest / production / Config.C
CommitLineData
ba2958c4 1// $Id$
b2b18d9e 2//
ba2958c4 3// Configuration for the Geant4 production 2010
4// By E. Sicking, CERN
b2b18d9e 5//
6
b2b18d9e 7#if !defined(__CINT__) || defined(__MAKECINT__)
8#include <Riostream.h>
9#include <TRandom.h>
10#include <TDatime.h>
11#include <TSystem.h>
12#include <TVirtualMC.h>
13#include <TGeant3TGeo.h>
14#include "STEER/AliRunLoader.h"
15#include "STEER/AliRun.h"
16#include "STEER/AliConfig.h"
17#include "PYTHIA6/AliDecayerPythia.h"
18#include "PYTHIA6/AliGenPythia.h"
19#include "TDPMjet/AliGenDPMjet.h"
20#include "STEER/AliMagFCheb.h"
21#include "STRUCT/AliBODY.h"
22#include "STRUCT/AliMAG.h"
23#include "STRUCT/AliABSOv3.h"
24#include "STRUCT/AliDIPOv3.h"
25#include "STRUCT/AliHALLv3.h"
26#include "STRUCT/AliFRAMEv2.h"
27#include "STRUCT/AliSHILv3.h"
28#include "STRUCT/AliPIPEv3.h"
29#include "ITS/AliITSv11Hybrid.h"
30#include "TPC/AliTPCv2.h"
31#include "TOF/AliTOFv6T0.h"
32#include "HMPID/AliHMPIDv3.h"
33#include "ZDC/AliZDCv3.h"
34#include "TRD/AliTRDv1.h"
35#include "TRD/AliTRDgeometry.h"
36#include "FMD/AliFMDv1.h"
37#include "MUON/AliMUONv1.h"
38#include "PHOS/AliPHOSv1.h"
39#include "PHOS/AliPHOSSimParam.h"
40#include "PMD/AliPMDv1.h"
41#include "T0/AliT0v1.h"
42#include "EMCAL/AliEMCALv2.h"
43#include "ACORDE/AliACORDEv1.h"
44#include "VZERO/AliVZEROv7.h"
45#endif
46
47
48enum PDC06Proc_t
49{
50 kPythia6, kPythia6D6T, kPythia6D6T_Flat, kPythia6ATLAS, kPythia6ATLAS_Flat, kPythia6_Perugia0, kPhojet, kRunMax
51};
52
53const char * pprRunName[] = {
54 "kPythia6", "kPythia6D6T", "kPythia6D6T_Flat", "kPythia6ATLAS", "kPythia6ATLAS_Flat", "kPythia6_Perugia0", "kPhojet"
55};
56
57enum Mag_t
58{
59 kNoField, k5kG, kFieldMax
60};
61
62const char * pprField[] = {
63 "kNoField", "k5kG"
64};
65
66enum PhysicsList_t
67{
68 QGSP_BERT_EMV,CHIPS,QGSP_BERT_CHIPS,QGSP_BERT_EMV_OPTICAL, CHIPS_OPTICAL, QGSP_BERT_CHIPS_OPTICAL, kListMax
69};
70
71const char * physicsListName[] = {
72 "QGSP_BERT_EMV", "CHIPS", "QGSP_BERT_CHIPS",
73 "QGSP_BERT_EMV_OPTICAL", "CHIPS_OPTICAL", "QGSP_BERT_CHIPS_OPTICAL"
74};
75
76//--- Functions ---
77class AliGenPythia;
78AliGenerator *MbPythia();
79AliGenerator *MbPythiaTuneD6T();
80AliGenerator *MbPhojet();
81void ProcessEnvironmentVars();
82
83// Geterator, field, beam energy
84static PDC06Proc_t proc = kPhojet;
85static Mag_t mag = k5kG;
86static Float_t energy = 10000; // energy in CMS
87static PhysicsList_t physicslist = QGSP_BERT_EMV;
88//========================//
89// Set Random Number seed //
90//========================//
91TDatime dt;
92static UInt_t seed = dt.Get();
93
94// Comment line
95static TString comment;
96
97void Config()
98{
99
100
101 // Get settings from environment variables
102 ProcessEnvironmentVars();
103
104 gRandom->SetSeed(seed);
105 cerr<<"Seed for random number generation= "<<seed<<endl;
106
107 // Libraries required by geant321
108#if defined(__CINT__)
109 gSystem->Load("liblhapdf"); // Parton density functions
110 gSystem->Load("libEGPythia6"); // TGenerator interface
111
112 if (proc == kPythia6 || proc == kPhojet) {
113 gSystem->Load("libpythia6"); // Pythia 6.2
114 } else {
115 gSystem->Load("libpythia6.4.21"); // Pythia 6.4
116 }
117
118 gSystem->Load("libAliPythia6"); // ALICE specific implementations
119 //gSystem->Load("libgeant321");
120#endif
121
122 // new TGeant3TGeo("C++ Interface to Geant3");
123
124 //=======================================================================
125 // Create the output file
126
127
128 AliRunLoader* rl=0x0;
129
130 cout<<"Config.C: Creating Run Loader ..."<<endl;
131 rl = AliRunLoader::Open("galice.root",
132 AliConfig::GetDefaultEventFolderName(),
133 "recreate");
134 if (rl == 0x0)
135 {
136 gAlice->Fatal("Config.C","Can not instatiate the Run Loader");
137 return;
138 }
139 rl->SetCompressionLevel(2);
140 rl->SetNumberOfEventsPerFile(1000);
141 gAlice->SetRunLoader(rl);
142 // gAlice->SetGeometryFromFile("geometry.root");
143 // gAlice->SetGeometryFromCDB();
144
145 // Set the trigger configuration: proton-proton
146 gAlice->SetTriggerDescriptor("p-p");
147 // AliSimulation::Instance()->SetTriggerConfig("p-p");
148
b2b18d9e 149 printf("\n \n Comment: %s \n \n", comment.Data());
fd03ab3d 150
b2b18d9e 151 rl->CdGAFile();
152
153 Int_t iABSO = 1;
154 Int_t iACORDE= 0;
155 Int_t iDIPO = 1;
156 Int_t iEMCAL = 1;
157 Int_t iFMD = 1;
158 Int_t iFRAME = 1;
159 Int_t iHALL = 1;
160 Int_t iITS = 1;
161 Int_t iMAG = 1;
162 Int_t iMUON = 1;
163 Int_t iPHOS = 1;
164 Int_t iPIPE = 1;
165 Int_t iPMD = 1;
166 Int_t iHMPID = 1;
167 Int_t iSHIL = 1;
168 Int_t iT0 = 1;
169 Int_t iTOF = 1;
170 Int_t iTPC = 1;
171 Int_t iTRD = 1;
172 Int_t iVZERO = 1;
173 Int_t iZDC = 1;
174
175
176 //=================== Alice BODY parameters =============================
177 AliBODY *BODY = new AliBODY("BODY", "Alice envelop");
178
179
180 if (iMAG)
181 {
182 //=================== MAG parameters ============================
183 // --- Start with Magnet since detector layouts may be depending ---
184 // --- on the selected Magnet dimensions ---
185 AliMAG *MAG = new AliMAG("MAG", "Magnet");
186 }
187
188
189 if (iABSO)
190 {
191 //=================== ABSO parameters ============================
192 AliABSO *ABSO = new AliABSOv3("ABSO", "Muon Absorber");
193 }
194
195 if (iDIPO)
196 {
197 //=================== DIPO parameters ============================
198
199 AliDIPO *DIPO = new AliDIPOv3("DIPO", "Dipole version 3");
200 }
201
202 if (iHALL)
203 {
204 //=================== HALL parameters ============================
205
206 AliHALL *HALL = new AliHALLv3("HALL", "Alice Hall");
207 }
208
209
210 if (iFRAME)
211 {
212 //=================== FRAME parameters ============================
213
214 AliFRAMEv2 *FRAME = new AliFRAMEv2("FRAME", "Space Frame");
215 FRAME->SetHoles(1);
216 }
217
218 if (iSHIL)
219 {
220 //=================== SHIL parameters ============================
221
222 AliSHIL *SHIL = new AliSHILv3("SHIL", "Shielding Version 3");
223 }
224
225
226 if (iPIPE)
227 {
228 //=================== PIPE parameters ============================
229
230 AliPIPE *PIPE = new AliPIPEv3("PIPE", "Beam Pipe");
231 }
232
233 if (iITS)
234 {
235 //=================== ITS parameters ============================
236
237 AliITS *ITS = new AliITSv11Hybrid("ITS","ITS v11Hybrid");
238 }
239
240 if (iTPC)
241 {
242 //============================ TPC parameters =====================
243
244 AliTPC *TPC = new AliTPCv2("TPC", "Default");
fd03ab3d 245 // TPC->SetPrimaryIonisation(); // not used with Geant3
b2b18d9e 246 }
247
248
249 if (iTOF) {
250 //=================== TOF parameters ============================
251
252 AliTOF *TOF = new AliTOFv6T0("TOF", "normal TOF");
253 }
254
255
256 if (iHMPID)
257 {
258 //=================== HMPID parameters ===========================
259
260 AliHMPID *HMPID = new AliHMPIDv3("HMPID", "normal HMPID");
261
262 }
263
264
265 if (iZDC)
266 {
267 //=================== ZDC parameters ============================
268
269 AliZDC *ZDC = new AliZDCv3("ZDC", "normal ZDC");
270 }
271
272 if (iTRD)
273 {
274 //=================== TRD parameters ============================
275
276 AliTRD *TRD = new AliTRDv1("TRD", "TRD slow simulator");
277 AliTRDgeometry *geoTRD = TRD->GetGeometry();
278 // Partial geometry: modules at 0,1,7,8,9,16,17
279 // starting at 3h in positive direction
280 geoTRD->SetSMstatus(2,0);
281 geoTRD->SetSMstatus(3,0);
282 geoTRD->SetSMstatus(4,0);
283 geoTRD->SetSMstatus(5,0);
284 geoTRD->SetSMstatus(6,0);
285 geoTRD->SetSMstatus(11,0);
286 geoTRD->SetSMstatus(12,0);
287 geoTRD->SetSMstatus(13,0);
288 geoTRD->SetSMstatus(14,0);
289 geoTRD->SetSMstatus(15,0);
290 geoTRD->SetSMstatus(16,0);
291 }
292
293 if (iFMD)
294 {
295 //=================== FMD parameters ============================
296
297 AliFMD *FMD = new AliFMDv1("FMD", "normal FMD");
298 }
299
300 if (iMUON)
301 {
302 //=================== MUON parameters ===========================
303 // New MUONv1 version (geometry defined via builders)
304
305 AliMUON *MUON = new AliMUONv1("MUON", "default");
306 }
307
308 if (iPHOS)
309 {
310 //=================== PHOS parameters ===========================
311
fd03ab3d 312 AliPHOS *PHOS = new AliPHOSv1("PHOS", "noCPV_Modules123");
b2b18d9e 313
b2b18d9e 314 }
315
b2b18d9e 316 if (iPMD)
317 {
318 //=================== PMD parameters ============================
319
320 AliPMD *PMD = new AliPMDv1("PMD", "normal PMD");
321 }
322
323 if (iT0)
324 {
325 //=================== T0 parameters ============================
326 AliT0 *T0 = new AliT0v1("T0", "T0 Detector");
327 }
328
329 if (iEMCAL)
330 {
331 //=================== EMCAL parameters ============================
332
333 AliEMCAL *EMCAL = new AliEMCALv2("EMCAL", "EMCAL_FIRSTYEAR");
334 }
335
336 if (iACORDE)
337 {
338 //=================== ACORDE parameters ============================
339
340 AliACORDE *ACORDE = new AliACORDEv1("ACORDE", "normal ACORDE");
341 }
342
343 if (iVZERO)
344 {
345 //=================== ACORDE parameters ============================
346
347 AliVZERO *VZERO = new AliVZEROv7("VZERO", "normal VZERO");
348 }
349
350
351
352 // Load Geant4 + Geant4 VMC libraries
353 //
354 if (gClassTable->GetID("TGeant4") == -1) {
ba2958c4 355 TString g4libsMacro = "$G4INSTALL/macro/g4libs.C";
356 //TString g4libsMacro = "$ALICE/geant4_vmc/examples/macro/g4libs.C";
b2b18d9e 357 // Load Geant4 libraries
ba2958c4 358 if (!gInterpreter->IsLoaded(g4libsMacro.Data())) {
359 gROOT->LoadMacro(g4libsMacro.Data());
b2b18d9e 360 gInterpreter->ProcessLine("g4libs()");
361 }
362 }
363
364
365 // Create Geant4 VMC
366 //
367 TGeant4 *geant4 = 0;
368 if ( ! gMC ) {
369 TG4RunConfiguration* runConfiguration=0x0;
370 for (Int_t iList = 0; iList < kListMax; iList++) {
371 if(iList<kListMax/2){
372 if(physicslist == iList){
373 runConfiguration =
374 new TG4RunConfiguration("geomRoot",
375 physicsListName[iList],
376 "specialCuts+stackPopper+stepLimiter",
377 true);
378 }
379 }
380 else if(iList>=kListMax/2){//add "optical" PL to HadronPhysicsList
381 if(physicslist == iList){
382 runConfiguration =
383 new TG4RunConfiguration("geomRoot",
384 Form("%s+optical",physicsListName[iList-kListMax/2]),
385 "specialCuts+stackPopper+stepLimiter",
386 true);
387 }
388 }
389 }
390 geant4 = new TGeant4("TGeant4", "The Geant4 Monte Carlo", runConfiguration);
391 cout << "Geant4 has been created." << endl;
392 }
393 else {
394 cout << "Monte Carlo has been already created." << endl;
395 }
396
397
398
399 // Customization of Geant4 VMC
400 //
401 geant4->ProcessGeantCommand("/mcVerbose/all 1");
402 geant4->ProcessGeantCommand("/mcVerbose/geometryManager 1");
403 geant4->ProcessGeantCommand("/mcVerbose/opGeometryManager 1");
404 geant4->ProcessGeantCommand("/mcTracking/loopVerbose 0");
405 geant4->ProcessGeantCommand("/mcPhysics/rangeCuts 0.01 mm");
406 geant4->ProcessGeantCommand("/mcPhysics/selectOpProcess Scintillation");
407 geant4->ProcessGeantCommand("/mcPhysics/setOpProcessActivation false");
408 geant4->ProcessGeantCommand("/mcVerbose/composedPhysicsList 2");
e6918581 409 geant4->ProcessGeantCommand("/mcTracking/skipNeutrino true");
b2b18d9e 410
fd03ab3d 411 /*
412 // Set PAI model for TPC (TPC_Ne-CO2-N-2)
413 geant4->ProcessGeantCommand("/mcPhysics/emModel/selectMedium 219");
414 geant4->ProcessGeantCommand("/mcPhysics/emModel/setElossModel PAI");
415 geant4->ProcessGeantCommand("/mcPhysics/emModel/setFluctModel PAI");
416 geant4->ProcessGeantCommand("/mcPhysics/emModel/setParticles all");
417
418 // Set PAI model for TRD (TRD_XeCO2)
419 geant4->ProcessGeantCommand("/mcPhysics/emModel/selectMedium 291");
420 geant4->ProcessGeantCommand("/mcPhysics/emModel/setElossModel PAI");
421 geant4->ProcessGeantCommand("/mcPhysics/emModel/setFluctModel PAI");
422 geant4->ProcessGeantCommand("/mcPhysics/emModel/setParticles all");
423 */
b2b18d9e 424
425
426 //
427 //=======================================================================
428 // ************* STEERING parameters FOR ALICE SIMULATION **************
429 // --- Specify event type to be tracked through the ALICE setup
430 // --- All positions are in cm, angles in degrees, and P and E in GeV
431
432
433 gMC->SetProcess("DCAY",1);
434 gMC->SetProcess("PAIR",1);
435 gMC->SetProcess("COMP",1);
436 gMC->SetProcess("PHOT",1);
437 gMC->SetProcess("PFIS",0);
438 gMC->SetProcess("DRAY",0);
439 gMC->SetProcess("ANNI",1);
440 gMC->SetProcess("BREM",1);
441 gMC->SetProcess("MUNU",1);
442 gMC->SetProcess("CKOV",1);
443 gMC->SetProcess("HADR",1);
444 gMC->SetProcess("LOSS",2);
445 gMC->SetProcess("MULS",1);
446 gMC->SetProcess("RAYL",1);
447
448 Float_t cut = 1.e-3; // 1MeV cut by default
449 Float_t tofmax = 1.e10;
450
451 gMC->SetCut("CUTGAM", cut);
452 gMC->SetCut("CUTELE", cut);
453 gMC->SetCut("CUTNEU", cut);
454 gMC->SetCut("CUTHAD", cut);
455 gMC->SetCut("CUTMUO", cut);
456 gMC->SetCut("BCUTE", cut);
457 gMC->SetCut("BCUTM", cut);
458 gMC->SetCut("DCUTE", cut);
459 gMC->SetCut("DCUTM", cut);
460 gMC->SetCut("PPCUTM", cut);
461 gMC->SetCut("TOFMAX", tofmax);
462
463
464
465
466 //======================//
467 // Set External decayer //
468 //======================//
469 TVirtualMCDecayer* decayer = new AliDecayerPythia();
470 decayer->SetForceDecay(kAll);
471 decayer->Init();
472 gMC->SetExternalDecayer(decayer);
473
474 //=========================//
475 // Generator Configuration //
476 //=========================//
477 AliGenerator* gener = 0x0;
478
479 if (proc == kPythia6) {
480 gener = MbPythia();
481 } else if (proc == kPythia6D6T) {
482 gener = MbPythiaTuneD6T();
483 } else if (proc == kPythia6D6T_Flat) {
484 gener = MbPythiaTuneD6T_Flat();
485 } else if (proc == kPythia6ATLAS) {
486 gener = MbPythiaTuneATLAS();
487 } else if (proc == kPythia6ATLAS_Flat) {
488 gener = MbPythiaTuneATLAS_Flat();
489 } else if (proc == kPhojet) {
490 gener = MbPhojet();
491 } else if (proc == kPythia6_Perugia0) {
492 gener = MbPythiaTunePerugia0();
493 }
494
495 //
496 // PRIMARY VERTEX
497 //
498 Double_t xv = 0;
499 Double_t yv = 0;
500 Double_t zv = 0;
501
502 Double_t vtxPos[3];
503 Double_t vtxErr[3];
504 for (Int_t i = 0; i < 3; i++) {
505 vtxPos[i] = 0.;
506 vtxErr[i] = 0.;
507 }
508
509 if(AliCDBManager::Instance()->IsDefaultStorageSet()){
510 AliCDBEntry* entry = AliCDBManager::Instance()->Get("GRP/Calib/MeanVertexSPD");
511 AliESDVertex* vertex = dynamic_cast<AliESDVertex*> (entry->GetObject());
512 vertex->GetXYZ(vtxPos);
513 vertex->GetSigmaXYZ(vtxErr);
514 }
515
516 printf("Vertex position from OCDB entry: x = %13.3f, y = %13.3f, z = %13.3f (sigma = %13.3f)\n",
517 vtxPos[0], vtxPos[1], vtxPos[2], vtxErr[2]);
518
519 gener->SetOrigin(vtxPos[0], vtxPos[1], vtxPos[2]); // vertex position
520 //
521 //
522 // Size of the interaction diamond
523 // Longitudinal
524 sigmaz = vtxErr[2];
525//
526 // Transverse
527 Float_t betast = 10; // beta* [m]
528 Float_t eps = 2.5e-6; // emittance [m]
529 Float_t gamma = energy / 2.0 / 0.938272; // relativistic gamma [1]
530 Float_t sigmaxy = TMath::Sqrt(eps * betast / gamma) / TMath::Sqrt(2.) * 100.; // [cm]
531
532 printf("\n \n Diamond size x-y: %10.3e z: %10.3e\n \n", sigmaxy, sigmaz);
533
534 gener->SetSigma(sigmaxy, sigmaxy, sigmaz); // Sigma in (X,Y,Z) (cm) on IP position
535 gener->SetVertexSmear(kPerEvent);
536 gener->Init();
537
538
539
540
541
542
543
544
545
546
547
548
549
550}
551//
552// PYTHIA
553//
554
555AliGenerator* MbPythia()
556{
557 comment = comment.Append(" pp: Pythia low-pt");
558//
559// Pythia
560 AliGenPythia* pythia = new AliGenPythia(-1);
561 pythia->SetMomentumRange(0, 999999.);
562 pythia->SetThetaRange(0., 180.);
563 pythia->SetYRange(-12.,12.);
564 pythia->SetPtRange(0,1000.);
565 pythia->SetProcess(kPyMb);
566 pythia->SetEnergyCMS(energy);
567
568 return pythia;
569}
570
571AliGenerator* MbPythiaTuneD6T()
572{
573 comment = comment.Append(" pp: Pythia low-pt");
574//
575// Pythia
576 AliGenPythia* pythia = new AliGenPythia(-1);
577 pythia->SetMomentumRange(0, 999999.);
578 pythia->SetThetaRange(0., 180.);
579 pythia->SetYRange(-12.,12.);
580 pythia->SetPtRange(0,1000.);
581 pythia->SetProcess(kPyMb);
582 pythia->SetEnergyCMS(energy);
583// Tune
584// 109 D6T : Rick Field's CDF Tune D6T (NB: needs CTEQ6L pdfs externally)
585 pythia->SetTune(109); // F I X
586 pythia->SetStrucFunc(kCTEQ6l);
587//
588 return pythia;
589}
590
591AliGenerator* MbPythiaTunePerugia0()
592{
593 comment = comment.Append(" pp: Pythia Minimum Bias Tune Perugia0");
594//
595// Pythia
596 AliGenPythia* pythia = new AliGenPythia(-1);
597 pythia->SetMomentumRange(0, 999999.);
598 pythia->SetThetaRange(0., 180.);
599 pythia->SetYRange(-12.,12.);
600 pythia->SetPtRange(0, 1000.);
601 pythia->SetProcess(kPyMb);
602 pythia->SetEnergyCMS(energy);
603// Tune
604 pythia->SetTune(320);
605//
606 return pythia;
607}
608
609AliGenerator* MbPythiaTuneATLAS()
610{
611 comment = comment.Append(" pp: Pythia low-pt");
612//
613// Pythia
614 AliGenPythia* pythia = new AliGenPythia(-1);
615 pythia->SetMomentumRange(0, 999999.);
616 pythia->SetThetaRange(0., 180.);
617 pythia->SetYRange(-12.,12.);
618 pythia->SetPtRange(0,1000.);
619 pythia->SetProcess(kPyMb);
620 pythia->SetEnergyCMS(energy);
621// Tune
622// C 306 ATLAS-CSC: Arthur Moraes' (new) ATLAS tune (needs CTEQ6L externally)
623 pythia->SetTune(306);
624 pythia->SetStrucFunc(kCTEQ6l);
625//
626 return pythia;
627}
628
629AliGenerator* MbPythiaTuneATLAS_Flat()
630{
631 AliGenPythia* pythia = MbPythiaTuneATLAS();
632
633 comment = comment.Append("; flat multiplicity distribution");
634
635 // set high multiplicity trigger
636 // this weight achieves a flat multiplicity distribution
637 TH1 *weight = new TH1D("weight","weight",201,-0.5,200.5);
638 weight->SetBinContent(1,5.49443);
639 weight->SetBinContent(2,8.770816);
640 weight->SetBinContent(6,0.4568624);
641 weight->SetBinContent(7,0.2919915);
642 weight->SetBinContent(8,0.6674189);
643 weight->SetBinContent(9,0.364737);
644 weight->SetBinContent(10,0.8818444);
645 weight->SetBinContent(11,0.531885);
646 weight->SetBinContent(12,1.035197);
647 weight->SetBinContent(13,0.9394057);
648 weight->SetBinContent(14,0.9643193);
649 weight->SetBinContent(15,0.94543);
650 weight->SetBinContent(16,0.9426507);
651 weight->SetBinContent(17,0.9423649);
652 weight->SetBinContent(18,0.789456);
653 weight->SetBinContent(19,1.149026);
654 weight->SetBinContent(20,1.100491);
655 weight->SetBinContent(21,0.6350525);
656 weight->SetBinContent(22,1.351941);
657 weight->SetBinContent(23,0.03233504);
658 weight->SetBinContent(24,0.9574557);
659 weight->SetBinContent(25,0.868133);
660 weight->SetBinContent(26,1.030998);
661 weight->SetBinContent(27,1.08897);
662 weight->SetBinContent(28,1.251382);
663 weight->SetBinContent(29,0.1391099);
664 weight->SetBinContent(30,1.192876);
665 weight->SetBinContent(31,0.448944);
666 weight->SetBinContent(32,1);
667 weight->SetBinContent(33,1);
668 weight->SetBinContent(34,1);
669 weight->SetBinContent(35,1);
670 weight->SetBinContent(36,0.9999997);
671 weight->SetBinContent(37,0.9999997);
672 weight->SetBinContent(38,0.9999996);
673 weight->SetBinContent(39,0.9999996);
674 weight->SetBinContent(40,0.9999995);
675 weight->SetBinContent(41,0.9999993);
676 weight->SetBinContent(42,1);
677 weight->SetBinContent(43,1);
678 weight->SetBinContent(44,1);
679 weight->SetBinContent(45,1);
680 weight->SetBinContent(46,1);
681 weight->SetBinContent(47,0.9999999);
682 weight->SetBinContent(48,0.9999998);
683 weight->SetBinContent(49,0.9999998);
684 weight->SetBinContent(50,0.9999999);
685 weight->SetBinContent(51,0.9999999);
686 weight->SetBinContent(52,0.9999999);
687 weight->SetBinContent(53,0.9999999);
688 weight->SetBinContent(54,0.9999998);
689 weight->SetBinContent(55,0.9999998);
690 weight->SetBinContent(56,0.9999998);
691 weight->SetBinContent(57,0.9999997);
692 weight->SetBinContent(58,0.9999996);
693 weight->SetBinContent(59,0.9999995);
694 weight->SetBinContent(60,1);
695 weight->SetBinContent(61,1);
696 weight->SetBinContent(62,1);
697 weight->SetBinContent(63,1);
698 weight->SetBinContent(64,1);
699 weight->SetBinContent(65,0.9999999);
700 weight->SetBinContent(66,0.9999998);
701 weight->SetBinContent(67,0.9999998);
702 weight->SetBinContent(68,0.9999999);
703 weight->SetBinContent(69,1);
704 weight->SetBinContent(70,1);
705 weight->SetBinContent(71,0.9999997);
706 weight->SetBinContent(72,0.9999995);
707 weight->SetBinContent(73,0.9999994);
708 weight->SetBinContent(74,1);
709 weight->SetBinContent(75,1);
710 weight->SetBinContent(76,1);
711 weight->SetBinContent(77,1);
712 weight->SetBinContent(78,0.9999999);
713 weight->SetBinContent(79,1);
714 weight->SetBinContent(80,1);
715 weight->SetEntries(526);
716
717 Int_t limit = weight->GetRandom();
718 pythia->SetTriggerChargedMultiplicity(limit, 1.4);
719
720 comment = comment.Append(Form("; multiplicity threshold set to %d in |eta| < 1.4", limit));
721
722 return pythia;
723}
724AliGenerator* MbPythiaTuneD6T_Flat()
725{
726 AliGenPythia* pythia = MbPythiaTuneD6T();
727
728 comment = comment.Append("; flat multiplicity distribution");
729
730 // set high multiplicity trigger
731 // this weight achieves a flat multiplicity distribution
732 TH1 *weight = new TH1D("weight","weight",201,-0.5,200.5);
733 weight->SetBinContent(1,5.49443);
734 weight->SetBinContent(2,8.770816);
735 weight->SetBinContent(6,0.4568624);
736 weight->SetBinContent(7,0.2919915);
737 weight->SetBinContent(8,0.6674189);
738 weight->SetBinContent(9,0.364737);
739 weight->SetBinContent(10,0.8818444);
740 weight->SetBinContent(11,0.531885);
741 weight->SetBinContent(12,1.035197);
742 weight->SetBinContent(13,0.9394057);
743 weight->SetBinContent(14,0.9643193);
744 weight->SetBinContent(15,0.94543);
745 weight->SetBinContent(16,0.9426507);
746 weight->SetBinContent(17,0.9423649);
747 weight->SetBinContent(18,0.789456);
748 weight->SetBinContent(19,1.149026);
749 weight->SetBinContent(20,1.100491);
750 weight->SetBinContent(21,0.6350525);
751 weight->SetBinContent(22,1.351941);
752 weight->SetBinContent(23,0.03233504);
753 weight->SetBinContent(24,0.9574557);
754 weight->SetBinContent(25,0.868133);
755 weight->SetBinContent(26,1.030998);
756 weight->SetBinContent(27,1.08897);
757 weight->SetBinContent(28,1.251382);
758 weight->SetBinContent(29,0.1391099);
759 weight->SetBinContent(30,1.192876);
760 weight->SetBinContent(31,0.448944);
761 weight->SetBinContent(32,1);
762 weight->SetBinContent(33,1);
763 weight->SetBinContent(34,1);
764 weight->SetBinContent(35,1);
765 weight->SetBinContent(36,0.9999997);
766 weight->SetBinContent(37,0.9999997);
767 weight->SetBinContent(38,0.9999996);
768 weight->SetBinContent(39,0.9999996);
769 weight->SetBinContent(40,0.9999995);
770 weight->SetBinContent(41,0.9999993);
771 weight->SetBinContent(42,1);
772 weight->SetBinContent(43,1);
773 weight->SetBinContent(44,1);
774 weight->SetBinContent(45,1);
775 weight->SetBinContent(46,1);
776 weight->SetBinContent(47,0.9999999);
777 weight->SetBinContent(48,0.9999998);
778 weight->SetBinContent(49,0.9999998);
779 weight->SetBinContent(50,0.9999999);
780 weight->SetBinContent(51,0.9999999);
781 weight->SetBinContent(52,0.9999999);
782 weight->SetBinContent(53,0.9999999);
783 weight->SetBinContent(54,0.9999998);
784 weight->SetBinContent(55,0.9999998);
785 weight->SetBinContent(56,0.9999998);
786 weight->SetBinContent(57,0.9999997);
787 weight->SetBinContent(58,0.9999996);
788 weight->SetBinContent(59,0.9999995);
789 weight->SetBinContent(60,1);
790 weight->SetBinContent(61,1);
791 weight->SetBinContent(62,1);
792 weight->SetBinContent(63,1);
793 weight->SetBinContent(64,1);
794 weight->SetBinContent(65,0.9999999);
795 weight->SetBinContent(66,0.9999998);
796 weight->SetBinContent(67,0.9999998);
797 weight->SetBinContent(68,0.9999999);
798 weight->SetBinContent(69,1);
799 weight->SetBinContent(70,1);
800 weight->SetBinContent(71,0.9999997);
801 weight->SetBinContent(72,0.9999995);
802 weight->SetBinContent(73,0.9999994);
803 weight->SetBinContent(74,1);
804 weight->SetBinContent(75,1);
805 weight->SetBinContent(76,1);
806 weight->SetBinContent(77,1);
807 weight->SetBinContent(78,0.9999999);
808 weight->SetBinContent(79,1);
809 weight->SetBinContent(80,1);
810 weight->SetEntries(526);
811
812 Int_t limit = weight->GetRandom();
813 pythia->SetTriggerChargedMultiplicity(limit, 1.4);
814
815 comment = comment.Append(Form("; multiplicity threshold set to %d in |eta| < 1.4", limit));
816
817 return pythia;
818}
819
820AliGenerator* MbPhojet()
821{
822 comment = comment.Append(" pp: Pythia low-pt");
823//
824// DPMJET
825#if defined(__CINT__)
826 gSystem->Load("libdpmjet"); // Parton density functions
827 gSystem->Load("libTDPMjet"); // Parton density functions
828#endif
829 AliGenDPMjet* dpmjet = new AliGenDPMjet(-1);
830 dpmjet->SetMomentumRange(0, 999999.);
831 dpmjet->SetThetaRange(0., 180.);
832 dpmjet->SetYRange(-12.,12.);
833 dpmjet->SetPtRange(0,1000.);
834 dpmjet->SetProcess(kDpmMb);
835 dpmjet->SetEnergyCMS(energy);
836
837 return dpmjet;
838}
839
840void ProcessEnvironmentVars()
841{
842 // Run type
843 if (gSystem->Getenv("CONFIG_RUN_TYPE")) {
844 for (Int_t iRun = 0; iRun < kRunMax; iRun++) {
845 if (strcmp(gSystem->Getenv("CONFIG_RUN_TYPE"), pprRunName[iRun])==0) {
846 proc = (PDC06Proc_t)iRun;
847 cout<<"Run type set to "<<pprRunName[iRun]<<endl;
848 }
849 }
850 }
851
852 // Field
853 if (gSystem->Getenv("CONFIG_FIELD")) {
854 for (Int_t iField = 0; iField < kFieldMax; iField++) {
855 if (strcmp(gSystem->Getenv("CONFIG_FIELD"), pprField[iField])==0) {
856 mag = (Mag_t)iField;
857 cout<<"Field set to "<<pprField[iField]<<endl;
858 }
859 }
860 }
861
862 // Energy
863 if (gSystem->Getenv("CONFIG_ENERGY")) {
864 energy = atoi(gSystem->Getenv("CONFIG_ENERGY"));
865 cout<<"Energy set to "<<energy<<" GeV"<<endl;
866 }
867
868 // Random Number seed
869 if (gSystem->Getenv("CONFIG_SEED")) {
870 seed = atoi(gSystem->Getenv("CONFIG_SEED"));
871 }
872
873 // Physics lists
874 if (gSystem->Getenv("CONFIG_PHYSICSLIST")) {
875 for (Int_t iList = 0; iList < kListMax; iList++) {
876 if (strcmp(gSystem->Getenv("CONFIG_PHYSICSLIST"), physicsListName[iList])==0) {
877 physicslist = (PhysicsList_t)iList;
878 cout<<"Physics list set to "<< physicsListName[iList]<<endl;
879 }
880 }
881 }
882
883}