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