]>
Commit | Line | Data |
---|---|---|
f7336fa3 | 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 | ||
afc51ac2 | 16 | /* $Id$ */ |
f7336fa3 | 17 | |
18 | /////////////////////////////////////////////////////////////////////////////// | |
19 | // // | |
20 | // TRD geometry class // | |
21 | // // | |
22 | /////////////////////////////////////////////////////////////////////////////// | |
23 | ||
793ff80c | 24 | |
bdbb05bb | 25 | #include <TError.h> |
b4a9cd27 | 26 | #include <TGeoManager.h> |
27 | #include <TGeoPhysicalNode.h> | |
28 | #include <TGeoMatrix.h> | |
29 | ||
bdbb05bb | 30 | |
31 | #include "AliRunLoader.h" | |
f7336fa3 | 32 | #include "AliTRDgeometry.h" |
a5cadd36 | 33 | #include "AliTRDpadPlane.h" |
f7336fa3 | 34 | |
b4a9cd27 | 35 | #include "AliAlignObj.h" |
36 | #include "AliAlignObjAngles.h" | |
37 | ||
ecb36af7 | 38 | #include "AliRun.h" |
39 | #include "AliTRD.h" | |
3551db50 | 40 | #include "AliTRDcalibDB.h" |
41 | #include "AliTRDCommonParam.h" | |
ecb36af7 | 42 | |
f7336fa3 | 43 | ClassImp(AliTRDgeometry) |
44 | ||
793ff80c | 45 | //_____________________________________________________________________________ |
46 | ||
47 | // | |
48 | // The geometry constants | |
49 | // | |
50 | const Int_t AliTRDgeometry::fgkNsect = kNsect; | |
51 | const Int_t AliTRDgeometry::fgkNplan = kNplan; | |
52 | const Int_t AliTRDgeometry::fgkNcham = kNcham; | |
53 | const Int_t AliTRDgeometry::fgkNdet = kNdet; | |
54 | ||
55 | // | |
56 | // Dimensions of the detector | |
57 | // | |
0a770ac9 | 58 | |
59 | // Inner and outer radius of the mother volumes | |
793ff80c | 60 | const Float_t AliTRDgeometry::fgkRmin = 294.0; |
61 | const Float_t AliTRDgeometry::fgkRmax = 368.0; | |
62 | ||
0a770ac9 | 63 | // Upper and lower length of the mother volumes |
793ff80c | 64 | const Float_t AliTRDgeometry::fgkZmax1 = 378.35; |
65 | const Float_t AliTRDgeometry::fgkZmax2 = 302.0; | |
66 | ||
0a770ac9 | 67 | // Parameter of the BTR mother volumes |
bd0f8685 | 68 | // CBL |
69 | //const Float_t AliTRDgeometry::fgkSheight = 74.0; | |
70 | const Float_t AliTRDgeometry::fgkSheight = 74.86; | |
71 | const Float_t AliTRDgeometry::fgkSwidth1 = 99.613; | |
72 | // CBL | |
73 | //const Float_t AliTRDgeometry::fgkSwidth2 = 125.707; | |
74 | const Float_t AliTRDgeometry::fgkSwidth2 = 126.012; | |
793ff80c | 75 | const Float_t AliTRDgeometry::fgkSlenTR1 = 751.0; |
76 | const Float_t AliTRDgeometry::fgkSlenTR2 = 313.5; | |
77 | const Float_t AliTRDgeometry::fgkSlenTR3 = 159.5; | |
78 | ||
73ae7b59 | 79 | // The super module side plates |
287c5d50 | 80 | const Float_t AliTRDgeometry::fgkSMpltT = 0.2; |
81 | const Float_t AliTRDgeometry::fgkSMgapT = 0.5; | |
73ae7b59 | 82 | |
0a770ac9 | 83 | // Height of different chamber parts |
84 | // Radiator | |
85 | const Float_t AliTRDgeometry::fgkCraH = 4.8; | |
86 | // Drift region | |
87 | const Float_t AliTRDgeometry::fgkCdrH = 3.0; | |
88 | // Amplification region | |
89 | const Float_t AliTRDgeometry::fgkCamH = 0.7; | |
90 | // Readout | |
73ae7b59 | 91 | const Float_t AliTRDgeometry::fgkCroH = 2.316; |
0a770ac9 | 92 | // Total height |
93 | const Float_t AliTRDgeometry::fgkCH = AliTRDgeometry::fgkCraH | |
94 | + AliTRDgeometry::fgkCdrH | |
95 | + AliTRDgeometry::fgkCamH | |
96 | + AliTRDgeometry::fgkCroH; | |
97 | ||
98 | // Vertical spacing of the chambers | |
73ae7b59 | 99 | const Float_t AliTRDgeometry::fgkVspace = 1.784; |
0a770ac9 | 100 | |
101 | // Horizontal spacing of the chambers | |
102 | const Float_t AliTRDgeometry::fgkHspace = 2.0; | |
103 | ||
104 | // Thicknesses of different parts of the chamber frame | |
105 | // Lower aluminum frame | |
106 | const Float_t AliTRDgeometry::fgkCalT = 0.3; | |
107 | // Lower G10 frame sides | |
108 | const Float_t AliTRDgeometry::fgkCclsT = 0.3; | |
109 | // Lower G10 frame front | |
110 | const Float_t AliTRDgeometry::fgkCclfT = 1.0; | |
111 | // Upper G10 frame | |
112 | const Float_t AliTRDgeometry::fgkCcuT = 0.9; | |
113 | // Upper Al frame | |
114 | const Float_t AliTRDgeometry::fgkCauT = 1.5; | |
115 | ||
116 | // Additional width of the readout chamber frames | |
117 | const Float_t AliTRDgeometry::fgkCroW = 0.9; | |
118 | ||
119 | // Difference of outer chamber width and pad plane width | |
73ae7b59 | 120 | //const Float_t AliTRDgeometry::fgkCpadW = 1.0; |
121 | const Float_t AliTRDgeometry::fgkCpadW = 0.0; | |
de6df1b1 | 122 | const Float_t AliTRDgeometry::fgkRpadW = 1.0; |
793ff80c | 123 | |
124 | // | |
125 | // Thickness of the the material layers | |
126 | // | |
db30bf0f | 127 | const Float_t AliTRDgeometry::fgkRaThick = 0.3646; |
793ff80c | 128 | const Float_t AliTRDgeometry::fgkMyThick = 0.005; |
0a770ac9 | 129 | const Float_t AliTRDgeometry::fgkDrThick = AliTRDgeometry::fgkCdrH; |
130 | const Float_t AliTRDgeometry::fgkAmThick = AliTRDgeometry::fgkCamH; | |
131 | const Float_t AliTRDgeometry::fgkXeThick = AliTRDgeometry::fgkDrThick | |
132 | + AliTRDgeometry::fgkAmThick; | |
793ff80c | 133 | const Float_t AliTRDgeometry::fgkCuThick = 0.001; |
134 | const Float_t AliTRDgeometry::fgkSuThick = 0.06; | |
135 | const Float_t AliTRDgeometry::fgkFeThick = 0.0044; | |
136 | const Float_t AliTRDgeometry::fgkCoThick = 0.02; | |
db30bf0f | 137 | const Float_t AliTRDgeometry::fgkWaThick = 0.02; |
793ff80c | 138 | |
139 | // | |
140 | // Position of the material layers | |
141 | // | |
0a770ac9 | 142 | const Float_t AliTRDgeometry::fgkRaZpos = -1.50; |
143 | const Float_t AliTRDgeometry::fgkMyZpos = 0.895; | |
144 | const Float_t AliTRDgeometry::fgkDrZpos = 2.4; | |
145 | const Float_t AliTRDgeometry::fgkAmZpos = 0.0; | |
146 | const Float_t AliTRDgeometry::fgkCuZpos = -0.9995; | |
793ff80c | 147 | const Float_t AliTRDgeometry::fgkSuZpos = 0.0000; |
0a770ac9 | 148 | const Float_t AliTRDgeometry::fgkFeZpos = 0.0322; |
149 | const Float_t AliTRDgeometry::fgkCoZpos = 0.97; | |
150 | const Float_t AliTRDgeometry::fgkWaZpos = 0.99; | |
3551db50 | 151 | |
152 | const Double_t AliTRDgeometry::fgkTime0Base = Rmin() + CraHght() + CdrHght() + CamHght()/2.; | |
153 | const Float_t AliTRDgeometry::fgkTime0[6] = { fgkTime0Base + 0 * (Cheight() + Cspace()), | |
154 | fgkTime0Base + 1 * (Cheight() + Cspace()), | |
155 | fgkTime0Base + 2 * (Cheight() + Cspace()), | |
156 | fgkTime0Base + 3 * (Cheight() + Cspace()), | |
157 | fgkTime0Base + 4 * (Cheight() + Cspace()), | |
158 | fgkTime0Base + 5 * (Cheight() + Cspace()) }; | |
793ff80c | 159 | |
f7336fa3 | 160 | //_____________________________________________________________________________ |
161 | AliTRDgeometry::AliTRDgeometry():AliGeometry() | |
162 | { | |
163 | // | |
164 | // AliTRDgeometry default constructor | |
165 | // | |
bd0f8685 | 166 | |
167 | fMatrixArray = 0; | |
168 | fMatrixCorrectionArray = 0; | |
169 | ||
f7336fa3 | 170 | Init(); |
bd0f8685 | 171 | |
f7336fa3 | 172 | } |
173 | ||
174 | //_____________________________________________________________________________ | |
175 | AliTRDgeometry::~AliTRDgeometry() | |
176 | { | |
8230f242 | 177 | // |
178 | // AliTRDgeometry destructor | |
179 | // | |
bd0f8685 | 180 | |
b4a9cd27 | 181 | delete fMatrixArray; |
182 | delete fMatrixCorrectionArray; | |
bd0f8685 | 183 | |
f7336fa3 | 184 | } |
185 | ||
186 | //_____________________________________________________________________________ | |
187 | void AliTRDgeometry::Init() | |
188 | { | |
189 | // | |
190 | // Initializes the geometry parameter | |
191 | // | |
f7336fa3 | 192 | // The maximum number of pads |
193 | // and the position of pad 0,0,0 | |
194 | // | |
195 | // chambers seen from the top: | |
196 | // +----------------------------+ | |
197 | // | | | |
793ff80c | 198 | // | | ^ |
199 | // | | rphi| | |
200 | // | | | | |
201 | // |0 | | | |
202 | // +----------------------------+ +------> | |
f7336fa3 | 203 | // z |
793ff80c | 204 | // chambers seen from the side: ^ |
205 | // +----------------------------+ drift| | |
206 | // |0 | | | |
207 | // | | | | |
208 | // +----------------------------+ +------> | |
f7336fa3 | 209 | // z |
210 | // | |
a2b90f83 | 211 | // IMPORTANT: time bin 0 is now the first one in the drift region |
212 | // closest to the readout !!! | |
793ff80c | 213 | // |
f7336fa3 | 214 | |
0a770ac9 | 215 | Int_t icham; |
216 | Int_t iplan; | |
217 | Int_t isect; | |
218 | ||
219 | // The outer width of the chambers | |
287c5d50 | 220 | fCwidth[0] = 90.4; |
e0d47c25 | 221 | fCwidth[1] = 94.8; |
222 | fCwidth[2] = 99.3; | |
223 | fCwidth[3] = 103.7; | |
224 | fCwidth[4] = 108.1; | |
225 | fCwidth[5] = 112.6; | |
0a770ac9 | 226 | |
227 | // The outer lengths of the chambers | |
73ae7b59 | 228 | // Includes the spacings between the chambers! |
8737e16f | 229 | Float_t length[kNplan][kNcham] = { { 124.0, 124.0, 110.0, 124.0, 124.0 } |
e0d47c25 | 230 | , { 124.0, 124.0, 110.0, 124.0, 124.0 } |
8737e16f | 231 | , { 131.0, 131.0, 110.0, 131.0, 131.0 } |
232 | , { 138.0, 138.0, 110.0, 138.0, 138.0 } | |
233 | , { 145.0, 145.0, 110.0, 145.0, 145.0 } | |
e0d47c25 | 234 | , { 147.0, 147.0, 110.0, 147.0, 147.0 } }; |
0a770ac9 | 235 | |
236 | for (icham = 0; icham < kNcham; icham++) { | |
237 | for (iplan = 0; iplan < kNplan; iplan++) { | |
238 | fClength[iplan][icham] = length[iplan][icham]; | |
0a770ac9 | 239 | } |
240 | } | |
241 | ||
793ff80c | 242 | // The rotation matrix elements |
243 | Float_t phi = 0; | |
244 | for (isect = 0; isect < fgkNsect; isect++) { | |
5443e65e | 245 | phi = -2.0 * TMath::Pi() / (Float_t) fgkNsect * ((Float_t) isect + 0.5); |
793ff80c | 246 | fRotA11[isect] = TMath::Cos(phi); |
247 | fRotA12[isect] = TMath::Sin(phi); | |
248 | fRotA21[isect] = TMath::Sin(phi); | |
249 | fRotA22[isect] = TMath::Cos(phi); | |
250 | phi = -1.0 * phi; | |
251 | fRotB11[isect] = TMath::Cos(phi); | |
252 | fRotB12[isect] = TMath::Sin(phi); | |
253 | fRotB21[isect] = TMath::Sin(phi); | |
254 | fRotB22[isect] = TMath::Cos(phi); | |
255 | } | |
bd0f8685 | 256 | |
257 | for (isect = 0; isect < fgkNsect; isect++) { | |
258 | SetSMstatus(isect,1); | |
259 | } | |
793ff80c | 260 | |
261 | } | |
262 | ||
f7336fa3 | 263 | //_____________________________________________________________________________ |
bd0f8685 | 264 | void AliTRDgeometry::CreateGeometry(Int_t *idtmed) |
265 | { | |
266 | // | |
267 | // Create the TRD geometry without hole | |
268 | // | |
269 | // | |
270 | // Names of the TRD volumina (xx = detector number): | |
271 | // | |
272 | // Volume (Air) wrapping the readout chamber components | |
273 | // UTxx includes: UAxx, UDxx, UFxx, UUxx | |
274 | // Obs: | |
275 | // UUxx the services volume has been reduced by 7.42 mm | |
276 | // in order to allow shifts in radial direction | |
277 | // | |
278 | // Lower part of the readout chambers (gas volume + radiator) | |
279 | // | |
280 | // UAxx Aluminum frames (Al) | |
281 | // UBxx G10 frames (C) | |
282 | // UCxx Inner volumes (Air) | |
283 | // | |
284 | // Upper part of the readout chambers (readout plane + fee) | |
285 | // | |
286 | // UDxx G10 frames (C) | |
287 | // UExx Inner volumes of the G10 (Air) | |
288 | // UFxx Aluminum frames (Al) | |
289 | // UGxx Inner volumes of the Al (Air) | |
290 | // | |
291 | // Inner material layers | |
292 | // | |
293 | // UHxx Radiator (Rohacell) | |
294 | // UIxx Entrance window (Mylar) | |
295 | // UJxx Drift volume (Xe/CO2) | |
296 | // UKxx Amplification volume (Xe/CO2) | |
297 | // ULxx Pad plane (Cu) | |
298 | // UMxx Support structure (Rohacell) | |
299 | // | |
300 | ||
301 | const Int_t kNparTrd = 4; | |
302 | const Int_t kNparCha = 3; | |
303 | ||
304 | Float_t xpos, ypos, zpos; | |
305 | ||
306 | Float_t parTrd[kNparTrd]; | |
307 | Float_t parCha[kNparCha]; | |
308 | ||
309 | Char_t cTagV[6]; | |
310 | Char_t cTagM[5]; | |
311 | ||
312 | // The TRD mother volume for one sector (Air), full length in z-direction | |
313 | // Provides material for side plates of super module | |
314 | parTrd[0] = fgkSwidth1/2.; | |
315 | parTrd[1] = fgkSwidth2/2.; | |
316 | parTrd[2] = fgkSlenTR1/2.; | |
317 | parTrd[3] = fgkSheight/2.; | |
318 | gMC->Gsvolu("UTR1","TRD1",idtmed[1302-1],parTrd,kNparTrd); | |
319 | ||
320 | // | |
321 | // The side plates of the super module (Al) | |
322 | parTrd[0] = fgkSwidth1/2. - fgkSMgapT; | |
323 | parTrd[1] = fgkSwidth2/2. - fgkSMgapT; | |
324 | parTrd[2] = fgkSlenTR1/2.; | |
325 | parTrd[3] = fgkSheight/2.; | |
326 | gMC->Gsvolu("UTS1","TRD1",idtmed[1301-1],parTrd,kNparTrd); | |
327 | ||
328 | // The inner part of the TRD mother volume for one sector (Air), | |
329 | // full length in z-direction | |
330 | parTrd[0] = fgkSwidth1/2. - fgkSMgapT - fgkSMpltT; | |
331 | parTrd[1] = fgkSwidth2/2. - fgkSMgapT - fgkSMpltT; | |
332 | parTrd[2] = fgkSlenTR1/2.; | |
333 | parTrd[3] = fgkSheight/2.; | |
334 | gMC->Gsvolu("UTI1","TRD1",idtmed[1302-1],parTrd,kNparTrd); | |
335 | ||
336 | for (Int_t icham = 0; icham < kNcham; icham++) { | |
337 | for (Int_t iplan = 0; iplan < kNplan; iplan++) { | |
338 | ||
339 | Int_t iDet = GetDetectorSec(iplan,icham); | |
340 | ||
341 | // The lower part of the readout chambers (gas volume + radiator) | |
342 | // The aluminum frames | |
343 | sprintf(cTagV,"UA%02d",iDet); | |
344 | parCha[0] = fCwidth[iplan]/2.; | |
345 | parCha[1] = fClength[iplan][icham]/2. - fgkHspace/2.; | |
346 | parCha[2] = fgkCraH/2. + fgkCdrH/2.; | |
347 | fChamberUAboxd[iDet][0] = parCha[0]; | |
348 | fChamberUAboxd[iDet][1] = parCha[1]; | |
349 | fChamberUAboxd[iDet][2] = parCha[2]; | |
350 | gMC->Gsvolu(cTagV,"BOX ",idtmed[1301-1],parCha,kNparCha); | |
351 | // The G10 frames | |
352 | sprintf(cTagV,"UB%02d",iDet); | |
353 | parCha[0] = fCwidth[iplan]/2. - fgkCalT; | |
354 | parCha[1] = -1.; | |
355 | parCha[2] = -1.; | |
356 | gMC->Gsvolu(cTagV,"BOX ",idtmed[1307-1],parCha,kNparCha); | |
357 | // The inner part (air) | |
358 | sprintf(cTagV,"UC%02d",iDet); | |
359 | parCha[0] = fCwidth[iplan]/2. - fgkCalT - fgkCclsT; | |
360 | parCha[1] = fClength[iplan][icham]/2. - fgkHspace/2.- fgkCclfT; | |
361 | parCha[2] = -1.; | |
362 | gMC->Gsvolu(cTagV,"BOX ",idtmed[1302-1],parCha,kNparCha); | |
363 | ||
364 | // The upper part of the readout chambers (readout plane) | |
365 | // The G10 frames | |
366 | sprintf(cTagV,"UD%02d",iDet); | |
367 | parCha[0] = fCwidth[iplan]/2. + fgkCroW; | |
368 | parCha[1] = fClength[iplan][icham]/2. - fgkHspace/2.; | |
369 | parCha[2] = fgkCamH/2.; | |
370 | fChamberUDboxd[iDet][0] = parCha[0]; | |
371 | fChamberUDboxd[iDet][1] = parCha[1]; | |
372 | fChamberUDboxd[iDet][2] = parCha[2]; | |
373 | gMC->Gsvolu(cTagV,"BOX ",idtmed[1307-1],parCha,kNparCha); | |
374 | // The inner part of the G10 frame (air) | |
375 | sprintf(cTagV,"UE%02d",iDet); | |
376 | parCha[0] = fCwidth[iplan]/2. + fgkCroW - fgkCcuT; | |
377 | parCha[1] = fClength[iplan][icham]/2. - fgkHspace/2.- fgkCcuT; | |
378 | parCha[2] = -1.; | |
379 | gMC->Gsvolu(cTagV,"BOX ",idtmed[1302-1],parCha,kNparCha); | |
380 | // The aluminum frames | |
381 | sprintf(cTagV,"UF%02d",iDet); | |
382 | parCha[0] = fCwidth[iplan]/2. + fgkCroW; | |
383 | parCha[1] = fClength[iplan][icham]/2. - fgkHspace/2.; | |
384 | parCha[2] = fgkCroH/2.; | |
385 | fChamberUFboxd[iDet][0] = parCha[0]; | |
386 | fChamberUFboxd[iDet][1] = parCha[1]; | |
387 | fChamberUFboxd[iDet][2] = parCha[2]; | |
388 | gMC->Gsvolu(cTagV,"BOX ",idtmed[1301-1],parCha,kNparCha); | |
389 | // The inner part of the aluminum frames | |
390 | sprintf(cTagV,"UG%02d",iDet); | |
391 | parCha[0] = fCwidth[iplan]/2. + fgkCroW - fgkCauT; | |
392 | parCha[1] = fClength[iplan][icham]/2. - fgkHspace/2.- fgkCauT; | |
393 | parCha[2] = -1.; | |
394 | gMC->Gsvolu(cTagV,"BOX ",idtmed[1302-1],parCha,kNparCha); | |
395 | ||
396 | // The material layers inside the chambers | |
397 | parCha[0] = -1.; | |
398 | parCha[1] = -1.; | |
399 | // Rohacell layer (radiator) | |
400 | parCha[2] = fgkRaThick/2; | |
401 | sprintf(cTagV,"UH%02d",iDet); | |
402 | gMC->Gsvolu(cTagV,"BOX ",idtmed[1315-1],parCha,kNparCha); | |
403 | // Mylar layer (entrance window + HV cathode) | |
404 | parCha[2] = fgkMyThick/2; | |
405 | sprintf(cTagV,"UI%02d",iDet); | |
406 | gMC->Gsvolu(cTagV,"BOX ",idtmed[1308-1],parCha,kNparCha); | |
407 | // Xe/Isobutane layer (drift volume) | |
408 | parCha[2] = fgkDrThick/2.; | |
409 | sprintf(cTagV,"UJ%02d",iDet); | |
410 | gMC->Gsvolu(cTagV,"BOX ",idtmed[1309-1],parCha,kNparCha); | |
411 | // Xe/Isobutane layer (amplification volume) | |
412 | parCha[2] = fgkAmThick/2.; | |
413 | sprintf(cTagV,"UK%02d",iDet); | |
414 | gMC->Gsvolu(cTagV,"BOX ",idtmed[1309-1],parCha,kNparCha); | |
415 | // Cu layer (pad plane) | |
416 | parCha[2] = fgkCuThick/2; | |
417 | sprintf(cTagV,"UL%02d",iDet); | |
418 | gMC->Gsvolu(cTagV,"BOX ",idtmed[1305-1],parCha,kNparCha); | |
419 | // G10 layer (support structure / honeycomb) | |
420 | parCha[2] = fgkSuThick/2; | |
421 | sprintf(cTagV,"UM%02d",iDet); | |
422 | gMC->Gsvolu(cTagV,"BOX ",idtmed[1313-1],parCha,kNparCha); | |
423 | ||
424 | // Position the layers in the chambers | |
425 | xpos = 0; | |
426 | ypos = 0; | |
427 | // Lower part | |
428 | // Rohacell layer (radiator) | |
429 | zpos = fgkRaZpos; | |
430 | sprintf(cTagV,"UH%02d",iDet); | |
431 | sprintf(cTagM,"UC%02d",iDet); | |
432 | gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY"); | |
433 | // Mylar layer (entrance window + HV cathode) | |
434 | zpos = fgkMyZpos; | |
435 | sprintf(cTagV,"UI%02d",iDet); | |
436 | sprintf(cTagM,"UC%02d",iDet); | |
437 | gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY"); | |
438 | // Xe/Isobutane layer (drift volume) | |
439 | zpos = fgkDrZpos; | |
440 | sprintf(cTagV,"UJ%02d",iDet); | |
441 | sprintf(cTagM,"UC%02d",iDet); | |
442 | gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY"); | |
443 | // Upper part | |
444 | // Xe/Isobutane layer (amplification volume) | |
445 | zpos = fgkAmZpos; | |
446 | sprintf(cTagV,"UK%02d",iDet); | |
447 | sprintf(cTagM,"UE%02d",iDet); | |
448 | gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY"); | |
449 | // Readout part | |
450 | // Cu layer (pad plane) | |
451 | zpos = fgkCuZpos; | |
452 | sprintf(cTagV,"UL%02d",iDet); | |
453 | sprintf(cTagM,"UG%02d",iDet); | |
454 | gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY"); | |
455 | // G10 layer (support structure) | |
456 | zpos = fgkSuZpos; | |
457 | sprintf(cTagV,"UM%02d",iDet); | |
458 | sprintf(cTagM,"UG%02d",iDet); | |
459 | gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY"); | |
460 | ||
461 | // Position the inner volumes of the chambers in the frames | |
462 | xpos = 0.0; | |
463 | ypos = 0.0; | |
464 | zpos = 0.0; | |
465 | // The inside of the lower G10 frame | |
466 | sprintf(cTagV,"UC%02d",iDet); | |
467 | sprintf(cTagM,"UB%02d",iDet); | |
468 | gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY"); | |
469 | // The lower G10 frame inside the aluminum frame | |
470 | sprintf(cTagV,"UB%02d",iDet); | |
471 | sprintf(cTagM,"UA%02d",iDet); | |
472 | gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY"); | |
473 | // The inside of the upper G10 frame | |
474 | sprintf(cTagV,"UE%02d",iDet); | |
475 | sprintf(cTagM,"UD%02d",iDet); | |
476 | gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY"); | |
477 | // The inside of the upper aluminum frame | |
478 | sprintf(cTagV,"UG%02d",iDet); | |
479 | sprintf(cTagM,"UF%02d",iDet); | |
480 | gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY"); | |
481 | ||
482 | // Position the frames of the chambers in the TRD mother volume | |
483 | xpos = 0.; | |
484 | ypos = - fClength[iplan][0] - fClength[iplan][1] - fClength[iplan][2]/2.; | |
485 | for (Int_t ic = 0; ic < icham; ic++) { | |
486 | ypos += fClength[iplan][ic]; | |
487 | } | |
488 | ypos += fClength[iplan][icham]/2.; | |
489 | zpos = fgkCraH/2. + fgkCdrH/2. - fgkSheight/2. + iplan * (fgkCH + fgkVspace); | |
490 | // The lower aluminum frame, radiator + drift region | |
491 | sprintf(cTagV,"UA%02d",iDet); | |
492 | fChamberUAorig[iDet][0] = xpos; | |
493 | fChamberUAorig[iDet][1] = ypos; | |
494 | fChamberUAorig[iDet][2] = zpos; | |
495 | // The upper G10 frame, amplification region | |
496 | sprintf(cTagV,"UD%02d",iDet); | |
497 | zpos += fgkCamH/2. + fgkCraH/2. + fgkCdrH/2.; | |
498 | fChamberUDorig[iDet][0] = xpos; | |
499 | fChamberUDorig[iDet][1] = ypos; | |
500 | fChamberUDorig[iDet][2] = zpos; | |
501 | // The upper aluminum frame | |
502 | sprintf(cTagV,"UF%02d",iDet); | |
503 | zpos += fgkCroH/2. + fgkCamH/2.; | |
504 | fChamberUForig[iDet][0] = xpos; | |
505 | fChamberUForig[iDet][1] = ypos; | |
506 | fChamberUForig[iDet][2] = zpos; | |
507 | ||
508 | } | |
509 | } | |
510 | ||
511 | // Create the volumes of the super module frame | |
512 | CreateFrame(idtmed); | |
513 | ||
514 | // Create the volumes of the services | |
515 | CreateServices(idtmed); | |
516 | ||
517 | for (Int_t icham = 0; icham < kNcham; icham++) { | |
518 | for (Int_t iplan = 0; iplan < kNplan; iplan++) { | |
519 | GroupChamber(iplan,icham,idtmed); | |
520 | } | |
521 | } | |
522 | ||
523 | xpos = 0.; | |
524 | ypos = 0.; | |
525 | zpos = 0.; | |
526 | gMC->Gspos("UTI1",1,"UTS1",xpos,ypos,zpos,0,"ONLY"); | |
527 | ||
528 | xpos = 0.; | |
529 | ypos = 0.; | |
530 | zpos = 0.; | |
531 | gMC->Gspos("UTS1",1,"UTR1",xpos,ypos,zpos,0,"ONLY"); | |
532 | ||
533 | // Put the TRD volumes into the space frame mother volumes | |
534 | // if enabled via status flag | |
535 | xpos = 0.; | |
536 | ypos = 0.; | |
537 | zpos = 0.; | |
538 | for (Int_t isect = 0; isect < kNsect; isect++) { | |
539 | if (fSMstatus[isect]) { | |
540 | sprintf(cTagV,"BTRD%d",isect); | |
541 | gMC->Gspos("UTR1",1,cTagV,xpos,ypos,zpos,0,"ONLY"); | |
542 | } | |
543 | } | |
544 | ||
545 | } | |
546 | ||
547 | //_____________________________________________________________________________ | |
548 | void AliTRDgeometry::CreateFrame(Int_t *idtmed) | |
549 | { | |
550 | // | |
551 | // Create the geometry of the frame of the supermodule | |
552 | // | |
553 | // Names of the TRD services volumina | |
554 | // | |
555 | // USRL Support rails for the chambers (Al) | |
556 | // USxx Support cross bars between the chambers (Al) | |
557 | // | |
558 | ||
559 | Int_t iplan = 0; | |
560 | ||
561 | Float_t xpos = 0.0; | |
562 | Float_t ypos = 0.0; | |
563 | Float_t zpos = 0.0; | |
564 | ||
565 | Char_t cTagV[5]; | |
566 | ||
567 | // | |
568 | // The chamber support rails | |
569 | // | |
570 | ||
571 | const Float_t kSRLwid = 2.0; | |
572 | const Float_t kSRLhgt = 2.3; | |
573 | const Float_t kSRLdst = 0.6; | |
574 | const Int_t kNparSRL = 3; | |
575 | Float_t parSRL[kNparSRL]; | |
576 | parSRL[0] = kSRLwid/2.; | |
577 | parSRL[1] = fgkSlenTR1/2.; | |
578 | parSRL[2] = kSRLhgt/2.; | |
579 | gMC->Gsvolu("USRL","BOX ",idtmed[1301-1],parSRL,kNparSRL); | |
580 | ||
581 | xpos = 0.0; | |
582 | ypos = 0.0; | |
583 | zpos = 0.0; | |
584 | for (iplan = 0; iplan < kNplan; iplan++) { | |
585 | ||
586 | xpos = fCwidth[iplan]/2. + kSRLwid/2. + kSRLdst; | |
587 | ypos = 0.0; | |
588 | zpos = fgkCraH + fgkCdrH - fgkSheight/2. - kSRLhgt/2. | |
589 | + iplan * (fgkCH + fgkVspace); | |
590 | gMC->Gspos("USRL",iplan+1 ,"UTI1", xpos,ypos,zpos,0,"ONLY"); | |
591 | gMC->Gspos("USRL",iplan+1+ kNplan,"UTI1",-xpos,ypos,zpos,0,"ONLY"); | |
592 | ||
593 | } | |
594 | ||
595 | // | |
596 | // The cross bars between the chambers | |
597 | // | |
598 | ||
599 | const Float_t kSCBwid = 1.0; | |
600 | const Int_t kNparSCB = 3; | |
601 | Float_t parSCB[kNparSCB]; | |
602 | parSCB[1] = kSCBwid/2.; | |
603 | parSCB[2] = fgkCH/2.; | |
604 | ||
605 | xpos = 0.0; | |
606 | ypos = 0.0; | |
607 | zpos = 0.0; | |
608 | for (iplan = 0; iplan < kNplan; iplan++) { | |
609 | ||
610 | parSCB[0] = fCwidth[iplan]/2. + kSRLdst/2.; | |
611 | ||
612 | sprintf(cTagV,"US0%01d",iplan); | |
613 | gMC->Gsvolu(cTagV,"BOX ",idtmed[1301-1],parSCB,kNparSCB); | |
614 | xpos = 0.0; | |
615 | ypos = fgkSlenTR1/2. - kSCBwid/2.; | |
616 | zpos = fgkCH/2. - fgkSheight/2. + iplan * (fgkCH + fgkVspace); | |
617 | gMC->Gspos(cTagV,1,"UTI1", xpos,ypos,zpos,0,"ONLY"); | |
618 | ||
619 | sprintf(cTagV,"US1%01d",iplan); | |
620 | gMC->Gsvolu(cTagV,"BOX ",idtmed[1301-1],parSCB,kNparSCB); | |
621 | xpos = 0.0; | |
622 | ypos = fClength[iplan][2]/2. + fClength[iplan][1]; | |
623 | zpos = fgkCH/2. - fgkSheight/2. + iplan * (fgkCH + fgkVspace); | |
624 | gMC->Gspos(cTagV,1,"UTI1", xpos,ypos,zpos,0,"ONLY"); | |
625 | ||
626 | sprintf(cTagV,"US2%01d",iplan); | |
627 | gMC->Gsvolu(cTagV,"BOX ",idtmed[1301-1],parSCB,kNparSCB); | |
628 | xpos = 0.0; | |
629 | ypos = fClength[iplan][2]/2.; | |
630 | zpos = fgkCH/2. - fgkSheight/2. + iplan * (fgkCH + fgkVspace); | |
631 | gMC->Gspos(cTagV,1,"UTI1", xpos,ypos,zpos,0,"ONLY"); | |
632 | ||
633 | sprintf(cTagV,"US3%01d",iplan); | |
634 | gMC->Gsvolu(cTagV,"BOX ",idtmed[1301-1],parSCB,kNparSCB); | |
635 | xpos = 0.0; | |
636 | ypos = - fClength[iplan][2]/2.; | |
637 | zpos = fgkCH/2. - fgkSheight/2. + iplan * (fgkCH + fgkVspace); | |
638 | gMC->Gspos(cTagV,1,"UTI1", xpos,ypos,zpos,0,"ONLY"); | |
639 | ||
640 | sprintf(cTagV,"US4%01d",iplan); | |
641 | gMC->Gsvolu(cTagV,"BOX ",idtmed[1301-1],parSCB,kNparSCB); | |
642 | xpos = 0.0; | |
643 | ypos = - fClength[iplan][2]/2. - fClength[iplan][1]; | |
644 | zpos = fgkCH/2. - fgkSheight/2. + iplan * (fgkCH + fgkVspace); | |
645 | gMC->Gspos(cTagV,1,"UTI1", xpos,ypos,zpos,0,"ONLY"); | |
646 | ||
647 | sprintf(cTagV,"US5%01d",iplan); | |
648 | gMC->Gsvolu(cTagV,"BOX ",idtmed[1301-1],parSCB,kNparSCB); | |
649 | xpos = 0.0; | |
650 | ypos = - fgkSlenTR1/2. + kSCBwid/2.; | |
651 | zpos = fgkCH/2. - fgkSheight/2. + iplan * (fgkCH + fgkVspace); | |
652 | gMC->Gspos(cTagV,1,"UTI1", xpos,ypos,zpos,0,"ONLY"); | |
653 | ||
654 | } | |
655 | ||
656 | } | |
657 | ||
658 | //_____________________________________________________________________________ | |
659 | void AliTRDgeometry::CreateServices(Int_t *idtmed) | |
660 | { | |
661 | // | |
662 | // Create the geometry of the services | |
663 | // | |
664 | // Names of the TRD services volumina | |
665 | // | |
666 | // UTCL Cooling arterias (Al) | |
667 | // UTCW Cooling arterias (Water) | |
668 | // UUxx Volumes for the services at the chambers (Air) | |
669 | // UTPW Power bars (Cu) | |
670 | // UTCP Cooling pipes (Al) | |
671 | // UTCH Cooling pipes (Water) | |
672 | // UTPL Power lines (Cu) | |
673 | // UMCM Readout MCMs (G10/Cu/Si) | |
674 | // | |
675 | ||
676 | Int_t iplan = 0; | |
677 | Int_t icham = 0; | |
678 | ||
679 | Float_t xpos = 0.0; | |
680 | Float_t ypos = 0.0; | |
681 | Float_t zpos = 0.0; | |
682 | ||
683 | Char_t cTagV[5]; | |
684 | ||
685 | // The rotation matrices | |
686 | const Int_t kNmatrix = 3; | |
687 | Int_t matrix[kNmatrix]; | |
688 | gMC->Matrix(matrix[0],100.0, 0.0, 90.0, 90.0, 10.0, 0.0); | |
689 | gMC->Matrix(matrix[1], 80.0, 0.0, 90.0, 90.0, 10.0,180.0); | |
690 | gMC->Matrix(matrix[2], 0.0, 0.0, 90.0, 90.0, 90.0, 0.0); | |
691 | ||
692 | AliTRDCommonParam* commonParam = AliTRDCommonParam::Instance(); | |
693 | if (!commonParam) | |
694 | { | |
695 | AliError("Could not get common params\n"); | |
696 | return; | |
697 | } | |
698 | ||
699 | // | |
700 | // The cooling arterias | |
701 | // | |
702 | ||
703 | // Width of the cooling arterias | |
704 | const Float_t kCOLwid = 0.5; | |
705 | // Height of the cooling arterias | |
706 | const Float_t kCOLhgt = 5.5; | |
707 | // Positioning of the cooling | |
708 | const Float_t kCOLposx = 1.6; | |
709 | const Float_t kCOLposz = -0.2; | |
710 | // Thickness of the walls of the cooling arterias | |
711 | const Float_t kCOLthk = 0.1; | |
712 | const Int_t kNparCOL = 3; | |
713 | Float_t parCOL[kNparCOL]; | |
714 | parCOL[0] = kCOLwid/2.; | |
715 | parCOL[1] = fgkSlenTR1/2.; | |
716 | parCOL[2] = kCOLhgt/2.; | |
717 | gMC->Gsvolu("UTCL","BOX ",idtmed[1324-1],parCOL,kNparCOL); | |
718 | parCOL[0] -= kCOLthk; | |
719 | parCOL[1] = fgkSlenTR1/2.; | |
720 | parCOL[2] -= kCOLthk; | |
721 | gMC->Gsvolu("UTCW","BOX ",idtmed[1314-1],parCOL,kNparCOL); | |
722 | ||
723 | xpos = 0.0; | |
724 | ypos = 0.0; | |
725 | zpos = 0.0; | |
726 | gMC->Gspos("UTCW",1,"UTCL", xpos,ypos,zpos,0,"ONLY"); | |
727 | ||
728 | for (iplan = 0; iplan < kNplan; iplan++) { // CHECK FOR OVERLAPS !!! | |
729 | //for (iplan = 1; iplan < kNplan; iplan++) { | |
730 | ||
731 | xpos = fCwidth[iplan]/2. + kCOLwid/2. + kCOLposx; | |
732 | ypos = 0.0; | |
733 | zpos = kCOLhgt/2. - fgkSheight/2. + kCOLposz + iplan * (fgkCH + fgkVspace); | |
734 | if (iplan == 0) zpos += 0.25; // To avoid overlaps ! | |
735 | gMC->Gspos("UTCL",iplan+1 ,"UTI1", xpos,ypos,zpos,matrix[0],"ONLY"); | |
736 | gMC->Gspos("UTCL",iplan+1+ kNplan,"UTI1",-xpos,ypos,zpos,matrix[1],"ONLY"); | |
737 | ||
738 | } | |
739 | ||
740 | // | |
741 | // The power bars | |
742 | // | |
743 | ||
744 | const Float_t kPWRwid = 0.6; | |
745 | const Float_t kPWRhgt = 4.5; | |
746 | const Float_t kPWRposx = 1.05; | |
747 | const Float_t kPWRposz = 0.9; | |
748 | const Int_t kNparPWR = 3; | |
749 | Float_t parPWR[kNparPWR]; | |
750 | parPWR[0] = kPWRwid/2.; | |
751 | parPWR[1] = fgkSlenTR1/2.; | |
752 | parPWR[2] = kPWRhgt/2.; | |
753 | gMC->Gsvolu("UTPW","BOX ",idtmed[1325-1],parPWR,kNparPWR); | |
754 | ||
755 | for (iplan = 0; iplan < kNplan; iplan++) { // CHECK FOR OVERLAPS !!! | |
756 | //for (iplan = 1; iplan < kNplan; iplan++) { | |
757 | ||
758 | xpos = fCwidth[iplan]/2. + kPWRwid/2. + kPWRposx; | |
759 | ypos = 0.0; | |
760 | zpos = kPWRhgt/2. - fgkSheight/2. + kPWRposz + iplan * (fgkCH + fgkVspace); | |
761 | gMC->Gspos("UTPW",iplan+1 ,"UTI1", xpos,ypos,zpos,matrix[0],"ONLY"); | |
762 | gMC->Gspos("UTPW",iplan+1+ kNplan,"UTI1",-xpos,ypos,zpos,matrix[1],"ONLY"); | |
763 | ||
764 | } | |
765 | ||
766 | // | |
767 | // The volumes for the services at the chambers | |
768 | // | |
769 | ||
770 | const Int_t kNparServ = 3; | |
771 | Float_t parServ[kNparServ]; | |
772 | ||
773 | for (icham = 0; icham < kNcham; icham++) { | |
774 | for (iplan = 0; iplan < kNplan; iplan++) { | |
775 | // Take out upper plane until TRD mothervolume is adjusted | |
776 | //for (iplan = 0; iplan < kNplan-1; iplan++) { | |
777 | ||
778 | Int_t iDet = GetDetectorSec(iplan,icham); | |
779 | ||
780 | sprintf(cTagV,"UU%02d",iDet); | |
781 | parServ[0] = fCwidth[iplan]/2.; | |
782 | parServ[1] = fClength[iplan][icham]/2. - fgkHspace/2.; | |
783 | parServ[2] = fgkVspace/2. - 0.742/2.; | |
784 | fChamberUUboxd[iDet][0] = parServ[0]; | |
785 | fChamberUUboxd[iDet][1] = parServ[1]; | |
786 | fChamberUUboxd[iDet][2] = parServ[2]; | |
787 | ||
788 | gMC->Gsvolu(cTagV,"BOX",idtmed[1302-1],parServ,kNparServ); | |
789 | xpos = 0.; | |
790 | ypos = - fClength[iplan][0] - fClength[iplan][1] - fClength[iplan][2]/2.; | |
791 | for (Int_t ic = 0; ic < icham; ic++) { | |
792 | ypos += fClength[iplan][ic]; | |
793 | } | |
794 | ypos += fClength[iplan][icham]/2.; | |
795 | zpos = fgkCH + fgkVspace/2. - fgkSheight/2. + iplan * (fgkCH + fgkVspace); | |
796 | zpos -= 0.742/2.; | |
797 | fChamberUUorig[iDet][0] = xpos; | |
798 | fChamberUUorig[iDet][1] = ypos; | |
799 | fChamberUUorig[iDet][2] = zpos; | |
800 | ||
801 | } | |
802 | } | |
803 | ||
804 | // | |
805 | // The cooling pipes inside the service volumes | |
806 | // | |
807 | ||
808 | const Int_t kNparTube = 3; | |
809 | Float_t parTube[kNparTube]; | |
810 | // The aluminum pipe for the cooling | |
811 | parTube[0] = 0.0; | |
812 | parTube[1] = 0.0; | |
813 | parTube[2] = 0.0; | |
814 | gMC->Gsvolu("UTCP","TUBE",idtmed[1324-1],parTube,0); | |
815 | // The cooling water | |
816 | parTube[0] = 0.0; | |
817 | parTube[1] = 0.2/2.; | |
818 | parTube[2] = -1.; | |
819 | gMC->Gsvolu("UTCH","TUBE",idtmed[1314-1],parTube,kNparTube); | |
820 | // Water inside the cooling pipe | |
821 | xpos = 0.0; | |
822 | ypos = 0.0; | |
823 | zpos = 0.0; | |
824 | gMC->Gspos("UTCH",1,"UTCP",xpos,ypos,zpos,0,"ONLY"); | |
825 | ||
826 | // Position the cooling pipes in the mother volume | |
827 | const Int_t kNpar = 3; | |
828 | Float_t par[kNpar]; | |
829 | for (icham = 0; icham < kNcham; icham++) { | |
830 | for (iplan = 0; iplan < kNplan; iplan++) { | |
831 | // Take out upper plane until TRD mothervolume is adjusted | |
832 | //for (iplan = 0; iplan < kNplan-1; iplan++) { | |
833 | Int_t iDet = GetDetectorSec(iplan,icham); | |
834 | Int_t iCopy = GetDetector(iplan,icham,0) * 100; | |
835 | Int_t nMCMrow = commonParam->GetRowMax(iplan,icham,0); | |
836 | Float_t ySize = (GetChamberLength(iplan,icham) - 2.*fgkRpadW) | |
837 | / ((Float_t) nMCMrow); | |
838 | sprintf(cTagV,"UU%02d",iDet); | |
839 | for (Int_t iMCMrow = 0; iMCMrow < nMCMrow; iMCMrow++) { | |
840 | xpos = 0.0; | |
841 | ypos = (0.5 + iMCMrow) * ySize - 1.9 | |
842 | - fClength[iplan][icham]/2. + fgkHspace/2.; | |
843 | zpos = 0.0 + 0.742/2.; | |
844 | par[0] = 0.0; | |
845 | par[1] = 0.3/2.; // Thickness of the cooling pipes | |
846 | par[2] = fCwidth[iplan]/2.; | |
847 | gMC->Gsposp("UTCP",iCopy+iMCMrow,cTagV,xpos,ypos,zpos | |
848 | ,matrix[2],"ONLY",par,kNpar); | |
849 | } | |
850 | } | |
851 | } | |
852 | ||
853 | // | |
854 | // The power lines | |
855 | // | |
856 | ||
857 | // The copper power lines | |
858 | parTube[0] = 0.0; | |
859 | parTube[1] = 0.0; | |
860 | parTube[2] = 0.0; | |
861 | gMC->Gsvolu("UTPL","TUBE",idtmed[1305-1],parTube,0); | |
862 | ||
863 | // Position the power lines in the mother volume | |
864 | for (icham = 0; icham < kNcham; icham++) { | |
865 | for (iplan = 0; iplan < kNplan; iplan++) { | |
866 | // Take out upper plane until TRD mothervolume is adjusted | |
867 | //for (iplan = 0; iplan < kNplan-1; iplan++) { | |
868 | Int_t iDet = GetDetectorSec(iplan,icham); | |
869 | Int_t iCopy = GetDetector(iplan,icham,0) * 100; | |
870 | Int_t nMCMrow = commonParam->GetRowMax(iplan,icham,0); | |
871 | Float_t ySize = (GetChamberLength(iplan,icham) - 2.*fgkRpadW) | |
872 | / ((Float_t) nMCMrow); | |
873 | sprintf(cTagV,"UU%02d",iDet); | |
874 | for (Int_t iMCMrow = 0; iMCMrow < nMCMrow; iMCMrow++) { | |
875 | xpos = 0.0; | |
876 | ypos = (0.5 + iMCMrow) * ySize - 1.0 | |
877 | - fClength[iplan][icham]/2. + fgkHspace/2.; | |
878 | zpos = -0.4 + 0.742/2.; | |
879 | par[0] = 0.0; | |
880 | par[1] = 0.2/2.; // Thickness of the power lines | |
881 | par[2] = fCwidth[iplan]/2.; | |
882 | gMC->Gsposp("UTPL",iCopy+iMCMrow,cTagV,xpos,ypos,zpos | |
883 | ,matrix[2],"ONLY",par,kNpar); | |
884 | } | |
885 | } | |
886 | } | |
887 | ||
888 | // | |
889 | // The MCMs | |
890 | // | |
891 | ||
892 | // The mother volume for the MCMs (air) | |
893 | const Int_t kNparMCM = 3; | |
894 | Float_t parMCM[kNparMCM]; | |
895 | parMCM[0] = 3.0/2.; | |
896 | parMCM[1] = 3.0/2.; | |
897 | parMCM[2] = 0.14/2.; | |
898 | gMC->Gsvolu("UMCM","BOX",idtmed[1302-1],parMCM,kNparMCM); | |
899 | ||
900 | // The MCM carrier G10 layer | |
901 | parMCM[0] = 3.0/2.; | |
902 | parMCM[1] = 3.0/2.; | |
903 | parMCM[2] = 0.1/2.; | |
904 | gMC->Gsvolu("UMC1","BOX",idtmed[1319-1],parMCM,kNparMCM); | |
905 | // The MCM carrier Cu layer | |
906 | parMCM[0] = 3.0/2.; | |
907 | parMCM[1] = 3.0/2.; | |
908 | parMCM[2] = 0.0162/2.; | |
909 | gMC->Gsvolu("UMC2","BOX",idtmed[1318-1],parMCM,kNparMCM); | |
910 | // The silicon of the chips | |
911 | parMCM[0] = 3.0/2.; | |
912 | parMCM[1] = 3.0/2.; | |
913 | parMCM[2] = 0.003/2.; | |
914 | gMC->Gsvolu("UMC3","BOX",idtmed[1320-1],parMCM,kNparMCM); | |
915 | ||
916 | // Put the MCM material inside the MCM mother volume | |
917 | xpos = 0.0; | |
918 | ypos = 0.0; | |
919 | zpos = -0.07 + 0.1/2.; | |
920 | gMC->Gspos("UMC1",1,"UMCM",xpos,ypos,zpos,0,"ONLY"); | |
921 | zpos += 0.1/2. + 0.0162/2.; | |
922 | gMC->Gspos("UMC2",1,"UMCM",xpos,ypos,zpos,0,"ONLY"); | |
923 | zpos += 0.00162/2 + 0.003/2.; | |
924 | gMC->Gspos("UMC3",1,"UMCM",xpos,ypos,zpos,0,"ONLY"); | |
925 | ||
926 | // Position the MCMs in the mother volume | |
927 | for (icham = 0; icham < kNcham; icham++) { | |
928 | for (iplan = 0; iplan < kNplan; iplan++) { | |
929 | // Take out upper plane until TRD mothervolume is adjusted | |
930 | //for (iplan = 0; iplan < kNplan-1; iplan++) { | |
931 | Int_t iDet = GetDetectorSec(iplan,icham); | |
932 | Int_t iCopy = GetDetector(iplan,icham,0) * 1000; | |
933 | Int_t nMCMrow = commonParam->GetRowMax(iplan,icham,0); | |
934 | Float_t ySize = (GetChamberLength(iplan,icham) - 2.*fgkRpadW) | |
935 | / ((Float_t) nMCMrow); | |
936 | Int_t nMCMcol = 8; | |
937 | Float_t xSize = (GetChamberWidth(iplan) - 2.* fgkCpadW) | |
938 | / ((Float_t) nMCMcol); | |
939 | sprintf(cTagV,"UU%02d",iDet); | |
940 | for (Int_t iMCMrow = 0; iMCMrow < nMCMrow; iMCMrow++) { | |
941 | for (Int_t iMCMcol = 0; iMCMcol < nMCMcol; iMCMcol++) { | |
942 | xpos = (0.5 + iMCMcol) * xSize + 1.0 | |
943 | - fCwidth[iplan]/2.; | |
944 | ypos = (0.5 + iMCMrow) * ySize + 1.0 | |
945 | - fClength[iplan][icham]/2. + fgkHspace/2.; | |
946 | zpos = -0.4 + 0.742/2.; | |
947 | par[0] = 0.0; | |
948 | par[1] = 0.2/2.; // Thickness of the power lines | |
949 | par[2] = fCwidth[iplan]/2.; | |
950 | gMC->Gspos("UMCM",iCopy+iMCMrow*10+iMCMcol,cTagV | |
951 | ,xpos,ypos,zpos,0,"ONLY"); | |
952 | } | |
953 | } | |
954 | ||
955 | } | |
956 | } | |
957 | ||
958 | } | |
959 | ||
960 | //_____________________________________________________________________________ | |
961 | void AliTRDgeometry::GroupChamber(Int_t iplan, Int_t icham, Int_t *idtmed) | |
f7336fa3 | 962 | { |
963 | // | |
bd0f8685 | 964 | // Group volumes UA, UD, UF, UU in a single chamber (Air) |
965 | // UA, UD, UF, UU are boxes | |
966 | // UT will be a box | |
0a770ac9 | 967 | // |
bd0f8685 | 968 | // ... for the moment there are no services (UU) for the upper plane ! |
969 | // | |
970 | ||
971 | const Int_t kNparCha = 3; | |
972 | ||
973 | Int_t iDet = GetDetectorSec(iplan,icham); | |
974 | ||
975 | Float_t xyzMin[3]; | |
976 | Float_t xyzMax[3]; | |
977 | Float_t xyzOrig[3]; | |
978 | Float_t xyzBoxd[3]; | |
979 | ||
980 | Char_t cTagV[5]; | |
981 | Char_t cTagM[5]; | |
982 | ||
983 | for (Int_t i = 0; i < 3; i++) { | |
984 | xyzMin[i] = +9999; xyzMax[i] = -9999; | |
985 | } | |
986 | ||
987 | for (Int_t i = 0; i < 3; i++) { | |
988 | ||
989 | xyzMin[i] = TMath::Min(xyzMin[i],fChamberUAorig[iDet][i]-fChamberUAboxd[iDet][i]); | |
990 | xyzMax[i] = TMath::Max(xyzMax[i],fChamberUAorig[iDet][i]+fChamberUAboxd[iDet][i]); | |
991 | ||
992 | xyzMin[i] = TMath::Min(xyzMin[i],fChamberUDorig[iDet][i]-fChamberUDboxd[iDet][i]); | |
993 | xyzMax[i] = TMath::Max(xyzMax[i],fChamberUDorig[iDet][i]+fChamberUDboxd[iDet][i]); | |
994 | ||
995 | xyzMin[i] = TMath::Min(xyzMin[i],fChamberUForig[iDet][i]-fChamberUFboxd[iDet][i]); | |
996 | xyzMax[i] = TMath::Max(xyzMax[i],fChamberUForig[iDet][i]+fChamberUFboxd[iDet][i]); | |
997 | ||
998 | // CBL | |
999 | //if (iplan < (kNplan-1)) { | |
1000 | xyzMin[i] = TMath::Min(xyzMin[i],fChamberUUorig[iDet][i]-fChamberUUboxd[iDet][i]); | |
1001 | xyzMax[i] = TMath::Max(xyzMax[i],fChamberUUorig[iDet][i]+fChamberUUboxd[iDet][i]); | |
1002 | //} | |
1003 | ||
1004 | xyzOrig[i] = 0.5*(xyzMax[i]+xyzMin[i]); | |
1005 | xyzBoxd[i] = 0.5*(xyzMax[i]-xyzMin[i]); | |
1006 | ||
1007 | } | |
1008 | ||
1009 | sprintf(cTagM,"UT%02d",iDet); | |
1010 | ||
1011 | gMC->Gsvolu(cTagM,"BOX ",idtmed[1302-1],xyzBoxd,kNparCha); | |
1012 | ||
1013 | sprintf(cTagV,"UA%02d",iDet); | |
1014 | gMC->Gspos(cTagV,1,cTagM, | |
1015 | fChamberUAorig[iDet][0]-xyzOrig[0], | |
1016 | fChamberUAorig[iDet][1]-xyzOrig[1], | |
1017 | fChamberUAorig[iDet][2]-xyzOrig[2], | |
1018 | 0,"ONLY"); | |
1019 | ||
1020 | sprintf(cTagV,"UD%02d",iDet); | |
1021 | gMC->Gspos(cTagV,1,cTagM, | |
1022 | fChamberUDorig[iDet][0]-xyzOrig[0], | |
1023 | fChamberUDorig[iDet][1]-xyzOrig[1], | |
1024 | fChamberUDorig[iDet][2]-xyzOrig[2], | |
1025 | 0,"ONLY"); | |
1026 | ||
1027 | sprintf(cTagV,"UF%02d",iDet); | |
1028 | gMC->Gspos(cTagV,1,cTagM, | |
1029 | fChamberUForig[iDet][0]-xyzOrig[0], | |
1030 | fChamberUForig[iDet][1]-xyzOrig[1], | |
1031 | fChamberUForig[iDet][2]-xyzOrig[2], | |
1032 | 0,"ONLY"); | |
1033 | ||
1034 | // CBL | |
1035 | //if (iplan < (kNplan-1)) { | |
1036 | sprintf(cTagV,"UU%02d",iDet); | |
1037 | gMC->Gspos(cTagV,1,cTagM, | |
1038 | fChamberUUorig[iDet][0]-xyzOrig[0], | |
1039 | fChamberUUorig[iDet][1]-xyzOrig[1], | |
1040 | fChamberUUorig[iDet][2]-xyzOrig[2], | |
1041 | 0,"ONLY"); | |
1042 | ||
1043 | // } | |
1044 | ||
1045 | sprintf(cTagV,"UT%02d",iDet); | |
1046 | gMC->Gspos(cTagV,1,"UTI1",xyzOrig[0],xyzOrig[1],xyzOrig[2],0,"ONLY"); | |
f7336fa3 | 1047 | |
1048 | } | |
1049 | ||
1050 | //_____________________________________________________________________________ | |
a5cadd36 | 1051 | Bool_t AliTRDgeometry::Local2Global(Int_t idet, Double_t *local |
dde59437 | 1052 | , Double_t *global) const |
f7336fa3 | 1053 | { |
1054 | // | |
1055 | // Converts local pad-coordinates (row,col,time) into | |
1056 | // global ALICE reference frame coordinates (x,y,z) | |
1057 | // | |
1058 | ||
793ff80c | 1059 | Int_t icham = GetChamber(idet); // Chamber info (0-4) |
1060 | Int_t isect = GetSector(idet); // Sector info (0-17) | |
1061 | Int_t iplan = GetPlane(idet); // Plane info (0-5) | |
f7336fa3 | 1062 | |
dde59437 | 1063 | return Local2Global(iplan,icham,isect,local,global); |
f7336fa3 | 1064 | |
1065 | } | |
1066 | ||
1067 | //_____________________________________________________________________________ | |
1068 | Bool_t AliTRDgeometry::Local2Global(Int_t iplan, Int_t icham, Int_t isect | |
dde59437 | 1069 | , Double_t *local, Double_t *global) const |
f7336fa3 | 1070 | { |
1071 | // | |
1072 | // Converts local pad-coordinates (row,col,time) into | |
1073 | // global ALICE reference frame coordinates (x,y,z) | |
1074 | // | |
1075 | ||
3551db50 | 1076 | AliTRDCommonParam* commonParam = AliTRDCommonParam::Instance(); |
1077 | if (!commonParam) | |
1078 | return kFALSE; | |
5443e65e | 1079 | |
3551db50 | 1080 | AliTRDcalibDB* calibration = AliTRDcalibDB::Instance(); |
1081 | if (!calibration) | |
1082 | return kFALSE; | |
1083 | ||
1084 | AliTRDpadPlane *padPlane = commonParam->GetPadPlane(iplan,icham); | |
f7336fa3 | 1085 | |
a5cadd36 | 1086 | // calculate (x,y,z) position in rotated chamber |
1087 | Int_t row = ((Int_t) local[0]); | |
1088 | Int_t col = ((Int_t) local[1]); | |
1089 | Float_t timeSlice = local[2] + 0.5; | |
3551db50 | 1090 | Float_t time0 = GetTime0(iplan); |
f7336fa3 | 1091 | |
7754cd1f | 1092 | Int_t idet = GetDetector(iplan, icham, isect); |
1093 | ||
a5cadd36 | 1094 | Double_t rot[3]; |
7754cd1f | 1095 | rot[0] = time0 - (timeSlice - calibration->GetT0(idet, col, row)) |
1096 | * calibration->GetVdrift(idet, col, row)/calibration->GetSamplingFrequency(); | |
a5cadd36 | 1097 | rot[1] = padPlane->GetColPos(col) - 0.5 * padPlane->GetColSize(col); |
1098 | rot[2] = padPlane->GetRowPos(row) - 0.5 * padPlane->GetRowSize(row); | |
f7336fa3 | 1099 | |
1100 | // Rotate back to original position | |
1101 | return RotateBack(idet,rot,global); | |
1102 | ||
1103 | } | |
1104 | ||
3d7b6a24 | 1105 | //_____________________________________________________________________________ |
a5cadd36 | 1106 | Bool_t AliTRDgeometry::Global2Local(Int_t mode, Double_t *local, Double_t *global |
bd0f8685 | 1107 | , Int_t* index) const |
3d7b6a24 | 1108 | { |
1109 | // | |
1110 | // Converts local pad-coordinates (row,col,time) into | |
1111 | // global ALICE reference frame coordinates (x,y,z) | |
1112 | // | |
e0d47c25 | 1113 | // index[0] = plane number |
1114 | // index[1] = chamber number | |
1115 | // index[2] = sector number | |
3d7b6a24 | 1116 | // |
1117 | // mode=0 - local coordinate in y, z, x - rotated global | |
1118 | // mode=2 - local coordinate in pad, and pad row, x - rotated global | |
1119 | // | |
e0d47c25 | 1120 | |
bd0f8685 | 1121 | Int_t idet = GetDetector(index[0],index[1],index[2]); // Detector number |
b4a9cd27 | 1122 | RotateBack(idet,global,local); |
bd0f8685 | 1123 | if (mode == 0) return kTRUE; |
3d7b6a24 | 1124 | |
1125 | return kTRUE; | |
1126 | ||
1127 | } | |
1128 | ||
a5cadd36 | 1129 | //_____________________________________________________________________________ |
3551db50 | 1130 | Bool_t AliTRDgeometry::Global2Detector(Double_t global[3], Int_t index[3]) |
3d7b6a24 | 1131 | { |
1132 | // | |
b4a9cd27 | 1133 | // Find detector for given global point - Ideal geometry |
1134 | // | |
1135 | // | |
e0d47c25 | 1136 | // input = global position |
1137 | // output = index | |
1138 | // index[0] = plane number | |
1139 | // index[1] = chamber number | |
1140 | // index[2] = sector number | |
3d7b6a24 | 1141 | // |
bd0f8685 | 1142 | |
3d7b6a24 | 1143 | // |
b4a9cd27 | 1144 | // Find sector |
1145 | // | |
1146 | Float_t fi = TMath::ATan2(global[1],global[0]); | |
bd0f8685 | 1147 | if (fi < 0) { |
1148 | fi += 2*TMath::Pi(); | |
1149 | } | |
1150 | index[2] = fgkNsect - 1 - TMath::Nint((fi - GetAlpha()/2.)/GetAlpha()); | |
1151 | ||
3d7b6a24 | 1152 | // |
b4a9cd27 | 1153 | // Find plane |
3d7b6a24 | 1154 | // |
1155 | Float_t locx = global[0] * fRotA11[index[2]] + global[1] * fRotA12[index[2]]; | |
1156 | index[0] = 0; | |
3551db50 | 1157 | Float_t max = locx - GetTime0(0); |
3d7b6a24 | 1158 | for (Int_t iplane=1; iplane<fgkNplan;iplane++){ |
3551db50 | 1159 | Float_t dist = TMath::Abs(locx - GetTime0(iplane)); |
3d7b6a24 | 1160 | if (dist < max){ |
1161 | index[0] = iplane; | |
1162 | max = dist; | |
1163 | } | |
1164 | } | |
bd0f8685 | 1165 | |
b4a9cd27 | 1166 | // |
1167 | // Find chamber | |
1168 | // | |
1169 | if (TMath::Abs(global[2]) < 0.5*GetChamberLength(index[0],2)){ | |
1170 | index[1]=2; | |
bd0f8685 | 1171 | } |
1172 | else{ | |
b4a9cd27 | 1173 | Double_t localZ = global[2]; |
bd0f8685 | 1174 | if (global[2] > 0){ |
b4a9cd27 | 1175 | localZ -= 0.5*(GetChamberLength(index[0],2)+GetChamberLength(index[0],1)); |
1176 | index[1] = (TMath::Abs(localZ) < 0.5*GetChamberLength(index[0],3)) ? 1:0; | |
1177 | } | |
1178 | else{ | |
1179 | localZ += 0.5*(GetChamberLength(index[0],2)+GetChamberLength(index[0],3)); | |
1180 | index[1] = (TMath::Abs(localZ) < 0.5*GetChamberLength(index[0],1)) ? 3:4; | |
1181 | } | |
1182 | } | |
bd0f8685 | 1183 | |
3d7b6a24 | 1184 | return kTRUE; |
3d7b6a24 | 1185 | |
bd0f8685 | 1186 | } |
3d7b6a24 | 1187 | |
f7336fa3 | 1188 | //_____________________________________________________________________________ |
a5cadd36 | 1189 | Bool_t AliTRDgeometry::Rotate(Int_t d, Double_t *pos, Double_t *rot) const |
f7336fa3 | 1190 | { |
1191 | // | |
1192 | // Rotates all chambers in the position of sector 0 and transforms | |
1193 | // the coordinates in the ALICE restframe <pos> into the | |
1194 | // corresponding local frame <rot>. | |
1195 | // | |
1196 | ||
793ff80c | 1197 | Int_t sector = GetSector(d); |
f7336fa3 | 1198 | |
793ff80c | 1199 | rot[0] = pos[0] * fRotA11[sector] + pos[1] * fRotA12[sector]; |
1200 | rot[1] = -pos[0] * fRotA21[sector] + pos[1] * fRotA22[sector]; | |
f7336fa3 | 1201 | rot[2] = pos[2]; |
1202 | ||
1203 | return kTRUE; | |
1204 | ||
1205 | } | |
1206 | ||
1207 | //_____________________________________________________________________________ | |
a5cadd36 | 1208 | Bool_t AliTRDgeometry::RotateBack(Int_t d, Double_t *rot, Double_t *pos) const |
f7336fa3 | 1209 | { |
1210 | // | |
1211 | // Rotates a chambers from the position of sector 0 into its | |
1212 | // original position and transforms the corresponding local frame | |
1213 | // coordinates <rot> into the coordinates of the ALICE restframe <pos>. | |
1214 | // | |
1215 | ||
793ff80c | 1216 | Int_t sector = GetSector(d); |
f7336fa3 | 1217 | |
793ff80c | 1218 | pos[0] = rot[0] * fRotB11[sector] + rot[1] * fRotB12[sector]; |
1219 | pos[1] = -rot[0] * fRotB21[sector] + rot[1] * fRotB22[sector]; | |
6f1e466d | 1220 | pos[2] = rot[2]; |
f7336fa3 | 1221 | |
1222 | return kTRUE; | |
1223 | ||
1224 | } | |
1225 | ||
1226 | //_____________________________________________________________________________ | |
3551db50 | 1227 | Int_t AliTRDgeometry::GetDetectorSec(Int_t p, Int_t c) |
0a770ac9 | 1228 | { |
1229 | // | |
1230 | // Convert plane / chamber into detector number for one single sector | |
1231 | // | |
1232 | ||
1233 | return (p + c * fgkNplan); | |
1234 | ||
1235 | } | |
1236 | ||
1237 | //_____________________________________________________________________________ | |
3551db50 | 1238 | Int_t AliTRDgeometry::GetDetector(Int_t p, Int_t c, Int_t s) |
f7336fa3 | 1239 | { |
1240 | // | |
1241 | // Convert plane / chamber / sector into detector number | |
1242 | // | |
1243 | ||
793ff80c | 1244 | return (p + c * fgkNplan + s * fgkNplan * fgkNcham); |
f7336fa3 | 1245 | |
1246 | } | |
1247 | ||
1248 | //_____________________________________________________________________________ | |
afc51ac2 | 1249 | Int_t AliTRDgeometry::GetPlane(Int_t d) const |
f7336fa3 | 1250 | { |
1251 | // | |
1252 | // Reconstruct the plane number from the detector number | |
1253 | // | |
1254 | ||
793ff80c | 1255 | return ((Int_t) (d % fgkNplan)); |
f7336fa3 | 1256 | |
1257 | } | |
1258 | ||
1259 | //_____________________________________________________________________________ | |
afc51ac2 | 1260 | Int_t AliTRDgeometry::GetChamber(Int_t d) const |
f7336fa3 | 1261 | { |
1262 | // | |
1263 | // Reconstruct the chamber number from the detector number | |
1264 | // | |
1265 | ||
793ff80c | 1266 | return ((Int_t) (d % (fgkNplan * fgkNcham)) / fgkNplan); |
f7336fa3 | 1267 | |
1268 | } | |
1269 | ||
1270 | //_____________________________________________________________________________ | |
afc51ac2 | 1271 | Int_t AliTRDgeometry::GetSector(Int_t d) const |
f7336fa3 | 1272 | { |
1273 | // | |
1274 | // Reconstruct the sector number from the detector number | |
1275 | // | |
1276 | ||
793ff80c | 1277 | return ((Int_t) (d / (fgkNplan * fgkNcham))); |
f7336fa3 | 1278 | |
1279 | } | |
1280 | ||
bdbb05bb | 1281 | //_____________________________________________________________________________ |
1282 | AliTRDgeometry* AliTRDgeometry::GetGeometry(AliRunLoader* runLoader) | |
1283 | { | |
1284 | // | |
1285 | // load the geometry from the galice file | |
1286 | // | |
1287 | ||
1288 | if (!runLoader) runLoader = AliRunLoader::GetRunLoader(); | |
1289 | if (!runLoader) { | |
1290 | ::Error("AliTRDgeometry::GetGeometry", "No run loader"); | |
1291 | return NULL; | |
1292 | } | |
1293 | ||
1294 | TDirectory* saveDir = gDirectory; | |
1295 | runLoader->CdGAFile(); | |
1296 | ||
ecb36af7 | 1297 | // Try from the galice.root file |
bdbb05bb | 1298 | AliTRDgeometry* geom = (AliTRDgeometry*) gDirectory->Get("TRDgeometry"); |
ecb36af7 | 1299 | |
1300 | if (!geom) { | |
1301 | // It is not in the file, try to get it from gAlice, | |
1302 | // which corresponds to the run loader | |
1303 | AliTRD * trd = (AliTRD*)runLoader->GetAliRun()->GetDetector("TRD"); | |
1304 | geom = trd->GetGeometry(); | |
1305 | } | |
bdbb05bb | 1306 | if (!geom) ::Error("AliTRDgeometry::GetGeometry", "Geometry not found"); |
1307 | ||
1308 | saveDir->cd(); | |
1309 | return geom; | |
b4a9cd27 | 1310 | |
bd0f8685 | 1311 | } |
b4a9cd27 | 1312 | |
1313 | //_____________________________________________________________________________ | |
bd0f8685 | 1314 | Bool_t AliTRDgeometry::ReadGeoMatrices() |
1315 | { | |
b4a9cd27 | 1316 | // |
1317 | // Read geo matrices from current gGeoManager for each TRD sector | |
1318 | // | |
1319 | ||
b4a9cd27 | 1320 | if (!gGeoManager) return kFALSE; |
1321 | fMatrixArray = new TObjArray(kNdet); | |
1322 | fMatrixCorrectionArray = new TObjArray(kNdet); | |
1323 | fMatrixGeo = new TObjArray(kNdet); | |
1324 | AliAlignObjAngles o; | |
bd0f8685 | 1325 | |
b4a9cd27 | 1326 | for (Int_t iLayer = AliAlignObj::kTRD1; iLayer <= AliAlignObj::kTRD6; iLayer++) { |
1327 | for (Int_t iModule = 0; iModule < AliAlignObj::LayerSize(iLayer); iModule++) { | |
1328 | UShort_t volid = AliAlignObj::LayerToVolUID(iLayer,iModule); | |
1329 | const char *path = AliAlignObj::GetVolPath(volid); | |
1330 | if (!gGeoManager->cd(path)) return kFALSE; | |
1331 | TGeoHMatrix* m = gGeoManager->GetCurrentMatrix(); | |
1332 | Int_t iLayerTRD = iLayer-AliAlignObj::kTRD1; | |
1333 | Int_t isector = Nsect()-1-(iModule/Ncham()); | |
1334 | Int_t ichamber = Ncham()-1-(iModule%Ncham()); | |
1335 | Int_t lid = GetDetector(iLayerTRD,ichamber,isector); | |
bd0f8685 | 1336 | |
b4a9cd27 | 1337 | // |
1338 | // local geo system z-x-y to x-y--z | |
1339 | // | |
1340 | fMatrixGeo->AddAt(new TGeoHMatrix(*m),lid); | |
1341 | ||
1342 | TGeoRotation mchange; | |
1343 | mchange.RotateY(90); mchange.RotateX(90); | |
bd0f8685 | 1344 | |
b4a9cd27 | 1345 | TGeoHMatrix gMatrix(mchange.Inverse()); |
1346 | gMatrix.MultiplyLeft(m); | |
1347 | fMatrixArray->AddAt(new TGeoHMatrix(gMatrix),lid); | |
bd0f8685 | 1348 | |
b4a9cd27 | 1349 | // |
1350 | // Cluster transformation matrix | |
1351 | // | |
1352 | TGeoHMatrix rotMatrix(mchange.Inverse()); | |
1353 | rotMatrix.MultiplyLeft(m); | |
1354 | Double_t sectorAngle = 20.*(isector%18)+10; | |
1355 | TGeoHMatrix rotSector; | |
1356 | rotSector.RotateZ(sectorAngle); | |
1357 | rotMatrix.MultiplyLeft(&rotSector); | |
bd0f8685 | 1358 | |
b4a9cd27 | 1359 | fMatrixCorrectionArray->AddAt(new TGeoHMatrix(rotMatrix),lid); |
bd0f8685 | 1360 | |
b4a9cd27 | 1361 | } |
1362 | } | |
bd0f8685 | 1363 | |
b4a9cd27 | 1364 | return kTRUE; |
b4a9cd27 | 1365 | |
bd0f8685 | 1366 | } |
b4a9cd27 | 1367 |