dcf2b8912749f8bcfe3e8d17062553bd8179fb27
[u/mrichter/AliRoot.git] / TRD / AliTRDconfig.C
1 void Config()
2 {
3
4   Int_t iField = 0;
5
6   new AliGeant3("C++ Interface to Geant3");
7
8   // Create the output file
9   TFile *rootfile = new TFile("TRD_test.root","recreate");
10   rootfile->SetCompressionLevel(2);
11
12   // Define the monte carlo
13   TGeant3 *geant3 = (TGeant3*) gMC;
14
15   // Set external decayer
16   AliDecayer* decayer = new AliDecayerPythia();
17   decayer->SetForceDecay(kAll);
18   //decayer->SetForceDecay(kAll);
19   decayer->Init();
20   gMC->SetExternalDecayer(decayer);
21
22   // ******* GEANT STEERING parameters FOR ALICE SIMULATION *******
23   geant3->SetTRIG(1); // Number of events to be processed 
24   geant3->SetSWIT(4,10);
25   geant3->SetDEBU(0,0,1);
26   //geant3->SetSWIT(2,2);
27   geant3->SetDCAY(1);
28   geant3->SetPAIR(1);
29   geant3->SetCOMP(1);
30   geant3->SetPHOT(1);
31   geant3->SetPFIS(0);
32   geant3->SetDRAY(0);
33   geant3->SetANNI(1);
34   geant3->SetBREM(1);
35   geant3->SetMUNU(1);
36   geant3->SetCKOV(1);
37   geant3->SetHADR(1); // Select pure GEANH (HADR 1) or GEANH/NUCRIN (HADR 3)
38   geant3->SetLOSS(2);
39   geant3->SetMULS(1);
40   geant3->SetRAYL(1);
41   geant3->SetAUTO(1); // Select automatic STMIN etc... calc. (AUTO 1) or manual (AUTO 0)
42   geant3->SetABAN(0); // Restore 3.16 behaviour for abandoned tracks
43   geant3->SetOPTI(2); // Select optimisation level for GEANT geometry searches (0,1,2)
44   geant3->SetERAN(5.e-7);
45
46   Float_t cut    = 1.e-3; // 1MeV cut by default
47   Float_t tofmax = 1.e10;
48   //             GAM ELEC NHAD CHAD MUON EBREM MUHAB EDEL MUDEL MUPA TOFMAX
49   geant3->SetCUTS(cut,cut, cut, cut, cut, cut,  cut,  cut, cut,  cut, tofmax);
50
51   // ************* STEERING parameters FOR ALICE SIMULATION **************
52   // --- Specify event type to be tracked through the ALICE setup
53   // --- All positions are in cm, angles in degrees, and P and E in GeV
54
55   AliGenCocktail *gener = new AliGenCocktail();
56
57   AliGenBox *genEl = new AliGenBox(100);
58   genEl->SetOrigin(0,0,0);        // Vertex position
59   genEl->SetSigma(0,0,0);         // Sigma in (X,Y,Z) (cm) on IP position
60   genEl->SetPart(11);             // Only electrons 
61
62   AliGenBox *genPi = new AliGenBox(100);
63   genPi->SetOrigin(0,0,0);        // Vertex position
64   genPi->SetSigma(0,0,0);         // Sigma in (X,Y,Z) (cm) on IP position
65   genPi->SetPart(-211);           // Only pions 
66
67   gener->AddGenerator(genEl,"Electrons",1);
68   gener->AddGenerator(genPi,"Pions"    ,1);
69
70   if (iField) {
71
72     // With magnetic field on
73     AliGenerator *gg = gener->FirstGenerator()->Generator();
74     gg->SetMomentumRange(3.00,3.01);
75     gg->SetPhiRange(76.0,92.0);
76     gg->SetThetaRange(83.0,97.0);
77     gg = gener->NextGenerator()->Generator();
78     gg->SetMomentumRange(0.560,0.561);
79     gg->SetPhiRange(62.0,78.0);
80     gg->SetThetaRange(83.0,97.0);
81
82     gener->Init();
83
84     // Specify maximum magnetic field in Tesla (neg. ==> default field)
85     // 0.4 T
86     gAlice->SetField(-999,2,2.0);    
87
88   }
89   else {
90
91     // With magnetic field off
92     AliGenerator *gg = gener->FirstGenerator()->Generator();
93     gg->SetMomentumRange(3.00,3.01);
94     gg->SetPhiRange(82.0,98.0);
95     gg->SetThetaRange(83.0,97.0);
96     gg = gener->NextGenerator()->Generator();
97     gg->SetMomentumRange(0.560,0.561);
98     gg->SetPhiRange(82.0,98.0);
99     gg->SetThetaRange(83.0,97.0);
100
101     gener->Init();
102
103     // Specify maximum magnetic field in Tesla (neg. ==> default field)
104     // No field
105     gAlice->SetField(0);    
106
107   }
108
109   Int_t iMAG   = 1;
110   Int_t iITS   = 1;
111   Int_t iTPC   = 1;
112   Int_t iTRD   = 1;
113   Int_t iABSO  = 1;
114   Int_t iDIPO  = 1;
115   Int_t iHALL  = 1;
116   Int_t iFRAME = 1;
117   Int_t iSHIL  = 1;
118   Int_t iPIPE  = 1;
119
120   //=================== Alice BODY parameters =============================
121   AliBODY *BODY = new AliBODY("BODY","Alice envelop");
122
123   if (iMAG) {
124     //=================== MAG parameters ============================
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
130   if (iABSO) { 
131     //=================== ABSO parameters ============================
132     AliABSO *ABSO  = new AliABSOv0("ABSO","Muon Absorber");
133   }
134
135   if (iDIPO) {
136     //=================== DIPO parameters ============================
137     AliDIPO *DIPO  = new AliDIPOv2("DIPO","Dipole version 2");
138   }
139
140   if (iHALL) {
141     //=================== HALL parameters ============================
142     AliHALL *HALL  = new AliHALL("HALL","Alice Hall");
143   }
144
145   if (iFRAME) {
146     //=================== FRAME parameters ============================
147     AliFRAME *FRAME  = new AliFRAMEv1("FRAME","Space Frame");
148   }
149
150   if (iSHIL) {
151     //=================== SHIL parameters ============================
152     AliSHIL *SHIL  = new AliSHILv0("SHIL","Shielding");
153   }
154
155   if (iPIPE) {
156     //=================== PIPE parameters ============================
157     AliPIPE *PIPE  = new AliPIPEv0("PIPE","Beam Pipe");
158   }
159
160   if (iITS) {
161     //=================== ITS parameters ============================
162     //
163     // EUCLID is a flag to output (=1) both geometry and media to two ASCII files 
164     // (called by default ITSgeometry.euc and ITSgeometry.tme) in a format
165     // understandable to the CAD system EUCLID. The default (=0) means that you 
166     // dont want to use this facility.
167     //
168
169     AliITS *ITS  = new AliITSv5asymm("ITS","Updates ITS TDR detailed version with asymmetric services");
170     ITS->SetEUCLID(0);
171
172   }
173
174   if (iTPC) {
175     //============================ TPC parameters ================================
176     // --- This allows the user to specify sectors for the SLOW (TPC geometry 2)
177     // --- Simulator. SecAL (SecAU) <0 means that ALL lower (upper)
178     // --- sectors are specified, any value other than that requires at least one 
179     // --- sector (lower or upper)to be specified!
180     // --- Reminder: sectors 1-24 are lower sectors (1-12 -> z>0, 13-24 -> z<0)
181     // ---           sectors 25-72 are the upper ones (25-48 -> z>0, 49-72 -> z<0)
182     // --- SecLows - number of lower sectors specified (up to 6)
183     // --- SecUps - number of upper sectors specified (up to 12)
184     // --- Sens - sensitive strips for the Slow Simulator !!!
185     // --- This does NOT work if all S or L-sectors are specified, i.e.
186     // --- if SecAL or SecAU < 0
187     //
188     //
189     //-----------------------------------------------------------------------------
190
191     AliTPC *TPC  = new AliTPCv2("TPC","Default");
192     // All sectors included 
193     TPC->SetSecAL(-1);
194     TPC->SetSecAU(-1);
195
196   }
197
198   if (iTRD) {
199     //=================== TRD parameters ============================
200   
201     AliTRD *TRD  = new AliTRDv1("TRD","TRD slow simulator");
202     TRD->SetSensChamber(2);
203     TRD->SetSensSector(13);
204   
205     // Set to detailed display
206     TRD->SetDisplayType(1);
207
208     // Draw TR photons
209     TRD->SetDrawTR(1);
210
211     // Switch on TR
212     AliTRDsim *TRDsim = TRD->CreateTR();
213
214   }
215         
216 }