ZDC/AliZDCv3.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 ///////////////////////////////////////////////////////////////////////
  18 //                                                                   //
  19 //              AliZDCv3 --- new ZDC geometry                        //
  20 //          with both ZDC arms geometry implemented                  //
  21 //                                                                   //
  22 ///////////////////////////////////////////////////////////////////////
  23
  24 // --- Standard libraries
  25 #include "stdio.h"
  26
  27 // --- ROOT system
  28 #include <TMath.h>
  29 #include <TRandom.h>
  30 #include <TSystem.h>
  31 #include <TTree.h>
  32 #include <TVirtualMC.h>
  33 #include <TGeoManager.h>
  34 #include <TGeoMatrix.h>
  35 #include <TGeoTube.h>
  36 #include <TGeoCone.h>
  37 #include <TGeoShape.h>
  38 #include <TGeoScaledShape.h>
  39 #include <TGeoCompositeShape.h>
  40 #include <TParticle.h>
  41
  42 // --- AliRoot classes
  43 #include "AliLog.h"
  44 #include "AliConst.h"
  45 #include "AliMagF.h"
  46 #include "AliRun.h"
  47 #include "AliZDCv3.h"
  48 #include "AliMC.h"
  49
  50 class  AliZDCHit;
  51 class  AliPDG;
  52 class  AliDetector;
  53
  54
  55 ClassImp(AliZDCv3)
  56
  57 //_____________________________________________________________________________
  58 AliZDCv3::AliZDCv3() :
  59   AliZDC(),
  60   fMedSensF1(0),
  61   fMedSensF2(0),
  62   fMedSensZP(0),
  63   fMedSensZN(0),
  64   fMedSensZEM(0),
  65   fMedSensGR(0),
  66   fMedSensPI(0),
  67   fMedSensTDI(0),
  68   fMedSensVColl(0),
  69   fMedSensLumi(0),
  70   fNalfan(0),
  71   fNalfap(0),
  72   fNben(0),
  73   fNbep(0),
  74   fZEMLength(0),
  75   fpLostITC(0),
  76   fpLostD1C(0),
  77   fpcVCollC(0),
  78   fpDetectedC(0),
  79   fnDetectedC(0),
  80   fpLostITA(0),
  81   fpLostD1A(0),
  82   fpLostTDI(0),
  83   fpcVCollA(0),
  84   fpDetectedA(0),
  85   fnDetectedA(0),
  86   fVCollSideCAperture(7./2.),
  87   fVCollSideCCentreY(0.),
  88   fVCollSideAAperture(7./2.),
  89   fVCollSideACentreY(0.),
  90   fLumiLength(15.)
  91 {
  92   //
  93   // Default constructor for Zero Degree Calorimeter
  94   //
  95
  96 }
  97
  98 //_____________________________________________________________________________
  99 AliZDCv3::AliZDCv3(const char *name, const char *title) :
 100   AliZDC(name,title),
 101   fMedSensF1(0),
 102   fMedSensF2(0),
 103   fMedSensZP(0),
 104   fMedSensZN(0),
 105   fMedSensZEM(0),
 106   fMedSensGR(0),
 107   fMedSensPI(0),
 108   fMedSensTDI(0),
 109   fMedSensVColl(0),
 110   fMedSensLumi(0),
 111   fNalfan(90),
 112   fNalfap(90),
 113   fNben(18),
 114   fNbep(28),
 115   fZEMLength(0),
 116   fpLostITC(0),
 117   fpLostD1C(0),
 118   fpcVCollC(0),
 119   fpDetectedC(0),
 120   fnDetectedC(0),
 121   fpLostITA(0),
 122   fpLostD1A(0),
 123   fpLostTDI(0),
 124   fpcVCollA(0),
 125   fpDetectedA(0),
 126   fnDetectedA(0),
 127   fVCollSideCAperture(7./2.),
 128   fVCollSideCCentreY(0.),
 129   fVCollSideAAperture(7./2.),
 130   fVCollSideACentreY(0.),
 131   fLumiLength(15.)
 132 {
 133   //
 134   // Standard constructor for Zero Degree Calorimeter
 135   //
 136   //
 137   // Check that DIPO, ABSO, DIPO and SHIL is there (otherwise tracking is wrong!!!)
 138
 139   AliModule* pipe=gAlice->GetModule("PIPE");
 140   AliModule* abso=gAlice->GetModule("ABSO");
 141   AliModule* dipo=gAlice->GetModule("DIPO");
 142   AliModule* shil=gAlice->GetModule("SHIL");
 143   if((!pipe) || (!abso) || (!dipo) || (!shil)) {
 144     Error("Constructor","ZDC needs PIPE, ABSO, DIPO and SHIL!!!\n");
 145     exit(1);
 146   }
 147   //
 148   Int_t ip,jp,kp;
 149   for(ip=0; ip<4; ip++){
 150      for(kp=0; kp<fNalfap; kp++){
 151         for(jp=0; jp<fNbep; jp++){
 152            fTablep[ip][kp][jp] = 0;
 153         }
 154      }
 155   }
 156   Int_t in,jn,kn;
 157   for(in=0; in<4; in++){
 158      for(kn=0; kn<fNalfan; kn++){
 159         for(jn=0; jn<fNben; jn++){
 160            fTablen[in][kn][jn] = 0;
 161         }
 162      }
 163   }
 164   //
 165   // Parameters for hadronic calorimeters geometry
 166   // Positions updated after post-installation measurements
 167   fDimZN[0] = 3.52;
 168   fDimZN[1] = 3.52;
 169   fDimZN[2] = 50.;
 170   fDimZP[0] = 11.2;
 171   fDimZP[1] = 6.;
 172   fDimZP[2] = 75.;
 173   fPosZNC[0] = 0.;
 174   fPosZNC[1] = 0.;
 175   fPosZNC[2] = -11397.3;
 176   fPosZPC[0] = 24.35;
 177   fPosZPC[1] = 0.;
 178   fPosZPC[2] = -11389.3;
 179   fPosZNA[0] = 0.;
 180   fPosZNA[1] = 0.;
 181   fPosZNA[2] = 11395.8;
 182   fPosZPA[0] = 24.35;
 183   fPosZPA[1] = 0.;
 184   fPosZPA[2] = 11387.8;
 185   fFibZN[0] = 0.;
 186   fFibZN[1] = 0.01825;
 187   fFibZN[2] = 50.;
 188   fFibZP[0] = 0.;
 189   fFibZP[1] = 0.0275;
 190   fFibZP[2] = 75.;
 191   // Parameters for EM calorimeter geometry
 192   fPosZEM[0] = 8.5;
 193   fPosZEM[1] = 0.;
 194   fPosZEM[2] = 735.;
 195   Float_t kDimZEMPb  = 0.15*(TMath::Sqrt(2.));  // z-dimension of the Pb slice
 196   Float_t kDimZEMAir = 0.001;                   // scotch
 197   Float_t kFibRadZEM = 0.0315;                  // External fiber radius (including cladding)
 198   Int_t   kDivZEM[3] = {92, 0, 20};             // Divisions for EM detector
 199   Float_t kDimZEM0 = 2*kDivZEM[2]*(kDimZEMPb+kDimZEMAir+kFibRadZEM*(TMath::Sqrt(2.)));
 200   fZEMLength = kDimZEM0;
 201
 202 }
 203
 204 //_____________________________________________________________________________
 205 void AliZDCv3::CreateGeometry()
 206 {
 207   //
 208   // Create the geometry for the Zero Degree Calorimeter version 2
 209   //* Initialize COMMON block ZDC_CGEOM
 210   //*
 211
 212   CreateBeamLine();
 213   CreateZDC();
 214 }
 215
 216 //_____________________________________________________________________________
 217 void AliZDCv3::CreateBeamLine()
 218 {
 219   //
 220   // Create the beam line elements
 221   //
 222
 223   Double_t zd1, zd2, zCorrDip, zInnTrip, zD1, zD2;
 224   Double_t conpar[9], tubpar[3], tubspar[5], boxpar[3];
 225
 226   //-- rotation matrices for the legs
 227   Int_t irotpipe1, irotpipe2;
 228   gMC->Matrix(irotpipe1,90.-1.0027,0.,90.,90.,1.0027,180.);
 229   gMC->Matrix(irotpipe2,90.+1.0027,0.,90.,90.,1.0027,0.);
 230
 231   //
 232   Int_t *idtmed = fIdtmed->GetArray();
 233
 234   ////////////////////////////////////////////////////////////////
 235   //                                                            //
 236   //                SIDE C - RB26 (dimuon side)                 //
 237   //                                                            //
 238   ///////////////////////////////////////////////////////////////
 239
 240
 241   // -- Mother of the ZDCs (Vacuum PCON)
 242   zd1 = 1921.6;
 243
 244   conpar[0] = 0.;
 245   conpar[1] = 360.;
 246   conpar[2] = 2.;
 247   conpar[3] = -13500.;
 248   conpar[4] = 0.;
 249   conpar[5] = 55.;
 250   conpar[6] = -zd1;
 251   conpar[7] = 0.;
 252   conpar[8] = 55.;
 253   gMC->Gsvolu("ZDCC", "PCON", idtmed[10], conpar, 9);
 254   gMC->Gspos("ZDCC", 1, "ALIC", 0., 0., 0., 0, "ONLY");
 255
 256
 257   // -- BEAM PIPE from compensator dipole to the beginning of D1)
 258   tubpar[0] = 6.3/2.;
 259   tubpar[1] = 6.7/2.;
 260   // From beginning of ZDC volumes to beginning of D1
 261   tubpar[2] = (5838.3-zd1)/2.;
 262   gMC->Gsvolu("QT01", "TUBE", idtmed[7], tubpar, 3);
 263   gMC->Gspos("QT01", 1, "ZDCC", 0., 0., -tubpar[2]-zd1, 0, "ONLY");
 264   // Ch.debug
 265   //printf("    QT01 TUBE pipe from z = %1.2f to z= %1.2f (D1 beg.)\n",-zd1,-2*tubpar[2]-zd1);
 266
 267   //-- BEAM PIPE from the end of D1 to the beginning of D2)
 268
 269   //-- FROM MAGNETIC BEGINNING OF D1 TO MAGNETIC END OF D1
 270   //--  Cylindrical pipe (r = 3.47) + conical flare
 271   // -> Beginning of D1
 272   zd1 += 2.*tubpar[2];
 273
 274   tubpar[0] = 6.94/2.;
 275   tubpar[1] = 7.34/2.;
 276   tubpar[2] = (6909.8-zd1)/2.;
 277   gMC->Gsvolu("QT02", "TUBE", idtmed[7], tubpar, 3);
 278   gMC->Gspos("QT02", 1, "ZDCC", 0., 0., -tubpar[2]-zd1, 0, "ONLY");
 279   // Ch.debug
 280   //printf("    QT02 TUBE pipe from z = %1.2f to z= %1.2f\n",-zd1,-2*tubpar[2]-zd1);
 281
 282   zd1 += 2.*tubpar[2];
 283
 284   tubpar[0] = 9./2.;
 285   tubpar[1] = 9.6/2.;
 286   tubpar[2] = (7022.8-zd1)/2.;
 287   gMC->Gsvolu("QT03", "TUBE", idtmed[7], tubpar, 3);
 288   gMC->Gspos("QT03", 1, "ZDCC", 0., 0., -tubpar[2]-zd1, 0, "ONLY");
 289   // Ch.debug
 290   //printf("    QT03 TUBE pipe from z = %1.2f to z= %1.2f\n",-zd1,-2*tubpar[2]-zd1);
 291
 292   zd1 += 2.*tubpar[2];
 293
 294   conpar[0] = 39.2/2.;
 295   conpar[1] = 18./2.;
 296   conpar[2] = 18.6/2.;
 297   conpar[3] = 9./2.;
 298   conpar[4] = 9.6/2.;
 299   gMC->Gsvolu("QC01", "CONE", idtmed[7], conpar, 5);
 300   gMC->Gspos("QC01", 1, "ZDCC", 0., 0., -conpar[0]-zd1, 0, "ONLY");
 301   // Ch.debug
 302   //printf("    QC01 CONE pipe from z = %1.2f to z= %1.2f\n",-zd1,-2*conpar[0]-zd1);
 303
 304   zd1 += conpar[0] * 2.;
 305
 306   // ******************************************************
 307   // N.B.-> according to last vacuum layout
 308   // private communication by D. Macina, mail 27/1/2009
 309   // ******************************************************
 310   // 2nd section of VCTCQ+VAMTF+TCTVB+VAMTF+TCLIA+VAMTF+1st part of VCTCP
 311   Float_t totLength1 = 160.8 + 78. + 148. + 78. + 148. + 78. + 9.3;
 312   //
 313   tubpar[0] = 18.6/2.;
 314   tubpar[1] = 7.6/2.;
 315   tubpar[2] = totLength1/2.;
 316 //  gMC->Gsvolu("QE01", "ELTU", idtmed[7], tubpar, 3);
 317   // temporary replace with a scaled tube (AG)
 318   TGeoTube *tubeQE01 = new TGeoTube(0.,tubpar[0],tubpar[2]);
 319   TGeoScale *scaleQE01 = new TGeoScale(1., tubpar[1]/tubpar[0], 1.);
 320   TGeoScaledShape *sshapeQE01 = new TGeoScaledShape(tubeQE01, scaleQE01);
 321   new TGeoVolume("QE01", sshapeQE01, gGeoManager->GetMedium(idtmed[7]));
 322
 323   tubpar[0] = 18.0/2.;
 324   tubpar[1] = 7.0/2.;
 325   tubpar[2] = totLength1/2.;
 326 //  gMC->Gsvolu("QE02", "ELTU", idtmed[10], tubpar, 3);
 327   // temporary replace with a scaled tube (AG)
 328   TGeoTube *tubeQE02 = new TGeoTube(0.,tubpar[0],tubpar[2]);
 329   TGeoScale *scaleQE02 = new TGeoScale(1., tubpar[1]/tubpar[0], 1.);
 330   TGeoScaledShape *sshapeQE02 = new TGeoScaledShape(tubeQE02, scaleQE02);
 331   new TGeoVolume("QE02", sshapeQE02, gGeoManager->GetMedium(idtmed[10]));
 332
 333   gMC->Gspos("QE01", 1, "ZDCC", 0., 0., -tubpar[2]-zd1, 0, "ONLY");
 334   gMC->Gspos("QE02", 1, "QE01", 0., 0., 0., 0, "ONLY");
 335   // Ch.debug
 336   //printf("    QE01 ELTU from z = %1.2f to z= %1.2f\n",-zd1,-2*tubpar[2]-zd1);
 337
 338   // Vertical collimator jaws (defined ONLY if fVCollAperture<3.5!)
 339   if(fVCollSideCAperture<3.5){
 340     boxpar[0] = 5.4/2.;
 341     boxpar[1] = (3.5-fVCollSideCAperture-fVCollSideCCentreY-0.7)/2.;
 342     if(boxpar[1]<0.) boxpar[1]=0.;
 343     boxpar[2] = 124.4/2.;
 344     printf("\n  AliZDCv3 -> Setting SideC VCollimator jaw: aperture %1.2f center %1.2f mod.thickness %1.3f\n\n",
 345         2*fVCollSideCAperture,fVCollSideCCentreY,2*boxpar[1]);
 346     gMC->Gsvolu("QCVC" , "BOX ", idtmed[13], boxpar, 3);
 347     gMC->Gspos("QCVC", 1, "QE02", -boxpar[0],  fVCollSideCAperture+fVCollSideCCentreY+boxpar[1], -totLength1/2.+160.8+78.+148./2., 0, "ONLY");
 348     gMC->Gspos("QCVC", 2, "QE02", -boxpar[0], -fVCollSideCAperture+fVCollSideCCentreY-boxpar[1], -totLength1/2.+160.8+78.+148./2., 0, "ONLY");
 349   }
 350
 351   zd1 += tubpar[2] * 2.;
 352
 353   // 2nd part of VCTCP
 354   conpar[0] = 31.5/2.;
 355   conpar[1] = 21.27/2.;
 356   conpar[2] = 21.87/2.;
 357   conpar[3] = 18.0/2.;
 358   conpar[4] = 18.6/2.;
 359   gMC->Gsvolu("QC02", "CONE", idtmed[7], conpar, 5);
 360   gMC->Gspos("QC02", 1, "ZDCC", 0., 0., -conpar[0]-zd1, 0, "ONLY");
 361   // Ch.debug
 362   //printf("    QC02 CONE pipe from z = %1.2f to z= %1.2f\n",-zd1,-2*conpar[0]-zd1);
 363
 364   zd1 += conpar[0] * 2.;
 365
 366   // 3rd section of VCTCP+VCDWC+VMLGB
 367   Float_t totLenght2 = 9.2 + 530.5+40.;
 368   tubpar[0] = 21.2/2.;
 369   tubpar[1] = 21.9/2.;
 370   tubpar[2] = totLenght2/2.;
 371   gMC->Gsvolu("QT04", "TUBE", idtmed[7], tubpar, 3);
 372   gMC->Gspos("QT04", 1, "ZDCC", 0., 0., -tubpar[2]-zd1, 0, "ONLY");
 373   // Ch.debug
 374   //printf("    QT04 TUBE pipe from z = %1.2f to z= %1.2f\n",-zd1,-2*tubpar[2]-zd1);
 375
 376   zd1 += tubpar[2] * 2.;
 377
 378   // First part of VCTCD
 379   // skewed transition cone from ID=212.7 mm to ID=797 mm
 380   conpar[0] = 121./2.;
 381   conpar[1] = 79.7/2.;
 382   conpar[2] = 81.3/2.;
 383   conpar[3] = 21.27/2.;
 384   conpar[4] = 21.87/2.;
 385   gMC->Gsvolu("QC03", "CONE", idtmed[7], conpar, 5);
 386   gMC->Gspos("QC03", 1, "ZDCC", 0., 0., -conpar[0]-zd1, 0, "ONLY");
 387   // Ch.debug
 388   //printf("    QC03 CONE pipe from z = %1.2f to z= %1.2f\n",-zd1,-2*conpar[0]-zd1);
 389
 390   zd1 += 2.*conpar[0];
 391
 392   // VCDGB + 1st part of VCTCH
 393   tubpar[0] = 79.7/2.;
 394   tubpar[1] = 81.3/2.;
 395   tubpar[2] = (5*475.2+97.)/2.;
 396   gMC->Gsvolu("QT05", "TUBE", idtmed[7], tubpar, 3);
 397   gMC->Gspos("QT05", 1, "ZDCC", 0., 0., -tubpar[2]-zd1, 0, "ONLY");
 398   // Ch.debug
 399   //printf("    QT05 TUBE pipe from z = %1.2f to z= %1.2f\n",-zd1,-2*tubpar[2]-zd1);
 400
 401   zd1 += 2.*tubpar[2];
 402
 403   // 2nd part of VCTCH
 404   // Transition from ID=797 mm to ID=196 mm:
 405   // in order to simulate the thin window opened in the transition cone
 406   // we divide the transition cone in three cones:
 407   // (1) 8 mm thick (2) 3 mm thick (3) the third 8 mm thick
 408
 409   // (1) 8 mm thick
 410   conpar[0] = 9.09/2.; // 15 degree
 411   conpar[1] = 74.82868/2.;
 412   conpar[2] = 76.42868/2.; // thickness 8 mm
 413   conpar[3] = 79.7/2.;
 414   conpar[4] = 81.3/2.; // thickness 8 mm
 415   gMC->Gsvolu("QC04", "CONE", idtmed[7], conpar, 5);
 416   gMC->Gspos("QC04", 1, "ZDCC", 0., 0., -conpar[0]-zd1, 0, "ONLY");
 417   // Ch.debug
 418   //printf("    QC04 CONE pipe from z = %1.2f to z= %1.2f\n",-zd1,-2*conpar[0]-zd1);
 419
 420   zd1 += 2.*conpar[0];
 421
 422   // (2) 3 mm thick
 423   conpar[0] = 96.2/2.; // 15 degree
 424   conpar[1] = 23.19588/2.;
 425   conpar[2] = 23.79588/2.; // thickness 3 mm
 426   conpar[3] = 74.82868/2.;
 427   conpar[4] = 75.42868/2.; // thickness 3 mm
 428   gMC->Gsvolu("QC05", "CONE", idtmed[7], conpar, 5);
 429   gMC->Gspos("QC05", 1, "ZDCC", 0., 0., -conpar[0]-zd1, 0, "ONLY");
 430   // Ch.debug
 431   //printf("    QC05 CONE pipe from z = %1.2f to z= %1.2f\n",-zd1,-2*conpar[0]-zd1);
 432
 433   zd1 += 2.*conpar[0];
 434
 435   // (3) 8 mm thick
 436   conpar[0] = 6.71/2.; // 15 degree
 437   conpar[1] = 19.6/2.;
 438   conpar[2] = 21.2/2.;// thickness 8 mm
 439   conpar[3] = 23.19588/2.;
 440   conpar[4] = 24.79588/2.;// thickness 8 mm
 441   gMC->Gsvolu("QC06", "CONE", idtmed[7], conpar, 5);
 442   gMC->Gspos("QC06", 1, "ZDCC", 0., 0., -conpar[0]-zd1, 0, "ONLY");
 443   // Ch.debug
 444   //printf("    QC06 CONE pipe from z = %1.2f to z= %1.2f\n",-zd1,-2*conpar[0]-zd1);
 445
 446   zd1 += 2.*conpar[0];
 447
 448   // VMZAR (5 volumes)
 449   tubpar[0] = 20.2/2.;
 450   tubpar[1] = 20.6/2.;
 451   tubpar[2] = 2.15/2.;
 452   gMC->Gsvolu("QT06", "TUBE", idtmed[7], tubpar, 3);
 453   gMC->Gspos("QT06", 1, "ZDCC", 0., 0., -tubpar[2]-zd1, 0, "ONLY");
 454   // Ch.debug
 455   //printf("    QT06 TUBE pipe from z = %1.2f to z= %1.2f\n",-zd1,-2*tubpar[2]-zd1);
 456
 457   zd1 += 2.*tubpar[2];
 458
 459   conpar[0] = 6.9/2.;
 460   conpar[1] = 23.9/2.;
 461   conpar[2] = 24.3/2.;
 462   conpar[3] = 20.2/2.;
 463   conpar[4] = 20.6/2.;
 464   gMC->Gsvolu("QC07", "CONE", idtmed[7], conpar, 5);
 465   gMC->Gspos("QC07", 1, "ZDCC", 0., 0., -conpar[0]-zd1, 0, "ONLY");
 466   // Ch.debug
 467   //printf("    QC07 CONE pipe from z = %1.2f to z= %1.2f\n",-zd1,-2*conpar[0]-zd1);
 468
 469   zd1 += 2.*conpar[0];
 470
 471   tubpar[0] = 23.9/2.;
 472   tubpar[1] = 25.5/2.;
 473   tubpar[2] = 17.0/2.;
 474   gMC->Gsvolu("QT07", "TUBE", idtmed[7], tubpar, 3);
 475   gMC->Gspos("QT07", 1, "ZDCC", 0., 0., -tubpar[2]-zd1, 0, "ONLY");
 476   // Ch.debug
 477   //printf("    QT07 TUBE pipe from z = %1.2f to z= %1.2f\n",-zd1,-2*tubpar[2]-zd1);
 478
 479   zd1 += 2.*tubpar[2];
 480
 481   conpar[0] = 6.9/2.;
 482   conpar[1] = 20.2/2.;
 483   conpar[2] = 20.6/2.;
 484   conpar[3] = 23.9/2.;
 485   conpar[4] = 24.3/2.;
 486   gMC->Gsvolu("QC08", "CONE", idtmed[7], conpar, 5);
 487   gMC->Gspos("QC08", 1, "ZDCC", 0., 0., -conpar[0]-zd1, 0, "ONLY");
 488   // Ch.debug
 489   //printf("    QC08 CONE pipe from z = %1.2f to z= %1.2f\n",-zd1,-2*conpar[0]-zd1);
 490
 491   zd1 += 2.*conpar[0];
 492
 493   tubpar[0] = 20.2/2.;
 494   tubpar[1] = 20.6/2.;
 495   tubpar[2] = 2.15/2.;
 496   gMC->Gsvolu("QT08", "TUBE", idtmed[7], tubpar, 3);
 497   gMC->Gspos("QT08", 1, "ZDCC", 0., 0., -tubpar[2]-zd1, 0, "ONLY");
 498   // Ch.debug
 499   //printf("    QT08 TUBE pipe from z = %1.2f to z= %1.2f\n",-zd1,-2*tubpar[2]-zd1);
 500
 501   zd1 += 2.*tubpar[2];
 502
 503   // Flange (ID=196 mm)(last part of VMZAR and first part of VCTYB)
 504   tubpar[0] = 19.6/2.;
 505   tubpar[1] = 25.3/2.;
 506   tubpar[2] = 4.9/2.;
 507   gMC->Gsvolu("QT09", "TUBE", idtmed[7], tubpar, 3);
 508   gMC->Gspos("QT09", 1, "ZDCC", 0., 0., -tubpar[2]-zd1, 0, "ONLY");
 509   // Ch.debug
 510   //printf("    QT09 TUBE pipe from z = %1.2f to z= %1.2f\n",-zd1,-2*tubpar[2]-zd1);
 511
 512   zd1 += 2.*tubpar[2];
 513   // Ch.debug
 514   //printf("    Beginning of VCTYB volume @ z = %1.2f \n",-zd1);
 515
 516   // simulation of the trousers (VCTYB)
 517   tubpar[0] = 19.6/2.;
 518   tubpar[1] = 20.0/2.;
 519   tubpar[2] = 3.9/2.;
 520   gMC->Gsvolu("QT10", "TUBE", idtmed[7], tubpar, 3);
 521   gMC->Gspos("QT10", 1, "ZDCC", 0., 0., -tubpar[2]-zd1, 0, "ONLY");
 522   // Ch.debug
 523   //printf("    QT10 TUBE pipe from z = %1.2f to z= %1.2f\n",-zd1,-2*tubpar[2]-zd1);
 524
 525   zd1 += 2.*tubpar[2];
 526
 527   // transition cone from ID=196. to ID=216.6
 528   conpar[0] = 32.55/2.;
 529   conpar[1] = 21.66/2.;
 530   conpar[2] = 22.06/2.;
 531   conpar[3] = 19.6/2.;
 532   conpar[4] = 20.0/2.;
 533   gMC->Gsvolu("QC09", "CONE", idtmed[7], conpar, 5);
 534   gMC->Gspos("QC09", 1, "ZDCC", 0., 0., -conpar[0]-zd1, 0, "ONLY");
 535   // Ch.debug
 536   //printf("    QC09 CONE pipe from z = %1.2f to z= %1.2f\n",-zd1,-2*conpar[0]-zd1);
 537
 538   zd1 += 2.*conpar[0];
 539
 540   // tube
 541   tubpar[0] = 21.66/2.;
 542   tubpar[1] = 22.06/2.;
 543   tubpar[2] = 28.6/2.;
 544   gMC->Gsvolu("QT11", "TUBE", idtmed[7], tubpar, 3);
 545   gMC->Gspos("QT11", 1, "ZDCC", 0., 0., -tubpar[2]-zd1, 0, "ONLY");
 546   // Ch.debug
 547   //printf("    QT11 TUBE pipe from z = %1.2f to z= %1.2f\n",-zd1,-2*tubpar[2]-zd1);
 548
 549   zd1 += 2.*tubpar[2];
 550   // Ch.debug
 551   //printf("    Beginning of recombination chamber @ z = %f \n",-zd1);
 552
 553   // --------------------------------------------------------
 554   // RECOMBINATION CHAMBER IMPLEMENTED USING TGeo CLASSES!!!!
 555   // author: Chiara (August 2008)
 556   // --------------------------------------------------------
 557   // TRANSFORMATION MATRICES
 558   // Combi transformation:
 559   Double_t dx = -3.970000;
 560   Double_t dy = 0.000000;
 561   Double_t dz = 0.0;
 562   // Rotation:
 563   Double_t thx = 84.989100;   Double_t phx = 180.000000;
 564   Double_t thy = 90.000000;   Double_t phy = 90.000000;
 565   Double_t thz = 185.010900;  Double_t phz = 0.000000;
 566   TGeoRotation *rotMatrix1c = new TGeoRotation("c",thx,phx,thy,phy,thz,phz);
 567   // Combi transformation:
 568   dx = -3.970000;
 569   dy = 0.000000;
 570   dz = 0.0;
 571   TGeoCombiTrans *rotMatrix2c = new TGeoCombiTrans("ZDCC_c1", dx,dy,dz,rotMatrix1c);
 572   rotMatrix2c->RegisterYourself();
 573   // Combi transformation:
 574   dx = 3.970000;
 575   dy = 0.000000;
 576   dz = 0.0;
 577   // Rotation:
 578   thx = 95.010900;   phx = 180.000000;
 579   thy = 90.000000;   phy = 90.000000;
 580   thz = 180.-5.010900;    phz = 0.000000;
 581   TGeoRotation *rotMatrix3c = new TGeoRotation("",thx,phx,thy,phy,thz,phz);
 582   TGeoCombiTrans *rotMatrix4c = new TGeoCombiTrans("ZDCC_c2", dx,dy,dz,rotMatrix3c);
 583   rotMatrix4c->RegisterYourself();
 584
 585   // VOLUMES DEFINITION
 586   // Volume: ZDCC
 587   TGeoVolume *pZDCC = gGeoManager->GetVolume("ZDCC");
 588
 589   conpar[0] = (90.1-0.95-0.26-0.0085)/2.;
 590   conpar[1] = 0.0/2.;
 591   conpar[2] = 21.6/2.;
 592   conpar[3] = 0.0/2.;
 593   conpar[4] = 5.8/2.;
 594   new TGeoCone("QCLext", conpar[0],conpar[1],conpar[2],conpar[3],conpar[4]);
 595
 596   conpar[0] = (90.1-0.95-0.26-0.0085)/2.;
 597   conpar[1] = 0.0/2.;
 598   conpar[2] = 21.2/2.;
 599   conpar[3] = 0.0/2.;
 600   conpar[4] = 5.4/2.;
 601   new TGeoCone("QCLint", conpar[0],conpar[1],conpar[2],conpar[3],conpar[4]);
 602
 603   // Outer trousers
 604   TGeoCompositeShape *pOutTrousersC = new TGeoCompositeShape("outTrousersC", "QCLext:ZDCC_c1+QCLext:ZDCC_c2");
 605
 606   // Volume: QCLext
 607   TGeoMedium *medZDCFe = gGeoManager->GetMedium("ZDC_ZIRONT");
 608   TGeoVolume *pQCLext = new TGeoVolume("QCLext",pOutTrousersC, medZDCFe);
 609   pQCLext->SetLineColor(kGreen);
 610   pQCLext->SetVisLeaves(kTRUE);
 611   //
 612   TGeoTranslation *tr1c = new TGeoTranslation(0., 0., (Double_t) -conpar[0]-0.95-zd1);
 613   //printf("    Recombination chamber from z = %1.2f to z= %1.2f\n",-zd1,-2*conpar[0]-0.95-zd1);
 614   //
 615   pZDCC->AddNode(pQCLext, 1, tr1c);
 616   // Inner trousers
 617   TGeoCompositeShape *pIntTrousersC = new TGeoCompositeShape("intTrousersC", "QCLint:ZDCC_c1+QCLint:ZDCC_c2");
 618   // Volume: QCLint
 619   TGeoMedium *medZDCvoid = gGeoManager->GetMedium("ZDC_ZVOID");
 620   TGeoVolume *pQCLint = new TGeoVolume("QCLint",pIntTrousersC, medZDCvoid);
 621   pQCLint->SetLineColor(kTeal);
 622   pQCLint->SetVisLeaves(kTRUE);
 623   pQCLext->AddNode(pQCLint, 1);
 624
 625   zd1 += 90.1;
 626   Double_t offset = 0.5;
 627   zd1 = zd1+offset;
 628
 629   //  second section : 2 tubes (ID = 54. OD = 58.)
 630   tubpar[0] = 5.4/2.;
 631   tubpar[1] = 5.8/2.;
 632   tubpar[2] = 40.0/2.;
 633   gMC->Gsvolu("QT12", "TUBE", idtmed[7], tubpar, 3);
 634   gMC->Gspos("QT12", 1, "ZDCC", -15.8/2., 0., -tubpar[2]-zd1, 0, "ONLY");
 635   gMC->Gspos("QT12", 2, "ZDCC",  15.8/2., 0., -tubpar[2]-zd1, 0, "ONLY");
 636   // Ch.debug
 637   //printf("    QT12 TUBE from z = %1.2f to z= %1.2f\n",-zd1,-2*tubpar[2]-zd1);
 638
 639   zd1 += 2.*tubpar[2];
 640
 641   // transition x2zdc to recombination chamber : skewed cone
 642   conpar[0] = (10.-0.2-offset)/2.;
 643   conpar[1] = 5.4/2.;
 644   conpar[2] = 5.8/2.;
 645   conpar[3] = 6.3/2.;
 646   conpar[4] = 7.0/2.;
 647   gMC->Gsvolu("QC10", "CONE", idtmed[7], conpar, 5);
 648   gMC->Gspos("QC10", 1, "ZDCC", -7.9-0.175, 0., -conpar[0]-0.1-zd1, irotpipe1, "ONLY");
 649   gMC->Gspos("QC10", 2, "ZDCC", 7.9+0.175, 0., -conpar[0]-0.1-zd1, irotpipe2, "ONLY");
 650   //printf("    QC10 CONE from z = %1.2f to z= %1.2f\n",-zd1,-2*conpar[0]-0.2-zd1);
 651
 652   zd1 += 2.*conpar[0]+0.2;
 653
 654   // 2 tubes (ID = 63 mm OD=70 mm)
 655   tubpar[0] = 6.3/2.;
 656   tubpar[1] = 7.0/2.;
 657   tubpar[2] = 639.8/2.;
 658   gMC->Gsvolu("QT13", "TUBE", idtmed[7], tubpar, 3);
 659   gMC->Gspos("QT13", 1, "ZDCC", -16.5/2., 0., -tubpar[2]-zd1, 0, "ONLY");
 660   gMC->Gspos("QT13", 2, "ZDCC",  16.5/2., 0., -tubpar[2]-zd1, 0, "ONLY");
 661   //printf("    QT13 TUBE from z = %1.2f to z= %1.2f\n",-zd1,-2*tubpar[2]-zd1);
 662
 663   zd1 += 2.*tubpar[2];
 664   //printf("    END OF SIDE C BEAM PIPE DEFINITION @ z = %f\n",-zd1);
 665
 666
 667   // -- Luminometer (Cu box) in front of ZN - side C
 668   boxpar[0] = 8.0/2.;
 669   boxpar[1] = 8.0/2.;
 670   boxpar[2] = fLumiLength/2.;
 671   gMC->Gsvolu("QLUC", "BOX ", idtmed[9], boxpar, 3);
 672   gMC->Gspos("QLUC", 1, "ZDCC", 0., 0.,  fPosZNC[2]+66.+boxpar[2], 0, "ONLY");
 673   //printf("    QLUC LUMINOMETER from z = %1.2f to z= %1.2f\n",  fPosZNC[2]+66., fPosZNC[2]+66.+2*boxpar[2]);
 674
 675   // --  END OF BEAM PIPE VOLUME DEFINITION FOR SIDE C (RB26 SIDE)
 676   // ----------------------------------------------------------------
 677
 678   ////////////////////////////////////////////////////////////////
 679   //                                                            //
 680   //                SIDE A - RB24                               //
 681   //                                                            //
 682   ///////////////////////////////////////////////////////////////
 683
 684   // Rotation Matrices definition
 685   Int_t irotpipe3, irotpipe4, irotpipe5;
 686   //-- rotation matrices for the tilted cone after the TDI to recenter vacuum chamber
 687   gMC->Matrix(irotpipe3,90.-1.8934,0.,90.,90.,1.8934,180.);
 688   //-- rotation matrices for the tilted tube before and after the TDI
 689   gMC->Matrix(irotpipe4,90.-3.8,0.,90.,90.,3.8,180.);
 690   //-- rotation matrix for the tilted cone after the TDI
 691   gMC->Matrix(irotpipe5,90.+9.8,0.,90.,90.,9.8,0.);
 692
 693   // -- Mother of the ZDCs (Vacuum PCON)
 694   zd2 = 1910.22;// zd2 initial value
 695
 696   conpar[0] = 0.;
 697   conpar[1] = 360.;
 698   conpar[2] = 2.;
 699   conpar[3] = zd2;
 700   conpar[4] = 0.;
 701   conpar[5] = 55.;
 702   conpar[6] = 13500.;
 703   conpar[7] = 0.;
 704   conpar[8] = 55.;
 705   gMC->Gsvolu("ZDCA", "PCON", idtmed[10], conpar, 9);
 706   gMC->Gspos("ZDCA", 1, "ALIC", 0., 0., 0., 0, "ONLY");
 707
 708   // To avoid overlaps 1 micron are left between certain volumes!
 709   Double_t dxNoOverlap = 0.0;
 710   //zd2 += dxNoOverlap;
 711
 712   // BEAM PIPE from 19.10 m to inner triplet beginning (22.965 m)
 713   tubpar[0] = 6.0/2.;
 714   tubpar[1] = 6.4/2.;
 715   tubpar[2] = 386.28/2. - dxNoOverlap;
 716   gMC->Gsvolu("QA01", "TUBE", idtmed[7], tubpar, 3);
 717   gMC->Gspos("QA01", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
 718   // Ch.debug
 719   //printf("    QA01 TUBE centred in %f from z = %1.2f to z= %1.2f (IT begin)\n",tubpar[2]+zd2,zd2,2*tubpar[2]+zd2);
 720
 721   zd2 += 2.*tubpar[2];
 722
 723   // -- FIRST SECTION OF THE BEAM PIPE (from beginning of inner triplet to
 724   //    beginning of D1)
 725   tubpar[0] = 6.3/2.;
 726   tubpar[1] = 6.7/2.;
 727   tubpar[2] = 3541.8/2. - dxNoOverlap;
 728   gMC->Gsvolu("QA02", "TUBE", idtmed[7], tubpar, 3);
 729   gMC->Gspos("QA02", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
 730   // Ch.debug
 731   //printf("    QA02 TUBE from z = %1.2f to z= %1.2f (D1 begin)\n",zd2,2*tubpar[2]+zd2);
 732
 733   zd2 += 2.*tubpar[2];
 734
 735
 736   // -- SECOND SECTION OF THE BEAM PIPE (from the beginning of D1 to the beginning of D2)
 737   //
 738   //  FROM (MAGNETIC) BEGINNING OF D1 TO THE (MAGNETIC) END OF D1 + 126.5 cm
 739   //  CYLINDRICAL PIPE of diameter increasing from 6.75 cm up to 8.0 cm
 740   //  from magnetic end :
 741   //  1) 80.1 cm still with ID = 6.75 radial beam screen
 742   //  2) 2.5 cm conical section from ID = 6.75 to ID = 8.0 cm
 743   //  3) 43.9 cm straight section (tube) with ID = 8.0 cm
 744
 745   tubpar[0] = 6.75/2.;
 746   tubpar[1] = 7.15/2.;
 747   tubpar[2] = (945.0+80.1)/2.;
 748   gMC->Gsvolu("QA03", "TUBE", idtmed[7], tubpar, 3);
 749   gMC->Gspos("QA03", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
 750   // Ch.debug
 751   //printf("    QA03 TUBE from z = %1.2f to z= %1.2f\n",zd2,2*tubpar[2]+zd2);
 752
 753   zd2 += 2.*tubpar[2];
 754
 755   // Transition Cone from ID=67.5 mm  to ID=80 mm
 756   conpar[0] = 2.5/2.;
 757   conpar[1] = 6.75/2.;
 758   conpar[2] = 7.15/2.;
 759   conpar[3] = 8.0/2.;
 760   conpar[4] = 8.4/2.;
 761   gMC->Gsvolu("QA04", "CONE", idtmed[7], conpar, 5);
 762   gMC->Gspos("QA04", 1, "ZDCA", 0., 0., conpar[0]+zd2, 0, "ONLY");
 763   //printf("    QA04 CONE from z = %1.2f to z= %1.2f\n",zd2,2*conpar[0]+zd2);
 764
 765   zd2 += 2.*conpar[0];
 766
 767   tubpar[0] = 8.0/2.;
 768   tubpar[1] = 8.4/2.;
 769   tubpar[2] = (43.9+20.+28.5+28.5)/2.;
 770   gMC->Gsvolu("QA05", "TUBE", idtmed[7], tubpar, 3);
 771   gMC->Gspos("QA05", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
 772   // Ch.debug
 773   //printf("    QA05 TUBE from z = %1.2f to z= %1.2f\n",zd2,2*tubpar[2]+zd2);
 774
 775   zd2 += 2.*tubpar[2];
 776
 777   // Second section of VAEHI (transition cone from ID=80mm to ID=98mm)
 778   conpar[0] = 4.0/2.;
 779   conpar[1] = 8.0/2.;
 780   conpar[2] = 8.4/2.;
 781   conpar[3] = 9.8/2.;
 782   conpar[4] = 10.2/2.;
 783   gMC->Gsvolu("QAV1", "CONE", idtmed[7], conpar, 5);
 784   gMC->Gspos("QAV1", 1, "ZDCA", 0., 0., conpar[0]+zd2, 0, "ONLY");
 785   //printf("    QAV1 CONE from z = %1.2f to z= %1.2f\n",zd2,2*conpar[0]+zd2);
 786
 787   zd2 += 2.*conpar[0];
 788
 789   //Third section of VAEHI (transition cone from ID=98mm to ID=90mm)
 790   conpar[0] = 1.0/2.;
 791   conpar[1] = 9.8/2.;
 792   conpar[2] = 10.2/2.;
 793   conpar[3] = 9.0/2.;
 794   conpar[4] = 9.4/2.;
 795   gMC->Gsvolu("QAV2", "CONE", idtmed[7], conpar, 5);
 796   gMC->Gspos("QAV2", 1, "ZDCA", 0., 0., conpar[0]+zd2, 0, "ONLY");
 797   //printf("    QAV2 CONE from z = %1.2f to z= %1.2f\n",zd2,2*conpar[0]+zd2);
 798
 799   zd2 += 2.*conpar[0];
 800
 801   // Fourth section of VAEHI (tube ID=90mm)
 802   tubpar[0] = 9.0/2.;
 803   tubpar[1] = 9.4/2.;
 804   tubpar[2] = 31.0/2.;
 805   gMC->Gsvolu("QAV3", "TUBE", idtmed[7], tubpar, 3);
 806   gMC->Gspos("QAV3", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
 807   // Ch.debug
 808   //printf("    QAV3 TUBE from z = %1.2f to z= %1.2f\n",zd2,2*tubpar[2]+zd2);
 809
 810   zd2 += 2.*tubpar[2];
 811
 812   //---------------------------- TCDD beginning ----------------------------------
 813   // space for the insertion of the collimator TCDD (2 m)
 814   // TCDD ZONE - 1st volume
 815   conpar[0] = 1.3/2.;
 816   conpar[1] = 9.0/2.;
 817   conpar[2] = 13.0/2.;
 818   conpar[3] = 9.6/2.;
 819   conpar[4] = 13.0/2.;
 820   gMC->Gsvolu("Q01T", "CONE", idtmed[7], conpar, 5);
 821   gMC->Gspos("Q01T", 1, "ZDCA", 0., 0., conpar[0]+zd2, 0, "ONLY");
 822   //printf("    Q01T CONE from z = %1.2f to z= %1.2f\n",zd2,2*conpar[0]+zd2);
 823
 824   zd2 += 2.*conpar[0];
 825
 826   // TCDD ZONE - 2nd volume
 827   tubpar[0] = 9.6/2.;
 828   tubpar[1] = 10.0/2.;
 829   tubpar[2] = 1.0/2.;
 830   gMC->Gsvolu("Q02T", "TUBE", idtmed[7], tubpar, 3);
 831   gMC->Gspos("Q02T", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
 832   // Ch.debug
 833   //printf("    Q02T TUBE from z = %1.2f to z= %1.2f\n",zd2,2*tubpar[2]+zd2);
 834
 835   zd2 += 2.*tubpar[2];
 836
 837   // TCDD ZONE - third volume
 838   conpar[0] = 9.04/2.;
 839   conpar[1] = 9.6/2.;
 840   conpar[2] = 10.0/2.;
 841   conpar[3] = 13.8/2.;
 842   conpar[4] = 14.2/2.;
 843   gMC->Gsvolu("Q03T", "CONE", idtmed[7], conpar, 5);
 844   gMC->Gspos("Q03T", 1, "ZDCA", 0., 0., conpar[0]+zd2, 0, "ONLY");
 845   //printf("    Q03T CONE from z = %1.2f to z= %1.2f\n",zd2,2*conpar[0]+zd2);
 846
 847   zd2 += 2.*conpar[0];
 848
 849   // TCDD ZONE - 4th volume
 850   tubpar[0] = 13.8/2.;
 851   tubpar[1] = 14.2/2.;
 852   tubpar[2] = 38.6/2.;
 853   gMC->Gsvolu("Q04T", "TUBE", idtmed[7], tubpar, 3);
 854   gMC->Gspos("Q04T", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
 855   // Ch.debug
 856   //printf("    Q04T TUBE from z = %1.2f to z= %1.2f\n",zd2,2*tubpar[2]+zd2);
 857
 858   zd2 += 2.*tubpar[2];
 859
 860   // TCDD ZONE - 5th volume
 861   tubpar[0] = 21.0/2.;
 862   tubpar[1] = 21.4/2.;
 863   tubpar[2] = 100.12/2.;
 864   gMC->Gsvolu("Q05T", "TUBE", idtmed[7], tubpar, 3);
 865   gMC->Gspos("Q05T", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
 866   // Ch.debug
 867   //printf("    Q05T TUBE from z = %1.2f to z= %1.2f\n",zd2,2*tubpar[2]+zd2);
 868
 869   zd2 += 2.*tubpar[2];
 870
 871   // TCDD ZONE - 6th volume
 872   tubpar[0] = 13.8/2.;
 873   tubpar[1] = 14.2/2.;
 874   tubpar[2] = 38.6/2.;
 875   gMC->Gsvolu("Q06T", "TUBE", idtmed[7], tubpar, 3);
 876   gMC->Gspos("Q06T", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
 877   // Ch.debug
 878   //printf("    Q06T TUBE from z = %1.2f to z= %1.2f\n",zd2,2*tubpar[2]+zd2);
 879
 880   zd2 += 2.*tubpar[2];
 881
 882   // TCDD ZONE - 7th volume
 883   conpar[0] = 11.34/2.;
 884   conpar[1] = 13.8/2.;
 885   conpar[2] = 14.2/2.;
 886   conpar[3] = 18.0/2.;
 887   conpar[4] = 18.4/2.;
 888   gMC->Gsvolu("Q07T", "CONE", idtmed[7], conpar, 5);
 889   gMC->Gspos("Q07T", 1, "ZDCA", 0., 0., conpar[0]+zd2, 0, "ONLY");
 890   //printf("    Q07T CONE from z = %1.2f to z= %1.2f\n",zd2,2*conpar[0]+zd2);
 891
 892   zd2 += 2.*conpar[0];
 893
 894   // Upper section : one single phi segment of a tube
 895   //  5 parameters for tubs: inner radius = 0.,
 896   //    outer radius = 7. cm, half length = 50 cm
 897   //    phi1 = 0., phi2 = 180.
 898   tubspar[0] = 0.0/2.;
 899   tubspar[1] = 14.0/2.;
 900   tubspar[2] = 100.0/2.;
 901   tubspar[3] = 0.;
 902   tubspar[4] = 180.;
 903   gMC->Gsvolu("Q08T", "TUBS", idtmed[7], tubspar, 5);
 904   // Ch.debug
 905   //printf("    upper part : one single phi segment of a tube (Q08T)\n");
 906
 907   // rectangular beam pipe inside TCDD upper section (Vacuum)
 908   boxpar[0] = 7.0/2.;
 909   boxpar[1] = 2.2/2.;
 910   boxpar[2] = 100./2.;
 911   gMC->Gsvolu("Q09T", "BOX ", idtmed[10], boxpar, 3);
 912   // positioning vacuum box in the upper section of TCDD
 913   gMC->Gspos("Q09T", 1, "Q08T", 0., 1.1,  0., 0, "ONLY");
 914
 915   // lower section : one single phi segment of a tube
 916   tubspar[0] = 0.0/2.;
 917   tubspar[1] = 14.0/2.;
 918   tubspar[2] = 100.0/2.;
 919   tubspar[3] = 180.;
 920   tubspar[4] = 360.;
 921   gMC->Gsvolu("Q10T", "TUBS", idtmed[7], tubspar, 5);
 922   // rectangular beam pipe inside TCDD lower section (Vacuum)
 923   boxpar[0] = 7.0/2.;
 924   boxpar[1] = 2.2/2.;
 925   boxpar[2] = 100./2.;
 926   gMC->Gsvolu("Q11T", "BOX ", idtmed[10], boxpar, 3);
 927   // positioning vacuum box in the lower section of TCDD
 928   gMC->Gspos("Q11T", 1, "Q10T", 0., -1.1,  0., 0, "ONLY");
 929
 930   // positioning  TCDD elements in ZDCA, (inside TCDD volume)
 931   gMC->Gspos("Q08T", 1, "ZDCA", 0., 2., -100.+zd2, 0, "ONLY");
 932   gMC->Gspos("Q10T", 1, "ZDCA", 0., -2., -100.+zd2, 0, "ONLY");
 933
 934   // RF screen
 935   boxpar[0] = 0.2/2.;
 936   boxpar[1] = 4.0/2.;
 937   boxpar[2] = 100./2.;
 938   gMC->Gsvolu("Q12T", "BOX ", idtmed[7], boxpar, 3);
 939   // positioning RF screen at both sides of TCDD
 940   gMC->Gspos("Q12T", 1, "ZDCA", tubspar[1]+boxpar[0], 0., -100.+zd2, 0, "ONLY");
 941   gMC->Gspos("Q12T", 2, "ZDCA", -tubspar[1]-boxpar[0], 0., -100.+zd2, 0, "ONLY");
 942   //---------------------------- TCDD end ---------------------------------------
 943
 944   // The following elliptical tube 180 mm x 70 mm
 945   // (obtained positioning the void QA09 in QA08)
 946   // represents VMTSA (780 mm) + space reserved to the TCTVB (1480 mm)+
 947   //            VMTSA (780 mm) + first part of VCTCP (93 mm)
 948
 949   tubpar[0] = 18.4/2.;
 950   tubpar[1] = 7.4/2.;
 951   tubpar[2] = 313.3/2.;
 952 //  gMC->Gsvolu("QA06", "ELTU", idtmed[7], tubpar, 3);
 953   // temporary replace with a scaled tube (AG)
 954   TGeoTube *tubeQA06 = new TGeoTube(0.,tubpar[0],tubpar[2]);
 955   TGeoScale *scaleQA06 = new TGeoScale(1., tubpar[1]/tubpar[0], 1.);
 956   TGeoScaledShape *sshapeQA06 = new TGeoScaledShape(tubeQA06, scaleQA06);
 957   new TGeoVolume("QA06", sshapeQA06, gGeoManager->GetMedium(idtmed[7]));
 958   //printf("    QA06 TUBE from z = %1.2f to z= %1.2f\n",zd2,2*tubpar[2]+zd2);
 959
 960   tubpar[0] = 18.0/2.;
 961   tubpar[1] = 7.0/2.;
 962   tubpar[2] = 313.3/2.;
 963 //  gMC->Gsvolu("QA07", "ELTU", idtmed[10], tubpar, 3);
 964   // temporary replace with a scaled tube (AG)
 965   TGeoTube *tubeQA07 = new TGeoTube(0.,tubpar[0],tubpar[2]);
 966   TGeoScale *scaleQA07 = new TGeoScale(1., tubpar[1]/tubpar[0], 1.);
 967   TGeoScaledShape *sshapeQA07 = new TGeoScaledShape(tubeQA07, scaleQA07);
 968   new TGeoVolume("QA07", sshapeQA07, gGeoManager->GetMedium(idtmed[7]));
 969   //printf("    QA07 TUBE from z = %1.2f to z= %1.2f\n",zd2,2*tubpar[2]+zd2);
 970   gMC->Gspos("QA06", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
 971   gMC->Gspos("QA07", 1, "QA06", 0., 0., 0., 0, "ONLY");
 972
 973   // Vertical collimator jaws (defined ONLY if fVCollAperture<3.5!)
 974   if(fVCollSideAAperture<3.5){
 975     boxpar[0] = 5.4/2.;
 976     boxpar[1] = (3.5-fVCollSideAAperture-fVCollSideACentreY-0.7)/2.;
 977     if(boxpar[1]<0.) boxpar[1]=0.;
 978     boxpar[2] = 124.4/2.;
 979     gMC->Gsvolu("QCVA" , "BOX ", idtmed[13], boxpar, 3);
 980     gMC->Gspos("QCVA", 1, "QA07", -boxpar[0], fVCollSideAAperture+fVCollSideACentreY+boxpar[1], -313.3/2.+78.+148./2., 0, "ONLY");
 981     gMC->Gspos("QCVA", 2, "QA07", -boxpar[0], -fVCollSideAAperture+fVCollSideACentreY-boxpar[1], -313.3/2.+78.+148./2., 0, "ONLY");
 982   }
 983
 984   zd2 += 2.*tubpar[2];
 985
 986   // VCTCP second part: transition cone from ID=180 to ID=212.7
 987   conpar[0] = 31.5/2.;
 988   conpar[1] = 18.0/2.;
 989   conpar[2] = 18.6/2.;
 990   conpar[3] = 21.27/2.;
 991   conpar[4] = 21.87/2.;
 992   gMC->Gsvolu("QA08", "CONE", idtmed[7], conpar, 5);
 993   gMC->Gspos("QA08", 1, "ZDCA", 0., 0., conpar[0]+zd2, 0, "ONLY");
 994   // Ch.debug
 995   //printf("    QA08 CONE from z = %Third part of VCTCR: tube (ID=196 mm) f to z = %f\n",zd2,2*conpar[0]+zd2);
 996
 997   zd2 += 2.*conpar[0];
 998
 999   // Tube ID 212.7 mm
1000   // Represents VCTCP third part (92 mm) + VCDWB (765 mm) + VMBGA (400 mm) +
1001   //            VCDWE (300 mm) + VMBGA (400 mm)
1002   tubpar[0] = 21.27/2.;
1003   tubpar[1] = 21.87/2.;
1004   tubpar[2] = 195.7/2.;
1005   gMC->Gsvolu("QA09", "TUBE", idtmed[7], tubpar, 3);
1006   gMC->Gspos("QA09", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
1007   //printf("    QA09 TUBE from z = %1.2f to z= %1.2f\n",zd2,2*tubpar[2]+zd2);
1008
1009   zd2 += 2.*tubpar[2];
1010
1011   // skewed transition piece (ID=212.7 mm to 332 mm) (before TDI)
1012   conpar[0] = (50.0-0.73-1.13)/2.;
1013   conpar[1] = 21.27/2.;
1014   conpar[2] = 21.87/2.;
1015   conpar[3] = 33.2/2.;
1016   conpar[4] = 33.8/2.;
1017   gMC->Gsvolu("QA10", "CONE", idtmed[7], conpar, 5);
1018   gMC->Gspos("QA10", 1, "ZDCA", -1.66, 0., conpar[0]+0.73+zd2, irotpipe4, "ONLY");
1019   // Ch.debug
1020   //printf("    QA10 skewed CONE from z = %1.2f to z= %1.2f\n",zd2,2*conpar[0]+0.73+1.13+zd2);
1021
1022   zd2 += 2.*conpar[0]+0.73+1.13;
1023
1024   // Vacuum chamber containing TDI
1025   tubpar[0] = 0.;
1026   tubpar[1] = 54.6/2.;
1027   tubpar[2] = 540.0/2.;
1028   gMC->Gsvolu("Q13TM", "TUBE", idtmed[10], tubpar, 3);
1029   gMC->Gspos("Q13TM", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
1030   tubpar[0] = 54.0/2.;
1031   tubpar[1] = 54.6/2.;
1032   tubpar[2] = 540.0/2.;
1033   gMC->Gsvolu("Q13T", "TUBE", idtmed[7], tubpar, 3);
1034   gMC->Gspos("Q13T", 1, "Q13TM", 0., 0., 0., 0, "ONLY");
1035   // Ch.debug
1036   //printf("    Q13T TUBE from z = %1.2f to z= %1.2f\n",zd2,2*tubpar[2]+zd2);
1037
1038   zd2 += 2.*tubpar[2];
1039
1040   //---------------- INSERT TDI INSIDE Q13T -----------------------------------
1041   boxpar[0] = 11.0/2.;
1042   boxpar[1] = 9.0/2.;
1043   boxpar[2] = 540.0/2.;
1044   gMC->Gsvolu("QTD1", "BOX ", idtmed[7], boxpar, 3);
1045   gMC->Gspos("QTD1", 1, "Q13TM", -3.8, 10.5,  0., 0, "ONLY");
1046   boxpar[0] = 11.0/2.;
1047   boxpar[1] = 9.0/2.;
1048   boxpar[2] = 540.0/2.;
1049   gMC->Gsvolu("QTD2", "BOX ", idtmed[7], boxpar, 3);
1050   gMC->Gspos("QTD2", 1, "Q13TM", -3.8, -10.5,  0., 0, "ONLY");
1051   boxpar[0] = 5.1/2.;
1052   boxpar[1] = 0.2/2.;
1053   boxpar[2] = 540.0/2.;
1054   gMC->Gsvolu("QTD3", "BOX ", idtmed[7], boxpar, 3);
1055   gMC->Gspos("QTD3", 1, "Q13TM", -3.8+5.5+boxpar[0], 6.1,  0., 0, "ONLY");
1056   gMC->Gspos("QTD3", 2, "Q13TM", -3.8+5.5+boxpar[0], -6.1,  0., 0, "ONLY");
1057   gMC->Gspos("QTD3", 3, "Q13TM", -3.8-5.5-boxpar[0], 6.1,  0., 0, "ONLY");
1058   gMC->Gspos("QTD3", 4, "Q13TM", -3.8-5.5-boxpar[0], -6.1,  0., 0, "ONLY");
1059   //
1060   tubspar[0] = 12.0/2.;
1061   tubspar[1] = 12.4/2.;
1062   tubspar[2] = 540.0/2.;
1063   tubspar[3] = 90.;
1064   tubspar[4] = 270.;
1065   gMC->Gsvolu("QTD4", "TUBS", idtmed[7], tubspar, 5);
1066   gMC->Gspos("QTD4", 1, "Q13TM", -3.8-10.6, 0.,  0., 0, "ONLY");
1067   tubspar[0] = 12.0/2.;
1068   tubspar[1] = 12.4/2.;
1069   tubspar[2] = 540.0/2.;
1070   tubspar[3] = -90.;
1071   tubspar[4] = 90.;
1072   gMC->Gsvolu("QTD5", "TUBS", idtmed[7], tubspar, 5);
1073   gMC->Gspos("QTD5", 1, "Q13TM", -3.8+10.6, 0.,  0., 0, "ONLY");
1074   //---------------- END DEFINING TDI INSIDE Q13T -------------------------------
1075
1076   // VCTCG skewed transition piece (ID=332 mm to 212.7 mm) (after TDI)
1077   conpar[0] = (50.0-2.92-1.89)/2.;
1078   conpar[1] = 33.2/2.;
1079   conpar[2] = 33.8/2.;
1080   conpar[3] = 21.27/2.;
1081   conpar[4] = 21.87/2.;
1082   gMC->Gsvolu("QA11", "CONE", idtmed[7], conpar, 5);
1083   gMC->Gspos("QA11", 1, "ZDCA", 4.32-3.8, 0., conpar[0]+2.92+zd2, irotpipe5, "ONLY");
1084   // Ch.debug
1085   //printf("    QA11 skewed CONE from z = %f to z =%f\n",zd2,2*conpar[0]+2.92+1.89+zd2);
1086
1087   zd2 += 2.*conpar[0]+2.92+1.89;
1088
1089   // The following tube ID 212.7 mm
1090   // represents VMBGA (400 mm) + VCDWE (300 mm) + VMBGA (400 mm) +
1091   //            BTVTS (600 mm) + VMLGB (400 mm)
1092   tubpar[0] = 21.27/2.;
1093   tubpar[1] = 21.87/2.;
1094   tubpar[2] = 210.0/2.;
1095   gMC->Gsvolu("QA12", "TUBE", idtmed[7], tubpar, 3);
1096   gMC->Gspos("QA12", 1, "ZDCA", 4., 0., tubpar[2]+zd2, 0, "ONLY");
1097   // Ch.debug
1098   //printf("    QA12 TUBE from z = %1.2f to z= %1.2f\n",zd2,2*tubpar[2]+zd2);
1099
1100   zd2 += 2.*tubpar[2];
1101
1102   // First part of VCTCC
1103   // skewed transition cone from ID=212.7 mm to ID=797 mm
1104   conpar[0] = (121.0-0.37-1.35)/2.;
1105   conpar[1] = 21.27/2.;
1106   conpar[2] = 21.87/2.;
1107   conpar[3] = 79.7/2.;
1108   conpar[4] = 81.3/2.;
1109   gMC->Gsvolu("QA13", "CONE", idtmed[7], conpar, 5);
1110   gMC->Gspos("QA13", 1, "ZDCA", 4.-2., 0., conpar[0]+0.37+zd2, irotpipe3, "ONLY");
1111   // Ch.debug
1112   //printf("    QA13 CONE from z = %1.2f to z= %1.2f\n",zd2,2*conpar[0]+0.37+1.35+zd2);
1113
1114   zd2 += 2.*conpar[0]+0.37+1.35;
1115
1116   // The following tube ID 797 mm  --- (volume QA16)
1117   // represents the second part of VCTCC (4272 mm) +
1118   //            4 x VCDGA (4 x 4272 mm) +
1119   //            the first part of VCTCR (850 mm)
1120   tubpar[0] = 79.7/2.;
1121   tubpar[1] = 81.3/2.;
1122   tubpar[2] = 2221./2.;
1123   gMC->Gsvolu("QA14", "TUBE", idtmed[7], tubpar, 3);
1124   gMC->Gspos("QA14", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
1125   // Ch.debug
1126   //printf("    QA14 TUBE from z = %1.2f to z= %1.2f\n",zd2,2*tubpar[2]+zd2);
1127
1128   zd2 += 2.*tubpar[2];
1129
1130   // Second part of VCTCR
1131   // Transition from ID=797 mm to ID=196 mm:
1132   // in order to simulate the thin window opened in the transition cone
1133   // we divide the transition cone in three cones:
1134   // (1) 8 mm thick (2) 3 mm thick (3) the third 8 mm thick
1135
1136   // (1) 8 mm thick
1137   conpar[0] = 9.09/2.; // 15 degree
1138   conpar[1] = 79.7/2.;
1139   conpar[2] = 81.3/2.; // thickness 8 mm
1140   conpar[3] = 74.82868/2.;
1141   conpar[4] = 76.42868/2.; // thickness 8 mm
1142   gMC->Gsvolu("QA15", "CONE", idtmed[7], conpar, 5);
1143   gMC->Gspos("QA15", 1, "ZDCA", 0., 0., conpar[0]+zd2, 0, "ONLY");
1144   //printf("    QA15 CONE from z = %1.2f to z= %1.2f\n",zd2,2*conpar[0]+zd2);
1145
1146   zd2 += 2.*conpar[0];
1147
1148   // (2) 3 mm thick
1149   conpar[0] = 96.2/2.; // 15 degree
1150   conpar[1] = 74.82868/2.;
1151   conpar[2] = 75.42868/2.; // thickness 3 mm
1152   conpar[3] = 23.19588/2.;
1153   conpar[4] = 23.79588/2.; // thickness 3 mm
1154   gMC->Gsvolu("QA16", "CONE", idtmed[7], conpar, 5);
1155   gMC->Gspos("QA16", 1, "ZDCA", 0., 0., conpar[0]+zd2, 0, "ONLY");
1156   //printf("    QA16 CONE from z = %1.2f to z= %1.2f\n",zd2,2*conpar[0]+zd2);
1157
1158   zd2 += 2.*conpar[0];
1159
1160   // (3) 8 mm thick
1161   conpar[0] = 6.71/2.; // 15 degree
1162   conpar[1] = 23.19588/2.;
1163   conpar[2] = 24.79588/2.;// thickness 8 mm
1164   conpar[3] = 19.6/2.;
1165   conpar[4] = 21.2/2.;// thickness 8 mm
1166   gMC->Gsvolu("QA17", "CONE", idtmed[7], conpar, 5);
1167   gMC->Gspos("QA17", 1, "ZDCA", 0., 0., conpar[0]+zd2, 0, "ONLY");
1168   //printf("    QA19 CONE from z = %1.2f to z= %1.2f\n",zd2,2*conpar[0]+zd2);
1169
1170   zd2 += 2.*conpar[0];
1171
1172   // Third part of VCTCR: tube (ID=196 mm)
1173   tubpar[0] = 19.6/2.;
1174   tubpar[1] = 21.2/2.;
1175   tubpar[2] = 9.55/2.;
1176   gMC->Gsvolu("QA18", "TUBE", idtmed[7], tubpar, 3);
1177   gMC->Gspos("QA18", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
1178   // Ch.debug
1179   //printf("    QA18 TUBE from z = %1.2f to z= %1.2f\n",zd2,2*tubpar[2]+zd2);
1180
1181   zd2 += 2.*tubpar[2];
1182
1183   // Flange (ID=196 mm) (last part of VCTCR and first part of VMZAR)
1184   tubpar[0] = 19.6/2.;
1185   tubpar[1] = 25.3/2.;
1186   tubpar[2] = 4.9/2.;
1187   gMC->Gsvolu("QF01", "TUBE", idtmed[7], tubpar, 3);
1188   gMC->Gspos("QF01", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
1189   // Ch.debug
1190   //printf("    QF01  TUBE from z = %1.2f to z= %1.2f\n",zd2,2*tubpar[2]+zd2);
1191
1192   zd2 += 2.*tubpar[2];
1193
1194   // VMZAR (5 volumes)
1195   tubpar[0] = 20.2/2.;
1196   tubpar[1] = 20.6/2.;
1197   tubpar[2] = 2.15/2.;
1198   gMC->Gsvolu("QA19", "TUBE", idtmed[7], tubpar, 3);
1199   gMC->Gspos("QA19", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
1200   // Ch.debug
1201   //printf("    QA19  TUBE from z = %1.2f to z= %1.2f\n",zd2,2*tubpar[2]+zd2);
1202
1203   zd2 += 2.*tubpar[2];
1204
1205   conpar[0] = 6.9/2.;
1206   conpar[1] = 20.2/2.;
1207   conpar[2] = 20.6/2.;
1208   conpar[3] = 23.9/2.;
1209   conpar[4] = 24.3/2.;
1210   gMC->Gsvolu("QA20", "CONE", idtmed[7], conpar, 5);
1211   gMC->Gspos("QA20", 1, "ZDCA", 0., 0., conpar[0]+zd2, 0, "ONLY");
1212   // Ch.debug
1213   //printf("    QA20 CONE from z = %1.2f to z= %1.2f\n",zd2,2*conpar[0]+zd2);
1214
1215   zd2 += 2.*conpar[0];
1216
1217   tubpar[0] = 23.9/2.;
1218   tubpar[1] = 25.5/2.;
1219   tubpar[2] = 17.0/2.;
1220   gMC->Gsvolu("QA21", "TUBE", idtmed[7], tubpar, 3);
1221   gMC->Gspos("QA21", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
1222   // Ch.debug
1223   //printf("    QA21  TUBE from z = %1.2f to z= %1.2f\n",zd2,2*tubpar[2]+zd2);
1224
1225   zd2 += 2.*tubpar[2];
1226
1227   conpar[0] = 6.9/2.;
1228   conpar[1] = 23.9/2.;
1229   conpar[2] = 24.3/2.;
1230   conpar[3] = 20.2/2.;
1231   conpar[4] = 20.6/2.;
1232   gMC->Gsvolu("QA22", "CONE", idtmed[7], conpar, 5);
1233   gMC->Gspos("QA22", 1, "ZDCA", 0., 0., conpar[0]+zd2, 0, "ONLY");
1234   // Ch.debug
1235   //printf("    QA22 CONE from z = %1.2f to z= %1.2f\n",zd2,2*conpar[0]+zd2);
1236
1237   zd2 += 2.*conpar[0];
1238
1239   tubpar[0] = 20.2/2.;
1240   tubpar[1] = 20.6/2.;
1241   tubpar[2] = 2.15/2.;
1242   gMC->Gsvolu("QA23", "TUBE", idtmed[7], tubpar, 3);
1243   gMC->Gspos("QA23", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
1244   // Ch.debug
1245   //printf("    QA23  TUBE from z = %1.2f to z= %1.2f\n",zd2,2*tubpar[2]+zd2);
1246
1247   zd2 += 2.*tubpar[2];
1248
1249   // Flange (ID=196 mm)(last part of VMZAR and first part of VCTYD)
1250   tubpar[0] = 19.6/2.;
1251   tubpar[1] = 25.3/2.;
1252   tubpar[2] = 4.9/2.;
1253   gMC->Gsvolu("QF02", "TUBE", idtmed[7], tubpar, 3);
1254   gMC->Gspos("QF02", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
1255   // Ch.debug
1256   //printf("    QF02 TUBE from z = %1.2f to z= %1.2f\n",zd2,2*tubpar[2]+zd2);
1257
1258   zd2 += 2.*tubpar[2];
1259
1260   // simulation of the trousers (VCTYB)
1261   tubpar[0] = 19.6/2.;
1262   tubpar[1] = 20.0/2.;
1263   tubpar[2] = 3.9/2.;
1264   gMC->Gsvolu("QA24", "TUBE", idtmed[7], tubpar, 3);
1265   gMC->Gspos("QA24", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
1266   // Ch.debug
1267   //printf("    QA24  TUBE from z = %1.2f to z= %1.2f\n",zd2,2*tubpar[2]+zd2);
1268
1269   zd2 += 2.*tubpar[2];
1270
1271   // transition cone from ID=196. to ID=216.6
1272   conpar[0] = 32.55/2.;
1273   conpar[1] = 19.6/2.;
1274   conpar[2] = 20.0/2.;
1275   conpar[3] = 21.66/2.;
1276   conpar[4] = 22.06/2.;
1277   gMC->Gsvolu("QA25", "CONE", idtmed[7], conpar, 5);
1278   gMC->Gspos("QA25", 1, "ZDCA", 0., 0., conpar[0]+zd2, 0, "ONLY");
1279   // Ch.debug
1280   //printf("    QA25 CONE from z = %1.2f to z= %1.2f\n",zd2,2*conpar[0]+zd2);
1281
1282   zd2 += 2.*conpar[0];
1283
1284   // tube
1285   tubpar[0] = 21.66/2.;
1286   tubpar[1] = 22.06/2.;
1287   tubpar[2] = 28.6/2.;
1288   gMC->Gsvolu("QA26", "TUBE", idtmed[7], tubpar, 3);
1289   gMC->Gspos("QA26", 1, "ZDCA", 0., 0., tubpar[2]+zd2, 0, "ONLY");
1290   // Ch.debug
1291   //printf("    QA26  TUBE from z = %1.2f to z= %1.2f\n",zd2,2*tubpar[2]+zd2);
1292
1293   zd2 += 2.*tubpar[2];
1294
1295   // --------------------------------------------------------
1296   // RECOMBINATION CHAMBER IMPLEMENTED USING TGeo CLASSES!!!!
1297   // author: Chiara (June 2008)
1298   // --------------------------------------------------------
1299   // TRANSFORMATION MATRICES
1300   // Combi transformation:
1301   dx = -3.970000;
1302   dy = 0.000000;
1303   dz = 0.0;
1304   // Rotation:
1305   thx = 84.989100;   phx = 0.000000;
1306   thy = 90.000000;   phy = 90.000000;
1307   thz = 5.010900;    phz = 180.000000;
1308   TGeoRotation *rotMatrix1 = new TGeoRotation("",thx,phx,thy,phy,thz,phz);
1309   // Combi transformation:
1310   dx = -3.970000;
1311   dy = 0.000000;
1312   dz = 0.0;
1313   TGeoCombiTrans *rotMatrix2 = new TGeoCombiTrans("ZDC_c1", dx,dy,dz,rotMatrix1);
1314   rotMatrix2->RegisterYourself();
1315   // Combi transformation:
1316   dx = 3.970000;
1317   dy = 0.000000;
1318   dz = 0.0;
1319   // Rotation:
1320   thx = 95.010900;   phx = 0.000000;
1321   thy = 90.000000;   phy = 90.000000;
1322   thz = 5.010900;    phz = 0.000000;
1323   TGeoRotation *rotMatrix3 = new TGeoRotation("",thx,phx,thy,phy,thz,phz);
1324   TGeoCombiTrans *rotMatrix4 = new TGeoCombiTrans("ZDC_c2", dx,dy,dz,rotMatrix3);
1325   rotMatrix4->RegisterYourself();
1326
1327
1328   // VOLUMES DEFINITION
1329   // Volume: ZDCA
1330   TGeoVolume *pZDCA = gGeoManager->GetVolume("ZDCA");
1331
1332   conpar[0] = (90.1-0.95-0.26)/2.;
1333   conpar[1] = 0.0/2.;
1334   conpar[2] = 21.6/2.;
1335   conpar[3] = 0.0/2.;
1336   conpar[4] = 5.8/2.;
1337   new TGeoCone("QALext", conpar[0],conpar[1],conpar[2],conpar[3],conpar[4]);
1338
1339   conpar[0] = (90.1-0.95-0.26)/2.;
1340   conpar[1] = 0.0/2.;
1341   conpar[2] = 21.2/2.;
1342   conpar[3] = 0.0/2.;
1343   conpar[4] = 5.4/2.;
1344   new TGeoCone("QALint", conpar[0],conpar[1],conpar[2],conpar[3],conpar[4]);
1345
1346   // Outer trousers
1347   TGeoCompositeShape *pOutTrousers = new TGeoCompositeShape("outTrousers", "QALext:ZDC_c1+QALext:ZDC_c2");
1348
1349   // Volume: QALext
1350   //TGeoMedium *medZDCFe = gGeoManager->GetMedium("ZDC_ZIRON");
1351   TGeoVolume *pQALext = new TGeoVolume("QALext",pOutTrousers, medZDCFe);
1352   pQALext->SetLineColor(kBlue);
1353   pQALext->SetVisLeaves(kTRUE);
1354   //
1355   TGeoTranslation *tr1 = new TGeoTranslation(0., 0., (Double_t) conpar[0]+0.95+zd2);
1356   pZDCA->AddNode(pQALext, 1, tr1);
1357   // Inner trousers
1358   TGeoCompositeShape *pIntTrousers = new TGeoCompositeShape("intTrousers", "QALint:ZDC_c1+QALint:ZDC_c2");
1359   // Volume: QALint
1360   //TGeoMedium *medZDCvoid = gGeoManager->GetMedium("ZDC_ZVOID");
1361   TGeoVolume *pQALint = new TGeoVolume("QALint",pIntTrousers, medZDCvoid);
1362   pQALint->SetLineColor(kAzure);
1363   pQALint->SetVisLeaves(kTRUE);
1364   pQALext->AddNode(pQALint, 1);
1365
1366   zd2 += 90.1;
1367
1368   //  second section : 2 tubes (ID = 54. OD = 58.)
1369   tubpar[0] = 5.4/2.;
1370   tubpar[1] = 5.8/2.;
1371   tubpar[2] = 40.0/2.;
1372   gMC->Gsvolu("QA27", "TUBE", idtmed[7], tubpar, 3);
1373   gMC->Gspos("QA27", 1, "ZDCA", -15.8/2., 0., tubpar[2]+zd2, 0, "ONLY");
1374   gMC->Gspos("QA27", 2, "ZDCA",  15.8/2., 0., tubpar[2]+zd2, 0, "ONLY");
1375   // Ch.debug
1376   //printf("    QA27 TUBE from z = %1.2f to z= %1.2f\n",zd2,2*tubpar[2]+zd2);
1377
1378   zd2 += 2.*tubpar[2];
1379
1380   // transition x2zdc to recombination chamber : skewed cone
1381   conpar[0] = (10.-1.)/2.;
1382   conpar[1] = 5.4/2.;
1383   conpar[2] = 5.8/2.;
1384   conpar[3] = 6.3/2.;
1385   conpar[4] = 7.0/2.;
1386   gMC->Gsvolu("QA28", "CONE", idtmed[7], conpar, 5);
1387   gMC->Gspos("QA28", 1, "ZDCA", -7.9-0.175, 0., conpar[0]+0.5+zd2, irotpipe1, "ONLY");
1388   gMC->Gspos("QA28", 2, "ZDCA", 7.9+0.175, 0., conpar[0]+0.5+zd2, irotpipe2, "ONLY");
1389   //printf("    QA28 CONE from z = %1.2f to z= %1.2f\n",zd2,2*conpar[0]+0.2+zd2);
1390
1391   zd2 += 2.*conpar[0]+1.;
1392
1393   // 2 tubes (ID = 63 mm OD=70 mm)
1394   tubpar[0] = 6.3/2.;
1395   tubpar[1] = 7.0/2.;
1396   tubpar[2] = (342.5+498.3)/2.;
1397   gMC->Gsvolu("QA29", "TUBE", idtmed[7], tubpar, 3);
1398   gMC->Gspos("QA29", 1, "ZDCA", -16.5/2., 0., tubpar[2]+zd2, 0, "ONLY");
1399   gMC->Gspos("QA29", 2, "ZDCA",  16.5/2., 0., tubpar[2]+zd2, 0, "ONLY");
1400   //printf("    QA29 TUBE from z = %1.2f to z= %1.2f\n",zd2,2*tubpar[2]+zd2);
1401
1402   zd2 += 2.*tubpar[2];
1403
1404   // -- Luminometer (Cu box) in front of ZN - side A
1405   boxpar[0] = 8.0/2.;
1406   boxpar[1] = 8.0/2.;
1407   boxpar[2] = fLumiLength/2.;
1408   gMC->Gsvolu("QLUA", "BOX ", idtmed[9], boxpar, 3);
1409   gMC->Gspos("QLUA", 1, "ZDCA", 0., 0.,  fPosZNA[2]-66.-boxpar[2], 0, "ONLY");
1410   //printf("    QLUA LUMINOMETER from z = %1.2f to z= %1.2f\n\n",  fPosZNA[2]-66., fPosZNA[2]-66.-2*boxpar[2]);
1411
1412   //printf("    END OF BEAM PIPE VOLUME DEFINITION AT z = %f\n",zd2);
1413
1414
1415   // ----------------------------------------------------------------
1416   // --  MAGNET DEFINITION  -> LHC OPTICS 6.5
1417   // ----------------------------------------------------------------
1418   // ***************************************************************
1419   //            SIDE C - RB26  (dimuon side)
1420   // ***************************************************************
1421   // --  COMPENSATOR DIPOLE (MBXW)
1422   zCorrDip = 1972.5;
1423
1424   // --  GAP (VACUUM WITH MAGNETIC FIELD)
1425   tubpar[0] = 0.;
1426   tubpar[1] = 3.14;
1427   tubpar[2] = 153./2.;
1428   gMC->Gsvolu("MBXW", "TUBE", idtmed[11], tubpar, 3);
1429
1430   // --  YOKE
1431   tubpar[0] = 4.5;
1432   tubpar[1] = 55.;
1433   tubpar[2] = 153./2.;
1434   gMC->Gsvolu("YMBX", "TUBE", idtmed[7], tubpar, 3);
1435
1436   gMC->Gspos("MBXW", 1, "ZDCC", 0., 0., -tubpar[2]-zCorrDip, 0, "ONLY");
1437   gMC->Gspos("YMBX", 1, "ZDCC", 0., 0., -tubpar[2]-zCorrDip, 0, "ONLY");
1438
1439
1440   // -- INNER TRIPLET
1441   zInnTrip = 2296.5;
1442
1443   // -- DEFINE MQXL AND MQX QUADRUPOLE ELEMENT
1444   // --  MQXL
1445   // --  GAP (VACUUM WITH MAGNETIC FIELD)
1446   tubpar[0] = 0.;
1447   tubpar[1] = 3.14;
1448   tubpar[2] = 637./2.;
1449   gMC->Gsvolu("MQXL", "TUBE", idtmed[11], tubpar, 3);
1450
1451   // --  YOKE
1452   tubpar[0] = 3.5;
1453   tubpar[1] = 22.;
1454   tubpar[2] = 637./2.;
1455   gMC->Gsvolu("YMQL", "TUBE", idtmed[7], tubpar, 3);
1456
1457   gMC->Gspos("MQXL", 1, "ZDCC", 0., 0., -tubpar[2]-zInnTrip, 0, "ONLY");
1458   gMC->Gspos("YMQL", 1, "ZDCC", 0., 0., -tubpar[2]-zInnTrip, 0, "ONLY");
1459
1460   gMC->Gspos("MQXL", 2, "ZDCC", 0., 0., -tubpar[2]-zInnTrip-2400., 0, "ONLY");
1461   gMC->Gspos("YMQL", 2, "ZDCC", 0., 0., -tubpar[2]-zInnTrip-2400., 0, "ONLY");
1462
1463   // --  MQX
1464   // --  GAP (VACUUM WITH MAGNETIC FIELD)
1465   tubpar[0] = 0.;
1466   tubpar[1] = 3.14;
1467   tubpar[2] = 550./2.;
1468   gMC->Gsvolu("MQX ", "TUBE", idtmed[11], tubpar, 3);
1469
1470   // --  YOKE
1471   tubpar[0] = 3.5;
1472   tubpar[1] = 22.;
1473   tubpar[2] = 550./2.;
1474   gMC->Gsvolu("YMQ ", "TUBE", idtmed[7], tubpar, 3);
1475
1476   gMC->Gspos("MQX ", 1, "ZDCC", 0., 0., -tubpar[2]-zInnTrip-908.5,  0, "ONLY");
1477   gMC->Gspos("YMQ ", 1, "ZDCC", 0., 0., -tubpar[2]-zInnTrip-908.5,  0, "ONLY");
1478
1479   gMC->Gspos("MQX ", 2, "ZDCC", 0., 0., -tubpar[2]-zInnTrip-1558.5, 0, "ONLY");
1480   gMC->Gspos("YMQ ", 2, "ZDCC", 0., 0., -tubpar[2]-zInnTrip-1558.5, 0, "ONLY");
1481
1482   // -- SEPARATOR DIPOLE D1
1483   zD1 = 5838.3001;
1484
1485   // --  GAP (VACUUM WITH MAGNETIC FIELD)
1486   tubpar[0] = 0.;
1487   tubpar[1] = 3.46;
1488   tubpar[2] = 945./2.;
1489   gMC->Gsvolu("MD1 ", "TUBE", idtmed[11], tubpar, 3);
1490
1491   // --  Insert horizontal Cu plates inside D1
1492   // --   (to simulate the vacuum chamber)
1493   boxpar[0] = TMath::Sqrt(tubpar[1]*tubpar[1]-(2.98+0.2)*(2.98+0.2)) - 0.05;
1494   boxpar[1] = 0.2/2.;
1495   boxpar[2] = 945./2.;
1496   gMC->Gsvolu("MD1V", "BOX ", idtmed[6], boxpar, 3);
1497   gMC->Gspos("MD1V", 1, "MD1 ", 0., 2.98+boxpar[1], 0., 0, "ONLY");
1498   gMC->Gspos("MD1V", 2, "MD1 ", 0., -2.98-boxpar[1], 0., 0, "ONLY");
1499
1500   // --  YOKE
1501   tubpar[0] = 3.68;
1502   tubpar[1] = 110./2.;
1503   tubpar[2] = 945./2.;
1504   gMC->Gsvolu("YD1 ", "TUBE", idtmed[7], tubpar, 3);
1505
1506   gMC->Gspos("YD1 ", 1, "ZDCC", 0., 0., -tubpar[2]-zD1, 0, "ONLY");
1507   gMC->Gspos("MD1 ", 1, "ZDCC", 0., 0., -tubpar[2]-zD1, 0, "ONLY");
1508   // Ch debug
1509   //printf("    MD1 from z = %1.2f to z= %1.2f cm\n",-zD1, -zD1-2*tubpar[2]);
1510
1511   // -- DIPOLE D2
1512   zD2 = 12167.8;
1513   // --  GAP (VACUUM WITH MAGNETIC FIELD)
1514   tubpar[0] = 0.;
1515   tubpar[1] = 7.5/2.;
1516   tubpar[2] = 945./2.;
1517   gMC->Gsvolu("MD2 ", "TUBE", idtmed[11], tubpar, 3);
1518
1519   // --  YOKE
1520   tubpar[0] = 0.;
1521   tubpar[1] = 55.;
1522   tubpar[2] = 945./2.;
1523   gMC->Gsvolu("YD2 ", "TUBE", idtmed[7], tubpar, 3);
1524
1525   gMC->Gspos("YD2 ", 1, "ZDCC", 0., 0., -tubpar[2]-zD2, 0, "ONLY");
1526   // Ch debug
1527   //printf("    YD2 from z = %1.2f to z= %1.2f cm\n",-zD2, -zD2-2*tubpar[2]);
1528
1529   gMC->Gspos("MD2 ", 1, "YD2 ", -9.4, 0., 0., 0, "ONLY");
1530   gMC->Gspos("MD2 ", 2, "YD2 ",  9.4, 0., 0., 0, "ONLY");
1531
1532   // ***************************************************************
1533   //            SIDE A - RB24
1534   // ***************************************************************
1535
1536   // COMPENSATOR DIPOLE (MCBWA) (2nd compensator)
1537   // --  GAP (VACUUM WITH MAGNETIC FIELD)
1538   tubpar[0] = 0.;
1539   tubpar[1] = 3.;
1540   tubpar[2] = 153./2.;
1541   gMC->Gsvolu("MCBW", "TUBE", idtmed[11], tubpar, 3);
1542   gMC->Gspos("MCBW", 1, "ZDCA", 0., 0., tubpar[2]+zCorrDip, 0, "ONLY");
1543
1544    // --  YOKE
1545   tubpar[0] = 4.5;
1546   tubpar[1] = 55.;
1547   tubpar[2] = 153./2.;
1548   gMC->Gsvolu("YMCB", "TUBE", idtmed[7], tubpar, 3);
1549   gMC->Gspos("YMCB", 1, "ZDCA", 0., 0., tubpar[2]+zCorrDip, 0, "ONLY");
1550
1551    // -- INNER TRIPLET
1552   // -- DEFINE MQX1 AND MQX2 QUADRUPOLE ELEMENT
1553   // --  MQX1
1554   // --  GAP (VACUUM WITH MAGNETIC FIELD)
1555   tubpar[0] = 0.;
1556   tubpar[1] = 3.14;
1557   tubpar[2] = 637./2.;
1558   gMC->Gsvolu("MQX1", "TUBE", idtmed[11], tubpar, 3);
1559   gMC->Gsvolu("MQX4", "TUBE", idtmed[11], tubpar, 3);
1560
1561   // --  YOKE
1562   tubpar[0] = 3.5;
1563   tubpar[1] = 22.;
1564   tubpar[2] = 637./2.;
1565   gMC->Gsvolu("YMQ1", "TUBE", idtmed[7], tubpar, 3);
1566
1567   // -- Q1
1568   gMC->Gspos("MQX1", 1, "ZDCA", 0., 0., tubpar[2]+zInnTrip, 0, "ONLY");
1569   gMC->Gspos("YMQ1", 1, "ZDCA", 0., 0., tubpar[2]+zInnTrip, 0, "ONLY");
1570
1571    // -- BEAM SCREEN FOR Q1
1572    tubpar[0] = 4.78/2.;
1573    tubpar[1] = 5.18/2.;
1574    tubpar[2] = 637./2.;
1575    gMC->Gsvolu("QBS1", "TUBE", idtmed[6], tubpar, 3);
1576    gMC->Gspos("QBS1", 1, "MQX1", 0., 0., 0., 0, "ONLY");
1577    // INSERT VERTICAL PLATE INSIDE Q1
1578    boxpar[0] = 0.2/2.0;
1579    boxpar[1] = TMath::Sqrt(tubpar[0]*tubpar[0]-(1.9+0.2)*(1.9+0.2));
1580    boxpar[2] =637./2.;
1581    gMC->Gsvolu("QBS2", "BOX ", idtmed[6], boxpar, 3);
1582    gMC->Gspos("QBS2", 1, "MQX1", 1.9+boxpar[0], 0., 0., 0, "ONLY");
1583    gMC->Gspos("QBS2", 2, "MQX1", -1.9-boxpar[0], 0., 0., 0, "ONLY");
1584
1585    // -- Q3
1586    gMC->Gspos("MQX4", 1, "ZDCA", 0., 0., tubpar[2]+zInnTrip+2400., 0, "ONLY");
1587    gMC->Gspos("YMQ1", 2, "ZDCA", 0., 0., tubpar[2]+zInnTrip+2400., 0, "ONLY");
1588
1589    // -- BEAM SCREEN FOR Q3
1590    tubpar[0] = 5.79/2.;
1591    tubpar[1] = 6.14/2.;
1592    tubpar[2] = 637./2.;
1593    gMC->Gsvolu("QBS3", "TUBE", idtmed[6], tubpar, 3);
1594    gMC->Gspos("QBS3", 1, "MQX4", 0., 0., 0., 0, "ONLY");
1595    // INSERT VERTICAL PLATE INSIDE Q3
1596    boxpar[0] = 0.2/2.0;
1597    boxpar[1] = TMath::Sqrt(tubpar[0]*tubpar[0]-(2.405+0.2)*(2.405+0.2));
1598    boxpar[2] =637./2.;
1599    gMC->Gsvolu("QBS4", "BOX ", idtmed[6], boxpar, 3);
1600    gMC->Gspos("QBS4", 1, "MQX4", 2.405+boxpar[0], 0., 0., 0, "ONLY");
1601    gMC->Gspos("QBS4", 2, "MQX4", -2.405-boxpar[0], 0., 0., 0, "ONLY");
1602
1603
1604
1605   // --  MQX2
1606   // --  GAP (VACUUM WITH MAGNETIC FIELD)
1607   tubpar[0] = 0.;
1608   tubpar[1] = 3.14;
1609   tubpar[2] = 550./2.;
1610   gMC->Gsvolu("MQX2", "TUBE", idtmed[11], tubpar, 3);
1611   gMC->Gsvolu("MQX3", "TUBE", idtmed[11], tubpar, 3);
1612
1613   // --  YOKE
1614   tubpar[0] = 3.5;
1615   tubpar[1] = 22.;
1616   tubpar[2] = 550./2.;
1617   gMC->Gsvolu("YMQ2", "TUBE", idtmed[7], tubpar, 3);
1618
1619    // -- BEAM SCREEN FOR Q2
1620    tubpar[0] = 5.79/2.;
1621    tubpar[1] = 6.14/2.;
1622    tubpar[2] = 550./2.;
1623    gMC->Gsvolu("QBS5", "TUBE", idtmed[6], tubpar, 3);
1624    //    VERTICAL PLATE INSIDE Q2
1625    boxpar[0] = 0.2/2.0;
1626    boxpar[1] = TMath::Sqrt(tubpar[0]*tubpar[0]-(2.405+0.2)*(2.405+0.2));
1627    boxpar[2] =550./2.;
1628    gMC->Gsvolu("QBS6", "BOX ", idtmed[6], boxpar, 3);
1629
1630   // -- Q2A
1631   gMC->Gspos("MQX2", 1, "ZDCA", 0., 0., tubpar[2]+zInnTrip+908.5,  0, "ONLY");
1632   gMC->Gspos("QBS5", 1, "MQX2", 0., 0., 0., 0, "ONLY");
1633   gMC->Gspos("QBS6", 1, "MQX2", 2.405+boxpar[0], 0., 0., 0, "ONLY");
1634   gMC->Gspos("QBS6", 2, "MQX2", -2.405-boxpar[0], 0., 0., 0, "ONLY");
1635   gMC->Gspos("YMQ2", 1, "ZDCA", 0., 0., tubpar[2]+zInnTrip+908.5,  0, "ONLY");
1636
1637
1638   // -- Q2B
1639   gMC->Gspos("MQX3", 1, "ZDCA", 0., 0., tubpar[2]+zInnTrip+1558.5, 0, "ONLY");
1640   gMC->Gspos("QBS5", 2, "MQX3", 0., 0., 0., 0, "ONLY");
1641   gMC->Gspos("QBS6", 3, "MQX3", 2.405+boxpar[0], 0., 0., 0, "ONLY");
1642   gMC->Gspos("QBS6", 4, "MQX3", -2.405-boxpar[0], 0., 0., 0, "ONLY");
1643   gMC->Gspos("YMQ2", 2, "ZDCA", 0., 0., tubpar[2]+zInnTrip+1558.5, 0, "ONLY");
1644
1645   // -- SEPARATOR DIPOLE D1
1646   // --  GAP (VACUUM WITH MAGNETIC FIELD)
1647   tubpar[0] = 0.;
1648   tubpar[1] = 6.75/2.;//3.375
1649   tubpar[2] = 945./2.;
1650   gMC->Gsvolu("MD1L", "TUBE", idtmed[11], tubpar, 3);
1651
1652   // --  The beam screen tube is provided by the beam pipe in D1 (QA03 volume)
1653   // --  Insert the beam screen horizontal Cu plates inside D1
1654   // --   (to simulate the vacuum chamber)
1655   boxpar[0] = TMath::Sqrt(tubpar[1]*tubpar[1]-(2.885+0.2)*(2.885+0.2));
1656   boxpar[1] = 0.2/2.;
1657   boxpar[2] =945./2.;
1658   gMC->Gsvolu("QBS7", "BOX ", idtmed[6], boxpar, 3);
1659   gMC->Gspos("QBS7", 1, "MD1L", 0., 2.885+boxpar[1],0., 0, "ONLY");
1660   gMC->Gspos("QBS7", 2, "MD1L", 0., -2.885-boxpar[1],0., 0, "ONLY");
1661
1662   // --  YOKE
1663   tubpar[0] = 3.68;
1664   tubpar[1] = 110./2;
1665   tubpar[2] = 945./2.;
1666   gMC->Gsvolu("YD1L", "TUBE", idtmed[7], tubpar, 3);
1667
1668   gMC->Gspos("YD1L", 1, "ZDCA", 0., 0., tubpar[2]+zD1, 0, "ONLY");
1669   gMC->Gspos("MD1L", 1, "ZDCA", 0., 0., tubpar[2]+zD1, 0, "ONLY");
1670
1671   // -- DIPOLE D2
1672   // --  GAP (VACUUM WITH MAGNETIC FIELD)
1673   tubpar[0] = 0.;
1674   tubpar[1] = 7.5/2.; // this has to be checked
1675   tubpar[2] = 945./2.;
1676   gMC->Gsvolu("MD2L", "TUBE", idtmed[11], tubpar, 3);
1677
1678   // --  YOKE
1679   tubpar[0] = 0.;
1680   tubpar[1] = 55.;
1681   tubpar[2] = 945./2.;
1682   gMC->Gsvolu("YD2L", "TUBE", idtmed[7], tubpar, 3);
1683
1684   gMC->Gspos("YD2L", 1, "ZDCA", 0., 0., tubpar[2]+zD2, 0, "ONLY");
1685
1686   gMC->Gspos("MD2L", 1, "YD2L", -9.4, 0., 0., 0, "ONLY");
1687   gMC->Gspos("MD2L", 2, "YD2L",  9.4, 0., 0., 0, "ONLY");
1688
1689   // -- END OF MAGNET DEFINITION
1690 }
1691
1692 //_____________________________________________________________________________
1693 void AliZDCv3::CreateZDC()
1694 {
1695  //
1696  // Create the various ZDCs (ZN + ZP)
1697  //
1698
1699   Float_t dimPb[6], dimVoid[6];
1700
1701   Int_t *idtmed = fIdtmed->GetArray();
1702
1703   // Parameters for hadronic calorimeters geometry
1704   // NB -> parameters used ONLY in CreateZDC()
1705   Float_t fGrvZN[3] = {0.03, 0.03, 50.};  // Grooves for neutron detector
1706   Float_t fGrvZP[3] = {0.04, 0.04, 75.};  // Grooves for proton detector
1707   Int_t   fDivZN[3] = {11, 11, 0};        // Division for neutron detector
1708   Int_t   fDivZP[3] = {7, 15, 0};         // Division for proton detector
1709   Int_t   fTowZN[2] = {2, 2};             // Tower for neutron detector
1710   Int_t   fTowZP[2] = {4, 1};             // Tower for proton detector
1711
1712   // Parameters for EM calorimeter geometry
1713   // NB -> parameters used ONLY in CreateZDC()
1714   Float_t kDimZEMPb  = 0.15*(TMath::Sqrt(2.));  // z-dimension of the Pb slice
1715   Float_t kFibRadZEM = 0.0315;                  // External fiber radius (including cladding)
1716   Int_t   fDivZEM[3] = {92, 0, 20};             // Divisions for EM detector
1717   Float_t fDimZEM[6] = {fZEMLength, 3.5, 3.5, 45., 0., 0.}; // Dimensions of EM detector
1718   Float_t fFibZEM2 = fDimZEM[2]/TMath::Sin(fDimZEM[3]*kDegrad)-kFibRadZEM;
1719   Float_t fFibZEM[3] = {0., 0.0275, fFibZEM2};  // Fibers for EM calorimeter
1720
1721
1722   //-- Create calorimeters geometry
1723
1724   // -------------------------------------------------------------------------------
1725   //--> Neutron calorimeter (ZN)
1726
1727   gMC->Gsvolu("ZNEU", "BOX ", idtmed[1], fDimZN, 3); // Passive material
1728   gMC->Gsvolu("ZNF1", "TUBE", idtmed[3], fFibZN, 3); // Active material
1729   gMC->Gsvolu("ZNF2", "TUBE", idtmed[4], fFibZN, 3);
1730   gMC->Gsvolu("ZNF3", "TUBE", idtmed[4], fFibZN, 3);
1731   gMC->Gsvolu("ZNF4", "TUBE", idtmed[3], fFibZN, 3);
1732   gMC->Gsvolu("ZNG1", "BOX ", idtmed[12], fGrvZN, 3); // Empty grooves
1733   gMC->Gsvolu("ZNG2", "BOX ", idtmed[12], fGrvZN, 3);
1734   gMC->Gsvolu("ZNG3", "BOX ", idtmed[12], fGrvZN, 3);
1735   gMC->Gsvolu("ZNG4", "BOX ", idtmed[12], fGrvZN, 3);
1736
1737   // Divide ZNEU in towers (for hits purposes)
1738
1739   gMC->Gsdvn("ZNTX", "ZNEU", fTowZN[0], 1); // x-tower
1740   gMC->Gsdvn("ZN1 ", "ZNTX", fTowZN[1], 2); // y-tower
1741
1742   //-- Divide ZN1 in minitowers
1743   //  fDivZN[0]= NUMBER OF FIBERS PER TOWER ALONG X-AXIS,
1744   //  fDivZN[1]= NUMBER OF FIBERS PER TOWER ALONG Y-AXIS
1745   //  (4 fibres per minitower)
1746
1747   gMC->Gsdvn("ZNSL", "ZN1 ", fDivZN[1], 2); // Slices
1748   gMC->Gsdvn("ZNST", "ZNSL", fDivZN[0], 1); // Sticks
1749
1750   // --- Position the empty grooves in the sticks (4 grooves per stick)
1751   Float_t dx = fDimZN[0] / fDivZN[0] / 4.;
1752   Float_t dy = fDimZN[1] / fDivZN[1] / 4.;
1753
1754   gMC->Gspos("ZNG1", 1, "ZNST", 0.-dx, 0.+dy, 0., 0, "ONLY");
1755   gMC->Gspos("ZNG2", 1, "ZNST", 0.+dx, 0.+dy, 0., 0, "ONLY");
1756   gMC->Gspos("ZNG3", 1, "ZNST", 0.-dx, 0.-dy, 0., 0, "ONLY");
1757   gMC->Gspos("ZNG4", 1, "ZNST", 0.+dx, 0.-dy, 0., 0, "ONLY");
1758
1759   // --- Position the fibers in the grooves
1760   gMC->Gspos("ZNF1", 1, "ZNG1", 0., 0., 0., 0, "ONLY");
1761   gMC->Gspos("ZNF2", 1, "ZNG2", 0., 0., 0., 0, "ONLY");
1762   gMC->Gspos("ZNF3", 1, "ZNG3", 0., 0., 0., 0, "ONLY");
1763   gMC->Gspos("ZNF4", 1, "ZNG4", 0., 0., 0., 0, "ONLY");
1764
1765   // --- Position the neutron calorimeter in ZDC
1766   // -- Rotation of ZDCs
1767   Int_t irotzdc;
1768   gMC->Matrix(irotzdc, 90., 180., 90., 90., 180., 0.);
1769   //
1770   gMC->Gspos("ZNEU", 1, "ZDCC", fPosZNC[0], fPosZNC[1], fPosZNC[2]-fDimZN[2], irotzdc, "ONLY");
1771   //Ch debug
1772   //printf("\n ZN -> %f < z < %f cm\n",fPosZN[2],fPosZN[2]-2*fDimZN[2]);
1773
1774   // --- Position the neutron calorimeter in ZDC2 (left line)
1775   // -- No Rotation of ZDCs
1776   gMC->Gspos("ZNEU", 2, "ZDCA", fPosZNA[0], fPosZNA[1], fPosZNA[2]+fDimZN[2], 0, "ONLY");
1777   //Ch debug
1778   //printf("\n ZN left -> %f < z < %f cm\n",fPosZNl[2],fPosZNl[2]+2*fDimZN[2]);
1779
1780
1781   // -------------------------------------------------------------------------------
1782   //--> Proton calorimeter (ZP)
1783
1784   gMC->Gsvolu("ZPRO", "BOX ", idtmed[2], fDimZP, 3); // Passive material
1785   gMC->Gsvolu("ZPF1", "TUBE", idtmed[3], fFibZP, 3); // Active material
1786   gMC->Gsvolu("ZPF2", "TUBE", idtmed[4], fFibZP, 3);
1787   gMC->Gsvolu("ZPF3", "TUBE", idtmed[4], fFibZP, 3);
1788   gMC->Gsvolu("ZPF4", "TUBE", idtmed[3], fFibZP, 3);
1789   gMC->Gsvolu("ZPG1", "BOX ", idtmed[12], fGrvZP, 3); // Empty grooves
1790   gMC->Gsvolu("ZPG2", "BOX ", idtmed[12], fGrvZP, 3);
1791   gMC->Gsvolu("ZPG3", "BOX ", idtmed[12], fGrvZP, 3);
1792   gMC->Gsvolu("ZPG4", "BOX ", idtmed[12], fGrvZP, 3);
1793
1794   //-- Divide ZPRO in towers(for hits purposes)
1795
1796   gMC->Gsdvn("ZPTX", "ZPRO", fTowZP[0], 1); // x-tower
1797   gMC->Gsdvn("ZP1 ", "ZPTX", fTowZP[1], 2); // y-tower
1798
1799
1800   //-- Divide ZP1 in minitowers
1801   //  fDivZP[0]= NUMBER OF FIBERS ALONG X-AXIS PER MINITOWER,
1802   //  fDivZP[1]= NUMBER OF FIBERS ALONG Y-AXIS PER MINITOWER
1803   //  (4 fiber per minitower)
1804
1805   gMC->Gsdvn("ZPSL", "ZP1 ", fDivZP[1], 2); // Slices
1806   gMC->Gsdvn("ZPST", "ZPSL", fDivZP[0], 1); // Sticks
1807
1808   // --- Position the empty grooves in the sticks (4 grooves per stick)
1809   dx = fDimZP[0] / fTowZP[0] / fDivZP[0] / 2.;
1810   dy = fDimZP[1] / fTowZP[1] / fDivZP[1] / 2.;
1811
1812   gMC->Gspos("ZPG1", 1, "ZPST", 0.-dx, 0.+dy, 0., 0, "ONLY");
1813   gMC->Gspos("ZPG2", 1, "ZPST", 0.+dx, 0.+dy, 0., 0, "ONLY");
1814   gMC->Gspos("ZPG3", 1, "ZPST", 0.-dx, 0.-dy, 0., 0, "ONLY");
1815   gMC->Gspos("ZPG4", 1, "ZPST", 0.+dx, 0.-dy, 0., 0, "ONLY");
1816
1817   // --- Position the fibers in the grooves
1818   gMC->Gspos("ZPF1", 1, "ZPG1", 0., 0., 0., 0, "ONLY");
1819   gMC->Gspos("ZPF2", 1, "ZPG2", 0., 0., 0., 0, "ONLY");
1820   gMC->Gspos("ZPF3", 1, "ZPG3", 0., 0., 0., 0, "ONLY");
1821   gMC->Gspos("ZPF4", 1, "ZPG4", 0., 0., 0., 0, "ONLY");
1822
1823
1824   // --- Position the proton calorimeter in ZDCC
1825   gMC->Gspos("ZPRO", 1, "ZDCC", fPosZPC[0], fPosZPC[1], fPosZPC[2]-fDimZP[2], irotzdc, "ONLY");
1826   //Ch debug
1827   //printf("\n ZP -> %f < z < %f cm\n",fPosZP[2],fPosZP[2]-2*fDimZP[2]);
1828
1829   // --- Position the proton calorimeter in ZDCA
1830   // --- No rotation
1831   gMC->Gspos("ZPRO", 2, "ZDCA", fPosZPA[0], fPosZPA[1], fPosZPA[2]+fDimZP[2], 0, "ONLY");
1832   //Ch debug
1833   //printf("\n ZP left -> %f < z < %f cm\n",fPosZPl[2],fPosZPl[2]+2*fDimZP[2]);
1834
1835
1836   // -------------------------------------------------------------------------------
1837   // -> EM calorimeter (ZEM)
1838
1839   gMC->Gsvolu("ZEM ", "PARA", idtmed[10], fDimZEM, 6);
1840
1841   Int_t irot1, irot2;
1842   gMC->Matrix(irot1,0.,0.,90.,90.,-90.,0.);                    // Rotation matrix 1
1843   gMC->Matrix(irot2,180.,0.,90.,fDimZEM[3]+90.,90.,fDimZEM[3]);// Rotation matrix 2
1844   //printf("irot1 = %d, irot2 = %d \n", irot1, irot2);
1845
1846   gMC->Gsvolu("ZEMF", "TUBE", idtmed[3], fFibZEM, 3);   // Active material
1847
1848   gMC->Gsdvn("ZETR", "ZEM ", fDivZEM[2], 1);            // Tranches
1849
1850   dimPb[0] = kDimZEMPb;                                 // Lead slices
1851   dimPb[1] = fDimZEM[2];
1852   dimPb[2] = fDimZEM[1];
1853   //dimPb[3] = fDimZEM[3]; //controllare
1854   dimPb[3] = 90.-fDimZEM[3]; //originale
1855   dimPb[4] = 0.;
1856   dimPb[5] = 0.;
1857   gMC->Gsvolu("ZEL0", "PARA", idtmed[5], dimPb, 6);
1858   gMC->Gsvolu("ZEL1", "PARA", idtmed[5], dimPb, 6);
1859   gMC->Gsvolu("ZEL2", "PARA", idtmed[5], dimPb, 6);
1860
1861   // --- Position the lead slices in the tranche
1862   Float_t zTran = fDimZEM[0]/fDivZEM[2];
1863   Float_t zTrPb = -zTran+kDimZEMPb;
1864   gMC->Gspos("ZEL0", 1, "ZETR", zTrPb, 0., 0., 0, "ONLY");
1865   gMC->Gspos("ZEL1", 1, "ZETR", kDimZEMPb, 0., 0., 0, "ONLY");
1866
1867   // --- Vacuum zone (to be filled with fibres)
1868   dimVoid[0] = (zTran-2*kDimZEMPb)/2.;
1869   dimVoid[1] = fDimZEM[2];
1870   dimVoid[2] = fDimZEM[1];
1871   dimVoid[3] = 90.-fDimZEM[3];
1872   dimVoid[4] = 0.;
1873   dimVoid[5] = 0.;
1874   gMC->Gsvolu("ZEV0", "PARA", idtmed[10], dimVoid,6);
1875   gMC->Gsvolu("ZEV1", "PARA", idtmed[10], dimVoid,6);
1876
1877   // --- Divide the vacuum slice into sticks along x axis
1878   gMC->Gsdvn("ZES0", "ZEV0", fDivZEM[0], 3);
1879   gMC->Gsdvn("ZES1", "ZEV1", fDivZEM[0], 3);
1880
1881   // --- Positioning the fibers into the sticks
1882   gMC->Gspos("ZEMF", 1,"ZES0", 0., 0., 0., irot2, "ONLY");
1883   gMC->Gspos("ZEMF", 1,"ZES1", 0., 0., 0., irot2, "ONLY");
1884
1885   // --- Positioning the vacuum slice into the tranche
1886   //Float_t displFib = fDimZEM[1]/fDivZEM[0];
1887   gMC->Gspos("ZEV0", 1,"ZETR", -dimVoid[0], 0., 0., 0, "ONLY");
1888   gMC->Gspos("ZEV1", 1,"ZETR", -dimVoid[0]+zTran, 0., 0., 0, "ONLY");
1889
1890   // --- Positioning the ZEM into the ZDC - rotation for 90 degrees
1891   // NB -> ZEM is positioned in ALIC (instead of in ZDC) volume
1892   gMC->Gspos("ZEM ", 1,"ALIC", -fPosZEM[0], fPosZEM[1], fPosZEM[2]+fDimZEM[0], irot1, "ONLY");
1893
1894   // Second EM ZDC (same side w.r.t. IP, just on the other side w.r.t. beam pipe)
1895   gMC->Gspos("ZEM ", 2,"ALIC", fPosZEM[0], fPosZEM[1], fPosZEM[2]+fDimZEM[0], irot1, "ONLY");
1896
1897   // --- Adding last slice at the end of the EM calorimeter
1898   Float_t zLastSlice = fPosZEM[2]+kDimZEMPb+2*fDimZEM[0];
1899   gMC->Gspos("ZEL2", 1,"ALIC", fPosZEM[0], fPosZEM[1], zLastSlice, irot1, "ONLY");
1900   //Ch debug
1901   //printf("\n ZEM lenght = %f cm\n",2*fZEMLength);
1902   //printf("\n ZEM -> %f < z < %f cm\n",fPosZEM[2],fPosZEM[2]+2*fZEMLength+zLastSlice+kDimZEMPb);
1903
1904 }
1905
1906 //_____________________________________________________________________________
1907 void AliZDCv3::DrawModule() const
1908 {
1909   //
1910   // Draw a shaded view of the Zero Degree Calorimeter version 1
1911   //
1912
1913   // Set everything unseen
1914   gMC->Gsatt("*", "seen", -1);
1915   //
1916   // Set ALIC mother transparent
1917   gMC->Gsatt("ALIC","SEEN",0);
1918   //
1919   // Set the volumes visible
1920   gMC->Gsatt("ZDCC","SEEN",0);
1921   gMC->Gsatt("QT01","SEEN",1);
1922   gMC->Gsatt("QT02","SEEN",1);
1923   gMC->Gsatt("QT03","SEEN",1);
1924   gMC->Gsatt("QT04","SEEN",1);
1925   gMC->Gsatt("QT05","SEEN",1);
1926   gMC->Gsatt("QT06","SEEN",1);
1927   gMC->Gsatt("QT07","SEEN",1);
1928   gMC->Gsatt("QT08","SEEN",1);
1929   gMC->Gsatt("QT09","SEEN",1);
1930   gMC->Gsatt("QT10","SEEN",1);
1931   gMC->Gsatt("QT11","SEEN",1);
1932   gMC->Gsatt("QT12","SEEN",1);
1933   gMC->Gsatt("QT13","SEEN",1);
1934   gMC->Gsatt("QC01","SEEN",1);
1935   gMC->Gsatt("QC02","SEEN",1);
1936   gMC->Gsatt("QC03","SEEN",1);
1937   gMC->Gsatt("QC04","SEEN",1);
1938   gMC->Gsatt("QC05","SEEN",1);
1939   gMC->Gsatt("QC06","SEEN",1);
1940   gMC->Gsatt("QC07","SEEN",1);
1941   gMC->Gsatt("QC08","SEEN",1);
1942   gMC->Gsatt("QC09","SEEN",1);
1943   gMC->Gsatt("QC10","SEEN",1);
1944   gMC->Gsatt("MQXL","SEEN",1);
1945   gMC->Gsatt("YMQL","SEEN",1);
1946   gMC->Gsatt("MQX ","SEEN",1);
1947   gMC->Gsatt("YMQ ","SEEN",1);
1948   gMC->Gsatt("ZQYX","SEEN",1);
1949   gMC->Gsatt("MD1 ","SEEN",1);
1950   gMC->Gsatt("MD1V","SEEN",1);
1951   gMC->Gsatt("YD1 ","SEEN",1);
1952   gMC->Gsatt("MD2 ","SEEN",1);
1953   gMC->Gsatt("YD2 ","SEEN",1);
1954   gMC->Gsatt("ZNEU","SEEN",0);
1955   gMC->Gsatt("ZNF1","SEEN",0);
1956   gMC->Gsatt("ZNF2","SEEN",0);
1957   gMC->Gsatt("ZNF3","SEEN",0);
1958   gMC->Gsatt("ZNF4","SEEN",0);
1959   gMC->Gsatt("ZNG1","SEEN",0);
1960   gMC->Gsatt("ZNG2","SEEN",0);
1961   gMC->Gsatt("ZNG3","SEEN",0);
1962   gMC->Gsatt("ZNG4","SEEN",0);
1963   gMC->Gsatt("ZNTX","SEEN",0);
1964   gMC->Gsatt("ZN1 ","COLO",4);
1965   gMC->Gsatt("ZN1 ","SEEN",1);
1966   gMC->Gsatt("ZNSL","SEEN",0);
1967   gMC->Gsatt("ZNST","SEEN",0);
1968   gMC->Gsatt("ZPRO","SEEN",0);
1969   gMC->Gsatt("ZPF1","SEEN",0);
1970   gMC->Gsatt("ZPF2","SEEN",0);
1971   gMC->Gsatt("ZPF3","SEEN",0);
1972   gMC->Gsatt("ZPF4","SEEN",0);
1973   gMC->Gsatt("ZPG1","SEEN",0);
1974   gMC->Gsatt("ZPG2","SEEN",0);
1975   gMC->Gsatt("ZPG3","SEEN",0);
1976   gMC->Gsatt("ZPG4","SEEN",0);
1977   gMC->Gsatt("ZPTX","SEEN",0);
1978   gMC->Gsatt("ZP1 ","COLO",6);
1979   gMC->Gsatt("ZP1 ","SEEN",1);
1980   gMC->Gsatt("ZPSL","SEEN",0);
1981   gMC->Gsatt("ZPST","SEEN",0);
1982   gMC->Gsatt("ZEM ","COLO",7);
1983   gMC->Gsatt("ZEM ","SEEN",1);
1984   gMC->Gsatt("ZEMF","SEEN",0);
1985   gMC->Gsatt("ZETR","SEEN",0);
1986   gMC->Gsatt("ZEL0","SEEN",0);
1987   gMC->Gsatt("ZEL1","SEEN",0);
1988   gMC->Gsatt("ZEL2","SEEN",0);
1989   gMC->Gsatt("ZEV0","SEEN",0);
1990   gMC->Gsatt("ZEV1","SEEN",0);
1991   gMC->Gsatt("ZES0","SEEN",0);
1992   gMC->Gsatt("ZES1","SEEN",0);
1993   //
1994   gMC->Gdopt("hide", "on");
1995   gMC->Gdopt("shad", "on");
1996   gMC->Gsatt("*", "fill", 7);
1997   gMC->SetClipBox(".");
1998   gMC->SetClipBox("*", 0, 100, -100, 100, 12000, 16000);
1999   gMC->DefaultRange();
2000   gMC->Gdraw("alic", 40, 30, 0, 488, 220, .07, .07);
2001   gMC->Gdhead(1111, "Zero Degree Calorimeter Version 3");
2002   gMC->Gdman(18, 4, "MAN");
2003 }
2004
2005 //_____________________________________________________________________________
2006 void AliZDCv3::CreateMaterials()
2007 {
2008   //
2009   // Create Materials for the Zero Degree Calorimeter
2010   //
2011   Float_t dens, ubuf[1], wmat[3], a[3], z[3];
2012
2013   // --- W alloy -> ZN passive material
2014   dens = 17.6;
2015   a[0] = 183.85;
2016   a[1] = 55.85;
2017   a[2] = 58.71;
2018   z[0] = 74.;
2019   z[1] = 26.;
2020   z[2] = 28.;
2021   wmat[0] = .93;
2022   wmat[1] = .03;
2023   wmat[2] = .04;
2024   AliMixture(1, "WALL", a, z, dens, 3, wmat);
2025
2026   // --- Brass (CuZn)  -> ZP passive material
2027   dens = 8.48;
2028   a[0] = 63.546;
2029   a[1] = 65.39;
2030   z[0] = 29.;
2031   z[1] = 30.;
2032   wmat[0] = .63;
2033   wmat[1] = .37;
2034   AliMixture(2, "BRASS", a, z, dens, 2, wmat);
2035
2036   // --- SiO2
2037   dens = 2.64;
2038   a[0] = 28.086;
2039   a[1] = 15.9994;
2040   z[0] = 14.;
2041   z[1] = 8.;
2042   wmat[0] = 1.;
2043   wmat[1] = 2.;
2044   AliMixture(3, "SIO2", a, z, dens, -2, wmat);
2045
2046   // --- Lead
2047   ubuf[0] = 1.12;
2048   AliMaterial(5, "LEAD", 207.19, 82., 11.35, .56, 0., ubuf, 1);
2049
2050   // --- Copper (energy loss taken into account)
2051   ubuf[0] = 1.10;
2052   AliMaterial(6, "COPP0", 63.54, 29., 8.96, 1.4, 0., ubuf, 1);
2053
2054   // --- Copper
2055   ubuf[0] = 1.10;
2056   AliMaterial(9, "COPP1", 63.54, 29., 8.96, 1.4, 0., ubuf, 1);
2057
2058   // --- Iron (energy loss taken into account)
2059   ubuf[0] = 1.1;
2060   AliMaterial(7, "IRON0", 55.85, 26., 7.87, 1.76, 0., ubuf, 1);
2061
2062   // --- Iron (no energy loss)
2063   ubuf[0] = 1.1;
2064   AliMaterial(8, "IRON1", 55.85, 26., 7.87, 1.76, 0., ubuf, 1);
2065
2066   // --- Tatalum
2067   ubuf[0] = 1.1;
2068   AliMaterial(13, "TANT", 183.84, 74., 19.3, 0.35, 0., ubuf, 1);
2069
2070   // ---------------------------------------------------------
2071   Float_t aResGas[3]={1.008,12.0107,15.9994};
2072   Float_t zResGas[3]={1.,6.,8.};
2073   Float_t wResGas[3]={0.28,0.28,0.44};
2074   Float_t dResGas = 3.2E-14;
2075
2076   // --- Vacuum (no magnetic field)
2077   AliMixture(10, "VOID", aResGas, zResGas, dResGas, 3, wResGas);
2078
2079   // --- Vacuum (with magnetic field)
2080   AliMixture(11, "VOIM", aResGas, zResGas, dResGas, 3, wResGas);
2081
2082   // --- Air (no magnetic field)
2083   Float_t aAir[4]={12.0107,14.0067,15.9994,39.948};
2084   Float_t zAir[4]={6.,7.,8.,18.};
2085   Float_t wAir[4]={0.000124,0.755267,0.231781,0.012827};
2086   Float_t dAir = 1.20479E-3;
2087   //
2088   AliMixture(12, "Air    $", aAir, zAir, dAir, 4, wAir);
2089
2090   // ---  Definition of tracking media:
2091
2092   // --- Tantalum = 1 ;
2093   // --- Brass = 2 ;
2094   // --- Fibers (SiO2) = 3 ;
2095   // --- Fibers (SiO2) = 4 ;
2096   // --- Lead = 5 ;
2097   // --- Copper (with high thr.)= 6 ;
2098   // --- Copper (with low thr.)=  9;
2099   // --- Iron (with energy loss) = 7 ;
2100   // --- Iron (without energy loss) = 8 ;
2101   // --- Vacuum (no field) = 10
2102   // --- Vacuum (with field) = 11
2103   // --- Air (no field) = 12
2104
2105   // ****************************************************
2106   //     Tracking media parameters
2107   //
2108   Float_t epsil  = 0.01;   // Tracking precision,
2109   Float_t stmin  = 0.01;   // Min. value 4 max. step (cm)
2110   Float_t stemax = 1.;     // Max. step permitted (cm)
2111   Float_t tmaxfd = 0.;     // Maximum angle due to field (degrees)
2112   Float_t tmaxfdv = 0.1;   // Maximum angle due to field (degrees)
2113   Float_t deemax = -1.;    // Maximum fractional energy loss
2114   Float_t nofieldm = 0.;   // Max. field value (no field)
2115   Float_t fieldm = 45.;    // Max. field value (with field)
2116   Int_t isvol = 0;         // ISVOL =0 -> not sensitive volume
2117   Int_t isvolActive = 1;   // ISVOL =1 -> sensitive volume
2118   Int_t inofld = 0;        // IFIELD=0 -> no magnetic field
2119   Int_t ifield =2;         // IFIELD=2 -> magnetic field defined in AliMagFC.h
2120   // *****************************************************
2121
2122   AliMedium(1, "ZWALL", 1, isvolActive, inofld, nofieldm, tmaxfd, stemax, deemax, epsil, stmin);
2123   AliMedium(2, "ZBRASS",2, isvolActive, inofld, nofieldm, tmaxfd, stemax, deemax, epsil, stmin);
2124   AliMedium(3, "ZSIO2", 3, isvolActive, inofld, nofieldm, tmaxfd, stemax, deemax, epsil, stmin);
2125   AliMedium(4, "ZQUAR", 3, isvolActive, inofld, nofieldm, tmaxfd, stemax, deemax, epsil, stmin);
2126   AliMedium(5, "ZLEAD", 5, isvolActive, inofld, nofieldm, tmaxfd, stemax, deemax, epsil, stmin);
2127   AliMedium(6, "ZCOPP", 6, isvol, inofld, nofieldm, tmaxfd, stemax, deemax, epsil, stmin);
2128   AliMedium(7, "ZIRON", 7, isvol, inofld, nofieldm, tmaxfd, stemax, deemax, epsil, stmin);
2129   AliMedium(8, "ZIRONN",8, isvol, inofld, nofieldm, tmaxfd, stemax, deemax, epsil, stmin);
2130   AliMedium(9, "ZCOPL", 6, isvol, inofld, nofieldm, tmaxfd, stemax, deemax, epsil, stmin);
2131   AliMedium(10,"ZVOID",10, isvol, inofld, nofieldm, tmaxfd, stemax, deemax, epsil, stmin);
2132   AliMedium(11,"ZVOIM",11, isvol, ifield, fieldm, tmaxfdv, stemax, deemax, epsil, stmin);
2133   AliMedium(12,"ZAIR", 12, isvolActive, inofld, nofieldm, tmaxfd, stemax, deemax, epsil, stmin);
2134   AliMedium(13,"ZTANT",13, isvolActive, inofld, nofieldm, tmaxfd, stemax, deemax, epsil, stmin);
2135   AliMedium(14, "ZIRONT", 7, isvol, inofld, nofieldm, tmaxfd, stemax, deemax, epsil, stmin);
2136
2137 }
2138
2139 //_____________________________________________________________________________
2140 void AliZDCv3::AddAlignableVolumes() const
2141 {
2142  //
2143  // Create entries for alignable volumes associating the symbolic volume
2144  // name with the corresponding volume path. Needs to be syncronized with
2145  // eventual changes in the geometry.
2146  //
2147  TString volpath1 = "ALIC_1/ZDCC_1/ZNEU_1";
2148  TString volpath2 = "ALIC_1/ZDCC_1/ZPRO_1";
2149  TString volpath3 = "ALIC_1/ZDCA_1/ZNEU_2";
2150  TString volpath4 = "ALIC_1/ZDCA_1/ZPRO_2";
2151
2152  TString symname1="ZDC/NeutronZDC_C";
2153  TString symname2="ZDC/ProtonZDC_C";
2154  TString symname3="ZDC/NeutronZDC_A";
2155  TString symname4="ZDC/ProtonZDC_A";
2156
2157  if(!gGeoManager->SetAlignableEntry(symname1.Data(),volpath1.Data()))
2158      AliFatal(Form("Alignable entry %s not created. Volume path %s not valid",   symname1.Data(),volpath1.Data()));
2159
2160  if(!gGeoManager->SetAlignableEntry(symname2.Data(),volpath2.Data()))
2161      AliFatal(Form("Alignable entry %s not created. Volume path %s not valid",   symname2.Data(),volpath2.Data()));
2162
2163  if(!gGeoManager->SetAlignableEntry(symname3.Data(),volpath3.Data()))
2164      AliFatal(Form("Alignable entry %s not created. Volume path %s not valid",   symname1.Data(),volpath1.Data()));
2165
2166  if(!gGeoManager->SetAlignableEntry(symname4.Data(),volpath4.Data()))
2167      AliFatal(Form("Alignable entry %s not created. Volume path %s not valid",   symname2.Data(),volpath2.Data()));
2168
2169 }
2170
2171
2172 //_____________________________________________________________________________
2173 void AliZDCv3::Init()
2174 {
2175  InitTables();
2176   Int_t *idtmed = fIdtmed->GetArray();
2177   //
2178   fMedSensZN     = idtmed[1];  // Sensitive volume: ZN passive material
2179   fMedSensZP     = idtmed[2];  // Sensitive volume: ZP passive material
2180   fMedSensF1     = idtmed[3];  // Sensitive volume: fibres type 1
2181   fMedSensF2     = idtmed[4];  // Sensitive volume: fibres type 2
2182   fMedSensZEM    = idtmed[5];  // Sensitive volume: ZEM passive material
2183   fMedSensTDI    = idtmed[6];  // Sensitive volume: TDI Cu shield
2184   fMedSensPI     = idtmed[7];  // Sensitive volume: beam pipes
2185   fMedSensLumi   = idtmed[9];  // Sensitive volume: luminometer
2186   fMedSensGR     = idtmed[12]; // Sensitive volume: air into the grooves
2187   fMedSensVColl  = idtmed[13]; // Sensitive volume: collimator jaws
2188 }
2189
2190 //_____________________________________________________________________________
2191 void AliZDCv3::InitTables()
2192 {
2193  //
2194  // Read light tables for Cerenkov light production parameterization
2195  //
2196
2197   Int_t k, j;
2198
2199   char *lightfName1,*lightfName2,*lightfName3,*lightfName4,
2200        *lightfName5,*lightfName6,*lightfName7,*lightfName8;
2201   FILE *fp1, *fp2, *fp3, *fp4, *fp5, *fp6, *fp7, *fp8;
2202
2203   //  --- Reading light tables for ZN
2204   lightfName1 = gSystem->ExpandPathName("$ALICE_ROOT/ZDC/light22620362207s");
2205   if((fp1 = fopen(lightfName1,"r")) == NULL){
2206      printf("Cannot open file fp1 \n");
2207      return;
2208   }
2209   lightfName2 = gSystem->ExpandPathName("$ALICE_ROOT/ZDC/light22620362208s");
2210   if((fp2 = fopen(lightfName2,"r")) == NULL){
2211      printf("Cannot open file fp2 \n");
2212      return;
2213   }
2214   lightfName3 = gSystem->ExpandPathName("$ALICE_ROOT/ZDC/light22620362209s");
2215   if((fp3 = fopen(lightfName3,"r")) == NULL){
2216      printf("Cannot open file fp3 \n");
2217      return;
2218   }
2219   lightfName4 = gSystem->ExpandPathName("$ALICE_ROOT/ZDC/light22620362210s");
2220   if((fp4 = fopen(lightfName4,"r")) == NULL){
2221      printf("Cannot open file fp4 \n");
2222      return;
2223   }
2224
2225   int read=1;
2226   for(k=0; k<fNalfan; k++){
2227      for(j=0; j<fNben; j++){
2228        read = fscanf(fp1,"%f",&fTablen[0][k][j]);
2229        if(read==0) AliDebug(3, " Error in reading light table 1");
2230        read = fscanf(fp2,"%f",&fTablen[1][k][j]);
2231        if(read==0) AliDebug(3, " Error in reading light table 2");
2232        read = fscanf(fp3,"%f",&fTablen[2][k][j]);
2233        if(read==0) AliDebug(3, " Error in reading light table 3");
2234        read = fscanf(fp4,"%f",&fTablen[3][k][j]);
2235        if(read==0) AliDebug(3, " Error in reading light table 4");
2236      }
2237   }
2238   fclose(fp1);
2239   fclose(fp2);
2240   fclose(fp3);
2241   fclose(fp4);
2242
2243   //  --- Reading light tables for ZP and ZEM
2244   lightfName5 = gSystem->ExpandPathName("$ALICE_ROOT/ZDC/light22620552207s");
2245   if((fp5 = fopen(lightfName5,"r")) == NULL){
2246      printf("Cannot open file fp5 \n");
2247      return;
2248   }
2249   lightfName6 = gSystem->ExpandPathName("$ALICE_ROOT/ZDC/light22620552208s");
2250   if((fp6 = fopen(lightfName6,"r")) == NULL){
2251      printf("Cannot open file fp6 \n");
2252      return;
2253   }
2254   lightfName7 = gSystem->ExpandPathName("$ALICE_ROOT/ZDC/light22620552209s");
2255   if((fp7 = fopen(lightfName7,"r")) == NULL){
2256      printf("Cannot open file fp7 \n");
2257      return;
2258   }
2259   lightfName8 = gSystem->ExpandPathName("$ALICE_ROOT/ZDC/light22620552210s");
2260   if((fp8 = fopen(lightfName8,"r")) == NULL){
2261      printf("Cannot open file fp8 \n");
2262      return;
2263   }
2264
2265   for(k=0; k<fNalfap; k++){
2266      for(j=0; j<fNbep; j++){
2267        read = fscanf(fp5,"%f",&fTablep[0][k][j]);
2268        if(read==0) AliDebug(3, " Error in reading light table 5");
2269        read = fscanf(fp6,"%f",&fTablep[1][k][j]);
2270        if(read==0) AliDebug(3, " Error in reading light table 6");
2271        read = fscanf(fp7,"%f",&fTablep[2][k][j]);
2272        if(read==0) AliDebug(3, " Error in reading light table 7");
2273        read = fscanf(fp8,"%f",&fTablep[3][k][j]);
2274        if(read==0) AliDebug(3, " Error in reading light table 8");
2275      }
2276   }
2277   fclose(fp5);
2278   fclose(fp6);
2279   fclose(fp7);
2280   fclose(fp8);
2281 }
2282 //_____________________________________________________________________________
2283 void AliZDCv3::StepManager()
2284 {
2285   //
2286   // Routine called at every step in the Zero Degree Calorimeters
2287   //
2288   Int_t   j, vol[2]={0,0}, ibeta=0, ialfa=0, ibe=0, nphe=0;
2289   Float_t hits[13], x[3], xdet[3], um[3], ud[3];
2290   Float_t destep=0., be=0., out=0.;
2291   Double_t s[3], p[4];
2292   const char *knamed;
2293   //
2294   for(j=0;j<13;j++) hits[j]=-999.;
2295   //
2296   // --- This part is for no shower developement in beam pipe, TDI, VColl
2297   // If particle interacts with beam pipe, TDI, VColl -> return
2298   if(fNoShower==1 && ((gMC->CurrentMedium() == fMedSensPI) || (gMC->CurrentMedium() == fMedSensTDI) ||
2299      (gMC->CurrentMedium() == fMedSensVColl || (gMC->CurrentMedium() == fMedSensLumi)))){
2300
2301     // If option NoShower is set -> StopTrack
2302
2303     Int_t ipr = 0;
2304       gMC->TrackPosition(s[0],s[1],s[2]);
2305       if(gMC->CurrentMedium() == fMedSensPI){
2306         knamed = gMC->CurrentVolName();
2307         if(!strncmp(knamed,"YMQ",3)){
2308           if(s[2]<0) fpLostITC += 1;
2309           else fpLostITA += 1;
2310           ipr=1;
2311         }
2312         else if(!strncmp(knamed,"YD1",3)){
2313           if(s[2]<0) fpLostD1C += 1;
2314           else fpLostD1A += 1;
2315           ipr=1;
2316         }
2317       }
2318       else if(gMC->CurrentMedium() == fMedSensTDI){
2319         knamed = gMC->CurrentVolName();
2320         if(!strncmp(knamed,"MD1",3)){
2321           if(s[2]<0) fpLostD1C += 1;
2322           else  fpLostD1A += 1;
2323           ipr=1;
2324         }
2325         else if(!strncmp(knamed,"QTD",3)) fpLostTDI += 1;
2326       }
2327       else if(gMC->CurrentMedium() == fMedSensVColl){
2328         knamed = gMC->CurrentVolName();
2329         if(!strncmp(knamed,"QCVC",4)) fpcVCollC++;
2330         else if(!strncmp(knamed,"QCVA",4))  fpcVCollA++;
2331         ipr=1;
2332       }
2333       //
2334       //gMC->TrackMomentum(p[0], p[1], p[2], p[3]);
2335       //printf("\t Particle: mass = %1.3f, E = %1.3f GeV, pz = %1.2f GeV -> stopped in volume %s\n",
2336       //     gMC->TrackMass(), p[3], p[2], gMC->CurrentVolName());
2337       //
2338       if(ipr!=0){
2339         printf("\n\t **********************************\n");
2340         printf("\t ********** Side C **********\n");
2341         printf("\t # of particles in IT = %d\n",fpLostITC);
2342         printf("\t # of particles in D1 = %d\n",fpLostD1C);
2343         printf("\t # of particles in VColl = %d\n",fpcVCollC);
2344         printf("\t ********** Side A **********\n");
2345         printf("\t # of particles in IT = %d\n",fpLostITA);
2346         printf("\t # of particles in D1 = %d\n",fpLostD1A);
2347         printf("\t # of particles in TDI = %d\n",fpLostTDI);
2348         printf("\t # of particles in VColl = %d\n",fpcVCollA);
2349         printf("\t **********************************\n");
2350       }
2351       gMC->StopTrack();
2352       return;
2353   }
2354
2355   if((gMC->CurrentMedium() == fMedSensZN) || (gMC->CurrentMedium() == fMedSensZP) ||
2356      (gMC->CurrentMedium() == fMedSensGR) || (gMC->CurrentMedium() == fMedSensF1) ||
2357      (gMC->CurrentMedium() == fMedSensF2) || (gMC->CurrentMedium() == fMedSensZEM)){
2358
2359
2360   //Particle coordinates
2361     gMC->TrackPosition(s[0],s[1],s[2]);
2362     for(j=0; j<=2; j++) x[j] = s[j];
2363     hits[0] = x[0];
2364     hits[1] = x[1];
2365     hits[2] = x[2];
2366
2367   // Determine in which ZDC the particle is
2368     knamed = gMC->CurrentVolName();
2369     if(!strncmp(knamed,"ZN",2)){
2370           if(x[2]<0.) vol[0]=1; // ZNC (dimuon side)
2371           else if(x[2]>0.) vol[0]=4; //ZNA
2372     }
2373     else if(!strncmp(knamed,"ZP",2)){
2374           if(x[2]<0.) vol[0]=2; //ZPC (dimuon side)
2375           else if(x[2]>0.) vol[0]=5; //ZPA
2376     }
2377     else if(!strncmp(knamed,"ZE",2)) vol[0]=3; //ZEM
2378
2379   // Determine in which quadrant the particle is
2380     if(vol[0]==1){      //Quadrant in ZNC
2381       // Calculating particle coordinates inside ZNC
2382       xdet[0] = x[0]-fPosZNC[0];
2383       xdet[1] = x[1]-fPosZNC[1];
2384       // Calculating quadrant in ZN
2385       if(xdet[0]<=0.){
2386         if(xdet[1]<=0.) vol[1]=1;
2387         else vol[1]=3;
2388       }
2389       else if(xdet[0]>0.){
2390         if(xdet[1]<=0.) vol[1]=2;
2391         else vol[1]=4;
2392       }
2393     }
2394
2395     else if(vol[0]==2){ //Quadrant in ZPC
2396       // Calculating particle coordinates inside ZPC
2397       xdet[0] = x[0]-fPosZPC[0];
2398       xdet[1] = x[1]-fPosZPC[1];
2399       if(xdet[0]>=fDimZP[0])  xdet[0]=fDimZP[0]-0.01;
2400       if(xdet[0]<=-fDimZP[0]) xdet[0]=-fDimZP[0]+0.01;
2401       // Calculating tower in ZP
2402       Float_t xqZP = xdet[0]/(fDimZP[0]/2.);
2403       for(int i=1; i<=4; i++){
2404          if(xqZP>=(i-3) && xqZP<(i-2)){
2405            vol[1] = i;
2406            break;
2407          }
2408       }
2409     }
2410     //
2411     // Quadrant in ZEM: vol[1] = 1 -> particle in 1st ZEM (placed at x = 8.5 cm)
2412     //                  vol[1] = 2 -> particle in 2nd ZEM (placed at x = -8.5 cm)
2413     else if(vol[0] == 3){
2414       if(x[0]>0.){
2415         vol[1] = 1;
2416         // Particle x-coordinate inside ZEM1
2417         xdet[0] = x[0]-fPosZEM[0];
2418       }
2419       else{
2420         vol[1] = 2;
2421         // Particle x-coordinate inside ZEM2
2422         xdet[0] = x[0]+fPosZEM[0];
2423       }
2424       xdet[1] = x[1]-fPosZEM[1];
2425     }
2426     //
2427     else if(vol[0]==4){ //Quadrant in ZNA
2428       // Calculating particle coordinates inside ZNA
2429       xdet[0] = x[0]-fPosZNA[0];
2430       xdet[1] = x[1]-fPosZNA[1];
2431       // Calculating quadrant in ZNA
2432       if(xdet[0]>=0.){
2433         if(xdet[1]<=0.) vol[1]=1;
2434         else vol[1]=3;
2435       }
2436       else if(xdet[0]<0.){
2437         if(xdet[1]<=0.) vol[1]=2;
2438         else vol[1]=4;
2439       }
2440     }
2441     //
2442     else if(vol[0]==5){ //Quadrant in ZPA
2443       // Calculating particle coordinates inside ZPA
2444       xdet[0] = x[0]-fPosZPA[0];
2445       xdet[1] = x[1]-fPosZPA[1];
2446       if(xdet[0]>=fDimZP[0])  xdet[0]=fDimZP[0]-0.01;
2447       if(xdet[0]<=-fDimZP[0]) xdet[0]=-fDimZP[0]+0.01;
2448       // Calculating tower in ZP
2449       Float_t xqZP = -xdet[0]/(fDimZP[0]/2.);
2450       for(int i=1; i<=4; i++){
2451          if(xqZP>=(i-3) && xqZP<(i-2)){
2452            vol[1] = i;
2453            break;
2454          }
2455       }
2456     }
2457     if((vol[1]!=1) && (vol[1]!=2) && (vol[1]!=3) && (vol[1]!=4))
2458       AliError(Form(" WRONG tower for det %d: tow %d with xdet=(%f, %f)\n",
2459                 vol[0], vol[1], xdet[0], xdet[1]));
2460     // Ch. debug
2461     //printf("\t *** det %d vol %d xdet(%f, %f)\n",vol[0], vol[1], xdet[0], xdet[1]);
2462
2463
2464     // Store impact point and kinetic energy of the ENTERING particle
2465
2466     if(gMC->IsTrackEntering()){
2467       //Particle energy
2468       gMC->TrackMomentum(p[0],p[1],p[2],p[3]);
2469       hits[3] = p[3];
2470
2471       // Impact point on ZDC
2472       // X takes into account the LHC x-axis sign
2473       // which is opposite to positive x on detector front face
2474       // for side A detectors (ZNA and ZPA)
2475       if(vol[0]==4 || vol[0]==5){
2476         hits[4] = -xdet[0];
2477       }
2478       else{
2479         hits[4] = xdet[0];
2480       }
2481       hits[5] = xdet[1];
2482       hits[6] = 0;
2483       hits[7] = 0;
2484       hits[8] = 0;
2485       hits[9] = 0;
2486       //
2487       Int_t curTrackN = gAlice->GetMCApp()->GetCurrentTrackNumber();
2488       TParticle *part = gAlice->GetMCApp()->Particle(curTrackN);
2489       hits[10] = part->GetPdgCode();
2490       //printf("\t PDGCode = %d\n", part->GetPdgCode());
2491       //
2492       Int_t imo = part->GetFirstMother();
2493       if(imo>0){
2494         TParticle * pmot = gAlice->GetMCApp()->Particle(imo);
2495         hits[11] = pmot->GetPdgCode();
2496       }
2497       else hits[11]=0;
2498       //
2499       hits[12] = 1.0e09*gMC->TrackTime(); // in ns!
2500       //printf("\t TrackTime = %f\n", hits[12]);
2501
2502       AddHit(curTrackN, vol, hits);
2503
2504       if(fNoShower==1){
2505         if(vol[0]==1){
2506           fnDetectedC += 1;
2507           if(fnDetectedC==1) printf("   ### Particle in ZNC\n\n");
2508         }
2509         else if(vol[0]==2){
2510           fpDetectedC += 1;
2511           if(fpDetectedC==1) printf("   ### Particle in ZPC\n\n");
2512         }
2513         else if(vol[0]==4){
2514           fnDetectedA += 1;
2515           if(fnDetectedA==1) printf("   ### Particle in ZNA\n\n");
2516         }
2517         else if(vol[0]==5){
2518           fpDetectedA += 1;
2519           if(fpDetectedA==1) printf("   ### Particle in ZPA\n\n");
2520         }
2521         //
2522         //printf("\t Pc: x %1.2f y %1.2f z %1.2f  E %1.2f GeV pz = %1.2f GeV in volume %s\n",
2523         //   x[0],x[1],x[3],p[3],p[2],gMC->CurrentVolName());
2524         //
2525         gMC->StopTrack();
2526         return;
2527       }
2528     }
2529
2530     // Particle energy loss
2531     if(gMC->Edep() != 0){
2532       hits[9] = gMC->Edep();
2533       hits[7] = 0.;
2534       hits[8] = 0.;
2535       AddHit(gAlice->GetMCApp()->GetCurrentTrackNumber(), vol, hits);
2536     }
2537   }
2538
2539
2540   // *** Light production in fibres
2541   if((gMC->CurrentMedium() == fMedSensF1) || (gMC->CurrentMedium() == fMedSensF2)){
2542
2543      //Select charged particles
2544      if((destep=gMC->Edep())){
2545
2546        // Particle velocity
2547        Float_t beta = 0.;
2548        gMC->TrackMomentum(p[0],p[1],p[2],p[3]);
2549        Float_t ptot=TMath::Sqrt(p[0]*p[0]+p[1]*p[1]+p[2]*p[2]);
2550        if(p[3] > 0.00001) beta =  ptot/p[3];
2551        else return;
2552        if(beta<0.67)return;
2553        else if((beta>=0.67) && (beta<=0.75)) ibeta = 0;
2554        else if((beta>0.75)  && (beta<=0.85)) ibeta = 1;
2555        else if((beta>0.85)  && (beta<=0.95)) ibeta = 2;
2556        else if(beta>0.95) ibeta = 3;
2557
2558        // Angle between particle trajectory and fibre axis
2559        // 1 -> Momentum directions
2560        um[0] = p[0]/ptot;
2561        um[1] = p[1]/ptot;
2562        um[2] = p[2]/ptot;
2563        gMC->Gmtod(um,ud,2);
2564        // 2 -> Angle < limit angle
2565        Double_t alfar = TMath::ACos(ud[2]);
2566        Double_t alfa = alfar*kRaddeg;
2567        if(alfa>=110.) return;
2568        //
2569        ialfa = Int_t(1.+alfa/2.);
2570
2571        // Distance between particle trajectory and fibre axis
2572        gMC->TrackPosition(s[0],s[1],s[2]);
2573        for(j=0; j<=2; j++){
2574           x[j] = s[j];
2575        }
2576        gMC->Gmtod(x,xdet,1);
2577        if(TMath::Abs(ud[0])>0.00001){
2578          Float_t dcoeff = ud[1]/ud[0];
2579          be = TMath::Abs((xdet[1]-dcoeff*xdet[0])/TMath::Sqrt(dcoeff*dcoeff+1.));
2580        }
2581        else{
2582          be = TMath::Abs(ud[0]);
2583        }
2584
2585        ibe = Int_t(be*1000.+1);
2586
2587        //Looking into the light tables
2588        Float_t charge = gMC->TrackCharge();
2589
2590        if(vol[0]==1 || vol[0]==4) {     // (1)  ZN fibres
2591          if(ibe>fNben) ibe=fNben;
2592          out =  charge*charge*fTablen[ibeta][ialfa][ibe];
2593          nphe = gRandom->Poisson(out);
2594          // Ch. debug
2595          //if(ibeta==3) printf("\t %f \t %f \t %f\n",alfa, be, out);
2596          //printf("\t ibeta = %d, ialfa = %d, ibe = %d -> nphe = %d\n\n",ibeta,ialfa,ibe,nphe);
2597          if(gMC->CurrentMedium() == fMedSensF1){
2598            hits[7] = nphe;      //fLightPMQ
2599            hits[8] = 0;
2600            hits[9] = 0;
2601            AddHit(gAlice->GetMCApp()->GetCurrentTrackNumber(), vol, hits);
2602          }
2603          else{
2604            hits[7] = 0;
2605            hits[8] = nphe;      //fLightPMC
2606            hits[9] = 0;
2607            AddHit(gAlice->GetMCApp()->GetCurrentTrackNumber(), vol, hits);
2608          }
2609        }
2610        else if(vol[0]==2 || vol[0]==5) {// (2) ZP fibres
2611          if(ibe>fNbep) ibe=fNbep;
2612          out =  charge*charge*fTablep[ibeta][ialfa][ibe];
2613          nphe = gRandom->Poisson(out);
2614          if(gMC->CurrentMedium() == fMedSensF1){
2615            hits[7] = nphe;      //fLightPMQ
2616            hits[8] = 0;
2617            hits[9] = 0;
2618            AddHit(gAlice->GetMCApp()->GetCurrentTrackNumber(), vol, hits);
2619          }
2620          else{
2621            hits[7] = 0;
2622            hits[8] = nphe;      //fLightPMC
2623            hits[9] = 0;
2624            AddHit(gAlice->GetMCApp()->GetCurrentTrackNumber(), vol, hits);
2625          }
2626        }
2627        else if((vol[0]==3)) {   // (3) ZEM fibres
2628          if(ibe>fNbep) ibe=fNbep;
2629          out =  charge*charge*fTablep[ibeta][ialfa][ibe];
2630          gMC->TrackPosition(s[0],s[1],s[2]);
2631          Float_t xalic[3];
2632          for(j=0; j<3; j++){
2633             xalic[j] = s[j];
2634          }
2635          // z-coordinate from ZEM front face
2636          // NB-> fPosZEM[2]+fZEMLength = -1000.+2*10.3 = 979.69 cm
2637          Float_t z = -xalic[2]+fPosZEM[2]+2*fZEMLength-xalic[1];
2638          //z = xalic[2]-fPosZEM[2]-fZEMLength-xalic[1]*(TMath::Tan(45.*kDegrad));
2639          //printf("     fPosZEM[2]+2*fZEMLength = %f", fPosZEM[2]+2*fZEMLength);
2640          //
2641          // Parametrization for light guide uniformity
2642          // NEW!!! Light guide tilted @ 51 degrees
2643          Float_t guiPar[4]={0.31,-0.0006305,0.01337,0.8895};
2644          Float_t guiEff = guiPar[0]*(guiPar[1]*z*z+guiPar[2]*z+guiPar[3]);
2645          out = out*guiEff;
2646          nphe = gRandom->Poisson(out);
2647          //printf("     out*guiEff = %f nphe = %d", out, nphe);
2648          if(vol[1] == 1){
2649            hits[7] = 0;
2650            hits[8] = nphe;      //fLightPMC (ZEM1)
2651            hits[9] = 0;
2652            AddHit(gAlice->GetMCApp()->GetCurrentTrackNumber(), vol, hits);
2653          }
2654          else{
2655            hits[7] = nphe;      //fLightPMQ (ZEM2)
2656            hits[8] = 0;
2657            hits[9] = 0;
2658            AddHit(gAlice->GetMCApp()->GetCurrentTrackNumber(), vol, hits);
2659          }
2660        }
2661      }
2662    }
2663 }