3 // Configuration for the Geant4 production 2010
7 #if !defined(__CINT__) || defined(__MAKECINT__)
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"
50 kPythia6, kPythia6D6T, kPythia6D6T_Flat, kPythia6ATLAS, kPythia6ATLAS_Flat, kPythia6_Perugia0, kPhojet, kRunMax
53 const char * pprRunName[] = {
54 "kPythia6", "kPythia6D6T", "kPythia6D6T_Flat", "kPythia6ATLAS", "kPythia6ATLAS_Flat", "kPythia6_Perugia0", "kPhojet"
59 kNoField, k5kG, kFieldMax
62 const char * pprField[] = {
68 QGSP_BERT_EMV,CHIPS,QGSP_BERT_CHIPS,QGSP_BERT_EMV_OPTICAL, CHIPS_OPTICAL, QGSP_BERT_CHIPS_OPTICAL, kListMax
71 const char * physicsListName[] = {
72 "QGSP_BERT_EMV", "CHIPS", "QGSP_BERT_CHIPS",
73 "QGSP_BERT_EMV_OPTICAL", "CHIPS_OPTICAL", "QGSP_BERT_CHIPS_OPTICAL"
78 AliGenerator *MbPythia();
79 AliGenerator *MbPythiaTuneD6T();
80 AliGenerator *MbPhojet();
81 void ProcessEnvironmentVars();
83 // Geterator, field, beam energy
84 static PDC06Proc_t proc = kPhojet;
85 static Mag_t mag = k5kG;
86 static Float_t energy = 10000; // energy in CMS
87 static PhysicsList_t physicslist = QGSP_BERT_EMV;
88 //========================//
89 // Set Random Number seed //
90 //========================//
92 static UInt_t seed = dt.Get();
95 static TString comment;
101 // Get settings from environment variables
102 ProcessEnvironmentVars();
104 gRandom->SetSeed(seed);
105 cerr<<"Seed for random number generation= "<<seed<<endl;
107 // Libraries required by geant321
108 #if defined(__CINT__)
109 gSystem->Load("liblhapdf"); // Parton density functions
110 gSystem->Load("libEGPythia6"); // TGenerator interface
112 if (proc == kPythia6 || proc == kPhojet) {
113 gSystem->Load("libpythia6"); // Pythia 6.2
115 gSystem->Load("libpythia6.4.21"); // Pythia 6.4
118 gSystem->Load("libAliPythia6"); // ALICE specific implementations
119 //gSystem->Load("libgeant321");
122 // new TGeant3TGeo("C++ Interface to Geant3");
124 //=======================================================================
125 // Create the output file
128 AliRunLoader* rl=0x0;
130 cout<<"Config.C: Creating Run Loader ..."<<endl;
131 rl = AliRunLoader::Open("galice.root",
132 AliConfig::GetDefaultEventFolderName(),
136 gAlice->Fatal("Config.C","Can not instatiate the Run Loader");
139 rl->SetCompressionLevel(2);
140 rl->SetNumberOfEventsPerFile(1000);
141 gAlice->SetRunLoader(rl);
142 // gAlice->SetGeometryFromFile("geometry.root");
143 // gAlice->SetGeometryFromCDB();
145 // Set the trigger configuration: proton-proton
146 gAlice->SetTriggerDescriptor("p-p");
147 // AliSimulation::Instance()->SetTriggerConfig("p-p");
149 printf("\n \n Comment: %s \n \n", comment.Data());
176 //=================== Alice BODY parameters =============================
177 AliBODY *BODY = new AliBODY("BODY", "Alice envelop");
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");
191 //=================== ABSO parameters ============================
192 AliABSO *ABSO = new AliABSOv3("ABSO", "Muon Absorber");
197 //=================== DIPO parameters ============================
199 AliDIPO *DIPO = new AliDIPOv3("DIPO", "Dipole version 3");
204 //=================== HALL parameters ============================
206 AliHALL *HALL = new AliHALLv3("HALL", "Alice Hall");
212 //=================== FRAME parameters ============================
214 AliFRAMEv2 *FRAME = new AliFRAMEv2("FRAME", "Space Frame");
220 //=================== SHIL parameters ============================
222 AliSHIL *SHIL = new AliSHILv3("SHIL", "Shielding Version 3");
228 //=================== PIPE parameters ============================
230 AliPIPE *PIPE = new AliPIPEv3("PIPE", "Beam Pipe");
235 //=================== ITS parameters ============================
237 AliITS *ITS = new AliITSv11Hybrid("ITS","ITS v11Hybrid");
242 //============================ TPC parameters =====================
244 AliTPC *TPC = new AliTPCv2("TPC", "Default");
245 // TPC->SetPrimaryIonisation(); // not used with Geant3
250 //=================== TOF parameters ============================
252 AliTOF *TOF = new AliTOFv6T0("TOF", "normal TOF");
258 //=================== HMPID parameters ===========================
260 AliHMPID *HMPID = new AliHMPIDv3("HMPID", "normal HMPID");
267 //=================== ZDC parameters ============================
269 AliZDC *ZDC = new AliZDCv3("ZDC", "normal ZDC");
274 //=================== TRD parameters ============================
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);
295 //=================== FMD parameters ============================
297 AliFMD *FMD = new AliFMDv1("FMD", "normal FMD");
302 //=================== MUON parameters ===========================
303 // New MUONv1 version (geometry defined via builders)
305 AliMUON *MUON = new AliMUONv1("MUON", "default");
310 //=================== PHOS parameters ===========================
312 AliPHOS *PHOS = new AliPHOSv1("PHOS", "noCPV_Modules123");
318 //=================== PMD parameters ============================
320 AliPMD *PMD = new AliPMDv1("PMD", "normal PMD");
325 //=================== T0 parameters ============================
326 AliT0 *T0 = new AliT0v1("T0", "T0 Detector");
331 //=================== EMCAL parameters ============================
333 AliEMCAL *EMCAL = new AliEMCALv2("EMCAL", "EMCAL_FIRSTYEAR");
338 //=================== ACORDE parameters ============================
340 AliACORDE *ACORDE = new AliACORDEv1("ACORDE", "normal ACORDE");
345 //=================== ACORDE parameters ============================
347 AliVZERO *VZERO = new AliVZEROv7("VZERO", "normal VZERO");
352 // Load Geant4 + Geant4 VMC libraries
354 if (gClassTable->GetID("TGeant4") == -1) {
355 TString g4libsMacro = "$G4INSTALL/macro/g4libs.C";
356 //TString g4libsMacro = "$ALICE/geant4_vmc/examples/macro/g4libs.C";
357 // Load Geant4 libraries
358 if (!gInterpreter->IsLoaded(g4libsMacro.Data())) {
359 gROOT->LoadMacro(g4libsMacro.Data());
360 gInterpreter->ProcessLine("g4libs()");
369 TG4RunConfiguration* runConfiguration=0x0;
370 for (Int_t iList = 0; iList < kListMax; iList++) {
371 if(iList<kListMax/2){
372 if(physicslist == iList){
374 new TG4RunConfiguration("geomRoot",
375 physicsListName[iList],
376 "specialCuts+stackPopper+stepLimiter",
380 else if(iList>=kListMax/2){//add "optical" PL to HadronPhysicsList
381 if(physicslist == iList){
383 new TG4RunConfiguration("geomRoot",
384 Form("%s+optical",physicsListName[iList-kListMax/2]),
385 "specialCuts+stackPopper+stepLimiter",
390 geant4 = new TGeant4("TGeant4", "The Geant4 Monte Carlo", runConfiguration);
391 cout << "Geant4 has been created." << endl;
394 cout << "Monte Carlo has been already created." << endl;
399 // Customization of Geant4 VMC
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");
409 geant4->ProcessGeantCommand("/mcTracking/skipNeutrino true");
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");
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");
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
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);
448 Float_t cut = 1.e-3; // 1MeV cut by default
449 Float_t tofmax = 1.e10;
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);
466 //======================//
467 // Set External decayer //
468 //======================//
469 TVirtualMCDecayer* decayer = new AliDecayerPythia();
470 decayer->SetForceDecay(kAll);
472 gMC->SetExternalDecayer(decayer);
474 //=========================//
475 // Generator Configuration //
476 //=========================//
477 AliGenerator* gener = 0x0;
479 if (proc == kPythia6) {
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) {
491 } else if (proc == kPythia6_Perugia0) {
492 gener = MbPythiaTunePerugia0();
504 for (Int_t i = 0; i < 3; i++) {
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);
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]);
519 gener->SetOrigin(vtxPos[0], vtxPos[1], vtxPos[2]); // vertex position
522 // Size of the interaction diamond
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]
532 printf("\n \n Diamond size x-y: %10.3e z: %10.3e\n \n", sigmaxy, sigmaz);
534 gener->SetSigma(sigmaxy, sigmaxy, sigmaz); // Sigma in (X,Y,Z) (cm) on IP position
535 gener->SetVertexSmear(kPerEvent);
555 AliGenerator* MbPythia()
557 comment = comment.Append(" pp: Pythia low-pt");
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);
571 AliGenerator* MbPythiaTuneD6T()
573 comment = comment.Append(" pp: Pythia low-pt");
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);
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);
591 AliGenerator* MbPythiaTunePerugia0()
593 comment = comment.Append(" pp: Pythia Minimum Bias Tune Perugia0");
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);
604 pythia->SetTune(320);
609 AliGenerator* MbPythiaTuneATLAS()
611 comment = comment.Append(" pp: Pythia low-pt");
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);
622 // C 306 ATLAS-CSC: Arthur Moraes' (new) ATLAS tune (needs CTEQ6L externally)
623 pythia->SetTune(306);
624 pythia->SetStrucFunc(kCTEQ6l);
629 AliGenerator* MbPythiaTuneATLAS_Flat()
631 AliGenPythia* pythia = MbPythiaTuneATLAS();
633 comment = comment.Append("; flat multiplicity distribution");
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);
717 Int_t limit = weight->GetRandom();
718 pythia->SetTriggerChargedMultiplicity(limit, 1.4);
720 comment = comment.Append(Form("; multiplicity threshold set to %d in |eta| < 1.4", limit));
724 AliGenerator* MbPythiaTuneD6T_Flat()
726 AliGenPythia* pythia = MbPythiaTuneD6T();
728 comment = comment.Append("; flat multiplicity distribution");
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);
812 Int_t limit = weight->GetRandom();
813 pythia->SetTriggerChargedMultiplicity(limit, 1.4);
815 comment = comment.Append(Form("; multiplicity threshold set to %d in |eta| < 1.4", limit));
820 AliGenerator* MbPhojet()
822 comment = comment.Append(" pp: Pythia low-pt");
825 #if defined(__CINT__)
826 gSystem->Load("libdpmjet"); // Parton density functions
827 gSystem->Load("libTDPMjet"); // Parton density functions
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);
840 void ProcessEnvironmentVars()
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;
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) {
857 cout<<"Field set to "<<pprField[iField]<<endl;
863 if (gSystem->Getenv("CONFIG_ENERGY")) {
864 energy = atoi(gSystem->Getenv("CONFIG_ENERGY"));
865 cout<<"Energy set to "<<energy<<" GeV"<<endl;
868 // Random Number seed
869 if (gSystem->Getenv("CONFIG_SEED")) {
870 seed = atoi(gSystem->Getenv("CONFIG_SEED"));
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;