a9e2aefa |
1 | /************************************************************************** |
2 | * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * |
3 | * * |
4 | * Author: The ALICE Off-line Project. * |
5 | * Contributors are mentioned in the code where appropriate. * |
6 | * * |
7 | * Permission to use, copy, modify and distribute this software and its * |
8 | * documentation strictly for non-commercial purposes is hereby granted * |
9 | * without fee, provided that the above copyright notice appears in all * |
10 | * copies and that both the copyright notice and this permission notice * |
11 | * appear in the supporting documentation. The authors make no claims * |
12 | * about the suitability of this software for any purpose. It is * |
13 | * provided "as is" without express or implied warranty. * |
14 | **************************************************************************/ |
15 | |
16 | /* |
17 | $Log$ |
18 | Revision 1.1.2.14 2000/06/14 14:37:25 morsch |
19 | Initialization of TriggerCircuit added (PC) |
20 | |
21 | Revision 1.1.2.13 2000/06/09 21:55:47 morsch |
22 | Most coding rule violations corrected. |
23 | |
24 | Revision 1.1.2.12 2000/05/05 11:34:29 morsch |
25 | Log inside comments. |
26 | |
27 | Revision 1.1.2.11 2000/05/05 10:06:48 morsch |
28 | Coding Rule violations regarding trigger section corrected (CP) |
29 | Log messages included. |
30 | */ |
31 | |
32 | ///////////////////////////////////////////////////////// |
33 | // Manager and hits classes for set:MUON version 0 // |
34 | ///////////////////////////////////////////////////////// |
35 | |
36 | #include <TTUBE.h> |
37 | #include <TNode.h> |
38 | #include <TRandom.h> |
39 | #include <TLorentzVector.h> |
40 | #include <iostream.h> |
41 | |
42 | #include "AliMUONv1.h" |
43 | #include "AliRun.h" |
44 | #include "AliMC.h" |
45 | #include "AliCallf77.h" |
46 | #include "AliConst.h" |
47 | #include "AliMUONChamber.h" |
48 | #include "AliMUONHit.h" |
49 | #include "AliMUONPadHit.h" |
50 | |
51 | ClassImp(AliMUONv1) |
52 | |
53 | //___________________________________________ |
54 | AliMUONv1::AliMUONv1() : AliMUON() |
55 | { |
56 | // Constructor |
57 | fChambers = 0; |
58 | } |
59 | |
60 | //___________________________________________ |
61 | AliMUONv1::AliMUONv1(const char *name, const char *title) |
62 | : AliMUON(name,title) |
63 | { |
64 | // Constructor |
65 | } |
66 | |
67 | //___________________________________________ |
68 | void AliMUONv1::CreateGeometry() |
69 | { |
70 | // |
71 | // Note: all chambers have the same structure, which could be |
72 | // easily parameterised. This was intentionally not done in order |
73 | // to give a starting point for the implementation of the actual |
74 | // design of each station. |
75 | Int_t *idtmed = fIdtmed->GetArray()-1099; |
76 | |
77 | // Distance between Stations |
78 | // |
79 | Float_t bpar[3]; |
80 | Float_t tpar[3]; |
81 | Float_t pgpar[10]; |
82 | Float_t zpos1, zpos2, zfpos; |
83 | Float_t dframep=.001; // Value for station 3 should be 6 ... |
84 | Float_t dframep1=.001; |
85 | // Bool_t frames=kTRUE; |
86 | Bool_t frames=kFALSE; |
87 | |
88 | Float_t dframez=0.9; |
89 | Float_t dr; |
90 | Float_t dstation; |
91 | |
92 | // |
93 | // Rotation matrices in the x-y plane |
94 | Int_t idrotm[1199]; |
95 | // phi= 0 deg |
96 | AliMatrix(idrotm[1100], 90., 0., 90., 90., 0., 0.); |
97 | // phi= 90 deg |
98 | AliMatrix(idrotm[1101], 90., 90., 90., 180., 0., 0.); |
99 | // phi= 180 deg |
100 | AliMatrix(idrotm[1102], 90., 180., 90., 270., 0., 0.); |
101 | // phi= 270 deg |
102 | AliMatrix(idrotm[1103], 90., 270., 90., 0., 0., 0.); |
103 | // |
104 | Float_t phi=2*TMath::Pi()/12/2; |
105 | |
106 | // |
107 | // pointer to the current chamber |
108 | // pointer to the current chamber |
109 | Int_t idAlu1=idtmed[1103]; |
110 | Int_t idAlu2=idtmed[1104]; |
111 | // Int_t idAlu1=idtmed[1100]; |
112 | // Int_t idAlu2=idtmed[1100]; |
113 | Int_t idAir=idtmed[1100]; |
114 | Int_t idGas=idtmed[1105]; |
115 | |
116 | |
117 | AliMUONChamber *iChamber, *iChamber1, *iChamber2; |
118 | //******************************************************************** |
119 | // Station 1 ** |
120 | //******************************************************************** |
121 | // CONCENTRIC |
122 | // indices 1 and 2 for first and second chambers in the station |
123 | // iChamber (first chamber) kept for other quanties than Z, |
124 | // assumed to be the same in both chambers |
125 | iChamber1 = iChamber = (AliMUONChamber*) (*fChambers)[0]; |
126 | iChamber2 =(AliMUONChamber*) (*fChambers)[1]; |
127 | zpos1=iChamber1->Z(); |
128 | zpos2=iChamber2->Z(); |
129 | dstation = zpos2 - zpos1; |
130 | zfpos=-(iChamber->DGas()+dframez+iChamber->DAlu())/2; |
131 | |
132 | // |
133 | // Mother volume |
134 | tpar[0] = iChamber->RInner()-dframep1; |
135 | tpar[1] = (iChamber->ROuter()+dframep1)/TMath::Cos(phi); |
136 | tpar[2] = dstation/4; |
137 | |
138 | gMC->Gsvolu("C01M", "TUBE", idAir, tpar, 3); |
139 | gMC->Gsvolu("C02M", "TUBE", idAir, tpar, 3); |
140 | gMC->Gspos("C01M", 1, "ALIC", 0., 0., zpos1 , 0, "ONLY"); |
141 | gMC->Gspos("C02M", 1, "ALIC", 0., 0., zpos2 , 0, "ONLY"); |
142 | // Aluminium frames |
143 | // Outer frames |
144 | pgpar[0] = 360/12/2; |
145 | pgpar[1] = 360.; |
146 | pgpar[2] = 12.; |
147 | pgpar[3] = 2; |
148 | pgpar[4] = -dframez/2; |
149 | pgpar[5] = iChamber->ROuter(); |
150 | pgpar[6] = pgpar[5]+dframep1; |
151 | pgpar[7] = +dframez/2; |
152 | pgpar[8] = pgpar[5]; |
153 | pgpar[9] = pgpar[6]; |
154 | gMC->Gsvolu("C01O", "PGON", idAlu1, pgpar, 10); |
155 | gMC->Gsvolu("C02O", "PGON", idAlu1, pgpar, 10); |
156 | gMC->Gspos("C01O",1,"C01M", 0.,0.,-zfpos, 0,"ONLY"); |
157 | gMC->Gspos("C01O",2,"C01M", 0.,0.,+zfpos, 0,"ONLY"); |
158 | gMC->Gspos("C02O",1,"C02M", 0.,0.,-zfpos, 0,"ONLY"); |
159 | gMC->Gspos("C02O",2,"C02M", 0.,0.,+zfpos, 0,"ONLY"); |
160 | // |
161 | // Inner frame |
162 | tpar[0]= iChamber->RInner()-dframep1; |
163 | tpar[1]= iChamber->RInner(); |
164 | tpar[2]= dframez/2; |
165 | gMC->Gsvolu("C01I", "TUBE", idAlu1, tpar, 3); |
166 | gMC->Gsvolu("C02I", "TUBE", idAlu1, tpar, 3); |
167 | |
168 | gMC->Gspos("C01I",1,"C01M", 0.,0.,-zfpos, 0,"ONLY"); |
169 | gMC->Gspos("C01I",2,"C01M", 0.,0.,+zfpos, 0,"ONLY"); |
170 | gMC->Gspos("C02I",1,"C02M", 0.,0.,-zfpos, 0,"ONLY"); |
171 | gMC->Gspos("C02I",2,"C02M", 0.,0.,+zfpos, 0,"ONLY"); |
172 | // |
173 | // Frame Crosses |
174 | if (frames) { |
175 | |
176 | bpar[0] = (iChamber->ROuter() - iChamber->RInner())/2; |
177 | bpar[1] = dframep1/2; |
178 | bpar[2] = dframez/2; |
179 | gMC->Gsvolu("C01B", "BOX", idAlu1, bpar, 3); |
180 | gMC->Gsvolu("C02B", "BOX", idAlu1, bpar, 3); |
181 | |
182 | gMC->Gspos("C01B",1,"C01M", +iChamber->RInner()+bpar[0] , 0,-zfpos, |
183 | idrotm[1100],"ONLY"); |
184 | gMC->Gspos("C01B",2,"C01M", -iChamber->RInner()-bpar[0] , 0,-zfpos, |
185 | idrotm[1100],"ONLY"); |
186 | gMC->Gspos("C01B",3,"C01M", 0, +iChamber->RInner()+bpar[0] ,-zfpos, |
187 | idrotm[1101],"ONLY"); |
188 | gMC->Gspos("C01B",4,"C01M", 0, -iChamber->RInner()-bpar[0] ,-zfpos, |
189 | idrotm[1101],"ONLY"); |
190 | gMC->Gspos("C01B",5,"C01M", +iChamber->RInner()+bpar[0] , 0,+zfpos, |
191 | idrotm[1100],"ONLY"); |
192 | gMC->Gspos("C01B",6,"C01M", -iChamber->RInner()-bpar[0] , 0,+zfpos, |
193 | idrotm[1100],"ONLY"); |
194 | gMC->Gspos("C01B",7,"C01M", 0, +iChamber->RInner()+bpar[0] ,+zfpos, |
195 | idrotm[1101],"ONLY"); |
196 | gMC->Gspos("C01B",8,"C01M", 0, -iChamber->RInner()-bpar[0] ,+zfpos, |
197 | idrotm[1101],"ONLY"); |
198 | |
199 | gMC->Gspos("C02B",1,"C02M", +iChamber->RInner()+bpar[0] , 0,-zfpos, |
200 | idrotm[1100],"ONLY"); |
201 | gMC->Gspos("C02B",2,"C02M", -iChamber->RInner()-bpar[0] , 0,-zfpos, |
202 | idrotm[1100],"ONLY"); |
203 | gMC->Gspos("C02B",3,"C02M", 0, +iChamber->RInner()+bpar[0] ,-zfpos, |
204 | idrotm[1101],"ONLY"); |
205 | gMC->Gspos("C02B",4,"C02M", 0, -iChamber->RInner()-bpar[0] ,-zfpos, |
206 | idrotm[1101],"ONLY"); |
207 | gMC->Gspos("C02B",5,"C02M", +iChamber->RInner()+bpar[0] , 0,+zfpos, |
208 | idrotm[1100],"ONLY"); |
209 | gMC->Gspos("C02B",6,"C02M", -iChamber->RInner()-bpar[0] , 0,+zfpos, |
210 | idrotm[1100],"ONLY"); |
211 | gMC->Gspos("C02B",7,"C02M", 0, +iChamber->RInner()+bpar[0] ,+zfpos, |
212 | idrotm[1101],"ONLY"); |
213 | gMC->Gspos("C02B",8,"C02M", 0, -iChamber->RInner()-bpar[0] ,+zfpos, |
214 | idrotm[1101],"ONLY"); |
215 | } |
216 | // |
217 | // Chamber Material represented by Alu sheet |
218 | tpar[0]= iChamber->RInner(); |
219 | tpar[1]= iChamber->ROuter(); |
220 | tpar[2] = (iChamber->DGas()+iChamber->DAlu())/2; |
221 | gMC->Gsvolu("C01A", "TUBE", idAlu2, tpar, 3); |
222 | gMC->Gsvolu("C02A", "TUBE",idAlu2, tpar, 3); |
223 | gMC->Gspos("C01A", 1, "C01M", 0., 0., 0., 0, "ONLY"); |
224 | gMC->Gspos("C02A", 1, "C02M", 0., 0., 0., 0, "ONLY"); |
225 | // |
226 | // Sensitive volumes |
227 | // tpar[2] = iChamber->DGas(); |
228 | tpar[2] = iChamber->DGas()/2; |
229 | gMC->Gsvolu("C01G", "TUBE", idtmed[1108], tpar, 3); |
230 | gMC->Gsvolu("C02G", "TUBE", idtmed[1108], tpar, 3); |
231 | gMC->Gspos("C01G", 1, "C01A", 0., 0., 0., 0, "ONLY"); |
232 | gMC->Gspos("C02G", 1, "C02A", 0., 0., 0., 0, "ONLY"); |
233 | // |
234 | // Frame Crosses to be placed inside gas |
235 | if (frames) { |
236 | |
237 | dr = (iChamber->ROuter() - iChamber->RInner()); |
238 | bpar[0] = TMath::Sqrt(dr*dr-dframep1*dframep1/4)/2; |
239 | bpar[1] = dframep1/2; |
240 | bpar[2] = iChamber->DGas()/2; |
241 | gMC->Gsvolu("C01F", "BOX", idAlu1, bpar, 3); |
242 | gMC->Gsvolu("C02F", "BOX", idAlu1, bpar, 3); |
243 | |
244 | gMC->Gspos("C01F",1,"C01G", +iChamber->RInner()+bpar[0] , 0, 0, |
245 | idrotm[1100],"ONLY"); |
246 | gMC->Gspos("C01F",2,"C01G", -iChamber->RInner()-bpar[0] , 0, 0, |
247 | idrotm[1100],"ONLY"); |
248 | gMC->Gspos("C01F",3,"C01G", 0, +iChamber->RInner()+bpar[0] , 0, |
249 | idrotm[1101],"ONLY"); |
250 | gMC->Gspos("C01F",4,"C01G", 0, -iChamber->RInner()-bpar[0] , 0, |
251 | idrotm[1101],"ONLY"); |
252 | |
253 | gMC->Gspos("C02F",1,"C02G", +iChamber->RInner()+bpar[0] , 0, 0, |
254 | idrotm[1100],"ONLY"); |
255 | gMC->Gspos("C02F",2,"C02G", -iChamber->RInner()-bpar[0] , 0, 0, |
256 | idrotm[1100],"ONLY"); |
257 | gMC->Gspos("C02F",3,"C02G", 0, +iChamber->RInner()+bpar[0] , 0, |
258 | idrotm[1101],"ONLY"); |
259 | gMC->Gspos("C02F",4,"C02G", 0, -iChamber->RInner()-bpar[0] , 0, |
260 | idrotm[1101],"ONLY"); |
261 | } |
262 | |
263 | // |
264 | // |
265 | //******************************************************************** |
266 | // Station 2 ** |
267 | //******************************************************************** |
268 | // indices 1 and 2 for first and second chambers in the station |
269 | // iChamber (first chamber) kept for other quanties than Z, |
270 | // assumed to be the same in both chambers |
271 | iChamber1 = iChamber = (AliMUONChamber*) (*fChambers)[2]; |
272 | iChamber2 =(AliMUONChamber*) (*fChambers)[3]; |
273 | zpos1=iChamber1->Z(); |
274 | zpos2=iChamber2->Z(); |
275 | dstation = zpos2 - zpos1; |
276 | zfpos=-(iChamber->DGas()+dframez+iChamber->DAlu())/2; |
277 | |
278 | // |
279 | // Mother volume |
280 | tpar[0] = iChamber->RInner()-dframep; |
281 | tpar[1] = (iChamber->ROuter()+dframep)/TMath::Cos(phi); |
282 | tpar[2] = dstation/4; |
283 | |
284 | gMC->Gsvolu("C03M", "TUBE", idAir, tpar, 3); |
285 | gMC->Gsvolu("C04M", "TUBE", idAir, tpar, 3); |
286 | gMC->Gspos("C03M", 1, "ALIC", 0., 0., zpos1 , 0, "ONLY"); |
287 | gMC->Gspos("C04M", 1, "ALIC", 0., 0., zpos2 , 0, "ONLY"); |
288 | // Aluminium frames |
289 | // Outer frames |
290 | pgpar[0] = 360/12/2; |
291 | pgpar[1] = 360.; |
292 | pgpar[2] = 12.; |
293 | pgpar[3] = 2; |
294 | pgpar[4] = -dframez/2; |
295 | pgpar[5] = iChamber->ROuter(); |
296 | pgpar[6] = pgpar[5]+dframep; |
297 | pgpar[7] = +dframez/2; |
298 | pgpar[8] = pgpar[5]; |
299 | pgpar[9] = pgpar[6]; |
300 | gMC->Gsvolu("C03O", "PGON", idAlu1, pgpar, 10); |
301 | gMC->Gsvolu("C04O", "PGON", idAlu1, pgpar, 10); |
302 | gMC->Gspos("C03O",1,"C03M", 0.,0.,-zfpos, 0,"ONLY"); |
303 | gMC->Gspos("C03O",2,"C03M", 0.,0.,+zfpos, 0,"ONLY"); |
304 | gMC->Gspos("C04O",1,"C04M", 0.,0.,-zfpos, 0,"ONLY"); |
305 | gMC->Gspos("C04O",2,"C04M", 0.,0.,+zfpos, 0,"ONLY"); |
306 | // |
307 | // Inner frame |
308 | tpar[0]= iChamber->RInner()-dframep; |
309 | tpar[1]= iChamber->RInner(); |
310 | tpar[2]= dframez/2; |
311 | gMC->Gsvolu("C03I", "TUBE", idAlu1, tpar, 3); |
312 | gMC->Gsvolu("C04I", "TUBE", idAlu1, tpar, 3); |
313 | |
314 | gMC->Gspos("C03I",1,"C03M", 0.,0.,-zfpos, 0,"ONLY"); |
315 | gMC->Gspos("C03I",2,"C03M", 0.,0.,+zfpos, 0,"ONLY"); |
316 | gMC->Gspos("C04I",1,"C04M", 0.,0.,-zfpos, 0,"ONLY"); |
317 | gMC->Gspos("C04I",2,"C04M", 0.,0.,+zfpos, 0,"ONLY"); |
318 | // |
319 | // Frame Crosses |
320 | if (frames) { |
321 | |
322 | bpar[0] = (iChamber->ROuter() - iChamber->RInner())/2; |
323 | bpar[1] = dframep/2; |
324 | bpar[2] = dframez/2; |
325 | gMC->Gsvolu("C03B", "BOX", idAlu1, bpar, 3); |
326 | gMC->Gsvolu("C04B", "BOX", idAlu1, bpar, 3); |
327 | |
328 | gMC->Gspos("C03B",1,"C03M", +iChamber->RInner()+bpar[0] , 0,-zfpos, |
329 | idrotm[1100],"ONLY"); |
330 | gMC->Gspos("C03B",2,"C03M", -iChamber->RInner()-bpar[0] , 0,-zfpos, |
331 | idrotm[1100],"ONLY"); |
332 | gMC->Gspos("C03B",3,"C03M", 0, +iChamber->RInner()+bpar[0] ,-zfpos, |
333 | idrotm[1101],"ONLY"); |
334 | gMC->Gspos("C03B",4,"C03M", 0, -iChamber->RInner()-bpar[0] ,-zfpos, |
335 | idrotm[1101],"ONLY"); |
336 | gMC->Gspos("C03B",5,"C03M", +iChamber->RInner()+bpar[0] , 0,+zfpos, |
337 | idrotm[1100],"ONLY"); |
338 | gMC->Gspos("C03B",6,"C03M", -iChamber->RInner()-bpar[0] , 0,+zfpos, |
339 | idrotm[1100],"ONLY"); |
340 | gMC->Gspos("C03B",7,"C03M", 0, +iChamber->RInner()+bpar[0] ,+zfpos, |
341 | idrotm[1101],"ONLY"); |
342 | gMC->Gspos("C03B",8,"C03M", 0, -iChamber->RInner()-bpar[0] ,+zfpos, |
343 | idrotm[1101],"ONLY"); |
344 | |
345 | gMC->Gspos("C04B",1,"C04M", +iChamber->RInner()+bpar[0] , 0,-zfpos, |
346 | idrotm[1100],"ONLY"); |
347 | gMC->Gspos("C04B",2,"C04M", -iChamber->RInner()-bpar[0] , 0,-zfpos, |
348 | idrotm[1100],"ONLY"); |
349 | gMC->Gspos("C04B",3,"C04M", 0, +iChamber->RInner()+bpar[0] ,-zfpos, |
350 | idrotm[1101],"ONLY"); |
351 | gMC->Gspos("C04B",4,"C04M", 0, -iChamber->RInner()-bpar[0] ,-zfpos, |
352 | idrotm[1101],"ONLY"); |
353 | gMC->Gspos("C04B",5,"C04M", +iChamber->RInner()+bpar[0] , 0,+zfpos, |
354 | idrotm[1100],"ONLY"); |
355 | gMC->Gspos("C04B",6,"C04M", -iChamber->RInner()-bpar[0] , 0,+zfpos, |
356 | idrotm[1100],"ONLY"); |
357 | gMC->Gspos("C04B",7,"C04M", 0, +iChamber->RInner()+bpar[0] ,+zfpos, |
358 | idrotm[1101],"ONLY"); |
359 | gMC->Gspos("C04B",8,"C04M", 0, -iChamber->RInner()-bpar[0] ,+zfpos, |
360 | idrotm[1101],"ONLY"); |
361 | } |
362 | // |
363 | // Chamber Material represented by Alu sheet |
364 | tpar[0]= iChamber->RInner(); |
365 | tpar[1]= iChamber->ROuter(); |
366 | tpar[2] = (iChamber->DGas()+iChamber->DAlu())/2; |
367 | gMC->Gsvolu("C03A", "TUBE", idAlu2, tpar, 3); |
368 | gMC->Gsvolu("C04A", "TUBE", idAlu2, tpar, 3); |
369 | gMC->Gspos("C03A", 1, "C03M", 0., 0., 0., 0, "ONLY"); |
370 | gMC->Gspos("C04A", 1, "C04M", 0., 0., 0., 0, "ONLY"); |
371 | // |
372 | // Sensitive volumes |
373 | // tpar[2] = iChamber->DGas(); |
374 | tpar[2] = iChamber->DGas()/2; |
375 | gMC->Gsvolu("C03G", "TUBE", idGas, tpar, 3); |
376 | gMC->Gsvolu("C04G", "TUBE", idGas, tpar, 3); |
377 | gMC->Gspos("C03G", 1, "C03A", 0., 0., 0., 0, "ONLY"); |
378 | gMC->Gspos("C04G", 1, "C04A", 0., 0., 0., 0, "ONLY"); |
379 | |
380 | if (frames) { |
381 | // |
382 | // Frame Crosses to be placed inside gas |
383 | dr = (iChamber->ROuter() - iChamber->RInner()); |
384 | bpar[0] = TMath::Sqrt(dr*dr-dframep*dframep/4)/2; |
385 | bpar[1] = dframep/2; |
386 | bpar[2] = iChamber->DGas()/2; |
387 | gMC->Gsvolu("C03F", "BOX", idAlu1, bpar, 3); |
388 | gMC->Gsvolu("C04F", "BOX", idAlu1, bpar, 3); |
389 | |
390 | gMC->Gspos("C03F",1,"C03G", +iChamber->RInner()+bpar[0] , 0, 0, |
391 | idrotm[1100],"ONLY"); |
392 | gMC->Gspos("C03F",2,"C03G", -iChamber->RInner()-bpar[0] , 0, 0, |
393 | idrotm[1100],"ONLY"); |
394 | gMC->Gspos("C03F",3,"C03G", 0, +iChamber->RInner()+bpar[0] , 0, |
395 | idrotm[1101],"ONLY"); |
396 | gMC->Gspos("C03F",4,"C03G", 0, -iChamber->RInner()-bpar[0] , 0, |
397 | idrotm[1101],"ONLY"); |
398 | |
399 | gMC->Gspos("C04F",1,"C04G", +iChamber->RInner()+bpar[0] , 0, 0, |
400 | idrotm[1100],"ONLY"); |
401 | gMC->Gspos("C04F",2,"C04G", -iChamber->RInner()-bpar[0] , 0, 0, |
402 | idrotm[1100],"ONLY"); |
403 | gMC->Gspos("C04F",3,"C04G", 0, +iChamber->RInner()+bpar[0] , 0, |
404 | idrotm[1101],"ONLY"); |
405 | gMC->Gspos("C04F",4,"C04G", 0, -iChamber->RInner()-bpar[0] , 0, |
406 | idrotm[1101],"ONLY"); |
407 | } |
408 | |
409 | //******************************************************************** |
410 | // Station 3 ** |
411 | //******************************************************************** |
412 | // indices 1 and 2 for first and second chambers in the station |
413 | // iChamber (first chamber) kept for other quanties than Z, |
414 | // assumed to be the same in both chambers |
415 | iChamber1 = iChamber = (AliMUONChamber*) (*fChambers)[4]; |
416 | iChamber2 =(AliMUONChamber*) (*fChambers)[5]; |
417 | zpos1=iChamber1->Z(); |
418 | zpos2=iChamber2->Z(); |
419 | dstation = zpos2 - zpos1; |
420 | |
421 | zfpos=-(iChamber->DGas()+dframez+iChamber->DAlu())/2; |
422 | // |
423 | // Mother volume |
424 | tpar[0] = iChamber->RInner()-dframep; |
425 | tpar[1] = (iChamber->ROuter()+dframep)/TMath::Cos(phi); |
426 | tpar[2] = dstation/4; |
427 | gMC->Gsvolu("C05M", "TUBE", idAir, tpar, 3); |
428 | gMC->Gsvolu("C06M", "TUBE", idAir, tpar, 3); |
429 | gMC->Gspos("C05M", 1, "ALIC", 0., 0., zpos1 , 0, "ONLY"); |
430 | gMC->Gspos("C06M", 1, "ALIC", 0., 0., zpos2 , 0, "ONLY"); |
431 | // Aluminium frames |
432 | // Outer frames |
433 | pgpar[0] = 360/12/2; |
434 | pgpar[1] = 360.; |
435 | pgpar[2] = 12.; |
436 | pgpar[3] = 2; |
437 | pgpar[4] = -dframez/2; |
438 | pgpar[5] = iChamber->ROuter(); |
439 | pgpar[6] = pgpar[5]+dframep; |
440 | pgpar[7] = +dframez/2; |
441 | pgpar[8] = pgpar[5]; |
442 | pgpar[9] = pgpar[6]; |
443 | gMC->Gsvolu("C05O", "PGON", idAlu1, pgpar, 10); |
444 | gMC->Gsvolu("C06O", "PGON", idAlu1, pgpar, 10); |
445 | gMC->Gspos("C05O",1,"C05M", 0.,0.,-zfpos, 0,"ONLY"); |
446 | gMC->Gspos("C05O",2,"C05M", 0.,0.,+zfpos, 0,"ONLY"); |
447 | gMC->Gspos("C06O",1,"C06M", 0.,0.,-zfpos, 0,"ONLY"); |
448 | gMC->Gspos("C06O",2,"C06M", 0.,0.,+zfpos, 0,"ONLY"); |
449 | // |
450 | // Inner frame |
451 | tpar[0]= iChamber->RInner()-dframep; |
452 | tpar[1]= iChamber->RInner(); |
453 | tpar[2]= dframez/2; |
454 | gMC->Gsvolu("C05I", "TUBE", idAlu1, tpar, 3); |
455 | gMC->Gsvolu("C06I", "TUBE", idAlu1, tpar, 3); |
456 | |
457 | gMC->Gspos("C05I",1,"C05M", 0.,0.,-zfpos, 0,"ONLY"); |
458 | gMC->Gspos("C05I",2,"C05M", 0.,0.,+zfpos, 0,"ONLY"); |
459 | gMC->Gspos("C06I",1,"C06M", 0.,0.,-zfpos, 0,"ONLY"); |
460 | gMC->Gspos("C06I",2,"C06M", 0.,0.,+zfpos, 0,"ONLY"); |
461 | // |
462 | // Frame Crosses |
463 | if (frames) { |
464 | bpar[0] = (iChamber->ROuter() - iChamber->RInner())/2; |
465 | bpar[1] = dframep/2; |
466 | bpar[2] = dframez/2; |
467 | gMC->Gsvolu("C05B", "BOX", idAlu1, bpar, 3); |
468 | gMC->Gsvolu("C06B", "BOX", idAlu1, bpar, 3); |
469 | |
470 | gMC->Gspos("C05B",1,"C05M", +iChamber->RInner()+bpar[0] , 0,-zfpos, |
471 | idrotm[1100],"ONLY"); |
472 | gMC->Gspos("C05B",2,"C05M", -iChamber->RInner()-bpar[0] , 0,-zfpos, |
473 | idrotm[1100],"ONLY"); |
474 | gMC->Gspos("C05B",3,"C05M", 0, +iChamber->RInner()+bpar[0] ,-zfpos, |
475 | idrotm[1101],"ONLY"); |
476 | gMC->Gspos("C05B",4,"C05M", 0, -iChamber->RInner()-bpar[0] ,-zfpos, |
477 | idrotm[1101],"ONLY"); |
478 | gMC->Gspos("C05B",5,"C05M", +iChamber->RInner()+bpar[0] , 0,+zfpos, |
479 | idrotm[1100],"ONLY"); |
480 | gMC->Gspos("C05B",6,"C05M", -iChamber->RInner()-bpar[0] , 0,+zfpos, |
481 | idrotm[1100],"ONLY"); |
482 | gMC->Gspos("C05B",7,"C05M", 0, +iChamber->RInner()+bpar[0] ,+zfpos, |
483 | idrotm[1101],"ONLY"); |
484 | gMC->Gspos("C05B",8,"C05M", 0, -iChamber->RInner()-bpar[0] ,+zfpos, |
485 | idrotm[1101],"ONLY"); |
486 | |
487 | gMC->Gspos("C06B",1,"C06M", +iChamber->RInner()+bpar[0] , 0,-zfpos, |
488 | idrotm[1100],"ONLY"); |
489 | gMC->Gspos("C06B",2,"C06M", -iChamber->RInner()-bpar[0] , 0,-zfpos, |
490 | idrotm[1100],"ONLY"); |
491 | gMC->Gspos("C06B",3,"C06M", 0, +iChamber->RInner()+bpar[0] ,-zfpos, |
492 | idrotm[1101],"ONLY"); |
493 | gMC->Gspos("C06B",4,"C06M", 0, -iChamber->RInner()-bpar[0] ,-zfpos, |
494 | idrotm[1101],"ONLY"); |
495 | gMC->Gspos("C06B",5,"C06M", +iChamber->RInner()+bpar[0] , 0,+zfpos, |
496 | idrotm[1100],"ONLY"); |
497 | gMC->Gspos("C06B",6,"C06M", -iChamber->RInner()-bpar[0] , 0,+zfpos, |
498 | idrotm[1100],"ONLY"); |
499 | gMC->Gspos("C06B",7,"C06M", 0, +iChamber->RInner()+bpar[0] ,+zfpos, |
500 | idrotm[1101],"ONLY"); |
501 | gMC->Gspos("C06B",8,"C06M", 0, -iChamber->RInner()-bpar[0] ,+zfpos, |
502 | idrotm[1101],"ONLY"); |
503 | } |
504 | |
505 | |
506 | // |
507 | // Chamber Material represented by Alu sheet |
508 | tpar[0]= iChamber->RInner(); |
509 | tpar[1]= iChamber->ROuter(); |
510 | tpar[2] = (iChamber->DGas()+iChamber->DAlu())/2; |
511 | gMC->Gsvolu("C05A", "TUBE", idAlu2, tpar, 3); |
512 | gMC->Gsvolu("C06A", "TUBE", idAlu2, tpar, 3); |
513 | gMC->Gspos("C05A", 1, "C05M", 0., 0., 0., 0, "ONLY"); |
514 | gMC->Gspos("C06A", 1, "C06M", 0., 0., 0., 0, "ONLY"); |
515 | // |
516 | // Sensitive volumes |
517 | tpar[2] = iChamber->DGas()/2.; |
518 | gMC->Gsvolu("C05G", "TUBE", idGas, tpar, 3); |
519 | gMC->Gsvolu("C06G", "TUBE", idGas, tpar, 3); |
520 | gMC->Gspos("C05G", 1, "C05A", 0., 0., 0., 0, "ONLY"); |
521 | gMC->Gspos("C06G", 1, "C06A", 0., 0., 0., 0, "ONLY"); |
522 | // |
523 | // Frame Crosses to be placed inside gas |
524 | if (frames) { |
525 | dr = (iChamber->ROuter() - iChamber->RInner()); |
526 | bpar[0] = TMath::Sqrt(dr*dr-dframep*dframep/4)/2; |
527 | bpar[1] = dframep/2; |
528 | bpar[2] = iChamber->DGas()/2; |
529 | gMC->Gsvolu("C05F", "BOX", idAlu1, bpar, 3); |
530 | gMC->Gsvolu("C06F", "BOX", idAlu1, bpar, 3); |
531 | |
532 | gMC->Gspos("C05F",1,"C05G", +iChamber->RInner()+bpar[0] , 0, 0, |
533 | idrotm[1100],"ONLY"); |
534 | gMC->Gspos("C05F",2,"C05G", -iChamber->RInner()-bpar[0] , 0, 0, |
535 | idrotm[1100],"ONLY"); |
536 | gMC->Gspos("C05F",3,"C05G", 0, +iChamber->RInner()+bpar[0] , 0, |
537 | idrotm[1101],"ONLY"); |
538 | gMC->Gspos("C05F",4,"C05G", 0, -iChamber->RInner()-bpar[0] , 0, |
539 | idrotm[1101],"ONLY"); |
540 | |
541 | gMC->Gspos("C06F",1,"C06G", +iChamber->RInner()+bpar[0] , 0, 0, |
542 | idrotm[1100],"ONLY"); |
543 | gMC->Gspos("C06F",2,"C06G", -iChamber->RInner()-bpar[0] , 0, 0, |
544 | idrotm[1100],"ONLY"); |
545 | gMC->Gspos("C06F",3,"C06G", 0, +iChamber->RInner()+bpar[0] , 0, |
546 | idrotm[1101],"ONLY"); |
547 | gMC->Gspos("C06F",4,"C06G", 0, -iChamber->RInner()-bpar[0] , 0, |
548 | idrotm[1101],"ONLY"); |
549 | } |
550 | |
551 | //******************************************************************** |
552 | // Station 4 ** |
553 | //******************************************************************** |
554 | // indices 1 and 2 for first and second chambers in the station |
555 | // iChamber (first chamber) kept for other quanties than Z, |
556 | // assumed to be the same in both chambers |
557 | iChamber1 = iChamber = (AliMUONChamber*) (*fChambers)[6]; |
558 | iChamber2 =(AliMUONChamber*) (*fChambers)[7]; |
559 | zpos1=iChamber1->Z(); |
560 | zpos2=iChamber2->Z(); |
561 | dstation = zpos2 - zpos1; |
562 | zfpos=-(iChamber->DGas()+dframez+iChamber->DAlu())/2; |
563 | |
564 | // |
565 | // Mother volume |
566 | tpar[0] = iChamber->RInner()-dframep; |
567 | tpar[1] = (iChamber->ROuter()+dframep)/TMath::Cos(phi); |
568 | tpar[2] = dstation/4; |
569 | |
570 | gMC->Gsvolu("C07M", "TUBE", idAir, tpar, 3); |
571 | gMC->Gsvolu("C08M", "TUBE", idAir, tpar, 3); |
572 | gMC->Gspos("C07M", 1, "ALIC", 0., 0., zpos1 , 0, "ONLY"); |
573 | gMC->Gspos("C08M", 1, "ALIC", 0., 0., zpos2 , 0, "ONLY"); |
574 | // Aluminium frames |
575 | // Outer frames |
576 | pgpar[0] = 360/12/2; |
577 | pgpar[1] = 360.; |
578 | pgpar[2] = 12.; |
579 | pgpar[3] = 2; |
580 | pgpar[4] = -dframez/2; |
581 | pgpar[5] = iChamber->ROuter(); |
582 | pgpar[6] = pgpar[5]+dframep; |
583 | pgpar[7] = +dframez/2; |
584 | pgpar[8] = pgpar[5]; |
585 | pgpar[9] = pgpar[6]; |
586 | gMC->Gsvolu("C07O", "PGON", idAlu1, pgpar, 10); |
587 | gMC->Gsvolu("C08O", "PGON", idAlu1, pgpar, 10); |
588 | gMC->Gspos("C07O",1,"C07M", 0.,0.,-zfpos, 0,"ONLY"); |
589 | gMC->Gspos("C07O",2,"C07M", 0.,0.,+zfpos, 0,"ONLY"); |
590 | gMC->Gspos("C08O",1,"C08M", 0.,0.,-zfpos, 0,"ONLY"); |
591 | gMC->Gspos("C08O",2,"C08M", 0.,0.,+zfpos, 0,"ONLY"); |
592 | // |
593 | // Inner frame |
594 | tpar[0]= iChamber->RInner()-dframep; |
595 | tpar[1]= iChamber->RInner(); |
596 | tpar[2]= dframez/2; |
597 | gMC->Gsvolu("C07I", "TUBE", idAlu1, tpar, 3); |
598 | gMC->Gsvolu("C08I", "TUBE", idAlu1, tpar, 3); |
599 | |
600 | gMC->Gspos("C07I",1,"C07M", 0.,0.,-zfpos, 0,"ONLY"); |
601 | gMC->Gspos("C07I",2,"C07M", 0.,0.,+zfpos, 0,"ONLY"); |
602 | gMC->Gspos("C08I",1,"C08M", 0.,0.,-zfpos, 0,"ONLY"); |
603 | gMC->Gspos("C08I",2,"C08M", 0.,0.,+zfpos, 0,"ONLY"); |
604 | // |
605 | // Frame Crosses |
606 | if (frames) { |
607 | bpar[0] = (iChamber->ROuter() - iChamber->RInner())/2; |
608 | bpar[1] = dframep/2; |
609 | bpar[2] = dframez/2; |
610 | gMC->Gsvolu("C07B", "BOX", idAlu1, bpar, 3); |
611 | gMC->Gsvolu("C08B", "BOX", idAlu1, bpar, 3); |
612 | |
613 | gMC->Gspos("C07B",1,"C07M", +iChamber->RInner()+bpar[0] , 0,-zfpos, |
614 | idrotm[1100],"ONLY"); |
615 | gMC->Gspos("C07B",2,"C07M", -iChamber->RInner()-bpar[0] , 0,-zfpos, |
616 | idrotm[1100],"ONLY"); |
617 | gMC->Gspos("C07B",3,"C07M", 0, +iChamber->RInner()+bpar[0] ,-zfpos, |
618 | idrotm[1101],"ONLY"); |
619 | gMC->Gspos("C07B",4,"C07M", 0, -iChamber->RInner()-bpar[0] ,-zfpos, |
620 | idrotm[1101],"ONLY"); |
621 | gMC->Gspos("C07B",5,"C07M", +iChamber->RInner()+bpar[0] , 0,+zfpos, |
622 | idrotm[1100],"ONLY"); |
623 | gMC->Gspos("C07B",6,"C07M", -iChamber->RInner()-bpar[0] , 0,+zfpos, |
624 | idrotm[1100],"ONLY"); |
625 | gMC->Gspos("C07B",7,"C07M", 0, +iChamber->RInner()+bpar[0] ,+zfpos, |
626 | idrotm[1101],"ONLY"); |
627 | gMC->Gspos("C07B",8,"C07M", 0, -iChamber->RInner()-bpar[0] ,+zfpos, |
628 | idrotm[1101],"ONLY"); |
629 | |
630 | gMC->Gspos("C08B",1,"C08M", +iChamber->RInner()+bpar[0] , 0,-zfpos, |
631 | idrotm[1100],"ONLY"); |
632 | gMC->Gspos("C08B",2,"C08M", -iChamber->RInner()-bpar[0] , 0,-zfpos, |
633 | idrotm[1100],"ONLY"); |
634 | gMC->Gspos("C08B",3,"C08M", 0, +iChamber->RInner()+bpar[0] ,-zfpos, |
635 | idrotm[1101],"ONLY"); |
636 | gMC->Gspos("C08B",4,"C08M", 0, -iChamber->RInner()-bpar[0] ,-zfpos, |
637 | idrotm[1101],"ONLY"); |
638 | gMC->Gspos("C08B",5,"C08M", +iChamber->RInner()+bpar[0] , 0,+zfpos, |
639 | idrotm[1100],"ONLY"); |
640 | gMC->Gspos("C08B",6,"C08M", -iChamber->RInner()-bpar[0] , 0,+zfpos, |
641 | idrotm[1100],"ONLY"); |
642 | gMC->Gspos("C08B",7,"C08M", 0, +iChamber->RInner()+bpar[0] ,+zfpos, |
643 | idrotm[1101],"ONLY"); |
644 | gMC->Gspos("C08B",8,"C08M", 0, -iChamber->RInner()-bpar[0] ,+zfpos, |
645 | idrotm[1101],"ONLY"); |
646 | } |
647 | |
648 | |
649 | // |
650 | // Chamber Material represented by Alu sheet |
651 | tpar[0]= iChamber->RInner(); |
652 | tpar[1]= iChamber->ROuter(); |
653 | tpar[2] = (iChamber->DGas()+iChamber->DAlu())/2; |
654 | gMC->Gsvolu("C07A", "TUBE", idAlu2, tpar, 3); |
655 | gMC->Gsvolu("C08A", "TUBE", idAlu2, tpar, 3); |
656 | gMC->Gspos("C07A", 1, "C07M", 0., 0., 0., 0, "ONLY"); |
657 | gMC->Gspos("C08A", 1, "C08M", 0., 0., 0., 0, "ONLY"); |
658 | // |
659 | // Sensitive volumes |
660 | // tpar[2] = iChamber->DGas(); |
661 | tpar[2] = iChamber->DGas()/2; |
662 | gMC->Gsvolu("C07G", "TUBE", idGas, tpar, 3); |
663 | gMC->Gsvolu("C08G", "TUBE", idGas, tpar, 3); |
664 | gMC->Gspos("C07G", 1, "C07A", 0., 0., 0., 0, "ONLY"); |
665 | gMC->Gspos("C08G", 1, "C08A", 0., 0., 0., 0, "ONLY"); |
666 | // |
667 | // Frame Crosses to be placed inside gas |
668 | if (frames) { |
669 | dr = (iChamber->ROuter() - iChamber->RInner()); |
670 | bpar[0] = TMath::Sqrt(dr*dr-dframep*dframep/4)/2; |
671 | bpar[1] = dframep/2; |
672 | bpar[2] = iChamber->DGas()/2; |
673 | gMC->Gsvolu("C07F", "BOX", idAlu1, bpar, 3); |
674 | gMC->Gsvolu("C08F", "BOX", idAlu1, bpar, 3); |
675 | |
676 | gMC->Gspos("C07F",1,"C07G", +iChamber->RInner()+bpar[0] , 0, 0, |
677 | idrotm[1100],"ONLY"); |
678 | gMC->Gspos("C07F",2,"C07G", -iChamber->RInner()-bpar[0] , 0, 0, |
679 | idrotm[1100],"ONLY"); |
680 | gMC->Gspos("C07F",3,"C07G", 0, +iChamber->RInner()+bpar[0] , 0, |
681 | idrotm[1101],"ONLY"); |
682 | gMC->Gspos("C07F",4,"C07G", 0, -iChamber->RInner()-bpar[0] , 0, |
683 | idrotm[1101],"ONLY"); |
684 | |
685 | gMC->Gspos("C08F",1,"C08G", +iChamber->RInner()+bpar[0] , 0, 0, |
686 | idrotm[1100],"ONLY"); |
687 | gMC->Gspos("C08F",2,"C08G", -iChamber->RInner()-bpar[0] , 0, 0, |
688 | idrotm[1100],"ONLY"); |
689 | gMC->Gspos("C08F",3,"C08G", 0, +iChamber->RInner()+bpar[0] , 0, |
690 | idrotm[1101],"ONLY"); |
691 | gMC->Gspos("C08F",4,"C08G", 0, -iChamber->RInner()-bpar[0] , 0, |
692 | idrotm[1101],"ONLY"); |
693 | } |
694 | //******************************************************************** |
695 | // Station 5 ** |
696 | //******************************************************************** |
697 | // indices 1 and 2 for first and second chambers in the station |
698 | // iChamber (first chamber) kept for other quanties than Z, |
699 | // assumed to be the same in both chambers |
700 | iChamber1 = iChamber = (AliMUONChamber*) (*fChambers)[8]; |
701 | iChamber2 =(AliMUONChamber*) (*fChambers)[9]; |
702 | zpos1=iChamber1->Z(); |
703 | zpos2=iChamber2->Z(); |
704 | dstation = zpos2 - zpos1; |
705 | zfpos=-(iChamber->DGas()+dframez+iChamber->DAlu())/2; |
706 | |
707 | // |
708 | // Mother volume |
709 | tpar[0] = iChamber->RInner()-dframep; |
710 | tpar[1] = (iChamber->ROuter()+dframep)/TMath::Cos(phi); |
711 | tpar[2] = dstation/4; |
712 | |
713 | gMC->Gsvolu("C09M", "TUBE", idAir, tpar, 3); |
714 | gMC->Gsvolu("C10M", "TUBE", idAir, tpar, 3); |
715 | gMC->Gspos("C09M", 1, "ALIC", 0., 0., zpos1 , 0, "ONLY"); |
716 | gMC->Gspos("C10M", 1, "ALIC", 0., 0., zpos2 , 0, "ONLY"); |
717 | // Aluminium frames |
718 | // Outer frames |
719 | pgpar[0] = 360/12/2; |
720 | pgpar[1] = 360.; |
721 | pgpar[2] = 12.; |
722 | pgpar[3] = 2; |
723 | pgpar[4] = -dframez/2; |
724 | pgpar[5] = iChamber->ROuter(); |
725 | pgpar[6] = pgpar[5]+dframep; |
726 | pgpar[7] = +dframez/2; |
727 | pgpar[8] = pgpar[5]; |
728 | pgpar[9] = pgpar[6]; |
729 | gMC->Gsvolu("C09O", "PGON", idAlu1, pgpar, 10); |
730 | gMC->Gsvolu("C10O", "PGON", idAlu1, pgpar, 10); |
731 | gMC->Gspos("C09O",1,"C09M", 0.,0.,-zfpos, 0,"ONLY"); |
732 | gMC->Gspos("C09O",2,"C09M", 0.,0.,+zfpos, 0,"ONLY"); |
733 | gMC->Gspos("C10O",1,"C10M", 0.,0.,-zfpos, 0,"ONLY"); |
734 | gMC->Gspos("C10O",2,"C10M", 0.,0.,+zfpos, 0,"ONLY"); |
735 | // |
736 | // Inner frame |
737 | tpar[0]= iChamber->RInner()-dframep; |
738 | tpar[1]= iChamber->RInner(); |
739 | tpar[2]= dframez/2; |
740 | gMC->Gsvolu("C09I", "TUBE", idAlu1, tpar, 3); |
741 | gMC->Gsvolu("C10I", "TUBE", idAlu1, tpar, 3); |
742 | |
743 | gMC->Gspos("C09I",1,"C09M", 0.,0.,-zfpos, 0,"ONLY"); |
744 | gMC->Gspos("C09I",2,"C09M", 0.,0.,+zfpos, 0,"ONLY"); |
745 | gMC->Gspos("C10I",1,"C10M", 0.,0.,-zfpos, 0,"ONLY"); |
746 | gMC->Gspos("C10I",2,"C10M", 0.,0.,+zfpos, 0,"ONLY"); |
747 | |
748 | if (frames) { |
749 | // |
750 | // Frame Crosses |
751 | |
752 | bpar[0] = (iChamber->ROuter() - iChamber->RInner())/2; |
753 | bpar[1] = dframep/2; |
754 | bpar[2] = dframez/2; |
755 | gMC->Gsvolu("C09B", "BOX", idAlu1, bpar, 3); |
756 | gMC->Gsvolu("C10B", "BOX", idAlu1, bpar, 3); |
757 | |
758 | gMC->Gspos("C09B",1,"C09M", +iChamber->RInner()+bpar[0] , 0,-zfpos, |
759 | idrotm[1100],"ONLY"); |
760 | gMC->Gspos("C09B",2,"C09M", -iChamber->RInner()-bpar[0] , 0,-zfpos, |
761 | idrotm[1100],"ONLY"); |
762 | gMC->Gspos("C09B",3,"C09M", 0, +iChamber->RInner()+bpar[0] ,-zfpos, |
763 | idrotm[1101],"ONLY"); |
764 | gMC->Gspos("C09B",4,"C09M", 0, -iChamber->RInner()-bpar[0] ,-zfpos, |
765 | idrotm[1101],"ONLY"); |
766 | gMC->Gspos("C09B",5,"C09M", +iChamber->RInner()+bpar[0] , 0,+zfpos, |
767 | idrotm[1100],"ONLY"); |
768 | gMC->Gspos("C09B",6,"C09M", -iChamber->RInner()-bpar[0] , 0,+zfpos, |
769 | idrotm[1100],"ONLY"); |
770 | gMC->Gspos("C09B",7,"C09M", 0, +iChamber->RInner()+bpar[0] ,+zfpos, |
771 | idrotm[1101],"ONLY"); |
772 | gMC->Gspos("C09B",8,"C09M", 0, -iChamber->RInner()-bpar[0] ,+zfpos, |
773 | idrotm[1101],"ONLY"); |
774 | |
775 | gMC->Gspos("C10B",1,"C10M", +iChamber->RInner()+bpar[0] , 0,-zfpos, |
776 | idrotm[1100],"ONLY"); |
777 | gMC->Gspos("C10B",2,"C10M", -iChamber->RInner()-bpar[0] , 0,-zfpos, |
778 | idrotm[1100],"ONLY"); |
779 | gMC->Gspos("C10B",3,"C10M", 0, +iChamber->RInner()+bpar[0] ,-zfpos, |
780 | idrotm[1101],"ONLY"); |
781 | gMC->Gspos("C10B",4,"C10M", 0, -iChamber->RInner()-bpar[0] ,-zfpos, |
782 | idrotm[1101],"ONLY"); |
783 | gMC->Gspos("C10B",5,"C10M", +iChamber->RInner()+bpar[0] , 0,+zfpos, |
784 | idrotm[1100],"ONLY"); |
785 | gMC->Gspos("C10B",6,"C10M", -iChamber->RInner()-bpar[0] , 0,+zfpos, |
786 | idrotm[1100],"ONLY"); |
787 | gMC->Gspos("C10B",7,"C10M", 0, +iChamber->RInner()+bpar[0] ,+zfpos, |
788 | idrotm[1101],"ONLY"); |
789 | gMC->Gspos("C10B",8,"C10M", 0, -iChamber->RInner()-bpar[0] ,+zfpos, |
790 | idrotm[1101],"ONLY"); |
791 | } |
792 | |
793 | |
794 | // |
795 | // Chamber Material represented by Alu sheet |
796 | tpar[0]= iChamber->RInner(); |
797 | tpar[1]= iChamber->ROuter(); |
798 | tpar[2] = (iChamber->DGas()+iChamber->DAlu())/2; |
799 | gMC->Gsvolu("C09A", "TUBE", idAlu2, tpar, 3); |
800 | gMC->Gsvolu("C10A", "TUBE", idAlu2, tpar, 3); |
801 | gMC->Gspos("C09A", 1, "C09M", 0., 0., 0., 0, "ONLY"); |
802 | gMC->Gspos("C10A", 1, "C10M", 0., 0., 0., 0, "ONLY"); |
803 | // |
804 | // Sensitive volumes |
805 | // tpar[2] = iChamber->DGas(); |
806 | tpar[2] = iChamber->DGas()/2; |
807 | gMC->Gsvolu("C09G", "TUBE", idGas, tpar, 3); |
808 | gMC->Gsvolu("C10G", "TUBE", idGas, tpar, 3); |
809 | gMC->Gspos("C09G", 1, "C09A", 0., 0., 0., 0, "ONLY"); |
810 | gMC->Gspos("C10G", 1, "C10A", 0., 0., 0., 0, "ONLY"); |
811 | // |
812 | // Frame Crosses to be placed inside gas |
813 | if (frames) { |
814 | dr = (iChamber->ROuter() - iChamber->RInner()); |
815 | bpar[0] = TMath::Sqrt(dr*dr-dframep*dframep/4)/2; |
816 | bpar[1] = dframep/2; |
817 | bpar[2] = iChamber->DGas()/2; |
818 | gMC->Gsvolu("C09F", "BOX", idAlu1, bpar, 3); |
819 | gMC->Gsvolu("C10F", "BOX", idAlu1, bpar, 3); |
820 | |
821 | gMC->Gspos("C09F",1,"C09G", +iChamber->RInner()+bpar[0] , 0, 0, |
822 | idrotm[1100],"ONLY"); |
823 | gMC->Gspos("C09F",2,"C09G", -iChamber->RInner()-bpar[0] , 0, 0, |
824 | idrotm[1100],"ONLY"); |
825 | gMC->Gspos("C09F",3,"C09G", 0, +iChamber->RInner()+bpar[0] , 0, |
826 | idrotm[1101],"ONLY"); |
827 | gMC->Gspos("C09F",4,"C09G", 0, -iChamber->RInner()-bpar[0] , 0, |
828 | idrotm[1101],"ONLY"); |
829 | |
830 | gMC->Gspos("C10F",1,"C10G", +iChamber->RInner()+bpar[0] , 0, 0, |
831 | idrotm[1100],"ONLY"); |
832 | gMC->Gspos("C10F",2,"C10G", -iChamber->RInner()-bpar[0] , 0, 0, |
833 | idrotm[1100],"ONLY"); |
834 | gMC->Gspos("C10F",3,"C10G", 0, +iChamber->RInner()+bpar[0] , 0, |
835 | idrotm[1101],"ONLY"); |
836 | gMC->Gspos("C10F",4,"C10G", 0, -iChamber->RInner()-bpar[0] , 0, |
837 | idrotm[1101],"ONLY"); |
838 | } |
839 | |
840 | /////////////////////////////////////// |
841 | // GEOMETRY FOR THE TRIGGER CHAMBERS // |
842 | /////////////////////////////////////// |
843 | |
844 | // 03/00 P. Dupieux : introduce a slighly more realistic |
845 | // geom. of the trigger readout planes with |
846 | // 2 Zpos per trigger plane (alternate |
847 | // between left and right of the trigger) |
848 | |
849 | // Parameters of the Trigger Chambers |
850 | |
851 | |
852 | const Float_t kXMC1MIN=34.; |
853 | const Float_t kXMC1MED=51.; |
854 | const Float_t kXMC1MAX=272.; |
855 | const Float_t kYMC1MIN=34.; |
856 | const Float_t kYMC1MAX=51.; |
857 | const Float_t kRMIN1=50.; |
858 | const Float_t kRMAX1=62.; |
859 | const Float_t kRMIN2=50.; |
860 | const Float_t kRMAX2=66.; |
861 | |
862 | // zposition of the middle of the gas gap in mother vol |
863 | const Float_t kZMCm=-3.6; |
864 | const Float_t kZMCp=+3.6; |
865 | |
866 | |
867 | // TRIGGER STATION 1 - TRIGGER STATION 1 - TRIGGER STATION 1 |
868 | |
869 | // iChamber 1 and 2 for first and second chambers in the station |
870 | // iChamber (first chamber) kept for other quanties than Z, |
871 | // assumed to be the same in both chambers |
872 | iChamber1 = iChamber = (AliMUONChamber*) (*fChambers)[10]; |
873 | iChamber2 =(AliMUONChamber*) (*fChambers)[11]; |
874 | |
875 | // 03/00 |
876 | // zpos1 and zpos2 are now the middle of the first and second |
877 | // plane of station 1 : |
878 | // zpos1=(16075+15995)/2=16035 mm, thick/2=40 mm |
879 | // zpos2=(16225+16145)/2=16185 mm, thick/2=40 mm |
880 | // |
881 | // zpos1m=15999 mm , zpos1p=16071 mm (middles of gas gaps) |
882 | // zpos2m=16149 mm , zpos2p=16221 mm (middles of gas gaps) |
883 | // rem : the total thickness accounts for 1 mm of al on both |
884 | // side of the RPCs (see zpos1 and zpos2), as previously |
885 | |
886 | zpos1=iChamber1->Z(); |
887 | zpos2=iChamber2->Z(); |
888 | |
889 | |
890 | // Mother volume definition |
891 | tpar[0] = iChamber->RInner(); |
892 | tpar[1] = iChamber->ROuter(); |
893 | tpar[2] = 4.0; |
894 | gMC->Gsvolu("CM11", "TUBE", idAir, tpar, 3); |
895 | gMC->Gsvolu("CM12", "TUBE", idAir, tpar, 3); |
896 | |
897 | // Definition of the flange between the beam shielding and the RPC |
898 | tpar[0]= kRMIN1; |
899 | tpar[1]= kRMAX1; |
900 | tpar[2]= 4.0; |
901 | |
902 | gMC->Gsvolu("CF1A", "TUBE", idAlu1, tpar, 3); //Al |
903 | gMC->Gspos("CF1A", 1, "CM11", 0., 0., 0., 0, "MANY"); |
904 | gMC->Gspos("CF1A", 2, "CM12", 0., 0., 0., 0, "MANY"); |
905 | |
906 | |
907 | // FIRST PLANE OF STATION 1 |
908 | |
909 | // ratios of zpos1m/zpos1p and inverse for first plane |
910 | Float_t zmp=(zpos1-3.6)/(zpos1+3.6); |
911 | Float_t zpm=1./zmp; |
912 | |
913 | |
914 | // Definition of prototype for chambers in the first plane |
915 | |
916 | tpar[0]= 0.; |
917 | tpar[1]= 0.; |
918 | tpar[2]= 0.; |
919 | |
920 | gMC->Gsvolu("CC1A", "BOX ", idAlu1, tpar, 0); //Al |
921 | gMC->Gsvolu("CB1A", "BOX ", idtmed[1107], tpar, 0); //Bakelite |
922 | gMC->Gsvolu("CG1A", "BOX ", idtmed[1106], tpar, 0); //Gas streamer |
923 | |
924 | // chamber type A |
925 | tpar[0] = -1.; |
926 | tpar[1] = -1.; |
927 | |
928 | const Float_t kXMC1A=kXMC1MED+(kXMC1MAX-kXMC1MED)/2.; |
929 | const Float_t kYMC1Am=0.; |
930 | const Float_t kYMC1Ap=0.; |
931 | |
932 | tpar[2] = 0.1; |
933 | gMC->Gsposp("CG1A", 1, "CB1A", 0., 0., 0., 0, "ONLY",tpar,3); |
934 | tpar[2] = 0.3; |
935 | gMC->Gsposp("CB1A", 1, "CC1A", 0., 0., 0., 0, "ONLY",tpar,3); |
936 | |
937 | tpar[2] = 0.4; |
938 | tpar[0] = (kXMC1MAX-kXMC1MED)/2.; |
939 | tpar[1] = kYMC1MIN; |
940 | |
941 | gMC->Gsposp("CC1A", 1, "CM11",kXMC1A,kYMC1Am,kZMCm, 0, "ONLY", tpar, 3); |
942 | gMC->Gsposp("CC1A", 2, "CM11",-kXMC1A,kYMC1Ap,kZMCp, 0, "ONLY", tpar, 3); |
943 | |
944 | // chamber type B |
945 | Float_t tpar1save=tpar[1]; |
946 | Float_t y1msave=kYMC1Am; |
947 | Float_t y1psave=kYMC1Ap; |
948 | |
949 | tpar[0] = (kXMC1MAX-kXMC1MIN)/2.; |
950 | tpar[1] = (kYMC1MAX-kYMC1MIN)/2.; |
951 | |
952 | const Float_t kXMC1B=kXMC1MIN+tpar[0]; |
953 | const Float_t kYMC1Bp=(y1msave+tpar1save)*zpm+tpar[1]; |
954 | const Float_t kYMC1Bm=(y1psave+tpar1save)*zmp+tpar[1]; |
955 | |
956 | gMC->Gsposp("CC1A", 3, "CM11",kXMC1B,kYMC1Bp,kZMCp, 0, "ONLY", tpar, 3); |
957 | gMC->Gsposp("CC1A", 4, "CM11",-kXMC1B,kYMC1Bm,kZMCm, 0, "ONLY", tpar, 3); |
958 | gMC->Gsposp("CC1A", 5, "CM11",kXMC1B,-kYMC1Bp,kZMCp, 0, "ONLY", tpar, 3); |
959 | gMC->Gsposp("CC1A", 6, "CM11",-kXMC1B,-kYMC1Bm,kZMCm, 0, "ONLY", tpar, 3); |
960 | |
961 | // chamber type C (end of type B !!) |
962 | tpar1save=tpar[1]; |
963 | y1msave=kYMC1Bm; |
964 | y1psave=kYMC1Bp; |
965 | |
966 | tpar[0] = kXMC1MAX/2; |
967 | tpar[1] = kYMC1MAX/2; |
968 | |
969 | const Float_t kXMC1C=tpar[0]; |
970 | // warning : same Z than type B |
971 | const Float_t kYMC1Cp=(y1psave+tpar1save)*1.+tpar[1]; |
972 | const Float_t kYMC1Cm=(y1msave+tpar1save)*1.+tpar[1]; |
973 | |
974 | gMC->Gsposp("CC1A", 7, "CM11",kXMC1C,kYMC1Cp,kZMCp, 0, "ONLY", tpar, 3); |
975 | gMC->Gsposp("CC1A", 8, "CM11",-kXMC1C,kYMC1Cm,kZMCm, 0, "ONLY", tpar, 3); |
976 | gMC->Gsposp("CC1A", 9, "CM11",kXMC1C,-kYMC1Cp,kZMCp, 0, "ONLY", tpar, 3); |
977 | gMC->Gsposp("CC1A", 10, "CM11",-kXMC1C,-kYMC1Cm,kZMCm, 0, "ONLY", tpar, 3); |
978 | |
979 | // chamber type D, E and F (same size) |
980 | tpar1save=tpar[1]; |
981 | y1msave=kYMC1Cm; |
982 | y1psave=kYMC1Cp; |
983 | |
984 | tpar[0] = kXMC1MAX/2.; |
985 | tpar[1] = kYMC1MIN; |
986 | |
987 | const Float_t kXMC1D=tpar[0]; |
988 | const Float_t kYMC1Dp=(y1msave+tpar1save)*zpm+tpar[1]; |
989 | const Float_t kYMC1Dm=(y1psave+tpar1save)*zmp+tpar[1]; |
990 | |
991 | gMC->Gsposp("CC1A", 11, "CM11",kXMC1D,kYMC1Dm,kZMCm, 0, "ONLY", tpar, 3); |
992 | gMC->Gsposp("CC1A", 12, "CM11",-kXMC1D,kYMC1Dp,kZMCp, 0, "ONLY", tpar, 3); |
993 | gMC->Gsposp("CC1A", 13, "CM11",kXMC1D,-kYMC1Dm,kZMCm, 0, "ONLY", tpar, 3); |
994 | gMC->Gsposp("CC1A", 14, "CM11",-kXMC1D,-kYMC1Dp,kZMCp, 0, "ONLY", tpar, 3); |
995 | |
996 | |
997 | tpar1save=tpar[1]; |
998 | y1msave=kYMC1Dm; |
999 | y1psave=kYMC1Dp; |
1000 | const Float_t kYMC1Ep=(y1msave+tpar1save)*zpm+tpar[1]; |
1001 | const Float_t kYMC1Em=(y1psave+tpar1save)*zmp+tpar[1]; |
1002 | |
1003 | gMC->Gsposp("CC1A", 15, "CM11",kXMC1D,kYMC1Ep,kZMCp, 0, "ONLY", tpar, 3); |
1004 | gMC->Gsposp("CC1A", 16, "CM11",-kXMC1D,kYMC1Em,kZMCm, 0, "ONLY", tpar, 3); |
1005 | gMC->Gsposp("CC1A", 17, "CM11",kXMC1D,-kYMC1Ep,kZMCp, 0, "ONLY", tpar, 3); |
1006 | gMC->Gsposp("CC1A", 18, "CM11",-kXMC1D,-kYMC1Em,kZMCm, 0, "ONLY", tpar, 3); |
1007 | |
1008 | tpar1save=tpar[1]; |
1009 | y1msave=kYMC1Em; |
1010 | y1psave=kYMC1Ep; |
1011 | const Float_t kYMC1Fp=(y1msave+tpar1save)*zpm+tpar[1]; |
1012 | const Float_t kYMC1Fm=(y1psave+tpar1save)*zmp+tpar[1]; |
1013 | |
1014 | gMC->Gsposp("CC1A", 19, "CM11",kXMC1D,kYMC1Fm,kZMCm, 0, "ONLY", tpar, 3); |
1015 | gMC->Gsposp("CC1A", 20, "CM11",-kXMC1D,kYMC1Fp,kZMCp, 0, "ONLY", tpar, 3); |
1016 | gMC->Gsposp("CC1A", 21, "CM11",kXMC1D,-kYMC1Fm,kZMCm, 0, "ONLY", tpar, 3); |
1017 | gMC->Gsposp("CC1A", 22, "CM11",-kXMC1D,-kYMC1Fp,kZMCp, 0, "ONLY", tpar, 3); |
1018 | |
1019 | // Positioning first plane in ALICE |
1020 | gMC->Gspos("CM11", 1, "ALIC", 0., 0., zpos1, 0, "ONLY"); |
1021 | |
1022 | // End of geometry definition for the first plane of station 1 |
1023 | |
1024 | |
1025 | |
1026 | // SECOND PLANE OF STATION 1 : proj ratio = zpos2/zpos1 |
1027 | |
1028 | const Float_t kZ12=zpos2/zpos1; |
1029 | |
1030 | // Definition of prototype for chambers in the second plane of station 1 |
1031 | |
1032 | tpar[0]= 0.; |
1033 | tpar[1]= 0.; |
1034 | tpar[2]= 0.; |
1035 | |
1036 | gMC->Gsvolu("CC2A", "BOX ", idAlu1, tpar, 0); //Al |
1037 | gMC->Gsvolu("CB2A", "BOX ", idtmed[1107], tpar, 0); //Bakelite |
1038 | gMC->Gsvolu("CG2A", "BOX ", idtmed[1106], tpar, 0); //Gas streamer |
1039 | |
1040 | // chamber type A |
1041 | tpar[0] = -1.; |
1042 | tpar[1] = -1.; |
1043 | |
1044 | const Float_t kXMC2A=kXMC1A*kZ12; |
1045 | const Float_t kYMC2Am=0.; |
1046 | const Float_t kYMC2Ap=0.; |
1047 | |
1048 | tpar[2] = 0.1; |
1049 | gMC->Gsposp("CG2A", 1, "CB2A", 0., 0., 0., 0, "ONLY",tpar,3); |
1050 | tpar[2] = 0.3; |
1051 | gMC->Gsposp("CB2A", 1, "CC2A", 0., 0., 0., 0, "ONLY",tpar,3); |
1052 | |
1053 | tpar[2] = 0.4; |
1054 | tpar[0] = ((kXMC1MAX-kXMC1MED)/2.)*kZ12; |
1055 | tpar[1] = kYMC1MIN*kZ12; |
1056 | |
1057 | gMC->Gsposp("CC2A", 1, "CM12",kXMC2A,kYMC2Am,kZMCm, 0, "ONLY", tpar, 3); |
1058 | gMC->Gsposp("CC2A", 2, "CM12",-kXMC2A,kYMC2Ap,kZMCp, 0, "ONLY", tpar, 3); |
1059 | |
1060 | |
1061 | // chamber type B |
1062 | |
1063 | tpar[0] = ((kXMC1MAX-kXMC1MIN)/2.)*kZ12; |
1064 | tpar[1] = ((kYMC1MAX-kYMC1MIN)/2.)*kZ12; |
1065 | |
1066 | const Float_t kXMC2B=kXMC1B*kZ12; |
1067 | const Float_t kYMC2Bp=kYMC1Bp*kZ12; |
1068 | const Float_t kYMC2Bm=kYMC1Bm*kZ12; |
1069 | gMC->Gsposp("CC2A", 3, "CM12",kXMC2B,kYMC2Bp,kZMCp, 0, "ONLY", tpar, 3); |
1070 | gMC->Gsposp("CC2A", 4, "CM12",-kXMC2B,kYMC2Bm,kZMCm, 0, "ONLY", tpar, 3); |
1071 | gMC->Gsposp("CC2A", 5, "CM12",kXMC2B,-kYMC2Bp,kZMCp, 0, "ONLY", tpar, 3); |
1072 | gMC->Gsposp("CC2A", 6, "CM12",-kXMC2B,-kYMC2Bm,kZMCm, 0, "ONLY", tpar, 3); |
1073 | |
1074 | |
1075 | // chamber type C (end of type B !!) |
1076 | |
1077 | tpar[0] = (kXMC1MAX/2)*kZ12; |
1078 | tpar[1] = (kYMC1MAX/2)*kZ12; |
1079 | |
1080 | const Float_t kXMC2C=kXMC1C*kZ12; |
1081 | const Float_t kYMC2Cp=kYMC1Cp*kZ12; |
1082 | const Float_t kYMC2Cm=kYMC1Cm*kZ12; |
1083 | gMC->Gsposp("CC2A", 7, "CM12",kXMC2C,kYMC2Cp,kZMCp, 0, "ONLY", tpar, 3); |
1084 | gMC->Gsposp("CC2A", 8, "CM12",-kXMC2C,kYMC2Cm,kZMCm, 0, "ONLY", tpar, 3); |
1085 | gMC->Gsposp("CC2A", 9, "CM12",kXMC2C,-kYMC2Cp,kZMCp, 0, "ONLY", tpar, 3); |
1086 | gMC->Gsposp("CC2A", 10, "CM12",-kXMC2C,-kYMC2Cm,kZMCm, 0, "ONLY", tpar, 3); |
1087 | |
1088 | // chamber type D, E and F (same size) |
1089 | |
1090 | tpar[0] = (kXMC1MAX/2.)*kZ12; |
1091 | tpar[1] = kYMC1MIN*kZ12; |
1092 | |
1093 | const Float_t kXMC2D=kXMC1D*kZ12; |
1094 | const Float_t kYMC2Dp=kYMC1Dp*kZ12; |
1095 | const Float_t kYMC2Dm=kYMC1Dm*kZ12; |
1096 | gMC->Gsposp("CC2A", 11, "CM12",kXMC2D,kYMC2Dm,kZMCm, 0, "ONLY", tpar, 3); |
1097 | gMC->Gsposp("CC2A", 12, "CM12",-kXMC2D,kYMC2Dp,kZMCp, 0, "ONLY", tpar, 3); |
1098 | gMC->Gsposp("CC2A", 13, "CM12",kXMC2D,-kYMC2Dm,kZMCm, 0, "ONLY", tpar, 3); |
1099 | gMC->Gsposp("CC2A", 14, "CM12",-kXMC2D,-kYMC2Dp,kZMCp, 0, "ONLY", tpar, 3); |
1100 | |
1101 | const Float_t kYMC2Ep=kYMC1Ep*kZ12; |
1102 | const Float_t kYMC2Em=kYMC1Em*kZ12; |
1103 | gMC->Gsposp("CC2A", 15, "CM12",kXMC2D,kYMC2Ep,kZMCp, 0, "ONLY", tpar, 3); |
1104 | gMC->Gsposp("CC2A", 16, "CM12",-kXMC2D,kYMC2Em,kZMCm, 0, "ONLY", tpar, 3); |
1105 | gMC->Gsposp("CC2A", 17, "CM12",kXMC2D,-kYMC2Ep,kZMCp, 0, "ONLY", tpar, 3); |
1106 | gMC->Gsposp("CC2A", 18, "CM12",-kXMC2D,-kYMC2Em,kZMCm, 0, "ONLY", tpar, 3); |
1107 | |
1108 | |
1109 | const Float_t kYMC2Fp=kYMC1Fp*kZ12; |
1110 | const Float_t kYMC2Fm=kYMC1Fm*kZ12; |
1111 | gMC->Gsposp("CC2A", 19, "CM12",kXMC2D,kYMC2Fm,kZMCm, 0, "ONLY", tpar, 3); |
1112 | gMC->Gsposp("CC2A", 20, "CM12",-kXMC2D,kYMC2Fp,kZMCp, 0, "ONLY", tpar, 3); |
1113 | gMC->Gsposp("CC2A", 21, "CM12",kXMC2D,-kYMC2Fm,kZMCm, 0, "ONLY", tpar, 3); |
1114 | gMC->Gsposp("CC2A", 22, "CM12",-kXMC2D,-kYMC2Fp,kZMCp, 0, "ONLY", tpar, 3); |
1115 | |
1116 | // Positioning second plane of station 1 in ALICE |
1117 | |
1118 | gMC->Gspos("CM12", 1, "ALIC", 0., 0., zpos2, 0, "ONLY"); |
1119 | |
1120 | // End of geometry definition for the second plane of station 1 |
1121 | |
1122 | |
1123 | |
1124 | // TRIGGER STATION 2 - TRIGGER STATION 2 - TRIGGER STATION 2 |
1125 | |
1126 | // 03/00 |
1127 | // zpos3 and zpos4 are now the middle of the first and second |
1128 | // plane of station 2 : |
1129 | // zpos3=(17075+16995)/2=17035 mm, thick/2=40 mm |
1130 | // zpos4=(17225+17145)/2=17185 mm, thick/2=40 mm |
1131 | // |
1132 | // zpos3m=16999 mm , zpos3p=17071 mm (middles of gas gaps) |
1133 | // zpos4m=17149 mm , zpos4p=17221 mm (middles of gas gaps) |
1134 | // rem : the total thickness accounts for 1 mm of al on both |
1135 | // side of the RPCs (see zpos3 and zpos4), as previously |
1136 | iChamber1 = iChamber = (AliMUONChamber*) (*fChambers)[12]; |
1137 | iChamber2 =(AliMUONChamber*) (*fChambers)[13]; |
1138 | Float_t zpos3=iChamber1->Z(); |
1139 | Float_t zpos4=iChamber2->Z(); |
1140 | |
1141 | |
1142 | // Mother volume definition |
1143 | tpar[0] = iChamber->RInner(); |
1144 | tpar[1] = iChamber->ROuter(); |
1145 | tpar[2] = 4.0; |
1146 | |
1147 | gMC->Gsvolu("CM21", "TUBE", idAir, tpar, 3); |
1148 | gMC->Gsvolu("CM22", "TUBE", idAir, tpar, 3); |
1149 | |
1150 | // Definition of the flange between the beam shielding and the RPC |
1151 | // ???? interface shielding |
1152 | |
1153 | tpar[0]= kRMIN2; |
1154 | tpar[1]= kRMAX2; |
1155 | tpar[2]= 4.0; |
1156 | |
1157 | gMC->Gsvolu("CF2A", "TUBE", idAlu1, tpar, 3); //Al |
1158 | gMC->Gspos("CF2A", 1, "CM21", 0., 0., 0., 0, "MANY"); |
1159 | gMC->Gspos("CF2A", 2, "CM22", 0., 0., 0., 0, "MANY"); |
1160 | |
1161 | |
1162 | |
1163 | // FIRST PLANE OF STATION 2 : proj ratio = zpos3/zpos1 |
1164 | |
1165 | const Float_t kZ13=zpos3/zpos1; |
1166 | |
1167 | // Definition of prototype for chambers in the first plane of station 2 |
1168 | tpar[0]= 0.; |
1169 | tpar[1]= 0.; |
1170 | tpar[2]= 0.; |
1171 | |
1172 | gMC->Gsvolu("CC3A", "BOX ", idAlu1, tpar, 0); //Al |
1173 | gMC->Gsvolu("CB3A", "BOX ", idtmed[1107], tpar, 0); //Bakelite |
1174 | gMC->Gsvolu("CG3A", "BOX ", idtmed[1106], tpar, 0); //Gas streamer |
1175 | |
1176 | |
1177 | // chamber type A |
1178 | tpar[0] = -1.; |
1179 | tpar[1] = -1.; |
1180 | |
1181 | const Float_t kXMC3A=kXMC1A*kZ13; |
1182 | const Float_t kYMC3Am=0.; |
1183 | const Float_t kYMC3Ap=0.; |
1184 | |
1185 | tpar[2] = 0.1; |
1186 | gMC->Gsposp("CG3A", 1, "CB3A", 0., 0., 0., 0, "ONLY",tpar,3); |
1187 | tpar[2] = 0.3; |
1188 | gMC->Gsposp("CB3A", 1, "CC3A", 0., 0., 0., 0, "ONLY",tpar,3); |
1189 | |
1190 | tpar[2] = 0.4; |
1191 | tpar[0] = ((kXMC1MAX-kXMC1MED)/2.)*kZ13; |
1192 | tpar[1] = kYMC1MIN*kZ13; |
1193 | gMC->Gsposp("CC3A", 1, "CM21",kXMC3A,kYMC3Am,kZMCm, 0, "ONLY", tpar, 3); |
1194 | gMC->Gsposp("CC3A", 2, "CM21",-kXMC3A,kYMC3Ap,kZMCp, 0, "ONLY", tpar, 3); |
1195 | |
1196 | |
1197 | // chamber type B |
1198 | tpar[0] = ((kXMC1MAX-kXMC1MIN)/2.)*kZ13; |
1199 | tpar[1] = ((kYMC1MAX-kYMC1MIN)/2.)*kZ13; |
1200 | |
1201 | const Float_t kXMC3B=kXMC1B*kZ13; |
1202 | const Float_t kYMC3Bp=kYMC1Bp*kZ13; |
1203 | const Float_t kYMC3Bm=kYMC1Bm*kZ13; |
1204 | gMC->Gsposp("CC3A", 3, "CM21",kXMC3B,kYMC3Bp,kZMCp, 0, "ONLY", tpar, 3); |
1205 | gMC->Gsposp("CC3A", 4, "CM21",-kXMC3B,kYMC3Bm,kZMCm, 0, "ONLY", tpar, 3); |
1206 | gMC->Gsposp("CC3A", 5, "CM21",kXMC3B,-kYMC3Bp,kZMCp, 0, "ONLY", tpar, 3); |
1207 | gMC->Gsposp("CC3A", 6, "CM21",-kXMC3B,-kYMC3Bm,kZMCm, 0, "ONLY", tpar, 3); |
1208 | |
1209 | |
1210 | // chamber type C (end of type B !!) |
1211 | tpar[0] = (kXMC1MAX/2)*kZ13; |
1212 | tpar[1] = (kYMC1MAX/2)*kZ13; |
1213 | |
1214 | const Float_t kXMC3C=kXMC1C*kZ13; |
1215 | const Float_t kYMC3Cp=kYMC1Cp*kZ13; |
1216 | const Float_t kYMC3Cm=kYMC1Cm*kZ13; |
1217 | gMC->Gsposp("CC3A", 7, "CM21",kXMC3C,kYMC3Cp,kZMCp, 0, "ONLY", tpar, 3); |
1218 | gMC->Gsposp("CC3A", 8, "CM21",-kXMC3C,kYMC3Cm,kZMCm, 0, "ONLY", tpar, 3); |
1219 | gMC->Gsposp("CC3A", 9, "CM21",kXMC3C,-kYMC3Cp,kZMCp, 0, "ONLY", tpar, 3); |
1220 | gMC->Gsposp("CC3A", 10, "CM21",-kXMC3C,-kYMC3Cm,kZMCm, 0, "ONLY", tpar, 3); |
1221 | |
1222 | |
1223 | // chamber type D, E and F (same size) |
1224 | |
1225 | tpar[0] = (kXMC1MAX/2.)*kZ13; |
1226 | tpar[1] = kYMC1MIN*kZ13; |
1227 | |
1228 | const Float_t kXMC3D=kXMC1D*kZ13; |
1229 | const Float_t kYMC3Dp=kYMC1Dp*kZ13; |
1230 | const Float_t kYMC3Dm=kYMC1Dm*kZ13; |
1231 | gMC->Gsposp("CC3A", 11, "CM21",kXMC3D,kYMC3Dm,kZMCm, 0, "ONLY", tpar, 3); |
1232 | gMC->Gsposp("CC3A", 12, "CM21",-kXMC3D,kYMC3Dp,kZMCp, 0, "ONLY", tpar, 3); |
1233 | gMC->Gsposp("CC3A", 13, "CM21",kXMC3D,-kYMC3Dm,kZMCm, 0, "ONLY", tpar, 3); |
1234 | gMC->Gsposp("CC3A", 14, "CM21",-kXMC3D,-kYMC3Dp,kZMCp, 0, "ONLY", tpar, 3); |
1235 | |
1236 | const Float_t kYMC3Ep=kYMC1Ep*kZ13; |
1237 | const Float_t kYMC3Em=kYMC1Em*kZ13; |
1238 | gMC->Gsposp("CC3A", 15, "CM21",kXMC3D,kYMC3Ep,kZMCp, 0, "ONLY", tpar, 3); |
1239 | gMC->Gsposp("CC3A", 16, "CM21",-kXMC3D,kYMC3Em,kZMCm, 0, "ONLY", tpar, 3); |
1240 | gMC->Gsposp("CC3A", 17, "CM21",kXMC3D,-kYMC3Ep,kZMCp, 0, "ONLY", tpar, 3); |
1241 | gMC->Gsposp("CC3A", 18, "CM21",-kXMC3D,-kYMC3Em,kZMCm, 0, "ONLY", tpar, 3); |
1242 | |
1243 | const Float_t kYMC3Fp=kYMC1Fp*kZ13; |
1244 | const Float_t kYMC3Fm=kYMC1Fm*kZ13; |
1245 | gMC->Gsposp("CC3A", 19, "CM21",kXMC3D,kYMC3Fm,kZMCm, 0, "ONLY", tpar, 3); |
1246 | gMC->Gsposp("CC3A", 20, "CM21",-kXMC3D,kYMC3Fp,kZMCp, 0, "ONLY", tpar, 3); |
1247 | gMC->Gsposp("CC3A", 21, "CM21",kXMC3D,-kYMC3Fm,kZMCm, 0, "ONLY", tpar, 3); |
1248 | gMC->Gsposp("CC3A", 22, "CM21",-kXMC3D,-kYMC3Fp,kZMCp, 0, "ONLY", tpar, 3); |
1249 | |
1250 | |
1251 | // Positioning first plane of station 2 in ALICE |
1252 | |
1253 | gMC->Gspos("CM21", 1, "ALIC", 0., 0., zpos3, 0, "ONLY"); |
1254 | |
1255 | // End of geometry definition for the first plane of station 2 |
1256 | |
1257 | |
1258 | |
1259 | |
1260 | // SECOND PLANE OF STATION 2 : proj ratio = zpos4/zpos1 |
1261 | |
1262 | const Float_t kZ14=zpos4/zpos1; |
1263 | |
1264 | // Definition of prototype for chambers in the second plane of station 2 |
1265 | |
1266 | tpar[0]= 0.; |
1267 | tpar[1]= 0.; |
1268 | tpar[2]= 0.; |
1269 | |
1270 | gMC->Gsvolu("CC4A", "BOX ", idAlu1, tpar, 0); //Al |
1271 | gMC->Gsvolu("CB4A", "BOX ", idtmed[1107], tpar, 0); //Bakelite |
1272 | gMC->Gsvolu("CG4A", "BOX ", idtmed[1106], tpar, 0); //Gas streamer |
1273 | |
1274 | // chamber type A |
1275 | tpar[0] = -1.; |
1276 | tpar[1] = -1.; |
1277 | |
1278 | const Float_t kXMC4A=kXMC1A*kZ14; |
1279 | const Float_t kYMC4Am=0.; |
1280 | const Float_t kYMC4Ap=0.; |
1281 | |
1282 | tpar[2] = 0.1; |
1283 | gMC->Gsposp("CG4A", 1, "CB4A", 0., 0., 0., 0, "ONLY",tpar,3); |
1284 | tpar[2] = 0.3; |
1285 | gMC->Gsposp("CB4A", 1, "CC4A", 0., 0., 0., 0, "ONLY",tpar,3); |
1286 | |
1287 | tpar[2] = 0.4; |
1288 | tpar[0] = ((kXMC1MAX-kXMC1MED)/2.)*kZ14; |
1289 | tpar[1] = kYMC1MIN*kZ14; |
1290 | gMC->Gsposp("CC4A", 1, "CM22",kXMC4A,kYMC4Am,kZMCm, 0, "ONLY", tpar, 3); |
1291 | gMC->Gsposp("CC4A", 2, "CM22",-kXMC4A,kYMC4Ap,kZMCp, 0, "ONLY", tpar, 3); |
1292 | |
1293 | |
1294 | // chamber type B |
1295 | tpar[0] = ((kXMC1MAX-kXMC1MIN)/2.)*kZ14; |
1296 | tpar[1] = ((kYMC1MAX-kYMC1MIN)/2.)*kZ14; |
1297 | |
1298 | const Float_t kXMC4B=kXMC1B*kZ14; |
1299 | const Float_t kYMC4Bp=kYMC1Bp*kZ14; |
1300 | const Float_t kYMC4Bm=kYMC1Bm*kZ14; |
1301 | gMC->Gsposp("CC4A", 3, "CM22",kXMC4B,kYMC4Bp,kZMCp, 0, "ONLY", tpar, 3); |
1302 | gMC->Gsposp("CC4A", 4, "CM22",-kXMC4B,kYMC4Bm,kZMCm, 0, "ONLY", tpar, 3); |
1303 | gMC->Gsposp("CC4A", 5, "CM22",kXMC4B,-kYMC4Bp,kZMCp, 0, "ONLY", tpar, 3); |
1304 | gMC->Gsposp("CC4A", 6, "CM22",-kXMC4B,-kYMC4Bm,kZMCm, 0, "ONLY", tpar, 3); |
1305 | |
1306 | |
1307 | // chamber type C (end of type B !!) |
1308 | tpar[0] =(kXMC1MAX/2)*kZ14; |
1309 | tpar[1] = (kYMC1MAX/2)*kZ14; |
1310 | |
1311 | const Float_t kXMC4C=kXMC1C*kZ14; |
1312 | const Float_t kYMC4Cp=kYMC1Cp*kZ14; |
1313 | const Float_t kYMC4Cm=kYMC1Cm*kZ14; |
1314 | gMC->Gsposp("CC4A", 7, "CM22",kXMC4C,kYMC4Cp,kZMCp, 0, "ONLY", tpar, 3); |
1315 | gMC->Gsposp("CC4A", 8, "CM22",-kXMC4C,kYMC4Cm,kZMCm, 0, "ONLY", tpar, 3); |
1316 | gMC->Gsposp("CC4A", 9, "CM22",kXMC4C,-kYMC4Cp,kZMCp, 0, "ONLY", tpar, 3); |
1317 | gMC->Gsposp("CC4A", 10, "CM22",-kXMC4C,-kYMC4Cm,kZMCm, 0, "ONLY", tpar, 3); |
1318 | |
1319 | |
1320 | // chamber type D, E and F (same size) |
1321 | tpar[0] = (kXMC1MAX/2.)*kZ14; |
1322 | tpar[1] = kYMC1MIN*kZ14; |
1323 | |
1324 | const Float_t kXMC4D=kXMC1D*kZ14; |
1325 | const Float_t kYMC4Dp=kYMC1Dp*kZ14; |
1326 | const Float_t kYMC4Dm=kYMC1Dm*kZ14; |
1327 | gMC->Gsposp("CC4A", 11, "CM22",kXMC4D,kYMC4Dm,kZMCm, 0, "ONLY", tpar, 3); |
1328 | gMC->Gsposp("CC4A", 12, "CM22",-kXMC4D,kYMC4Dp,kZMCp, 0, "ONLY", tpar, 3); |
1329 | gMC->Gsposp("CC4A", 13, "CM22",kXMC4D,-kYMC4Dm,kZMCm, 0, "ONLY", tpar, 3); |
1330 | gMC->Gsposp("CC4A", 14, "CM22",-kXMC4D,-kYMC4Dp,kZMCp, 0, "ONLY", tpar, 3); |
1331 | |
1332 | const Float_t kYMC4Ep=kYMC1Ep*kZ14; |
1333 | const Float_t kYMC4Em=kYMC1Em*kZ14; |
1334 | gMC->Gsposp("CC4A", 15, "CM22",kXMC4D,kYMC4Ep,kZMCp, 0, "ONLY", tpar, 3); |
1335 | gMC->Gsposp("CC4A", 16, "CM22",-kXMC4D,kYMC4Em,kZMCm, 0, "ONLY", tpar, 3); |
1336 | gMC->Gsposp("CC4A", 17, "CM22",kXMC4D,-kYMC4Ep,kZMCp, 0, "ONLY", tpar, 3); |
1337 | gMC->Gsposp("CC4A", 18, "CM22",-kXMC4D,-kYMC4Em,kZMCm, 0, "ONLY", tpar, 3); |
1338 | |
1339 | const Float_t kYMC4Fp=kYMC1Fp*kZ14; |
1340 | const Float_t kYMC4Fm=kYMC1Fm*kZ14; |
1341 | gMC->Gsposp("CC4A", 19, "CM22",kXMC4D,kYMC4Fm,kZMCm, 0, "ONLY", tpar, 3); |
1342 | gMC->Gsposp("CC4A", 20, "CM22",-kXMC4D,kYMC4Fp,kZMCp, 0, "ONLY", tpar, 3); |
1343 | gMC->Gsposp("CC4A", 21, "CM22",kXMC4D,-kYMC4Fm,kZMCm, 0, "ONLY", tpar, 3); |
1344 | gMC->Gsposp("CC4A", 22, "CM22",-kXMC4D,-kYMC4Fp,kZMCp, 0, "ONLY", tpar, 3); |
1345 | |
1346 | |
1347 | // Positioning second plane of station 2 in ALICE |
1348 | |
1349 | gMC->Gspos("CM22", 1, "ALIC", 0., 0., zpos4, 0, "ONLY"); |
1350 | |
1351 | // End of geometry definition for the second plane of station 2 |
1352 | |
1353 | // End of trigger geometry definition |
1354 | |
1355 | } |
1356 | |
1357 | |
1358 | |
1359 | //___________________________________________ |
1360 | void AliMUONv1::CreateMaterials() |
1361 | { |
1362 | // *** DEFINITION OF AVAILABLE MUON MATERIALS *** |
1363 | // |
1364 | // Ar-CO2 gas |
1365 | Float_t ag1[3] = { 39.95,12.01,16. }; |
1366 | Float_t zg1[3] = { 18.,6.,8. }; |
1367 | Float_t wg1[3] = { .8,.0667,.13333 }; |
1368 | Float_t dg1 = .001821; |
1369 | // |
1370 | // Ar-buthane-freon gas -- trigger chambers |
1371 | Float_t atr1[4] = { 39.95,12.01,1.01,19. }; |
1372 | Float_t ztr1[4] = { 18.,6.,1.,9. }; |
1373 | Float_t wtr1[4] = { .56,.1262857,.2857143,.028 }; |
1374 | Float_t dtr1 = .002599; |
1375 | // |
1376 | // Ar-CO2 gas |
1377 | Float_t agas[3] = { 39.95,12.01,16. }; |
1378 | Float_t zgas[3] = { 18.,6.,8. }; |
1379 | Float_t wgas[3] = { .74,.086684,.173316 }; |
1380 | Float_t dgas = .0018327; |
1381 | // |
1382 | // Ar-Isobutane gas (80%+20%) -- tracking |
1383 | Float_t ag[3] = { 39.95,12.01,1.01 }; |
1384 | Float_t zg[3] = { 18.,6.,1. }; |
1385 | Float_t wg[3] = { .8,.057,.143 }; |
1386 | Float_t dg = .0019596; |
1387 | // |
1388 | // Ar-Isobutane-Forane-SF6 gas (49%+7%+40%+4%) -- trigger |
1389 | Float_t atrig[5] = { 39.95,12.01,1.01,19.,32.066 }; |
1390 | Float_t ztrig[5] = { 18.,6.,1.,9.,16. }; |
1391 | Float_t wtrig[5] = { .49,1.08,1.5,1.84,0.04 }; |
1392 | Float_t dtrig = .0031463; |
1393 | // |
1394 | // bakelite |
1395 | |
1396 | Float_t abak[3] = {12.01 , 1.01 , 16.}; |
1397 | Float_t zbak[3] = {6. , 1. , 8.}; |
1398 | Float_t wbak[3] = {6. , 6. , 1.}; |
1399 | Float_t dbak = 1.4; |
1400 | |
1401 | Float_t epsil, stmin, deemax, tmaxfd, stemax; |
1402 | |
1403 | Int_t iSXFLD = gAlice->Field()->Integ(); |
1404 | Float_t sXMGMX = gAlice->Field()->Max(); |
1405 | // |
1406 | // --- Define the various materials for GEANT --- |
1407 | AliMaterial(9, "ALUMINIUM$", 26.98, 13., 2.7, 8.9, 37.2); |
1408 | AliMaterial(10, "ALUMINIUM$", 26.98, 13., 2.7, 8.9, 37.2); |
1409 | AliMaterial(15, "AIR$ ", 14.61, 7.3, .001205, 30423.24, 67500); |
1410 | AliMixture(19, "Bakelite$", abak, zbak, dbak, -3, wbak); |
1411 | AliMixture(20, "ArC4H10 GAS$", ag, zg, dg, 3, wg); |
1412 | AliMixture(21, "TRIG GAS$", atrig, ztrig, dtrig, -5, wtrig); |
1413 | AliMixture(22, "ArCO2 80%$", ag1, zg1, dg1, 3, wg1); |
1414 | AliMixture(23, "Ar-freon $", atr1, ztr1, dtr1, 4, wtr1); |
1415 | AliMixture(24, "ArCO2 GAS$", agas, zgas, dgas, 3, wgas); |
1416 | |
1417 | epsil = .001; // Tracking precision, |
1418 | stemax = -1.; // Maximum displacement for multiple scat |
1419 | tmaxfd = -20.; // Maximum angle due to field deflection |
1420 | deemax = -.3; // Maximum fractional energy loss, DLS |
1421 | stmin = -.8; |
1422 | // |
1423 | // Air |
1424 | AliMedium(1, "AIR_CH_US ", 15, 1, iSXFLD, sXMGMX, tmaxfd, stemax, deemax, epsil, stmin); |
1425 | // |
1426 | // Aluminum |
1427 | |
1428 | AliMedium(4, "ALU_CH_US ", 9, 0, iSXFLD, sXMGMX, tmaxfd, fMaxStepAlu, |
1429 | fMaxDestepAlu, epsil, stmin); |
1430 | AliMedium(5, "ALU_CH_US ", 10, 0, iSXFLD, sXMGMX, tmaxfd, fMaxStepAlu, |
1431 | fMaxDestepAlu, epsil, stmin); |
1432 | // |
1433 | // Ar-isoC4H10 gas |
1434 | |
1435 | AliMedium(6, "AR_CH_US ", 20, 1, iSXFLD, sXMGMX, tmaxfd, fMaxStepGas, |
1436 | fMaxDestepGas, epsil, stmin); |
1437 | // |
1438 | // Ar-Isobuthane-Forane-SF6 gas |
1439 | |
1440 | AliMedium(7, "GAS_CH_TRIGGER ", 21, 1, iSXFLD, sXMGMX, tmaxfd, stemax, deemax, epsil, stmin); |
1441 | |
1442 | AliMedium(8, "BAKE_CH_TRIGGER ", 19, 0, iSXFLD, sXMGMX, tmaxfd, fMaxStepAlu, |
1443 | fMaxDestepAlu, epsil, stmin); |
1444 | |
1445 | AliMedium(9, "ARG_CO2 ", 22, 1, iSXFLD, sXMGMX, tmaxfd, fMaxStepGas, |
1446 | fMaxDestepAlu, epsil, stmin); |
1447 | } |
1448 | |
1449 | //___________________________________________ |
1450 | |
1451 | void AliMUONv1::Init() |
1452 | { |
1453 | // |
1454 | // Initialize Tracking Chambers |
1455 | // |
1456 | |
1457 | printf("\n\n\n Start Init for version 1 - CPC chamber type\n\n\n"); |
1458 | |
1459 | for (Int_t i=0; i<kNCH; i++) { |
1460 | ( (AliMUONChamber*) (*fChambers)[i])->Init(); |
1461 | } |
1462 | |
1463 | // |
1464 | // Set the chamber (sensitive region) GEANT identifier |
1465 | AliMC* gMC = AliMC::GetMC(); |
1466 | ((AliMUONChamber*)(*fChambers)[0])->SetGid(gMC->VolId("C01G")); |
1467 | ((AliMUONChamber*)(*fChambers)[1])->SetGid(gMC->VolId("C02G")); |
1468 | ((AliMUONChamber*)(*fChambers)[2])->SetGid(gMC->VolId("C03G")); |
1469 | ((AliMUONChamber*)(*fChambers)[3])->SetGid(gMC->VolId("C04G")); |
1470 | ((AliMUONChamber*)(*fChambers)[4])->SetGid(gMC->VolId("C05G")); |
1471 | ((AliMUONChamber*)(*fChambers)[5])->SetGid(gMC->VolId("C06G")); |
1472 | ((AliMUONChamber*)(*fChambers)[6])->SetGid(gMC->VolId("C07G")); |
1473 | ((AliMUONChamber*)(*fChambers)[7])->SetGid(gMC->VolId("C08G")); |
1474 | ((AliMUONChamber*)(*fChambers)[8])->SetGid(gMC->VolId("C09G")); |
1475 | ((AliMUONChamber*)(*fChambers)[9])->SetGid(gMC->VolId("C10G")); |
1476 | ((AliMUONChamber*)(*fChambers)[10])->SetGid(gMC->VolId("CG1A")); |
1477 | ((AliMUONChamber*)(*fChambers)[11])->SetGid(gMC->VolId("CG2A")); |
1478 | ((AliMUONChamber*)(*fChambers)[12])->SetGid(gMC->VolId("CG3A")); |
1479 | ((AliMUONChamber*)(*fChambers)[13])->SetGid(gMC->VolId("CG4A")); |
1480 | |
1481 | printf("\n\n\n Finished Init for version 0 - CPC chamber type\n\n\n"); |
1482 | |
1483 | //cp |
1484 | printf("\n\n\n Start Init for Trigger Circuits\n\n\n"); |
1485 | for (Int_t i=0; i<kNTriggerCircuit; i++) { |
1486 | ( (AliMUONTriggerCircuit*) (*fTriggerCircuits)[i])->Init(i); |
1487 | } |
1488 | printf(" Finished Init for Trigger Circuits\n\n\n"); |
1489 | //cp |
1490 | |
1491 | } |
1492 | |
1493 | //___________________________________________ |
1494 | void AliMUONv1::StepManager() |
1495 | { |
1496 | Int_t copy, id; |
1497 | static Int_t idvol; |
1498 | static Int_t vol[2]; |
1499 | Int_t ipart; |
1500 | TLorentzVector pos; |
1501 | TLorentzVector mom; |
1502 | Float_t theta,phi; |
1503 | Float_t destep, step; |
1504 | |
1505 | static Float_t eloss, eloss2, xhit, yhit, tof, tlength; |
1506 | const Float_t kBig=1.e10; |
1507 | |
1508 | // modifs perso |
1509 | static Float_t hits[15]; |
1510 | |
1511 | TClonesArray &lhits = *fHits; |
1512 | |
1513 | // |
1514 | // Set maximum step size for gas |
1515 | // numed=gMC->GetMedium(); |
1516 | // |
1517 | // Only charged tracks |
1518 | if( !(gMC->TrackCharge()) ) return; |
1519 | // |
1520 | // Only gas gap inside chamber |
1521 | // Tag chambers and record hits when track enters |
1522 | idvol=-1; |
1523 | id=gMC->CurrentVolID(copy); |
1524 | |
1525 | for (Int_t i=1; i<=kNCH; i++) { |
1526 | if(id==((AliMUONChamber*)(*fChambers)[i-1])->GetGid()){ |
1527 | vol[0]=i; |
1528 | idvol=i-1; |
1529 | } |
1530 | } |
1531 | if (idvol == -1) return; |
1532 | // |
1533 | // Get current particle id (ipart), track position (pos) and momentum (mom) |
1534 | gMC->TrackPosition(pos); |
1535 | gMC->TrackMomentum(mom); |
1536 | |
1537 | ipart = gMC->TrackPid(); |
1538 | //Int_t ipart1 = gMC->IdFromPDG(ipart); |
1539 | //printf("ich, ipart %d %d \n",vol[0],ipart1); |
1540 | |
1541 | // |
1542 | // momentum loss and steplength in last step |
1543 | destep = gMC->Edep(); |
1544 | step = gMC->TrackStep(); |
1545 | |
1546 | // |
1547 | // record hits when track enters ... |
1548 | if( gMC->IsTrackEntering()) { |
1549 | gMC->SetMaxStep(fMaxStepGas); |
1550 | Double_t tc = mom[0]*mom[0]+mom[1]*mom[1]; |
1551 | Double_t rt = TMath::Sqrt(tc); |
1552 | Double_t pmom = TMath::Sqrt(tc+mom[2]*mom[2]); |
1553 | Double_t tx=mom[0]/pmom; |
1554 | Double_t ty=mom[1]/pmom; |
1555 | Double_t tz=mom[2]/pmom; |
1556 | Double_t s=((AliMUONChamber*)(*fChambers)[idvol]) |
1557 | ->ResponseModel() |
1558 | ->Pitch()/tz; |
1559 | theta = Float_t(TMath::ATan2(rt,Double_t(mom[2])))*kRaddeg; |
1560 | phi = Float_t(TMath::ATan2(Double_t(mom[1]),Double_t(mom[0])))*kRaddeg; |
1561 | hits[0] = Float_t(ipart); // Geant3 particle type |
1562 | hits[1] = pos[0]+s*tx; // X-position for hit |
1563 | hits[2] = pos[1]+s*ty; // Y-position for hit |
1564 | hits[3] = pos[2]+s*tz; // Z-position for hit |
1565 | hits[4] = theta; // theta angle of incidence |
1566 | hits[5] = phi; // phi angle of incidence |
1567 | hits[8] = (Float_t) fNPadHits; // first padhit |
1568 | hits[9] = -1; // last pad hit |
1569 | |
1570 | // modifs perso |
1571 | hits[10] = mom[3]; // hit momentum P |
1572 | hits[11] = mom[0]; // Px/P |
1573 | hits[12] = mom[1]; // Py/P |
1574 | hits[13] = mom[2]; // Pz/P |
1575 | // fin modifs perso |
1576 | tof=gMC->TrackTime(); |
1577 | hits[14] = tof; // Time of flight |
1578 | // phi angle of incidence |
1579 | tlength = 0; |
1580 | eloss = 0; |
1581 | eloss2 = 0; |
1582 | xhit = pos[0]; |
1583 | yhit = pos[1]; |
1584 | // Only if not trigger chamber |
1585 | if(idvol<10) { |
1586 | // |
1587 | // Initialize hit position (cursor) in the segmentation model |
1588 | ((AliMUONChamber*) (*fChambers)[idvol]) |
1589 | ->SigGenInit(pos[0], pos[1], pos[2]); |
1590 | } else { |
1591 | //geant3->Gpcxyz(); |
1592 | //printf("In the Trigger Chamber #%d\n",idvol-9); |
1593 | } |
1594 | } |
1595 | eloss2+=destep; |
1596 | |
1597 | // |
1598 | // Calculate the charge induced on a pad (disintegration) in case |
1599 | // |
1600 | // Mip left chamber ... |
1601 | if( gMC->IsTrackExiting() || gMC->IsTrackStop() || gMC->IsTrackDisappeared()){ |
1602 | gMC->SetMaxStep(kBig); |
1603 | eloss += destep; |
1604 | tlength += step; |
1605 | |
1606 | Float_t x0,y0; |
1607 | |
1608 | if(idvol<10) { |
1609 | // tracking chambers |
1610 | x0 = 0.5*(xhit+pos[0]); |
1611 | y0 = 0.5*(yhit+pos[1]); |
1612 | } else { |
1613 | // trigger chambers |
1614 | x0=xhit; |
1615 | y0=yhit; |
1616 | } |
1617 | |
1618 | |
1619 | if (eloss >0) MakePadHits(x0,y0,eloss,tof,idvol); |
1620 | |
1621 | |
1622 | hits[6]=tlength; |
1623 | hits[7]=eloss2; |
1624 | if (fNPadHits > (Int_t)hits[8]) { |
1625 | hits[8]= hits[8]+1; |
1626 | hits[9]= (Float_t) fNPadHits; |
1627 | } |
1628 | |
1629 | new(lhits[fNhits++]) |
1630 | AliMUONHit(fIshunt,gAlice->CurrentTrack(),vol,hits); |
1631 | eloss = 0; |
1632 | // |
1633 | // Check additional signal generation conditions |
1634 | // defined by the segmentation |
1635 | // model (boundary crossing conditions) |
1636 | } else if |
1637 | (((AliMUONChamber*) (*fChambers)[idvol]) |
1638 | ->SigGenCond(pos[0], pos[1], pos[2])) |
1639 | { |
1640 | ((AliMUONChamber*) (*fChambers)[idvol]) |
1641 | ->SigGenInit(pos[0], pos[1], pos[2]); |
1642 | // printf("\n-> MakePadHits, reason special %d",ipart); |
1643 | if (eloss > 0 && idvol < 10) |
1644 | MakePadHits(0.5*(xhit+pos[0]),0.5*(yhit+pos[1]),eloss,tof,idvol); |
1645 | xhit = pos[0]; |
1646 | yhit = pos[1]; |
1647 | eloss = destep; |
1648 | tlength += step ; |
1649 | // |
1650 | // nothing special happened, add up energy loss |
1651 | } else { |
1652 | eloss += destep; |
1653 | tlength += step ; |
1654 | } |
1655 | } |
1656 | |
1657 | |