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.19  2002/03/25 12:34:40  morsch
  19 Obsolete support rollers removed.
  20
  21 Revision 1.18  2002/02/01 18:02:41  morsch
  22 Material of beam pipe between inner Be piece and forward detectors
  23 can be set by SetPipeMaterial(mat), mat = kInox, kAlu, kBe ...
  24
  25 Revision 1.17  2001/09/24 13:11:50  morsch
  26 Ion pump and bellows moved out by 15 cm to make space for forward
  27 detectors.
  28
  29 Revision 1.16  2001/05/16 14:57:22  alibrary
  30 New files for folders and Stack
  31
  32 Revision 1.15  2001/05/02 11:50:18  morsch
  33 New layout of the non-absorber side provided by Y. Viyogi. Not the final design
  34 but the prsent most realistic.
  35
  36 Revision 1.14  2001/01/20 16:56:33  morsch
  37 Put air in connecting tubes and flanges of vacuum pump.
  38
  39 Revision 1.13  2001/01/20 16:35:27  morsch
  40 Increase mother volume for bellows.
  41
  42 Revision 1.12  2000/12/21 16:41:06  morsch
  43 Coding convention clean-up (RS3)
  44
  45 Revision 1.11  2000/11/28 16:06:57  morsch
  46 Undulated beam-pipe replaced by Al-Be (40,60) pipe 1.5 mm thick.
  47
  48 Revision 1.10  2000/11/24 13:00:37  morsch
  49 - Geometry and materials imported from euclid output
  50 - include comments
  51 - better struturing of volume tree
  52 - improved version of flange close to front absorber
  53 - more realistic pump materials
  54 - undulated beam pipe imported from v3.
  55
  56 Revision 1.9  2000/10/02 21:28:15  fca
  57 Removal of useless dependecies via forward declarations
  58
  59 Revision 1.8  2000/06/11 12:37:01  morsch
  60 Coding rule violations corrected
  61
  62 Revision 1.7  2000/02/23 16:25:24  fca
  63 AliVMC and AliGeant3 classes introduced
  64 ReadEuclid moved from AliRun to AliModule
  65
  66 Revision 1.6  1999/09/29 09:24:30  fca
  67 Introduction of the Copyright and cvs Log
  68
  69 */
  70
  71 ////////////////////////////////////////////////
  72 //  Beam pipe class                            /
  73 ////////////////////////////////////////////////
  74
  75 #include "AliPIPEv0.h"
  76 #include "AliRun.h"
  77 #include "AliConst.h"
  78 #include "AliMagF.h"
  79 #include "AliMC.h"
  80 #include "TSystem.h"
  81
  82 #include <iostream.h>
  83
  84 ClassImp(AliPIPEv0)
  85
  86 //_____________________________________________________________________________
  87 AliPIPEv0::AliPIPEv0()
  88 {
  89 // Constructor
  90     SetPipeMaterial();
  91 }
  92
  93 //_____________________________________________________________________________
  94 AliPIPEv0::AliPIPEv0(const char *name, const char *title)
  95   : AliPIPE(name,title)
  96 {
  97 // Constructor
  98     SetPipeMaterial();
  99 }
 100
 101
 102 //___________________________________________
 103 void AliPIPEv0::CreateGeometry()
 104 {
 105 //Begin_Html
 106 /*
 107 <img src="picts/pipe.gif">
 108 */
 109 //End_Html
 110
 111
 112 //Begin_Html
 113 /*
 114 <img src="picts/tree_pipe.gif">
 115 */
 116 //End_Html
 117
 118     if(fDebug) printf("%s: Create PIPEv0 geometry \n",ClassName());
 119
 120
 121     Int_t *idtmed = fIdtmed->GetArray();
 122     Float_t ppcon[84], ptube[3], pbox[3];
 123     Int_t i=0;
 124
 125
 126     Int_t   idrotm[2099];
 127     AliMatrix(idrotm[2001],90.,240.,  0.,  0., 90.,150.);
 128     AliMatrix(idrotm[2002],90.,  0.,  0.,  0., 90.,270.);
 129     AliMatrix(idrotm[2003],90.,120.,  0.,  0., 90., 30.);
 130     AliMatrix(idrotm[2004],90.,315., 90., 45.,  0.,  0.);
 131     AliMatrix(idrotm[2005],90.,270., 90.,  0.,  0.,  0.);
 132     AliMatrix(idrotm[2006],90.,225., 90.,315.,  0.,  0.);
 133     AliMatrix(idrotm[2007],90.,180., 90.,270.,  0.,  0.);
 134     AliMatrix(idrotm[2008],90.,135., 90.,225.,  0.,  0.);
 135     AliMatrix(idrotm[2009],90., 90., 90.,180.,  0.,  0.);
 136     AliMatrix(idrotm[2010],90., 45., 90.,135.,  0.,  0.);
 137     idrotm[2011] = 0;
 138     AliMatrix(idrotm[2012],90.,180., 90., 90.,180.,  0.);
 139     AliMatrix(idrotm[2013],90.,  0., 90., 90.,180.,  0.);
 140 //
 141 //  Bellow
 142 //
 143 //  distance between bellows
 144     const Float_t dzbb =  9.0;
 145 //  total size of bellow section
 146     const Float_t dzb  = 15.0;
 147 //  size of undulated region
 148     const Float_t dzub =  2.0;
 149
 150 // half-lengths of various beam pipe sections
 151 // central Be-Pipe
 152     Float_t hlenQbbe1 = 40.;
 153     Float_t hlenQbbe2 = 36.5;
 154     Float_t hlenQbbe  = (hlenQbbe1+hlenQbbe2)/2.;
 155 //
 156 //
 157 //    Float_t hlenQbt1 = 5.5/2.;
 158 //
 159 //  Pipe outside central region (non-absober side)
 160     Float_t hlenQbab = 157.5;
 161 //
 162 //  Flange non-absorber side
 163     Float_t hlenQb29 = 11.5/2.+1.75 + 5.0;
 164 //
 165 //  Bellow element
 166     Float_t hlenQbe0 = dzb;
 167 //
 168 //  Inox pipe between Be and Bellow (absorber side)
 169     Float_t hlenQb24[3] = {10.5/2., 1.8, 3.3};
 170 //
 171 //
 172     Float_t hlenQb28 = (800.-hlenQbbe1-2.*hlenQbab-4.*hlenQb29-2.*hlenQbe0)/2.;
 173 //
 174 //  Position of the pump
 175     Float_t zPump = hlenQbbe1+2.*hlenQbab+2.*hlenQb29+dzb;
 176 //
 177 //  Inner beam pipe radius
 178     Float_t RinBe = 2.9;
 179     Float_t RinSt = 2.92;
 180 //
 181 //
 182     Float_t RoutBe = 2.98;
 183     Float_t RoutSt = 3.00;
 184
 185
 186 //
 187     Float_t dz;
 188
 189 //
 190 // The peam pipe up to the Front Absorber
 191 //
 192 // Mother Volume QBPM
 193     ppcon[0]  =   0;
 194     ppcon[1]  = 360;
 195     ppcon[2]  =  18;
 196 //  1:
 197     ppcon[3]  = -90.;
 198     ppcon[4]  =   0.;
 199     ppcon[5]  =   4.4;
 200 //  2
 201     ppcon[6]  = -90+2.*hlenQb24[2]+2.8+2.*hlenQb24[1];
 202     ppcon[7]  =   0.;
 203     ppcon[8]  =   4.4;
 204 //  3
 205     ppcon[9]  = ppcon[6];
 206     ppcon[10] =   0.;
 207     ppcon[11] =   4.1;
 208 //  4
 209     ppcon[12] = ppcon[9] + 2. + 2.*dzub+0.2;
 210     ppcon[13] =   0.;
 211     ppcon[14] =   4.1;
 212 //  5
 213     ppcon[15] = ppcon[12];
 214     ppcon[16] =   0.;
 215     ppcon[17] =   3.2;
 216 //  6
 217     ppcon[18] = ppcon[15] + 2.* dzbb-0.4;
 218     ppcon[19] =   0.;
 219     ppcon[20] =   3.2;
 220 //  7
 221     ppcon[21] = ppcon[18];
 222     ppcon[22] =   0.;
 223     ppcon[23] =   4.1;
 224 //  8
 225     ppcon[24] = -44.;
 226     ppcon[25] =   0.;
 227     ppcon[26] =   4.1;
 228 //  9
 229     ppcon[27] = -44.;
 230     ppcon[28] =    0;
 231     ppcon[29] =    3.0;
 232 //  10
 233     ppcon[30] =  38.;
 234     ppcon[31] =    0;
 235     ppcon[32] =    3.0;
 236 //  11
 237     ppcon[33] =  38.;
 238     ppcon[34] =    0;
 239     ppcon[35] =    3.6;
 240 //  12
 241     ppcon[36] = hlenQbbe1+2.*hlenQbab-0.1;
 242     ppcon[37] =    0.;
 243     ppcon[38] =    3.6;
 244 //  13
 245     ppcon[39] = ppcon[36];
 246     ppcon[40] =    0.;
 247     ppcon[41] =    3.6;
 248 //  14
 249     ppcon[42] = ppcon[39]+2.*hlenQb29-5.;
 250     ppcon[43] =    0.;
 251     ppcon[44] =    3.6;
 252 //  15
 253     ppcon[45] = ppcon[42];
 254     ppcon[46] =    0.;
 255     ppcon[47] =   56.;
 256 //  16
 257     ppcon[48] = ppcon[45]+2.*dzb+10.;
 258     ppcon[49] =    0.;
 259     ppcon[50] =   56.;
 260 //  17
 261     ppcon[51] =   ppcon[48];
 262     ppcon[52] =    0.;
 263     ppcon[53] =    3.6;
 264 //  18
 265     ppcon[54] =  800.;
 266     ppcon[55] =    0.;
 267     ppcon[56] =    3.6;
 268
 269     gMC->Gsvolu("QBPM", "PCON", idtmed[kAir], ppcon,57);
 270
 271
 272 //
 273 // volume definitions of various sections
 274 //
 275
 276 //
 277 // The Vacuum
 278     gMC->Gsvolu("QBVA","TUBE", idtmed[kVac], ptube, 0);
 279     ptube[0] =   0.0;
 280     ptube[1] =   RinSt;
 281     ptube[2] =   (90.-hlenQbbe2)/2.;
 282     dz = -90. + ptube[2];
 283     gMC->Gsposp ("QBVA", 1, "QBPM", 0., 0., dz , 0, "ONLY", ptube, 3);
 284     dz = dz + ptube[2];
 285
 286     ptube[1] =   RinBe;
 287     ptube[2] =   hlenQbbe+hlenQbab;
 288     dz = dz + ptube[2];
 289     gMC->Gsposp ("QBVA", 2, "QBPM", 0., 0., dz , 0, "ONLY", ptube, 3);
 290     dz = dz + ptube[2];
 291
 292     ptube[1] =   RinSt;
 293     ptube[2] =   (800.-hlenQbbe1-2.*hlenQbab)/2.;
 294     dz = dz + ptube[2];
 295     gMC->Gsposp ("QBVA", 3, "QBPM", 0., 0., dz , 0, "ONLY", ptube, 3);
 296
 297 //
 298 // Be Pipe in central Alice
 299     ptube[0] = RinBe;
 300     ptube[1] = RoutBe;
 301     ptube[2] = hlenQbbe;
 302
 303     gMC->Gsvolu("QBBE","TUBE", idtmed[kBe], ptube, 3);
 304
 305 //
 306 //  Support Ring
 307 //
 308     //  Mother
 309     ppcon[0]  =   0;
 310     ppcon[1]  = 360;
 311     ppcon[2]  =   6;
 312 //  1:
 313     ppcon[3]  = -1.;
 314     ppcon[4]  = 3.0;
 315     ppcon[5]  = 3.4;
 316 //  2
 317     ppcon[6]  = -0.8;
 318     ppcon[7]  = 3.0;
 319     ppcon[8]  = 3.4;
 320 //  3
 321     ppcon[9]  = -0.8;
 322     ppcon[10] = 3.0;
 323     ppcon[11] = 3.2;
 324 //  4
 325     ppcon[12] = 0.8;
 326     ppcon[13] = 3.0;
 327     ppcon[14] = 3.2;
 328 //  5
 329     ppcon[15] = 0.8;
 330     ppcon[16] = 3.0;
 331     ppcon[17] = 3.4;
 332 //  6
 333     ppcon[18] = 1.0;
 334     ppcon[19] = 3.0;
 335     ppcon[20] = 3.4;
 336
 337
 338     gMC->Gsvolu("QBSR", "PCON", idtmed[kC], ppcon,21);
 339     gMC->Gspos("QBSR", 1, "QBPM", 0.0, 0.0,  40., 0, "ONLY");
 340     gMC->Gspos("QBSR", 2, "QBPM", 0.0, 0.0, 150., 0, "ONLY");
 341     gMC->Gspos("QBSR", 3, "QBPM", 0.0, 0.0, 260., 0, "ONLY");
 342     gMC->Gspos("QBSR", 4, "QBPM", 0.0, 0.0,- 46., 0, "ONLY");
 343 //
 344 // Flange and Fixed Point: non absorber side
 345 //
 346 // ---------->
 347 //
 348 //  Mother
 349     ppcon[0]  =   0;
 350     ppcon[1]  = 360;
 351     ppcon[2]  =   4;
 352 //  1:
 353     ppcon[3]  = -hlenQb29;
 354     ppcon[4]  = RinSt;
 355     ppcon[5]  = 5.8;
 356 //  2
 357     ppcon[6]  = ppcon[3]+3.6;
 358     ppcon[7]  = RinSt;
 359     ppcon[8]  = 5.8;
 360 //  3
 361     ppcon[9]  = ppcon[6];
 362     ppcon[10] = RinSt;
 363     ppcon[11] = 3.6;
 364 //  4
 365     ppcon[12] = hlenQb29;
 366     ppcon[13] = RinSt;
 367     ppcon[14] = 3.6;
 368
 369     gMC->Gsvolu("QB29", "PCON", idtmed[kAir], ppcon,15);
 370
 371
 372 //    Flange
 373     ptube[0] = RinSt;
 374     ptube[1] = 5.7;
 375     ptube[2] = 1.75;
 376     gMC->Gsvolu("QF29","TUBE", idtmed[kInox], ptube, 3);
 377     gMC->Gspos("QF29", 1, "QB29", 0.0, 0.0, -hlenQb29+1.75, 0, "ONLY");
 378 //    Pipe
 379     ptube[0] = RinSt;
 380     ptube[1] = 3.0;
 381     ptube[2] = hlenQb29;
 382     gMC->Gsvolu("QS29","TUBE", idtmed[kInox], ptube, 3);
 383     gMC->Gspos("QS29", 1, "QB29", 0.0, 0.0, 0., 0, "ONLY");
 384 //    Fixed point
 385     ptube[0] = RinSt;
 386     ptube[1] = 3.5;
 387     ptube[2] = 0.3;
 388     gMC->Gsvolu("QP29","TUBE", idtmed[kInox], ptube, 3);
 389     gMC->Gspos("QP29", 1, "QB29", 0.0, 0.0, -hlenQb29+9.75+3., 0, "ONLY");
 390
 391 //
 392 //
 393 // Inox beam pipe: final section on non-absorber side
 394
 395     ptube[0] =   RinSt;
 396     ptube[1] =   RoutSt;
 397     ptube[2] =   hlenQb28;
 398
 399     gMC->Gsvolu("QB28","TUBE", idtmed[kInox], ptube, 3);
 400
 401
 402 //  Al-Be (40-60 wgt%, rho=2.7 g/cm**3) beam pipe
 403 //
 404 //  This section is under study (A.M. 1/2/2002)
 405 //
 406
 407     ptube[0] = RinBe;
 408     if (fPipeMaterial == kAlu) {
 409         ptube[1] = 3.0;
 410     } else if (fPipeMaterial == kBe) {
 411         ptube[1] = RoutBe;
 412     } else if (fPipeMaterial == kInox){
 413         ptube[1] = RoutSt;
 414     }
 415     ptube[2] =   hlenQbab;
 416
 417     gMC->Gsvolu("QBAB","TUBE", idtmed[fPipeMaterial], ptube, 3);
 418
 419 // 2.5 mm thick SS tube for hanging pump
 420 /*
 421     ptube[0] = Rin;
 422     ptube[1] = 3.15;
 423     ptube[2] = hlenQb26;
 424
 425     gMC->Gsvolu("QB26","TUBE", idtmed[kInox], ptube, 3);
 426 */
 427 //
 428 // Bellows
 429 //
 430 //
 431 // Mother Volume
 432     Float_t pconQBE0[33];
 433     pconQBE0[ 0]= 0;
 434     pconQBE0[ 1]= 360;
 435     pconQBE0[ 2]= 10;
 436 //  1
 437     pconQBE0[ 3]= -dzb;
 438     pconQBE0[ 4]= RinSt;
 439     pconQBE0[ 5]= RoutSt;
 440 //  2
 441     pconQBE0[ 6]= -dzb+2.;
 442     pconQBE0[ 7]= RinSt;
 443     pconQBE0[ 8]= RoutSt;
 444 //  3
 445     pconQBE0[ 9]= -dzb+2.;
 446     pconQBE0[10]= RinSt;
 447     pconQBE0[11]= 4.00;
 448 //  4
 449     pconQBE0[12]= -dzb+2.+2.*dzub;
 450     pconQBE0[13]= RinSt;
 451     pconQBE0[14]= 4.00;
 452 //  5
 453     pconQBE0[15]= -dzb+2.+2.*dzub;
 454     pconQBE0[16]= RinSt;
 455     pconQBE0[17]= RoutSt;
 456 //  6
 457     pconQBE0[18]= -dzb+2.+2.*dzub+2.*dzbb;
 458     pconQBE0[19]= RinSt;
 459     pconQBE0[20]= RoutSt;
 460 //  7
 461     pconQBE0[21]= -dzb+2.+2.*dzub+2.*dzbb;
 462     pconQBE0[22]= RinSt;
 463     pconQBE0[23]= 4.00;
 464 //  8
 465     pconQBE0[24]= -dzb+2.+4.*dzub+2.*dzbb;
 466     pconQBE0[25]= RinSt;
 467     pconQBE0[26]= 4.00;
 468 //  9
 469     pconQBE0[27]= -dzb+2.+4.*dzub+2.*dzbb;
 470     pconQBE0[28]= RinSt;
 471     pconQBE0[29]= RoutSt;
 472 //  10
 473     pconQBE0[30]= +dzb;
 474     pconQBE0[31]= RinSt;
 475     pconQBE0[32]= RoutSt;
 476     gMC->Gsvolu("QBE0", "PCON", idtmed[kAir], pconQBE0, 33);
 477 //
 478 //  Undulated piece mother
 479     ptube[0] =  RinSt;
 480     ptube[1] =  4.00;
 481     ptube[2] =  dzub;
 482     gMC->Gsvolu("QBEM","TUBE", idtmed[kAir], ptube, 3);
 483     dz = -dzb+2.+dzub;
 484     gMC->Gspos("QBEM", 2 ,"QBE0", 0.0, 0.0,   dz, 0 , "ONLY");
 485     gMC->Gspos("QBEM", 1 ,"QBE0", 0.0, 0.0,  -dz, idrotm[2012], "ONLY");
 486 //
 487     Float_t pund[30];
 488     Float_t uw = 0.02;
 489     Float_t ur = 2.*dzub-12.*uw;
 490     Float_t uz = ur/13.;
 491     Float_t ut = uz+uw;
 492
 493     pund[ 0] = 0;
 494     pund[ 1] = 360;
 495     pund[ 2] =  9;
 496
 497     pund[ 3] = -ut;
 498     pund[ 4] = RinSt;
 499     pund[ 5] = RinSt+uw;
 500
 501     pund[ 6] = -ut+uz;
 502     pund[ 7] = pund[4];
 503     pund[ 8] = pund[5];
 504
 505     pund[ 9] = pund[6];
 506     pund[10] = pund[4];
 507     pund[11] = 4.0;
 508
 509     pund[12] = pund[9]+uw;
 510     pund[13] = pund[10];
 511     pund[14] = pund[11];
 512
 513     pund[15] = pund[12];
 514     pund[16] = 4.0-uw;
 515     pund[17] = pund[14];
 516
 517     pund[18] = pund[12]+uz;
 518     pund[19] = pund[16];
 519     pund[20] = pund[17];
 520
 521     pund[21] = pund[18];
 522     pund[22] = RinSt;
 523     pund[23] = pund[20];
 524
 525     pund[24] = pund[21]+uw;
 526     pund[25] = pund[22];
 527     pund[26] = pund[23];
 528
 529     pund[27] = pund[24];
 530     pund[28] = pund[25];
 531     pund[29] = RinSt+uw;
 532
 533     gMC->Gsvolu("QBEU", "PCON", idtmed[kInox], pund, 30);
 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