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