]> git.uio.no Git - u/mrichter/AliRoot.git/blob - STRUCT/AliABSOv0.cxx
MAG geometry with resolved MANY (I.Hrivnacova)
[u/mrichter/AliRoot.git] / STRUCT / AliABSOv0.cxx
1 /**************************************************************************
2  * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
3  *                                                                        *
4  * Author: The ALICE Off-line Project.                                    *
5  * Contributors are mentioned in the code where appropriate.              *
6  *                                                                        *
7  * Permission to use, copy, modify and distribute this software and its   *
8  * documentation strictly for non-commercial purposes is hereby granted   *
9  * without fee, provided that the above copyright notice appears in all   *
10  * copies and that both the copyright notice and this permission notice   *
11  * appear in the supporting documentation. The authors make no claims     *
12  * about the suitability of this software for any purpose. It is          *
13  * provided "as is" without express or implied warranty.                  *
14  **************************************************************************/
15
16 /*
17 $Log$
18 Revision 1.11.6.1  2002/07/24 10:08:27  alibrary
19 Updating VirtualMC
20
21 Revision 1.12  2002/07/15 07:56:41  morsch
22 Steel rear and  Al ring.
23
24 Revision 1.11  2001/11/29 14:16:51  morsch
25 - truncated inner W-cone
26 - new support structure
27
28 Revision 1.10  2001/10/25 08:50:57  morsch
29 New beamshield geometry with increased tolerances and insulation thickness.
30
31 Revision 1.9  2001/05/16 14:57:22  alibrary
32 New files for folders and Stack
33
34 Revision 1.8  2001/01/12 13:16:09  morsch
35 Store absorber composition information in fMLayers and fZLayers
36 Rear 25 cm Fe + 35 cm Cu
37
38 Revision 1.7  2000/10/02 21:28:15  fca
39 Removal of useless dependecies via forward declarations
40
41 Revision 1.6  2000/06/15 09:40:31  morsch
42 Obsolete typedef keyword removed
43
44 Revision 1.5  2000/06/12 19:39:01  morsch
45 New structure of beam pipe and heating jacket.
46
47 Revision 1.4  2000/04/03 08:13:40  fca
48 Introduce extra scope for non ANSI compliant C++ compilers
49
50 Revision 1.3  2000/01/18 17:49:56  morsch
51 Serious overlap of ABSM with shield corrected
52 Small error in ARPB parameters corrected
53
54 Revision 1.2  2000/01/13 11:23:59  morsch
55 Last layer of Pb outer angle corrected
56
57 Revision 1.1  2000/01/12 15:39:30  morsch
58 Standard version of ABSO
59
60 */
61
62 ///////////////////////////////////////////////////////////////////////////////
63 //                                                                           //
64 //  Muon ABSOrber                                                            //
65 //  This class contains the description of the muon absorber geometry        //
66 //                                                                           //
67 //Begin_Html
68 /*
69 <img src="picts/AliABSOClass.gif">
70 </pre>
71 <br clear=left>
72 <font size=+2 color=red>
73 <p>The responsible person for this module is
74 <a href="mailto:andreas.morsch@cern.ch">Andreas Morsch</a>.
75 </font>
76 <pre>
77 */
78 //End_Html
79 //                                                                           //
80 //                                                                           //
81 ///////////////////////////////////////////////////////////////////////////////
82
83 #include "AliABSOv0.h"
84 #include "AliRun.h"
85 #include "AliMC.h"
86 #include "AliConst.h"
87
88 ClassImp(AliABSOv0)
89  
90 //_____________________________________________________________________________
91 AliABSOv0::AliABSOv0()
92 {
93   //
94   // Default constructor
95   //
96 }
97  
98 //_____________________________________________________________________________
99 AliABSOv0::AliABSOv0(const char *name, const char *title)
100        : AliABSO(name,title)
101 {
102   //
103   // Standard constructor
104   //
105   SetMarkerColor(7);
106   SetMarkerStyle(2);
107   SetMarkerSize(0.4);
108 }
109  
110 //_____________________________________________________________________________
111 void AliABSOv0::CreateGeometry()
112 {
113     //
114     // Creation of the geometry of the muon absorber
115     //
116     //Begin_Html
117     /*
118       <img src="picts/AliABSOv0Tree.gif">
119     */
120     //End_Html
121     //Begin_Html
122     /*
123       <img src="picts/AliABSOv0.gif">
124     */
125     //End_Html
126     
127     //
128     //
129
130     enum {kC=1605, kAl=1608, kFe=1609, kCu=1610, kW=1611, kPb=1612,
131           kNiCuW=1620, kVacuum=1615, kAir=1614, kConcrete=1616,
132           kPolyCH2=1617, kSteel=1609, kInsulation=1613, kPolyCc=1619};    
133     
134     Int_t *idtmed = fIdtmed->GetArray()-1599;
135     
136     Float_t par[24], cpar[5], cpar0[5], pcpar[12], tpar[3], tpar0[3]; 
137     Float_t dz;
138     Int_t idrotm[1699];
139 #include "ABSOSHILConst.h"
140 #include "ABSOConst.h"
141 //
142 // Structure of Tracking Region
143 //
144   Float_t dzFe = 25.;
145
146 // 3 < theta < 9
147     fNLayers[0] = 5; 
148     fMLayers[0][0]  = kAir;              fZLayers[0][0] = zAbsStart;
149     fMLayers[0][1]  = kC;                fZLayers[0][1] = zAbsCc;             
150     fMLayers[0][2]  = kConcrete;         fZLayers[0][2] = zRear-dRear-dzFe;
151     fMLayers[0][3]  = kSteel;            fZLayers[0][3] = zRear-dRear;
152     fMLayers[0][4]  = kSteel;            fZLayers[0][4] = zRear;
153 // 2 < theta < 3
154     fNLayers[1] = 5; 
155     fMLayers[1][0] = fMLayers[0][0];      fZLayers[1][0] = fZLayers[0][0];
156     fMLayers[1][1] = fMLayers[0][1];      fZLayers[1][1] = fZLayers[0][1];
157     fMLayers[1][2] = fMLayers[0][2];      fZLayers[1][2] = fZLayers[0][2];
158     fMLayers[1][3] = fMLayers[0][3];      fZLayers[1][3] = fZLayers[0][3];
159     fMLayers[1][4] = kNiCuW;              fZLayers[1][4] = fZLayers[0][4];
160 //    
161
162     Float_t dTube=0.1;                     // tube thickness
163     Float_t dInsu=0.5;                     // insulation thickness
164     Float_t dEnve=0.1;                     // protective envelope thickness
165     Float_t dFree=0.5;                     // clearance thickness
166
167
168 // Mother volume and outer shielding: Pb
169   par[0]  = 0.;
170   par[1]  = 360.;
171   par[2]  = 7.;
172     
173   par[3]  = -(zRear-zAbsStart)/2.;
174   par[4]  = rAbs;
175   par[5]  = zAbsStart * TMath::Tan(theta1);
176
177   par[6]  = par[3]+(zNose-zAbsStart);
178   par[7]  = rAbs;
179   par[8]  = zNose * TMath::Tan(theta1);
180
181   par[9]  = par[3]+(zConeTPC-zAbsStart);
182   par[10] = rAbs;
183   par[11] = par[8] + (par[9] - par[6]) * TMath::Tan(theta2);
184
185   par[12]  = par[3]+(zOpen-zAbsStart);
186   par[13] = rAbs;
187   par[14] = par[11] + (par[12] - par[9]) * TMath::Tan(accMax);
188
189   par[15] = par[3]+(zRear-dRear-zAbsStart);
190   par[16] = rAbs   + (par[15] - par[12]) * TMath::Tan(thetaOpen1) ;
191   par[17] = par[14] + (par[15] - par[12]) * TMath::Tan(accMax);
192
193   par[18] = par[3]+(zRear-dRear-zAbsStart);
194   par[19] = (zRear-dRear) * TMath::Tan(accMin);
195   par[20] = par[14] + (par[18] - par[12]) * TMath::Tan(accMax);
196
197   par[21] = -par[3];
198   par[22] =  zRear* TMath::Tan(accMin);
199   par[23] = par[20] + (par[21] - par[18]) * TMath::Tan(accMax);
200   gMC->Gsvolu("ABSS", "PCON", idtmed[kPb], par, 24);
201   { // Begin local scope for i
202       for (Int_t i=4; i<18; i+=3) par[i]  = 0;
203   } // End local scope for i
204   gMC->Gsvolu("ABSM", "PCON", idtmed[kVacuum+40], par, 24);
205   gMC->Gspos("ABSS", 1, "ABSM", 0., 0., 0., 0, "ONLY");
206
207 //
208 // Steel envelope
209 //
210   par[4] = par[5] -dSteel;
211   par[7] = par[8] -dSteel;
212   par[10]= par[11]-dSteel;  
213   par[13]= par[14]-dSteel;  
214   par[16]= par[17]-dSteel;  
215   par[19]= par[20]-dSteel;  
216   par[22]= par[23]-dSteel;  
217   gMC->Gsvolu("ABST", "PCON", idtmed[kSteel], par, 24);
218   gMC->Gspos("ABST", 1, "ABSS", 0., 0., 0., 0, "ONLY");
219 //
220 // Polyethylene shield
221 // 
222   cpar[0] = (zRear - zConeTPC) / 2.;
223   cpar[1] = zConeTPC * TMath::Tan(accMax);
224   cpar[2] = cpar[1] + dPoly;
225   cpar[3] = zRear * TMath::Tan(accMax);
226   cpar[4] = cpar[3] + dPoly;
227   gMC->Gsvolu("APOL", "CONE", idtmed[kPolyCH2+40], cpar, 5);
228   dz = (zRear-zAbsStart)/2.-cpar[0];
229   gMC->Gspos("APOL", 1, "ABSS", 0., 0., dz, 0, "ONLY");
230
231 //
232 // Tungsten nose to protect TPC
233 // 
234   cpar[0] = (zNose - zAbsStart) / 2.;
235   cpar[1] = zAbsStart * TMath::Tan(accMax);
236   cpar[2] = zAbsStart * TMath::Tan(theta1)-dSteel;
237   cpar[3] = zNose * TMath::Tan(accMax);
238   cpar[4] = zNose * TMath::Tan(theta1)-dSteel;
239   gMC->Gsvolu("ANOS", "CONE", idtmed[kW], cpar, 5);
240   //
241   dz = -(zRear-zAbsStart)/2.+cpar[0];
242   gMC->Gspos("ANOS", 1, "ABSS", 0., 0., dz, 0, "ONLY");
243   //
244   // Tungsten inner shield
245   //
246   Float_t zW = zTwoDeg+.1;
247   Float_t dZ = zW+(zRear-dRear-zW)/2.;
248   //
249   pcpar[0]  = 0.;
250   pcpar[1]  = 360.;
251   pcpar[2]  = 3.;
252   pcpar[3]  = zW-dZ;
253   pcpar[4]  = rAbs;
254   pcpar[5]  = zW * TMath::Tan(accMin);
255   pcpar[6]  = zOpen-dZ;
256   pcpar[7]  = rAbs;
257   pcpar[8]  = zOpen * TMath::Tan(accMin);
258   pcpar[9]  = zRear-dRear-dZ;
259   pcpar[10] = rAbs+(zRear-dRear-zOpen) * TMath::Tan(thetaOpen1);
260   pcpar[11] = (zRear-dRear) * TMath::Tan(accMin);
261   
262   gMC->Gsvolu("AWIN", "PCON", idtmed[kNiCuW+40], pcpar, 12);
263   dz=(zW+zRear-dRear)/2-(zAbsStart+zRear)/2.;
264   gMC->Gspos("AWIN", 1, "ABSS", 0., 0., dz, 0, "ONLY");
265 //
266 // First part replaced by Carbon  
267 //
268   cpar[0] = (200.-zW)/2.;
269   cpar[1] = rAbs;
270   cpar[2] = pcpar[5];
271   cpar[3] = rAbs;
272   cpar[4] = 200. * TMath::Tan(accMin);
273   gMC->Gsvolu("ACNO", "CONE", idtmed[kC], cpar, 5);
274   dz = zW-dZ+cpar[0];
275   gMC->Gspos("ACNO", 1, "AWIN", 0., 0., dz, 0, "ONLY");
276
277   Float_t zWW = 383.5;
278 /*  
279   cpar[0] = (zRear-dRear-zWW)/2.;
280   cpar[1] = rAbs + (zWW-zOpen) *  TMath::Tan(thetaOpen1);
281   cpar[2] =  zWW * TMath::Tan(accMin);
282   cpar[3] = pcpar[10];
283   cpar[4] = pcpar[11];
284   gMC->Gsvolu("AWNO", "CONE", idtmed[kCu+40], cpar, 5);
285   dz = zWW-dZ+cpar[0];
286   
287   gMC->Gspos("AWNO", 1, "AWIN", 0., 0., dz, 0, "ONLY");
288 */
289   //
290   //     Inner tracking region
291   //
292   //     mother volume: Cu
293   //
294   //
295   pcpar[0]  = 0.;
296   pcpar[1]  = 360.;
297   pcpar[2]  = 3.;
298   pcpar[3]  = -(zRear-zAbsStart)/2.;
299   pcpar[4]  = rAbs;
300   pcpar[5]  = zAbsStart * TMath::Tan(accMax);
301   pcpar[6]  = pcpar[3]+(zTwoDeg-zAbsStart);
302   pcpar[7]  = rAbs;
303   pcpar[8]  = zTwoDeg * TMath::Tan(accMax);
304   pcpar[9]  = -pcpar[3];
305   pcpar[10] = zRear * TMath::Tan(accMin);
306   pcpar[11] = zRear * TMath::Tan(accMax);
307   gMC->Gsvolu("AITR", "PCON", idtmed[fMLayers[0][4]], pcpar, 12);
308   //
309   // special Pb medium for last 5 cm of Pb
310   Float_t zr=zRear-2.-0.001;
311   cpar[0] = 1.0;
312   cpar[1] = zr * TMath::Tan(thetaR);
313   cpar[2] = zr * TMath::Tan(accMax);
314   cpar[3] = cpar[1] + TMath::Tan(thetaR) * 2;
315   cpar[4] = cpar[2] + TMath::Tan(accMax) * 2;
316   gMC->Gsvolu("ARPB", "CONE", idtmed[fMLayers[0][4]], cpar, 5);
317   dz=(zRear-zAbsStart)/2.-cpar[0]-0.001;
318   gMC->Gspos("ARPB", 1, "AITR", 0., 0., dz, 0, "ONLY");
319   //
320   //     concrete cone: concrete 
321   //
322   pcpar[9]  = pcpar[3]+(zRear-dRear-zAbsStart);
323   pcpar[10] = (zRear-dRear) * TMath::Tan(accMin);
324   pcpar[11] = (zRear-dRear) * TMath::Tan(accMax);
325   gMC->Gsvolu("ACON", "PCON", idtmed[fMLayers[0][2]+40], pcpar, 12);
326   gMC->Gspos("ACON", 1, "AITR", 0., 0., 0., 0, "ONLY");
327 //
328 //    Fe Cone 
329 //
330   zr = zRear-dRear-dzFe;
331   cpar[0]  = dzFe/2.;
332   cpar[1] = zr * TMath::Tan(accMin);
333   cpar[2] = zr * TMath::Tan(accMax);
334   cpar[3] = cpar[1] + TMath::Tan(accMin) * dzFe;
335   cpar[4] = cpar[2] + TMath::Tan(accMax) * dzFe;
336   gMC->Gsvolu("ACFE", "CONE",idtmed[fMLayers[0][3]], cpar, 5);
337
338   dz = (zRear-zAbsStart)/2.-dRear-dzFe/2.;
339
340   gMC->Gspos("ACFE", 1, "ACON", 0., 0., dz, 0, "ONLY");
341
342   
343   //
344   //
345   //     carbon cone: carbon
346   //
347   pcpar[9]  = pcpar[3]+(zAbsCc-zAbsStart);
348   pcpar[10]  = zAbsCc * TMath::Tan(accMin);
349   pcpar[11]  = zAbsCc * TMath::Tan(accMax);
350   gMC->Gsvolu("ACAR", "PCON", idtmed[fMLayers[0][1]+40], pcpar, 12);
351   gMC->Gspos("ACAR", 1, "ACON", 0., 0., 0., 0, "ONLY");
352  //
353  //     carbon cone outer region
354  //
355   cpar[0]  = 10.;
356   cpar[1]  = rAbs;
357   cpar[2]  = zAbsStart* TMath::Tan(accMax);
358   cpar[3]  = rAbs;
359   cpar[4]  = cpar[2]+2. * cpar[0] * TMath::Tan(accMax);
360
361   gMC->Gsvolu("ACAO", "CONE", idtmed[fMLayers[0][1]], cpar, 5);
362   dz=-(zRear-zAbsStart)/2.+cpar[0];
363   gMC->Gspos("ACAO", 1, "ACAR", 0., 0., dz, 0, "ONLY");
364   //
365   //     inner W shield
366   Float_t epsi=0.;
367   Float_t repsi=1.;
368   
369   zr=zRear-(dRear-epsi);
370   cpar[0] = (dRear-epsi)/2.;
371   cpar[1] = zr * TMath::Tan(accMin);
372   cpar[2] = zr * TMath::Tan(thetaR*repsi);
373   cpar[3] = cpar[1] + TMath::Tan(accMin) * (dRear-epsi);
374   cpar[4] = cpar[2] + TMath::Tan(thetaR*repsi) * (dRear-epsi);
375   gMC->Gsvolu("ARW0", "CONE", idtmed[fMLayers[1][4]+40], cpar, 5);
376   dz=(zRear-zAbsStart)/2.-cpar[0];
377   gMC->Gspos("ARW0", 1, "AITR", 0., 0., dz, 0, "ONLY");
378   //
379   // special W medium for last 5 cm of W
380   zr=zRear-5;
381   cpar[0] = 2.5;
382   cpar[1] = zr * TMath::Tan(accMin);
383   cpar[2] = zr * TMath::Tan(thetaR*repsi);
384   cpar[3] = cpar[1] + TMath::Tan(accMin) * 5.;
385   cpar[4] = cpar[2] + TMath::Tan(thetaR*repsi) * 5.;
386   gMC->Gsvolu("ARW1", "CONE", idtmed[fMLayers[1][4]+20], cpar, 5);
387   dz=(dRear-epsi)/2.-cpar[0];
388   gMC->Gspos("ARW1", 1, "ARW0", 0., 0., dz, 0, "ONLY");
389   //
390   // Cu
391   Float_t drMin=TMath::Tan(thetaR) * 5;
392   Float_t drMax=TMath::Tan(accMax) * 5;
393   gMC->Gsvolu("ARPE", "CONE", idtmed[fMLayers[0][4]], cpar, 0);
394   cpar[0]=2.5;
395   { // Begin local scope for i
396       for (Int_t i=0; i<3; i++) {
397           zr=zRear-dRear+5+i*10.;
398           cpar[1] = zr * TMath::Tan(thetaR);
399           cpar[2] = zr * TMath::Tan(accMax);
400           cpar[3] = cpar[1] + drMin;
401           cpar[4] = cpar[2] + drMax;
402           dz=(zRear-zAbsStart)/2.-cpar[0]-5.-(2-i)*10;
403           gMC->Gsposp("ARPE", i+1, "AITR", 0., 0., dz, 0, "ONLY",cpar,5);
404       }
405   } // End local scope for i
406   gMC->Gspos("AITR", 1, "ABSS", 0., 0., 0., 0, "ONLY"); 
407   dz = (zRear-zAbsStart)/2.+zAbsStart;
408   gMC->Gspos("ABSM", 1, "ALIC", 0., 0., dz, 0, "ONLY"); 
409 //
410 //
411 // vacuum system
412 //
413 // pipe and heating jackets
414 //
415 //
416 // cylindrical piece
417   tpar0[2]=(zOpen-zAbsStart)/2;
418   tpar0[0]=rVacu;
419   tpar0[1]=rVacu+dTube+dInsu+dEnve;
420   gMC->Gsvolu("AV11", "TUBE", idtmed[kSteel+40], tpar0, 3);
421 //
422 // insulation
423
424   tpar[2]=tpar0[2];
425   tpar[0]=rVacu+dTube;
426   tpar[1]=tpar[0]+dInsu;
427   gMC->Gsvolu("AI11", "TUBE", idtmed[kInsulation+40], tpar, 3);
428   gMC->Gspos("AI11", 1, "AV11", 0., 0., 0., 0, "ONLY"); 
429 //
430   dz=-(zRear-zAbsStart)/2.+tpar0[2];
431   gMC->Gspos("AV11", 1, "ABSM", 0., 0., dz, 0, "ONLY"); 
432 //
433 // conical piece
434
435   cpar0[0]=(zRear-dRear-zOpen)/2;
436   cpar0[1]= rVacu-0.05;
437   cpar0[2]= rVacu+dTube+dInsu+dEnve;
438   Float_t dR=2.*cpar0[0]*TMath::Tan(thetaOpen1);
439   cpar0[3]=cpar0[1]+dR;
440   cpar0[4]=cpar0[2]+dR;
441   gMC->Gsvolu("AV21", "CONE", idtmed[kSteel+40], cpar0, 5);
442   dTube+=0.05;
443
444 //
445 // insulation
446   cpar[0]=cpar0[0];
447   cpar[1]=cpar0[1]+dTube;
448   cpar[2]=cpar0[1]+dTube+dInsu;
449   cpar[3]=cpar0[3]+dTube;
450   cpar[4]=cpar0[3]+dTube+dInsu;
451   gMC->Gsvolu("AI21", "CONE", idtmed[kInsulation+40], cpar, 5);
452   gMC->Gspos("AI21", 1, "AV21", 0., 0., 0., 0, "ONLY"); 
453   
454   dz=(zRear-zAbsStart)/2.-cpar0[0]-dRear;
455   gMC->Gspos("AV21", 1, "ABSM", 0., 0., dz, 0, "ONLY"); 
456 //
457 // Support cone 
458
459   par[0]  =  22.5;
460   par[1]  = 360.0;
461   par[2]  =   8.0;
462   par[3]  =   4.0;
463     
464   par[4]  = zRear;
465   par[5]  = 100.;
466   par[6]  = 180.;
467   
468   par[7]  = zRear+20.;
469   par[8]  = 100.;
470   par[9]  = 180.;
471
472   par[10] = zRear+20.;
473   par[11] = 178.;
474   par[12] = 180.;
475
476   par[13] = 600.;
477   par[14] = 178.;
478   par[15] = 180.;
479   
480
481   gMC->Gsvolu("ASSS", "PGON", idtmed[kAl], par, 16);
482   gMC->Gspos("ASSS", 1, "ALIC", 0., 0., 0., 0, "ONLY");
483
484   Float_t trap[11];
485   trap[ 0] = (530.-170.)/2.;
486   trap[ 2] = 0.;
487   trap[ 3] = 2.; 
488   trap[ 4] = (600.-(zRear+2.))/2.;;
489   trap[ 5] = trap[4];
490   trap[ 6] = 0.;
491   trap[ 7] = 2.;
492   trap[ 8] = 5.;
493   trap[ 9] = 5.;
494   trap[10] = 0.;
495   trap[ 1] = -TMath::ATan((trap[4]-trap[8])/2./trap[0])*180./TMath::Pi();
496   AliMatrix(idrotm[1600], 180., 0., 90., 0., 90., 90.);
497   AliMatrix(idrotm[1601], 180., 0., 90., 0., 90., 270.);
498   gMC->Gsvolu("ASST", "TRAP", idtmed[kSteel], trap, 11);
499   dz = (600.+zRear+2.)/2.+(trap[4]-trap[8])/2.;
500   Float_t dy =  170.+trap[0];
501   
502 //  gMC->Gspos("ASST", 1, "ALIC", 0.,  dy, dz, idrotm[1600], "ONLY");
503 //  gMC->Gspos("ASST", 2, "ALIC", 0., -dy, dz, idrotm[1601], "ONLY");
504 }
505
506 //_____________________________________________________________________________
507
508 void AliABSOv0::Init()
509 {
510   //
511   // Initialisation of the muon absorber after it has been built
512   Int_t i;
513   //
514   if(fDebug) {
515     printf("\n%s: ",ClassName());
516     for(i=0;i<35;i++) printf("*");
517     printf(" ABSOv0_INIT ");
518     for(i=0;i<35;i++) printf("*");
519     printf("\n%s: ",ClassName());
520     //
521     for(i=0;i<80;i++) printf("*");
522     printf("\n");
523   }
524 }
525  
526
527
528
529
530
531
532
533
534