The present commit corresponds to an important change in the way the
[u/mrichter/AliRoot.git] / doc / aliroot-primer / scripts / Config.C
1 // Function converting pseudorapidity
2 // interval to polar angle interval. It is used to set 
3 // the limits in the generator
4 Float_t EtaToTheta(Float_t arg){
5   return (180./TMath::Pi())*2.*atan(exp(-arg));
6 }
7
8 // Set Random Number seed using the current time
9 TDatime dat;
10 static UInt_t sseed = dat.Get();
11
12 void Config()
13 {
14   gRandom->SetSeed(sseed);
15   cout<<"Seed for random number generation= "<<gRandom->GetSeed()<<endl; 
16     
17   // Load GEANT 3 library. It has to be in LD_LIBRARY_PATH
18   gSystem->Load("libgeant321");
19
20   // Instantiation of the particle transport package. gMC is set internaly
21   new TGeant3TGeo("C++ Interface to Geant3");
22
23   // Create run loader and set some properties
24   AliRunLoader* rl =  AliRunLoader::Open("galice.root",
25                                          AliConfig::GetDefaultEventFolderName(),
26                                          "recreate");
27   if (!rl) Fatal("Config.C","Can not instatiate the Run Loader");
28   rl->SetCompressionLevel(2);
29   rl->SetNumberOfEventsPerFile(3);
30
31   // Register the run loader in gAlice
32   gAlice->SetRunLoader(rl);
33
34   // Set external decayer
35   TVirtualMCDecayer *decayer = new AliDecayerPythia();
36   decayer->SetForceDecay(kAll); // kAll means no specific decay is forced
37   decayer->Init();
38
39   // Register the external decayer in the transport package
40   gMC->SetExternalDecayer(decayer);
41
42   // STEERING parameters FOR ALICE SIMULATION
43   // Specify event type to be transported through the ALICE setup
44   // All positions are in cm, angles in degrees, and P and E in GeV
45   // For the details see the GEANT 3 manual
46
47   // Switch on/off the physics processes (global)
48   // Please consult the file data/galice.cuts for detector
49   // specific settings, i.e. DRAY
50   gMC->SetProcess("DCAY",1); // Particle decay
51   gMC->SetProcess("PAIR",1); // Pair production
52   gMC->SetProcess("COMP",1); // Compton scattering
53   gMC->SetProcess("PHOT",1); // Photo effect
54   gMC->SetProcess("PFIS",0); // Photo fission
55   gMC->SetProcess("DRAY",0); // Delta rays
56   gMC->SetProcess("ANNI",1); // Positron annihilation
57   gMC->SetProcess("BREM",1); // Bremstrahlung
58   gMC->SetProcess("MUNU",1); // Muon nuclear interactions
59   gMC->SetProcess("CKOV",1); // Cerenkov production
60   gMC->SetProcess("HADR",1); // Hadronic interactions
61   gMC->SetProcess("LOSS",2); // Energy loss (2=complete fluct.)
62   gMC->SetProcess("MULS",1); // Multiple scattering
63   gMC->SetProcess("RAYL",1); // Rayleigh scattering
64   
65   // Set the transport package cuts
66   Float_t cut = 1.e-3;        // 1MeV cut by default
67   Float_t tofmax = 1.e10;
68
69   gMC->SetCut("CUTGAM", cut); // Cut for gammas
70   gMC->SetCut("CUTELE", cut); // Cut for electrons
71   gMC->SetCut("CUTNEU", cut); // Cut for neutral hadrons
72   gMC->SetCut("CUTHAD", cut); // Cut for charged hadrons
73   gMC->SetCut("CUTMUO", cut); // Cut for muons
74   gMC->SetCut("BCUTE",  cut); // Cut for electron brems.
75   gMC->SetCut("BCUTM",  cut); // Cut for muon brems.
76   gMC->SetCut("DCUTE",  cut); // Cut for electron delta-rays
77   gMC->SetCut("DCUTM",  cut); // Cut for muon delta-rays
78   gMC->SetCut("PPCUTM", cut); // Cut for e+e- pairs by muons
79   gMC->SetCut("TOFMAX", tofmax); // Time of flight cut
80   
81   // Set up the particle generation
82
83   // AliGenCocktail permits to combine several different generators
84   AliGenCocktail *gener = new AliGenCocktail();
85
86   // The phi range is always inside 0-360
87   gener->SetPhiRange(0, 360);
88
89   // Set pseudorapidity range from -8 to 8.
90   Float_t thmin = EtaToTheta(8);   // theta min. <---> eta max
91   Float_t thmax = EtaToTheta(-8);  // theta max. <---> eta min 
92   gener->SetThetaRange(thmin,thmax);
93
94   gener->SetOrigin(0, 0, 0);  // vertex position
95   gener->SetSigma(0, 0, 5.3);   // Sigma in (X,Y,Z) (cm) on IP position
96   gener->SetCutVertexZ(1.);     // Truncate at 1 sigma
97   gener->SetVertexSmear(kPerEvent); 
98
99   // First cocktail component: 100 ``background'' particles 
100   AliGenHIJINGpara *hijingparam = new AliGenHIJINGpara(100);
101   hijingparam->SetMomentumRange(0.2, 999);
102   gener->AddGenerator(hijingparam,"HIJING PARAM",1);
103
104   // Second cocktail component: one gamma in PHOS direction
105   AliGenBox *genbox = new AliGenBox(1);
106   genbox->SetMomentumRange(10,11.);
107   genbox->SetPhiRange(270.5,270.7);
108   genbox->SetThetaRange(90.5,90.7);
109   genbox->SetPart(kGamma);
110   gener->AddGenerator(genbox,"GENBOX GAMMA for PHOS",1);
111
112   gener->Init(); // Initialization of the coctail generator
113
114   // Field (the last parameter is 1 => L3 0.4 T)
115   TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", 2, 1., 1., 10., AliMagF::k5kG));
116
117   // Make sure the current ROOT directory is in galice.root 
118   rl->CdGAFile();
119
120   // Build the setup and set some detector parameters
121
122   // ALICE BODY parameters. BODY is always present
123   AliBODY *BODY = new AliBODY("BODY", "ALICE envelop");
124
125   // Start with Magnet since detector layouts may be depending
126   // on the selected Magnet dimensions
127   AliMAG *MAG = new AliMAG("MAG", "Magnet");
128
129   AliABSO *ABSO = new AliABSOv0("ABSO", "Muon Absorber");       // Absorber
130
131   AliDIPO *DIPO = new AliDIPOv2("DIPO", "Dipole version 2");    // Dipole magnet
132
133   AliHALL *HALL = new AliHALL("HALL", "ALICE Hall");            // Hall
134     
135   AliFRAMEv2 *FRAME = new AliFRAMEv2("FRAME", "Space Frame");   // Space frame
136
137   AliSHIL *SHIL = new AliSHILv2("SHIL", "Shielding Version 2"); // Shielding
138
139   AliPIPE *PIPE = new AliPIPEv0("PIPE", "Beam Pipe");           // Beam pipe
140     
141   // ITS parameters
142   AliITSvPPRasymmFMD *ITS  = new AliITSvPPRasymmFMD("ITS",
143                "ITS PPR detailed version with asymmetric services");
144   ITS->SetMinorVersion(2);      // don't change it if you're not an ITS developer
145   ITS->SetReadDet(kFALSE);      // don't change it if you're not an ITS developer
146   ITS->SetThicknessDet1(200.);  // detector thickness on layer 1:[100,300] mkm
147   ITS->SetThicknessDet2(200.);  // detector thickness on layer 2:[100,300] mkm
148   ITS->SetThicknessChip1(150.); // chip thickness on layer 1: [150,300] mkm
149   ITS->SetThicknessChip2(150.); // chip thickness on layer 2: [150,300]
150   ITS->SetRails(0);             // 1 --> rails in ; 0 --> rails out
151   ITS->SetCoolingFluid(1);      // 1 --> water ; 0 --> freon
152   ITS->SetEUCLID(0);            // no output for the EUCLID CAD system 
153
154   
155   AliTPC *TPC = new AliTPCv2("TPC", "Default");                 // TPC
156
157   AliTOF *TOF = new AliTOFv5T0("TOF", "normal TOF");            // TOF
158
159   AliHMPID *HMPID = new AliHMPIDv1("HMPID", "normal HMPID");         // HMPID
160
161   AliZDC *ZDC = new AliZDCv2("ZDC", "normal ZDC");              // ZDC
162
163   AliTRD *TRD = new AliTRDv1("TRD", "TRD slow simulator");      // TRD
164
165   AliFMD *FMD = new AliFMDv1("FMD", "normal FMD");              // FMD
166
167   AliMUON *MUON = new AliMUONv1("MUON", "default");             // MUON
168
169   AliPHOS *PHOS = new AliPHOSv1("PHOS", "IHEP");                // PHOS
170
171   AliPMD *PMD = new AliPMDv1("PMD", "normal PMD");              // PMD
172
173   AliT0 *T0 = new AliT0v1("T0", "T0 Detector");  // T0
174
175   // EMCAL
176   AliEMCAL *EMCAL = new AliEMCALv2("EMCAL", "EMCAL_COMPLETE");
177
178   AliVZERO *VZERO = new AliVZEROv7("VZERO", "normal VZERO");    // VZERO
179 }