1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
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 **************************************************************************/
18 //_________________________________________________________________________
19 // Geometry class for PHOS : EMCA (Electromagnetic Calorimeter)
20 // Its data members provide geometry parametrization of EMCA
21 // which can be changed in the constructor only.
22 // Author : Yves Schutz (SUBATECH)
23 // Modified : Yuri Kharlov (IHEP, Protvino)
25 // Modified : Dmitri Peressounko (RRC "Kurchatov Institute")
28 // --- AliRoot header files ---
30 #include "AliPHOSEMCAGeometry.h"
32 ClassImp(AliPHOSEMCAGeometry)
34 //____________________________________________________________________________
35 AliPHOSEMCAGeometry::AliPHOSEMCAGeometry():
37 fStripWallWidthOut(0.f),
38 fStripWallWidthIn(0.f),
40 fInnerThermoWidthX(0.f),
41 fInnerThermoWidthY(0.f),
42 fInnerThermoWidthZ(0.f),
49 fAlCoverThickness(0.f),
50 fOuterThermoWidthXUp(0.f),
51 fOuterThermoWidthXLow(0.f),
52 fOuterThermoWidthY(0.f),
53 fOuterThermoWidthZ(0.f),
56 fFiberGlassSup2X(0.f),
57 fFiberGlassSup1X(0.f),
63 fWarmUpperThickness(0.f),
64 fWarmBottomThickness(0.f),
65 fWarmAlCoverWidthX(0.f),
66 fWarmAlCoverWidthY(0.f),
67 fWarmAlCoverWidthZ(0.f),
68 fFiberGlassSup1Y(0.f),
69 fFiberGlassSup2Y(0.f),
71 fTSupport1Thickness(0.f),
72 fTSupport2Thickness(0.f),
75 fIPtoOuterCoverDistance(0.f),
76 fIPtoCrystalSurface(0.f),
77 fSupportPlateThickness(0.f),
88 // Initializes the EMC parameters
89 // Coordinate system chosen: x across beam, z along beam, y out of beam.
90 // Reference point for all volumes incide module is
91 // center of module in x,z on the upper surface of support beam
93 //Distance from IP to surface of the crystals
94 fIPtoCrystalSurface = 460.0 ;
99 fCrystalHalfSize[0] = 2.2 /2 ; //Half-Sizes of crystall
100 fCrystalHalfSize[1] = 18.0 /2 ;
101 fCrystalHalfSize[2] = 2.2 /2 ;
105 //fPinDiodeSize[0] = 1.71 ; //Values of ame PIN diode
106 //fPinDiodeSize[1] = 0.0280 ; // OHO 0.0280 is the depth of active layer
107 //fPinDiodeSize[2] = 1.61 ;
109 fPinDiodeHalfSize[0] = 0.5000 /2 ; // APD 5 mm side
110 fPinDiodeHalfSize[1] = 0.0100 /2 ; // APD bulk thickness
111 fPinDiodeHalfSize[2] = 0.5000 /2 ; // APD 5 mm side
113 fPreampHalfSize[0] = 1.5 / 2 ; // Preamplifier
114 fPreampHalfSize[1] = 0.5 / 2 ;
115 fPreampHalfSize[2] = 1.5 / 2 ;
117 //Strip unit (8x2 crystals)
119 fNCellsXInStrip = 8 ; //Number of crystals in strip unit along x-axis
120 fNCellsZInStrip = 2 ; //Number of crystals in strip unit along z-axis
121 fNStripX = 8 ; //Number of strip units across along x-axis
122 fNStripZ = 28 ; //Number of strips along z-axis
124 fStripWallWidthOut = 0.01 ; // Side to another strip
125 fStripWallWidthIn = 0.02 ; // Side betveen crystals in one strip
127 fTyvecThickness = 0.0175 ; //Thickness of the tyvec
129 fAirGapLed = 1.5 - 2 * fPreampHalfSize[1] - 2 * fPinDiodeHalfSize[1] ; // Air gap before crystalls for LED system
130 // Note, that Cell in Strip 1.5 longer then crystall
132 //---Now calculate thechnical sizes for GEANT implementation
134 fWrappedHalfSize[0] = (2*fTyvecThickness + 2*fCrystalHalfSize[0])/2 ; //This will be size of crystall
135 fWrappedHalfSize[1] = fCrystalHalfSize[1] ; //wrapped into tyvec
136 fWrappedHalfSize[2] = (2*fTyvecThickness + 2*fCrystalHalfSize[2])/2 ; //
138 fAirCellHalfSize[0] = fWrappedHalfSize[0] + 0.01;
139 fAirCellHalfSize[1] = (fAirGapLed + 2*fPreampHalfSize[1] +
140 2*fPinDiodeHalfSize[1] + 2*fWrappedHalfSize[1])/2 ; //in strip
141 fAirCellHalfSize[2] = fWrappedHalfSize[2] + 0.01;
143 // fSupportPlateHalfSize[0] = ( (fNCellsXInStrip-1)*fStripWallWidthIn + 2*fStripWallWidthOut +
144 // fNCellsXInStrip * (2*fTyvecThickness + 2*fCrystalHalfSize[0]) )/2 ;
145 fSupportPlateHalfSize[0] = 18.04 /2 ;
146 fSupportPlateHalfSize[1] = 6.0 /2 ;
147 // fSupportPlateHalfSize[2] = ( (fNCellsZInStrip-1)*fStripWallWidthIn + 2*fStripWallWidthOut +
148 // fNCellsZInStrip * (2*fTyvecThickness + 2*fCrystalHalfSize[2]) )/2;
149 fSupportPlateHalfSize[2] = 4.51 /2 ;
150 fSupportPlateThickness = 0.3 ;
151 fSupportPlateInHalfSize[0] = fSupportPlateHalfSize[0] ; //Half-sizes of the air
152 fSupportPlateInHalfSize[1] = fSupportPlateHalfSize[1]-fSupportPlateThickness ; //box in the support plate
153 fSupportPlateInHalfSize[2] = fSupportPlateHalfSize[2]-fSupportPlateThickness/2 ;
155 fStripHalfSize[0]= fSupportPlateHalfSize[0] ;
156 fStripHalfSize[1]= ( 2*fSupportPlateHalfSize[1] + 2*fAirCellHalfSize[1] )/2;
157 fStripHalfSize[2]= fSupportPlateHalfSize[2] ;
159 // ------- Inner hermoinsulation ---------------
160 fInnerThermoWidthX = 2.0 ; // Width of the innerthermoinsulation across the beam
161 fInnerThermoWidthY = 2.0 ; // Width of the upper cover of innerthermoinsulation
162 fInnerThermoWidthZ = 2.0 ; // Width of the innerthermoinsulation along the beam
164 fInnerThermoHalfSize[0] = (2 * fStripHalfSize[0] * fNStripX + 2 * fInnerThermoWidthX ) /2 ;
165 fInnerThermoHalfSize[1] = (2 * fStripHalfSize[1] + fInnerThermoWidthY ) /2 ;
166 fInnerThermoHalfSize[2] = (2 * fStripHalfSize[2] * fNStripZ + 2 * fInnerThermoWidthZ ) /2 ;
168 // ------- Air gap between inner thermoinsulation and passive coller ---------
170 fAirGapWidthX = 0.2 ; // Width of the air gap across the beam
171 fAirGapWidthY = 0.2 ; // Width of the upper air gap
172 fAirGapWidthZ = 0.2 ; // Width of the air gap along the beam
174 fAirGapHalfSize[0] = (2 * fInnerThermoHalfSize[0] + 2 * fAirGapWidthX ) /2 ;
175 fAirGapHalfSize[1] = (2 * fInnerThermoHalfSize[1] + fAirGapWidthY ) /2 ;
176 fAirGapHalfSize[2] = (2 * fInnerThermoHalfSize[2] + 2 * fAirGapWidthZ ) /2 ;
178 // ------- Passive Cooler ------------------------
180 fCoolerWidthX = 2.0 ; // Width of the passive coller across the beam
181 fCoolerWidthY = 0.3 ; // Width of the upper cover of cooler
182 fCoolerWidthZ = 2.0 ; // Width of the passive cooler along the beam
184 fCoolerHalfSize[0] = (2 * fAirGapHalfSize[0] + 2 * fCoolerWidthX ) /2 ;
185 fCoolerHalfSize[1] = (2 * fAirGapHalfSize[1] + fCoolerWidthY ) /2 ;
186 fCoolerHalfSize[2] = (2 * fAirGapHalfSize[2] + 2 * fCoolerWidthZ ) /2 ;
188 // ------- Outer thermoinsulation and Al cover -------------------------------
190 fAlCoverThickness = 0.1 ; //Thickness of the Al cover of the module
192 fOuterThermoWidthXUp = 156.0 - fAlCoverThickness ;
193 //width of the upper surface of the PHOS module accross the beam
194 fOuterThermoWidthY = 6.0 ; // with of the upper cover of outer thermoinsulation
195 fOuterThermoWidthZ = 6.0 ; //width of the thermoinsulation along the beam
197 fAlFrontCoverX = 6.0 ; //Width of Al strip around fiberglass window: across
198 fAlFrontCoverZ = 6.0 ; //and along the beam
201 // Calculate distance from IP to upper cover
202 fIPtoOuterCoverDistance = fIPtoCrystalSurface - fAirGapLed - fInnerThermoWidthY - fAirGapWidthY -
203 fCoolerWidthY - fOuterThermoWidthY - fAlCoverThickness ;
205 Float_t tanA = fOuterThermoWidthXUp / (2.*fIPtoOuterCoverDistance) ;
206 // tan(a) where A = angle between IP to center and IP to side across beam
208 fOuterThermoWidthXLow = fOuterThermoWidthXUp +
209 2 * (2* fCoolerHalfSize[1] + fOuterThermoWidthY) * tanA
210 - fAlCoverThickness ;
211 //width of the lower surface of the COOL section accross the beam
214 fOuterThermoParams[0] = fOuterThermoWidthXUp / 2 ; // half-length along x at the z surface positioned at -DZ;
215 fOuterThermoParams[1] = fOuterThermoWidthXLow/ 2 ; // half-length along x at the z surface positioned at +DZ;
216 fOuterThermoParams[2] = ( 2 * fCoolerHalfSize[2] + 2 * fOuterThermoWidthZ ) / 2 ;
217 // `half-length along the y-axis' in out case this is z axis
218 fOuterThermoParams[3] = ( 2* fCoolerHalfSize[1] + fOuterThermoWidthY) /2 ;
219 // `half-length along the z-axis' in our case this is y axis
221 fAlCoverParams[0] = fOuterThermoParams[0] + fAlCoverThickness ;
222 fAlCoverParams[1] = fOuterThermoParams[1] + fAlCoverThickness ;
223 fAlCoverParams[2] = fOuterThermoParams[2] + fAlCoverThickness ;
224 fAlCoverParams[3] = fOuterThermoParams[3] + fAlCoverThickness /2 ;
227 fFiberGlassHalfSize[0] = fAlCoverParams[0] - fAlFrontCoverX ;
228 fFiberGlassHalfSize[1] = fAlCoverParams[2] - fAlFrontCoverZ ; //Note, here other ref. system
229 fFiberGlassHalfSize[2] = fAlCoverThickness / 2 ;
232 //============Now warm section======================
234 fWarmAlCoverWidthX = 2 * fAlCoverParams[1] ; //Across beam
235 fWarmAlCoverWidthY = 159.0 ; //along beam
238 fTSupport1Thickness = 3.5 ;
239 fTSupport2Thickness = 5.0 ;
240 fTSupport1Width = 10.6 ;
241 fTSupport2Width = 3.1 ;
242 fNTSupports = fNStripX + 1 ;
243 fTSupportDist = 7.48 ;
246 fAirSpaceFeeX = 148.6 ; //Across beam
247 fAirSpaceFeeY = 135.0 ; //along beam
248 fAirSpaceFeeZ = 19.0 ; //out of beam
251 fWarmBottomThickness = 4.0 ;
252 fWarmUpperThickness = 4.0 ;
255 fFrameThickness = 5.0 ;
256 fFrameHeight = 15.0 ;
259 fFiberGlassSup1X = 6.0 ;
260 fFiberGlassSup1Y = 3.9 + fWarmUpperThickness ;
262 fFiberGlassSup2X = 3.0 ;
263 fFiberGlassSup2Y = fFrameHeight ;
265 //Now calculate Half-sizes
267 fWarmAlCoverWidthZ = fAirSpaceFeeZ + fWarmBottomThickness + fWarmUpperThickness +
268 fTSupport1Thickness + fTSupport2Thickness ;
271 fWarmAlCoverHalfSize[0] = fWarmAlCoverWidthX / 2 ;
272 fWarmAlCoverHalfSize[1] = fWarmAlCoverWidthY / 2 ;
273 fWarmAlCoverHalfSize[2] = fWarmAlCoverWidthZ / 2 ;
276 fWarmThermoHalfSize[0] = fWarmAlCoverHalfSize[0] - fAlCoverThickness ;
277 fWarmThermoHalfSize[1] = fWarmAlCoverHalfSize[1] - fAlCoverThickness ;
278 fWarmThermoHalfSize[2] = fWarmAlCoverHalfSize[2] - fAlCoverThickness /2 ;
282 fTSupport1HalfSize[0] = fTSupport1Width /2 ; //Across beam
283 fTSupport1HalfSize[1] = (fAirSpaceFeeY + 2*fFiberGlassSup1X) /2 ; //along beam
284 fTSupport1HalfSize[2] = fTSupport1Thickness /2; //out of beam
286 fTSupport2HalfSize[0] = fTSupport2Width /2; //Across beam
287 fTSupport2HalfSize[1] = fTSupport1HalfSize[1] ; //along beam
288 fTSupport2HalfSize[2] = fTSupport2Thickness /2; //out of beam
291 fTCables1HalfSize[0] = (2*fTSupport1HalfSize[0]*fNTSupports + (fNTSupports-1)* fTSupportDist) / 2 ; //Across beam
292 fTCables1HalfSize[1] = fTSupport1HalfSize[1] ; //along beam
293 fTCables1HalfSize[2] = fTSupport1HalfSize[2] ; //out of beam
295 fTCables2HalfSize[0] = fTCables1HalfSize[0] ; //Across beam
296 fTCables2HalfSize[1] = fTSupport2HalfSize[1] ; //along beam
297 fTCables2HalfSize[2] = fTSupport2HalfSize[2] ; //out of beam
299 //frame: we define two frames along beam ...Z and across beam ...X
300 fFrameXHalfSize[0] = (fAirSpaceFeeX + 2 * fFiberGlassSup2X + 2* fFrameThickness) /2 ;
301 fFrameXHalfSize[1] = fFrameThickness /2 ;
302 fFrameXHalfSize[2] = fFrameHeight /2 ;
304 fFrameXPosition[0] = 0 ;
305 fFrameXPosition[1] = fAirSpaceFeeY /2 + fFiberGlassSup2X + fFrameXHalfSize[1] ;
306 fFrameXPosition[2] = fWarmThermoHalfSize[2] - fFrameHeight/ 2 - fWarmBottomThickness ;
308 fFrameZHalfSize[0] = fFrameThickness /2 ;
309 fFrameZHalfSize[1] = (fAirSpaceFeeY + 2 * fFiberGlassSup2X) /2 ;
310 fFrameZHalfSize[2] = fFrameHeight /2 ;
312 fFrameZPosition[0] = fAirSpaceFeeX /2 + fFiberGlassSup2X + fFrameZHalfSize[0] ;
313 fFrameZPosition[1] = 0 ;
314 fFrameZPosition[2] = fWarmThermoHalfSize[2] - fFrameHeight/ 2 - fWarmBottomThickness ;
316 //Fiberglass support define 4 fiber glass supports 2 along Z and 2 along X
318 fFGupXHalfSize[0] = fFrameXHalfSize[0] ;
319 fFGupXHalfSize[1] = fFiberGlassSup1X /2 ;
320 fFGupXHalfSize[2] = fFiberGlassSup1Y /2;
322 fFGupXPosition[0] = 0 ;
323 fFGupXPosition[1] = fAirSpaceFeeY /2 + fFGupXHalfSize[1] ;
324 fFGupXPosition[2] = fWarmThermoHalfSize[2] - fFrameHeight - fWarmBottomThickness - fFGupXHalfSize[2] ;
326 fFGupZHalfSize[0] = fFiberGlassSup1X /2 ;
327 fFGupZHalfSize[1] = fAirSpaceFeeY /2 ;
328 fFGupZHalfSize[2] = fFiberGlassSup1Y /2;
330 fFGupZPosition[0] = fAirSpaceFeeX /2 + fFGupZHalfSize[0] ;
331 fFGupZPosition[1] = 0 ;
332 fFGupZPosition[2] = fWarmThermoHalfSize[2] - fFrameHeight - fWarmBottomThickness - fFGupXHalfSize[2] ;
334 fFGlowXHalfSize[0] = fFrameXHalfSize[0] - 2*fFrameZHalfSize[0] ;
335 fFGlowXHalfSize[1] = fFiberGlassSup2X /2 ;
336 fFGlowXHalfSize[2] = fFrameXHalfSize[2] ;
338 fFGlowXPosition[0] = 0 ;
339 fFGlowXPosition[1] = fAirSpaceFeeY /2 + fFGlowXHalfSize[1] ;
340 fFGlowXPosition[2] = fWarmThermoHalfSize[2] - fWarmBottomThickness - fFGlowXHalfSize[2] ;
342 fFGlowZHalfSize[0] = fFiberGlassSup2X /2 ;
343 fFGlowZHalfSize[1] = fAirSpaceFeeY /2 ;
344 fFGlowZHalfSize[2] = fFrameZHalfSize[2] ;
346 fFGlowZPosition[0] = fAirSpaceFeeX /2 + fFGlowZHalfSize[0] ;
347 fFGlowZPosition[1] = 0 ;
348 fFGlowZPosition[2] = fWarmThermoHalfSize[2] - fWarmBottomThickness - fFGlowXHalfSize[2] ;
351 // --- Air Gap for FEE ----
353 fFEEAirHalfSize[0] = fAirSpaceFeeX /2 ;
354 fFEEAirHalfSize[1] = fAirSpaceFeeY /2;
355 fFEEAirHalfSize[2] = fAirSpaceFeeZ /2;
357 fFEEAirPosition[0] = 0 ;
358 fFEEAirPosition[1] = 0 ;
359 fFEEAirPosition[2] = fWarmThermoHalfSize[2] - fWarmBottomThickness - fFEEAirHalfSize[2] ;
361 // --- Calculate the oveol dimentions of the EMC module
363 fEMCParams[3] = fAlCoverParams[3] + fWarmAlCoverHalfSize[2] ; //Size out of beam
364 fEMCParams[0] = fAlCoverParams[0] ; //Upper size across the beam
365 fEMCParams[1] = (fAlCoverParams[1] - fAlCoverParams[0])*fEMCParams[3]/fAlCoverParams[3]
366 + fAlCoverParams[0] ; //Lower size across the beam
367 fEMCParams[2] = fWarmAlCoverHalfSize[1] ; // Size along the beam
369 fNPhi = fNStripX * fNCellsXInStrip ; //number of crystals across the beam
370 fNZ = fNStripZ * fNCellsZInStrip ; //number of crystals along the beam