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