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