VZERO/AliVZEROv7.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 /* $Id$ */
  17
  18 ///////////////////////////////////////////////////////////////////////
  19 //                                                                   //
  20 //  (V-zero) detector  version 7 as designed by the Lyon and         //
  21 //   Mexico groups and Carlos Perez Lara from Pontificia Universidad //
  22 //   Catolica del Peru                                               //
  23 //   All comments should be sent to Brigitte CHEYNIS:                //
  24 //                     b.cheynis@ipnl.in2p3.fr                       //
  25 //   Geometry of April 2006 done with ROOT geometrical modeler       //
  26 //   V0R (now V0C) sits between Z values  -89.5 and  -84.8 cm        //
  27 //   V0L (now V0A) sits between Z values +325.0 and +330.0 cm        //
  28 //   New coordinate system has been implemented in october 2003      //
  29 //   Revision of the V0A part by Lizardo Valencia  in July 2008      //
  30 //                                                                   //
  31 ///////////////////////////////////////////////////////////////////////
  32
  33 // --- Standard libraries ---
  34 #include <Riostream.h>
  35
  36 // --- ROOT libraries ---
  37 #include <TClonesArray.h>
  38 #include <TGeoGlobalMagField.h>
  39 #include <TMath.h>
  40 #include <TParticle.h>
  41 #include <TVirtualMC.h>
  42
  43 #include <TGeoManager.h>
  44 #include <TGeoMatrix.h>
  45 #include <TGeoMaterial.h>
  46 #include <TGeoMedium.h>
  47 #include <TGeoVolume.h>
  48 #include "TGeoTube.h"
  49 #include "TGeoArb8.h"
  50 #include "TGeoCompositeShape.h"
  51
  52 // --- AliRoot header files ---
  53 #include "AliRun.h"
  54 #include "AliMC.h"
  55 #include "AliMagF.h"
  56 #include "AliVZEROLoader.h"
  57 #include "AliVZEROdigit.h"
  58 #include "AliVZEROhit.h"
  59 #include "AliVZEROv7.h"
  60 #include "AliLog.h"
  61 #include "AliTrackReference.h"
  62
  63 ClassImp(AliVZEROv7)
  64
  65 //_____________________________________________________________________________
  66 AliVZEROv7:: AliVZEROv7():AliVZERO(),
  67    fCellId(0),
  68    fTrackPosition(),
  69    fTrackMomentum(),
  70    fV0CHeight1(2.5),
  71    fV0CHeight2(4.4),
  72    fV0CHeight3(7.4),
  73    fV0CHeight4(12.5),
  74    fV0CRMin(4.6),
  75    fV0CRBox(38.0),
  76    fV0CLidThickness(0.30),
  77    fV0CCellThickness(2.00),
  78    fV0CBoxThickness(4.70),
  79    fV0COffsetFibers(1.0),
  80    fV0CLightYield(93.75),
  81    fV0CLightAttenuation(0.05),
  82    fV0CnMeters(15.0),
  83    fV0CFibToPhot(0.3),
  84    fV0AR0(4.2),
  85    fV0AR1(7.6),
  86    fV0AR2(13.8),
  87    fV0AR3(22.7),
  88    fV0AR4(41.3),
  89    fV0AR5(43.3),
  90    fV0AR6(72.6),
  91    fV0AR7(92.0), // Distance from origin to outtermost intersection sector7 and sector8
  92    fV0ASciWd(2.5),
  93    fV0APlaWd(0.5),
  94    fV0APlaAl(0.06),
  95    fV0AOctWd(0.75),
  96    fV0AFraWd(0.2),
  97    fV0AOctH1(1.0),
  98    fV0AOctH2(2.0),
  99    fV0ABasHt(2.0),
 100    fV0AFibRd(0.1),
 101    fV0APlaEx(4.4),
 102    fV0APMBWd(24.6),
 103    fV0APMBHt(22.0),
 104    fV0APMBTh(7.1),
 105    fV0APMBWdW(0.3),
 106    fV0APMBHtW(1.0),
 107    fV0APMBAng(30.0),
 108    fV0APMBThW(0.3),
 109    fV0APMTR1(2.44),
 110    fV0APMTR2(2.54),
 111    fV0APMTR3(2.54),
 112    fV0APMTR4(2.70),
 113    fV0APMTH(10.0),
 114    fV0APMTB(1.0),
 115    fV0AnMeters(fV0AR6*0.01),
 116    fV0ALightYield(93.75),
 117    fV0ALightAttenuation(0.05),
 118    fV0AFibToPhot(0.3),
 119    fVersion(7)
 120 {
 121 // Standard default constructor
 122 }
 123
 124 //_____________________________________________________________________________
 125 AliVZEROv7::AliVZEROv7(const char *name, const char *title):AliVZERO(name,title),
 126    fCellId(0),
 127    fTrackPosition(),
 128    fTrackMomentum(),
 129    fV0CHeight1(2.5),
 130    fV0CHeight2(4.4),
 131    fV0CHeight3(7.4),
 132    fV0CHeight4(12.5),
 133    fV0CRMin(4.6),
 134    fV0CRBox(38.0),
 135    fV0CLidThickness(0.30),
 136    fV0CCellThickness(2.00),
 137    fV0CBoxThickness(4.70),
 138    fV0COffsetFibers(1.0),
 139    fV0CLightYield(93.75),
 140    fV0CLightAttenuation(0.05),
 141    fV0CnMeters(15.0),
 142    fV0CFibToPhot(0.3),
 143    fV0AR0(4.2),
 144    fV0AR1(7.6),
 145    fV0AR2(13.8),
 146    fV0AR3(22.7),
 147    fV0AR4(41.3),
 148    fV0AR5(43.3),
 149    fV0AR6(72.6),
 150    fV0AR7(92.0), // Distance from origin to outtermost intersection of sector7 and sector8
 151    fV0ASciWd(2.5),
 152    fV0APlaWd(0.5),
 153    fV0APlaAl(0.06),
 154    fV0AOctWd(0.75),
 155    fV0AFraWd(0.2),
 156    fV0AOctH1(1.0),
 157    fV0AOctH2(2.0),
 158    fV0ABasHt(2.0),
 159    fV0AFibRd(0.1),
 160    fV0APlaEx(4.4),
 161    fV0APMBWd(24.6),
 162    fV0APMBHt(22.0),
 163    fV0APMBTh(7.1),
 164    fV0APMBWdW(0.3),
 165    fV0APMBHtW(1.0),
 166    fV0APMBAng(30.0),
 167    fV0APMBThW(0.3),
 168    fV0APMTR1(2.44),
 169    fV0APMTR2(2.54),
 170    fV0APMTR3(2.54),
 171    fV0APMTR4(2.70),
 172    fV0APMTH(10.0),
 173    fV0APMTB(1.0),
 174    fV0AnMeters(fV0AR6*0.01),
 175    fV0ALightYield(93.75),
 176    fV0ALightAttenuation(0.05),
 177    fV0AFibToPhot(0.3),
 178    fVersion(7)
 179
 180
 181 {
 182 // Standard constructor for V-zero Detector  version 7
 183
 184   AliDebug(2,"Create VZERO object ");
 185
 186 //  fVersion            =     7;  // version number
 187
 188 //   // V0C Parameters related to geometry: All in cm
 189 //   fV0CHeight1         =    2.5; // height of cell 1
 190 //   fV0CHeight2         =    4.4; // height of cell 2
 191 //   fV0CHeight3         =    7.4; // height of cell 3
 192 //   fV0CHeight4         =   12.5; // height of cell 4
 193 //   fV0CRMin            =    4.6; // inner radius of box
 194 //   fV0CRBox            =   38.0; // outer radius of box
 195 //   fV0CLidThickness    =   0.30; // thickness of Carbon lid
 196 //   fV0CCellThickness   =   2.00; // thickness of elementary cell
 197 //   fV0CBoxThickness    =   4.70; // thickness of V0C Box
 198 //   fV0COffsetFibers    =    1.0; // offset to output fibers
 199 //   // V0C Parameters related to light output
 200 //   fV0CLightYield         =  93.75; // Light yield in BC408 (93.75 eV per photon)
 201 //   fV0CLightAttenuation   =   0.05; // Light attenuation in fiber (0.05 per meter)
 202 //   fV0CnMeters            =   15.0; // Number of meters of clear fibers to PM
 203 //   fV0CFibToPhot          =    0.3; // Attenuation at fiber-photocathode interface
 204 //
 205 //   // V0A Parameters related to geometry: All in cm
 206 //   fV0AR0     =  4.2;  // Radius of hole
 207 //   fV0AR1     =  7.6;  // Maximun radius of 1st cell
 208 //   fV0AR2     = 13.8; // Maximun radius of 2nd cell
 209 //   fV0AR3     = 22.7; // Maximun radius of 3rd cell
 210 //   fV0AR4     = 41.3; // Maximun radius of 4th cell
 211 //   fV0AR5     = 43.3; // Radius circunscrite to innermost octagon
 212 //   fV0AR6     = 68.0; // Radius circunscrite to outtermost octagon
 213 //   fV0ASciWd  =  2.5;  // Scintillator thickness
 214 //   fV0APlaWd  =  0.5;  // Plates thinckness
 215 //   fV0APlaAl  = 0.06; // Plates AlMg3 thinckness
 216 //   fV0AOctWd  = 0.75; // Innermost octagon thickness
 217 //   fV0AOctH1  =  1.0;  // Height of innermost octagon
 218 //   fV0AOctH2  =  2.0;  // Height of outtermost octagon
 219 //   fV0AFibRd  =  0.1;  // Radius of Fiber
 220 //   fV0AFraWd  =  0.2;  // Support Frame thickness
 221 //   fV0APMBWd  = 24.6;  // Width of PM Box
 222 //   fV0APMBHt  = 22.0;  // Height of PM Box
 223 //   fV0APMBTh  =  7.1;  // Thickness of PM Box
 224 //   fV0APMBWdW =  0.3;  // Thickness of PM Box Side1 Wall
 225 //   fV0APMBHtW =  1.0;  // Thickness of PM Box Side2 Wall
 226 //   fV0APMBThW =  0.3;  // Thickness of PM Box Top Wall
 227 //   fV0APMBAng = 30.0;  // Angle between PM Box and Support
 228 //   fV0APMTR1  = 2.44;  // PMT Glass
 229 //   fV0APMTR2  = 2.54;  // PMT Glass
 230 //   fV0APMTR3  = 2.54;  // PMT Cover
 231 //   fV0APMTR4  = 2.70;  // PMT Cover
 232 //   fV0APMTH   = 10.0;  // PMT Height
 233 //   fV0APMTB   =  1.0;  // PMT Basis
 234 //   fV0APlaEx  =  4.4;  // Plates Extension height
 235 //   fV0ABasHt  =  2.0;  // Basis Height
 236 //   // V0A Parameters related to light output
 237 //   fV0ALightYield         =  93.75;      // Light yield in BC404
 238 //   fV0ALightAttenuation   =   0.05;      // Light attenuation in WLS fiber, per meter
 239 //   fV0AnMeters            = fV0AR6*0.01; // Tentative value, in meters
 240 //   fV0AFibToPhot          =    0.3;      // Attenuation at fiber-photocathode interface
 241 }
 242
 243 //_____________________________________________________________________________
 244 void AliVZEROv7::CreateGeometry()
 245 {
 246 // Constructs TGeo geometry
 247
 248   AliDebug(2,"VZERO ConstructGeometry");
 249   TGeoVolume *top = gGeoManager->GetVolume("ALIC");
 250
 251   ///////////////////////////////////////////////////////////////////////////
 252   // Construct the geometry of V0C Detector. Brigitte CHEYNIS
 253
 254     const int kColorVZERO  = kGreen;
 255     TGeoMedium *medV0CAlu = gGeoManager->GetMedium("VZERO_V0CAlu");
 256     TGeoMedium *medV0CCar = gGeoManager->GetMedium("VZERO_V0CCar");
 257     TGeoMedium *medV0CSci = gGeoManager->GetMedium("VZERO_V0CSci");
 258     TGeoVolume *v0RI = new TGeoVolumeAssembly("V0RI");
 259     Float_t heightRight, r4Right;
 260     Float_t zdet = 90.0 - 0.5 - fV0CBoxThickness/2.0;
 261     heightRight  = fV0CHeight1 + fV0CHeight2 + fV0CHeight3 + fV0CHeight4;
 262     r4Right      = fV0CRMin + heightRight + 3.0*0.2; // 3 spacings of 2mm between rings
 263
 264     // Creation of  carbon lids (3.0 mm thick) to keep V0C box shut :
 265     Float_t   partube[3];
 266     partube[0] =   fV0CRMin;
 267     partube[1] =   fV0CRBox;
 268     partube[2] =   fV0CLidThickness/2.0;
 269     TGeoTube   *sV0CA = new TGeoTube("V0CA", partube[0], partube[1], partube[2]);
 270     TGeoVolume *v0CA  = new TGeoVolume("V0CA",sV0CA,medV0CCar);
 271     TGeoTranslation *tr2 = new TGeoTranslation(0.,0., fV0CBoxThickness/2.0-partube[2]);
 272     TGeoTranslation *tr3 = new TGeoTranslation(0.,0.,-fV0CBoxThickness/2.0+partube[2]);
 273     v0RI->AddNode(v0CA,1,tr2);
 274     v0RI->AddNode(v0CA,2,tr3);
 275     v0CA->SetLineColor(kYellow);
 276
 277     // Creation of aluminum rings 3.0 mm thick to maintain the v0RI pieces :
 278     partube[0] =   fV0CRMin - 0.3;
 279     partube[1] =   fV0CRMin;
 280     partube[2] =   fV0CBoxThickness/2.0;
 281     TGeoTube   *sV0IR = new TGeoTube("V0IR", partube[0], partube[1], partube[2]);
 282     TGeoVolume *v0IR  = new TGeoVolume("V0IR",sV0IR,medV0CAlu);
 283     v0RI->AddNode(v0IR,1,0);
 284     v0IR->SetLineColor(kYellow);
 285     partube[0] =   fV0CRBox;
 286     partube[1] =   fV0CRBox + 0.3;
 287     partube[2] =   fV0CBoxThickness/2.0;
 288     TGeoTube   *sV0ER = new TGeoTube("V0ER", partube[0], partube[1], partube[2]);
 289     TGeoVolume *v0ER  = new TGeoVolume("V0ER",sV0ER,medV0CAlu);
 290     v0RI->AddNode(v0ER,1,0);
 291     v0ER->SetLineColor(kYellow);
 292
 293     // Creation of assembly V0R0 of scintillator cells within one sector
 294     TGeoVolume *v0R0 = new TGeoVolumeAssembly("V0R0");
 295
 296     // Elementary cell of ring 1  - right part - :
 297     // (cells of ring 1 will be shifted by 2.0 cm backwards to output fibers)
 298     Float_t   r1Right =  fV0CRMin + fV0CHeight1;
 299     Float_t   offset  = fV0CBoxThickness/2.0 - fV0CLidThickness - fV0CCellThickness/2.0;
 300     Float_t   partubs[5];
 301     partubs[0]     =  fV0CRMin;
 302     partubs[1]     =  r1Right;
 303     partubs[2]     =  fV0CCellThickness/2.0;
 304     partubs[3]     =  90.0-22.5;
 305     partubs[4]     = 135.0-22.5;
 306     TGeoTubeSeg *sV0R1 = new TGeoTubeSeg("V0R1", partubs[0], partubs[1], partubs[2],
 307                                          partubs[3], partubs[4]);
 308     TGeoVolume  *v0R1  = new TGeoVolume("V0R1",sV0R1,medV0CSci);
 309     TGeoTranslation *tr4 = new TGeoTranslation(0.,0.,-offset);
 310     v0R0->AddNode(v0R1,1,tr4);
 311     v0R1->SetLineColor(kColorVZERO);
 312
 313     // Elementary cell of ring 2 - right part - :
 314     // (cells of ring 2 will be shifted by 1.0 cm backwards to output fibers)
 315     Float_t   r2Right  =  r1Right + fV0CHeight2;
 316     partubs[0]     =  r1Right;  //  must be equal to 7.1
 317     partubs[1]     =  r2Right;  //  must be equal to 11.5
 318     TGeoTubeSeg *sV0R2 = new TGeoTubeSeg("V0R2", partubs[0], partubs[1], partubs[2],
 319                                          partubs[3], partubs[4]);
 320     TGeoVolume  *v0R2  = new TGeoVolume("V0R2",sV0R2,medV0CSci);
 321     TGeoTranslation *tr5 = new TGeoTranslation(0.0,0.2,-offset + fV0COffsetFibers);
 322     v0R0->AddNode(v0R2,1,tr5);
 323     v0R2->SetLineColor(kColorVZERO);
 324
 325     // Ring 3 - right part -  :
 326     r2Right  =  r2Right + 0.2;
 327     Float_t   r3Right  =  r2Right + fV0CHeight3;
 328     partubs[0]     =  r2Right;  //  must be equal to 11.7
 329     partubs[1]     =  r3Right;  //  must be equal to 19.1
 330     partubs[3]     =  90.0-22.5;
 331     partubs[4]     = 112.5-22.5;
 332     TGeoTubeSeg *sV0R3 = new TGeoTubeSeg("V0R3", partubs[0], partubs[1], partubs[2],
 333                                          partubs[3], partubs[4]);
 334     TGeoVolume  *v0R3  = new TGeoVolume("V0R3",sV0R3,medV0CSci);
 335     TGeoTranslation *tr6 = new TGeoTranslation(0.,0.2,-offset + 2.0*fV0COffsetFibers);
 336     v0R0->AddNode(v0R3,1,tr6);
 337     v0R3->SetLineColor(kColorVZERO);
 338     partubs[3]     = 112.5-22.5;
 339     partubs[4]     = 135.0-22.5;
 340     TGeoTubeSeg *sV0R4 = new TGeoTubeSeg("V0R4", partubs[0], partubs[1], partubs[2],
 341                                          partubs[3], partubs[4]);
 342     TGeoVolume  *v0R4  = new TGeoVolume("V0R4",sV0R4,medV0CSci);
 343     v0R0->AddNode(v0R4,1,tr6);
 344     v0R4->SetLineColor(kColorVZERO);
 345
 346     // Ring 4 - right part -  :
 347     Float_t x = TMath::ATan(3.5/257.5) * ((180./TMath::Pi()));
 348     r3Right = r3Right + 0.2 + 0.2;   // + 0.2 because no shift in translation here !!
 349     partubs[0]     =  r3Right;  //  must be equal to 19.5
 350     partubs[1]     =  r4Right;  //  must be equal to 32.0
 351     partubs[3]     =  90.0-22.5+x;
 352     partubs[4]     = 112.5-22.5-x;
 353     TGeoTubeSeg *sV0R5 = new TGeoTubeSeg("V0R5", partubs[0], partubs[1], partubs[2],
 354                                          partubs[3], partubs[4]);
 355     TGeoVolume  *v0R5  = new TGeoVolume("V0R5",sV0R5,medV0CSci);
 356     TGeoTranslation *tr7 = new TGeoTranslation(0.,0.0,-offset + 2.0*fV0COffsetFibers);
 357     v0R0->AddNode(v0R5,1,tr7);
 358     v0R5->SetLineColor(kColorVZERO);
 359     partubs[3]     = 112.5-22.5+x;
 360     partubs[4]     = 135.0-22.5-x;
 361     TGeoTubeSeg *sV0R6 = new TGeoTubeSeg("V0R6", partubs[0], partubs[1], partubs[2],
 362                                          partubs[3], partubs[4]);
 363     TGeoVolume  *v0R6  = new TGeoVolume("V0R6",sV0R6,medV0CSci);
 364     v0R0->AddNode(v0R6,1,tr7);
 365     v0R6->SetLineColor(kColorVZERO);
 366     Float_t  phi;
 367     Float_t  phiDeg= 180./4.;
 368     Int_t    nsecR = 1;     // number of sectors in right part of V0
 369     for (phi = 22.5; phi < 360.0; phi = phi + phiDeg) {
 370       TGeoRotation  *rot1 = new TGeoRotation("rot1", 90.0, +phi, 90., 90.+phi, 0.0, 0.0 );
 371       v0RI->AddNode(v0R0,nsecR,rot1);
 372       nsecR++;
 373     }
 374
 375   ///////////////////////////////////////////////////////////////////////////
 376   // Construct the geometry of V0A Detector. Carlos PEREZ, PUCP
 377   // Revision by Lizardo VALENCIA, UNAM Mexico in July 2008
 378
 379     const int kV0AColorSci   = 5;
 380     const int kV0AColorPlaIn = 3;
 381     const int kV0AColorPlaOu = 41;
 382     const int kV0AColorOct   = 7;
 383     const int kV0AColorFra   = 6;
 384     const int kV0AColorFib   = 11;
 385     const int kV0AColorPMG   = 1;
 386     const int kV0AColorPMA   = 2;
 387     TGeoMedium *medV0ASci = gGeoManager->GetMedium("VZERO_V0ASci");
 388     TGeoMedium *medV0APlaIn = gGeoManager->GetMedium("VZERO_V0APlaIn");
 389     TGeoMedium *medV0APlaOu = gGeoManager->GetMedium("VZERO_V0APlaOu");
 390     TGeoMedium *medV0ASup = gGeoManager->GetMedium("VZERO_V0ALuc");
 391     TGeoMedium *medV0AFra = gGeoManager->GetMedium("VZERO_V0ALuc");
 392     TGeoMedium *medV0AFib = gGeoManager->GetMedium("VZERO_V0AFib");
 393     TGeoMedium *medV0APMGlass = gGeoManager->GetMedium("VZERO_V0APMG");
 394     TGeoMedium *medV0APMAlum = gGeoManager->GetMedium("VZERO_V0APMA");
 395     double pi = TMath::Pi();
 396     double sin225   = TMath::Sin(pi/8.);
 397     double cos225   = TMath::Cos(pi/8.);
 398     double sin45    = TMath::Sin(pi/4.); // lucky: Sin45=Cos45
 399     double cos45    = TMath::Cos(pi/4.);
 400     double v0APts[16];
 401     double sin654   = TMath::Sin(1.14);
 402     double cos654   = TMath::Cos(1.14);
 403
 404     ////////////////////////////
 405     /// Definition sector 1
 406     TGeoVolume *v0ASec = new TGeoVolumeAssembly("V0ASec");
 407
 408     /// For boolean sustraction
 409     double preShape = 0.2;
 410     for (int i=0;i<2;i++) {
 411       v0APts[0+8*i] = fV0AR0-fV0AFraWd/2.-preShape;  v0APts[1+8*i] = -preShape;
 412       v0APts[2+8*i] = fV0AR0-fV0AFraWd/2.-preShape;  v0APts[3+8*i] = fV0AFraWd/2.;
 413       v0APts[4+8*i] = fV0AR4+fV0AFraWd/2.+preShape;  v0APts[5+8*i] = fV0AFraWd/2.;
 414       v0APts[6+8*i] = fV0AR4+fV0AFraWd/2.+preShape;  v0APts[7+8*i] = -preShape;
 415     }
 416     new TGeoArb8("sV0ACha1",fV0ASciWd/1.5,v0APts);
 417     for (int i=0;i<2;i++) {
 418       v0APts[0+8*i] = fV0AR0*sin45-preShape;
 419       v0APts[1+8*i] = (fV0AR0-fV0AFraWd)*sin45-preShape;
 420       v0APts[2+8*i] = (fV0AR0-fV0AFraWd/2.)*sin45-preShape;
 421       v0APts[3+8*i] = (fV0AR0-fV0AFraWd/2.)*sin45;
 422       v0APts[4+8*i] = (fV0AR4+fV0AFraWd/2.)*sin45+preShape;
 423       v0APts[5+8*i] = (fV0AR4+fV0AFraWd/2.)*sin45+2.*preShape;
 424       v0APts[6+8*i] = (fV0AR4+fV0AFraWd)*sin45+preShape;
 425       v0APts[7+8*i] = fV0AR4*sin45+preShape;
 426     }
 427     new TGeoArb8("sV0ACha2", fV0ASciWd/2.+2.*preShape, v0APts);
 428     new TGeoCompositeShape("sV0ACha12","sV0ACha1+sV0ACha2");
 429     new TGeoTube("sV0ANail1SciHole", 0.0, 0.4, 1.65);
 430     TGeoTranslation *pos1 = new TGeoTranslation("pos1", 42.9, 0.51, 0.0);
 431     pos1->RegisterYourself();
 432     new TGeoTube("sV0ANail2SciHole", 0.0, 0.4, 1.65);
 433     TGeoTranslation *pos2 = new TGeoTranslation("pos2", 30.73,29.98,0.0);
 434     pos2->RegisterYourself();
 435     new TGeoCompositeShape("sV0ANailsSciHoles","sV0ANail1SciHole:pos1+sV0ANail2SciHole:pos2");
 436     new TGeoCompositeShape("sV0ACha","sV0ACha12+sV0ANailsSciHoles");
 437     new TGeoTubeSeg("sV0AFicR5", fV0AR4+fV0AFraWd/2., fV0AR4 + fV0AR0, fV0ASciWd/2.+2*preShape, 0, 45);
 438
 439     /// Frame
 440     TGeoVolume *v0AFra = new TGeoVolumeAssembly("V0AFra");
 441     for (int i=0;i<2;i++) {
 442       v0APts[0+8*i] = fV0AR0-fV0AFraWd/2.;  v0APts[1+8*i] = 0.;
 443       v0APts[2+8*i] = fV0AR0-fV0AFraWd/2.;  v0APts[3+8*i] = fV0AFraWd/2.;
 444       v0APts[4+8*i] = fV0AR4+fV0AFraWd/2.;  v0APts[5+8*i] = fV0AFraWd/2.;
 445       v0APts[6+8*i] = fV0AR4+fV0AFraWd/2.;  v0APts[7+8*i] = 0.;
 446     }
 447     TGeoArb8 *sV0AFraB1 = new TGeoArb8("sV0AFraB1",fV0ASciWd/2.,v0APts);
 448     TGeoVolume *v0AFraB1 = new TGeoVolume("V0AFraB1",sV0AFraB1,medV0AFra);
 449     for (int i=0;i<2;i++) {
 450       v0APts[0+8*i] = fV0AR0*sin45;
 451       v0APts[1+8*i] = (fV0AR0-fV0AFraWd)*sin45;
 452       v0APts[2+8*i] = (fV0AR0-fV0AFraWd/2.)*sin45;
 453       v0APts[3+8*i] = (fV0AR0-fV0AFraWd/2.)*sin45;
 454       v0APts[4+8*i] = (fV0AR4+fV0AFraWd/2.)*sin45;
 455       v0APts[5+8*i] = (fV0AR4+fV0AFraWd/2.)*sin45;
 456       v0APts[6+8*i] = (fV0AR4+fV0AFraWd)*sin45;
 457       v0APts[7+8*i] = fV0AR4*sin45;
 458     }
 459     TGeoArb8 *sV0AFraB2 = new TGeoArb8("sV0AFraB2", fV0ASciWd/2., v0APts);
 460     TGeoVolume *v0AFraB2 = new TGeoVolume("V0AFraB2",sV0AFraB2,medV0AFra);
 461     v0AFraB1->SetLineColor(kV0AColorFra); v0AFraB2->SetLineColor(kV0AColorFra);
 462     v0AFra->AddNode(v0AFraB1,1);
 463     v0AFra->AddNode(v0AFraB2,1);
 464     new TGeoTubeSeg( "sV0AFraR1b", fV0AR0-fV0AFraWd/2.,
 465                      fV0AR0+fV0AFraWd/2., fV0ASciWd/2., 0, 45);
 466     new TGeoTubeSeg( "sV0AFraR2b", fV0AR1-fV0AFraWd/2.,
 467                      fV0AR1+fV0AFraWd/2., fV0ASciWd/2., 0, 45);
 468     new TGeoTubeSeg( "sV0AFraR3b", fV0AR2-fV0AFraWd/2.,
 469                      fV0AR2+fV0AFraWd/2., fV0ASciWd/2., 0, 45);
 470     new TGeoTubeSeg( "sV0AFraR4b", fV0AR3-fV0AFraWd/2.,
 471                      fV0AR3+fV0AFraWd/2., fV0ASciWd/2., 0, 45);
 472     new TGeoTubeSeg( "sV0AFraR5b", fV0AR4-fV0AFraWd/2.,
 473                      fV0AR4+fV0AFraWd/2., fV0ASciWd/2., 0, 45);
 474     TGeoCompositeShape *sV0AFraR1 = new TGeoCompositeShape("sV0AFraR1","sV0AFraR1b-sV0ACha");
 475     TGeoCompositeShape *sV0AFraR2 = new TGeoCompositeShape("sV0AFraR2","sV0AFraR2b-sV0ACha");
 476     TGeoCompositeShape *sV0AFraR3 = new TGeoCompositeShape("sV0AFraR3","sV0AFraR3b-sV0ACha");
 477     TGeoCompositeShape *sV0AFraR4 = new TGeoCompositeShape("sV0AFraR4","sV0AFraR4b-sV0ACha");
 478     TGeoCompositeShape *sV0AFraR5 = new TGeoCompositeShape("sV0AFraR5","sV0AFraR5b-sV0ACha");
 479     TGeoVolume *v0AFraR1 = new TGeoVolume("V0AFraR1",sV0AFraR1,medV0AFra);
 480     TGeoVolume *v0AFraR2 = new TGeoVolume("V0AFraR2",sV0AFraR2,medV0AFra);
 481     TGeoVolume *v0AFraR3 = new TGeoVolume("V0AFraR3",sV0AFraR3,medV0AFra);
 482     TGeoVolume *v0AFraR4 = new TGeoVolume("V0AFraR4",sV0AFraR4,medV0AFra);
 483     TGeoVolume *v0AFraR5 = new TGeoVolume("V0AFraR5",sV0AFraR5,medV0AFra);
 484     v0AFraR1->SetLineColor(kV0AColorFra); v0AFraR2->SetLineColor(kV0AColorFra);
 485     v0AFraR3->SetLineColor(kV0AColorFra); v0AFraR4->SetLineColor(kV0AColorFra);
 486     v0AFraR5->SetLineColor(kV0AColorFra);
 487     v0AFra->AddNode(v0AFraR1,1);
 488     v0AFra->AddNode(v0AFraR2,1);
 489     v0AFra->AddNode(v0AFraR3,1);
 490     v0AFra->AddNode(v0AFraR4,1);
 491     v0AFra->AddNode(v0AFraR5,1);
 492     v0ASec->AddNode(v0AFra,1);
 493
 494     /// Sensitive scintilator
 495     TGeoVolume *v0ASci = new TGeoVolumeAssembly("V0ASci");
 496     new TGeoTubeSeg( "sV0AR1b", fV0AR0+fV0AFraWd/2.,
 497                      fV0AR1-fV0AFraWd/2., fV0ASciWd/2., 0, 45);
 498     new TGeoTubeSeg( "sV0AR2b", fV0AR1+fV0AFraWd/2.,
 499                      fV0AR2-fV0AFraWd/2., fV0ASciWd/2., 0, 45);
 500     new TGeoTubeSeg( "sV0AR3b", fV0AR2+fV0AFraWd/2.,
 501                      fV0AR3-fV0AFraWd/2., fV0ASciWd/2., 0, 45);
 502     new TGeoTubeSeg( "sV0AR4b", fV0AR3+fV0AFraWd/2.,
 503                      fV0AR4-fV0AFraWd/2., fV0ASciWd/2., 0, 45);
 504     TGeoCompositeShape *sV0AR1 = new TGeoCompositeShape("sV0AR1","sV0AR1b-sV0ACha");
 505     TGeoCompositeShape *sV0AR2 = new TGeoCompositeShape("sV0AR2","sV0AR2b-sV0ACha");
 506     TGeoCompositeShape *sV0AR3 = new TGeoCompositeShape("sV0AR3","sV0AR3b-sV0ACha");
 507     TGeoCompositeShape *sV0AR4 = new TGeoCompositeShape("sV0AR4","sV0AR4b-sV0ACha");
 508     TGeoVolume *v0L1 = new TGeoVolume("V0L1",sV0AR1,medV0ASci);
 509     TGeoVolume *v0L2 = new TGeoVolume("V0L2",sV0AR2,medV0ASci);
 510     TGeoVolume *v0L3 = new TGeoVolume("V0L3",sV0AR3,medV0ASci);
 511     TGeoVolume *v0L4 = new TGeoVolume("V0L4",sV0AR4,medV0ASci);
 512     v0L1->SetLineColor(kV0AColorSci); v0L2->SetLineColor(kV0AColorSci);
 513     v0L3->SetLineColor(kV0AColorSci); v0L4->SetLineColor(kV0AColorSci);
 514     v0ASec->AddNode(v0L1,1);
 515     v0ASec->AddNode(v0L2,1);
 516     v0ASec->AddNode(v0L3,1);
 517     v0ASec->AddNode(v0L4,1);
 518
 519     /// Segment of octagon
 520     for (int i=0;i<2;i++) {
 521     v0APts[0+8*i] =  fV0AR6-fV0AOctH2;  v0APts[1+8*i] = 0.;
 522       v0APts[2+8*i] = (fV0AR6-fV0AOctH2)*sin45;  v0APts[3+8*i] = (fV0AR6-fV0AOctH2)*sin45;
 523       v0APts[4+8*i] = fV0AR6*sin45;             v0APts[5+8*i] = fV0AR6*sin45;
 524       v0APts[6+8*i] = fV0AR6;                   v0APts[7+8*i] = 0.;
 525     }
 526     TGeoArb8 *sV0AOct2 = new TGeoArb8("sV0AOct2", (fV0ASciWd+2*fV0AOctWd)/2., v0APts);
 527     TGeoVolume *v0AOct2 = new TGeoVolume("V0AOct2", sV0AOct2,medV0ASup);
 528     v0AOct2->SetLineColor(kV0AColorOct);
 529     TGeoVolume *v0ASup = new TGeoVolumeAssembly("V0ASup");
 530     v0ASup->AddNode(v0AOct2,1);
 531     v0ASec->AddNode(v0ASup,1);
 532
 533     //Bunch of fibers
 534     v0APts[ 0] = v0APts[ 2] = -13.0;
 535     v0APts[ 1] = v0APts[ 7] = (fV0ASciWd+fV0AOctWd)/2.-0.01;
 536     v0APts[ 3] = v0APts[ 5] = (fV0ASciWd+fV0AOctWd)/2.+0.01;
 537     v0APts[ 4] = v0APts[ 6] = +13.0;
 538     v0APts[ 8] = v0APts[10] = -10.0;
 539     v0APts[ 9] = v0APts[15] = 0.;
 540     v0APts[11] = v0APts[13] = 0.25;
 541     v0APts[12] = v0APts[14] = +10.0;
 542     new TGeoArb8("sV0AFib1", 11.5, v0APts);
 543     TGeoRotation *rot = new TGeoRotation("rot");
 544     rot->RotateX(-90);
 545     rot->RotateZ(-90.+22.5);
 546     TGeoCombiTrans *fibpos1 = new TGeoCombiTrans("fibpos1", (fV0AR6-fV0AOctH2+fV0AR5)*cos225/2. - 3.8, (fV0AR6-fV0AOctH2+fV0AR5)*sin225/2.-1.8, 0, rot);
 547     fibpos1->RegisterYourself();
 548     TGeoCompositeShape *sV0AFib1Hole = new TGeoCompositeShape("sV0AFib1Hole","sV0AFib1:fibpos1-sV0AFicR5");
 549     TGeoVolume *v0AFib1Hole = new TGeoVolume("V0AFib1Hole",sV0AFib1Hole,medV0AFib);
 550     v0AFib1Hole->SetLineColor(kV0AColorFib);
 551     new TGeoArb8("sV0AFib2", 11.5, v0APts);
 552     rot = new TGeoRotation("rot");
 553     rot->RotateX(-90);
 554     rot->RotateY(180);
 555     rot->RotateZ(-90.+22.5);
 556     TGeoCombiTrans *fibpos2 = new TGeoCombiTrans("fibpos2", (fV0AR6-fV0AOctH2+fV0AR5)*cos225/2. - 3.8, (fV0AR6-fV0AOctH2+fV0AR5)*sin225/2.-1.8, 0, rot);
 557     fibpos2->RegisterYourself();
 558     TGeoCompositeShape *sV0AFib2Hole = new TGeoCompositeShape("sV0AFib2Hole","sV0AFib2:fibpos2-sV0AFicR5");
 559     TGeoVolume *v0AFib2Hole = new TGeoVolume("V0AFib2Hole",sV0AFib2Hole,medV0AFib);
 560     v0AFib2Hole->SetLineColor(kV0AColorFib);
 561     TGeoVolume *v0AFib = new TGeoVolumeAssembly("V0AFib");
 562     v0AFib->AddNode(v0AFib1Hole,1);
 563     v0AFib->AddNode(v0AFib2Hole,1);
 564     v0ASec->AddNode(v0AFib,1);
 565
 566      /// Plates
 567     new TGeoTube("sV0ANail1PlaInHole", 0.0, 0.4, (fV0APlaWd-2*fV0APlaAl)/2.);
 568     new TGeoTube("sV0ANail2PlaInHole", 0.0, 0.4, (fV0APlaWd-2*fV0APlaAl)/2.);
 569     new TGeoCompositeShape("sV0ANailsPlaInHoles","sV0ANail1PlaInHole:pos1+sV0ANail2PlaInHole:pos2");
 570     new TGeoTube("sV0ANail1PlaOuHole", 0.0, 0.4, (fV0APlaAl)/2.);
 571     new TGeoTube("sV0ANail2PlaOuHole", 0.0, 0.4, (fV0APlaAl)/2.);
 572     new TGeoCompositeShape("sV0ANailsPlaOuHoles","sV0ANail1PlaOuHole:pos1+sV0ANail2PlaOuHole:pos2");
 573     for (int i=0;i<2;i++) {
 574       v0APts[0+8*i] = fV0AR0;                   v0APts[1+8*i] = 0.;
 575       v0APts[2+8*i] = fV0AR0*sin45;             v0APts[3+8*i] = fV0AR0*sin45;
 576       v0APts[4+8*i] = fV0AR6 * sin45;   v0APts[5+8*i] = fV0AR6*sin45;
 577       v0APts[6+8*i] = fV0AR6;           v0APts[7+8*i] = 0.;
 578     }
 579     new TGeoArb8("sV0APlaIn", (fV0APlaWd-2*fV0APlaAl)/2., v0APts);
 580     TGeoCompositeShape *sV0APlaInNailsHoles = new TGeoCompositeShape("sV0APlaInNailsHoles","sV0APlaIn-sV0ANailsPlaInHoles");
 581     TGeoVolume *v0APlaInNailsHoles = new TGeoVolume("V0APlaInNailsHoles", sV0APlaInNailsHoles, medV0APlaIn);
 582     new TGeoArb8("sV0APlaOu", fV0APlaAl/2., v0APts);
 583     TGeoCompositeShape *sV0APlaOuNailsHoles = new TGeoCompositeShape("sV0APlaOuNailsHoles","sV0APlaOu-sV0ANailsPlaOuHoles");
 584     TGeoVolume *v0APlaOuNailsHoles = new TGeoVolume("V0APlaOuNailsHoles", sV0APlaOuNailsHoles, medV0APlaOu);
 585     v0APlaInNailsHoles->SetLineColor(kV0AColorPlaIn); v0APlaOuNailsHoles->SetLineColor(kV0AColorPlaOu);
 586     TGeoVolume *v0APla = new TGeoVolumeAssembly("V0APla");
 587     v0APla->AddNode(v0APlaInNailsHoles,1);
 588     v0APla->AddNode(v0APlaOuNailsHoles,1,new TGeoTranslation(0,0,(fV0APlaWd-fV0APlaAl)/2.));
 589     v0APla->AddNode(v0APlaOuNailsHoles,2,new TGeoTranslation(0,0,-(fV0APlaWd-fV0APlaAl)/2.));
 590     v0ASec->AddNode(v0APla,1,new TGeoTranslation(0,0,(fV0ASciWd+2*fV0AOctWd+fV0APlaWd)/2.));
 591     v0ASec->AddNode(v0APla,2,new TGeoTranslation(0,0,-(fV0ASciWd+2*fV0AOctWd+fV0APlaWd)/2.));
 592
 593      /// Non-sensitive scintilator
 594     new TGeoTubeSeg("sV0AR5S2", fV0AR4+fV0AFraWd/2., fV0AR4 + fV0AR0, fV0ASciWd/2.+2*preShape, 0, 45);
 595     TGeoCompositeShape *sV0AR5 = new TGeoCompositeShape("V0AR5","sV0AR5S2 - sV0ACha");
 596     TGeoVolume *v0AR5 = new TGeoVolume("V0AR5",sV0AR5,medV0ASci);
 597     v0AR5->SetLineColor(kV0AColorSci);
 598     v0ASci->AddNode(v0AR5,1);
 599     v0ASec->AddNode(v0ASci,1);
 600
 601     /// PMBox
 602     TGeoVolume* v0APM = new TGeoVolumeAssembly("V0APM");
 603     new TGeoBBox("sV0APMB1", fV0APMBWd/2., fV0APMBHt/2., fV0APMBTh/2.);
 604     new TGeoBBox("sV0APMB2", fV0APMBWd/2.-fV0APMBWdW, fV0APMBHt/2.-fV0APMBHtW, fV0APMBTh/2.-fV0APMBThW);
 605     TGeoCompositeShape *sV0APMB = new TGeoCompositeShape("sV0APMB","sV0APMB1-sV0APMB2");
 606     TGeoVolume *v0APMB = new TGeoVolume("V0APMB",sV0APMB, medV0APMAlum);
 607     v0APMB->SetLineColor(kV0AColorPMA);
 608     v0APM->AddNode(v0APMB,1);
 609
 610     /// PMTubes
 611     TGeoTube *sV0APMT1 = new TGeoTube("sV0APMT1", fV0APMTR1, fV0APMTR2, fV0APMTH/2.);
 612     TGeoVolume *v0APMT1 = new TGeoVolume("V0APMT1", sV0APMT1, medV0APMGlass);
 613     TGeoTube *sV0APMT2 = new TGeoTube("sV0APMT2", fV0APMTR3, fV0APMTR4, fV0APMTH/2.);
 614     TGeoVolume *v0APMT2 = new TGeoVolume("V0APMT2", sV0APMT2, medV0APMAlum);
 615     TGeoVolume *v0APMT = new TGeoVolumeAssembly("V0APMT");
 616     TGeoTube *sV0APMTT = new TGeoTube("sV0APMTT", 0., fV0APMTR4, fV0APMTB/2.);
 617     TGeoVolume *v0APMTT = new TGeoVolume("V0APMTT", sV0APMTT, medV0APMAlum);
 618     v0APMT1->SetLineColor(kV0AColorPMG);
 619     v0APMT2->SetLineColor(kV0AColorPMA);
 620     v0APMTT->SetLineColor(kV0AColorPMA);
 621     rot = new TGeoRotation("rot", 90, 0, 180, 0, 90, 90);
 622     v0APMT->AddNode(v0APMT1,1,rot);
 623     v0APMT->AddNode(v0APMT2,1,rot);
 624     v0APMT->AddNode(v0APMTT,1,new TGeoCombiTrans(0,-(fV0APMTH+fV0APMTB)/2.,0,rot));
 625     double autoShift = (fV0APMBWd-2*fV0APMBWdW)/4.;
 626     v0APM->AddNode(v0APMT, 1, new TGeoTranslation(-1.5*autoShift, 0, 0));
 627     v0APM->AddNode(v0APMT, 2, new TGeoTranslation(-0.5*autoShift, 0, 0));
 628     v0APM->AddNode(v0APMT, 3, new TGeoTranslation(+0.5*autoShift, 0, 0));
 629     v0APM->AddNode(v0APMT, 4, new TGeoTranslation(+1.5*autoShift, 0, 0));
 630
 631     /// PM
 632     rot = new TGeoRotation("rot");
 633     rot->RotateX(90-fV0APMBAng);
 634     rot->RotateZ(-90.+22.5);
 635     double cosAngPMB = TMath::Cos(fV0APMBAng*TMath::DegToRad());
 636     double sinAngPMB = TMath::Sin(fV0APMBAng*TMath::DegToRad());
 637     double shiftZ = fV0APMBHt/2. * cosAngPMB
 638       -   ( fV0ASciWd + 2 * fV0AOctWd + 2 * fV0APlaWd )/2.   -   fV0APMBTh/2. * sinAngPMB;
 639     double shiftR = fV0AR6  +  fV0AOctH2 + fV0APlaAl;
 640     v0ASec->AddNode(v0APM,1, new TGeoCombiTrans( shiftR*cos225+1.07, shiftR*sin225, shiftZ, rot));
 641
 642     // Aluminium nails
 643     TGeoTube *sV0ANail1 = new TGeoTube("sV0ANail1", 0.0, 0.4, 5.09/2.);
 644     TGeoVolume *v0ANail1 = new TGeoVolume("V0ANail1", sV0ANail1, medV0APMAlum);
 645     v0ANail1->SetLineColor(kV0AColorPMA);// this is the color for aluminium
 646     v0ASec->AddNode(v0ANail1,1,new TGeoTranslation(42.9, 0.51, 0.0));
 647     TGeoTube *sV0ANail2 = new TGeoTube("sV0ANail2", 0.0, 0.4, 5.09/2.);
 648     TGeoVolume *v0ANail2 = new TGeoVolume("V0ANail2", sV0ANail2, medV0APMAlum);
 649     v0ANail2->SetLineColor(kV0AColorPMA);
 650     v0ASec->AddNode(v0ANail2,1,new TGeoTranslation(30.73,29.98,0.0));
 651
 652     /// Replicate sectors and adding sector to v0LE volume
 653     TGeoVolume *v0LE = new TGeoVolumeAssembly("V0LE");
 654     for(int i=0; i<4; i++) {
 655        TGeoRotation *rotation = new TGeoRotation("rotation", 90., i*45., 90., 90.+i*45., 0., 0.);
 656        v0LE->AddNode(v0ASec,i+1,rotation);
 657     }
 658
 659     //Upper supports
 660     for (int i=0;i<2;i++){
 661     v0APts[0+8*i] = 0.0;            v0APts[1+8*i] = 45.5;
 662     v0APts[2+8*i] = 0.0;            v0APts[3+8*i] = 70.4;
 663     v0APts[4+8*i] = 4.0;            v0APts[5+8*i] = 68.9;
 664     v0APts[6+8*i] = 4.0;            v0APts[7+8*i] = 45.5;
 665     }
 666     TGeoArb8 *sV0ASuppur = new TGeoArb8("sV0ASuppur", 1.65, v0APts);
 667     TGeoVolume *v0ASuppur = new TGeoVolume("V0ASuppur", sV0ASuppur, medV0ASup);
 668     v0ASuppur->SetLineColor(kV0AColorOct);
 669     v0LE->AddNode(v0ASuppur,1);
 670     for (int i=0;i<2;i++){
 671     v0APts[0+8*i] = -0.0;           v0APts[1+8*i] = 45.5;
 672     v0APts[2+8*i] = -0.0;           v0APts[3+8*i] = 70.4;
 673     v0APts[4+8*i] = -4.0;           v0APts[5+8*i] = 68.9;
 674     v0APts[6+8*i] = -4.0;           v0APts[7+8*i] = 45.5;
 675     }
 676     TGeoArb8 *sV0ASuppul = new TGeoArb8("sV0ASuppul", 1.65, v0APts);
 677     TGeoVolume *v0ASuppul = new TGeoVolume("V0ASuppul", sV0ASuppul, medV0ASup);
 678     v0ASuppul->SetLineColor(kV0AColorOct);
 679     v0LE->AddNode(v0ASuppul,1);
 680
 681
 682     //Definition of sector 5
 683     TGeoVolume *v0ASec5 = new TGeoVolumeAssembly("V0ASec5");
 684
 685     /// For boolean sustraction
 686     double preShape5 = 0.2;
 687     for (int i=0;i<2;i++) {
 688       v0APts[0+8*i] = -fV0AR0+fV0AFraWd/2.-preShape5;  v0APts[1+8*i] = -preShape5;
 689       v0APts[2+8*i] = -fV0AR0+fV0AFraWd/2.-preShape5;  v0APts[3+8*i] = fV0AFraWd/2.;
 690       v0APts[4+8*i] = -fV0AR4-fV0AFraWd/2.+preShape5;  v0APts[5+8*i] = fV0AFraWd/2.;
 691       v0APts[6+8*i] = -fV0AR4-fV0AFraWd/2.+preShape5;  v0APts[7+8*i] = -preShape5;
 692     }
 693     new TGeoArb8("sV0ACha15",fV0ASciWd/1.5,v0APts);
 694     for (int i=0;i<2;i++) {
 695       v0APts[0+8*i] = -fV0AR0*cos45+preShape;
 696       v0APts[1+8*i] = -(fV0AR0-fV0AFraWd)*sin45+preShape;
 697       v0APts[2+8*i] = -(fV0AR0-fV0AFraWd/2.)*cos45+preShape;
 698       v0APts[3+8*i] = -(fV0AR0-fV0AFraWd/2.)*sin45;
 699       v0APts[4+8*i] = -(fV0AR4+fV0AFraWd/2.)*cos45-preShape;
 700       v0APts[5+8*i] = -(fV0AR4+fV0AFraWd/2.)*sin45-2.*preShape;
 701       v0APts[6+8*i] = -(fV0AR4+fV0AFraWd)*cos45-preShape;
 702       v0APts[7+8*i] = -fV0AR4*sin45-preShape;
 703     }
 704     new TGeoArb8("sV0ACha25", fV0ASciWd/2.+2.*preShape5, v0APts);
 705     new TGeoCompositeShape("sV0ACha125","sV0ACha15+sV0ACha25");
 706     new TGeoTube("sV0ANail15Hole", 0.0, 0.4, 1.65);
 707     TGeoTranslation *pos15 = new TGeoTranslation("pos15", -42.9, -0.51, 0.0);
 708     pos15->RegisterYourself();
 709     new TGeoTube("sV0ANail25Hole", 0.0, 0.4, 1.65);
 710     TGeoTranslation *pos25 = new TGeoTranslation("pos25",-30.8,-30.04,0.0);
 711     pos25->RegisterYourself();
 712     TGeoTranslation *pos35 = new TGeoTranslation("pos35",-30.05,-30.79,0.0);
 713     pos35->RegisterYourself();
 714     new TGeoCompositeShape("sV0ANailsHoles5","sV0ANail15Hole:pos15+sV0ANail25Hole:pos25");
 715     new TGeoCompositeShape("sV0ACha5","sV0ACha125+sV0ANailsHoles5");
 716     new TGeoTubeSeg("sV0AFicR55", fV0AR4+fV0AFraWd/2., fV0AR4 + fV0AR0, fV0ASciWd/2. +2*preShape, 180.0, 225.0);
 717     new TGeoTube("sV0ANail1PlaInHole5", 0.0, 0.4, (fV0APlaWd-2*fV0APlaAl)/2.);
 718     new TGeoTube("sV0ANail2PlaInHole5", 0.0, 0.4, (fV0APlaWd-2*fV0APlaAl)/2.);
 719     new TGeoTube("sV0ANail3PlaInHole5", 0.0, 0.4, (fV0APlaWd-2*fV0APlaAl)/2.);
 720     new TGeoCompositeShape("sV0ANailsPlaInHoles5","sV0ANail1PlaInHole5:pos15+sV0ANail2PlaInHole5:pos25+sV0ANail3PlaInHole5:pos35");
 721     new TGeoTube("sV0ANail1PlaOuHole5", 0.0, 0.4, (fV0APlaAl)/2.);
 722     new TGeoTube("sV0ANail2PlaOuHole5", 0.0, 0.4, (fV0APlaAl)/2.);
 723     new TGeoTube("sV0ANail3PlaOuHole5", 0.0, 0.4, (fV0APlaAl)/2.);
 724     new TGeoCompositeShape("sV0ANailsPlaOuHoles5","sV0ANail1PlaOuHole5:pos15+sV0ANail2PlaOuHole5:pos25+sV0ANail3PlaOuHole5:pos35");
 725
 726     /// Frame
 727     TGeoVolume *v0AFra5 = new TGeoVolumeAssembly("V0AFra5");
 728     for (int i=0;i<2;i++) {
 729       v0APts[0+8*i] = -fV0AR0+fV0AFraWd/2.;  v0APts[1+8*i] = 0.0;
 730       v0APts[2+8*i] = -fV0AR0+fV0AFraWd/2.;  v0APts[3+8*i] = fV0AFraWd/2.;
 731       v0APts[4+8*i] = -fV0AR4-fV0AFraWd/2.;  v0APts[5+8*i] = fV0AFraWd/2.;
 732       v0APts[6+8*i] = -fV0AR4-fV0AFraWd/2.;  v0APts[7+8*i] = 0.0;
 733     }
 734     TGeoArb8 *sV0AFraB15 = new TGeoArb8("sV0AFraB15",fV0ASciWd/2.,v0APts);
 735     TGeoVolume *v0AFraB15 = new TGeoVolume("V0AFraB15",sV0AFraB15,medV0AFra);
 736     for (int i=0;i<2;i++) {
 737       v0APts[0+8*i] = -fV0AR0*cos45;
 738       v0APts[1+8*i] = -(fV0AR0-fV0AFraWd)*sin45;
 739       v0APts[2+8*i] = -(fV0AR0-fV0AFraWd/2.)*cos45;
 740       v0APts[3+8*i] = -(fV0AR0-fV0AFraWd/2.)*sin45;
 741       v0APts[4+8*i] = -(fV0AR4+fV0AFraWd/2.)*cos45;
 742       v0APts[5+8*i] = -(fV0AR4+fV0AFraWd/2.)*sin45;
 743       v0APts[6+8*i] = -(fV0AR4+fV0AFraWd)*cos45;
 744       v0APts[7+8*i] = -fV0AR4*sin45;
 745     }
 746     TGeoArb8 *sV0AFraB25 = new TGeoArb8("sV0AFraB25", fV0ASciWd/2., v0APts);
 747     TGeoVolume *v0AFraB25 = new TGeoVolume("V0AFraB25",sV0AFraB25,medV0AFra);
 748     v0AFraB15->SetLineColor(kV0AColorFra); v0AFraB25->SetLineColor(kV0AColorFra);
 749     v0AFra5->AddNode(v0AFraB15,1);
 750     v0AFra5->AddNode(v0AFraB25,1);  // Prefer 2 GeoObjects insted of 3 GeoMovements
 751     new TGeoTubeSeg( "sV0AFraR1b5", fV0AR0-fV0AFraWd/2.,
 752                      fV0AR0+fV0AFraWd/2., fV0ASciWd/2., 180.0, 225.0);
 753     new TGeoTubeSeg( "sV0AFraR2b5", fV0AR1-fV0AFraWd/2.,
 754                      fV0AR1+fV0AFraWd/2., fV0ASciWd/2., 180.0, 225.0);
 755     new TGeoTubeSeg( "sV0AFraR3b5", fV0AR2-fV0AFraWd/2.,
 756                      fV0AR2+fV0AFraWd/2., fV0ASciWd/2., 180.0, 225.0);
 757     new TGeoTubeSeg( "sV0AFraR4b5", fV0AR3-fV0AFraWd/2.,
 758                      fV0AR3+fV0AFraWd/2., fV0ASciWd/2., 180.0, 225.0);
 759     new TGeoTubeSeg( "sV0AFraR5b5", fV0AR4-fV0AFraWd/2.,
 760                      fV0AR4+fV0AFraWd/2., fV0ASciWd/2., 180.0, 225.0);
 761     TGeoCompositeShape *sV0AFraR15 = new TGeoCompositeShape("sV0AFraR15","sV0AFraR1b5-sV0ACha5");
 762     TGeoCompositeShape *sV0AFraR25 = new TGeoCompositeShape("sV0AFraR25","sV0AFraR2b5-sV0ACha5");
 763     TGeoCompositeShape *sV0AFraR35 = new TGeoCompositeShape("sV0AFraR35","sV0AFraR3b5-sV0ACha5");
 764     TGeoCompositeShape *sV0AFraR45 = new TGeoCompositeShape("sV0AFraR45","sV0AFraR4b5-sV0ACha5");
 765     TGeoCompositeShape *sV0AFraR55 = new TGeoCompositeShape("sV0AFraR55","sV0AFraR5b5-sV0ACha5");
 766     TGeoVolume *v0AFraR15 = new TGeoVolume("V0AFraR15",sV0AFraR15,medV0AFra);
 767     TGeoVolume *v0AFraR25 = new TGeoVolume("V0AFraR25",sV0AFraR25,medV0AFra);
 768     TGeoVolume *v0AFraR35 = new TGeoVolume("V0AFraR35",sV0AFraR35,medV0AFra);
 769     TGeoVolume *v0AFraR45 = new TGeoVolume("V0AFraR45",sV0AFraR45,medV0AFra);
 770     TGeoVolume *v0AFraR55 = new TGeoVolume("V0AFraR55",sV0AFraR55,medV0AFra);
 771     v0AFraR15->SetLineColor(kV0AColorFra); v0AFraR25->SetLineColor(kV0AColorFra);
 772     v0AFraR35->SetLineColor(kV0AColorFra); v0AFraR45->SetLineColor(kV0AColorFra);
 773     v0AFraR55->SetLineColor(kV0AColorFra);
 774     v0AFra5->AddNode(v0AFraR15,1);
 775     v0AFra5->AddNode(v0AFraR25,1);
 776     v0AFra5->AddNode(v0AFraR35,1);
 777     v0AFra5->AddNode(v0AFraR45,1);
 778     v0AFra5->AddNode(v0AFraR55,1);
 779     v0ASec5->AddNode(v0AFra5,1);
 780
 781     /// Sensitive scintilator
 782     TGeoVolume *v0ASci5 = new TGeoVolumeAssembly("V0ASci5");
 783     new TGeoTubeSeg( "sV0AR1b5", fV0AR0+fV0AFraWd/2.,
 784                      fV0AR1-fV0AFraWd/2., fV0ASciWd/2., 180.0, 225.0);
 785     new TGeoTubeSeg( "sV0AR2b5", fV0AR1+fV0AFraWd/2.,
 786                      fV0AR2-fV0AFraWd/2., fV0ASciWd/2., 180.0, 225.0);
 787     new TGeoTubeSeg( "sV0AR3b5", fV0AR2+fV0AFraWd/2.,
 788                      fV0AR3-fV0AFraWd/2., fV0ASciWd/2., 180.0, 225.0);
 789     new TGeoTubeSeg( "sV0AR4b5", fV0AR3+fV0AFraWd/2.,
 790                      fV0AR4-fV0AFraWd/2., fV0ASciWd/2., 180.0, 225.0);
 791     TGeoCompositeShape *sV0AR15 = new TGeoCompositeShape("sV0AR15","sV0AR1b5-sV0ACha5");
 792     TGeoCompositeShape *sV0AR25 = new TGeoCompositeShape("sV0AR25","sV0AR2b5-sV0ACha5");
 793     TGeoCompositeShape *sV0AR35 = new TGeoCompositeShape("sV0AR35","sV0AR3b5-sV0ACha5");
 794     TGeoCompositeShape *sV0AR45 = new TGeoCompositeShape("sV0AR45","sV0AR4b5-sV0ACha5");
 795     TGeoVolume *v0L15 = new TGeoVolume("V0L15",sV0AR15,medV0ASci);
 796     TGeoVolume *v0L25 = new TGeoVolume("V0L25",sV0AR25,medV0ASci);
 797     TGeoVolume *v0L35 = new TGeoVolume("V0L35",sV0AR35,medV0ASci);
 798     TGeoVolume *v0L45 = new TGeoVolume("V0L45",sV0AR45,medV0ASci);
 799     v0L15->SetLineColor(kV0AColorSci); v0L25->SetLineColor(kV0AColorSci);
 800     v0L35->SetLineColor(kV0AColorSci); v0L45->SetLineColor(kV0AColorSci);
 801     v0ASci5->AddNode(v0L15,1);
 802     v0ASci5->AddNode(v0L25,1);
 803     v0ASci5->AddNode(v0L35,1);
 804     v0ASci5->AddNode(v0L45,1);
 805
 806      /// Segment of octagon
 807     for (int i=0;i<2;i++) {
 808     v0APts[0+8*i] = -fV0AR6+fV0AOctH2;           v0APts[1+8*i] = 0.;
 809     v0APts[2+8*i] = -(fV0AR7-fV0AOctH2)*cos654;   v0APts[3+8*i] = -(fV0AR7-fV0AOctH2)*sin654;
 810     v0APts[4+8*i] = -fV0AR7*cos654;              v0APts[5+8*i] = -fV0AR7*sin654;
 811     v0APts[6+8*i] = -fV0AR6;                     v0APts[7+8*i] = 0.;
 812     }
 813     TGeoArb8 *sV0AOct25 = new TGeoArb8("sV0AOct25", (fV0ASciWd+2*fV0AOctWd)/2., v0APts);
 814     TGeoVolume *v0AOct25 = new TGeoVolume("V0AOct25", sV0AOct25,medV0ASup);
 815     v0AOct25->SetLineColor(kV0AColorOct);
 816     TGeoVolume *v0ASup5 = new TGeoVolumeAssembly("V0ASup5");
 817     v0ASup5->AddNode(v0AOct25,1);
 818     v0ASec5->AddNode(v0ASup5,1);
 819
 820     //Bunch of fibers
 821     v0APts[ 0] = v0APts[ 2] = -14.0;
 822     v0APts[ 1] = v0APts[ 7] = (fV0ASciWd+fV0AOctWd)/2.-0.01;
 823     v0APts[ 3] = v0APts[ 5] = (fV0ASciWd+fV0AOctWd)/2.+0.01;
 824     v0APts[ 4] = v0APts[ 6] = +14.0;
 825     v0APts[ 8] = v0APts[10] = -10.0;
 826     v0APts[ 9] = v0APts[15] = 0.;
 827     v0APts[11] = v0APts[13] = 0.25;
 828     v0APts[12] = v0APts[14] = +10.0;
 829     new TGeoArb8("sV0AFib15", 11.8, v0APts);
 830     rot = new TGeoRotation("rot");
 831     rot->RotateX(-90);
 832     rot->RotateZ(90+22.5);
 833     TGeoCombiTrans *fib15pos = new TGeoCombiTrans("fib15pos", -(fV0AR6-fV0AOctH2+fV0AR5)*cos225/2. +
 834     3.3, -(fV0AR6-fV0AOctH2+fV0AR5)*sin225/2. + 1.5, 0, rot);
 835     fib15pos->RegisterYourself();
 836     TGeoCompositeShape *sV0AFib15Hole = new TGeoCompositeShape("sV0AFib15Hole", "sV0AFib15:fib15pos-sV0AFicR55");
 837     TGeoVolume *v0AFib15Hole = new TGeoVolume("V0AFib15",sV0AFib15Hole,medV0AFib);
 838     v0AFib15Hole->SetLineColor(kV0AColorFib);
 839     new TGeoArb8("sV0AFib25", 11.8, v0APts);
 840     rot = new TGeoRotation("rot");
 841     rot->RotateX(-90);
 842     rot->RotateY(180);
 843     rot->RotateZ(90+22.5);
 844     TGeoCombiTrans *fib25pos = new TGeoCombiTrans("fib25pos", -(fV0AR6-fV0AOctH2+fV0AR5)*cos225/2. +
 845     3.3, -(fV0AR6-fV0AOctH2+fV0AR5)*sin225/2. + 1.5, 0, rot);
 846     fib25pos->RegisterYourself();
 847     TGeoCompositeShape *sV0AFib25Hole = new TGeoCompositeShape("sV0AFib25Hole","sV0AFib25:fib25pos-sV0AFicR55");
 848     TGeoVolume *v0AFib25Hole = new TGeoVolume("V0AFib25Hole",sV0AFib25Hole,medV0AFib);
 849     v0AFib25Hole->SetLineColor(kV0AColorFib);
 850     TGeoVolume *v0AFib5 = new TGeoVolumeAssembly("V0AFib5");
 851     v0AFib5->AddNode(v0AFib15Hole,1);
 852     v0AFib5->AddNode(v0AFib25Hole,1);
 853     v0ASec5->AddNode(v0AFib5,1);
 854
 855     /// Non-sensitive scintilator
 856     new TGeoTubeSeg("sV0AR5S25", fV0AR4+fV0AFraWd/2., fV0AR4 + fV0AR0, fV0ASciWd/2.+2*preShape, 180.0, 225.0);
 857     TGeoCompositeShape *sV0AR55 = new TGeoCompositeShape("V0AR55","sV0AR5S25 - sV0ACha5");
 858     TGeoVolume *v0AR55 = new TGeoVolume("V0AR55",sV0AR55,medV0ASci);
 859     v0AR55->SetLineColor(kV0AColorSci);
 860     v0ASci5->AddNode(v0AR55,1);
 861     v0ASec5->AddNode(v0ASci5,1);
 862
 863     /// Plates
 864     for (int i=0;i<2;i++) {
 865       v0APts[0+8*i] = -fV0AR0;                  v0APts[1+8*i] = 0.;
 866       v0APts[2+8*i] = -fV0AR0*cos654;           v0APts[3+8*i] = -fV0AR0*sin654;
 867       v0APts[4+8*i] = -fV0AR7*cos654;           v0APts[5+8*i] = -fV0AR7*sin654;
 868       v0APts[6+8*i] = -fV0AR6;                  v0APts[7+8*i] = 0.;
 869     }
 870     new TGeoArb8("sV0APlaIn5", (fV0APlaWd-2*fV0APlaAl)/2., v0APts);
 871     TGeoCompositeShape *sV0APlaInNailsHoles5 = new TGeoCompositeShape("sV0APlaInNailsHoles5","sV0APlaIn5-sV0ANailsPlaInHoles5");
 872     TGeoVolume *v0APlaInNailsHoles5 = new TGeoVolume("V0APlaInNailsHoles5", sV0APlaInNailsHoles5, medV0APlaIn);
 873     new TGeoArb8("sV0APlaOu5", fV0APlaAl/2., v0APts);
 874     TGeoCompositeShape *sV0APlaOuNailsHoles5 = new TGeoCompositeShape("sV0APlaOuNailsHoles5","sV0APlaOu5-sV0ANailsPlaOuHoles5");
 875     TGeoVolume *v0APlaOuNailsHoles5 = new TGeoVolume("V0APlaOuNailsHoles5", sV0APlaOuNailsHoles5, medV0APlaOu);
 876     v0APlaInNailsHoles5->SetLineColor(kV0AColorPlaIn); v0APlaOuNailsHoles5->SetLineColor(kV0AColorPlaOu);
 877     TGeoVolume *v0APla5 = new TGeoVolumeAssembly("V0APla5");
 878     v0APla5->AddNode(v0APlaInNailsHoles5,1);
 879     v0APla5->AddNode(v0APlaOuNailsHoles5,1,new TGeoTranslation(0,0,(fV0APlaWd-fV0APlaAl)/2.));
 880     v0APla5->AddNode(v0APlaOuNailsHoles5,2,new TGeoTranslation(0,0,-(fV0APlaWd-fV0APlaAl)/2.));
 881     v0ASec5->AddNode(v0APla5,1,new TGeoTranslation(0,0,(fV0ASciWd+2*fV0AOctWd+fV0APlaWd)/2.));
 882     v0ASec5->AddNode(v0APla5,2,new TGeoTranslation(0,0,-(fV0ASciWd+2*fV0AOctWd+fV0APlaWd)/2.));
 883
 884     /// PMBox
 885     TGeoVolume* v0APM5 = new TGeoVolumeAssembly("V0APM5");
 886     new TGeoBBox("sV0APMB15", fV0APMBWd/2., fV0APMBHt/2., fV0APMBTh/2.);
 887     new TGeoBBox("sV0APMB25", fV0APMBWd/2.-fV0APMBWdW, fV0APMBHt/2.-fV0APMBHtW, fV0APMBTh/2.-fV0APMBThW);
 888     TGeoCompositeShape *sV0APMB5 = new TGeoCompositeShape("sV0APMB5","sV0APMB15-sV0APMB25");
 889     TGeoVolume *v0APMB5 = new TGeoVolume("V0APMB5",sV0APMB5, medV0APMAlum);
 890     v0APMB5->SetLineColor(kV0AColorPMA);
 891     v0APM5->AddNode(v0APMB5,1);
 892
 893     /// PMTubes
 894     TGeoTube *sV0APMT15 = new TGeoTube("sV0APMT15", fV0APMTR1, fV0APMTR2, fV0APMTH/2.);
 895     TGeoVolume *v0APMT15 = new TGeoVolume("V0APMT15", sV0APMT15, medV0APMGlass);
 896     TGeoTube *sV0APMT25 = new TGeoTube("sV0APMT25", fV0APMTR3, fV0APMTR4, fV0APMTH/2.);
 897     TGeoVolume *v0APMT25 = new TGeoVolume("V0APMT25", sV0APMT25, medV0APMAlum);
 898     TGeoVolume *v0APMT5 = new TGeoVolumeAssembly("V0APMT5");
 899     TGeoTube *sV0APMTT5 = new TGeoTube("sV0APMTT5", 0., fV0APMTR4, fV0APMTB/2.);
 900     TGeoVolume *v0APMTT5 = new TGeoVolume("V0APMTT5", sV0APMTT5, medV0APMAlum);
 901     v0APMT5->SetLineColor(kV0AColorPMG);
 902     v0APMT25->SetLineColor(kV0AColorPMA);
 903     v0APMTT5->SetLineColor(kV0AColorPMA);
 904     rot = new TGeoRotation("rot", 90, 0, 180, 0, 90, 90);
 905     v0APMT5->AddNode(v0APMT15,1,rot);
 906     v0APMT5->AddNode(v0APMT25,1,rot);
 907     v0APMT5->AddNode(v0APMTT5,1,new TGeoCombiTrans(0,(fV0APMTH+fV0APMTB)/2.,0,rot));
 908     double autoShift5 = (fV0APMBWd-2*fV0APMBWdW)/4.;
 909     v0APM5->AddNode(v0APMT5, 1, new TGeoTranslation(-1.5*autoShift5, 0, 0));
 910     v0APM5->AddNode(v0APMT5, 2, new TGeoTranslation(-0.5*autoShift5, 0, 0));
 911     v0APM5->AddNode(v0APMT5, 3, new TGeoTranslation(+0.5*autoShift5, 0, 0));
 912     v0APM5->AddNode(v0APMT5, 4, new TGeoTranslation(+1.5*autoShift5, 0, 0));
 913
 914     /// PM
 915     rot = new TGeoRotation("rot");
 916     rot->RotateX(-90+30);
 917     rot->RotateY(0);
 918     rot->RotateZ(-65-3);
 919     double cosAngPMB5 = TMath::Cos(fV0APMBAng*TMath::DegToRad());
 920     double sinAngPMB5 = TMath::Sin(fV0APMBAng*TMath::DegToRad());
 921     double shiftZ5 = fV0APMBHt/2. * cosAngPMB5
 922       -   ( fV0ASciWd + 2 * fV0AOctWd + 2 * fV0APlaWd )/2.   -   fV0APMBTh/2. * sinAngPMB5;
 923     double shiftR5 = fV0AR6  +  fV0AOctH2 + fV0APlaAl;
 924     v0ASec5->AddNode(v0APM5,1, new TGeoCombiTrans( -shiftR5*cos225-1.3, -shiftR5*sin225, shiftZ5, rot));
 925
 926     // Aluminium nails
 927     TGeoTube *sV0ANail51 = new TGeoTube("sV0ANail51", 0.0, 0.4, 5.09/2.);
 928     TGeoVolume *v0ANail51 = new TGeoVolume("V0ANail51", sV0ANail51, medV0APMAlum);
 929     v0ANail51->SetLineColor(kV0AColorPMA);// this is the color for aluminium
 930     v0ASec5->AddNode(v0ANail51,1,new TGeoTranslation(-42.9,-0.51,0.0));
 931     TGeoTube *sV0ANail52 = new TGeoTube("sV0ANail52", 0.0, 0.4, 5.09/2.);
 932     TGeoVolume *v0ANail52 = new TGeoVolume("V0ANail52", sV0ANail52, medV0APMAlum);
 933     v0ANail52->SetLineColor(kV0AColorPMA);
 934     v0ASec5->AddNode(v0ANail52,1,new TGeoTranslation(-30.8,-30.04,0.0));
 935
 936     // Adding sector to v0LE volume
 937     v0LE->AddNode(v0ASec5, 1);
 938
 939
 940     //Definition of  sector 6
 941     TGeoVolume *v0ASec6 = new TGeoVolumeAssembly("V0ASec6");
 942
 943     /// For boolean sustraction
 944     double preShape6 = 0.2;
 945     for (int i=0;i<2;i++) {
 946     v0APts[0+8*i] = -preShape6;                      v0APts[1+8*i] = -fV0AR0+fV0AFraWd/2.-preShape6;
 947     v0APts[2+8*i] = -fV0AFraWd/2.;                    v0APts[3+8*i] = -fV0AR0+fV0AFraWd/2.-preShape6;
 948     v0APts[4+8*i] = -fV0AFraWd/2.;                    v0APts[5+8*i] = -fV0AR4-fV0AFraWd/2.+preShape6;
 949     v0APts[6+8*i] = -preShape6;                      v0APts[7+8*i] = -fV0AR4-fV0AFraWd/2.+preShape6;
 950     }
 951     new TGeoArb8("sV0ACha16",fV0ASciWd/1.5,v0APts);
 952     for (int i=0;i<2;i++) {
 953     v0APts[0+8*i] = -fV0AR0*cos45+preShape;
 954     v0APts[1+8*i] = -(fV0AR0-fV0AFraWd)*sin45+preShape;
 955     v0APts[2+8*i] = -(fV0AR0-fV0AFraWd/2.)*cos45+preShape;
 956     v0APts[3+8*i] = -(fV0AR0-fV0AFraWd/2.)*sin45;
 957     v0APts[4+8*i] = -(fV0AR4+fV0AFraWd/2.)*cos45-preShape;
 958     v0APts[5+8*i] = -(fV0AR4+fV0AFraWd/2.)*sin45-preShape;
 959     v0APts[6+8*i] = -(fV0AR4+fV0AFraWd)*cos45-preShape;
 960     v0APts[7+8*i] = -fV0AR4*sin45-preShape;
 961     }
 962     new TGeoArb8("sV0ACha26", fV0ASciWd/2.+2.*preShape, v0APts);
 963     new TGeoCompositeShape("sV0ACha126","sV0ACha16+sV0ACha26");
 964     new TGeoTube("sV0ANail16Hole", 0.0, 0.4, 1.65);
 965     TGeoTranslation *pos16 = new TGeoTranslation("pos16",-0.51,-42.9,0.0);
 966     pos16->RegisterYourself();
 967     new TGeoTube("sV0ANail26Hole", 0.0, 0.4, 1.65);
 968     TGeoTranslation *pos26 = new TGeoTranslation("pos26",-30.05,-30.79,0.0);
 969     pos26->RegisterYourself();
 970     new TGeoCompositeShape("sV0ANailsHoles6","sV0ANail16Hole:pos16+sV0ANail26Hole:pos26");
 971     new TGeoCompositeShape("sV0ACha6","sV0ACha126+sV0ANailsHoles6");
 972     new TGeoTubeSeg("sV0AFicR56", fV0AR4+fV0AFraWd/2., fV0AR4 + fV0AR0, fV0ASciWd/2.+2*preShape, 225, 270.0);
 973     new TGeoTube("sV0ANail1PlaInHole6", 0.0, 0.4, (fV0APlaWd-2*fV0APlaAl)/2.);
 974     new TGeoTube("sV0ANail1PlaOuHole6", 0.0, 0.4, (fV0APlaAl)/2.);
 975
 976     /// Frame
 977     TGeoVolume *v0AFra6 = new TGeoVolumeAssembly("V0AFra6");
 978     for (int i=0;i<2;i++) {
 979     v0APts[0+8*i] = 0.;                              v0APts[1+8*i] = -fV0AR0+fV0AFraWd/2.;
 980     v0APts[2+8*i] = -fV0AFraWd/2.;                    v0APts[3+8*i] = -fV0AR0+fV0AFraWd/2.;
 981     v0APts[4+8*i] = -fV0AFraWd/2.;                    v0APts[5+8*i] = -fV0AR4-fV0AFraWd/2.;
 982     v0APts[6+8*i] = 0.;                              v0APts[7+8*i] = -fV0AR4-fV0AFraWd/2.;
 983     }
 984     TGeoArb8 *sV0AFraB16 = new TGeoArb8("sV0AFraB16",fV0ASciWd/2.,v0APts);
 985     TGeoVolume *v0AFraB16 = new TGeoVolume("V0AFraB16",sV0AFraB16,medV0AFra);
 986     for (int i=0;i<2;i++) {
 987     v0APts[0+8*i] = -fV0AR0*cos45;
 988     v0APts[1+8*i] = -(fV0AR0+fV0AFraWd)*sin45;
 989     v0APts[2+8*i] = -(fV0AR0+fV0AFraWd/2.)*cos45;
 990     v0APts[3+8*i] = -(fV0AR0+fV0AFraWd/2.)*sin45;
 991     v0APts[4+8*i] = -(fV0AR4-fV0AFraWd/2.)*cos45;
 992     v0APts[5+8*i] = -(fV0AR4-fV0AFraWd/2.)*sin45;
 993     v0APts[6+8*i] = -(fV0AR4-fV0AFraWd)*cos45;
 994     v0APts[7+8*i] = -fV0AR4*sin45;
 995     }
 996     TGeoArb8 *sV0AFraB26 = new TGeoArb8("sV0AFraB26", fV0ASciWd/2., v0APts);
 997     TGeoVolume *v0AFraB26 = new TGeoVolume("V0AFraB26",sV0AFraB26,medV0AFra);
 998     v0AFraB16->SetLineColor(kV0AColorFra); v0AFraB26->SetLineColor(kV0AColorFra);
 999     v0AFra6->AddNode(v0AFraB16,1);
1000     v0AFra6->AddNode(v0AFraB26,1);  // Prefer 2 GeoObjects insted of 3 GeoMovements
1001     new TGeoTubeSeg( "sV0AFraR1b6", fV0AR0-fV0AFraWd/2.,
1002              fV0AR0+fV0AFraWd/2., fV0ASciWd/2., 225.0, 270.0);
1003     new TGeoTubeSeg( "sV0AFraR2b6", fV0AR1-fV0AFraWd/2.,
1004              fV0AR1+fV0AFraWd/2., fV0ASciWd/2., 225.0, 270.0);
1005     new TGeoTubeSeg( "sV0AFraR3b6", fV0AR2-fV0AFraWd/2.,
1006              fV0AR2+fV0AFraWd/2., fV0ASciWd/2., 225.0, 270.0);
1007     new TGeoTubeSeg( "sV0AFraR4b6", fV0AR3-fV0AFraWd/2.,
1008              fV0AR3+fV0AFraWd/2., fV0ASciWd/2., 225.0, 270.0);
1009     new TGeoTubeSeg( "sV0AFraR5b6", fV0AR4-fV0AFraWd/2.,
1010              fV0AR4+fV0AFraWd/2., fV0ASciWd/2., 225.0, 270.0);
1011     TGeoCompositeShape *sV0AFraR16 = new TGeoCompositeShape("sV0AFraR16","sV0AFraR1b6-sV0ACha6");
1012     TGeoCompositeShape *sV0AFraR26 = new TGeoCompositeShape("sV0AFraR26","sV0AFraR2b6-sV0ACha6");
1013     TGeoCompositeShape *sV0AFraR36 = new TGeoCompositeShape("sV0AFraR36","sV0AFraR3b6-sV0ACha6");
1014     TGeoCompositeShape *sV0AFraR46 = new TGeoCompositeShape("sV0AFraR46","sV0AFraR4b6-sV0ACha6");
1015     TGeoCompositeShape *sV0AFraR56 = new TGeoCompositeShape("sV0AFraR56","sV0AFraR5b6-sV0ACha6");
1016     TGeoVolume *v0AFraR16 = new TGeoVolume("V0AFraR16",sV0AFraR16,medV0AFra);
1017     TGeoVolume *v0AFraR26 = new TGeoVolume("V0AFraR26",sV0AFraR26,medV0AFra);
1018     TGeoVolume *v0AFraR36 = new TGeoVolume("V0AFraR36",sV0AFraR36,medV0AFra);
1019     TGeoVolume *v0AFraR46 = new TGeoVolume("V0AFraR46",sV0AFraR46,medV0AFra);
1020     TGeoVolume *v0AFraR56 = new TGeoVolume("V0AFraR56",sV0AFraR56,medV0AFra);
1021     v0AFraR16->SetLineColor(kV0AColorFra); v0AFraR26->SetLineColor(kV0AColorFra);
1022     v0AFraR36->SetLineColor(kV0AColorFra); v0AFraR46->SetLineColor(kV0AColorFra);
1023     v0AFraR56->SetLineColor(kV0AColorFra);
1024     v0AFra6->AddNode(v0AFraR16,1);
1025     v0AFra6->AddNode(v0AFraR26,1);
1026     v0AFra6->AddNode(v0AFraR36,1);
1027     v0AFra6->AddNode(v0AFraR46,1);
1028     v0AFra6->AddNode(v0AFraR56,1);
1029     v0ASec6->AddNode(v0AFra6,1);
1030
1031     /// Sensitive scintilator
1032     TGeoVolume *v0ASci6 = new TGeoVolumeAssembly("V0ASci6");
1033     new TGeoTubeSeg( "sV0AR1b6", fV0AR0+fV0AFraWd/2.,
1034                      fV0AR1-fV0AFraWd/2., fV0ASciWd/2., 225.0, 270.0);
1035     new TGeoTubeSeg( "sV0AR2b6", fV0AR1+fV0AFraWd/2.,
1036              fV0AR2-fV0AFraWd/2., fV0ASciWd/2., 225.0, 270.0);
1037     new TGeoTubeSeg( "sV0AR3b6", fV0AR2+fV0AFraWd/2.,
1038              fV0AR3-fV0AFraWd/2., fV0ASciWd/2., 225.0, 270.0);
1039     new TGeoTubeSeg( "sV0AR4b6", fV0AR3+fV0AFraWd/2.,
1040              fV0AR4-fV0AFraWd/2., fV0ASciWd/2., 225.0, 270.0);
1041     TGeoCompositeShape *sV0AR16 = new TGeoCompositeShape("sV0AR16","sV0AR1b6-sV0ACha6");
1042     TGeoCompositeShape *sV0AR26 = new TGeoCompositeShape("sV0AR26","sV0AR2b6-sV0ACha6");
1043     TGeoCompositeShape *sV0AR36 = new TGeoCompositeShape("sV0AR36","sV0AR3b6-sV0ACha6");
1044     TGeoCompositeShape *sV0AR46 = new TGeoCompositeShape("sV0AR46","sV0AR4b6-sV0ACha6");
1045     TGeoVolume *v0L16 = new TGeoVolume("V0L16",sV0AR16,medV0ASci);
1046     TGeoVolume *v0L26 = new TGeoVolume("V0L26",sV0AR26,medV0ASci);
1047     TGeoVolume *v0L36 = new TGeoVolume("V0L36",sV0AR36,medV0ASci);
1048     TGeoVolume *v0L46 = new TGeoVolume("V0L46",sV0AR46,medV0ASci);
1049     v0L16->SetLineColor(kV0AColorSci); v0L26->SetLineColor(kV0AColorSci);
1050     v0L36->SetLineColor(kV0AColorSci); v0L46->SetLineColor(kV0AColorSci);
1051     v0ASci6->AddNode(v0L16,1);
1052     v0ASci6->AddNode(v0L26,1);
1053     v0ASci6->AddNode(v0L36,1);
1054     v0ASci6->AddNode(v0L46,1);
1055
1056     // Bunch of fibers
1057     for (int i=0;i<2;i++) {
1058     v0APts[0+8*i] = -10.0;          v0APts[1+8*i] = 13.1;
1059     v0APts[2+8*i] = 10.0;           v0APts[3+8*i] = 13.1;
1060     v0APts[4+8*i] = 8.0;            v0APts[5+8*i] = -29.0;
1061     v0APts[6+8*i] = -12.0;          v0APts[7+8*i] = -12.0;
1062     }
1063     new TGeoArb8("sV0AFib16", 0.01, v0APts);
1064     rot = new TGeoRotation("rot");
1065     rot->RotateX(2.0);
1066     rot->RotateY(180.0);
1067     rot->RotateZ(90+22.5);
1068     TGeoCombiTrans *fib16pos = new TGeoCombiTrans("fib16pos", -40.0 + 3.3, -50.0 + 1.5, 0.5, rot);
1069     fib16pos->RegisterYourself();
1070     TGeoCompositeShape *sV0AFib16Hole = new TGeoCompositeShape("sV0AFib16Hole", "sV0AFib16:fib16pos-sV0AFicR56");
1071     TGeoVolume *v0AFib16Hole = new TGeoVolume("V0AFib16Hole",sV0AFib16Hole,medV0AFib);
1072     v0AFib16Hole->SetLineColor(kV0AColorFib);
1073     new TGeoArb8("sV0AFib26", 0.01, v0APts);
1074     rot = new TGeoRotation("rot");
1075     rot->RotateX(-2.0);
1076     rot->RotateY(180.0);
1077     rot->RotateZ(90+22.5);
1078     TGeoCombiTrans *fib26pos = new TGeoCombiTrans("fib26pos", -40.0 + 3.3, -50.0 + 1.5, -0.5, rot);
1079     fib26pos->RegisterYourself();
1080     TGeoCompositeShape *sV0AFib26Hole = new TGeoCompositeShape("sV0AFib26Hole", "sV0AFib26:fib26pos-sV0AFicR56");
1081     TGeoVolume *v0AFib26Hole = new TGeoVolume("V0AFib26Hole",sV0AFib26Hole,medV0AFib);
1082     v0AFib26Hole->SetLineColor(kV0AColorFib);
1083     TGeoVolume *v0AFib6 = new TGeoVolumeAssembly("V0AFib6");
1084     v0AFib6->AddNode(v0AFib16Hole,1);
1085     v0AFib6->AddNode(v0AFib26Hole,1);
1086     v0ASec6->AddNode(v0AFib6,1);
1087
1088     /// Non-sensitive scintilator
1089     new TGeoTubeSeg("sV0AR5S26", fV0AR4+fV0AFraWd/2., fV0AR4 + fV0AR0, fV0ASciWd/2.+2*preShape, 225.0, 270.0);
1090     TGeoCompositeShape *sV0AR56 = new TGeoCompositeShape("V0AR56","sV0AR5S26 - sV0ACha6");
1091     TGeoVolume *v0AR56 = new TGeoVolume("V0AR56",sV0AR56,medV0ASci);
1092     v0AR56->SetLineColor(kV0AColorSci);
1093     v0ASci6->AddNode(v0AR56,1);
1094     v0ASec6->AddNode(v0ASci6,1);
1095
1096     /// Segment of octagon
1097     for (int i=0;i<2;i++) {
1098     v0APts[0+8*i] = 0.;                          v0APts[1+8*i] = -(fV0AR7-fV0AOctH2)*sin654;
1099     v0APts[2+8*i] = 0.;                          v0APts[3+8*i] = -fV0AR7*sin654;
1100     v0APts[4+8*i] = -fV0AR7*cos654;              v0APts[5+8*i] = -fV0AR7*sin654;
1101     v0APts[6+8*i] = -(fV0AR7-fV0AOctH2)*cos654;  v0APts[7+8*i] = -(fV0AR7-fV0AOctH2)*sin654;
1102     }
1103     TGeoArb8 *sV0AOct26 = new TGeoArb8("sV0AOct26", (fV0ASciWd+2*fV0AOctWd)/2., v0APts);
1104     TGeoVolume *v0AOct26 = new TGeoVolume("V0AOct26", sV0AOct26,medV0ASup);
1105     v0AOct26->SetLineColor(kV0AColorOct);
1106     TGeoVolume *v0ASup6 = new TGeoVolumeAssembly("V0ASup6");
1107     v0ASup6->AddNode(v0AOct26,1);
1108     v0ASec6->AddNode(v0ASup6,1);
1109
1110     /// Plates
1111     for (int i=0;i<2;i++) {
1112     v0APts[0+8*i] = 0.;                 v0APts[1+8*i] = -fV0AR0;
1113     v0APts[2+8*i] = 0.;                 v0APts[3+8*i] = -fV0AR7*sin654;
1114     v0APts[4+8*i] = -fV0AR7*cos654;     v0APts[5+8*i] = -fV0AR7*sin654;
1115     v0APts[6+8*i] = -fV0AR0*cos654;             v0APts[7+8*i] = -fV0AR0*sin654;
1116     }
1117     new TGeoArb8("sV0APlaIn6", (fV0APlaWd-2*fV0APlaAl)/2., v0APts);
1118     TGeoCompositeShape *sV0APlaInNailsHoles6 = new TGeoCompositeShape("sV0APlaInNailsHoles6","sV0APlaIn6-sV0ANail1PlaInHole6:pos16");
1119     TGeoVolume *v0APlaInNailsHoles6 = new TGeoVolume("V0APlaInNailsHoles6", sV0APlaInNailsHoles6, medV0APlaIn);
1120     new TGeoArb8("sV0APlaOu6", fV0APlaAl/2., v0APts);
1121     TGeoCompositeShape *sV0APlaOuNailsHoles6 = new TGeoCompositeShape("sV0APlaOuNailsHoles6","sV0APlaOu6-sV0ANail1PlaOuHole6:pos16");
1122     TGeoVolume *v0APlaOuNailsHoles6 = new TGeoVolume("V0APlaOuNailsHoles6", sV0APlaOuNailsHoles6, medV0APlaOu);
1123     v0APlaInNailsHoles6->SetLineColor(kV0AColorPlaIn); v0APlaOuNailsHoles6->SetLineColor(kV0AColorPlaOu);
1124     TGeoVolume *v0APla6 = new TGeoVolumeAssembly("V0APla6");
1125     v0APla6->AddNode(v0APlaInNailsHoles6,1);
1126     v0APla6->AddNode(v0APlaOuNailsHoles6,1,new TGeoTranslation(0,0,(fV0APlaWd-fV0APlaAl)/2.));
1127     v0APla6->AddNode(v0APlaOuNailsHoles6,2,new TGeoTranslation(0,0,-(fV0APlaWd-fV0APlaAl)/2.));
1128     v0ASec6->AddNode(v0APla6,1,new TGeoTranslation(0,0,(fV0ASciWd+2*fV0AOctWd+fV0APlaWd)/2.));
1129     v0ASec6->AddNode(v0APla6,2,new TGeoTranslation(0,0,-(fV0ASciWd+2*fV0AOctWd+fV0APlaWd)/2.));
1130
1131     /// PMBox
1132     TGeoVolume* v0APM6 = new TGeoVolumeAssembly("V0APM6");
1133     new TGeoBBox("sV0APMB16", fV0APMBWd/2., fV0APMBHt/2., fV0APMBTh/2.);
1134     new TGeoBBox("sV0APMB26", fV0APMBWd/2.-fV0APMBWdW, fV0APMBHt/2.-fV0APMBHtW, fV0APMBTh/2.-fV0APMBThW);
1135     TGeoCompositeShape *sV0APMB6 = new TGeoCompositeShape("sV0APMB6","sV0APMB16-sV0APMB26");
1136     TGeoVolume *v0APMB6 = new TGeoVolume("V0APMB6",sV0APMB6, medV0APMAlum);
1137     v0APMB6->SetLineColor(kV0AColorPMA);
1138     v0APM6->AddNode(v0APMB6,1);
1139
1140     /// PMTubes
1141     TGeoTube *sV0APMT16 = new TGeoTube("sV0APMT16", fV0APMTR1, fV0APMTR2, fV0APMTH/2.);
1142     TGeoVolume *v0APMT16 = new TGeoVolume("V0APMT16", sV0APMT16, medV0APMGlass);
1143     TGeoTube *sV0APMT26 = new TGeoTube("sV0APMT26", fV0APMTR3, fV0APMTR4, fV0APMTH/2.);
1144     TGeoVolume *v0APMT26 = new TGeoVolume("V0APMT26", sV0APMT26, medV0APMAlum);
1145     TGeoVolume *v0APMT6 = new TGeoVolumeAssembly("V0APMT6");
1146     TGeoTube *sV0APMTT6 = new TGeoTube("sV0APMTT6", 0., fV0APMTR4, fV0APMTB/2.);
1147     TGeoVolume *v0APMTT6 = new TGeoVolume("V0APMTT6", sV0APMTT6, medV0APMAlum);
1148     v0APMT6->SetLineColor(kV0AColorPMG);
1149     v0APMT26->SetLineColor(kV0AColorPMA);
1150     v0APMTT6->SetLineColor(kV0AColorPMA);
1151     rot = new TGeoRotation("rot", 90, 0, 180, 0, 90, 90);
1152     v0APMT6->AddNode(v0APMT16,1,rot);
1153     v0APMT6->AddNode(v0APMT26,1,rot);
1154     v0APMT6->AddNode(v0APMTT6,1,new TGeoCombiTrans(0,(fV0APMTH+fV0APMTB)/2.,0,rot));
1155     double autoShift6 = (fV0APMBWd-2*fV0APMBWdW)/4.;
1156     v0APM6->AddNode(v0APMT6, 1, new TGeoTranslation(-1.5*autoShift6, 0, 0));
1157     v0APM6->AddNode(v0APMT6, 2, new TGeoTranslation(-0.5*autoShift6, 0, 0));
1158     v0APM6->AddNode(v0APMT6, 3, new TGeoTranslation(+0.5*autoShift6, 0, 0));
1159     v0APM6->AddNode(v0APMT6, 4, new TGeoTranslation(+1.5*autoShift6, 0, 0));
1160
1161     /// PM
1162     rot = new TGeoRotation("rot");
1163     rot->RotateX(-90+30);
1164     rot->RotateY(0);
1165     rot->RotateZ(-65-3);
1166     double cosAngPMB6 = TMath::Cos(fV0APMBAng*TMath::DegToRad());
1167     double sinAngPMB6 = TMath::Sin(fV0APMBAng*TMath::DegToRad());
1168     double shiftZ6 = fV0APMBHt/2. * cosAngPMB6
1169       -   ( fV0ASciWd + 2 * fV0AOctWd + 2 * fV0APlaWd )/2.   -   fV0APMBTh/2. * sinAngPMB6;
1170     double shiftR6 = fV0AR6  + fV0AR1 + fV0AOctWd + fV0APlaAl/3.;
1171     v0ASec6->AddNode(v0APM6,1, new TGeoCombiTrans( -shiftR6*cos45-1.3, -shiftR6*sin45, shiftZ6, rot));
1172
1173     /// Support
1174     TGeoBBox *sV0ASuppbl = new TGeoBBox("sV0ASuppbl", 2.0, 18.13, fV0ASciWd/2.);
1175     TGeoVolume *v0ASuppbl = new TGeoVolume("V0ASuppbl", sV0ASuppbl, medV0ASup);
1176     v0ASuppbl->SetLineColor(kV0AColorOct);
1177     v0ASec6->AddNode(v0ASuppbl,1,new TGeoTranslation(-2.0,-63.64,0.0));
1178
1179     // Aluminium nail
1180     TGeoTube *sV0ANail61 = new TGeoTube("sV0ANail61", 0.0, 0.4, 5.09/2.);
1181     TGeoVolume *v0ANail61 = new TGeoVolume("V0ANail61", sV0ANail61, medV0APMAlum);
1182     v0ANail61->SetLineColor(kV0AColorPMA);// this is the color for aluminium
1183     v0ASec6->AddNode(v0ANail61,1,new TGeoTranslation(-0.51,-42.9,0.0));
1184     TGeoTube *sV0ANail62 = new TGeoTube("sV0ANail62", 0.0, 0.4, 5.09/2.);
1185     TGeoVolume *v0ANail62 = new TGeoVolume("V0ANail62", sV0ANail62, medV0APMAlum);
1186     v0ANail62->SetLineColor(kV0AColorPMA);
1187     v0ASec6->AddNode(v0ANail62,1,new TGeoTranslation(-30.05,-30.79,0.0));
1188
1189     // Adding sector to v0LE volume
1190     v0LE->AddNode(v0ASec6, 1);
1191
1192
1193      //Definition of sector 7
1194     TGeoVolume *v0ASec7 = new TGeoVolumeAssembly("V0ASec7");
1195
1196     /// For boolean sustraction
1197     double preShape7 = 0.2;
1198     for (int i=0;i<2;i++) {
1199     v0APts[0+8*i] = 0.0;                      v0APts[1+8*i] = -fV0AR0+fV0AFraWd/2.-preShape7;
1200     v0APts[2+8*i] = fV0AFraWd/2.;                    v0APts[3+8*i] = -fV0AR0+fV0AFraWd/2.-preShape7;
1201     v0APts[4+8*i] = fV0AFraWd/2.;                    v0APts[5+8*i] = -fV0AR4-fV0AFraWd/2.+preShape7;
1202     v0APts[6+8*i] = 0.0;                      v0APts[7+8*i] = -fV0AR4-fV0AFraWd/2.+preShape7;
1203     }
1204     new TGeoArb8("sV0ACha17",fV0ASciWd/1.5,v0APts);
1205     for (int i=0;i<2;i++) {
1206     v0APts[0+8*i] = fV0AR0*cos45-preShape7;
1207     v0APts[1+8*i] = -(fV0AR0-fV0AFraWd)*sin45+preShape7;
1208     v0APts[2+8*i] = (fV0AR0-fV0AFraWd/2.)*cos45-preShape7;
1209     v0APts[3+8*i] = -(fV0AR0-fV0AFraWd/2.)*sin45;
1210     v0APts[4+8*i] = (fV0AR4+fV0AFraWd/2.)*cos45+preShape7;
1211     v0APts[5+8*i] = -(fV0AR4+fV0AFraWd/2.)*sin45-2.*preShape7;
1212     v0APts[6+8*i] = (fV0AR4+fV0AFraWd)*cos45+preShape7;
1213     v0APts[7+8*i] = -fV0AR4*sin45-preShape7;
1214     }
1215     new TGeoArb8("sV0ACha27", fV0ASciWd/2.+2.*preShape, v0APts);
1216     new TGeoCompositeShape("sV0ACha127","sV0ACha17+sV0ACha27");
1217     new TGeoTube("sV0ANail17Hole", 0.0, 0.4, 1.65);
1218     TGeoTranslation *pos17 = new TGeoTranslation("pos17",0.51,-42.9,0.0);
1219     pos17->RegisterYourself();
1220     new TGeoTube("sV0ANail27Hole", 0.0, 0.4, 1.65);
1221     TGeoTranslation *pos27 = new TGeoTranslation("pos27",30.05,-30.79,0.0);
1222     pos27->RegisterYourself();
1223     new TGeoCompositeShape("sV0ANailsHoles7","sV0ANail17Hole:pos17+sV0ANail27Hole:pos27");
1224     new TGeoCompositeShape("sV0ACha7","sV0ACha127+sV0ANailsHoles7");
1225     new TGeoTubeSeg("sV0AFicR57", fV0AR4+fV0AFraWd/2., fV0AR4 + fV0AR0, fV0ASciWd/2.+2*preShape, 270.0, 315.0);
1226     new TGeoTube("sV0ANail1PlaInHole7", 0.0, 0.4, (fV0APlaWd-2*fV0APlaAl)/2.);
1227     new TGeoTube("sV0ANail1PlaOuHole7", 0.0, 0.4, (fV0APlaAl)/2.);
1228
1229     /// Frame
1230     TGeoVolume *v0AFra7 = new TGeoVolumeAssembly("V0AFra7");
1231     for (int i=0;i<2;i++) {
1232     v0APts[0+8*i] = 0.;                              v0APts[1+8*i] = -fV0AR0-fV0AFraWd/2.;
1233     v0APts[2+8*i] = fV0AFraWd/2.;                    v0APts[3+8*i] = -fV0AR0-fV0AFraWd/2.;
1234     v0APts[4+8*i] = fV0AFraWd/2.;                    v0APts[5+8*i] = -fV0AR4+fV0AFraWd/2.;
1235     v0APts[6+8*i] = 0.;                              v0APts[7+8*i] = -fV0AR4+fV0AFraWd/2.;
1236     }
1237     TGeoArb8 *sV0AFraB17 = new TGeoArb8("sV0AFraB17",fV0ASciWd/2.,v0APts);
1238     TGeoVolume *v0AFraB17 = new TGeoVolume("V0AFraB17",sV0AFraB17,medV0AFra);
1239     for (int i=0;i<2;i++) {
1240     v0APts[0+8*i] = fV0AR0*cos45;
1241     v0APts[1+8*i] = -(fV0AR0-fV0AFraWd)*sin45;
1242     v0APts[2+8*i] = (fV0AR0-fV0AFraWd/2.)*cos45;
1243     v0APts[3+8*i] = -(fV0AR0-fV0AFraWd/2.)*sin45;
1244     v0APts[4+8*i] = (fV0AR4+fV0AFraWd/2.)*cos45/2.;
1245     v0APts[5+8*i] = -(fV0AR4+fV0AFraWd/2.)*sin45/2.;
1246     v0APts[6+8*i] = (fV0AR4+fV0AFraWd)*cos45/2.;
1247     v0APts[7+8*i] = -fV0AR4*sin45/2.;
1248     }
1249     TGeoArb8 *sV0AFraB27 = new TGeoArb8("sV0AFraB27", fV0ASciWd/2., v0APts);
1250     TGeoVolume *v0AFraB27 = new TGeoVolume("V0AFraB27",sV0AFraB27,medV0AFra);
1251     v0AFraB17->SetLineColor(kV0AColorFra); v0AFraB27->SetLineColor(kV0AColorFra);
1252     v0AFra7->AddNode(v0AFraB17,1);
1253     v0AFra7->AddNode(v0AFraB27,1);  // Prefer 2 GeoObjects insted of 3 GeoMovements
1254     new TGeoTubeSeg( "sV0AFraR1b7", fV0AR0-fV0AFraWd/2.,
1255              fV0AR0+fV0AFraWd/2., fV0ASciWd/2., 270.0, 315.0);
1256     new TGeoTubeSeg( "sV0AFraR2b7", fV0AR1-fV0AFraWd/2.,
1257              fV0AR1+fV0AFraWd/2., fV0ASciWd/2., 270.0, 315.0);
1258     new TGeoTubeSeg( "sV0AFraR3b7", fV0AR2-fV0AFraWd/2.,
1259              fV0AR2+fV0AFraWd/2., fV0ASciWd/2., 270.0, 315.0);
1260     new TGeoTubeSeg( "sV0AFraR4b7", fV0AR3-fV0AFraWd/2.,
1261              fV0AR3+fV0AFraWd/2., fV0ASciWd/2., 270.0, 315.0);
1262     new TGeoTubeSeg( "sV0AFraR5b7", fV0AR4-fV0AFraWd/2.,
1263              fV0AR4+fV0AFraWd/2., fV0ASciWd/2., 270.0, 315.0);
1264     TGeoCompositeShape *sV0AFraR17 = new TGeoCompositeShape("sV0AFraR17","sV0AFraR1b7-sV0ACha7");
1265     TGeoCompositeShape *sV0AFraR27 = new TGeoCompositeShape("sV0AFraR27","sV0AFraR2b7-sV0ACha7");
1266     TGeoCompositeShape *sV0AFraR37 = new TGeoCompositeShape("sV0AFraR37","sV0AFraR3b7-sV0ACha7");
1267     TGeoCompositeShape *sV0AFraR47 = new TGeoCompositeShape("sV0AFraR47","sV0AFraR4b7-sV0ACha7");
1268     TGeoCompositeShape *sV0AFraR57 = new TGeoCompositeShape("sV0AFraR57","sV0AFraR5b7-sV0ACha7");
1269     TGeoVolume *v0AFraR17 = new TGeoVolume("V0AFraR17",sV0AFraR17,medV0AFra);
1270     TGeoVolume *v0AFraR27 = new TGeoVolume("V0AFraR27",sV0AFraR27,medV0AFra);
1271     TGeoVolume *v0AFraR37 = new TGeoVolume("V0AFraR37",sV0AFraR37,medV0AFra);
1272     TGeoVolume *v0AFraR47 = new TGeoVolume("V0AFraR47",sV0AFraR47,medV0AFra);
1273     TGeoVolume *v0AFraR57 = new TGeoVolume("V0AFraR57",sV0AFraR57,medV0AFra);
1274     v0AFraR17->SetLineColor(kV0AColorFra); v0AFraR27->SetLineColor(kV0AColorFra);
1275     v0AFraR37->SetLineColor(kV0AColorFra); v0AFraR47->SetLineColor(kV0AColorFra);
1276     v0AFraR57->SetLineColor(kV0AColorFra);
1277     v0AFra7->AddNode(v0AFraR17,1);
1278     v0AFra7->AddNode(v0AFraR27,1);
1279     v0AFra7->AddNode(v0AFraR37,1);
1280     v0AFra7->AddNode(v0AFraR47,1);
1281     v0AFra7->AddNode(v0AFraR57,1);
1282     v0ASec7->AddNode(v0AFra7,1);
1283
1284     /// Sensitive scintilator
1285     TGeoVolume *v0ASci7 = new TGeoVolumeAssembly("V0ASci7");
1286     new TGeoTubeSeg( "sV0AR1b7", fV0AR0+fV0AFraWd/2.,
1287              fV0AR1-fV0AFraWd/2., fV0ASciWd/2., 270.0, 315.0);
1288     new TGeoTubeSeg( "sV0AR2b7", fV0AR1+fV0AFraWd/2.,
1289              fV0AR2-fV0AFraWd/2., fV0ASciWd/2., 270.0, 315.0);
1290     new TGeoTubeSeg( "sV0AR3b7", fV0AR2+fV0AFraWd/2.,
1291              fV0AR3-fV0AFraWd/2., fV0ASciWd/2., 270.0, 315.0);
1292     new TGeoTubeSeg( "sV0AR4b7", fV0AR3+fV0AFraWd/2.,
1293              fV0AR4-fV0AFraWd/2., fV0ASciWd/2., 270.0, 315.0);
1294     TGeoCompositeShape *sV0AR17 = new TGeoCompositeShape("sV0AR17","sV0AR1b7-sV0ACha7");
1295     TGeoCompositeShape *sV0AR27 = new TGeoCompositeShape("sV0AR27","sV0AR2b7-sV0ACha7");
1296     TGeoCompositeShape *sV0AR37 = new TGeoCompositeShape("sV0AR37","sV0AR3b7-sV0ACha7");
1297     TGeoCompositeShape *sV0AR47 = new TGeoCompositeShape("sV0AR47","sV0AR4b7-sV0ACha7");
1298     TGeoVolume *v0L17 = new TGeoVolume("V0L17",sV0AR17,medV0ASci);
1299     TGeoVolume *v0L27 = new TGeoVolume("V0L27",sV0AR27,medV0ASci);
1300     TGeoVolume *v0L37 = new TGeoVolume("V0L37",sV0AR37,medV0ASci);
1301     TGeoVolume *v0L47 = new TGeoVolume("V0L47",sV0AR47,medV0ASci);
1302     v0L17->SetLineColor(kV0AColorSci); v0L27->SetLineColor(kV0AColorSci);
1303     v0L37->SetLineColor(kV0AColorSci); v0L47->SetLineColor(kV0AColorSci);
1304     v0ASci7->AddNode(v0L17,1);
1305     v0ASci7->AddNode(v0L27,1);
1306     v0ASci7->AddNode(v0L37,1);
1307     v0ASci7->AddNode(v0L47,1);
1308
1309     // Bunch of fibers
1310     for (int i=0;i<2;i++) {
1311     v0APts[0+8*i] = -10.0;            v0APts[1+8*i] = 13.1;
1312     v0APts[2+8*i] = 10.0;            v0APts[3+8*i] = 13.1;
1313     v0APts[4+8*i] = 8.0;            v0APts[5+8*i] = -29.0;
1314     v0APts[6+8*i] = -12.0;            v0APts[7+8*i] = -12.0;
1315     }
1316     new TGeoArb8("sV0AFib17", 0.01, v0APts);
1317     rot = new TGeoRotation("rot");
1318     rot->RotateX(-2.0);
1319     rot->RotateY(0.0);
1320     rot->RotateZ(248.0);
1321     TGeoCombiTrans *fib17pos = new TGeoCombiTrans("fib17pos", 40.0 - 3.3, -50.0 + 1.5, 0.5, rot);
1322     fib17pos->RegisterYourself();
1323     TGeoCompositeShape *sV0AFib17Hole = new TGeoCompositeShape("sV0AFib17Hole", "sV0AFib17:fib17pos-sV0AFicR57");
1324     TGeoVolume *v0AFib17Hole = new TGeoVolume("V0AFib17Hole",sV0AFib17Hole,medV0AFib);
1325     v0AFib17Hole->SetLineColor(kV0AColorFib);
1326     new TGeoArb8("sV0AFib27", 0.01, v0APts);
1327     rot = new TGeoRotation("rot");
1328     rot->RotateX(2.0);
1329     rot->RotateY(0.0);
1330     rot->RotateZ(248.0);
1331     TGeoCombiTrans *fib27pos = new TGeoCombiTrans("fib27pos", 40.0 - 3.3, -50.0 + 1.5, -0.5, rot);
1332     fib27pos->RegisterYourself();
1333     TGeoCompositeShape *sV0AFib27Hole = new TGeoCompositeShape("sV0AFib27Hole", "sV0AFib27:fib27pos-sV0AFicR57");
1334     TGeoVolume *v0AFib27Hole = new TGeoVolume("V0AFib27Hole",sV0AFib27Hole,medV0AFib);
1335     v0AFib27Hole->SetLineColor(kV0AColorFib);
1336     TGeoVolume *v0AFib7 = new TGeoVolumeAssembly("V0AFib7");
1337     v0AFib7->AddNode(v0AFib17Hole,1);
1338     v0AFib7->AddNode(v0AFib27Hole,1);
1339     v0ASec7->AddNode(v0AFib7,1);
1340
1341     /// Non-sensitive scintilator
1342     new TGeoTubeSeg("sV0AR5S27", fV0AR4+fV0AFraWd/2., fV0AR4 + fV0AR0, fV0ASciWd/2.+2*preShape, 270.0, 315.0);
1343     TGeoCompositeShape *sV0AR57 = new TGeoCompositeShape("V0AR57","sV0AR5S27 - sV0ACha7");
1344     TGeoVolume *v0AR57 = new TGeoVolume("V0AR57",sV0AR57,medV0ASci);
1345     v0AR57->SetLineColor(kV0AColorSci);
1346     v0ASci7->AddNode(v0AR57,1);
1347     v0ASec7->AddNode(v0ASci7,1);
1348
1349     /// Segment of octagon
1350     for (int i=0;i<2;i++) {
1351     v0APts[0+8*i] = 0.;                          v0APts[1+8*i] = -(fV0AR7-fV0AOctH2)*sin654;
1352     v0APts[2+8*i] = 0.;                          v0APts[3+8*i] = -fV0AR7*sin654;
1353     v0APts[4+8*i] = fV0AR7*cos654;               v0APts[5+8*i] = -fV0AR7*sin654;
1354     v0APts[6+8*i] = (fV0AR7-fV0AOctH2)*cos654;   v0APts[7+8*i] = -(fV0AR7-fV0AOctH2)*sin654;
1355     }
1356     TGeoArb8 *sV0AOct27 = new TGeoArb8("sV0AOct27", (fV0ASciWd+2*fV0AOctWd)/2., v0APts);
1357     TGeoVolume *v0AOct27 = new TGeoVolume("V0AOct27", sV0AOct27,medV0ASup);
1358     v0AOct27->SetLineColor(kV0AColorOct);
1359     TGeoVolume *v0ASup7 = new TGeoVolumeAssembly("V0ASup7");
1360     v0ASup7->AddNode(v0AOct27,1);
1361     v0ASec7->AddNode(v0ASup7,1);
1362
1363     /// Plates
1364     for (int i=0;i<2;i++) {
1365     v0APts[0+8*i] = 0.;                 v0APts[1+8*i] = -fV0AR0;
1366     v0APts[2+8*i] = 0.;                 v0APts[3+8*i] = -fV0AR7*sin654;
1367     v0APts[4+8*i] = fV0AR7*cos654;      v0APts[5+8*i] = -fV0AR7*sin654;
1368     v0APts[6+8*i] = fV0AR0*cos654;              v0APts[7+8*i] = -fV0AR0*sin654;
1369     }
1370     new TGeoArb8("sV0APlaIn7", (fV0APlaWd-2*fV0APlaAl)/2., v0APts);
1371     TGeoCompositeShape *sV0APlaInNailsHoles7 = new TGeoCompositeShape("sV0APlaInNailsHoles7","sV0APlaIn7-sV0ANail1PlaInHole7:pos17");
1372     TGeoVolume *v0APlaInNailsHoles7 = new TGeoVolume("V0APlaInNailsHoles7", sV0APlaInNailsHoles7, medV0APlaIn);
1373     new TGeoArb8("sV0APlaOu7", fV0APlaAl/2., v0APts);
1374     TGeoCompositeShape *sV0APlaOuNailsHoles7 = new TGeoCompositeShape("sV0APlaOuNailsHoles7","sV0APlaOu7-sV0ANail1PlaOuHole7:pos17");
1375     TGeoVolume *v0APlaOuNailsHoles7 = new TGeoVolume("V0APlaOuNailsHoles7", sV0APlaOuNailsHoles7, medV0APlaOu);
1376     v0APlaInNailsHoles7->SetLineColor(kV0AColorPlaIn); v0APlaOuNailsHoles7->SetLineColor(kV0AColorPlaOu);
1377     TGeoVolume *v0APla7 = new TGeoVolumeAssembly("V0APla7");
1378     v0APla7->AddNode(v0APlaInNailsHoles7,1);
1379     v0APla7->AddNode(v0APlaOuNailsHoles7,1,new TGeoTranslation(0,0,(fV0APlaWd-fV0APlaAl)/2.));
1380     v0APla7->AddNode(v0APlaOuNailsHoles7,2,new TGeoTranslation(0,0,-(fV0APlaWd-fV0APlaAl)/2.));
1381     v0ASec7->AddNode(v0APla7,1,new TGeoTranslation(0,0,(fV0ASciWd+2*fV0AOctWd+fV0APlaWd)/2.));
1382     v0ASec7->AddNode(v0APla7,2,new TGeoTranslation(0,0,-(fV0ASciWd+2*fV0AOctWd+fV0APlaWd)/2.));
1383
1384     /// PMBox
1385     TGeoVolume* v0APM7 = new TGeoVolumeAssembly("V0APM7");
1386     new TGeoBBox("sV0APMB17", fV0APMBWd/2., fV0APMBHt/2., fV0APMBTh/2.);
1387     new TGeoBBox("sV0APMB27", fV0APMBWd/2.-fV0APMBWdW, fV0APMBHt/2.-fV0APMBHtW, fV0APMBTh/2.-fV0APMBThW);
1388     TGeoCompositeShape *sV0APMB7 = new TGeoCompositeShape("sV0APMB7","sV0APMB17-sV0APMB27");
1389     TGeoVolume *v0APMB7 = new TGeoVolume("V0APMB7",sV0APMB7, medV0APMAlum);
1390     v0APMB7->SetLineColor(kV0AColorPMA);
1391     v0APM7->AddNode(v0APMB7,1);
1392
1393     /// PMTubes
1394     TGeoTube *sV0APMT17 = new TGeoTube("sV0APMT17", fV0APMTR1, fV0APMTR2, fV0APMTH/2.);
1395     TGeoVolume *v0APMT17 = new TGeoVolume("V0APMT17", sV0APMT17, medV0APMGlass);
1396     TGeoTube *sV0APMT27 = new TGeoTube("sV0APMT27", fV0APMTR3, fV0APMTR4, fV0APMTH/2.);
1397     TGeoVolume *v0APMT27 = new TGeoVolume("V0APMT27", sV0APMT27, medV0APMAlum);
1398     TGeoVolume *v0APMT7 = new TGeoVolumeAssembly("V0APMT7"); // pk si no choca con la 752 o con la 794
1399     TGeoTube *sV0APMTT7 = new TGeoTube("sV0APMTT7", 0., fV0APMTR4, fV0APMTB/2.);
1400     TGeoVolume *v0APMTT7 = new TGeoVolume("V0APMTT7", sV0APMTT7, medV0APMAlum);
1401     v0APMT7->SetLineColor(kV0AColorPMG);
1402     v0APMT27->SetLineColor(kV0AColorPMA);
1403     v0APMTT7->SetLineColor(kV0AColorPMA);
1404     rot = new TGeoRotation("rot", 90, 0, 180, 0, 90, 90);
1405     v0APMT7->AddNode(v0APMT17,1,rot);
1406     v0APMT7->AddNode(v0APMT27,1,rot);
1407     v0APMT7->AddNode(v0APMTT7,1,new TGeoCombiTrans(0,(fV0APMTH+fV0APMTB)/2.,0,rot));
1408     v0APM7->AddNode(v0APMT7, 1, new TGeoTranslation(-1.5*autoShift, 0, 0));
1409     v0APM7->AddNode(v0APMT7, 2, new TGeoTranslation(-0.5*autoShift, 0, 0));
1410     v0APM7->AddNode(v0APMT7, 3, new TGeoTranslation(+0.5*autoShift, 0, 0));
1411     v0APM7->AddNode(v0APMT7, 4, new TGeoTranslation(+1.5*autoShift, 0, 0));
1412
1413     /// PM
1414     rot = new TGeoRotation("rot");
1415     rot->RotateX(-90+30);
1416     rot->RotateY(0);
1417     rot->RotateZ(65+3);
1418     double shiftZ7 = fV0APMBHt/2. * cosAngPMB
1419       -   ( fV0ASciWd + 2 * fV0AOctWd + 2 * fV0APlaWd )/2.   -   fV0APMBTh/2. * sinAngPMB;
1420     double shiftR7 = fV0AR6  + fV0AR1 + fV0AOctWd + fV0APlaAl/3.;
1421     v0ASec7->AddNode(v0APM7,1, new TGeoCombiTrans( shiftR7*cos45+1.3, -shiftR7*sin45, shiftZ7, rot));
1422
1423     // Aluminium nail
1424     TGeoTube *sV0ANail71 = new TGeoTube("sV0ANail71", 0.0, 0.4, 5.09/2.);
1425     TGeoVolume *v0ANail71 = new TGeoVolume("V0ANail71", sV0ANail71, medV0APMAlum);
1426     v0ANail71->SetLineColor(kV0AColorPMA);// this is the color for aluminium
1427     v0ASec7->AddNode(v0ANail71,1,new TGeoTranslation(0.51,-42.9,0.0));
1428     TGeoTube *sV0ANail72 = new TGeoTube("sV0ANail72", 0.0, 0.4, 5.09/2.);
1429     TGeoVolume *v0ANail72 = new TGeoVolume("V0ANail72", sV0ANail72, medV0APMAlum);
1430     v0ANail72->SetLineColor(kV0AColorPMA);
1431     v0ASec7->AddNode(v0ANail72,1,new TGeoTranslation(30.05,-30.79,0.0));
1432
1433     // Support
1434     TGeoBBox *sV0ASuppbr = new TGeoBBox("sV0ASuppbr", 2.0, 18.13, fV0ASciWd/2.);
1435     TGeoVolume *v0ASuppbr = new TGeoVolume("V0ASuppbr", sV0ASuppbr, medV0ASup);
1436     v0ASuppbr->SetLineColor(kV0AColorOct);
1437     v0ASec7->AddNode(v0ASuppbr,1,new TGeoTranslation(2.0,-63.64,0.0));
1438
1439     // Adding sector to v0LE volume
1440     v0LE->AddNode(v0ASec7,1);
1441
1442
1443    //Definition of sector 8
1444    TGeoVolume *v0ASec8 = new TGeoVolumeAssembly("V0ASec8");
1445
1446   /// For boolean sustraction
1447       for (int i=0;i<2;i++) {
1448       v0APts[0+8*i] = fV0AR0-fV0AFraWd/2.+preShape;  v0APts[1+8*i] = 0.0;
1449       v0APts[2+8*i] = fV0AR0-fV0AFraWd/2.+preShape;  v0APts[3+8*i] = -preShape;
1450       v0APts[4+8*i] = fV0AR4+fV0AFraWd/2.-preShape;  v0APts[5+8*i] = -preShape;
1451       v0APts[6+8*i] = fV0AR4+fV0AFraWd/2.-preShape;  v0APts[7+8*i] = 0.0;
1452     }
1453     new TGeoArb8("sV0ACha18",fV0ASciWd/1.5,v0APts);
1454     for (int i=0;i<2;i++) {
1455       v0APts[0+8*i] = fV0AR0*sin45-fV0AFraWd/2.-preShape;
1456       v0APts[1+8*i] = -(fV0AR0-fV0AFraWd)*sin45+preShape;
1457       v0APts[2+8*i] = (fV0AR0-fV0AFraWd/2.)*sin45-preShape;
1458       v0APts[3+8*i] = -(fV0AR0-fV0AFraWd/2.)*sin45;
1459       v0APts[4+8*i] = (fV0AR4+fV0AFraWd/2.)*sin45+preShape;
1460       v0APts[5+8*i] = -(fV0AR4+fV0AFraWd/2.)*sin45-2.*preShape;
1461       v0APts[6+8*i] = (fV0AR4+fV0AFraWd)*sin45+preShape;
1462       v0APts[7+8*i] = -fV0AR4*sin45-preShape;
1463     }
1464     new TGeoArb8("sV0ACha28", fV0ASciWd/2.+2.*preShape, v0APts);
1465     new TGeoCompositeShape("sV0ACha128","sV0ACha18+sV0ACha28");
1466     new TGeoTube("sV0ANail18Hole", 0.0, 0.4, 1.65);
1467     TGeoTranslation *pos18 = new TGeoTranslation("pos18",42.9,-.51,0.0);
1468     pos18->RegisterYourself();
1469     new TGeoTube("sV0ANail28Hole", 0.0, 0.4, 1.65);
1470     TGeoTranslation *pos28 = new TGeoTranslation("pos28",30.8,-30.04,0.0);
1471     pos28->RegisterYourself();
1472     TGeoTranslation *pos38 = new TGeoTranslation("pos38",30.05,-30.79,0.0);
1473     pos38->RegisterYourself();
1474     new TGeoCompositeShape("sV0ANailsHoles8","sV0ANail18Hole:pos18+sV0ANail28Hole:pos28");
1475     new TGeoCompositeShape("sV0ACha8","sV0ACha128+sV0ANailsHoles8");
1476     new TGeoTubeSeg("sV0AFicR58", fV0AR4+fV0AFraWd/2., fV0AR4 + fV0AR0, fV0ASciWd/2. +2*preShape, 315.0, 360.0);
1477     new TGeoTube("sV0ANail1PlaInHole8", 0.0, 0.4, (fV0APlaWd-2*fV0APlaAl)/2.);
1478     new TGeoTube("sV0ANail2PlaInHole8", 0.0, 0.4, (fV0APlaWd-2*fV0APlaAl)/2.);
1479     new TGeoTube("sV0ANail3PlaInHole8", 0.0, 0.4, (fV0APlaWd-2*fV0APlaAl)/2.);
1480     new TGeoCompositeShape("sV0ANailsPlaInHoles8","sV0ANail1PlaInHole8:pos18+sV0ANail2PlaInHole8:pos28+sV0ANail3PlaInHole8:pos38");
1481     new TGeoTube("sV0ANail1PlaOuHole8", 0.0, 0.4, (fV0APlaAl)/2.);
1482     new TGeoTube("sV0ANail2PlaOuHole8", 0.0, 0.4, (fV0APlaAl)/2.);
1483     new TGeoTube("sV0ANail3PlaOuHole8", 0.0, 0.4, (fV0APlaAl)/2.);
1484     new TGeoCompositeShape("sV0ANailsPlaOuHoles8","sV0ANail1PlaOuHole8:pos18+sV0ANail2PlaOuHole8:pos28+sV0ANail3PlaOuHole8:pos38");
1485
1486     /// Frame
1487     TGeoVolume *v0AFra8 = new TGeoVolumeAssembly("V0AFra8");
1488     for (int i=0;i<2;i++) {
1489       v0APts[0+8*i] = fV0AR0-fV0AFraWd/2.;  v0APts[1+8*i] = 0.0;
1490       v0APts[2+8*i] = fV0AR0-fV0AFraWd/2.;  v0APts[3+8*i] = -fV0AFraWd/8.;
1491       v0APts[4+8*i] = fV0AR4+fV0AFraWd/2.;  v0APts[5+8*i] = -fV0AFraWd/8.;
1492       v0APts[6+8*i] = fV0AR4+fV0AFraWd/2.;  v0APts[7+8*i] = 0.0;
1493     }
1494     TGeoArb8 *sV0AFraB18 = new TGeoArb8("sV0AFraB18",fV0ASciWd/2.,v0APts);
1495     TGeoVolume *v0AFraB18 = new TGeoVolume("V0AFraB18",sV0AFraB18,medV0AFra);
1496     for (int i=0;i<2;i++) {
1497       v0APts[0+8*i] = fV0AR0*sin45-fV0AFraWd/2.;
1498       v0APts[1+8*i] = -(fV0AR0-fV0AFraWd)*sin45;
1499       v0APts[2+8*i] = (fV0AR0-fV0AFraWd/2.)*sin45;
1500       v0APts[3+8*i] = -(fV0AR0-fV0AFraWd/2.)*sin45;
1501       v0APts[4+8*i] = (fV0AR4+fV0AFraWd/2.)*sin45/2.;
1502       v0APts[5+8*i] = -(fV0AR4+fV0AFraWd/2.)*sin45/2.;
1503       v0APts[6+8*i] = (fV0AR4+fV0AFraWd)*sin45/2.;
1504       v0APts[7+8*i] = -fV0AR4*sin45/2.;
1505     }
1506     TGeoArb8 *sV0AFraB28 = new TGeoArb8("sV0AFraB28", fV0ASciWd/2., v0APts);
1507     TGeoVolume *v0AFraB28 = new TGeoVolume("V0AFraB28",sV0AFraB28,medV0AFra);
1508     v0AFraB18->SetLineColor(kV0AColorFra); v0AFraB28->SetLineColor(kV0AColorFra);
1509     v0AFra8->AddNode(v0AFraB18,1);
1510     v0AFra8->AddNode(v0AFraB28,1);  // Prefer 2 GeoObjects insted of 3 GeoMovements
1511     new TGeoTubeSeg( "sV0AFraR1b8", fV0AR0-fV0AFraWd/2.,
1512                      fV0AR0+fV0AFraWd/2., fV0ASciWd/2., 315.0, 360.0);
1513     new TGeoTubeSeg( "sV0AFraR2b8", fV0AR1-fV0AFraWd/2.,
1514                      fV0AR1+fV0AFraWd/2., fV0ASciWd/2., 315.0, 360.0);
1515     new TGeoTubeSeg( "sV0AFraR3b8", fV0AR2-fV0AFraWd/2.,
1516                      fV0AR2+fV0AFraWd/2., fV0ASciWd/2., 315.0, 360.0);
1517     new TGeoTubeSeg( "sV0AFraR4b8", fV0AR3-fV0AFraWd/2.,
1518                      fV0AR3+fV0AFraWd/2., fV0ASciWd/2., 315.0, 360.0);
1519     new TGeoTubeSeg( "sV0AFraR5b8", fV0AR4-fV0AFraWd/2.,
1520                      fV0AR4+fV0AFraWd/2., fV0ASciWd/2., 315.0, 360.0);
1521     TGeoCompositeShape *sV0AFraR18 = new TGeoCompositeShape("sV0AFraR18","sV0AFraR1b8-sV0ACha8");
1522     TGeoCompositeShape *sV0AFraR28 = new TGeoCompositeShape("sV0AFraR28","sV0AFraR2b8-sV0ACha8");
1523     TGeoCompositeShape *sV0AFraR38 = new TGeoCompositeShape("sV0AFraR38","sV0AFraR3b8-sV0ACha8");
1524     TGeoCompositeShape *sV0AFraR48 = new TGeoCompositeShape("sV0AFraR48","sV0AFraR4b8-sV0ACha8");
1525     TGeoCompositeShape *sV0AFraR58 = new TGeoCompositeShape("sV0AFraR58","sV0AFraR5b8-sV0ACha8");
1526     TGeoVolume *v0AFraR18 = new TGeoVolume("V0AFraR18",sV0AFraR18,medV0AFra);
1527     TGeoVolume *v0AFraR28 = new TGeoVolume("V0AFraR28",sV0AFraR28,medV0AFra);
1528     TGeoVolume *v0AFraR38 = new TGeoVolume("V0AFraR38",sV0AFraR38,medV0AFra);
1529     TGeoVolume *v0AFraR48 = new TGeoVolume("V0AFraR48",sV0AFraR48,medV0AFra);
1530     TGeoVolume *v0AFraR58 = new TGeoVolume("V0AFraR58",sV0AFraR58,medV0AFra);
1531     v0AFraR18->SetLineColor(kV0AColorFra); v0AFraR28->SetLineColor(kV0AColorFra);
1532     v0AFraR38->SetLineColor(kV0AColorFra); v0AFraR48->SetLineColor(kV0AColorFra);
1533     v0AFraR58->SetLineColor(kV0AColorFra);
1534     v0AFra8->AddNode(v0AFraR18,1);
1535     v0AFra8->AddNode(v0AFraR28,1);
1536     v0AFra8->AddNode(v0AFraR38,1);
1537     v0AFra8->AddNode(v0AFraR48,1);
1538     v0AFra8->AddNode(v0AFraR58,1);
1539     v0ASec8->AddNode(v0AFra8,1);
1540
1541     /// Sensitive scintilator
1542     TGeoVolume *v0ASci8 = new TGeoVolumeAssembly("V0ASci8");
1543     new TGeoTubeSeg( "sV0AR1b8", fV0AR0+fV0AFraWd/2.,
1544                      fV0AR1-fV0AFraWd/2., fV0ASciWd/2., 315.0, 360.0);
1545     new TGeoTubeSeg( "sV0AR2b8", fV0AR1+fV0AFraWd/2.,
1546                      fV0AR2-fV0AFraWd/2., fV0ASciWd/2., 315.0, 360.0);
1547     new TGeoTubeSeg( "sV0AR3b8", fV0AR2+fV0AFraWd/2.,
1548                      fV0AR3-fV0AFraWd/2., fV0ASciWd/2., 315.0, 360.0);
1549     new TGeoTubeSeg( "sV0AR4b8", fV0AR3+fV0AFraWd/2.,
1550                      fV0AR4-fV0AFraWd/2., fV0ASciWd/2., 315.0, 360.0);
1551     TGeoCompositeShape *sV0AR18 = new TGeoCompositeShape("sV0AR18","sV0AR1b8-sV0ACha8");
1552     TGeoCompositeShape *sV0AR28 = new TGeoCompositeShape("sV0AR28","sV0AR2b8-sV0ACha8");
1553     TGeoCompositeShape *sV0AR38 = new TGeoCompositeShape("sV0AR38","sV0AR3b8-sV0ACha8");
1554     TGeoCompositeShape *sV0AR48 = new TGeoCompositeShape("sV0AR48","sV0AR4b8-sV0ACha8");
1555     TGeoVolume *v0L18 = new TGeoVolume("V0L18",sV0AR18,medV0ASci);
1556     TGeoVolume *v0L28 = new TGeoVolume("V0L28",sV0AR28,medV0ASci);
1557     TGeoVolume *v0L38 = new TGeoVolume("V0L38",sV0AR38,medV0ASci);
1558     TGeoVolume *v0L48 = new TGeoVolume("V0L48",sV0AR48,medV0ASci);
1559     v0L18->SetLineColor(kV0AColorSci); v0L28->SetLineColor(kV0AColorSci);
1560     v0L38->SetLineColor(kV0AColorSci); v0L48->SetLineColor(kV0AColorSci);
1561     v0ASci8->AddNode(v0L18,1);
1562     v0ASci8->AddNode(v0L28,1);
1563     v0ASci8->AddNode(v0L38,1);
1564     v0ASci8->AddNode(v0L48,1);
1565
1566     /// Segment of octagon
1567     for (int i=0;i<2;i++) {
1568     v0APts[0+8*i] = fV0AR6-fV0AOctH2;            v0APts[1+8*i] = 0.;
1569     v0APts[2+8*i] = (fV0AR7-fV0AOctH2)*cos654;   v0APts[3+8*i] = -(fV0AR7-fV0AOctH2)*sin654;
1570     v0APts[4+8*i] = fV0AR7*cos654;               v0APts[5+8*i] = -fV0AR7*sin654;
1571     v0APts[6+8*i] = fV0AR6;                      v0APts[7+8*i] = 0.;
1572     }
1573     TGeoArb8 *sV0AOct28 = new TGeoArb8("sV0AOct28", (fV0ASciWd+2*fV0AOctWd)/2., v0APts);
1574     TGeoVolume *v0AOct28 = new TGeoVolume("V0AOct28", sV0AOct28,medV0ASup);
1575     v0AOct28->SetLineColor(kV0AColorOct);
1576     TGeoVolume *v0ASup8 = new TGeoVolumeAssembly("V0ASup8");
1577     v0ASup8->AddNode(v0AOct28,1);
1578     v0ASec8->AddNode(v0ASup8,1);
1579
1580     //Bunch of fibers
1581     v0APts[ 0] = v0APts[ 2] = -14.0;
1582     v0APts[ 1] = v0APts[ 7] = (fV0ASciWd+fV0AOctWd)/2.-0.01;
1583     v0APts[ 3] = v0APts[ 5] = (fV0ASciWd+fV0AOctWd)/2.+0.01;
1584     v0APts[ 4] = v0APts[ 6] = +14.0;
1585     v0APts[ 8] = v0APts[10] = -10.0;
1586     v0APts[ 9] = v0APts[15] = 0.;
1587     v0APts[11] = v0APts[13] = 0.25;
1588     v0APts[12] = v0APts[14] = +10.0;
1589     new TGeoArb8("sV0AFib18", 11.8, v0APts);
1590     rot = new TGeoRotation("rot");
1591     rot->RotateX(-90);
1592     rot->RotateZ(-90-22.5);
1593     TGeoCombiTrans *fib18pos = new TGeoCombiTrans("fib18pos", (fV0AR6-fV0AOctH2+fV0AR5)*cos225/2. - 3.3, -(fV0AR6-fV0AOctH2+fV0AR5)*sin225/2. + 1.5, 0, rot);
1594     fib18pos->RegisterYourself();
1595     TGeoCompositeShape *sV0AFib18Hole = new TGeoCompositeShape("sV0AFib18Hole", "sV0AFib18:fib18pos-sV0AFicR58");
1596     TGeoVolume *v0AFib18Hole = new TGeoVolume("V0AFib18",sV0AFib18Hole,medV0AFib);
1597     v0AFib18Hole->SetLineColor(kV0AColorFib);
1598     new TGeoArb8("sV0AFib28", 11.8, v0APts);
1599     rot = new TGeoRotation("rot");
1600     rot->RotateX(-90);
1601     rot->RotateY(180);
1602     rot->RotateZ(-90-22.5);
1603     TGeoCombiTrans *fib28pos = new TGeoCombiTrans("fib28pos", (fV0AR6-fV0AOctH2+fV0AR5)*cos225/2. - 3.3, -(fV0AR6-fV0AOctH2+fV0AR5)*sin225/2. + 1.5, 0, rot);
1604     fib28pos->RegisterYourself();
1605     TGeoCompositeShape *sV0AFib28Hole = new TGeoCompositeShape("sV0AFib28Hole", "sV0AFib28:fib28pos-sV0AFicR58");
1606     TGeoVolume *v0AFib28Hole = new TGeoVolume("V0AFib28Hole",sV0AFib28Hole,medV0AFib);
1607     v0AFib28Hole->SetLineColor(kV0AColorFib);
1608     TGeoVolume *v0AFib8 = new TGeoVolumeAssembly("V0AFib8");
1609     v0AFib8->AddNode(v0AFib18Hole,1);
1610     v0AFib8->AddNode(v0AFib28Hole,1);
1611     v0ASec8->AddNode(v0AFib8,1);
1612
1613     /// Non-sensitive scintilator
1614     new TGeoTubeSeg("sV0AR5S28", fV0AR4+fV0AFraWd/2., fV0AR4 + fV0AR0, fV0ASciWd/2.+2*preShape, 315.0, 360.0);
1615     TGeoCompositeShape *sV0AR58 = new TGeoCompositeShape("V0AR58","sV0AR5S28 - sV0ACha8");
1616     TGeoVolume *v0AR58 = new TGeoVolume("V0AR58",sV0AR58,medV0ASci);
1617     v0AR58->SetLineColor(kV0AColorSci);
1618     v0ASci8->AddNode(v0AR58,1);
1619     v0ASec8->AddNode(v0ASci8,1);
1620
1621     /// Plates
1622     for (int i=0;i<2;i++) {
1623       v0APts[0+8*i] = fV0AR0;                   v0APts[1+8*i] = 0.;
1624       v0APts[2+8*i] = fV0AR0*cos654;            v0APts[3+8*i] = -fV0AR0*sin654;
1625       v0APts[4+8*i] = fV0AR7*cos654;            v0APts[5+8*i] = -fV0AR7*sin654;
1626       v0APts[6+8*i] = fV0AR6;                   v0APts[7+8*i] = 0.;
1627     }
1628     new TGeoArb8("sV0APlaIn8", (fV0APlaWd-2*fV0APlaAl)/2., v0APts);
1629     TGeoCompositeShape *sV0APlaInNailsHoles8 = new TGeoCompositeShape("sV0APlaInNailsHoles8","sV0APlaIn8-sV0ANailsPlaInHoles8");
1630     TGeoVolume *v0APlaInNailsHoles8 = new TGeoVolume("V0APlaInNailsHoles8", sV0APlaInNailsHoles8, medV0APlaIn);
1631     new TGeoArb8("sV0APlaOu8", fV0APlaAl/2., v0APts);
1632     TGeoCompositeShape *sV0APlaOuNailsHoles8 = new TGeoCompositeShape("sV0APlaOuNailsHoles8","sV0APlaOu8-sV0ANailsPlaOuHoles8");
1633     TGeoVolume *v0APlaOuNailsHoles8 = new TGeoVolume("V0APlaOuNailsHoles8", sV0APlaOuNailsHoles8, medV0APlaOu);
1634     v0APlaInNailsHoles8->SetLineColor(kV0AColorPlaIn); v0APlaOuNailsHoles8->SetLineColor(kV0AColorPlaOu);
1635     TGeoVolume *v0APla8 = new TGeoVolumeAssembly("V0APla8");
1636     v0APla8->AddNode(v0APlaInNailsHoles8,1);
1637     v0APla8->AddNode(v0APlaOuNailsHoles8,1,new TGeoTranslation(0,0,(fV0APlaWd-fV0APlaAl)/2.));
1638     v0APla8->AddNode(v0APlaOuNailsHoles8,2,new TGeoTranslation(0,0,-(fV0APlaWd-fV0APlaAl)/2.));
1639     v0ASec8->AddNode(v0APla8,1,new TGeoTranslation(0,0,(fV0ASciWd+2*fV0AOctWd+fV0APlaWd)/2.));
1640     v0ASec8->AddNode(v0APla8,2,new TGeoTranslation(0,0,-(fV0ASciWd+2*fV0AOctWd+fV0APlaWd)/2.));
1641
1642     /// PMBox
1643     TGeoVolume* v0APM8 = new TGeoVolumeAssembly("V0APM1");
1644     new TGeoBBox("sV0APMB18", fV0APMBWd/2., fV0APMBHt/2., fV0APMBTh/2.);
1645     new TGeoBBox("sV0APMB28", fV0APMBWd/2.-fV0APMBWdW, fV0APMBHt/2.-fV0APMBHtW, fV0APMBTh/2.-fV0APMBThW);
1646     TGeoCompositeShape *sV0APMB8 = new TGeoCompositeShape("sV0APMB8","sV0APMB18-sV0APMB28");
1647     TGeoVolume *v0APMB8 = new TGeoVolume("V0APMB8",sV0APMB8, medV0APMAlum);
1648     v0APMB8->SetLineColor(kV0AColorPMA);
1649     v0APM8->AddNode(v0APMB8,1);
1650
1651     /// PMTubes
1652     TGeoTube *sV0APMT18 = new TGeoTube("sV0APMT18", fV0APMTR1, fV0APMTR2, fV0APMTH/2.);
1653     TGeoVolume *v0APMT18 = new TGeoVolume("V0APMT18", sV0APMT18, medV0APMGlass);
1654     TGeoTube *sV0APMT28 = new TGeoTube("sV0APMT28", fV0APMTR3, fV0APMTR4, fV0APMTH/2.);
1655     TGeoVolume *v0APMT28 = new TGeoVolume("V0APMT28", sV0APMT28, medV0APMAlum);
1656     TGeoVolume *v0APMT8 = new TGeoVolumeAssembly("V0APMT8");
1657     TGeoTube *sV0APMTT8 = new TGeoTube("sV0APMTT8", 0., fV0APMTR4, fV0APMTB/2.);
1658     TGeoVolume *v0APMTT8 = new TGeoVolume("V0APMTT8", sV0APMTT8, medV0APMAlum);
1659     v0APMT8->SetLineColor(kV0AColorPMG);
1660     v0APMT28->SetLineColor(kV0AColorPMA);
1661     v0APMTT8->SetLineColor(kV0AColorPMA);
1662     rot = new TGeoRotation("rot", 90, 0, 180, 0, 90, 90);
1663     v0APMT8->AddNode(v0APMT18,1,rot);
1664     v0APMT8->AddNode(v0APMT28,1,rot);
1665     v0APMT8->AddNode(v0APMTT8,1,new TGeoCombiTrans(0,(fV0APMTH+fV0APMTB)/2.,0,rot));
1666     double autoShift8 = (fV0APMBWd-2*fV0APMBWdW)/4.;
1667     v0APM8->AddNode(v0APMT8, 1, new TGeoTranslation(-1.5*autoShift8, 0, 0));
1668     v0APM8->AddNode(v0APMT8, 2, new TGeoTranslation(-0.5*autoShift8, 0, 0));
1669     v0APM8->AddNode(v0APMT8, 3, new TGeoTranslation(+0.5*autoShift8, 0, 0));
1670     v0APM8->AddNode(v0APMT8, 4, new TGeoTranslation(+1.5*autoShift8, 0, 0));
1671
1672     /// PM
1673     rot = new TGeoRotation("rot");
1674     rot->RotateX(-90+30);
1675     rot->RotateY(0);
1676     rot->RotateZ(65+3);
1677     double shiftZ8 = fV0APMBHt/2. * cosAngPMB
1678       -   ( fV0ASciWd + 2 * fV0AOctWd + 2 * fV0APlaWd )/2.   -   fV0APMBTh/2. * sinAngPMB;
1679     double shiftR8 = fV0AR6  +  fV0AOctH2 + fV0APlaAl;
1680     v0ASec8->AddNode(v0APM8,1, new TGeoCombiTrans( shiftR8*cos225 + 1.3, -shiftR8*sin225, shiftZ8, rot));
1681
1682     // Aluminium nails
1683     TGeoTube *sV0ANail81 = new TGeoTube("sV0ANail81", 0.0, 0.4, 5.09/2.);
1684     TGeoVolume *v0ANail81 = new TGeoVolume("V0ANail81", sV0ANail81, medV0APMAlum);
1685     v0ANail81->SetLineColor(kV0AColorPMA);// this is the color for aluminium
1686     v0ASec8->AddNode(v0ANail81,1,new TGeoTranslation(42.9,-.51,0.0));
1687     TGeoTube *sV0ANail82 = new TGeoTube("sV0ANail82", 0.0, 0.4, 5.09/2.);
1688     TGeoVolume *v0ANail82 = new TGeoVolume("V0ANail82", sV0ANail82, medV0APMAlum);
1689     v0ANail82->SetLineColor(kV0AColorPMA);
1690     v0ASec8->AddNode(v0ANail82,1,new TGeoTranslation(30.8,-30.04,0.0));
1691
1692     // Adding sector to v0LE volume
1693     v0LE->AddNode(v0ASec8, 1);
1694
1695
1696     // Adding detectors to top volume
1697     TGeoVolume *vZERO = new TGeoVolumeAssembly("VZERO");
1698     vZERO->AddNode(v0RI,1,new TGeoTranslation(0, 0, -zdet));
1699     vZERO->AddNode(v0LE,1,new TGeoTranslation(0, 0, +327.5));
1700     top->AddNode(vZERO,1);
1701 }
1702
1703 //_____________________________________________________________________________
1704 void AliVZEROv7::AddAlignableVolumes() const
1705 {
1706   //
1707   // Create entries for alignable volumes associating the symbolic volume
1708   // name with the corresponding volume path. Needs to be syncronized with
1709   // eventual changes in the geometry.
1710   //
1711   TString vpC = "/ALIC_1/VZERO_1/V0RI_1";
1712   TString vpA = "/ALIC_1/VZERO_1/V0LE_1";
1713   TString snC = "VZERO/V0C";
1714   TString snA = "VZERO/V0A";
1715
1716   if(!gGeoManager->SetAlignableEntry(snC.Data(),vpC.Data()))
1717     AliFatal(Form("Alignable entry %s not created. Volume path %s not valid", snC.Data(),vpC.Data()));
1718   if(!gGeoManager->SetAlignableEntry(snA.Data(),vpA.Data()))
1719     AliFatal(Form("Alignable entry %s not created. Volume path %s not valid", snA.Data(),vpA.Data()));
1720
1721 }
1722
1723 //_____________________________________________________________________________
1724 void AliVZEROv7::CreateMaterials()
1725 {
1726
1727 // Creates materials used for geometry
1728
1729   AliDebug(2,"Create materials");
1730   // Parameters for simulation scope
1731   Int_t     fieldType       = ((AliMagF*)TGeoGlobalMagField::Instance()->GetField())->Integ();     // Field type
1732   Double_t  maxField        = ((AliMagF*)TGeoGlobalMagField::Instance()->GetField())->Max();       // Field max.
1733   Double_t  maxBending      = 10;    // Max Angle
1734   Double_t  maxStepSize     = 0.01;  // Max step size
1735   Double_t  maxEnergyLoss   = 1;     // Max Delta E
1736   Double_t  precision       = 0.003; // Precision
1737   Double_t  minStepSize     = 0.003; // Minimum step size
1738
1739   Int_t    id;
1740   Double_t a, z, radLength, absLength;
1741   Float_t density, as[4], zs[4], ws[4];
1742
1743 // Parameters  for V0CPrePlates: Aluminium
1744    a = 26.98;
1745    z = 13.00;
1746    density     = 2.7;
1747    radLength   = 8.9;
1748    absLength   = 37.2;
1749    id = 2;
1750    AliMaterial( id, "V0CAlu", a, z, density, radLength, absLength, 0, 0);
1751    AliMedium(id, "V0CAlu", id, 1, fieldType, maxField, maxBending, maxStepSize,
1752              maxEnergyLoss, precision, minStepSize);
1753
1754 // Parameters  for V0CPlates: Carbon
1755    a = 12.01;
1756    z =  6.00;
1757    density   = 2.265;
1758    radLength = 18.8;
1759    absLength = 49.9;
1760    id = 3;
1761    AliMaterial(id, "V0CCar",  a, z, density, radLength, absLength, 0, 0);
1762    AliMedium(id, "V0CCar", id, 1, fieldType, maxField, maxBending, maxStepSize,
1763              maxEnergyLoss, precision, minStepSize);
1764
1765 // Parameters  for V0Cscintillator: BC408
1766    as[0] = 1.00794;     as[1] = 12.011;
1767    zs[0] = 1.;          zs[1] = 6.;
1768    ws[0] = 1.;          ws[1] = 1.;
1769    density      = 1.032;
1770    id           = 4;
1771    AliMixture(id, "V0CSci", as, zs, density, -2, ws);
1772    AliMedium(id,"V0CSci", id, 1, fieldType, maxField, maxBending, maxStepSize,
1773              maxEnergyLoss, precision, minStepSize);
1774
1775 // Parameters for V0Ascintilator: BC404
1776    as[0] = 1.00794;     as[1] = 12.011;
1777    zs[0] = 1.;          zs[1] = 6.;
1778    ws[0] = 5.21;        ws[1] = 4.74;
1779    density      = 1.032;
1780    id           = 5;
1781    AliMixture(id, "V0ASci", as, zs, density, -2, ws);
1782    AliMedium(id,  "V0ASci", id, 1, fieldType, maxField, maxBending, maxStepSize,
1783              maxEnergyLoss, precision, minStepSize);
1784
1785 // Parameters for V0ALuc: Lucita but for the simulation BC404
1786    as[0] = 1.00794;     as[1] = 12.011;
1787    zs[0] = 1.;          zs[1] = 6.;
1788    ws[0] = 5.21;        ws[1] = 4.74;
1789    density      = 1.032;
1790    id           = 6;
1791    AliMixture(id, "V0ALuc", as, zs, density, -2, ws);
1792    AliMedium(id, "V0ALuc", id, 1, fieldType, maxField, maxBending, maxStepSize,
1793              maxEnergyLoss, precision, minStepSize);
1794
1795 // Parameters for V0Aplate: EuroComposite - EC-PI 626 PS - AlMg3
1796    as[0] = 26.982;      as[1] = 24.305;
1797    zs[0] = 13.;         zs[1] = 12.;
1798    ws[0] = 1.;          ws[1] = 3.;
1799    density      = 3.034;
1800    id           = 7;
1801    AliMixture(id, "V0APlaOu", as, zs, density, -2, ws);
1802    AliMedium(id, "V0APlaOu", id, 1, fieldType, maxField, maxBending, maxStepSize,
1803              maxEnergyLoss, precision, minStepSize);
1804
1805 // Parameters for V0Aplate: EuroComposite - EC-PI 626 PS - EC-PI 6.4-42
1806    as[0] = 1.00794;     as[1] = 12.011;
1807    zs[0] = 1.;          zs[1] = 6.;
1808    ws[0] = 5.21;        ws[1] = 4.74;
1809    density      = 0.042;
1810    id           = 8;
1811    AliMixture(id, "V0APlaIn", as, zs, density, -2, ws);
1812    AliMedium(id, "V0APlaIn", id, 1, fieldType, maxField, maxBending, maxStepSize,
1813              maxEnergyLoss, precision, minStepSize);
1814
1815 // Parameters for V0Afiber: BC9929AMC Plastic Scintillating Fiber from Saint-Gobain
1816    as[0] = 1.00794;     as[1] = 12.011;
1817    zs[0] = 1.;          zs[1] = 6.;
1818    ws[0] = 4.82;        ws[1] = 4.85;
1819    density      = 1.05;
1820    id           = 9;
1821    AliMixture(id, "V0AFib", as, zs, density, -2, ws);
1822    AliMedium(id, "V0AFib", id, 1, fieldType, maxField, maxBending, maxStepSize,
1823              maxEnergyLoss, precision, minStepSize);
1824
1825 // Parameters for V0APMA: Aluminium
1826    a = 26.98;
1827    z = 13.00;
1828    density     = 2.7;
1829    radLength   = 8.9;
1830    absLength   = 37.2;
1831    id = 10;
1832    AliMaterial(id, "V0APMA",  a, z, density, radLength, absLength, 0, 0);
1833    AliMedium(id, "V0APMA", id, 1, fieldType, maxField, maxBending, maxStepSize,
1834              maxEnergyLoss, precision, minStepSize);
1835
1836 // Parameters for V0APMG: Glass for the simulation Aluminium
1837    a = 26.98;
1838    z = 13.00;
1839    density   = 2.7;
1840    radLength = 8.9;
1841    absLength = 37.2;
1842    id = 11;
1843    AliMaterial(id, "V0APMG",  a, z, density, radLength, absLength, 0, 0);
1844    AliMedium(id, "V0APMG", id, 1, fieldType, maxField, maxBending, maxStepSize,
1845              maxEnergyLoss, precision, minStepSize);
1846 }
1847
1848 //_____________________________________________________________________________
1849 void AliVZEROv7::DrawModule() const
1850 {
1851 //  Drawing is done in DrawVZERO.C
1852
1853    AliDebug(2,"DrawModule");
1854 }
1855
1856
1857 //_____________________________________________________________________________
1858 void AliVZEROv7::DrawGeometry()
1859 {
1860 //  Drawing of V0 geometry done in DrawV0.C
1861
1862    AliDebug(2,"DrawGeometry");
1863 }
1864
1865 //_____________________________________________________________________________
1866 void AliVZEROv7::Init()
1867 {
1868 // Initialises version of the VZERO Detector given in Config
1869 // Just prints an information message
1870
1871 //   AliInfo(Form("VZERO version %d initialized \n",IsVersion()));
1872
1873    AliDebug(1,"VZERO version 7 initialized");
1874    AliVZERO::Init();
1875 }
1876
1877 //_____________________________________________________________________________
1878 void AliVZEROv7::StepManager()
1879 {
1880 // Step Manager, called at each step
1881
1882   Int_t     copy;
1883   static    Int_t   vol[4];
1884   static    Float_t hits[21];
1885   static    Float_t eloss, tlength;
1886   static    Int_t   nPhotonsInStep = 0;
1887   static    Int_t   nPhotons = 0;
1888   static    Int_t   numStep = 0;
1889   Int_t     ringNumber;
1890   Float_t   destep, step;
1891   numStep += 1;
1892
1893   //   We keep only charged tracks :
1894   if ( !gMC->TrackCharge() || !gMC->IsTrackAlive() ) return;
1895
1896   vol[0]    = gMC->CurrentVolOffID(1, vol[1]);
1897   vol[2]    = gMC->CurrentVolID(copy);
1898   vol[3]    = copy;
1899   static Int_t idV0R1 = gMC->VolId("V0R1");
1900   static Int_t idV0L1 = gMC->VolId("V0L1");
1901   static Int_t idV0L15 = gMC->VolId("V0L15");
1902   static Int_t idV0L16 = gMC->VolId("V0L16");
1903   static Int_t idV0L17 = gMC->VolId("V0L17");
1904   static Int_t idV0L18 = gMC->VolId("V0L18");
1905   static Int_t idV0R2 = gMC->VolId("V0R2");
1906   static Int_t idV0L2 = gMC->VolId("V0L2");
1907   static Int_t idV0L25 = gMC->VolId("V0L25");
1908   static Int_t idV0L26 = gMC->VolId("V0L26");
1909   static Int_t idV0L27 = gMC->VolId("V0L27");
1910   static Int_t idV0L28 = gMC->VolId("V0L28");
1911   static Int_t idV0R3 = gMC->VolId("V0R3");
1912   static Int_t idV0L3 = gMC->VolId("V0L3");
1913   static Int_t idV0L35 = gMC->VolId("V0L35");
1914   static Int_t idV0L36 = gMC->VolId("V0L36");
1915   static Int_t idV0L37 = gMC->VolId("V0L37");
1916   static Int_t idV0L38 = gMC->VolId("V0L38");
1917   static Int_t idV0R4 = gMC->VolId("V0R4");
1918   static Int_t idV0L4 = gMC->VolId("V0L4");
1919   static Int_t idV0L45 = gMC->VolId("V0L45");
1920   static Int_t idV0L46 = gMC->VolId("V0L46");
1921   static Int_t idV0L47 = gMC->VolId("V0L47");
1922   static Int_t idV0L48 = gMC->VolId("V0L48");
1923   static Int_t idV0R5 = gMC->VolId("V0R5");
1924   static Int_t idV0R6 = gMC->VolId("V0R6");
1925   bool   hitOnV0C = true;
1926   double lightYield;
1927   double lightAttenuation;
1928   double nMeters;
1929   double fibToPhot;
1930   if      ( gMC->CurrentVolID(copy) == idV0R1 || gMC->CurrentVolID(copy) == idV0L1 || gMC->CurrentVolID(copy) == idV0L15 ||  gMC->CurrentVolID(copy) == idV0L16  || gMC->CurrentVolID(copy) == idV0L17  || gMC->CurrentVolID(copy) == idV0L18  )
1931     ringNumber = 1;
1932   else if ( gMC->CurrentVolID(copy) == idV0R2 || gMC->CurrentVolID(copy) == idV0L2 || gMC->CurrentVolID(copy) == idV0L25 || gMC->CurrentVolID(copy) == idV0L26 || gMC->CurrentVolID(copy) == idV0L27 || gMC->CurrentVolID(copy) == idV0L28 )
1933     ringNumber = 2;
1934   else if ( gMC->CurrentVolID(copy) == idV0R3 || gMC->CurrentVolID(copy) == idV0R4
1935             || gMC->CurrentVolID(copy) == idV0L3 || gMC->CurrentVolID(copy) == idV0L35 || gMC->CurrentVolID(copy) == idV0L36 || gMC->CurrentVolID(copy) == idV0L37 || gMC->CurrentVolID(copy) == idV0L38 ) ringNumber = 3;
1936   else if ( gMC->CurrentVolID(copy) == idV0R5 || gMC->CurrentVolID(copy) == idV0R6
1937             || gMC->CurrentVolID(copy) == idV0L4 || gMC->CurrentVolID(copy) == idV0L45 || gMC->CurrentVolID(copy) == idV0L46 || gMC->CurrentVolID(copy) == idV0L47 || gMC->CurrentVolID(copy) == idV0L48 ) ringNumber = 4;
1938   else ringNumber = 0;
1939   if  (ringNumber) {
1940     if (gMC->CurrentVolID(copy) == idV0L1 || gMC->CurrentVolID(copy) == idV0L15  || gMC->CurrentVolID(copy) == idV0L16  || gMC->CurrentVolID(copy) == idV0L17 || gMC->CurrentVolID(copy) == idV0L18 || gMC->CurrentVolID(copy) == idV0L2 || gMC->CurrentVolID(copy) == idV0L25 || gMC->CurrentVolID(copy) == idV0L26  || gMC->CurrentVolID(copy) == idV0L27  || gMC->CurrentVolID(copy) == idV0L28 || gMC->CurrentVolID(copy) == idV0L3 || gMC->CurrentVolID(copy) == idV0L35 || gMC->CurrentVolID(copy) == idV0L36 || gMC->CurrentVolID(copy) == idV0L37 || gMC->CurrentVolID(copy) == idV0L38 || gMC->CurrentVolID(copy) == idV0L4 || gMC->CurrentVolID(copy) == idV0L45 || gMC->CurrentVolID(copy) == idV0L46 || gMC->CurrentVolID(copy) == idV0L47 || gMC->CurrentVolID(copy) == idV0L48)
1941       hitOnV0C = false;
1942     destep = gMC->Edep();
1943     step   = gMC->TrackStep();
1944     if (hitOnV0C) {
1945       lightYield = fV0CLightYield;
1946       lightAttenuation = fV0CLightAttenuation;
1947       nMeters = fV0CnMeters;
1948       fibToPhot = fV0CFibToPhot;
1949     } else {
1950       lightYield = fV0ALightYield;
1951       lightAttenuation = fV0ALightAttenuation;
1952       nMeters = fV0AnMeters;
1953       fibToPhot = fV0AFibToPhot;
1954     }
1955     nPhotonsInStep  = Int_t(destep / (lightYield *1e-9) );
1956     nPhotonsInStep  = gRandom->Poisson(nPhotonsInStep);
1957     eloss    += destep;
1958     tlength  += step;
1959     if ( gMC->IsTrackEntering() ) {
1960       nPhotons  =  nPhotonsInStep;
1961       gMC->TrackPosition(fTrackPosition);
1962       gMC->TrackMomentum(fTrackMomentum);
1963       Float_t pt  = TMath::Sqrt( fTrackMomentum.Px() * fTrackMomentum.Px()
1964                                  + fTrackMomentum.Py() * fTrackMomentum.Py() );
1965       TParticle *par = gAlice->GetMCApp()->Particle(gAlice->GetMCApp()->GetCurrentTrackNumber());
1966       hits[0]  = fTrackPosition.X();
1967       hits[1]  = fTrackPosition.Y();
1968       hits[2]  = fTrackPosition.Z();
1969       hits[3]  = Float_t (gMC->TrackPid());
1970       hits[4]  = gMC->TrackTime();
1971       hits[5]  = gMC->TrackCharge();
1972       hits[6]  = fTrackMomentum.Theta()*TMath::RadToDeg();
1973       hits[7]  = fTrackMomentum.Phi()*TMath::RadToDeg();
1974       hits[8]  = ringNumber;
1975       hits[9]  = pt;
1976       hits[10] = fTrackMomentum.P();
1977       hits[11] = fTrackMomentum.Px();
1978       hits[12] = fTrackMomentum.Py();
1979       hits[13] = fTrackMomentum.Pz();
1980       hits[14] = par->Vx();
1981       hits[15] = par->Vy();
1982       hits[16] = par->Vz();
1983       tlength  = 0.0;
1984       eloss    = 0.0;
1985
1986       //////////////////////////
1987       ///// Display V0A geometry
1988       //      if (!hitOnV0C) {
1989       //        FILE *of;
1990       //        of = fopen("V0A.out", "a");
1991       //        // x, y, z, ringnumber, cellid
1992       //        fprintf( of, "%f %f %f %f %d \n",  hits[0], hits[1], hits[2], hits[8], GetCellId (vol, hits) );
1993       //        fclose(of);
1994       //      }
1995       //////////////////////////
1996     }
1997     nPhotons  = nPhotons + nPhotonsInStep;
1998     if( gMC->IsTrackExiting() || gMC->IsTrackStop() || gMC->IsTrackDisappeared()){
1999       nPhotons = nPhotons - Int_t((Float_t(nPhotons) * lightAttenuation * nMeters));
2000       nPhotons = nPhotons - Int_t( Float_t(nPhotons) * fibToPhot);
2001       hits[17] = eloss;
2002       hits[18] = tlength;
2003       hits[19] = nPhotons;
2004       hits[20] = GetCellId (vol, hits);
2005       AddHit(gAlice->GetMCApp()->GetCurrentTrackNumber(), vol, hits);
2006       tlength         = 0.0;
2007       eloss           = 0.0;
2008       nPhotons        = 0;
2009       nPhotonsInStep  = 0;
2010       numStep         = 0;
2011     }
2012     if( gMC->IsTrackEntering() || gMC->IsTrackExiting() ) {
2013       AddTrackReference(gAlice->GetMCApp()->GetCurrentTrackNumber(), AliTrackReference::kVZERO);
2014     }
2015   }
2016 }
2017
2018 //_____________________________________________________________________________
2019 void AliVZEROv7::AddHit(Int_t track, Int_t *vol, Float_t *hits)
2020 {
2021 //  Adds a VZERO hit
2022
2023   TClonesArray &lhits = *fHits;
2024   new(lhits[fNhits++]) AliVZEROhit(fIshunt,track,vol,hits);
2025 }
2026
2027 //_____________________________________________________________________________
2028 void AliVZEROv7::AddDigits(Int_t *tracks, Int_t* digits)
2029 {
2030 //  Adds a VZERO digit
2031
2032    TClonesArray  &ldigits = *fDigits;
2033    new(ldigits[fNdigits++]) AliVZEROdigit(tracks, digits);
2034 }
2035
2036 //_____________________________________________________________________________
2037 void AliVZEROv7::MakeBranch(Option_t *option)
2038 {
2039 // Creates new branches in the current Root Tree
2040
2041   char branchname[10];
2042   sprintf(branchname,"%s",GetName());
2043   AliDebug(2,Form("fBufferSize = %d",fBufferSize));
2044   const char *cH = strstr(option,"H");
2045   if (fHits   && fLoader->TreeH() && cH) {
2046     fLoader->TreeH()->Branch(branchname,&fHits, fBufferSize);
2047     AliDebug(2,Form("Making Branch %s for hits",branchname));
2048   }
2049   const char *cD = strstr(option,"D");
2050   if (fDigits   && fLoader->TreeD() && cD) {
2051     fLoader->TreeD()->Branch(branchname,&fDigits, fBufferSize);
2052     AliDebug(2,Form("Making Branch %s for digits",branchname));
2053   }
2054 }
2055
2056 //_____________________________________________________________________________
2057 Int_t AliVZEROv7::GetCellId(Int_t *vol, Float_t *hits)
2058 {
2059   //   Returns Id of scintillator cell
2060   //   Right side from  0 to 47
2061   //   Left  side from 48 to 79
2062   //   hits[8] = ring number (1 to 4)
2063   //   vol[1]  = copy number (1 to 8)
2064
2065   Int_t index      = vol[1];
2066   Int_t ringNumber = Int_t(hits[8]);
2067   fCellId          = 0;
2068
2069   Float_t phi = Float_t(TMath::ATan2(Double_t(hits[1]),Double_t(hits[0])) );
2070   Float_t kRaddeg = 180.0/TMath::Pi();
2071   phi = kRaddeg * phi;
2072
2073   if (index < 7) index = index + 8;
2074
2075   if (hits[2] < 0.0) {
2076     if(ringNumber < 3) {
2077       index = (index - 7) + ( ( ringNumber - 1 ) * 8);
2078     } else if (ringNumber >= 3) {
2079       if ( gMC->CurrentVolID(vol[1]) == gMC->VolId("V0R3") || gMC->CurrentVolID(vol[1])
2080            == gMC->VolId("V0R5") )  index = (index*2-14)+((ringNumber-2)*16);
2081       if ( gMC->CurrentVolID(vol[1]) == gMC->VolId("V0R4") || gMC->CurrentVolID(vol[1])
2082            == gMC->VolId("V0R6") )  index = (index*2-13)+((ringNumber-2)*16);
2083     }
2084     fCellId   = index;
2085   } else if (hits[2] > 0.0) {
2086     //    cout << " vol[0] = " << vol[0] << " copy : " << vol[1]
2087     //   << " called " << gMC->VolName(vol[0]) << endl;
2088     // cout << " vol[2] = " << vol[2] << " copy : " << vol[3]
2089     //   << " called " << gMC->VolName(vol[2]) << endl;
2090     if (gMC->CurrentVolID(vol[2]) == gMC->VolId("V0L1")) fCellId = vol[1]+47;
2091     if (gMC->CurrentVolID(vol[2]) == gMC->VolId("V0L2")) fCellId = 8+vol[1]+47;
2092     if (gMC->CurrentVolID(vol[2]) == gMC->VolId("V0L3")) fCellId = 16+vol[1]+47;
2093     if (gMC->CurrentVolID(vol[2]) == gMC->VolId("V0L4")) fCellId = 24+vol[1]+47;
2094     if (gMC->CurrentVolID(vol[2]) == gMC->VolId("V0L15")) fCellId = 48+4;
2095     if (gMC->CurrentVolID(vol[2]) == gMC->VolId("V0L16")) fCellId = 48+5;
2096     if (gMC->CurrentVolID(vol[2]) == gMC->VolId("V0L17")) fCellId = 48+6;
2097     if (gMC->CurrentVolID(vol[2]) == gMC->VolId("V0L18")) fCellId = 48+7;
2098     if (gMC->CurrentVolID(vol[2]) == gMC->VolId("V0L25")) fCellId = 8+48+4;
2099     if (gMC->CurrentVolID(vol[2]) == gMC->VolId("V0L26")) fCellId = 8+48+5;
2100     if (gMC->CurrentVolID(vol[2]) == gMC->VolId("V0L27")) fCellId = 8+48+6;
2101     if (gMC->CurrentVolID(vol[2]) == gMC->VolId("V0L28")) fCellId = 8+48+7;
2102     if (gMC->CurrentVolID(vol[2]) == gMC->VolId("V0L35")) fCellId = 16+48+4;
2103     if (gMC->CurrentVolID(vol[2]) == gMC->VolId("V0L36")) fCellId = 16+48+5;
2104     if (gMC->CurrentVolID(vol[2]) == gMC->VolId("V0L37")) fCellId = 16+48+6;
2105     if (gMC->CurrentVolID(vol[2]) == gMC->VolId("V0L38")) fCellId = 16+48+7;
2106     if (gMC->CurrentVolID(vol[2]) == gMC->VolId("V0L45")) fCellId = 24+48+4;
2107     if (gMC->CurrentVolID(vol[2]) == gMC->VolId("V0L46")) fCellId = 24+48+5;
2108     if (gMC->CurrentVolID(vol[2]) == gMC->VolId("V0L47")) fCellId = 24+48+6;
2109     if (gMC->CurrentVolID(vol[2]) == gMC->VolId("V0L48")) fCellId = 24+48+7;
2110   }
2111
2112   return fCellId;
2113 }
2114