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