1 /// One can use the configuration macro in compiled mode by
2 // root [0] gSystem->Load("libgeant321");
3 // root [0] gSystem->SetIncludePath("-I$ROOTSYS/include -I$ALICE_ROOT/include\
4 // -I$ALICE_ROOT -I$ALICE/geant3/TGeant3");
5 // root [0] .x grun.C(1,"Config.C++"
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 "EVGEN/AliGenCocktail.h"
19 #include "EVGEN/AliGenHIJINGpara.h"
20 #include "EVGEN/AliGenFixed.h"
21 #include "EVGEN/AliGenBox.h"
22 #include "STEER/AliMagWrapCheb.h"
23 #include "STRUCT/AliBODY.h"
24 #include "STRUCT/AliMAG.h"
25 #include "STRUCT/AliABSOv3.h"
26 #include "STRUCT/AliDIPOv3.h"
27 #include "STRUCT/AliHALLv3.h"
28 #include "STRUCT/AliFRAMEv2.h"
29 #include "STRUCT/AliSHILv3.h"
30 #include "STRUCT/AliPIPEv3.h"
31 #include "ITS/AliITSv11.h"
32 #include "ITS/UPGRADE/AliITSUv11.h"
33 #include "TPC/AliTPCv2.h"
34 #include "TOF/AliTOFv6T0.h"
35 #include "HMPID/AliHMPIDv3.h"
36 #include "ZDC/AliZDCv3.h"
37 #include "TRD/AliTRDv1.h"
38 #include "TRD/AliTRDgeometry.h"
39 #include "FMD/AliFMDv1.h"
40 #include "MUON/AliMUONv1.h"
41 #include "PHOS/AliPHOSv1.h"
42 #include "PMD/AliPMDv1.h"
43 #include "T0/AliT0v1.h"
44 #include "EMCAL/AliEMCALv2.h"
45 #include "ACORDE/AliACORDEv1.h"
46 #include "VZERO/AliVZEROv7.h"
47 #include <TVirtualMagField.h>
50 /* $Id: Config.C 47147 2011-02-07 11:46:44Z amastros $ */
53 kDefaultPPTrig, kDefaultPbPbTrig
56 const char * pprTrigConfName[] = {
60 Float_t EtaToTheta(Float_t arg);
62 static PprTrigConf_t strig = kDefaultPPTrig;// default PP trigger configuration
66 // ThetaRange is (0., 180.). It was (0.28,179.72) 7/12/00 09:00
67 // Theta range given through pseudorapidity limits 22/6/2001
69 // Set Random Number seed
70 gRandom->SetSeed(0); // Set 0 to use the currecnt time
73 // libraries required by geant321
75 gSystem->Load("liblhapdf");
76 gSystem->Load("libEGPythia6");
77 gSystem->Load("libpythia6");
78 gSystem->Load("libAliPythia6");
79 gSystem->Load("libgeant321");
81 gSystem->Load("libITSUpgradeBase.so");
82 gSystem->Load("libITSUpgradeSim.so");
85 new TGeant3TGeo("C++ Interface to Geant3");
90 rl = AliRunLoader::Open("galice.root",
91 AliConfig::GetDefaultEventFolderName(),
95 gAlice->Fatal("Config.C","Can not instatiate the Run Loader");
98 rl->SetCompressionLevel(2);
99 rl->SetNumberOfEventsPerFile(2000);
100 gAlice->SetRunLoader(rl);
102 // Set the trigger configuration
103 // gAlice->SetTriggerDescriptor(pprTrigConfName[strig]);
104 //cout<<"Trigger configuration is set to "<<pprTrigConfName[strig]<<endl;
105 AliSimulation::Instance()->SetTriggerConfig(pprTrigConfName[strig]);
106 cout<<"Trigger configuration is set to pprTrigConfName[strig] "<<endl;
109 // Set External decayer
110 TVirtualMCDecayer *decayer = new AliDecayerPythia();
112 decayer->SetForceDecay(kAll);
114 gMC->SetExternalDecayer(decayer);
115 //=======================================================================
116 // ************* STEERING parameters FOR ALICE SIMULATION **************
117 // --- Specify event type to be tracked through the ALICE setup
118 // --- All positions are in cm, angles in degrees, and P and E in GeV
121 gMC->SetProcess("DCAY",1);
122 gMC->SetProcess("PAIR",1);
123 gMC->SetProcess("COMP",1);
124 gMC->SetProcess("PHOT",1);
125 gMC->SetProcess("PFIS",0);
126 gMC->SetProcess("DRAY",0);
127 gMC->SetProcess("ANNI",1);
128 gMC->SetProcess("BREM",1);
129 gMC->SetProcess("MUNU",1);
130 gMC->SetProcess("CKOV",1);
131 gMC->SetProcess("HADR",0);
132 gMC->SetProcess("LOSS",2);
133 gMC->SetProcess("MULS",1);
134 gMC->SetProcess("RAYL",1);
136 Float_t cut = 1.e-3; // 1MeV cut by default
137 Float_t tofmax = 1.e10;
139 gMC->SetCut("CUTGAM", cut);
140 gMC->SetCut("CUTELE", cut);
141 gMC->SetCut("CUTNEU", cut);
142 gMC->SetCut("CUTHAD", cut);
143 gMC->SetCut("CUTMUO", cut);
144 gMC->SetCut("BCUTE", cut);
145 gMC->SetCut("BCUTM", cut);
146 gMC->SetCut("DCUTE", cut);
147 gMC->SetCut("DCUTM", cut);
148 gMC->SetCut("PPCUTM", cut);
149 gMC->SetCut("TOFMAX", tofmax);
151 // Special generation for Valgrind tests
152 // Each detector is fired by few particles selected
153 // to cover specific cases
156 // The cocktail itself
158 AliGenCocktail *gener = new AliGenCocktail();
159 Int_t nparticles = 20;
160 gener->SetPtRange(0.1,2.0001); //
161 gener->SetPhiRange(0, 360);
162 Float_t thmin = EtaToTheta(1.5); // theta min. <---> eta max
163 Float_t thmax = EtaToTheta(-1.5); // theta max. <---> eta min
164 gener->SetThetaRange(thmin,thmax);
165 gener->SetOrigin(0., 0., 0); //vertex position
166 gener->SetSigma(0., 0., 0); //Sigma in (X,Y,Z) (cm) on IP position
168 AliGenBox *gbox1 = new AliGenBox(nparticles/2);
169 gbox1->SetPart(kPiMinus);
170 gener->AddGenerator(gbox1,"GENBOX PIONS for ITS",1);
171 AliGenBox *gbox2 = new AliGenBox(nparticles/2);
172 gbox2->SetPart(kPiPlus);
173 gener->AddGenerator(gbox2,"GENBOX PIONS for ITS",1);
179 int nParticles = 50;//14022;
180 AliGenHIJINGpara *gener = new AliGenHIJINGpara(nParticles);
181 gener->SetMomentumRange(0.01, 5);
182 gener->SetPhiRange(0., 360.);
183 Float_t thmin = EtaToTheta(1.5); // theta min. <---> eta max
184 Float_t thmax = EtaToTheta(-1.5); // theta max. <---> eta min
185 gener->SetThetaRange(thmin,thmax);
186 gener->SetOrigin(0, 0, 0); //vertex position
187 gener->SetSigma(0, 0, 0); //Sigma in (X,Y,Z) (cm) on IP position
192 // Activate this line if you want the vertex smearing to happen
195 //VertexSmear_t perTrack;
196 //gener->SetVertexSmear(perTrack);
198 //AliMagF* field = new AliMagF("map","map",2, -1.,1., 15, AliMagF::k5kGUniform);
199 //TGeoGlobalMagField::Instance()->SetField(field);
200 TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", -1., -1., AliMagF::k5kG));
224 //=================== Alice BODY parameters =============================
225 AliBODY *BODY = new AliBODY("BODY", "Alice envelop");
229 //=================== MAG parameters ============================
230 // --- Start with Magnet since detector layouts may be depending ---
231 // --- on the selected Magnet dimensions ---
232 AliMAG *MAG = new AliMAG("MAG", "Magnet");
238 //=================== ABSO parameters ============================
239 AliABSO *ABSO = new AliABSOv3("ABSO", "Muon Absorber");
244 //=================== DIPO parameters ============================
246 AliDIPO *DIPO = new AliDIPOv3("DIPO", "Dipole version 3");
251 //=================== HALL parameters ============================
253 AliHALL *HALL = new AliHALLv3("HALL", "Alice Hall");
259 //=================== FRAME parameters ============================
261 AliFRAMEv2 *FRAME = new AliFRAMEv2("FRAME", "Space Frame");
267 //=================== SHIL parameters ============================
269 AliSHIL *SHIL = new AliSHILv3("SHIL", "Shielding Version 3");
275 //=================== PIPE parameters ============================
277 AliPIPE *PIPE = new AliPIPEv3("PIPE", "Beam Pipe");
282 //=================== ITS parameters ============================
284 // create segmentations:
285 AliITSUSegmentationPix* seg0 = new AliITSUSegmentationPix(0, // segID (0:9)
286 5, // chips per module
287 1750, // ncols (total for module)
289 20., // default row pitch in microns
290 20., // default col pitch in microns
291 300. // sensor thickness in microns
292 ); // see AliITSUSegmentationPix.h for extra options
293 seg0->Store(AliITSUGeomTGeo::GetITSsegmentationFileName());
294 AliITSUSegmentationPix* seg1 = new AliITSUSegmentationPix(1, // segID (0:9)
295 5, // chips per module
296 1750, // ncols (total for module)
298 20., // default row pitch in microns
299 20., // default col pitch in microns
300 300. // sensor thickness in microns
301 ); // see AliITSUSegmentationPix.h for extra options
302 seg1->Store(AliITSUGeomTGeo::GetITSsegmentationFileName());
304 int nmod,nlad; // modules per ladded, n ladders
305 const double kUM2CM=1e-4;
307 double thickLr = 0.15;
308 AliITSUv11 *ITS = new AliITSUv11("ITS Upgrade",7);
311 ITS->DefineLayerTurbo(0, 2.2, nmod*seg0->Dz()*kUM2CM, nlad, nmod, seg0->Dx()*kUM2CM, tilt, thickLr, seg0->Dy()*kUM2CM, seg0->GetDetTypeID());
314 ITS->DefineLayerTurbo(1, 2.8, nmod*seg0->Dz()*kUM2CM, nlad, nmod, seg0->Dx()*kUM2CM, tilt, thickLr, seg0->Dy()*kUM2CM, seg0->GetDetTypeID());
317 ITS->DefineLayerTurbo(2, 3.6, nmod*seg0->Dz()*kUM2CM, nlad, nmod, seg0->Dx()*kUM2CM, tilt, thickLr, seg0->Dy()*kUM2CM, seg0->GetDetTypeID());
320 ITS->DefineLayerTurbo(3, 20.0, nmod*seg0->Dz()*kUM2CM, nlad, nmod, seg0->Dx()*kUM2CM, tilt, thickLr, seg1->Dy()*kUM2CM, seg1->GetDetTypeID());
323 ITS->DefineLayerTurbo(4, 22.0, nmod*seg0->Dz()*kUM2CM, nlad, nmod, seg0->Dx()*kUM2CM, tilt, thickLr, seg1->Dy()*kUM2CM, seg1->GetDetTypeID());
326 ITS->DefineLayerTurbo(5, 40.0, nmod*seg0->Dz()*kUM2CM, nlad, nmod, seg0->Dx()*kUM2CM, tilt, thickLr, seg1->Dy()*kUM2CM, seg1->GetDetTypeID()); //41 creates ovl!
329 ITS->DefineLayerTurbo(6, 43.0, nmod*seg0->Dz()*kUM2CM, nlad, nmod, seg0->Dx()*kUM2CM, tilt, thickLr, seg1->Dy()*kUM2CM, seg1->GetDetTypeID());
331 ITS->AddBeamPipe(2.0, 2.08, 100);
337 //============================ TPC parameters ===================
338 AliTPC *TPC = new AliTPCv2("TPC", "Default");
343 //=================== TOF parameters ============================
344 AliTOF *TOF = new AliTOFv6T0("TOF", "normal TOF");
350 //=================== HMPID parameters ===========================
351 AliHMPID *HMPID = new AliHMPIDv3("HMPID", "normal HMPID");
358 //=================== ZDC parameters ============================
360 AliZDC *ZDC = new AliZDCv3("ZDC", "normal ZDC");
365 //=================== TRD parameters ============================
367 AliTRD *TRD = new AliTRDv1("TRD", "TRD slow simulator");
368 AliTRDgeometry *geoTRD = TRD->GetGeometry();
369 // Partial geometry: modules at 2,3,4,6,11,12,14,15
370 // starting at 6h in positive direction
371 geoTRD->SetSMstatus(0,0);
372 geoTRD->SetSMstatus(1,0);
373 geoTRD->SetSMstatus(5,0);
374 geoTRD->SetSMstatus(7,0);
375 geoTRD->SetSMstatus(8,0);
376 geoTRD->SetSMstatus(9,0);
377 geoTRD->SetSMstatus(10,0);
378 geoTRD->SetSMstatus(13,0);
379 geoTRD->SetSMstatus(16,0);
380 geoTRD->SetSMstatus(17,0);
385 //=================== FMD parameters ============================
386 AliFMD *FMD = new AliFMDv1("FMD", "normal FMD");
391 //=================== MUON parameters ===========================
392 // New MUONv1 version (geometry defined via builders)
393 AliMUON *MUON = new AliMUONv1("MUON","default");
395 //=================== PHOS parameters ===========================
399 AliPHOS *PHOS = new AliPHOSv1("PHOS", "IHEP");
405 //=================== PMD parameters ============================
406 AliPMD *PMD = new AliPMDv1("PMD", "normal PMD");
411 //=================== T0 parameters ============================
412 AliT0 *T0 = new AliT0v1("T0", "T0 Detector");
417 //=================== EMCAL parameters ============================
418 AliEMCAL *EMCAL = new AliEMCALv2("EMCAL", "EMCAL_COMPLETEV1");
423 //=================== ACORDE parameters ============================
424 AliACORDE *ACORDE = new AliACORDEv1("ACORDE", "normal ACORDE");
429 //=================== VZERO parameters ============================
430 AliVZERO *VZERO = new AliVZEROv7("VZERO", "normal VZERO");
436 Float_t EtaToTheta(Float_t arg){
437 return (180./TMath::Pi())*2.*atan(exp(-arg));