]>
Commit | Line | Data |
---|---|---|
4c039060 | 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$ | |
b16a1b1e | 18 | Revision 1.23 2002/10/14 14:57:39 hristov |
19 | Merging the VirtualMC branch to the main development branch (HEAD) | |
20 | ||
b9d0a01d | 21 | Revision 1.19.4.2 2002/10/11 08:04:28 hristov |
22 | Updating VirtualMC to v3-09-02 | |
23 | ||
24 | Revision 1.22 2002/09/02 15:29:20 morsch | |
25 | Gsbool calls to resolve MANY added (I. Hrivnacova). | |
26 | ||
fddcf8f7 | 27 | Revision 1.19.4.1 2002/06/10 15:13:48 hristov |
28 | Merged with v3-08-02 | |
29 | ||
30 | Revision 1.21 2002/05/11 19:14:44 morsch | |
31 | PCONE QBEU had last z defined twice. | |
32 | ||
13a638e8 | 33 | Revision 1.20 2002/05/02 12:36:43 morsch |
34 | New version of beam-pipe geometry. Ref. AL502206PL | |
35 | (used in production readiness report). | |
36 | ||
614dff1c | 37 | Revision 1.19 2002/03/25 12:34:40 morsch |
38 | Obsolete support rollers removed. | |
39 | ||
4a7308f7 | 40 | Revision 1.18 2002/02/01 18:02:41 morsch |
41 | Material of beam pipe between inner Be piece and forward detectors | |
42 | can be set by SetPipeMaterial(mat), mat = kInox, kAlu, kBe ... | |
43 | ||
765e9e1a | 44 | Revision 1.17 2001/09/24 13:11:50 morsch |
45 | Ion pump and bellows moved out by 15 cm to make space for forward | |
46 | detectors. | |
47 | ||
22811e34 | 48 | Revision 1.16 2001/05/16 14:57:22 alibrary |
49 | New files for folders and Stack | |
50 | ||
9e1a0ddb | 51 | Revision 1.15 2001/05/02 11:50:18 morsch |
52 | New layout of the non-absorber side provided by Y. Viyogi. Not the final design | |
53 | but the prsent most realistic. | |
54 | ||
cd9c4f3a | 55 | Revision 1.14 2001/01/20 16:56:33 morsch |
56 | Put air in connecting tubes and flanges of vacuum pump. | |
57 | ||
14b75006 | 58 | Revision 1.13 2001/01/20 16:35:27 morsch |
59 | Increase mother volume for bellows. | |
60 | ||
1416b7f3 | 61 | Revision 1.12 2000/12/21 16:41:06 morsch |
62 | Coding convention clean-up (RS3) | |
63 | ||
13391c3e | 64 | Revision 1.11 2000/11/28 16:06:57 morsch |
65 | Undulated beam-pipe replaced by Al-Be (40,60) pipe 1.5 mm thick. | |
66 | ||
2196bc42 | 67 | Revision 1.10 2000/11/24 13:00:37 morsch |
68 | - Geometry and materials imported from euclid output | |
69 | - include comments | |
70 | - better struturing of volume tree | |
71 | - improved version of flange close to front absorber | |
72 | - more realistic pump materials | |
73 | - undulated beam pipe imported from v3. | |
74 | ||
7bd7faf8 | 75 | Revision 1.9 2000/10/02 21:28:15 fca |
76 | Removal of useless dependecies via forward declarations | |
77 | ||
94de3818 | 78 | Revision 1.8 2000/06/11 12:37:01 morsch |
79 | Coding rule violations corrected | |
80 | ||
b43eb0dc | 81 | Revision 1.7 2000/02/23 16:25:24 fca |
82 | AliVMC and AliGeant3 classes introduced | |
83 | ReadEuclid moved from AliRun to AliModule | |
84 | ||
b13db077 | 85 | Revision 1.6 1999/09/29 09:24:30 fca |
86 | Introduction of the Copyright and cvs Log | |
87 | ||
4c039060 | 88 | */ |
89 | ||
fe4da5cc | 90 | //////////////////////////////////////////////// |
91 | // Beam pipe class / | |
92 | //////////////////////////////////////////////// | |
93 | ||
fe4da5cc | 94 | #include "AliPIPEv0.h" |
95 | #include "AliRun.h" | |
7bd7faf8 | 96 | #include "AliConst.h" |
97 | #include "AliMagF.h" | |
94de3818 | 98 | #include "AliMC.h" |
fe4da5cc | 99 | #include "TSystem.h" |
22811e34 | 100 | |
b16a1b1e | 101 | #include <Riostream.h> |
fe4da5cc | 102 | |
103 | ClassImp(AliPIPEv0) | |
104 | ||
105 | //_____________________________________________________________________________ | |
b8032157 | 106 | AliPIPEv0::AliPIPEv0() |
fe4da5cc | 107 | { |
b43eb0dc | 108 | // Constructor |
765e9e1a | 109 | SetPipeMaterial(); |
fe4da5cc | 110 | } |
111 | ||
112 | //_____________________________________________________________________________ | |
113 | AliPIPEv0::AliPIPEv0(const char *name, const char *title) | |
114 | : AliPIPE(name,title) | |
115 | { | |
b43eb0dc | 116 | // Constructor |
765e9e1a | 117 | SetPipeMaterial(); |
fe4da5cc | 118 | } |
119 | ||
120 | ||
121 | //___________________________________________ | |
122 | void AliPIPEv0::CreateGeometry() | |
123 | { | |
fe4da5cc | 124 | //Begin_Html |
125 | /* | |
1439f98e | 126 | <img src="picts/pipe.gif"> |
fe4da5cc | 127 | */ |
128 | //End_Html | |
129 | ||
130 | ||
131 | //Begin_Html | |
132 | /* | |
1439f98e | 133 | <img src="picts/tree_pipe.gif"> |
fe4da5cc | 134 | */ |
135 | //End_Html | |
136 | ||
9e1a0ddb | 137 | if(fDebug) printf("%s: Create PIPEv0 geometry \n",ClassName()); |
fe4da5cc | 138 | |
fe4da5cc | 139 | |
7bd7faf8 | 140 | Int_t *idtmed = fIdtmed->GetArray(); |
22811e34 | 141 | Float_t ppcon[84], ptube[3], pbox[3]; |
7bd7faf8 | 142 | Int_t i=0; |
143 | ||
7bd7faf8 | 144 | |
145 | Int_t idrotm[2099]; | |
146 | AliMatrix(idrotm[2001],90.,240., 0., 0., 90.,150.); | |
147 | AliMatrix(idrotm[2002],90., 0., 0., 0., 90.,270.); | |
148 | AliMatrix(idrotm[2003],90.,120., 0., 0., 90., 30.); | |
149 | AliMatrix(idrotm[2004],90.,315., 90., 45., 0., 0.); | |
150 | AliMatrix(idrotm[2005],90.,270., 90., 0., 0., 0.); | |
151 | AliMatrix(idrotm[2006],90.,225., 90.,315., 0., 0.); | |
152 | AliMatrix(idrotm[2007],90.,180., 90.,270., 0., 0.); | |
153 | AliMatrix(idrotm[2008],90.,135., 90.,225., 0., 0.); | |
154 | AliMatrix(idrotm[2009],90., 90., 90.,180., 0., 0.); | |
155 | AliMatrix(idrotm[2010],90., 45., 90.,135., 0., 0.); | |
156 | idrotm[2011] = 0; | |
157 | AliMatrix(idrotm[2012],90.,180., 90., 90.,180., 0.); | |
158 | AliMatrix(idrotm[2013],90., 0., 90., 90.,180., 0.); | |
614dff1c | 159 | // |
160 | // Bellow | |
161 | // | |
22811e34 | 162 | // distance between bellows |
614dff1c | 163 | const Float_t dzbb = 9.0; |
164 | // total size of bellow section | |
165 | const Float_t dzb = 15.0; | |
166 | // size of undulated region | |
167 | const Float_t dzub = 2.0; | |
168 | ||
169 | // half-lengths of various beam pipe sections | |
170 | // central Be-Pipe | |
171 | Float_t hlenQbbe1 = 40.; | |
172 | Float_t hlenQbbe2 = 36.5; | |
173 | Float_t hlenQbbe = (hlenQbbe1+hlenQbbe2)/2.; | |
174 | // | |
175 | // | |
176 | // Float_t hlenQbt1 = 5.5/2.; | |
177 | // | |
178 | // Pipe outside central region (non-absober side) | |
179 | Float_t hlenQbab = 157.5; | |
180 | // | |
181 | // Flange non-absorber side | |
182 | Float_t hlenQb29 = 11.5/2.+1.75 + 5.0; | |
183 | // | |
184 | // Bellow element | |
185 | Float_t hlenQbe0 = dzb; | |
186 | // | |
187 | // Inox pipe between Be and Bellow (absorber side) | |
188 | Float_t hlenQb24[3] = {10.5/2., 1.8, 3.3}; | |
189 | // | |
190 | // | |
191 | Float_t hlenQb28 = (800.-hlenQbbe1-2.*hlenQbab-4.*hlenQb29-2.*hlenQbe0)/2.; | |
192 | // | |
193 | // Position of the pump | |
194 | Float_t zPump = hlenQbbe1+2.*hlenQbab+2.*hlenQb29+dzb; | |
195 | // | |
196 | // Inner beam pipe radius | |
197 | Float_t RinBe = 2.9; | |
198 | Float_t RinSt = 2.92; | |
199 | // | |
200 | // | |
201 | Float_t RoutBe = 2.98; | |
202 | Float_t RoutSt = 3.00; | |
22811e34 | 203 | |
614dff1c | 204 | |
205 | // | |
206 | Float_t dz; | |
7bd7faf8 | 207 | |
fe4da5cc | 208 | // |
209 | // The peam pipe up to the Front Absorber | |
fe4da5cc | 210 | // |
7bd7faf8 | 211 | // Mother Volume QBPM |
614dff1c | 212 | ppcon[0] = 0; |
213 | ppcon[1] = 360; | |
214 | ppcon[2] = 18; | |
7bd7faf8 | 215 | // 1: |
614dff1c | 216 | ppcon[3] = -90.; |
217 | ppcon[4] = 0.; | |
218 | ppcon[5] = 4.4; | |
7bd7faf8 | 219 | // 2 |
614dff1c | 220 | ppcon[6] = -90+2.*hlenQb24[2]+2.8+2.*hlenQb24[1]; |
221 | ppcon[7] = 0.; | |
222 | ppcon[8] = 4.4; | |
7bd7faf8 | 223 | // 3 |
614dff1c | 224 | ppcon[9] = ppcon[6]; |
225 | ppcon[10] = 0.; | |
226 | ppcon[11] = 4.1; | |
22811e34 | 227 | // 4 |
614dff1c | 228 | ppcon[12] = ppcon[9] + 2. + 2.*dzub+0.2; |
229 | ppcon[13] = 0.; | |
230 | ppcon[14] = 4.1; | |
22811e34 | 231 | // 5 |
614dff1c | 232 | ppcon[15] = ppcon[12]; |
233 | ppcon[16] = 0.; | |
234 | ppcon[17] = 3.2; | |
22811e34 | 235 | // 6 |
614dff1c | 236 | ppcon[18] = ppcon[15] + 2.* dzbb-0.4; |
237 | ppcon[19] = 0.; | |
238 | ppcon[20] = 3.2; | |
7bd7faf8 | 239 | // 7 |
614dff1c | 240 | ppcon[21] = ppcon[18]; |
241 | ppcon[22] = 0.; | |
242 | ppcon[23] = 4.1; | |
7bd7faf8 | 243 | // 8 |
614dff1c | 244 | ppcon[24] = -44.; |
245 | ppcon[25] = 0.; | |
246 | ppcon[26] = 4.1; | |
7bd7faf8 | 247 | // 9 |
614dff1c | 248 | ppcon[27] = -44.; |
7bd7faf8 | 249 | ppcon[28] = 0; |
614dff1c | 250 | ppcon[29] = 3.0; |
22811e34 | 251 | // 10 |
614dff1c | 252 | ppcon[30] = 38.; |
7bd7faf8 | 253 | ppcon[31] = 0; |
614dff1c | 254 | ppcon[32] = 3.0; |
22811e34 | 255 | // 11 |
614dff1c | 256 | ppcon[33] = 38.; |
7bd7faf8 | 257 | ppcon[34] = 0; |
614dff1c | 258 | ppcon[35] = 3.6; |
22811e34 | 259 | // 12 |
614dff1c | 260 | ppcon[36] = hlenQbbe1+2.*hlenQbab-0.1; |
261 | ppcon[37] = 0.; | |
262 | ppcon[38] = 3.6; | |
263 | // 13 | |
264 | ppcon[39] = ppcon[36]; | |
265 | ppcon[40] = 0.; | |
266 | ppcon[41] = 3.6; | |
22811e34 | 267 | // 14 |
614dff1c | 268 | ppcon[42] = ppcon[39]+2.*hlenQb29-5.; |
269 | ppcon[43] = 0.; | |
270 | ppcon[44] = 3.6; | |
271 | // 15 | |
272 | ppcon[45] = ppcon[42]; | |
273 | ppcon[46] = 0.; | |
274 | ppcon[47] = 56.; | |
22811e34 | 275 | // 16 |
614dff1c | 276 | ppcon[48] = ppcon[45]+2.*dzb+10.; |
277 | ppcon[49] = 0.; | |
278 | ppcon[50] = 56.; | |
22811e34 | 279 | // 17 |
614dff1c | 280 | ppcon[51] = ppcon[48]; |
22811e34 | 281 | ppcon[52] = 0.; |
614dff1c | 282 | ppcon[53] = 3.6; |
22811e34 | 283 | // 18 |
614dff1c | 284 | ppcon[54] = 800.; |
22811e34 | 285 | ppcon[55] = 0.; |
614dff1c | 286 | ppcon[56] = 3.6; |
7bd7faf8 | 287 | |
614dff1c | 288 | gMC->Gsvolu("QBPM", "PCON", idtmed[kAir], ppcon,57); |
289 | ||
290 | ||
7bd7faf8 | 291 | // |
cd9c4f3a | 292 | // volume definitions of various sections |
293 | // | |
cd9c4f3a | 294 | |
614dff1c | 295 | // |
296 | // The Vacuum | |
297 | gMC->Gsvolu("QBVA","TUBE", idtmed[kVac], ptube, 0); | |
298 | ptube[0] = 0.0; | |
299 | ptube[1] = RinSt; | |
300 | ptube[2] = (90.-hlenQbbe2)/2.; | |
301 | dz = -90. + ptube[2]; | |
302 | gMC->Gsposp ("QBVA", 1, "QBPM", 0., 0., dz , 0, "ONLY", ptube, 3); | |
303 | dz = dz + ptube[2]; | |
304 | ||
305 | ptube[1] = RinBe; | |
306 | ptube[2] = hlenQbbe+hlenQbab; | |
307 | dz = dz + ptube[2]; | |
308 | gMC->Gsposp ("QBVA", 2, "QBPM", 0., 0., dz , 0, "ONLY", ptube, 3); | |
309 | dz = dz + ptube[2]; | |
310 | ||
311 | ptube[1] = RinSt; | |
312 | ptube[2] = (800.-hlenQbbe1-2.*hlenQbab)/2.; | |
313 | dz = dz + ptube[2]; | |
314 | gMC->Gsposp ("QBVA", 3, "QBPM", 0., 0., dz , 0, "ONLY", ptube, 3); | |
cd9c4f3a | 315 | |
614dff1c | 316 | // |
317 | // Be Pipe in central Alice | |
318 | ptube[0] = RinBe; | |
319 | ptube[1] = RoutBe; | |
cd9c4f3a | 320 | ptube[2] = hlenQbbe; |
7bd7faf8 | 321 | |
322 | gMC->Gsvolu("QBBE","TUBE", idtmed[kBe], ptube, 3); | |
7bd7faf8 | 323 | |
324 | // | |
614dff1c | 325 | // Support Ring |
7bd7faf8 | 326 | // |
614dff1c | 327 | // Mother |
328 | ppcon[0] = 0; | |
329 | ppcon[1] = 360; | |
330 | ppcon[2] = 6; | |
331 | // 1: | |
332 | ppcon[3] = -1.; | |
333 | ppcon[4] = 3.0; | |
334 | ppcon[5] = 3.4; | |
335 | // 2 | |
336 | ppcon[6] = -0.8; | |
337 | ppcon[7] = 3.0; | |
338 | ppcon[8] = 3.4; | |
339 | // 3 | |
340 | ppcon[9] = -0.8; | |
341 | ppcon[10] = 3.0; | |
342 | ppcon[11] = 3.2; | |
343 | // 4 | |
344 | ppcon[12] = 0.8; | |
345 | ppcon[13] = 3.0; | |
346 | ppcon[14] = 3.2; | |
347 | // 5 | |
348 | ppcon[15] = 0.8; | |
349 | ppcon[16] = 3.0; | |
350 | ppcon[17] = 3.4; | |
351 | // 6 | |
352 | ppcon[18] = 1.0; | |
353 | ppcon[19] = 3.0; | |
354 | ppcon[20] = 3.4; | |
7bd7faf8 | 355 | |
7bd7faf8 | 356 | |
614dff1c | 357 | gMC->Gsvolu("QBSR", "PCON", idtmed[kC], ppcon,21); |
358 | gMC->Gspos("QBSR", 1, "QBPM", 0.0, 0.0, 40., 0, "ONLY"); | |
359 | gMC->Gspos("QBSR", 2, "QBPM", 0.0, 0.0, 150., 0, "ONLY"); | |
360 | gMC->Gspos("QBSR", 3, "QBPM", 0.0, 0.0, 260., 0, "ONLY"); | |
361 | gMC->Gspos("QBSR", 4, "QBPM", 0.0, 0.0,- 46., 0, "ONLY"); | |
7bd7faf8 | 362 | // |
614dff1c | 363 | // Flange and Fixed Point: non absorber side |
7bd7faf8 | 364 | // |
614dff1c | 365 | // ----------> |
366 | // | |
367 | // Mother | |
368 | ppcon[0] = 0; | |
369 | ppcon[1] = 360; | |
370 | ppcon[2] = 4; | |
371 | // 1: | |
372 | ppcon[3] = -hlenQb29; | |
373 | ppcon[4] = RinSt; | |
374 | ppcon[5] = 5.8; | |
375 | // 2 | |
376 | ppcon[6] = ppcon[3]+3.6; | |
377 | ppcon[7] = RinSt; | |
378 | ppcon[8] = 5.8; | |
379 | // 3 | |
380 | ppcon[9] = ppcon[6]; | |
381 | ppcon[10] = RinSt; | |
382 | ppcon[11] = 3.6; | |
383 | // 4 | |
384 | ppcon[12] = hlenQb29; | |
385 | ppcon[13] = RinSt; | |
386 | ppcon[14] = 3.6; | |
7bd7faf8 | 387 | |
614dff1c | 388 | gMC->Gsvolu("QB29", "PCON", idtmed[kAir], ppcon,15); |
7bd7faf8 | 389 | |
7bd7faf8 | 390 | |
614dff1c | 391 | // Flange |
392 | ptube[0] = RinSt; | |
393 | ptube[1] = 5.7; | |
394 | ptube[2] = 1.75; | |
395 | gMC->Gsvolu("QF29","TUBE", idtmed[kInox], ptube, 3); | |
396 | gMC->Gspos("QF29", 1, "QB29", 0.0, 0.0, -hlenQb29+1.75, 0, "ONLY"); | |
397 | // Pipe | |
398 | ptube[0] = RinSt; | |
399 | ptube[1] = 3.0; | |
cd9c4f3a | 400 | ptube[2] = hlenQb29; |
614dff1c | 401 | gMC->Gsvolu("QS29","TUBE", idtmed[kInox], ptube, 3); |
402 | gMC->Gspos("QS29", 1, "QB29", 0.0, 0.0, 0., 0, "ONLY"); | |
403 | // Fixed point | |
404 | ptube[0] = RinSt; | |
405 | ptube[1] = 3.5; | |
406 | ptube[2] = 0.3; | |
407 | gMC->Gsvolu("QP29","TUBE", idtmed[kInox], ptube, 3); | |
408 | gMC->Gspos("QP29", 1, "QB29", 0.0, 0.0, -hlenQb29+9.75+3., 0, "ONLY"); | |
22811e34 | 409 | |
614dff1c | 410 | // |
7bd7faf8 | 411 | // |
cd9c4f3a | 412 | // Inox beam pipe: final section on non-absorber side |
7bd7faf8 | 413 | |
614dff1c | 414 | ptube[0] = RinSt; |
415 | ptube[1] = RoutSt; | |
cd9c4f3a | 416 | ptube[2] = hlenQb28; |
7bd7faf8 | 417 | |
418 | gMC->Gsvolu("QB28","TUBE", idtmed[kInox], ptube, 3); | |
7bd7faf8 | 419 | |
2196bc42 | 420 | |
421 | // Al-Be (40-60 wgt%, rho=2.7 g/cm**3) beam pipe | |
422 | // | |
614dff1c | 423 | // This section is under study (A.M. 1/2/2002) |
765e9e1a | 424 | // |
425 | ||
614dff1c | 426 | ptube[0] = RinBe; |
765e9e1a | 427 | if (fPipeMaterial == kAlu) { |
428 | ptube[1] = 3.0; | |
429 | } else if (fPipeMaterial == kBe) { | |
614dff1c | 430 | ptube[1] = RoutBe; |
765e9e1a | 431 | } else if (fPipeMaterial == kInox){ |
614dff1c | 432 | ptube[1] = RoutSt; |
765e9e1a | 433 | } |
cd9c4f3a | 434 | ptube[2] = hlenQbab; |
2196bc42 | 435 | |
765e9e1a | 436 | gMC->Gsvolu("QBAB","TUBE", idtmed[fPipeMaterial], ptube, 3); |
2196bc42 | 437 | |
cd9c4f3a | 438 | // 2.5 mm thick SS tube for hanging pump |
614dff1c | 439 | /* |
440 | ptube[0] = Rin; | |
7bd7faf8 | 441 | ptube[1] = 3.15; |
cd9c4f3a | 442 | ptube[2] = hlenQb26; |
7bd7faf8 | 443 | |
cd9c4f3a | 444 | gMC->Gsvolu("QB26","TUBE", idtmed[kInox], ptube, 3); |
614dff1c | 445 | */ |
7bd7faf8 | 446 | // |
447 | // Bellows | |
fe4da5cc | 448 | // |
22811e34 | 449 | // |
7bd7faf8 | 450 | // Mother Volume |
614dff1c | 451 | Float_t pconQBE0[33]; |
452 | pconQBE0[ 0]= 0; | |
453 | pconQBE0[ 1]= 360; | |
454 | pconQBE0[ 2]= 10; | |
455 | // 1 | |
456 | pconQBE0[ 3]= -dzb; | |
457 | pconQBE0[ 4]= RinSt; | |
458 | pconQBE0[ 5]= RoutSt; | |
459 | // 2 | |
460 | pconQBE0[ 6]= -dzb+2.; | |
461 | pconQBE0[ 7]= RinSt; | |
462 | pconQBE0[ 8]= RoutSt; | |
463 | // 3 | |
464 | pconQBE0[ 9]= -dzb+2.; | |
465 | pconQBE0[10]= RinSt; | |
466 | pconQBE0[11]= 4.00; | |
467 | // 4 | |
468 | pconQBE0[12]= -dzb+2.+2.*dzub; | |
469 | pconQBE0[13]= RinSt; | |
470 | pconQBE0[14]= 4.00; | |
471 | // 5 | |
472 | pconQBE0[15]= -dzb+2.+2.*dzub; | |
473 | pconQBE0[16]= RinSt; | |
474 | pconQBE0[17]= RoutSt; | |
475 | // 6 | |
476 | pconQBE0[18]= -dzb+2.+2.*dzub+2.*dzbb; | |
477 | pconQBE0[19]= RinSt; | |
478 | pconQBE0[20]= RoutSt; | |
479 | // 7 | |
480 | pconQBE0[21]= -dzb+2.+2.*dzub+2.*dzbb; | |
481 | pconQBE0[22]= RinSt; | |
482 | pconQBE0[23]= 4.00; | |
483 | // 8 | |
484 | pconQBE0[24]= -dzb+2.+4.*dzub+2.*dzbb; | |
485 | pconQBE0[25]= RinSt; | |
486 | pconQBE0[26]= 4.00; | |
487 | // 9 | |
488 | pconQBE0[27]= -dzb+2.+4.*dzub+2.*dzbb; | |
489 | pconQBE0[28]= RinSt; | |
490 | pconQBE0[29]= RoutSt; | |
491 | // 10 | |
492 | pconQBE0[30]= +dzb; | |
493 | pconQBE0[31]= RinSt; | |
494 | pconQBE0[32]= RoutSt; | |
495 | gMC->Gsvolu("QBE0", "PCON", idtmed[kAir], pconQBE0, 33); | |
496 | // | |
497 | // Undulated piece mother | |
498 | ptube[0] = RinSt; | |
499 | ptube[1] = 4.00; | |
500 | ptube[2] = dzub; | |
7bd7faf8 | 501 | gMC->Gsvolu("QBEM","TUBE", idtmed[kAir], ptube, 3); |
614dff1c | 502 | dz = -dzb+2.+dzub; |
503 | gMC->Gspos("QBEM", 2 ,"QBE0", 0.0, 0.0, dz, 0 , "ONLY"); | |
504 | gMC->Gspos("QBEM", 1 ,"QBE0", 0.0, 0.0, -dz, idrotm[2012], "ONLY"); | |
505 | // | |
506 | Float_t pund[30]; | |
507 | Float_t uw = 0.02; | |
508 | Float_t ur = 2.*dzub-12.*uw; | |
509 | Float_t uz = ur/13.; | |
510 | Float_t ut = uz+uw; | |
511 | ||
512 | pund[ 0] = 0; | |
513 | pund[ 1] = 360; | |
13a638e8 | 514 | pund[ 2] = 8; |
614dff1c | 515 | |
516 | pund[ 3] = -ut; | |
517 | pund[ 4] = RinSt; | |
518 | pund[ 5] = RinSt+uw; | |
519 | ||
520 | pund[ 6] = -ut+uz; | |
521 | pund[ 7] = pund[4]; | |
522 | pund[ 8] = pund[5]; | |
523 | ||
524 | pund[ 9] = pund[6]; | |
525 | pund[10] = pund[4]; | |
526 | pund[11] = 4.0; | |
527 | ||
528 | pund[12] = pund[9]+uw; | |
529 | pund[13] = pund[10]; | |
530 | pund[14] = pund[11]; | |
531 | ||
532 | pund[15] = pund[12]; | |
533 | pund[16] = 4.0-uw; | |
534 | pund[17] = pund[14]; | |
535 | ||
536 | pund[18] = pund[12]+uz; | |
537 | pund[19] = pund[16]; | |
538 | pund[20] = pund[17]; | |
539 | ||
540 | pund[21] = pund[18]; | |
541 | pund[22] = RinSt; | |
542 | pund[23] = pund[20]; | |
543 | ||
544 | pund[24] = pund[21]+uw; | |
545 | pund[25] = pund[22]; | |
546 | pund[26] = pund[23]; | |
614dff1c | 547 | |
13a638e8 | 548 | gMC->Gsvolu("QBEU", "PCON", idtmed[kInox], pund, 27); |
614dff1c | 549 | |
550 | for (i = 0; i < 6; i++) | |
551 | { | |
552 | dz = -dzub+(1+2*i)*ut; | |
553 | gMC->Gspos("QBEU", i+1 ,"QBEM", 0.0, 0.0, dz, 0 , "ONLY"); | |
7bd7faf8 | 554 | } |
614dff1c | 555 | ptube[0] = RinSt; |
556 | ptube[1] = RinSt+uw; | |
557 | ptube[2] = uz; | |
558 | gMC->Gsvolu("QBEW","TUBE", idtmed[kInox], ptube, 3); | |
559 | gMC->Gspos("QBEW", 1 ,"QBEM", 0.0, 0.0, dzub-uz, 0 , "ONLY"); | |
7bd7faf8 | 560 | // |
614dff1c | 561 | // BeamPipe |
562 | gMC->Gsvolu("QBEP","TUBE", idtmed[kInox], ptube, 0); | |
563 | ptube[0] = RinSt; | |
564 | ptube[1] = RoutSt; | |
565 | ptube[2] = 1.; | |
566 | gMC->Gsposp("QBEP", 1 ,"QBE0", 0.0, 0.0, -dzb+1., 0 , "ONLY", ptube, 3); | |
567 | gMC->Gsposp("QBEP", 2 ,"QBE0", 0.0, 0.0, dzb-1., 0 , "ONLY", ptube, 3); | |
568 | ptube[2] = dzbb; | |
569 | gMC->Gsposp("QBEP", 3 ,"QBE0", 0.0, 0.0, 0., 0 , "ONLY", ptube, 3); | |
570 | // | |
571 | // End undulated part | |
572 | // | |
573 | // | |
574 | // ----> End Bellow | |
575 | // | |
576 | // **** Placement of various sections on non-absorber side **** | |
7bd7faf8 | 577 | // |
cd9c4f3a | 578 | // |
614dff1c | 579 | // first the beryllium section |
580 | Float_t zpos = -(hlenQbbe2-hlenQbbe1)/2; | |
cd9c4f3a | 581 | gMC->Gspos("QBBE", 1, "QBPM", 0., 0., zpos, 0, "ONLY"); |
582 | ||
583 | // next meta-metal transition QBT1 on on-absorber side | |
614dff1c | 584 | // zpos = zpos + hlenQbbe + hlenQbt1; |
585 | // gMC->Gspos("QBT1", 1, "QBPM", 0., 0., zpos, 0, "ONLY"); | |
cd9c4f3a | 586 | |
587 | // Aluminium OR Al-be alloy section | |
614dff1c | 588 | zpos = hlenQbbe1+hlenQbab; |
cd9c4f3a | 589 | gMC->Gspos("QBAB", 1, "QBPM", 0.0, 0.0, zpos, 0, "ONLY"); |
614dff1c | 590 | // |
cd9c4f3a | 591 | // inox flange at the start of bellow |
cd9c4f3a | 592 | zpos = zpos + hlenQbab + hlenQb29; |
614dff1c | 593 | gMC->Gspos("QB29", 1, "QBPM", 0.0, 0.0, zpos, idrotm[2012], "ONLY"); |
594 | // | |
cd9c4f3a | 595 | // bellow section |
cd9c4f3a | 596 | zpos = zpos + hlenQb29 + hlenQbe0; |
614dff1c | 597 | gMC->Gspos("QBE0", 2 ,"QBPM", 0.0, 0.0, zpos, 0, "ONLY"); |
598 | // | |
cd9c4f3a | 599 | // inox flange at the end of bellow and start of thick inox for pump |
cd9c4f3a | 600 | zpos = zpos + hlenQbe0 + hlenQb29; |
614dff1c | 601 | gMC->Gspos("QB29", 2, "QBPM", 0.0, 0.0, zpos, 0, "ONLY"); |
602 | // | |
cd9c4f3a | 603 | //last inox section till 800 cm |
cd9c4f3a | 604 | zpos = zpos + hlenQb29 + hlenQb28; |
605 | gMC->Gspos("QB28", 1, "QBPM", 0.0, 0.0, zpos, 0, "ONLY"); | |
614dff1c | 606 | |
cd9c4f3a | 607 | //******** end of placement on no-absorber side ********* |
608 | ||
609 | // | |
610 | // **** Absorber side ***** | |
611 | // | |
612 | // | |
614dff1c | 613 | /* |
cd9c4f3a | 614 | // metal-metal transition : Be-Alu on absorber side |
615 | // Mother Volume | |
616 | ptube[0] = 2.900; | |
617 | ptube[1] = 4.200; | |
618 | ptube[2] = 2.750; | |
619 | gMC->Gsvolu("QBT2","TUBE", idtmed[kAir], ptube, 3); | |
620 | // z = 43.3 - 48.8 | |
621 | gMC->Gspos("QBT2", 1, "QBPM", 0., 0., -hlenQbbe-ptube[2], idrotm[2012], "ONLY"); | |
622 | ||
623 | ptube[0] = 2.900; | |
624 | ptube[1] = 3.150; | |
625 | ptube[2] = 0.375; | |
626 | // | |
627 | // Be-part | |
628 | gMC->Gsvolu("QB02","TUBE", idtmed[kAlu], ptube, 3); | |
629 | ||
630 | ptube[1] = 3.000; | |
631 | gMC->Gsvolu("QBA2","TUBE", idtmed[kBe], ptube, 3); | |
632 | ||
633 | gMC->Gspos("QBA2", 1, "QB02", 0., 0., 0, 0, "ONLY"); | |
634 | // z = -2.75 -> -2.00 | |
635 | gMC->Gspos("QB02", 1, "QBT2", 0., 0.,-2.75+ptube[2], 0, "ONLY"); | |
636 | ||
637 | // Alu part | |
638 | ptube[0] = 2.900; | |
639 | ptube[1] = 3.150; | |
640 | ptube[2] = 2.375; | |
641 | // z = -2.00 -> 2.75 | |
642 | gMC->Gsvolu("QB04","TUBE", idtmed[kAlu], ptube, 3); | |
643 | gMC->Gspos("QB04", 1, "QBT2", 0., 0.,-2.+ptube[2], 0, "ONLY"); | |
644 | ||
645 | ||
646 | ptube[0] = 3.15; | |
647 | ptube[1] = 3.50; | |
648 | ptube[2] = 0.10; | |
649 | // z = 2.55 -> 2.75 | |
650 | gMC->Gsvolu("QB06","TUBE", idtmed[kAlu], ptube, 3); | |
651 | gMC->Gspos("QB06", 1, "QBT2", 0., 0., 2.55+ptube[2], 0, "ONLY"); | |
652 | ||
653 | ||
654 | // Fixation | |
655 | ptube[0] = 0.0; | |
656 | ptube[1] = 0.1; | |
657 | ptube[2] = 0.5; | |
658 | ||
659 | gMC->Gsvolu("QBA8","TUBE", idtmed[kInox], ptube, 3); | |
660 | gMC->Gspos("QBA8", 1 ,"QBT2", 0.000, 3.650, -1.25, idrotm[2002], "ONLY"); | |
661 | gMC->Gspos("QBA8", 2 ,"QBT2", 3.161, -1.825, -1.25, idrotm[2001], "ONLY"); | |
662 | gMC->Gspos("QBA8", 3 ,"QBT2", -3.161, -1.825, -1.25, idrotm[2003], "ONLY"); | |
663 | ||
664 | // Carbon ring | |
665 | ptube[0] = 3.15; | |
666 | ptube[1] = 4.10; | |
667 | ptube[2] = 0.55; | |
668 | ||
669 | gMC->Gsvolu("QB77","TUBE", idtmed[kC], ptube, 3); | |
670 | ||
671 | ptube[0] = 3.15; | |
672 | ptube[1] = 3.50; | |
673 | ptube[2] = 0.10; | |
674 | gMC->Gsvolu("QBB7","TUBE", idtmed[kInox], ptube, 3); | |
675 | gMC->Gspos("QBB7", 1, "QB77", 0.0, 0.0, 0.55-0.2, 0, "ONLY"); | |
676 | gMC->Gspos("QB77", 1, "QBT2", 0.0, 0.0, 2., 0, "ONLY"); | |
614dff1c | 677 | */ |
678 | // | |
679 | // Beam pipes between elements | |
680 | // | |
cd9c4f3a | 681 | |
614dff1c | 682 | gMC->Gsvolu("QB24","TUBE", idtmed[kInox], ptube, 0); |
683 | ptube[0] = RinSt; | |
684 | ptube[1] = RoutSt; | |
685 | ptube[2] = hlenQb24[0]; | |
686 | dz = hlenQbbe2 + ptube[2]; | |
687 | gMC->Gsposp("QB24", 1 ,"QBPM", 0.0, 0.0, -dz, 0, "ONLY", ptube, 3); | |
688 | // | |
689 | // Bellow on absorber side | |
690 | dz = dz+hlenQb24[0] + dzb; | |
691 | gMC->Gspos("QBE0", 1 ,"QBPM", 0.0, 0.0, -dz, 0, "ONLY"); | |
692 | // | |
693 | ptube[2] = hlenQb24[1]; | |
694 | dz = dz + dzb + ptube[2]; | |
695 | gMC->Gsposp("QB24", 2 ,"QBPM", 0.0, 0.0, -dz, 0, "ONLY", ptube, 3); | |
696 | dz = dz + ptube[2]; | |
7bd7faf8 | 697 | |
7bd7faf8 | 698 | // |
7bd7faf8 | 699 | // Flange |
700 | // | |
701 | // Mother Volume | |
614dff1c | 702 | ptube[0] = RinSt; |
7bd7faf8 | 703 | ptube[1] = 4.300; |
614dff1c | 704 | ptube[2] = 1.4; |
7bd7faf8 | 705 | |
614dff1c | 706 | gMC->Gsvolu("QFA0","TUBE", idtmed[kInox], ptube, 3); |
707 | dz = dz + ptube[2]; | |
708 | gMC->Gspos("QFA0", 1 ,"QBPM", 0.0, 0.0, -dz, 0, "ONLY"); | |
709 | dz = dz + ptube[2]; | |
7bd7faf8 | 710 | // |
614dff1c | 711 | // |
712 | ptube[0] = RinSt; | |
713 | ptube[1] = RoutSt; | |
714 | ptube[2] = hlenQb24[2]; | |
715 | dz = dz + ptube[2]; | |
716 | gMC->Gsposp("QB24", 3 ,"QBPM", 0.0, 0.0, -dz, 0, "ONLY", ptube, 3); | |
cd9c4f3a | 717 | |
718 | ||
fe4da5cc | 719 | // --- Place the PIPE ghost volume (QBPM) in its mother volume (ALIC) |
cd9c4f3a | 720 | // by rotating it to 180 deg. and make it invisible |
fe4da5cc | 721 | // |
7bd7faf8 | 722 | gMC->Gspos("QBPM",1,"ALIC",0,0,0,idrotm[2013], "ONLY"); |
fddcf8f7 | 723 | gMC->Gsbool("QBPM", "L3DX"); |
724 | gMC->Gsbool("QBPM", "L3O3"); | |
725 | gMC->Gsbool("QBPM", "L3O4"); | |
7bd7faf8 | 726 | |
cd9c4f3a | 727 | |
7bd7faf8 | 728 | // |
cd9c4f3a | 729 | // ******** Ion Pump volume description starts here ****** |
7bd7faf8 | 730 | // |
7bd7faf8 | 731 | // |
732 | // Getters -> | |
733 | pbox[0] = 6.50; | |
734 | pbox[1] = 6.75; | |
735 | pbox[2] = 15.60; | |
736 | gMC->Gsvolu("QI32","BOX", idtmed[kInox], pbox, 3); | |
737 | ||
738 | pbox[0] = 5.90; | |
739 | pbox[1] = 6.15; | |
740 | pbox[2] = 15.00; | |
741 | gMC->Gsvolu("QI42","BOX", idtmed[kGetter], pbox, 3); | |
742 | gMC->Gspos("QI42", 1, "QI32", 0.0, 0.0, 0.0, 0, "ONLY"); | |
743 | // <- | |
744 | ||
745 | ptube[0] = 0.0; | |
746 | ptube[1] = 19.0; | |
747 | ptube[2] = 2.5; | |
748 | gMC->Gsvolu("QI33","TUBE", idtmed[kInox], ptube, 3); | |
749 | ||
750 | ||
751 | ptube[0] = 0.0; | |
752 | ptube[1] = 15.0; | |
753 | ptube[2] = 2.5; | |
14b75006 | 754 | gMC->Gsvolu("QI43","TUBE", idtmed[kAir], ptube, 3); |
7bd7faf8 | 755 | gMC->Gspos("QI43", 1, "QI33", 0.0, 0.0, 0.0, 0, "ONLY"); |
756 | // | |
757 | // Connecting tube -> | |
758 | ptube[0] = 0.0; | |
759 | ptube[1] = 5.4; | |
cd9c4f3a | 760 | ptube[2] = 14.6; |
7bd7faf8 | 761 | gMC->Gsvolu("QI34","TUBE", idtmed[kInox], ptube, 3); |
762 | ||
763 | ptube[0] = 0.0; | |
764 | ptube[1] = 4.8; | |
cd9c4f3a | 765 | ptube[2] = 14.6; |
14b75006 | 766 | gMC->Gsvolu("QI44","TUBE", idtmed[kAir], ptube, 3); |
7bd7faf8 | 767 | gMC->Gspos("QI44", 1, "QI34", 0.0, 0.0, 0.0, 0, "ONLY"); |
768 | // <- | |
769 | ||
770 | // | |
771 | // Flange -> | |
cd9c4f3a | 772 | ptube[0] = 5.41; |
7bd7faf8 | 773 | ptube[1] = 7.30; |
774 | ptube[2] = 2.15; | |
775 | gMC->Gsvolu("QI35","TUBE", idtmed[kInox], ptube, 3); | |
7bd7faf8 | 776 | // <- |
cd9c4f3a | 777 | gMC->Gspos("QI32", 1, "QBPM", 0.0, -44.25, zPump, 0, "ONLY"); |
778 | gMC->Gspos("QI33", 1, "QBPM", 0.0, -35.00, zPump,idrotm[2002], "ONLY"); | |
779 | gMC->Gspos("QI34", 1, "QBPM", 0.0, -17.90, zPump,idrotm[2002], "ONLY"); | |
780 | gMC->Gspos("QI35", 1, "QBPM", 0.0, -24.35, zPump,idrotm[2002], "ONLY"); | |
fe4da5cc | 781 | |
22811e34 | 782 | gMC->Gsatt("QBPM", "SEEN", 1); |
783 | gMC->Gsatt("QBEM", "SEEN", 1); | |
fe4da5cc | 784 | } |
785 | ||
cd9c4f3a | 786 | |
fe4da5cc | 787 | |
788 | //___________________________________________ | |
789 | void AliPIPEv0::CreateMaterials() | |
790 | { | |
7bd7faf8 | 791 | // |
13391c3e | 792 | // Define materials for beam pipe |
7bd7faf8 | 793 | // |
13391c3e | 794 | |
9e1a0ddb | 795 | if(fDebug) printf("%s: Create PIPEv0 materials \n",ClassName()); |
7bd7faf8 | 796 | Int_t isxfld = gAlice->Field()->Integ(); |
797 | Float_t sxmgmx = gAlice->Field()->Max(); | |
2196bc42 | 798 | // Steel (Inox) |
7bd7faf8 | 799 | Float_t asteel[4] = { 55.847,51.9961,58.6934,28.0855 }; |
800 | Float_t zsteel[4] = { 26.,24.,28.,14. }; | |
801 | Float_t wsteel[4] = { .715,.18,.1,.005 }; | |
2196bc42 | 802 | // AlBe - alloy |
803 | Float_t aAlBe[2] = { 26.98, 9.01}; | |
804 | Float_t zAlBe[2] = { 13.00, 4.00}; | |
805 | Float_t wAlBe[2] = { 0.4, 0.6}; | |
806 | ||
7bd7faf8 | 807 | // |
808 | // Berillium | |
809 | AliMaterial(5, "BERILLIUM$", 9.01, 4., 1.848, 35.3, 36.7); | |
810 | // | |
811 | // Carbon | |
812 | AliMaterial(6, "CARBON$ ", 12.01, 6., 2.265, 18.8, 49.9); | |
813 | // | |
814 | // Aluminum | |
815 | AliMaterial(9, "ALUMINIUM$", 26.98, 13., 2.7, 8.9, 37.2); | |
816 | // | |
817 | // Air | |
818 | AliMaterial(15, "AIR$ ", 14.61, 7.3, .001205, 30423.24, 67500.); | |
819 | // | |
820 | // Vacuum | |
821 | AliMaterial(16, "VACUUM$ ", 1e-16, 1e-16, 1e-16, 1e16, 1e16); | |
822 | // | |
823 | // stainless Steel | |
824 | AliMixture(19, "STAINLESS STEEL$", asteel, zsteel, 7.88, 4, wsteel); | |
825 | // | |
826 | // reduced density steel to approximate pump getter material | |
827 | AliMixture(20, "GETTER$", asteel, zsteel, 1.00, 4, wsteel); | |
2196bc42 | 828 | // Al-Be alloy |
829 | // | |
830 | AliMixture(21, "AlBe$", aAlBe, zAlBe, 2.07, 2, wAlBe); | |
7bd7faf8 | 831 | // |
832 | // **************** | |
833 | // Defines tracking media parameters. | |
834 | // | |
835 | Float_t epsil = .001; // Tracking precision, | |
836 | Float_t stemax = -0.01; // Maximum displacement for multiple scat | |
837 | Float_t tmaxfd = -20.; // Maximum angle due to field deflection | |
838 | Float_t deemax = -.3; // Maximum fractional energy loss, DLS | |
839 | Float_t stmin = -.8; | |
840 | // *************** | |
841 | // | |
842 | // Beryllium | |
843 | ||
844 | AliMedium(5, "BE", 5, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin); | |
845 | ||
846 | // Carbon | |
847 | AliMedium(6, "C", 6, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin); | |
848 | // | |
849 | // Aluminum | |
850 | AliMedium(9, "ALU", 9, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin); | |
851 | // | |
852 | // Air | |
853 | AliMedium(15, "AIR", 15, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin); | |
854 | // | |
855 | // Vacuum | |
856 | AliMedium(16, "VACUUM", 16, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin); | |
857 | // | |
858 | // Steel | |
859 | AliMedium(19, "INOX", 19, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin); | |
860 | // | |
861 | // Getter | |
862 | AliMedium(20, "GETTER", 20, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin); | |
2196bc42 | 863 | // |
864 | // AlBe - Aloy | |
865 | AliMedium(21, "AlBe" , 21, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin); | |
866 | ||
fe4da5cc | 867 | } |
868 | ||
869 | ||
fe4da5cc | 870 | |
871 | ||
872 | ||
873 | ||
874 | ||
875 | ||
876 | ||
877 |