Changes required by the coding conventions
[u/mrichter/AliRoot.git] / MONITOR / galice.C
1 void galice()
2 {
3     AliRunLoader* runLoader = AliRunLoader::Open("galice.root", AliConfig::fgkDefaultEventFolderName, "recreate");
4     runLoader->SetCompressionLevel(2);
5     runLoader->SetNumberOfEventsPerFile(1);
6     gAlice->SetRunLoader(runLoader);
7
8     gSystem->Load("libgeant321");
9     new TGeant3("C++ Interface to Geant3");
10
11     // Field (L3 0.4 T)
12     AliMagFMaps* field = new AliMagFMaps("Maps","Maps", 2, 1., 10., AliMagFMaps::k4kG);
13     gAlice->SetField(field);    
14
15     Int_t   iABSO  =  0;
16     Int_t   iDIPO  =  0;
17     Int_t   iFMD   =  0;
18     Int_t   iFRAME =  0;
19     Int_t   iHALL  =  0;
20     Int_t   iITS   =  1;
21     Int_t   iMAG   =  0;
22     Int_t   iMUON  =  0;
23     Int_t   iPHOS  =  0;
24     Int_t   iPIPE  =  0;
25     Int_t   iPMD   =  0;
26     Int_t   iRICH  =  0;
27     Int_t   iSHIL  =  0;
28     Int_t   iSTART =  0;
29     Int_t   iTOF   =  0;
30     Int_t   iTPC   =  1;
31     Int_t   iTRD   =  0;
32     Int_t   iZDC   =  0;
33     Int_t   iEMCAL =  0;
34     Int_t   iCRT   =  0;
35     Int_t   iVZERO =  0;
36     runLoader->CdGAFile();
37     //=================== Alice BODY parameters =============================
38     AliBODY *BODY = new AliBODY("BODY", "Alice envelop");
39
40     if (iMAG)
41     {
42         //=================== MAG parameters ============================
43         // --- Start with Magnet since detector layouts may be depending ---
44         // --- on the selected Magnet dimensions ---
45         AliMAG *MAG = new AliMAG("MAG", "Magnet");
46     }
47
48
49     if (iABSO)
50     {
51         //=================== ABSO parameters ============================
52         AliABSO *ABSO = new AliABSOv0("ABSO", "Muon Absorber");
53     }
54
55     if (iDIPO)
56     {
57         //=================== DIPO parameters ============================
58
59         AliDIPO *DIPO = new AliDIPOv2("DIPO", "Dipole version 2");
60     }
61
62     if (iHALL)
63     {
64         //=================== HALL parameters ============================
65
66         AliHALL *HALL = new AliHALL("HALL", "Alice Hall");
67     }
68
69
70     if (iFRAME)
71     {
72         //=================== FRAME parameters ============================
73
74         AliFRAMEv2 *FRAME = new AliFRAMEv2("FRAME", "Space Frame");
75         if (geo == kHoles) {
76             FRAME->SetHoles(1);
77         } else {
78             FRAME->SetHoles(0);
79         }
80     }
81
82     if (iSHIL)
83     {
84         //=================== SHIL parameters ============================
85
86         AliSHIL *SHIL = new AliSHILv2("SHIL", "Shielding Version 2");
87     }
88
89
90     if (iPIPE)
91     {
92         //=================== PIPE parameters ============================
93
94         AliPIPE *PIPE = new AliPIPEv0("PIPE", "Beam Pipe");
95     }
96  
97     if(iITS) {
98
99     //=================== ITS parameters ============================
100     //
101     // As the innermost detector in ALICE, the Inner Tracking System "impacts" on
102     // almost all other detectors. This involves the fact that the ITS geometry
103     // still has several options to be followed in parallel in order to determine
104     // the best set-up which minimizes the induced background. All the geometries
105     // available to date are described in the following. Read carefully the comments
106     // and use the default version (the only one uncommented) unless you are making
107     // comparisons and you know what you are doing. In this case just uncomment the
108     // ITS geometry you want to use and run Aliroot.
109     //
110     // Detailed geometries:         
111     //
112     //
113     //AliITS *ITS  = new AliITSv5symm("ITS","Updated ITS TDR detailed version with symmetric services");
114     //
115     //AliITS *ITS  = new AliITSv5asymm("ITS","Updates ITS TDR detailed version with asymmetric services");
116     //
117         AliITSvPPRasymm *ITS  = new AliITSvPPRasymm("ITS","New ITS PPR detailed version with asymmetric services");
118         ITS->SetMinorVersion(2);                                         // don't touch this parameter if you're not an ITS developer
119         ITS->SetReadDet(kFALSE);                                         // don't touch this parameter if you're not an ITS developer
120     //    ITS->SetWriteDet("$ALICE_ROOT/ITS/ITSgeometry_vPPRasymm2.det");  // don't touch this parameter if you're not an ITS developer
121         ITS->SetThicknessDet1(200.);   // detector thickness on layer 1 must be in the range [100,300]
122         ITS->SetThicknessDet2(200.);   // detector thickness on layer 2 must be in the range [100,300]
123         ITS->SetThicknessChip1(200.);  // chip thickness on layer 1 must be in the range [150,300]
124         ITS->SetThicknessChip2(200.);  // chip thickness on layer 2 must be in the range [150,300]
125         ITS->SetRails(0);            // 1 --> rails in ; 0 --> rails out
126         ITS->SetCoolingFluid(1);   // 1 --> water ; 0 --> freon
127         //
128     //AliITSvPPRsymm *ITS  = new AliITSvPPRsymm("ITS","New ITS PPR detailed version with symmetric services");
129     //ITS->SetMinorVersion(2);                                       // don't touch this parameter if you're not an ITS developer
130     //ITS->SetReadDet(kFALSE);                                       // don't touch this parameter if you're not an ITS developer
131     //ITS->SetWriteDet("$ALICE_ROOT/ITS/ITSgeometry_vPPRsymm2.det"); // don't touch this parameter if you're not an ITS developer
132     //ITS->SetThicknessDet1(200.);   // detector thickness on layer 1 must be in the range [100,300]
133     //ITS->SetThicknessDet2(200.);   // detector thickness on layer 2 must be in the range [100,300]
134     //ITS->SetThicknessChip1(200.);  // chip thickness on layer 1 must be in the range [150,300]
135     //ITS->SetThicknessChip2(200.);  // chip thickness on layer 2 must be in the range [150,300]
136     //ITS->SetRails(0);              // 1 --> rails in ; 0 --> rails out
137     //ITS->SetCoolingFluid(1);       // 1 --> water ; 0 --> freon
138     //
139     //
140     // Coarse geometries (warning: no hits are produced with these coarse geometries and they unuseful 
141     // for reconstruction !):
142     //                                                     
143     //
144     //AliITSvPPRcoarseasymm *ITS  = new AliITSvPPRcoarseasymm("ITS","New ITS PPR coarse version with asymmetric services");
145     //ITS->SetRails(0);                // 1 --> rails in ; 0 --> rails out
146     //ITS->SetSupportMaterial(0);      // 0 --> Copper ; 1 --> Aluminum ; 2 --> Carbon
147     //
148     //AliITS *ITS  = new AliITSvPPRcoarsesymm("ITS","New ITS PPR coarse version with symmetric services");
149     //ITS->SetRails(0);                // 1 --> rails in ; 0 --> rails out
150     //ITS->SetSupportMaterial(0);      // 0 --> Copper ; 1 --> Aluminum ; 2 --> Carbon
151     //                      
152     //
153     //
154     // Geant3 <-> EUCLID conversion
155     // ============================
156     //
157     // SetEUCLID is a flag to output (=1) or not to output (=0) both geometry and
158     // media to two ASCII files (called by default ITSgeometry.euc and
159     // ITSgeometry.tme) in a format understandable to the CAD system EUCLID.
160     // The default (=0) means that you dont want to use this facility.
161     //
162      ITS->SetEUCLID(0);  
163     }
164
165     if (iTPC)
166     {
167         //============================ TPC parameters ================================
168         // --- This allows the user to specify sectors for the SLOW (TPC geometry 2)
169         // --- Simulator. SecAL (SecAU) <0 means that ALL lower (upper)
170         // --- sectors are specified, any value other than that requires at least one 
171         // --- sector (lower or upper)to be specified!
172         // --- Reminder: sectors 1-24 are lower sectors (1-12 -> z>0, 13-24 -> z<0)
173         // ---           sectors 25-72 are the upper ones (25-48 -> z>0, 49-72 -> z<0)
174         // --- SecLows - number of lower sectors specified (up to 6)
175         // --- SecUps - number of upper sectors specified (up to 12)
176         // --- Sens - sensitive strips for the Slow Simulator !!!
177         // --- This does NOT work if all S or L-sectors are specified, i.e.
178         // --- if SecAL or SecAU < 0
179         //
180         //
181         //-----------------------------------------------------------------------------
182
183         //  gROOT->LoadMacro("SetTPCParam.C");
184         //  AliTPCParam *param = SetTPCParam();
185         AliTPC *TPC = new AliTPCv2("TPC", "Default");
186
187         // All sectors included 
188         TPC->SetSecAU(-1);
189         TPC->SetSecAL(-1);
190     }
191
192
193     if (iTOF) {
194         if (geo == kHoles) {
195         //=================== TOF parameters ============================
196             AliTOF *TOF = new AliTOFv2FHoles("TOF", "TOF with Holes");
197         } else {
198             AliTOF *TOF = new AliTOFv4T0("TOF", "normal TOF");
199         }
200     }
201
202
203     if (iRICH)
204     {
205         //=================== RICH parameters ===========================
206         AliRICH *RICH = new AliRICHv3("RICH", "normal RICH");
207
208     }
209
210
211     if (iZDC)
212     {
213         //=================== ZDC parameters ============================
214
215         AliZDC *ZDC = new AliZDCv2("ZDC", "normal ZDC");
216     }
217
218     if (iTRD)
219     {
220         //=================== TRD parameters ============================
221
222         AliTRD *TRD = new AliTRDv1("TRD", "TRD slow simulator");
223
224         // Select the gas mixture (0: 97% Xe + 3% isobutane, 1: 90% Xe + 10% CO2)
225         TRD->SetGasMix(1);
226         if (geo == kHoles) {
227             // With hole in front of PHOS
228             TRD->SetPHOShole();
229             // With hole in front of RICH
230             TRD->SetRICHhole();
231         }
232             // Switch on TR
233             AliTRDsim *TRDsim = TRD->CreateTR();
234     }
235
236     if (iFMD)
237     {
238         //=================== FMD parameters ============================
239         AliFMD *FMD = new AliFMDv1("FMD", "normal FMD");
240         FMD->SetRingsSi1(256);
241         FMD->SetRingsSi2(128);
242         FMD->SetSectorsSi1(20);
243         FMD->SetSectorsSi2(40);      
244    }
245
246     if (iMUON)
247     {
248         //=================== MUON parameters ===========================
249
250         AliMUON *MUON = new AliMUONv1("MUON", "default");
251     }
252     //=================== PHOS parameters ===========================
253
254     if (iPHOS)
255     {
256         AliPHOS *PHOS = new AliPHOSv1("PHOS", "IHEP");
257     }
258
259
260     if (iPMD)
261     {
262         //=================== PMD parameters ============================
263         AliPMD *PMD = new AliPMDv1("PMD", "normal PMD");
264     }
265
266     if (iSTART)
267     {
268         //=================== START parameters ============================
269         AliSTART *START = new AliSTARTv1("START", "START Detector");
270     }
271
272     if (iEMCAL)
273     {
274         //=================== EMCAL parameters ============================
275         AliEMCAL *EMCAL = new AliEMCALv1("EMCAL", "EMCALArch1a");
276     }
277
278      if (iCRT)
279     {
280         //=================== CRT parameters ============================
281         AliCRT *CRT = new AliCRTv0("CRT", "normal ACORDE");
282     }
283
284      if (iVZERO)
285     {
286         //=================== CRT parameters ============================
287         AliVZERO *VZERO = new AliVZEROv2("VZERO", "normal VZERO");
288     }
289
290      AliConfig::Instance()->Add(gMC);
291      gMC->Init();
292      gAlice->InitLoaders();
293      gAlice->Write();
294      runLoader->MakeTree("E");
295      runLoader->Write();
296      delete runLoader;
297 }