]> git.uio.no Git - u/mrichter/AliRoot.git/blame - TRD/AliTRDgeometryHole.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Pyhtia comparison extended
[u/mrichter/AliRoot.git] / TRD / AliTRDgeometryHole.cxx
CommitLineData
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
73ae7b59 16/*
17$Log$
f2649dd2 18Revision 1.11 2003/09/18 09:06:07 cblume
19Geometry update, Removal of compiler warnings
20
73ae7b59 21Revision 1.9 2002/11/21 22:38:47 alibrary
22Removing AliMC and AliMCProcess
23
24Revision 1.8 2002/10/31 17:45:35 cblume
25New chamber geometry
26
27Revision 1.7 2002/02/11 14:21:16 cblume
28Update of the geometry. Get rid of MANY
29
30Revision 1.6 2001/05/11 07:56:12 hristov
31Consistent declarations needed on Alpha
32
33Revision 1.5 2001/02/14 18:22:26 cblume
34Change in the geometry of the padplane
35
36Revision 1.4 2000/11/01 14:53:21 cblume
37Merge with TRD-develop
38
39Revision 1.1.4.4 2000/10/15 23:40:01 cblume
40Remove AliTRDconst
41
42Revision 1.1.4.3 2000/10/06 16:49:46 cblume
43Made Getters const
44
45Revision 1.1.4.2 2000/10/04 16:34:58 cblume
46Replace include files by forward declarations
47
48Revision 1.1.4.1 2000/09/22 14:43:41 cblume
49Allow the pad/timebin-dimensions to be changed after initialization
50
51Revision 1.3 2000/10/02 21:28:19 fca
52Removal of useless dependecies via forward declarations
53
54Revision 1.2 2000/06/08 18:32:58 cblume
55Make code compliant to coding conventions
56
57Revision 1.1 2000/02/28 19:01:42 cblume
58Add new TRD classes
59
60*/
f7336fa3 61
62///////////////////////////////////////////////////////////////////////////////
63// //
64// TRD geometry with holes //
65// //
66///////////////////////////////////////////////////////////////////////////////
67
73ae7b59 68#include "TVirtualMC.h"
69
793ff80c 70#include "AliTRDgeometryHole.h"
71
f7336fa3 72ClassImp(AliTRDgeometryHole)
73
74//_____________________________________________________________________________
75AliTRDgeometryHole::AliTRDgeometryHole():AliTRDgeometry()
76{
77 //
78 // AliTRDgeometryHole default constructor
79 //
80
81 Init();
82
83}
84
85//_____________________________________________________________________________
86AliTRDgeometryHole::~AliTRDgeometryHole()
87{
8230f242 88 //
89 // AliTRDgeometryHole destructor
90 //
f7336fa3 91
92}
93
f2649dd2 94Bool_t AliTRDgeometryHole::IsHole(Int_t iplan, Int_t icham, Int_t isec) const
95{
96 // Position of Holes for PHOS (P) and RICH (R) starting at 6h
97 // P P P - - R R R - - - - - - - - P P
98 //Int_t cham[18] = {1, 1, 1, 0, 0, 2, 2, 2, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1};
99 Int_t cham[18] = {2, 2, 2, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 0, 0 }; // did we want this?
100
101 if ((cham[17-isec]==1) && (fClengthPH[iplan][icham]<0.00001)) return kTRUE;
102 if ((cham[17-isec]==2) &&(fClengthRH[iplan][icham]<0.000001)) return kTRUE;
103 return kFALSE;
104
105}
106
f7336fa3 107//_____________________________________________________________________________
108void AliTRDgeometryHole::Init()
109{
110 //
111 // Initializes the geometry parameter
112 //
113
114 Int_t iplan;
0a770ac9 115 Int_t icham;
116
117 // The outer lengths of the chambers for the sectors with holes for the PHOS
8737e16f 118 Float_t lengthPH[kNplan][kNcham] = { { 0.0, 0.0, 0.0, 117.0, 124.0 }
0a770ac9 119 , { 0.0, 0.0, 0.0, 124.0, 131.0 }
8737e16f 120 , { 0.0, 0.0, 0.0, 131.0, 138.0 }
121 , { 0.0, 0.0, 0.0, 138.0, 145.0 }
122 , { 0.0, 0.0, 0.0, 140.0, 147.0 }
123 , { 0.0, 0.0, 0.0, 140.0, 147.0 } };
0a770ac9 124
125 // The outer lengths of the chambers for the sectors with holes for the RICH
8737e16f 126 Float_t lengthRH[kNplan][kNcham] = { { 0.0, 0.0, 0.0, 0.0, 87.5 }
0a770ac9 127 , { 0.0, 0.0, 0.0, 0.0, 101.5 }
8737e16f 128 , { 0.0, 0.0, 0.0, 0.0, 115.5 }
129 , { 0.0, 0.0, 0.0, 0.0, 129.5 }
130 , { 0.0, 0.0, 0.0, 0.0, 133.5 }
131 , { 0.0, 0.0, 0.0, 0.0, 133.5 } };
0a770ac9 132
133 for (icham = 0; icham < kNcham; icham++) {
134 for (iplan = 0; iplan < kNplan; iplan++) {
135 fClengthPH[iplan][icham] = lengthPH[iplan][icham];
136 fClengthRH[iplan][icham] = lengthRH[iplan][icham];
f7336fa3 137 }
f7336fa3 138 }
139
140}
141
142//_____________________________________________________________________________
0a770ac9 143void AliTRDgeometryHole::CreateGeometry(Int_t *idtmed)
71d9fa7b 144{
145 //
0a770ac9 146 // Create the TRD geometry with holes
71d9fa7b 147 //
0a770ac9 148 // Names of the TRD volumina (xx = detector number):
149 //
150 // Lower part of the readout chambers (gas volume + radiator)
151 //
152 // UAxx Aluminum frames (Al)
153 // UBxx G10 frames (C)
154 // UCxx Inner volumes (Air)
155 //
156 // Upper part of the readout chambers (readout plane + fee)
f7336fa3 157 //
0a770ac9 158 // UDxx G10 frames (C)
159 // UExx Inner volumes of the G10 (Air)
160 // UFxx Aluminum frames (Al)
161 // UGxx Inner volumes of the Al (Air)
162 //
163 // Inner material layers
164 //
165 // UHxx Radiator (Rohacell)
166 // UIxx Entrance window (Mylar)
167 // UJxx Drift volume (Xe/CO2)
168 // UKxx Amplification volume (Xe/CO2)
169 // ULxx Pad plane (Cu)
170 // UMxx Support structure (Rohacell)
171 // UNxx FEE + signal lines (Cu)
172 // UOxx Cooling device (Al)
173 // UPxx Cooling device (Water)
f7336fa3 174 //
175
0a770ac9 176 const Int_t kNdet = kNplan * kNcham;
f7336fa3 177
8230f242 178 const Int_t kNparTrd = 4;
179 const Int_t kNparCha = 3;
0a770ac9 180
181 Float_t xpos, ypos, zpos;
f7336fa3 182
8230f242 183 Float_t parTrd[kNparTrd];
184 Float_t parCha[kNparCha];
f7336fa3 185
0a770ac9 186 Char_t cTagV[5];
187 Char_t cTagM[5];
f7336fa3 188
189 AliTRDgeometry::CreateGeometry(idtmed);
190
0a770ac9 191 // The TRD mother volume for one sector (Air), full length in z-direction
793ff80c 192 parTrd[0] = fgkSwidth1/2.;
193 parTrd[1] = fgkSwidth2/2.;
194 parTrd[2] = fgkSlenTR1/2.;
195 parTrd[3] = fgkSheight/2.;
0a770ac9 196 gMC->Gsvolu("UTR1","TRD1",idtmed[1302-1],parTrd,kNparTrd);
197
198 // The TRD mother volume for one sector (Air), leaving hole for PHOS
793ff80c 199 parTrd[0] = fgkSwidth1/2.;
200 parTrd[1] = fgkSwidth2/2.;
201 parTrd[2] = fgkSlenTR2/2.;
202 parTrd[3] = fgkSheight/2.;
0a770ac9 203 gMC->Gsvolu("UTR2","TRD1",idtmed[1302-1],parTrd,kNparTrd);
f7336fa3 204
0a770ac9 205 // The TRD mother volume for one sector (Air), leaving hole for RICH
793ff80c 206 parTrd[0] = fgkSwidth1/2.;
207 parTrd[1] = fgkSwidth2/2.;
208 parTrd[2] = fgkSlenTR3/2.;
209 parTrd[3] = fgkSheight/2.;
0a770ac9 210 gMC->Gsvolu("UTR3","TRD1",idtmed[1302-1],parTrd,kNparTrd);
211
212 for (Int_t icham = 0; icham < kNcham; icham++) {
213 for (Int_t iplan = 0; iplan < kNplan; iplan++) {
214
215 Int_t iDet = GetDetectorSec(iplan,icham);
216
217 // The lower part of the readout chambers (gas volume + radiator)
218 // The aluminum frames
219 sprintf(cTagV,"UA%02d",iDet);
220 parCha[0] = fCwidth[iplan]/2.;
221 parCha[1] = fClength[iplan][icham]/2. - fgkHspace/2.;
222 parCha[2] = fgkCraH/2. + fgkCdrH/2.;
223 gMC->Gsvolu(cTagV,"BOX ",idtmed[1301-1],parCha,kNparCha);
224 // The G10 frames
225 sprintf(cTagV,"UB%02d",iDet);
226 parCha[0] = fCwidth[iplan]/2. - fgkCalT;
227 parCha[1] = -1.;
228 parCha[2] = -1.;
229 gMC->Gsvolu(cTagV,"BOX ",idtmed[1307-1],parCha,kNparCha);
230 // The inner part (air)
231 sprintf(cTagV,"UC%02d",iDet);
232 parCha[0] = fCwidth[iplan]/2. - fgkCalT - fgkCclsT;
233 parCha[1] = fClength[iplan][icham]/2. - fgkHspace/2.- fgkCclfT;
234 parCha[2] = -1.;
235 gMC->Gsvolu(cTagV,"BOX ",idtmed[1302-1],parCha,kNparCha);
236 if (fClengthPH[iplan][icham] > 0.0) {
237 // The aluminum frames
238 sprintf(cTagV,"UA%02d",iDet+kNdet);
239 parCha[0] = fCwidth[iplan]/2.;
240 parCha[1] = fClengthPH[iplan][icham]/2. - fgkHspace/2.;
241 parCha[2] = fgkCraH/2. + fgkCdrH/2.;
242 gMC->Gsvolu(cTagV,"BOX ",idtmed[1301-1],parCha,kNparCha);
243 // The G10 frames
244 sprintf(cTagV,"UB%02d",iDet+kNdet);
245 parCha[0] = fCwidth[iplan]/2. - fgkCalT;
246 parCha[1] = -1.;
247 parCha[2] = -1.;
248 gMC->Gsvolu(cTagV,"BOX ",idtmed[1307-1],parCha,kNparCha);
249 // The inner part (air)
250 sprintf(cTagV,"UC%02d",iDet+kNdet);
251 parCha[0] = fCwidth[iplan]/2. - fgkCalT - fgkCclsT;
252 parCha[1] = fClengthPH[iplan][icham]/2. - fgkHspace/2.- fgkCclfT;
253 parCha[2] = -1.;
254 gMC->Gsvolu(cTagV,"BOX ",idtmed[1302-1],parCha,kNparCha);
255 }
256 if (fClengthRH[iplan][icham] > 0.0) {
257 // The aluminum frames
258 sprintf(cTagV,"UA%02d",iDet+2*kNdet);
259 parCha[0] = fCwidth[iplan]/2.;
260 parCha[1] = fClengthRH[iplan][icham]/2. - fgkHspace/2.;
261 parCha[2] = fgkCraH/2. + fgkCdrH/2.;
262 gMC->Gsvolu(cTagV,"BOX ",idtmed[1301-1],parCha,kNparCha);
263 // The G10 frames
264 sprintf(cTagV,"UB%02d",iDet+2*kNdet);
265 parCha[0] = fCwidth[iplan]/2. - fgkCalT;
266 parCha[1] = -1.;
267 parCha[2] = -1.;
268 gMC->Gsvolu(cTagV,"BOX ",idtmed[1307-1],parCha,kNparCha);
269 // The inner part (air)
270 sprintf(cTagV,"UC%02d",iDet+2*kNdet);
271 parCha[0] = fCwidth[iplan]/2. - fgkCalT - fgkCclsT;
272 parCha[1] = fClengthRH[iplan][icham]/2. - fgkHspace/2.- fgkCclfT;
273 parCha[2] = -1.;
274 gMC->Gsvolu(cTagV,"BOX ",idtmed[1302-1],parCha,kNparCha);
275 }
276
277 // The upper part of the readout chambers (readout plane + fee)
278 // The G10 frames
279 sprintf(cTagV,"UD%02d",iDet);
280 parCha[0] = fCwidth[iplan]/2. + fgkCroW;
281 parCha[1] = fClength[iplan][icham]/2. - fgkHspace/2.;
282 parCha[2] = fgkCamH/2.;
283 gMC->Gsvolu(cTagV,"BOX ",idtmed[1307-1],parCha,kNparCha);
284 // The inner part of the G10 frame (air)
285 sprintf(cTagV,"UE%02d",iDet);
286 parCha[0] = fCwidth[iplan]/2. + fgkCroW - fgkCcuT;
287 parCha[1] = fClength[iplan][icham]/2. - fgkHspace/2.- fgkCcuT;
288 parCha[2] = -1.;
289 gMC->Gsvolu(cTagV,"BOX ",idtmed[1302-1],parCha,kNparCha);
290 // The aluminum frames
291 sprintf(cTagV,"UF%02d",iDet);
292 parCha[0] = fCwidth[iplan]/2. + fgkCroW;
293 parCha[1] = fClength[iplan][icham]/2. - fgkHspace/2.;
294 parCha[2] = fgkCroH/2.;
295 gMC->Gsvolu(cTagV,"BOX ",idtmed[1301-1],parCha,kNparCha);
296 // The inner part of the aluminum frames
297 sprintf(cTagV,"UG%02d",iDet);
298 parCha[0] = fCwidth[iplan]/2. + fgkCroW - fgkCauT;
299 parCha[1] = fClength[iplan][icham]/2. - fgkHspace/2.- fgkCauT;
300 parCha[2] = -1.;
301 gMC->Gsvolu(cTagV,"BOX ",idtmed[1302-1],parCha,kNparCha);
302 if (fClengthPH[iplan][icham] > 0.0) {
303 sprintf(cTagV,"UD%02d",iDet+kNdet);
304 parCha[0] = fCwidth[iplan]/2. + fgkCroW;
305 parCha[1] = fClengthPH[iplan][icham]/2. - fgkHspace/2.;
306 parCha[2] = fgkCamH/2.;
307 gMC->Gsvolu(cTagV,"BOX ",idtmed[1307-1],parCha,kNparCha);
308 // The inner part of the G10 frame (air)
309 sprintf(cTagV,"UE%02d",iDet+kNdet);
310 parCha[0] = fCwidth[iplan]/2. + fgkCroW - fgkCcuT;
311 parCha[1] = fClengthPH[iplan][icham]/2. - fgkHspace/2.- fgkCcuT;
312 parCha[2] = -1.;
313 gMC->Gsvolu(cTagV,"BOX ",idtmed[1302-1],parCha,kNparCha);
314 // The aluminum frames
315 sprintf(cTagV,"UF%02d",iDet+kNdet);
316 parCha[0] = fCwidth[iplan]/2. + fgkCroW;
317 parCha[1] = fClengthPH[iplan][icham]/2. - fgkHspace/2.;
318 parCha[2] = fgkCroH/2.;
319 gMC->Gsvolu(cTagV,"BOX ",idtmed[1301-1],parCha,kNparCha);
320 // The inner part of the aluminum frames
321 sprintf(cTagV,"UG%02d",iDet+kNdet);
322 parCha[0] = fCwidth[iplan]/2. + fgkCroW - fgkCauT;
323 parCha[1] = fClengthPH[iplan][icham]/2. - fgkHspace/2.- fgkCauT;
324 parCha[2] = -1.;
325 gMC->Gsvolu(cTagV,"BOX ",idtmed[1302-1],parCha,kNparCha);
326 }
327 if (fClengthRH[iplan][icham] > 0.0) {
328 sprintf(cTagV,"UD%02d",iDet+2*kNdet);
329 parCha[0] = fCwidth[iplan]/2. + fgkCroW;
330 parCha[1] = fClengthRH[iplan][icham]/2. - fgkHspace/2.;
331 parCha[2] = fgkCamH/2.;
332 gMC->Gsvolu(cTagV,"BOX ",idtmed[1307-1],parCha,kNparCha);
333 // The inner part of the G10 frame (air)
334 sprintf(cTagV,"UE%02d",iDet+2*kNdet);
335 parCha[0] = fCwidth[iplan]/2. + fgkCroW - fgkCcuT;
336 parCha[1] = fClengthRH[iplan][icham]/2. - fgkHspace/2.- fgkCcuT;
337 parCha[2] = -1.;
338 gMC->Gsvolu(cTagV,"BOX ",idtmed[1302-1],parCha,kNparCha);
339 // The aluminum frames
340 sprintf(cTagV,"UF%02d",iDet+2*kNdet);
341 parCha[0] = fCwidth[iplan]/2. + fgkCroW;
342 parCha[1] = fClengthRH[iplan][icham]/2. - fgkHspace/2.;
343 parCha[2] = fgkCroH/2.;
344 gMC->Gsvolu(cTagV,"BOX ",idtmed[1301-1],parCha,kNparCha);
345 // The inner part of the aluminum frames
346 sprintf(cTagV,"UG%02d",iDet+2*kNdet);
347 parCha[0] = fCwidth[iplan]/2. + fgkCroW - fgkCauT;
348 parCha[1] = fClengthRH[iplan][icham]/2. - fgkHspace/2.- fgkCauT;
349 parCha[2] = -1.;
350 gMC->Gsvolu(cTagV,"BOX ",idtmed[1302-1],parCha,kNparCha);
351 }
f7336fa3 352
0a770ac9 353 // The material layers inside the chambers
354 parCha[0] = -1.;
355 parCha[1] = -1.;
356 // Rohacell layer (radiator)
357 parCha[2] = fgkRaThick/2;
358 sprintf(cTagV,"UH%02d",iDet);
359 gMC->Gsvolu(cTagV,"BOX ",idtmed[1315-1],parCha,kNparCha);
360 // Mylar layer (entrance window + HV cathode)
361 parCha[2] = fgkMyThick/2;
362 sprintf(cTagV,"UI%02d",iDet);
363 gMC->Gsvolu(cTagV,"BOX ",idtmed[1308-1],parCha,kNparCha);
364 // Xe/Isobutane layer (drift volume)
365 parCha[2] = fgkDrThick/2.;
366 sprintf(cTagV,"UJ%02d",iDet);
367 gMC->Gsvolu(cTagV,"BOX ",idtmed[1309-1],parCha,kNparCha);
368 // Xe/Isobutane layer (amplification volume)
369 parCha[2] = fgkAmThick/2.;
370 sprintf(cTagV,"UK%02d",iDet);
371 gMC->Gsvolu(cTagV,"BOX ",idtmed[1309-1],parCha,kNparCha);
372 // Cu layer (pad plane)
373 parCha[2] = fgkCuThick/2;
374 sprintf(cTagV,"UL%02d",iDet);
375 gMC->Gsvolu(cTagV,"BOX ",idtmed[1305-1],parCha,kNparCha);
376 // G10 layer (support structure / honeycomb)
377 parCha[2] = fgkSuThick/2;
378 sprintf(cTagV,"UM%02d",iDet);
379 gMC->Gsvolu(cTagV,"BOX ",idtmed[1313-1],parCha,kNparCha);
380 // Cu layer (FEE + signal lines)
381 parCha[2] = fgkFeThick/2;
382 sprintf(cTagV,"UN%02d",iDet);
383 gMC->Gsvolu(cTagV,"BOX ",idtmed[1305-1],parCha,kNparCha);
384 // Al layer (cooling devices)
385 parCha[2] = fgkCoThick/2;
386 sprintf(cTagV,"UO%02d",iDet);
387 gMC->Gsvolu(cTagV,"BOX ",idtmed[1301-1],parCha,kNparCha);
388 // Water layer (cooling)
389 parCha[2] = fgkWaThick/2;
390 sprintf(cTagV,"UP%02d",iDet);
391 gMC->Gsvolu(cTagV,"BOX ",idtmed[1314-1],parCha,kNparCha);
392 if (fClengthPH[iplan][icham] > 0.0) {
393 // Rohacell layer (radiator)
394 parCha[2] = fgkRaThick/2;
395 sprintf(cTagV,"UH%02d",iDet+kNdet);
396 gMC->Gsvolu(cTagV,"BOX ",idtmed[1315-1],parCha,kNparCha);
397 // Mylar layer (entrance window + HV cathode)
398 parCha[2] = fgkMyThick/2;
399 sprintf(cTagV,"UI%02d",iDet+kNdet);
400 gMC->Gsvolu(cTagV,"BOX ",idtmed[1308-1],parCha,kNparCha);
401 // Xe/Isobutane layer (drift volume)
402 parCha[2] = fgkDrThick/2.;
403 sprintf(cTagV,"UJ%02d",iDet+kNdet);
404 gMC->Gsvolu(cTagV,"BOX ",idtmed[1309-1],parCha,kNparCha);
405 // Xe/Isobutane layer (amplification volume)
406 parCha[2] = fgkAmThick/2.;
407 sprintf(cTagV,"UK%02d",iDet+kNdet);
408 gMC->Gsvolu(cTagV,"BOX ",idtmed[1309-1],parCha,kNparCha);
409 // Cu layer (pad plane)
410 parCha[2] = fgkCuThick/2;
411 sprintf(cTagV,"UL%02d",iDet+kNdet);
412 gMC->Gsvolu(cTagV,"BOX ",idtmed[1305-1],parCha,kNparCha);
413 // G10 layer (support structure / honeycomb)
414 parCha[2] = fgkSuThick/2;
415 sprintf(cTagV,"UM%02d",iDet+kNdet);
416 gMC->Gsvolu(cTagV,"BOX ",idtmed[1313-1],parCha,kNparCha);
417 // Cu layer (FEE + signal lines)
418 parCha[2] = fgkFeThick/2;
419 sprintf(cTagV,"UN%02d",iDet+kNdet);
420 gMC->Gsvolu(cTagV,"BOX ",idtmed[1305-1],parCha,kNparCha);
421 // Al layer (cooling devices)
422 parCha[2] = fgkCoThick/2;
423 sprintf(cTagV,"UO%02d",iDet+kNdet);
424 gMC->Gsvolu(cTagV,"BOX ",idtmed[1301-1],parCha,kNparCha);
425 // Water layer (cooling)
426 parCha[2] = fgkWaThick/2;
427 sprintf(cTagV,"UP%02d",iDet+kNdet);
428 gMC->Gsvolu(cTagV,"BOX ",idtmed[1314-1],parCha,kNparCha);
429 }
430 if (fClengthRH[iplan][icham] > 0.0) {
431 // Rohacell layer (radiator)
432 parCha[2] = fgkRaThick/2;
433 sprintf(cTagV,"UH%02d",iDet+2*kNdet);
434 gMC->Gsvolu(cTagV,"BOX ",idtmed[1315-1],parCha,kNparCha);
435 // Mylar layer (entrance window + HV cathode)
436 parCha[2] = fgkMyThick/2;
437 sprintf(cTagV,"UI%02d",iDet+2*kNdet);
438 gMC->Gsvolu(cTagV,"BOX ",idtmed[1308-1],parCha,kNparCha);
439 // Xe/Isobutane layer (drift volume)
440 parCha[2] = fgkDrThick/2.;
441 sprintf(cTagV,"UJ%02d",iDet+2*kNdet);
442 gMC->Gsvolu(cTagV,"BOX ",idtmed[1309-1],parCha,kNparCha);
443 // Xe/Isobutane layer (amplification volume)
444 parCha[2] = fgkAmThick/2.;
445 sprintf(cTagV,"UK%02d",iDet+2*kNdet);
446 gMC->Gsvolu(cTagV,"BOX ",idtmed[1309-1],parCha,kNparCha);
447 // Cu layer (pad plane)
448 parCha[2] = fgkCuThick/2;
449 sprintf(cTagV,"UL%02d",iDet+2*kNdet);
450 gMC->Gsvolu(cTagV,"BOX ",idtmed[1305-1],parCha,kNparCha);
451 // G10 layer (support structure / honeycomb)
452 parCha[2] = fgkSuThick/2;
453 sprintf(cTagV,"UM%02d",iDet+2*kNdet);
454 gMC->Gsvolu(cTagV,"BOX ",idtmed[1313-1],parCha,kNparCha);
455 // Cu layer (FEE + signal lines)
456 parCha[2] = fgkFeThick/2;
457 sprintf(cTagV,"UN%02d",iDet+2*kNdet);
458 gMC->Gsvolu(cTagV,"BOX ",idtmed[1305-1],parCha,kNparCha);
459 // Al layer (cooling devices)
460 parCha[2] = fgkCoThick/2.;
461 sprintf(cTagV,"UO%02d",iDet+2*kNdet);
462 gMC->Gsvolu(cTagV,"BOX ",idtmed[1301-1],parCha,kNparCha);
463 // Water layer (cooling)
464 parCha[2] = fgkWaThick/2;
465 sprintf(cTagV,"UP%02d",iDet+2*kNdet);
466 gMC->Gsvolu(cTagV,"BOX ",idtmed[1314-1],parCha,kNparCha);
467 }
468
469 // Position the layers in the chambers
470 xpos = 0;
471 ypos = 0;
472 // Lower part
473 // Rohacell layer (radiator)
474 zpos = fgkRaZpos;
475 sprintf(cTagV,"UH%02d",iDet);
476 sprintf(cTagM,"UC%02d",iDet);
477 gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
478 // Mylar layer (entrance window + HV cathode)
479 zpos = fgkMyZpos;
480 sprintf(cTagV,"UI%02d",iDet);
481 sprintf(cTagM,"UC%02d",iDet);
482 gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
483 // Xe/Isobutane layer (drift volume)
484 zpos = fgkDrZpos;
485 sprintf(cTagV,"UJ%02d",iDet);
486 sprintf(cTagM,"UC%02d",iDet);
487 gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
488 // Upper part
489 // Xe/Isobutane layer (amplification volume)
490 zpos = fgkAmZpos;
491 sprintf(cTagV,"UK%02d",iDet);
492 sprintf(cTagM,"UE%02d",iDet);
493 gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
494 // Readout part
495 // Cu layer (pad plane)
496 zpos = fgkCuZpos;
497 sprintf(cTagV,"UL%02d",iDet);
498 sprintf(cTagM,"UG%02d",iDet);
499 gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
500 // G10 layer (support structure)
501 zpos = fgkSuZpos;
502 sprintf(cTagV,"UM%02d",iDet);
503 sprintf(cTagM,"UG%02d",iDet);
504 gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
505 // Cu layer (FEE + signal lines)
506 zpos = fgkFeZpos;
507 sprintf(cTagV,"UN%02d",iDet);
508 sprintf(cTagM,"UG%02d",iDet);
509 gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
510 // Al layer (cooling devices)
511 zpos = fgkCoZpos;
512 sprintf(cTagV,"UO%02d",iDet);
513 sprintf(cTagM,"UG%02d",iDet);
514 gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
515 // Water layer (cooling)
516 zpos = fgkWaZpos;
517 sprintf(cTagV,"UP%02d",iDet);
518 sprintf(cTagM,"UG%02d",iDet);
519 gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
520 if (fClengthPH[iplan][icham] > 0.0) {
521 // Lower part
522 // Rohacell layer (radiator)
523 zpos = fgkRaZpos;
524 sprintf(cTagV,"UH%02d",iDet+kNdet);
525 sprintf(cTagM,"UC%02d",iDet+kNdet);
526 gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
527 // Mylar layer (entrance window + HV cathode)
528 zpos = fgkMyZpos;
529 sprintf(cTagV,"UI%02d",iDet+kNdet);
530 sprintf(cTagM,"UC%02d",iDet+kNdet);
531 gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
532 // Xe/Isobutane layer (drift volume)
533 zpos = fgkDrZpos;
534 sprintf(cTagV,"UJ%02d",iDet+kNdet);
535 sprintf(cTagM,"UC%02d",iDet+kNdet);
536 gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
537 // Upper part
538 // Xe/Isobutane layer (amplification volume)
539 zpos = fgkAmZpos;
540 sprintf(cTagV,"UK%02d",iDet+kNdet);
541 sprintf(cTagM,"UE%02d",iDet+kNdet);
542 gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
543 // Readout part
544 // Cu layer (pad plane)
545 zpos = fgkCuZpos;
546 sprintf(cTagV,"UL%02d",iDet+kNdet);
547 sprintf(cTagM,"UG%02d",iDet+kNdet);
548 gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
549 // G10 layer (support structure)
550 zpos = fgkSuZpos;
551 sprintf(cTagV,"UM%02d",iDet+kNdet);
552 sprintf(cTagM,"UG%02d",iDet+kNdet);
553 gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
554 // Cu layer (FEE + signal lines)
555 zpos = fgkFeZpos;
556 sprintf(cTagV,"UN%02d",iDet+kNdet);
557 sprintf(cTagM,"UG%02d",iDet+kNdet);
558 gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
559 // Al layer (cooling devices)
560 zpos = fgkCoZpos;
561 sprintf(cTagV,"UO%02d",iDet+kNdet);
562 sprintf(cTagM,"UG%02d",iDet+kNdet);
563 gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
564 // Water layer (cooling)
565 zpos = fgkWaZpos;
566 sprintf(cTagV,"UP%02d",iDet+kNdet);
567 sprintf(cTagM,"UG%02d",iDet+kNdet);
568 gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
569 }
570 if (fClengthRH[iplan][icham] > 0.0) {
571 // Lower part
572 // Rohacell layer (radiator)
573 zpos = fgkRaZpos;
574 sprintf(cTagV,"UH%02d",iDet+2*kNdet);
575 sprintf(cTagM,"UC%02d",iDet+2*kNdet);
576 gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
577 // Mylar layer (entrance window + HV cathode)
578 zpos = fgkMyZpos;
579 sprintf(cTagV,"UI%02d",iDet+2*kNdet);
580 sprintf(cTagM,"UC%02d",iDet+2*kNdet);
581 gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
582 // Xe/Isobutane layer (drift volume)
583 zpos = fgkDrZpos;
584 sprintf(cTagV,"UJ%02d",iDet+2*kNdet);
585 sprintf(cTagM,"UC%02d",iDet+2*kNdet);
586 gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
587 // Upper part
588 // Xe/Isobutane layer (amplification volume)
589 zpos = fgkAmZpos;
590 sprintf(cTagV,"UK%02d",iDet+2*kNdet);
591 sprintf(cTagM,"UE%02d",iDet+2*kNdet);
592 gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
593 // Readout part
594 // Cu layer (pad plane)
595 zpos = fgkCuZpos;
596 sprintf(cTagV,"UL%02d",iDet+2*kNdet);
597 sprintf(cTagM,"UG%02d",iDet+2*kNdet);
598 gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
599 // G10 layer (support structure)
600 zpos = fgkSuZpos;
601 sprintf(cTagV,"UM%02d",iDet+2*kNdet);
602 sprintf(cTagM,"UG%02d",iDet+2*kNdet);
603 gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
604 // Cu layer (FEE + signal lines)
605 zpos = fgkFeZpos;
606 sprintf(cTagV,"UN%02d",iDet+2*kNdet);
607 sprintf(cTagM,"UG%02d",iDet+2*kNdet);
608 gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
609 // Al layer (cooling devices)
610 zpos = fgkCoZpos;
611 sprintf(cTagV,"UO%02d",iDet+2*kNdet);
612 sprintf(cTagM,"UG%02d",iDet+2*kNdet);
613 gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
614 // Water layer (cooling)
615 zpos = fgkWaZpos;
616 sprintf(cTagV,"UP%02d",iDet+2*kNdet);
617 sprintf(cTagM,"UG%02d",iDet+2*kNdet);
618 gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
619 }
620
621 // Position the inner volumes of the chambers in the frames
622 xpos = 0.0;
623 ypos = 0.0;
624 zpos = 0.0;
625 // The inside of the lower G10 frame
626 sprintf(cTagV,"UC%02d",iDet);
627 sprintf(cTagM,"UB%02d",iDet);
628 gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
629 // The lower G10 frame inside the aluminum frame
630 sprintf(cTagV,"UB%02d",iDet);
631 sprintf(cTagM,"UA%02d",iDet);
632 gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
633 // The inside of the upper G10 frame
634 sprintf(cTagV,"UE%02d",iDet);
635 sprintf(cTagM,"UD%02d",iDet);
636 gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
637 // The inside of the upper aluminum frame
638 sprintf(cTagV,"UG%02d",iDet);
639 sprintf(cTagM,"UF%02d",iDet);
640 gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
641 if (fClengthPH[iplan][icham] > 0.0) {
642 // The inside of the lower G10 frame
643 sprintf(cTagV,"UC%02d",iDet+kNdet);
644 sprintf(cTagM,"UB%02d",iDet+kNdet);
645 gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
646 // The lower G10 frame inside the aluminum frame
647 sprintf(cTagV,"UB%02d",iDet+kNdet);
648 sprintf(cTagM,"UA%02d",iDet+kNdet);
649 gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
650 // The inside of the upper G10 frame
651 sprintf(cTagV,"UE%02d",iDet+kNdet);
652 sprintf(cTagM,"UD%02d",iDet+kNdet);
653 gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
654 // The inside of the upper aluminum frame
655 sprintf(cTagV,"UG%02d",iDet+kNdet);
656 sprintf(cTagM,"UF%02d",iDet+kNdet);
657 gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
658 }
659 if (fClengthRH[iplan][icham] > 0.0) {
660 // The inside of the lower G10 frame
661 sprintf(cTagV,"UC%02d",iDet+2*kNdet);
662 sprintf(cTagM,"UB%02d",iDet+2*kNdet);
663 gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
664 // The lower G10 frame inside the aluminum frame
665 sprintf(cTagV,"UB%02d",iDet+2*kNdet);
666 sprintf(cTagM,"UA%02d",iDet+2*kNdet);
667 gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
668 // The inside of the upper G10 frame
669 sprintf(cTagV,"UE%02d",iDet+2*kNdet);
670 sprintf(cTagM,"UD%02d",iDet+2*kNdet);
671 gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
672 // The inside of the upper aluminum frame
673 sprintf(cTagV,"UG%02d",iDet+2*kNdet);
674 sprintf(cTagM,"UF%02d",iDet+2*kNdet);
675 gMC->Gspos(cTagV,1,cTagM,xpos,ypos,zpos,0,"ONLY");
676 }
677
678 // Position the frames of the chambers in the TRD mother volume
679 xpos = 0.;
680 ypos = - fClength[iplan][0] - fClength[iplan][1] - fClength[iplan][2]/2.;
681 for (Int_t ic = 0; ic < icham; ic++) {
682 ypos += fClength[iplan][ic];
683 }
684 ypos += fClength[iplan][icham]/2.;
685 zpos = fgkCraH/2. + fgkCdrH/2. - fgkSheight/2. + iplan * (fgkCH + fgkVspace);
686 // The lower aluminum frame, radiator + drift region
687 sprintf(cTagV,"UA%02d",iDet);
688 gMC->Gspos(cTagV,1,"UTR1",xpos,ypos,zpos,0,"ONLY");
689 // The upper G10 frame, amplification region
690 sprintf(cTagV,"UD%02d",iDet);
691 zpos += fgkCamH/2. + fgkCraH/2. + fgkCdrH/2.;
692 gMC->Gspos(cTagV,1,"UTR1",xpos,ypos,zpos,0,"ONLY");
693 // The upper aluminum frame
694 sprintf(cTagV,"UF%02d",iDet);
695 zpos += fgkCroH/2. + fgkCamH/2.;
696 gMC->Gspos(cTagV,1,"UTR1",xpos,ypos,zpos,0,"ONLY");
697 if (fClengthPH[iplan][icham] > 0.0) {
698 xpos = 0.;
699 ypos = - fgkSlenTR2/2.;
700 for (Int_t ic = 0; ic < icham; ic++) {
701 ypos += fClengthPH[iplan][ic];
702 }
703 ypos += fClengthPH[iplan][icham]/2.;
704 zpos = fgkCraH/2. + fgkCdrH/2. - fgkSheight/2. + iplan * (fgkCH + fgkVspace);
705 // The lower aluminum frame, radiator + drift region
706 sprintf(cTagV,"UA%02d",iDet+kNdet);
707 gMC->Gspos(cTagV,1,"UTR2",xpos,ypos,zpos,0,"ONLY");
708 // The upper G10 frame, amplification region
709 sprintf(cTagV,"UD%02d",iDet+kNdet);
710 zpos += fgkCamH/2. + fgkCraH/2. + fgkCdrH/2.;
711 gMC->Gspos(cTagV,1,"UTR2",xpos,ypos,zpos,0,"ONLY");
712 // The upper aluminum frame
713 sprintf(cTagV,"UF%02d",iDet+kNdet);
714 zpos += fgkCroH/2. + fgkCamH/2.;
715 gMC->Gspos(cTagV,1,"UTR2",xpos,ypos,zpos,0,"ONLY");
716 }
717 if (fClengthRH[iplan][icham] > 0.0) {
718 xpos = 0.;
719 ypos = - fgkSlenTR3/2.;
720 for (Int_t ic = 0; ic < icham; ic++) {
721 ypos += fClengthRH[iplan][ic];
722 }
723 ypos += fClengthRH[iplan][icham]/2.;
724 zpos = fgkCraH/2. + fgkCdrH/2. - fgkSheight/2. + iplan * (fgkCH + fgkVspace);
725 // The lower aluminum frame, radiator + drift region
726 sprintf(cTagV,"UA%02d",iDet+2*kNdet);
727 gMC->Gspos(cTagV,1,"UTR3",xpos,ypos,zpos,0,"ONLY");
728 // The upper G10 frame, amplification region
729 sprintf(cTagV,"UD%02d",iDet+2*kNdet);
730 zpos += fgkCamH/2. + fgkCraH/2. + fgkCdrH/2.;
731 gMC->Gspos(cTagV,1,"UTR3",xpos,ypos,zpos,0,"ONLY");
732 // The upper aluminum frame
733 sprintf(cTagV,"UF%02d",iDet+2*kNdet);
734 zpos += fgkCroH/2. + fgkCamH/2.;
735 gMC->Gspos(cTagV,1,"UTR3",xpos,ypos,zpos,0,"ONLY");
736 }
737
738 }
f7336fa3 739 }
740
0a770ac9 741 xpos = 0.;
742 ypos = 0.;
743 zpos = 0.;
744 gMC->Gspos("UTR1",1,"BTR1",xpos,ypos,zpos,0,"ONLY");
745 gMC->Gspos("UTR2",2,"BTR2",xpos,ypos,zpos,0,"ONLY");
746 gMC->Gspos("UTR3",3,"BTR3",xpos,ypos,zpos,0,"ONLY");
f7336fa3 747
748}
8737e16f 749
750//_____________________________________________________________________________
751void AliTRDgeometryHole::SetOldGeometry()
752{
753 //
754 // Use the old chamber lengths
755 //
756
757 Int_t icham;
758 Int_t iplan;
759
760 AliTRDgeometry::SetOldGeometry();
761
762 // The outer lengths of the chambers for the sectors with holes for the PHOS
763 Float_t lengthPH[kNplan][kNcham] = { { 0.0, 0.0, 0.0, 116.5, 123.5 }
764 , { 0.0, 0.0, 0.0, 124.0, 131.0 }
765 , { 0.0, 0.0, 0.0, 131.5, 134.5 }
766 , { 0.0, 0.0, 0.0, 139.0, 142.0 }
767 , { 0.0, 0.0, 0.0, 146.0, 142.0 }
768 , { 0.0, 0.0, 0.0, 153.5, 134.5 } };
769
770 // The outer lengths of the chambers for the sectors with holes for the RICH
771 Float_t lengthRH[kNplan][kNcham] = { { 0.0, 0.0, 0.0, 0.0, 86.5 }
772 , { 0.0, 0.0, 0.0, 0.0, 101.5 }
773 , { 0.0, 0.0, 0.0, 0.0, 112.5 }
774 , { 0.0, 0.0, 0.0, 0.0, 127.5 }
775 , { 0.0, 0.0, 0.0, 0.0, 134.5 }
776 , { 0.0, 0.0, 0.0, 0.0, 134.5 } };
777
778 for (icham = 0; icham < kNcham; icham++) {
779 for (iplan = 0; iplan < kNplan; iplan++) {
780 fClengthPH[iplan][icham] = lengthPH[iplan][icham];
781 fClengthRH[iplan][icham] = lengthRH[iplan][icham];
782 }
783 }
784
785}
786
787
788