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.15  2002/10/29 09:53:40  morsch
  19 Constants start with k. Warnings corrected.
  20
  21 Revision 1.14  2002/10/21 09:09:33  alibrary
  22 Removing unused variables
  23
  24 Revision 1.13  2002/10/14 14:57:34  hristov
  25 Merging the VirtualMC branch to the main development branch (HEAD)
  26
  27 Revision 1.11.6.1  2002/07/24 10:08:27  alibrary
  28 Updating VirtualMC
  29
  30 Revision 1.12  2002/07/15 07:56:41  morsch
  31 Steel rear and  Al ring.
  32
  33 Revision 1.11  2001/11/29 14:16:51  morsch
  34 - truncated inner W-cone
  35 - new support structure
  36
  37 Revision 1.10  2001/10/25 08:50:57  morsch
  38 New beamshield geometry with increased tolerances and insulation thickness.
  39
  40 Revision 1.9  2001/05/16 14:57:22  alibrary
  41 New files for folders and Stack
  42
  43 Revision 1.8  2001/01/12 13:16:09  morsch
  44 Store absorber composition information in fMLayers and fZLayers
  45 Rear 25 cm Fe + 35 cm Cu
  46
  47 Revision 1.7  2000/10/02 21:28:15  fca
  48 Removal of useless dependecies via forward declarations
  49
  50 Revision 1.6  2000/06/15 09:40:31  morsch
  51 Obsolete typedef keyword removed
  52
  53 Revision 1.5  2000/06/12 19:39:01  morsch
  54 New structure of beam pipe and heating jacket.
  55
  56 Revision 1.4  2000/04/03 08:13:40  fca
  57 Introduce extra scope for non ANSI compliant C++ compilers
  58
  59 Revision 1.3  2000/01/18 17:49:56  morsch
  60 Serious overlap of ABSM with shield corrected
  61 Small error in ARPB parameters corrected
  62
  63 Revision 1.2  2000/01/13 11:23:59  morsch
  64 Last layer of Pb outer angle corrected
  65
  66 Revision 1.1  2000/01/12 15:39:30  morsch
  67 Standard version of ABSO
  68
  69 */
  70
  71 ///////////////////////////////////////////////////////////////////////////////
  72 //                                                                           //
  73 //  Muon ABSOrber                                                            //
  74 //  This class contains the description of the muon absorber geometry        //
  75 //                                                                           //
  76 //Begin_Html
  77 /*
  78 <img src="picts/AliABSOClass.gif">
  79 </pre>
  80 <br clear=left>
  81 <font size=+2 color=red>
  82 <p>The responsible person for this module is
  83 <a href="mailto:andreas.morsch@cern.ch">Andreas Morsch</a>.
  84 </font>
  85 <pre>
  86 */
  87 //End_Html
  88 //                                                                           //
  89 //                                                                           //
  90 ///////////////////////////////////////////////////////////////////////////////
  91
  92 #include "AliABSOv0.h"
  93 #include "AliRun.h"
  94 #include "AliConst.h"
  95
  96 ClassImp(AliABSOv0)
  97
  98 //_____________________________________________________________________________
  99 AliABSOv0::AliABSOv0()
 100 {
 101   //
 102   // Default constructor
 103   //
 104 }
 105
 106 //_____________________________________________________________________________
 107 AliABSOv0::AliABSOv0(const char *name, const char *title)
 108        : AliABSO(name,title)
 109 {
 110   //
 111   // Standard constructor
 112   //
 113   SetMarkerColor(7);
 114   SetMarkerStyle(2);
 115   SetMarkerSize(0.4);
 116 }
 117
 118 //_____________________________________________________________________________
 119 void AliABSOv0::CreateGeometry()
 120 {
 121     //
 122     // Creation of the geometry of the muon absorber
 123     //
 124     //Begin_Html
 125     /*
 126       <img src="picts/AliABSOv0Tree.gif">
 127     */
 128     //End_Html
 129     //Begin_Html
 130     /*
 131       <img src="picts/AliABSOv0.gif">
 132     */
 133     //End_Html
 134
 135     //
 136     //
 137
 138     enum {kC=1605, kAl=1608, kFe=1609, kCu=1610, kW=1611, kPb=1612,
 139           kNiCuW=1620, kVacuum=1615, kAir=1614, kConcrete=1616,
 140           kPolyCH2=1617, kSteel=1609, kInsulation=1613, kPolyCc=1619};
 141
 142     Int_t *idtmed = fIdtmed->GetArray()-1599;
 143
 144     Float_t par[24], cpar[5], cpar0[5], pcpar[12], tpar[3], tpar0[3];
 145     Float_t dz;
 146     Int_t idrotm[1699];
 147 #include "ABSOSHILConst.h"
 148 #include "ABSOConst.h"
 149 //
 150 // Structure of Tracking Region
 151 //
 152   Float_t dzFe = 25.;
 153
 154 // 3 < theta < 9
 155     fNLayers[0] = 5;
 156     fMLayers[0][0]  = kAir;              fZLayers[0][0] = kZAbsStart;
 157     fMLayers[0][1]  = kC;                fZLayers[0][1] = kZAbsCc;
 158     fMLayers[0][2]  = kConcrete;         fZLayers[0][2] = kZRear-kDRear-dzFe;
 159     fMLayers[0][3]  = kSteel;            fZLayers[0][3] = kZRear-kDRear;
 160     fMLayers[0][4]  = kSteel;            fZLayers[0][4] = kZRear;
 161 // 2 < theta < 3
 162     fNLayers[1] = 6;
 163
 164     fMLayers[1][0] = kAir          ;      fZLayers[1][0] = fZLayers[0][0]-10.;
 165     fMLayers[1][1] = kAl           ;      fZLayers[1][1] = fZLayers[0][0];
 166     fMLayers[1][2] = fMLayers[0][1];      fZLayers[1][2] = fZLayers[0][1];
 167     fMLayers[1][3] = fMLayers[0][2];      fZLayers[1][3] = fZLayers[0][2];
 168     fMLayers[1][4] = fMLayers[0][3];      fZLayers[1][4] = fZLayers[0][3];
 169     fMLayers[1][5] = kNiCuW;              fZLayers[1][5] = fZLayers[0][4];
 170 //
 171
 172     Float_t dTube=0.1;                     // tube thickness
 173     Float_t dInsu=0.5;                     // insulation thickness
 174     Float_t dEnve=0.1;                     // protective envelope thickness
 175     //    Float_t dFree=0.5;                     // clearance thickness
 176
 177
 178 // Mother volume and outer shielding: Pb
 179   par[0]  = 0.;
 180   par[1]  = 360.;
 181   par[2]  = 7.;
 182
 183   par[3]  = -(kZRear-kZAbsStart)/2.;
 184   par[4]  = kRAbs;
 185   par[5]  = kZAbsStart * TMath::Tan(kTheta1);
 186
 187   par[6]  = par[3]+(kZNose-kZAbsStart);
 188   par[7]  = kRAbs;
 189   par[8]  = kZNose * TMath::Tan(kTheta1);
 190
 191   par[9]  = par[3]+(kZConeTPC-kZAbsStart);
 192   par[10] = kRAbs;
 193   par[11] = par[8] + (par[9] - par[6]) * TMath::Tan(kTheta2);
 194
 195   par[12]  = par[3]+(kZOpen-kZAbsStart);
 196   par[13] = kRAbs;
 197   par[14] = par[11] + (par[12] - par[9]) * TMath::Tan(kAccMax);
 198
 199   par[15] = par[3]+(kZRear-kDRear-kZAbsStart);
 200   par[16] = kRAbs   + (par[15] - par[12]) * TMath::Tan(kThetaOpen1) ;
 201   par[17] = par[14] + (par[15] - par[12]) * TMath::Tan(kAccMax);
 202
 203   par[18] = par[3]+(kZRear-kDRear-kZAbsStart);
 204   par[19] = (kZRear-kDRear) * TMath::Tan(kAccMin);
 205   par[20] = par[14] + (par[18] - par[12]) * TMath::Tan(kAccMax);
 206
 207   par[21] = -par[3];
 208   par[22] =  kZRear* TMath::Tan(kAccMin);
 209   par[23] = par[20] + (par[21] - par[18]) * TMath::Tan(kAccMax);
 210   gMC->Gsvolu("ABSS", "PCON", idtmed[kPb], par, 24);
 211   { // Begin local scope for i
 212       for (Int_t i=4; i<18; i+=3) par[i]  = 0;
 213   } // End local scope for i
 214   gMC->Gsvolu("ABSM", "PCON", idtmed[kVacuum+40], par, 24);
 215   gMC->Gspos("ABSS", 1, "ABSM", 0., 0., 0., 0, "ONLY");
 216
 217 //
 218 // Steel envelope
 219 //
 220   par[4] = par[5] -kDSteel;
 221   par[7] = par[8] -kDSteel;
 222   par[10]= par[11]-kDSteel;
 223   par[13]= par[14]-kDSteel;
 224   par[16]= par[17]-kDSteel;
 225   par[19]= par[20]-kDSteel;
 226   par[22]= par[23]-kDSteel;
 227   gMC->Gsvolu("ABST", "PCON", idtmed[kSteel], par, 24);
 228   gMC->Gspos("ABST", 1, "ABSS", 0., 0., 0., 0, "ONLY");
 229 //
 230 // Polyethylene shield
 231 //
 232   cpar[0] = (kZRear - kZConeTPC) / 2.;
 233   cpar[1] = kZConeTPC * TMath::Tan(kAccMax);
 234   cpar[2] = cpar[1] + kDPoly;
 235   cpar[3] = kZRear * TMath::Tan(kAccMax);
 236   cpar[4] = cpar[3] + kDPoly;
 237   gMC->Gsvolu("APOL", "CONE", idtmed[kPolyCH2+40], cpar, 5);
 238   dz = (kZRear-kZAbsStart)/2.-cpar[0];
 239   gMC->Gspos("APOL", 1, "ABSS", 0., 0., dz, 0, "ONLY");
 240
 241 //
 242 // Tungsten nose to protect TPC
 243 //
 244   cpar[0] = (kZNose - kZAbsStart) / 2.;
 245   cpar[1] = kZAbsStart * TMath::Tan(kAccMax);
 246   cpar[2] = kZAbsStart * TMath::Tan(kTheta1)-kDSteel;
 247   cpar[3] = kZNose * TMath::Tan(kAccMax);
 248   cpar[4] = kZNose * TMath::Tan(kTheta1)-kDSteel;
 249   gMC->Gsvolu("ANOS", "CONE", idtmed[kW], cpar, 5);
 250   //
 251   dz = -(kZRear-kZAbsStart)/2.+cpar[0];
 252   gMC->Gspos("ANOS", 1, "ABSS", 0., 0., dz, 0, "ONLY");
 253   //
 254   // Tungsten inner shield
 255   //
 256   Float_t zW = kZTwoDeg+.1;
 257   Float_t dZ = zW+(kZRear-kDRear-zW)/2.;
 258   //
 259   pcpar[0]  = 0.;
 260   pcpar[1]  = 360.;
 261   pcpar[2]  = 3.;
 262   pcpar[3]  = zW-dZ;
 263   pcpar[4]  = kRAbs;
 264   pcpar[5]  = zW * TMath::Tan(kAccMin);
 265   pcpar[6]  = kZOpen-dZ;
 266   pcpar[7]  = kRAbs;
 267   pcpar[8]  = kZOpen * TMath::Tan(kAccMin);
 268   pcpar[9]  = kZRear-kDRear-dZ;
 269   pcpar[10] = kRAbs+(kZRear-kDRear-kZOpen) * TMath::Tan(kThetaOpen1);
 270   pcpar[11] = (kZRear-kDRear) * TMath::Tan(kAccMin);
 271
 272   gMC->Gsvolu("AWIN", "PCON", idtmed[kNiCuW+40], pcpar, 12);
 273   dz=(zW+kZRear-kDRear)/2-(kZAbsStart+kZRear)/2.;
 274   gMC->Gspos("AWIN", 1, "ABSS", 0., 0., dz, 0, "ONLY");
 275 //
 276 // First part replaced by Carbon
 277 //
 278   cpar[0] = (200.-zW)/2.;
 279   cpar[1] = kRAbs;
 280   cpar[2] = pcpar[5];
 281   cpar[3] = kRAbs;
 282   cpar[4] = 200. * TMath::Tan(kAccMin);
 283   gMC->Gsvolu("ACNO", "CONE", idtmed[kC], cpar, 5);
 284   dz = zW-dZ+cpar[0];
 285   gMC->Gspos("ACNO", 1, "AWIN", 0., 0., dz, 0, "ONLY");
 286
 287 /*
 288   Float_t zWW = 383.5;
 289   cpar[0] = (kZRear-kDRear-zWW)/2.;
 290   cpar[1] = kRAbs + (zWW-kZOpen) *  TMath::Tan(kThetaOpen1);
 291   cpar[2] =  zWW * TMath::Tan(kAccMin);
 292   cpar[3] = pcpar[10];
 293   cpar[4] = pcpar[11];
 294   gMC->Gsvolu("AWNO", "CONE", idtmed[kCu+40], cpar, 5);
 295   dz = zWW-dZ+cpar[0];
 296
 297   gMC->Gspos("AWNO", 1, "AWIN", 0., 0., dz, 0, "ONLY");
 298 */
 299   //
 300   //     Inner tracking region
 301   //
 302   //     mother volume: Cu
 303   //
 304   //
 305   pcpar[0]  = 0.;
 306   pcpar[1]  = 360.;
 307   pcpar[2]  = 3.;
 308   pcpar[3]  = -(kZRear-kZAbsStart)/2.;
 309   pcpar[4]  = kRAbs;
 310   pcpar[5]  = kZAbsStart * TMath::Tan(kAccMax);
 311   pcpar[6]  = pcpar[3]+(kZTwoDeg-kZAbsStart);
 312   pcpar[7]  = kRAbs;
 313   pcpar[8]  = kZTwoDeg * TMath::Tan(kAccMax);
 314   pcpar[9]  = -pcpar[3];
 315   pcpar[10] = kZRear * TMath::Tan(kAccMin);
 316   pcpar[11] = kZRear * TMath::Tan(kAccMax);
 317   gMC->Gsvolu("AITR", "PCON", idtmed[fMLayers[0][4]], pcpar, 12);
 318   //
 319   // special Pb medium for last 5 cm of Pb
 320   Float_t zr=kZRear-2.-0.001;
 321   cpar[0] = 1.0;
 322   cpar[1] = zr * TMath::Tan(kThetaR);
 323   cpar[2] = zr * TMath::Tan(kAccMax);
 324   cpar[3] = cpar[1] + TMath::Tan(kThetaR) * 2;
 325   cpar[4] = cpar[2] + TMath::Tan(kAccMax) * 2;
 326   gMC->Gsvolu("ARPB", "CONE", idtmed[fMLayers[0][4]], cpar, 5);
 327   dz=(kZRear-kZAbsStart)/2.-cpar[0]-0.001;
 328   gMC->Gspos("ARPB", 1, "AITR", 0., 0., dz, 0, "ONLY");
 329   //
 330   //     concrete cone: concrete
 331   //
 332   pcpar[9]  = pcpar[3]+(kZRear-kDRear-kZAbsStart);
 333   pcpar[10] = (kZRear-kDRear) * TMath::Tan(kAccMin);
 334   pcpar[11] = (kZRear-kDRear) * TMath::Tan(kAccMax);
 335   gMC->Gsvolu("ACON", "PCON", idtmed[fMLayers[0][2]+40], pcpar, 12);
 336   gMC->Gspos("ACON", 1, "AITR", 0., 0., 0., 0, "ONLY");
 337 //
 338 //    Fe Cone
 339 //
 340   zr = kZRear-kDRear-dzFe;
 341   cpar[0]  = dzFe/2.;
 342   cpar[1] = zr * TMath::Tan(kAccMin);
 343   cpar[2] = zr * TMath::Tan(kAccMax);
 344   cpar[3] = cpar[1] + TMath::Tan(kAccMin) * dzFe;
 345   cpar[4] = cpar[2] + TMath::Tan(kAccMax) * dzFe;
 346   gMC->Gsvolu("ACFE", "CONE",idtmed[fMLayers[0][3]], cpar, 5);
 347
 348   dz = (kZRear-kZAbsStart)/2.-kDRear-dzFe/2.;
 349
 350   gMC->Gspos("ACFE", 1, "ACON", 0., 0., dz, 0, "ONLY");
 351
 352
 353   //
 354   //
 355   //     carbon cone: carbon
 356   //
 357   pcpar[9]  = pcpar[3]+(kZAbsCc-kZAbsStart);
 358   pcpar[10]  = kZAbsCc * TMath::Tan(kAccMin);
 359   pcpar[11]  = kZAbsCc * TMath::Tan(kAccMax);
 360   gMC->Gsvolu("ACAR", "PCON", idtmed[fMLayers[0][1]+40], pcpar, 12);
 361   gMC->Gspos("ACAR", 1, "ACON", 0., 0., 0., 0, "ONLY");
 362  //
 363  //     carbon cone outer region
 364  //
 365   cpar[0]  = 10.;
 366   cpar[1]  = kRAbs;
 367   cpar[2]  = kZAbsStart* TMath::Tan(kAccMax);
 368   cpar[3]  = kRAbs;
 369   cpar[4]  = cpar[2]+2. * cpar[0] * TMath::Tan(kAccMax);
 370
 371   gMC->Gsvolu("ACAO", "CONE", idtmed[fMLayers[0][1]], cpar, 5);
 372   dz=-(kZRear-kZAbsStart)/2.+cpar[0];
 373   gMC->Gspos("ACAO", 1, "ACAR", 0., 0., dz, 0, "ONLY");
 374   //
 375   //     inner W shield
 376   Float_t epsi=0.;
 377   Float_t repsi=1.;
 378
 379   zr=kZRear-(kDRear-epsi);
 380   cpar[0] = (kDRear-epsi)/2.;
 381   cpar[1] = zr * TMath::Tan(kAccMin);
 382   cpar[2] = zr * TMath::Tan(kThetaR*repsi);
 383   cpar[3] = cpar[1] + TMath::Tan(kAccMin) * (kDRear-epsi);
 384   cpar[4] = cpar[2] + TMath::Tan(kThetaR*repsi) * (kDRear-epsi);
 385   gMC->Gsvolu("ARW0", "CONE", idtmed[fMLayers[1][4]+40], cpar, 5);
 386   dz=(kZRear-kZAbsStart)/2.-cpar[0];
 387   gMC->Gspos("ARW0", 1, "AITR", 0., 0., dz, 0, "ONLY");
 388   //
 389   // special W medium for last 5 cm of W
 390   zr=kZRear-5;
 391   cpar[0] = 2.5;
 392   cpar[1] = zr * TMath::Tan(kAccMin);
 393   cpar[2] = zr * TMath::Tan(kThetaR*repsi);
 394   cpar[3] = cpar[1] + TMath::Tan(kAccMin) * 5.;
 395   cpar[4] = cpar[2] + TMath::Tan(kThetaR*repsi) * 5.;
 396   gMC->Gsvolu("ARW1", "CONE", idtmed[fMLayers[1][4]+20], cpar, 5);
 397   dz=(kDRear-epsi)/2.-cpar[0];
 398   gMC->Gspos("ARW1", 1, "ARW0", 0., 0., dz, 0, "ONLY");
 399   //
 400   // Cu
 401   Float_t drMin=TMath::Tan(kThetaR) * 5;
 402   Float_t drMax=TMath::Tan(kAccMax) * 5;
 403   gMC->Gsvolu("ARPE", "CONE", idtmed[fMLayers[0][4]], cpar, 0);
 404   cpar[0]=2.5;
 405   { // Begin local scope for i
 406       for (Int_t i=0; i<3; i++) {
 407           zr=kZRear-kDRear+5+i*10.;
 408           cpar[1] = zr * TMath::Tan(kThetaR);
 409           cpar[2] = zr * TMath::Tan(kAccMax);
 410           cpar[3] = cpar[1] + drMin;
 411           cpar[4] = cpar[2] + drMax;
 412           dz=(kZRear-kZAbsStart)/2.-cpar[0]-5.-(2-i)*10;
 413           gMC->Gsposp("ARPE", i+1, "AITR", 0., 0., dz, 0, "ONLY",cpar,5);
 414       }
 415   } // End local scope for i
 416   gMC->Gspos("AITR", 1, "ABSS", 0., 0., 0., 0, "ONLY");
 417   dz = (kZRear-kZAbsStart)/2.+kZAbsStart;
 418   gMC->Gspos("ABSM", 1, "ALIC", 0., 0., dz, 0, "ONLY");
 419 //
 420 //
 421 // vacuum system
 422 //
 423 // pipe and heating jackets
 424 //
 425 //
 426 // cylindrical piece
 427   tpar0[2]=(kZOpen-kZAbsStart)/2;
 428   tpar0[0]=kRVacu;
 429   tpar0[1]=kRVacu+dTube+dInsu+dEnve;
 430   gMC->Gsvolu("AV11", "TUBE", idtmed[kSteel+40], tpar0, 3);
 431 //
 432 // insulation
 433
 434   tpar[2]=tpar0[2];
 435   tpar[0]=kRVacu+dTube;
 436   tpar[1]=tpar[0]+dInsu;
 437   gMC->Gsvolu("AI11", "TUBE", idtmed[kInsulation+40], tpar, 3);
 438   gMC->Gspos("AI11", 1, "AV11", 0., 0., 0., 0, "ONLY");
 439 //
 440   dz=-(kZRear-kZAbsStart)/2.+tpar0[2];
 441   gMC->Gspos("AV11", 1, "ABSM", 0., 0., dz, 0, "ONLY");
 442 //
 443 // conical piece
 444
 445   cpar0[0]=(kZRear-kDRear-kZOpen)/2;
 446   cpar0[1]= kRVacu-0.05;
 447   cpar0[2]= kRVacu+dTube+dInsu+dEnve;
 448   Float_t dR=2.*cpar0[0]*TMath::Tan(kThetaOpen1);
 449   cpar0[3]=cpar0[1]+dR;
 450   cpar0[4]=cpar0[2]+dR;
 451   gMC->Gsvolu("AV21", "CONE", idtmed[kSteel+40], cpar0, 5);
 452   dTube+=0.05;
 453
 454 //
 455 // insulation
 456   cpar[0]=cpar0[0];
 457   cpar[1]=cpar0[1]+dTube;
 458   cpar[2]=cpar0[1]+dTube+dInsu;
 459   cpar[3]=cpar0[3]+dTube;
 460   cpar[4]=cpar0[3]+dTube+dInsu;
 461   gMC->Gsvolu("AI21", "CONE", idtmed[kInsulation+40], cpar, 5);
 462   gMC->Gspos("AI21", 1, "AV21", 0., 0., 0., 0, "ONLY");
 463
 464   dz=(kZRear-kZAbsStart)/2.-cpar0[0]-kDRear;
 465   gMC->Gspos("AV21", 1, "ABSM", 0., 0., dz, 0, "ONLY");
 466 //
 467 // Support cone
 468
 469   par[0]  =  22.5;
 470   par[1]  = 360.0;
 471   par[2]  =   8.0;
 472   par[3]  =   4.0;
 473
 474   par[4]  = kZRear;
 475   par[5]  = 100.;
 476   par[6]  = 180.;
 477
 478   par[7]  = kZRear+20.;
 479   par[8]  = 100.;
 480   par[9]  = 180.;
 481
 482   par[10] = kZRear+20.;
 483   par[11] = 178.;
 484   par[12] = 180.;
 485
 486   par[13] = 600.;
 487   par[14] = 178.;
 488   par[15] = 180.;
 489
 490
 491   gMC->Gsvolu("ASSS", "PGON", idtmed[kAl], par, 16);
 492   gMC->Gspos("ASSS", 1, "ALIC", 0., 0., 0., 0, "ONLY");
 493
 494   Float_t trap[11];
 495   trap[ 0] = (530.-170.)/2.;
 496   trap[ 2] = 0.;
 497   trap[ 3] = 2.;
 498   trap[ 4] = (600.-(kZRear+2.))/2.;;
 499   trap[ 5] = trap[4];
 500   trap[ 6] = 0.;
 501   trap[ 7] = 2.;
 502   trap[ 8] = 5.;
 503   trap[ 9] = 5.;
 504   trap[10] = 0.;
 505   trap[ 1] = -TMath::ATan((trap[4]-trap[8])/2./trap[0])*180./TMath::Pi();
 506   AliMatrix(idrotm[1600], 180., 0., 90., 0., 90., 90.);
 507   AliMatrix(idrotm[1601], 180., 0., 90., 0., 90., 270.);
 508   gMC->Gsvolu("ASST", "TRAP", idtmed[kSteel], trap, 11);
 509   dz = (600.+kZRear+2.)/2.+(trap[4]-trap[8])/2.;
 510   //  Float_t dy =  170.+trap[0];
 511
 512 //  gMC->Gspos("ASST", 1, "ALIC", 0.,  dy, dz, idrotm[1600], "ONLY");
 513 //  gMC->Gspos("ASST", 2, "ALIC", 0., -dy, dz, idrotm[1601], "ONLY");
 514 }
 515
 516 //_____________________________________________________________________________
 517
 518 void AliABSOv0::Init()
 519 {
 520   //
 521   // Initialisation of the muon absorber after it has been built
 522   Int_t i;
 523   //
 524   if(fDebug) {
 525     printf("\n%s: ",ClassName());
 526     for(i=0;i<35;i++) printf("*");
 527     printf(" ABSOv0_INIT ");
 528     for(i=0;i<35;i++) printf("*");
 529     printf("\n%s: ",ClassName());
 530     //
 531     for(i=0;i<80;i++) printf("*");
 532     printf("\n");
 533   }
 534 }
 535
 536
 537
 538
 539
 540
 541
 542
 543
 544