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.25  2002/11/21 16:28:39  alibrary
  19 Removing AliMCProcess and AliMC
  20
  21 Revision 1.24  2002/10/22 15:02:47  alibrary
  22 Introducing Riostream.h
  23
  24 Revision 1.23  2002/10/14 14:57:39  hristov
  25 Merging the VirtualMC branch to the main development branch (HEAD)
  26
  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
  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
  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
  43 Revision 1.19  2002/03/25 12:34:40  morsch
  44 Obsolete support rollers removed.
  45
  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
  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
  54 Revision 1.16  2001/05/16 14:57:22  alibrary
  55 New files for folders and Stack
  56
  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
  61 Revision 1.14  2001/01/20 16:56:33  morsch
  62 Put air in connecting tubes and flanges of vacuum pump.
  63
  64 Revision 1.13  2001/01/20 16:35:27  morsch
  65 Increase mother volume for bellows.
  66
  67 Revision 1.12  2000/12/21 16:41:06  morsch
  68 Coding convention clean-up (RS3)
  69
  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
  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
  81 Revision 1.9  2000/10/02 21:28:15  fca
  82 Removal of useless dependecies via forward declarations
  83
  84 Revision 1.8  2000/06/11 12:37:01  morsch
  85 Coding rule violations corrected
  86
  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
  91 Revision 1.6  1999/09/29 09:24:30  fca
  92 Introduction of the Copyright and cvs Log
  93
  94 */
  95
  96 ////////////////////////////////////////////////
  97 //  Beam pipe class                            /
  98 ////////////////////////////////////////////////
  99
 100 #include "AliPIPEv0.h"
 101 #include "AliRun.h"
 102 #include "AliConst.h"
 103 #include "AliMagF.h"
 104 #include "TSystem.h"
 105
 106 #include <Riostream.h>
 107
 108 ClassImp(AliPIPEv0)
 109
 110 //_____________________________________________________________________________
 111 AliPIPEv0::AliPIPEv0()
 112 {
 113 // Constructor
 114     SetPipeMaterial();
 115 }
 116
 117 //_____________________________________________________________________________
 118 AliPIPEv0::AliPIPEv0(const char *name, const char *title)
 119   : AliPIPE(name,title)
 120 {
 121 // Constructor
 122     SetPipeMaterial();
 123 }
 124
 125
 126 //___________________________________________
 127 void AliPIPEv0::CreateGeometry()
 128 {
 129 //Begin_Html
 130 /*
 131 <img src="picts/pipe.gif">
 132 */
 133 //End_Html
 134
 135
 136 //Begin_Html
 137 /*
 138 <img src="picts/tree_pipe.gif">
 139 */
 140 //End_Html
 141
 142     if(fDebug) printf("%s: Create PIPEv0 geometry \n",ClassName());
 143
 144
 145     Int_t *idtmed = fIdtmed->GetArray();
 146     Float_t ppcon[90], ptube[3], pbox[3];
 147     Int_t i=0;
 148
 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.);
 164 //
 165 //  Bellow
 166 //
 167 //  distance between bellows
 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;
 208
 209
 210 //
 211     Float_t dz;
 212
 213 //
 214 // The peam pipe up to the Front Absorber
 215 //
 216 // Mother Volume QBPM
 217     ppcon[0]  =   0;
 218     ppcon[1]  = 360;
 219     ppcon[2]  =  20;
 220 //  1
 221     ppcon[3]  = -90.;
 222     ppcon[4]  =   0.;
 223     ppcon[5]  =   3.1;
 224 //  2
 225     ppcon[6]  = -84.;
 226     ppcon[7]  =   0.;
 227     ppcon[8]  =   3.1;
 228 //  3
 229     ppcon[9]  = -84.;
 230     ppcon[10] =   0.;
 231     ppcon[11] =   4.4;
 232 //  4
 233     ppcon[12] = -90+2.*hlenQb24[2]+2.8+2.*hlenQb24[1];
 234     ppcon[13] =   0.;
 235     ppcon[14] =   4.4;
 236 //  5
 237     ppcon[15]  = ppcon[12];
 238     ppcon[16] =   0.;
 239     ppcon[17] =   4.1;
 240 //  6
 241     ppcon[18] = ppcon[15] + 2. + 2.*dzub+0.2;
 242     ppcon[19] =   0.;
 243     ppcon[20] =   4.1;
 244 //  7
 245     ppcon[21] = ppcon[18];
 246     ppcon[22] =   0.;
 247     ppcon[23] =   3.2;
 248 //  8
 249     ppcon[24] = ppcon[21] + 2.* dzbb-0.4;
 250     ppcon[25] =   0.;
 251     ppcon[26] =   3.2;
 252 //  9
 253     ppcon[27] = ppcon[24];
 254     ppcon[28] =   0.;
 255     ppcon[29] =   4.1;
 256 //  10
 257     ppcon[30] = -44.;
 258     ppcon[31] =   0.;
 259     ppcon[32] =   4.1;
 260 //  11
 261     ppcon[33] = -44.;
 262     ppcon[34] =    0;
 263     ppcon[35] =    3.0;
 264 //  12
 265     ppcon[36] =  38.;
 266     ppcon[37] =    0;
 267     ppcon[38] =    3.0;
 268 //  13
 269     ppcon[39] =  38.;
 270     ppcon[40] =    0;
 271     ppcon[41] =    3.6;
 272 //  14
 273     ppcon[42] = hlenQbbe1+2.*hlenQbab-0.1;
 274     ppcon[43] =    0.;
 275     ppcon[44] =    3.6;
 276 //  15
 277     ppcon[45] = ppcon[42];
 278     ppcon[46] =    0.;
 279     ppcon[47] =    3.6;
 280 //  16
 281     ppcon[48] = ppcon[45]+2.*hlenQb29-5.;
 282     ppcon[49] =    0.;
 283     ppcon[50] =    3.6;
 284 //  17
 285     ppcon[51] = ppcon[48];
 286     ppcon[52] =    0.;
 287     ppcon[53] =   56.;
 288 //  18
 289     ppcon[54] = ppcon[51]+2.*dzb+10.;
 290     ppcon[55] =    0.;
 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;
 300
 301     gMC->Gsvolu("QBPM", "PCON", idtmed[kAir], ppcon,63);
 302
 303
 304 //
 305 // volume definitions of various sections
 306 //
 307
 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);
 328
 329 //
 330 // Be Pipe in central Alice
 331     ptube[0] = RinBe;
 332     ptube[1] = RoutBe;
 333     ptube[2] = hlenQbbe;
 334
 335     gMC->Gsvolu("QBBE","TUBE", idtmed[kBe], ptube, 3);
 336
 337 //
 338 //  Support Ring
 339 //
 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;
 368
 369
 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");
 375 //
 376 // Flange and Fixed Point: non absorber side
 377 //
 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;
 400
 401     gMC->Gsvolu("QB29", "PCON", idtmed[kAir], ppcon,15);
 402
 403
 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;
 413     ptube[2] = hlenQb29;
 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");
 422
 423 //
 424 //
 425 // Inox beam pipe: final section on non-absorber side
 426
 427     ptube[0] =   RinSt;
 428     ptube[1] =   RoutSt;
 429     ptube[2] =   hlenQb28;
 430
 431     gMC->Gsvolu("QB28","TUBE", idtmed[kInox], ptube, 3);
 432
 433
 434 //  Al-Be (40-60 wgt%, rho=2.7 g/cm**3) beam pipe
 435 //
 436 //  This section is under study (A.M. 1/2/2002)
 437 //
 438
 439     ptube[0] = RinBe;
 440     if (fPipeMaterial == kAlu) {
 441         ptube[1] = 3.0;
 442     } else if (fPipeMaterial == kBe) {
 443         ptube[1] = RoutBe;
 444     } else if (fPipeMaterial == kInox){
 445         ptube[1] = RoutSt;
 446     }
 447     ptube[2] =   hlenQbab;
 448
 449     gMC->Gsvolu("QBAB","TUBE", idtmed[fPipeMaterial], ptube, 3);
 450
 451 // 2.5 mm thick SS tube for hanging pump
 452 /*
 453     ptube[0] = Rin;
 454     ptube[1] = 3.15;
 455     ptube[2] = hlenQb26;
 456
 457     gMC->Gsvolu("QB26","TUBE", idtmed[kInox], ptube, 3);
 458 */
 459 //
 460 // Bellows
 461 //
 462 //
 463 // Mother Volume
 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;
 514     gMC->Gsvolu("QBEM","TUBE", idtmed[kAir], ptube, 3);
 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;
 527     pund[ 2] =  8;
 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];
 560
 561     gMC->Gsvolu("QBEU", "PCON", idtmed[kInox], pund, 27);
 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");
 567     }
 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");
 573 //
 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 ****
 590 //
 591     //
 592     // first the beryllium section
 593     Float_t zpos = -(hlenQbbe2-hlenQbbe1)/2;
 594     gMC->Gspos("QBBE", 1, "QBPM", 0., 0., zpos, 0, "ONLY");
 595
 596     // next meta-metal transition QBT1 on on-absorber side
 597 //    zpos = zpos + hlenQbbe + hlenQbt1;
 598 //    gMC->Gspos("QBT1", 1, "QBPM", 0., 0.,  zpos, 0, "ONLY");
 599
 600     // Aluminium OR Al-be alloy section
 601     zpos = hlenQbbe1+hlenQbab;
 602     gMC->Gspos("QBAB", 1, "QBPM", 0.0, 0.0, zpos, 0, "ONLY");
 603     //
 604     // inox flange at the start of bellow
 605     zpos = zpos + hlenQbab + hlenQb29;
 606     gMC->Gspos("QB29", 1, "QBPM", 0.0, 0.0, zpos, idrotm[2012], "ONLY");
 607     //
 608     // bellow section
 609     zpos = zpos + hlenQb29 + hlenQbe0;
 610     gMC->Gspos("QBE0", 2 ,"QBPM", 0.0, 0.0, zpos, 0, "ONLY");
 611     //
 612     // inox flange at the end of bellow and start of thick inox for pump
 613     zpos = zpos + hlenQbe0 + hlenQb29;
 614     gMC->Gspos("QB29", 2, "QBPM", 0.0, 0.0, zpos, 0, "ONLY");
 615     //
 616     //last inox section till 800 cm
 617     zpos = zpos + hlenQb29 + hlenQb28;
 618     gMC->Gspos("QB28", 1, "QBPM", 0.0, 0.0, zpos, 0, "ONLY");
 619
 620 //******** end of placement on no-absorber side *********
 621
 622     //
 623     // **** Absorber side *****
 624     //
 625     //
 626 /*
 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");
 690  */
 691 //
 692 // Beam pipes between elements
 693 //
 694
 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];
 710
 711 //
 712 // Flange
 713 //
 714 //  Mother Volume
 715     ptube[0] = RinSt;
 716     ptube[1] = 4.300;
 717     ptube[2] = 1.4;
 718
 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];
 723 //
 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);
 730
 731
 732 // --- Place the PIPE ghost volume (QBPM) in its mother volume (ALIC)
 733 //    by rotating it to 180 deg. and make it invisible
 734 //
 735     gMC->Gspos("QBPM",1,"ALIC",0,0,0,idrotm[2013], "ONLY");
 736     gMC->Gsbool("QBPM", "L3DX");
 737     gMC->Gsbool("QBPM", "L3O3");
 738     gMC->Gsbool("QBPM", "L3O4");
 739
 740
 741 //
 742 // ******** Ion Pump volume description starts here ******
 743 //
 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;
 767     gMC->Gsvolu("QI43","TUBE", idtmed[kAir], ptube, 3);
 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;
 773     ptube[2] = 14.6;
 774     gMC->Gsvolu("QI34","TUBE", idtmed[kInox], ptube, 3);
 775
 776     ptube[0] =  0.0;
 777     ptube[1] =  4.8;
 778     ptube[2] = 14.6;
 779     gMC->Gsvolu("QI44","TUBE", idtmed[kAir], ptube, 3);
 780     gMC->Gspos("QI44", 1, "QI34", 0.0, 0.0, 0.0, 0, "ONLY");
 781 // <-
 782
 783   //
 784   // Flange ->
 785     ptube[0] =  5.41;
 786     ptube[1] =  7.30;
 787     ptube[2] =  2.15;
 788     gMC->Gsvolu("QI35","TUBE", idtmed[kInox], ptube, 3);
 789 // <-
 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");
 794
 795     gMC->Gsatt("QBPM", "SEEN", 1);
 796     gMC->Gsatt("QBEM", "SEEN", 1);
 797 }
 798
 799
 800
 801 //___________________________________________
 802 void AliPIPEv0::CreateMaterials()
 803 {
 804   //
 805   // Define materials for beam pipe
 806   //
 807
 808   if(fDebug) printf("%s: Create PIPEv0 materials \n",ClassName());
 809   Int_t   isxfld = gAlice->Field()->Integ();
 810   Float_t sxmgmx = gAlice->Field()->Max();
 811   // Steel (Inox)
 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 };
 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
 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);
 841   //     Al-Be alloy
 842   //
 843   AliMixture(21, "AlBe$", aAlBe, zAlBe, 2.07, 2, wAlBe);
 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);
 876   //
 877   //   AlBe - Aloy
 878   AliMedium(21, "AlBe"  , 21, 0, isxfld, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
 879
 880 }
 881
 882
 883
 884
 885
 886
 887
 888
 889
 890