STRUCT/AliABSOv0.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.6  2000/06/15 09:40:31  morsch
  19 Obsolete typedef keyword removed
  20
  21 Revision 1.5  2000/06/12 19:39:01  morsch
  22 New structure of beam pipe and heating jacket.
  23
  24 Revision 1.4  2000/04/03 08:13:40  fca
  25 Introduce extra scope for non ANSI compliant C++ compilers
  26
  27 Revision 1.3  2000/01/18 17:49:56  morsch
  28 Serious overlap of ABSM with shield corrected
  29 Small error in ARPB parameters corrected
  30
  31 Revision 1.2  2000/01/13 11:23:59  morsch
  32 Last layer of Pb outer angle corrected
  33
  34 Revision 1.1  2000/01/12 15:39:30  morsch
  35 Standard version of ABSO
  36
  37 */
  38
  39 ///////////////////////////////////////////////////////////////////////////////
  40 //                                                                           //
  41 //  Muon ABSOrber                                                            //
  42 //  This class contains the description of the muon absorber geometry        //
  43 //                                                                           //
  44 //Begin_Html
  45 /*
  46 <img src="picts/AliABSOClass.gif">
  47 </pre>
  48 <br clear=left>
  49 <font size=+2 color=red>
  50 <p>The responsible person for this module is
  51 <a href="mailto:andreas.morsch@cern.ch">Andreas Morsch</a>.
  52 </font>
  53 <pre>
  54 */
  55 //End_Html
  56 //                                                                           //
  57 //                                                                           //
  58 ///////////////////////////////////////////////////////////////////////////////
  59
  60 #include "AliABSOv0.h"
  61 #include "AliRun.h"
  62 #include "AliMC.h"
  63 #include "AliConst.h"
  64
  65 ClassImp(AliABSOv0)
  66
  67 //_____________________________________________________________________________
  68 AliABSOv0::AliABSOv0()
  69 {
  70   //
  71   // Default constructor
  72   //
  73 }
  74
  75 //_____________________________________________________________________________
  76 AliABSOv0::AliABSOv0(const char *name, const char *title)
  77        : AliABSO(name,title)
  78 {
  79   //
  80   // Standard constructor
  81   //
  82   SetMarkerColor(7);
  83   SetMarkerStyle(2);
  84   SetMarkerSize(0.4);
  85 }
  86
  87 //_____________________________________________________________________________
  88 void AliABSOv0::CreateGeometry()
  89 {
  90     //
  91     // Creation of the geometry of the muon absorber
  92     //
  93     //Begin_Html
  94     /*
  95       <img src="picts/AliABSOv0Tree.gif">
  96     */
  97     //End_Html
  98     //Begin_Html
  99     /*
 100       <img src="picts/AliABSOv0.gif">
 101     */
 102     //End_Html
 103
 104     //
 105     //
 106
 107     enum {kC=1605, kAl=1608, kFe=1609, kCu=1610, kW=1611, kPb=1612,
 108                   kNiCuW=1620, kVacuum=1615, kAir=1614, kConcrete=1616,
 109                   kPolyCH2=1617, kSteel=1609, kInsulation=1613};
 110
 111     Int_t *idtmed = fIdtmed->GetArray()-1599;
 112
 113     Float_t par[24], cpar[5], cpar0[5], pcpar[12], tpar[3], tpar0[3];
 114     Float_t dz;
 115
 116 #include "ABSOSHILConst.h"
 117 #include "ABSOConst.h"
 118     Float_t dTube=0.1;                     // tube thickness
 119     Float_t dInsu=0.5;                     // insulation thickness
 120     Float_t dEnve=0.1;                     // protective envelope thickness
 121     Float_t dFree=0.5;                     // clearance thickness
 122
 123
 124 // Mother volume and outer shielding: Pb
 125   par[0]  = 0.;
 126   par[1]  = 360.;
 127   par[2]  = 7.;
 128
 129   par[3]  = -(zRear-zAbsStart)/2.;
 130   par[4]  = rAbs;
 131   par[5]  = zAbsStart * TMath::Tan(theta1);
 132
 133   par[6]  = par[3]+(zNose-zAbsStart);
 134   par[7]  = rAbs;
 135   par[8]  = zNose * TMath::Tan(theta1);
 136
 137   par[9]  = par[3]+(zConeTPC-zAbsStart);
 138   par[10] = rAbs;
 139   par[11] = par[8] + (par[9] - par[6]) * TMath::Tan(theta2);
 140
 141   par[12]  = par[3]+(zOpen-zAbsStart);
 142   par[13] = rAbs;
 143   par[14] = par[11] + (par[12] - par[9]) * TMath::Tan(accMax);
 144
 145   par[15] = par[3]+(zRear-dRear-zAbsStart);
 146   par[16] = rAbs   + (par[15] - par[12]) * TMath::Tan(thetaOpen1) ;
 147   par[17] = par[14] + (par[15] - par[12]) * TMath::Tan(accMax);
 148
 149   par[18] = par[3]+(zRear-dRear-zAbsStart);
 150   par[19] = (zRear-dRear) * TMath::Tan(accMin);
 151   par[20] = par[14] + (par[18] - par[12]) * TMath::Tan(accMax);
 152
 153   par[21] = -par[3];
 154   par[22] =  zRear* TMath::Tan(accMin);
 155   par[23] = par[20] + (par[21] - par[18]) * TMath::Tan(accMax);
 156   gMC->Gsvolu("ABSS", "PCON", idtmed[kPb], par, 24);
 157   { // Begin local scope for i
 158       for (Int_t i=4; i<18; i+=3) par[i]  = 0;
 159   } // End local scope for i
 160   gMC->Gsvolu("ABSM", "PCON", idtmed[kVacuum+40], par, 24);
 161   gMC->Gspos("ABSS", 1, "ABSM", 0., 0., 0., 0, "ONLY");
 162
 163 //
 164 // Steel envelope
 165 //
 166   par[4] = par[5] -dSteel;
 167   par[7] = par[8] -dSteel;
 168   par[10]= par[11]-dSteel;
 169   par[13]= par[14]-dSteel;
 170   par[16]= par[17]-dSteel;
 171   par[19]= par[20]-dSteel;
 172   par[22]= par[23]-dSteel;
 173   gMC->Gsvolu("ABST", "PCON", idtmed[kSteel], par, 24);
 174   gMC->Gspos("ABST", 1, "ABSS", 0., 0., 0., 0, "ONLY");
 175 //
 176 // Polyethylene shield
 177 //
 178   cpar[0] = (zRear - zConeTPC) / 2.;
 179   cpar[1] = zConeTPC * TMath::Tan(accMax);
 180   cpar[2] = cpar[1] + dPoly;
 181   cpar[3] = zRear * TMath::Tan(accMax);
 182   cpar[4] = cpar[3] + dPoly;
 183   gMC->Gsvolu("APOL", "CONE", idtmed[kPolyCH2+40], cpar, 5);
 184   dz = (zRear-zAbsStart)/2.-cpar[0];
 185   gMC->Gspos("APOL", 1, "ABSS", 0., 0., dz, 0, "ONLY");
 186
 187 //
 188 // Tungsten nose to protect TPC
 189 //
 190   cpar[0] = (zNose - zAbsStart) / 2.;
 191   cpar[1] = zAbsStart * TMath::Tan(accMax);
 192   cpar[2] = zAbsStart * TMath::Tan(theta1)-dSteel;
 193   cpar[3] = zNose * TMath::Tan(accMax);
 194   cpar[4] = zNose * TMath::Tan(theta1)-dSteel;
 195   gMC->Gsvolu("ANOS", "CONE", idtmed[kW], cpar, 5);
 196 //
 197   dz = -(zRear-zAbsStart)/2.+cpar[0];
 198   gMC->Gspos("ANOS", 1, "ABSS", 0., 0., dz, 0, "ONLY");
 199 //
 200 // Tungsten inner shield
 201 //
 202   Float_t zW=zTwoDeg+.1;
 203   Float_t dZ = zW+(zRear-dRear-zW)/2.;
 204   //
 205   pcpar[0]  = 0.;
 206   pcpar[1]  = 360.;
 207   pcpar[2]  = 3.;
 208   pcpar[3]  = zW-dZ;
 209   pcpar[4]  = rAbs;
 210   pcpar[5]  = zW * TMath::Tan(accMin);
 211   pcpar[6]  = zOpen-dZ;
 212   pcpar[7]  = rAbs;
 213   pcpar[8]  = zOpen * TMath::Tan(accMin);
 214   pcpar[9]  = zRear-dRear-dZ;
 215   pcpar[10] = rAbs+(zRear-dRear-zOpen) * TMath::Tan(thetaOpen1);
 216   pcpar[11] = (zRear-dRear) * TMath::Tan(accMin);
 217
 218   gMC->Gsvolu("AWIN", "PCON", idtmed[kNiCuW+40], pcpar, 12);
 219   //
 220   dz=(zW+zRear-dRear)/2-(zAbsStart+zRear)/2.;
 221   gMC->Gspos("AWIN", 1, "ABSS", 0., 0., dz, 0, "ONLY");
 222
 223   //     Inner tracking region
 224   //
 225   //     mother volume: Pb
 226   //
 227   pcpar[0]  = 0.;
 228   pcpar[1]  = 360.;
 229   pcpar[2]  = 3.;
 230   pcpar[3]  = -(zRear-zAbsStart)/2.;
 231   pcpar[4]  = rAbs;
 232   pcpar[5]  = zAbsStart * TMath::Tan(accMax);
 233   pcpar[6]  = pcpar[3]+(zTwoDeg-zAbsStart);
 234   pcpar[7]  = rAbs;
 235   pcpar[8]  = zTwoDeg * TMath::Tan(accMax);
 236   pcpar[9]  = -pcpar[3];
 237   pcpar[10] = zRear * TMath::Tan(accMin);
 238   pcpar[11] = zRear * TMath::Tan(accMax);
 239   gMC->Gsvolu("AITR", "PCON", idtmed[kPb], pcpar, 12);
 240   //
 241   // special Pb medium for last 5 cm of Pb
 242   Float_t zr=zRear-2.-0.001;
 243   cpar[0] = 1.0;
 244   cpar[1] = zr * TMath::Tan(thetaR);
 245   cpar[2] = zr * TMath::Tan(accMax);
 246   cpar[3] = cpar[1] + TMath::Tan(thetaR) * 2;
 247   cpar[4] = cpar[2] + TMath::Tan(accMax) * 2;
 248   gMC->Gsvolu("ARPB", "CONE", idtmed[kPb], cpar, 5);
 249   dz=(zRear-zAbsStart)/2.-cpar[0]-0.001;
 250   gMC->Gspos("ARPB", 1, "AITR", 0., 0., dz, 0, "ONLY");
 251   //
 252   //     concrete cone: concrete
 253   //
 254   pcpar[9]  = pcpar[3]+(zRear-dRear-zAbsStart);
 255   pcpar[10] = (zRear-dRear) * TMath::Tan(accMin);
 256   pcpar[11] = (zRear-dRear) * TMath::Tan(accMax);
 257   gMC->Gsvolu("ACON", "PCON", idtmed[kConcrete+40], pcpar, 12);
 258   gMC->Gspos("ACON", 1, "AITR", 0., 0., 0., 0, "ONLY");
 259   //
 260   //     carbon cone: carbon
 261   //
 262   pcpar[9]  = pcpar[3]+(zAbsCc-zAbsStart);
 263   pcpar[10]  = zAbsCc * TMath::Tan(accMin);
 264   pcpar[11]  = zAbsCc * TMath::Tan(accMax);
 265   gMC->Gsvolu("ACAR", "PCON", idtmed[kC+40], pcpar, 12);
 266   gMC->Gspos("ACAR", 1, "ACON", 0., 0., 0., 0, "ONLY");
 267  //
 268  //     carbon cone outer region
 269  //
 270   cpar[0]  = 10.;
 271   cpar[1]  = rAbs;
 272   cpar[2]  = zAbsStart* TMath::Tan(accMax);
 273   cpar[3]  = rAbs;
 274   cpar[4]  = cpar[2]+2. * cpar[0] * TMath::Tan(accMax);
 275
 276   gMC->Gsvolu("ACAO", "CONE", idtmed[kC], cpar, 5);
 277   dz=-(zRear-zAbsStart)/2.+cpar[0];
 278   gMC->Gspos("ACAO", 1, "ACAR", 0., 0., dz, 0, "ONLY");
 279   //
 280   //     inner W shield
 281   Float_t epsi=0.;
 282   Float_t repsi=1.;
 283
 284   zr=zRear-(dRear-epsi);
 285   cpar[0] = (dRear-epsi)/2.;
 286   cpar[1] = zr * TMath::Tan(accMin);
 287   cpar[2] = zr * TMath::Tan(thetaR*repsi);
 288   cpar[3] = cpar[1] + TMath::Tan(accMin) * (dRear-epsi);
 289   cpar[4] = cpar[2] + TMath::Tan(thetaR*repsi) * (dRear-epsi);
 290   gMC->Gsvolu("ARW0", "CONE", idtmed[kNiCuW+40], cpar, 5);
 291   dz=(zRear-zAbsStart)/2.-cpar[0];
 292   gMC->Gspos("ARW0", 1, "AITR", 0., 0., dz, 0, "ONLY");
 293   //
 294   // special W medium for last 5 cm of W
 295   zr=zRear-5;
 296   cpar[0] = 2.5;
 297   cpar[1] = zr * TMath::Tan(accMin);
 298   cpar[2] = zr * TMath::Tan(thetaR*repsi);
 299   cpar[3] = cpar[1] + TMath::Tan(accMin) * 5.;
 300   cpar[4] = cpar[2] + TMath::Tan(thetaR*repsi) * 5.;
 301   gMC->Gsvolu("ARW1", "CONE", idtmed[kNiCuW+20], cpar, 5);
 302   dz=(dRear-epsi)/2.-cpar[0];
 303   gMC->Gspos("ARW1", 1, "ARW0", 0., 0., dz, 0, "ONLY");
 304   //
 305   // PolyEthylene Layers
 306   Float_t drMin=TMath::Tan(thetaR) * 5;
 307   Float_t drMax=TMath::Tan(accMax) * 5;
 308   gMC->Gsvolu("ARPE", "CONE", idtmed[kPolyCH2], cpar, 0);
 309   cpar[0]=2.5;
 310   { // Begin local scope for i
 311       for (Int_t i=0; i<3; i++) {
 312           zr=zRear-dRear+5+i*10.;
 313           cpar[1] = zr * TMath::Tan(thetaR);
 314           cpar[2] = zr * TMath::Tan(accMax);
 315           cpar[3] = cpar[1] + drMin;
 316           cpar[4] = cpar[2] + drMax;
 317           dz=(zRear-zAbsStart)/2.-cpar[0]-5.-(2-i)*10;
 318           gMC->Gsposp("ARPE", i+1, "AITR", 0., 0., dz, 0, "ONLY",cpar,5);
 319       }
 320   } // End local scope for i
 321   gMC->Gspos("AITR", 1, "ABSS", 0., 0., 0., 0, "ONLY");
 322   dz = (zRear-zAbsStart)/2.+zAbsStart;
 323   gMC->Gspos("ABSM", 1, "ALIC", 0., 0., dz, 0, "ONLY");
 324 //
 325 //
 326 // vacuum system
 327 //
 328 // pipe and heating jackets
 329 //
 330 //
 331 // cylindrical piece
 332   tpar0[2]=(zOpen-zAbsStart)/2;
 333   tpar0[0]=rVacu;
 334   tpar0[1]=rAbs;
 335   gMC->Gsvolu("AV11", "TUBE", idtmed[kSteel+40], tpar0, 3);
 336 //
 337 // insulation
 338
 339   tpar[2]=tpar0[2];
 340   tpar[0]=rVacu+dTube;
 341   tpar[1]=tpar[0]+dInsu;
 342   gMC->Gsvolu("AI11", "TUBE", idtmed[kInsulation+40], tpar, 3);
 343   gMC->Gspos("AI11", 1, "AV11", 0., 0., 0., 0, "ONLY");
 344 //
 345 // clearance
 346   tpar[0]=tpar[1]+dEnve;
 347   tpar[1]=tpar[0]+dFree;
 348   gMC->Gsvolu("AP11", "TUBE", idtmed[kAir+40], tpar, 3);
 349   gMC->Gspos("AP11", 1, "AV11", 0., 0., 0., 0, "ONLY");
 350 //
 351   dz=-(zRear-zAbsStart)/2.+tpar0[2];
 352   gMC->Gspos("AV11", 1, "ABSM", 0., 0., dz, 0, "ONLY");
 353 //
 354 // conical piece
 355
 356   cpar0[0]=(zRear-dRear-zOpen)/2;
 357   cpar0[1]=rVacu-0.05;
 358   cpar0[2]=rAbs;
 359   Float_t dR=2.*cpar0[0]*TMath::Tan(thetaOpen1);
 360   cpar0[3]=cpar0[1]+dR;
 361   cpar0[4]=cpar0[2]+dR;
 362   gMC->Gsvolu("AV21", "CONE", idtmed[kSteel+40], cpar0, 5);
 363   dTube+=0.05;
 364
 365 //
 366 // insulation
 367   cpar[0]=cpar0[0];
 368   cpar[1]=cpar0[1]+dTube;
 369   cpar[2]=cpar0[1]+dTube+dInsu;
 370   cpar[3]=cpar0[3]+dTube;
 371   cpar[4]=cpar0[3]+dTube+dInsu;
 372   gMC->Gsvolu("AI21", "CONE", idtmed[kInsulation+40], cpar, 5);
 373   gMC->Gspos("AI21", 1, "AV21", 0., 0., 0., 0, "ONLY");
 374 //
 375 // clearance
 376   cpar[1]=cpar0[1]+dTube+dInsu+dEnve;
 377   cpar[2]=rAbs;
 378   cpar[3]=cpar0[1]+dTube+dInsu+dEnve+dR;
 379   cpar[4]=rAbs+dR;
 380
 381   gMC->Gsvolu("AP21", "CONE", idtmed[kAir+40], cpar, 5);
 382   gMC->Gspos("AP21", 1, "AV21", 0., 0., 0., 0, "ONLY");
 383
 384   dz=(zRear-zAbsStart)/2.-cpar0[0]-dRear;
 385   gMC->Gspos("AV21", 1, "ABSM", 0., 0., dz, 0, "ONLY");
 386 }
 387
 388 //_____________________________________________________________________________
 389
 390 void AliABSOv0::Init()
 391 {
 392   //
 393   // Initialisation of the muon absorber after it has been built
 394   Int_t i;
 395   //
 396   printf("\n");
 397   for(i=0;i<35;i++) printf("*");
 398   printf(" ABSOv0_INIT ");
 399   for(i=0;i<35;i++) printf("*");
 400   printf("\n");
 401   //
 402   for(i=0;i<80;i++) printf("*");
 403   printf("\n");
 404 }
 405
 406
 407
 408
 409
 410
 411
 412
 413
 414