800da27d2038a90bb2e658fe5146c8ac9d0f9cce
[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.3  2000/01/18 17:49:56  morsch
19 Serious overlap of ABSM with shield corrected
20 Small error in ARPB parameters corrected
21
22 Revision 1.2  2000/01/13 11:23:59  morsch
23 Last layer of Pb outer angle corrected
24
25 Revision 1.1  2000/01/12 15:39:30  morsch
26 Standar version of ABSO
27
28 */
29
30 ///////////////////////////////////////////////////////////////////////////////
31 //                                                                           //
32 //  Muon ABSOrber                                                            //
33 //  This class contains the description of the muon absorber geometry        //
34 //                                                                           //
35 //Begin_Html
36 /*
37 <img src="picts/AliABSOClass.gif">
38 </pre>
39 <br clear=left>
40 <font size=+2 color=red>
41 <p>The responsible person for this module is
42 <a href="mailto:andreas.morsch@cern.ch">Andreas Morsch</a>.
43 </font>
44 <pre>
45 */
46 //End_Html
47 //                                                                           //
48 //                                                                           //
49 ///////////////////////////////////////////////////////////////////////////////
50
51 #include "AliABSOv0.h"
52 #include "AliRun.h"
53 #include "AliConst.h"
54  
55 ClassImp(AliABSOv0)
56  
57 //_____________________________________________________________________________
58 AliABSOv0::AliABSOv0()
59 {
60   //
61   // Default constructor
62   //
63 }
64  
65 //_____________________________________________________________________________
66 AliABSOv0::AliABSOv0(const char *name, const char *title)
67        : AliABSO(name,title)
68 {
69   //
70   // Standard constructor
71   //
72   SetMarkerColor(7);
73   SetMarkerStyle(2);
74   SetMarkerSize(0.4);
75 }
76  
77 //_____________________________________________________________________________
78 void AliABSOv0::CreateGeometry()
79 {
80   //
81   // Creation of the geometry of the muon absorber
82   //
83   //Begin_Html
84   /*
85     <img src="picts/AliABSOv0Tree.gif">
86   */
87   //End_Html
88   //Begin_Html
89   /*
90     <img src="picts/AliABSOv0.gif">
91   */
92   //End_Html
93
94
95   Int_t *idtmed = fIdtmed->GetArray()-1599;
96   
97   Float_t par[24], cpar[5], cpar0[5], pcpar[12], tpar[3], tpar0[3]; 
98   Float_t dz;
99 #include "ShieldConst.h"
100 // Mother volume and outer shielding: Pb
101     
102   par[0]  = 0.;
103   par[1]  = 360.;
104   par[2]  = 7.;
105
106   par[3]  = -(abs_l-abs_d)/2.;
107   par[4]  = r_abs;
108   par[5]  = abs_d * TMath::Tan(theta1);
109
110   par[6]  = par[3]+(z_nose-abs_d);
111   par[7]  = r_abs;
112   par[8]  = z_nose * TMath::Tan(theta1);
113
114   par[9]  = par[3]+(z_cone-abs_d);
115   par[10] = r_abs;
116   par[11] = par[8] + (par[9] - par[6]) * TMath::Tan(theta2);
117
118   par[12]  = par[3]+(abs_c-abs_d);
119   par[13] = r_abs;
120   par[14] = par[11] + (par[12] - par[9]) * TMath::Tan(acc_max);
121
122   par[15] = par[3]+(abs_l-d_rear-abs_d);
123   par[16] = r_abs   + (par[15] - par[12]) * TMath::Tan(theta_open1) ;
124   par[17] = par[14] + (par[15] - par[12]) * TMath::Tan(acc_max);
125
126   par[18] = par[3]+(abs_l-d_rear-abs_d);
127   par[19] = (abs_l-d_rear) * TMath::Tan(acc_min);
128   par[20] = par[14] + (par[18] - par[12]) * TMath::Tan(acc_max);
129
130   par[21] = -par[3];
131   par[22] =  abs_l* TMath::Tan(acc_min);
132   par[23] = par[20] + (par[21] - par[18]) * TMath::Tan(acc_max);
133   gMC->Gsvolu("ABSS", "PCON", idtmed[1612], par, 24);
134   { // Begin local scope for i
135     for (Int_t i=4; i<18; i+=3) par[i]  = 0;
136   } // End local scope for i
137   gMC->Gsvolu("ABSM", "PCON", idtmed[1655], par, 24);
138   gMC->Gspos("ABSS", 1, "ABSM", 0., 0., 0., 0, "ONLY");
139
140 //
141 // Steel envelope
142 //
143   par[4] = par[5] -d_steel;
144   par[7] = par[8] -d_steel;
145   par[10]= par[11]-d_steel;  
146   par[13]= par[14]-d_steel;  
147   par[16]= par[17]-d_steel;  
148   par[19]= par[20]-d_steel;  
149   par[22]= par[23]-d_steel;  
150   gMC->Gsvolu("ABST", "PCON", idtmed[1618], par, 24);
151   gMC->Gspos("ABST", 1, "ABSS", 0., 0., 0., 0, "ONLY");
152 //
153 // Polyethylene shield
154 // 
155   cpar[0] = (abs_l - z_cone) / 2.;
156   cpar[1] = z_cone * TMath::Tan(acc_max);
157   cpar[2] = cpar[1] + d_poly;
158   cpar[3] = abs_l * TMath::Tan(acc_max);
159   cpar[4] = cpar[3] + d_poly;
160   gMC->Gsvolu("APOL", "CONE", idtmed[1657], cpar, 5);
161   dz = (abs_l-abs_d)/2.-cpar[0];
162   gMC->Gspos("APOL", 1, "ABSS", 0., 0., dz, 0, "ONLY");
163
164 //
165 // Tungsten nose to protect TPC
166 // 
167   cpar[0] = (z_nose - abs_d) / 2.;
168   cpar[1] = abs_d * TMath::Tan(acc_max);
169   cpar[2] = abs_d * TMath::Tan(theta1)-d_steel;
170   cpar[3] = z_nose * TMath::Tan(acc_max);
171   cpar[4] = z_nose * TMath::Tan(theta1)-d_steel;
172   gMC->Gsvolu("ANOS", "CONE", idtmed[1611], cpar, 5);
173   //
174   dz = -(abs_l-abs_d)/2.+cpar[0];
175   gMC->Gspos("ANOS", 1, "ABSS", 0., 0., dz, 0, "ONLY");
176 //
177 // Tungsten inner shield
178 //
179   cpar[0] = (abs_l-d_rear - abs_c)/ 2.;
180   cpar[1] = r_abs;
181   cpar[2] = abs_c * TMath::Tan(acc_min);
182   cpar[3] = r_abs + 2. * cpar[0] * TMath::Tan(theta_open1);
183   cpar[4] = (abs_l-d_rear)  * TMath::Tan(acc_min);
184   gMC->Gsvolu("AWIN", "CONE", idtmed[1651], cpar, 5);
185   //
186   dz = (abs_l-abs_d)/2.-cpar[0]-d_rear;
187   gMC->Gspos("AWIN", 1, "ABSS", 0., 0., dz, 0, "ONLY");
188
189   //     Inner tracking region
190   //
191   //     mother volume: Pb
192   //
193   pcpar[0]  = 0.;
194   pcpar[1]  = 360.;
195   pcpar[2]  = 3.;
196   pcpar[3]  = -(abs_l-abs_d)/2.;
197   pcpar[4]  = r_abs;
198   pcpar[5]  = abs_d * TMath::Tan(acc_max);
199   pcpar[6]  = pcpar[3]+(z_2deg-abs_d);
200   pcpar[7]  = r_abs;
201   pcpar[8]  = z_2deg * TMath::Tan(acc_max);
202   pcpar[9]  = -pcpar[3];
203   pcpar[10] = abs_l * TMath::Tan(acc_min);
204   pcpar[11] = abs_l * TMath::Tan(acc_max);
205   gMC->Gsvolu("AITR", "PCON", idtmed[1612], pcpar, 12);
206   //
207   // special Pb medium for last 5 cm of Pb
208   zr=abs_l-5;
209   cpar[0] = 2.5;
210   cpar[1] = zr * TMath::Tan(theta_r);
211   cpar[2] = zr * TMath::Tan(acc_max);
212   cpar[3] = cpar[1] + TMath::Tan(theta_r) * 5;
213   cpar[4] = cpar[2] + TMath::Tan(acc_max) * 5;
214   gMC->Gsvolu("ARPB", "CONE", idtmed[1632], cpar, 5);
215   dz=(abs_l-abs_d)/2.-cpar[0];
216   gMC->Gspos("ARPB", 1, "AITR", 0., 0., dz, 0, "ONLY");
217
218   //
219   //     concrete cone: concrete 
220   //
221   pcpar[9]  = par[3]+(abs_l-d_rear-abs_d);
222   pcpar[10] = (abs_l-d_rear) * TMath::Tan(acc_min);
223   pcpar[11] = (abs_l-d_rear) * TMath::Tan(acc_max);
224   gMC->Gsvolu("ACON", "PCON", idtmed[1616], pcpar, 12);
225   gMC->Gspos("ACON", 1, "AITR", 0., 0., 0., 0, "ONLY");
226   //
227   //     carbon cone: carbon
228   //
229   pcpar[9]  = pcpar[3]+(abs_cc-abs_d);
230   pcpar[10]  = abs_cc * TMath::Tan(acc_min);
231   pcpar[11]  = abs_cc * TMath::Tan(acc_max);
232   gMC->Gsvolu("ACAR", "PCON", idtmed[1605], pcpar, 12);
233   gMC->Gspos("ACAR", 1, "ACON", 0., 0., 0., 0, "ONLY");
234   //
235   //     inner W shield
236   zr=abs_l-d_rear;
237   cpar[0] = d_rear/2.;
238   cpar[1] = zr * TMath::Tan(acc_min);
239   cpar[2] = zr * TMath::Tan(theta_r);
240   cpar[3] = cpar[1] + TMath::Tan(acc_min) * d_rear;
241   cpar[4] = cpar[2] + TMath::Tan(theta_r) * d_rear;
242   gMC->Gsvolu("ARW0", "CONE", idtmed[1611], cpar, 5);
243   dz=(abs_l-abs_d)/2.-cpar[0];
244   gMC->Gspos("ARW0", 1, "AITR", 0., 0., dz, 0, "ONLY");
245   //
246   // special W medium for last 5 cm of W
247   zr=abs_l-5;
248   cpar[0] = 2.5;
249   cpar[1] = zr * TMath::Tan(acc_min);
250   cpar[2] = zr * TMath::Tan(theta_r);
251   cpar[3] = cpar[1] + TMath::Tan(acc_min) * 5.;
252   cpar[4] = cpar[2] + TMath::Tan(theta_r) * 5.;
253   gMC->Gsvolu("ARW1", "CONE", idtmed[1631], cpar, 5);
254   dz=d_rear/2.-cpar[0];
255   gMC->Gspos("ARW1", 1, "ARW0", 0., 0., dz, 0, "ONLY");
256   //
257   // PolyEthylene Layers
258   Float_t dr_min=TMath::Tan(theta_r) * 5;
259   Float_t dr_max=TMath::Tan(acc_max) * 5;
260   gMC->Gsvolu("ARPE", "CONE", idtmed[1617], cpar, 0);
261   cpar[0]=2.5;
262   { // Begin local scope for i
263     for (Int_t i=0; i<3; i++) {
264       zr=abs_l-d_rear+5+i*10.;
265       cpar[1] = zr * TMath::Tan(theta_r);
266       cpar[2] = zr * TMath::Tan(acc_max);
267       cpar[3] = cpar[1] + dr_min;
268       cpar[4] = cpar[2] + dr_max;
269       dz=(abs_l-abs_d)/2.-cpar[0]-5.-(2-i)*10;
270       gMC->Gsposp("ARPE", i+1, "AITR", 0., 0., dz, 0, "ONLY",cpar,5);
271     }
272   } // End local scope for i
273   gMC->Gspos("AITR", 1, "ABSS", 0., 0., 0., 0, "ONLY"); 
274   dz = (abs_l-abs_d)/2.+abs_d;
275   gMC->Gspos("ABSM", 1, "ALIC", 0., 0., dz, 0, "ONLY"); 
276 //
277 //
278 // vacuum system
279 //
280 // pipe and heating jackets
281 //
282 //
283 // cylindrical piece
284   tpar0[2]=(abs_c-abs_d)/2;
285   tpar0[0]=r_vacu;
286   tpar0[1]=r_abs;
287   gMC->Gsvolu("AV11", "TUBE", idtmed[1658], tpar0, 3);
288 //
289 // insulation
290   tpar[2]=tpar0[2];
291   tpar[0]=tpar0[0]+d_tube;
292   tpar[1]=tpar0[0]+d_tube+d_insu;
293   gMC->Gsvolu("AI11", "TUBE", idtmed[1653], tpar, 3);
294   gMC->Gspos("AI11", 1, "AV11", 0., 0., 0., 0, "ONLY"); 
295 //
296 // clearance
297   tpar[0]=tpar0[1]-d_prot-d_free;
298   tpar[1]=tpar0[1]-d_prot;
299   gMC->Gsvolu("AP11", "TUBE", idtmed[1655], tpar, 3);
300   gMC->Gspos("AP11", 1, "AV11", 0., 0., 0., 0, "ONLY"); 
301   
302   dz=-(abs_l-abs_d)/2.+tpar0[2];
303   gMC->Gspos("AV11", 1, "ABSM", 0., 0., dz, 0, "ONLY"); 
304   
305 //
306 // conical piece
307   cpar0[0]=(abs_l-d_rear-abs_c)/2;
308   cpar0[1]=r_vacu;
309   cpar0[2]=r_abs;
310   cpar0[3]=cpar0[1]+2.*cpar0[0]*TMath::Tan(theta_open1);
311   cpar0[4]=cpar0[2]+2.*cpar0[0]*TMath::Tan(theta_open1);
312   gMC->Gsvolu("AV21", "CONE", idtmed[1658], cpar0, 5);
313 //
314 // insulation
315   cpar[0]=cpar0[0];
316   cpar[1]=cpar0[1]+d_tube;
317   cpar[2]=cpar0[1]+d_tube+d_insu;
318   cpar[3]=cpar0[3]+d_tube;
319   cpar[4]=cpar0[3]+d_tube+d_insu;
320   gMC->Gsvolu("AI21", "CONE", idtmed[1653], cpar, 5);
321   gMC->Gspos("AI21", 1, "AV21", 0., 0., 0., 0, "ONLY"); 
322 //
323 // clearance
324   cpar[1]=cpar0[2]-d_prot-d_free;
325   cpar[2]=cpar0[2]-d_prot;
326   cpar[3]=cpar0[4]-d_prot-d_free;
327   cpar[4]=cpar0[4]-d_prot;
328   gMC->Gsvolu("AP21", "CONE", idtmed[1655], cpar, 5);
329   gMC->Gspos("AP21", 1, "AV21", 0., 0., 0., 0, "ONLY"); 
330   
331   dz=(abs_l-abs_d)/2.-cpar0[0]-d_rear;
332   gMC->Gspos("AV21", 1, "ABSM", 0., 0., dz, 0, "ONLY"); 
333
334 }
335
336 //_____________________________________________________________________________
337
338 void AliABSOv0::Init()
339 {
340   //
341   // Initialisation of the muon absorber after it has been built
342   Int_t i;
343   //
344   printf("\n");
345   for(i=0;i<35;i++) printf("*");
346   printf(" ABSOv0_INIT ");
347   for(i=0;i<35;i++) printf("*");
348   printf("\n");
349   //
350   for(i=0;i<80;i++) printf("*");
351   printf("\n");
352 }
353  
354
355
356
357
358
359
360
361
362