Shil.cxx

Go to the documentation of this file.

// Copyright 2019-2020 CERN and copyright holders of ALICE O2.
// See https://alice-o2.web.cern.ch/copyright for details of the copyright holders.
// All rights not expressly granted are reserved.
//
// This software is distributed under the terms of the GNU General Public
// License v3 (GPL Version 3), copied verbatim in the file "COPYING".
//
// In applying this license CERN does not waive the privileges and immunities
// granted to it by virtue of its status as an Intergovernmental Organization
// or submit itself to any jurisdiction.
 
#include <DetectorsBase/Detector.h>
#include <DetectorsPassive/Shil.h>
#include <DetectorsBase/MaterialManager.h>
#include <TGeoCompositeShape.h>
#include <TGeoCone.h>
#include <TGeoManager.h>
#include <TGeoMatrix.h>
#include <TGeoPcon.h>
#include <TGeoTube.h>
#include <TGeoVolume.h>
#ifdef NDEBUG
#undef NDEBUG
#endif
#include <cassert>
 
using namespace o2::passive;
 
#define kDegrad TMath::DegToRad()
 
Shil::~Shil() = default;
 
Shil::Shil() : PassiveBase("SHIL", "") {}
Shil::Shil(const char* name, const char* Title) : PassiveBase(name, Title) {}
Shil::Shil(const Shil& rhs) = default;
 
Shil& Shil::operator=(const Shil& rhs)
{
  // self assignment
  if (this == &rhs) {
    return *this;
  }
 
  // base class assignment
  PassiveBase::operator=(rhs);
 
  return *this;
}
 
void InvertPcon(TGeoPcon* pcon)
{
  //
  //  z -> -z
  //
  Int_t nz = pcon->GetNz();
  Double_t* z = new Double_t[nz];
  Double_t* rmin = new Double_t[nz];
  Double_t* rmax = new Double_t[nz];
 
  Double_t* z0 = pcon->GetZ();
  Double_t* rmin0 = pcon->GetRmin();
  Double_t* rmax0 = pcon->GetRmax();
 
  for (Int_t i = 0; i < nz; i++) {
    z[i] = z0[i];
    rmin[i] = rmin0[i];
    rmax[i] = rmax0[i];
  }
 
  for (Int_t i = 0; i < nz; i++) {
    Int_t j = nz - i - 1;
    pcon->DefineSection(i, -z[j], rmin[j], rmax[j]);
  }
 
  delete[] z;
  delete[] rmin;
  delete[] rmax;
}
 
namespace
{
TGeoPcon* MakeShapeFromTemplate(const TGeoPcon* pcon, Float_t drMin, Float_t drMax)
{
  //
  // Returns new shape based on a template changing
  // the inner radii by drMin and the outer radii by drMax.
  //
  Int_t nz = pcon->GetNz();
  TGeoPcon* cpcon = new TGeoPcon(0., 360., nz);
  for (Int_t i = 0; i < nz; i++) {
    cpcon->DefineSection(i, pcon->GetZ(i), pcon->GetRmin(i) + drMin, pcon->GetRmax(i) + drMax);
  }
  return cpcon;
}
} // namespace
 
void Shil::ConstructGeometry()
{
  createMaterials();
 
  //
  // The geometry of the small angle absorber "Beam Shield"
  //
 
  //
  // The top volume
  //
  TGeoVolume* top = gGeoManager->GetVolume("cave");
 
  Float_t dz, dr, z, rmax;
 
  //  Rotations
  TGeoRotation* rot000 = new TGeoRotation("rot000", 90., 0., 90., 90., 0., 0.);
  TGeoRotation* rot090 = new TGeoRotation("rot090", 90., 90., 90., 180., 0., 0.);
  TGeoRotation* rot180 = new TGeoRotation("rot180", 90., 180., 90., 270., 0., 0.);
  TGeoRotation* rot270 = new TGeoRotation("rot270", 90., 270., 90., 0., 0., 0.);
  Float_t alhc = 0.794;
  TGeoRotation* rotxzlhc = new TGeoRotation("rotxzlhc", 0., -alhc, 0.);
  TGeoRotation* rotlhc = new TGeoRotation("rotlhc", 0., alhc, 0.);
 
  //
  // Media
  //
  auto& matmgr = o2::base::MaterialManager::Instance();
  auto kMedNiW = matmgr.getTGeoMedium("SHIL_Ni/W0");
  auto kMedNiWsh = matmgr.getTGeoMedium("SHIL_Ni/W3");
  //
  auto kMedSteel = matmgr.getTGeoMedium("SHIL_ST_C0");
  auto kMedSteelSh = matmgr.getTGeoMedium("SHIL_ST_C3");
  //
  auto kMedAir = matmgr.getTGeoMedium("SHIL_AIR_C0");
  auto kMedAirMu = matmgr.getTGeoMedium("SHIL_AIR_MUON");
  //
  auto kMedPb = matmgr.getTGeoMedium("SHIL_PB_C0");
  auto kMedPbSh = matmgr.getTGeoMedium("SHIL_PB_C2");
  //
  auto kMedConcSh = matmgr.getTGeoMedium("SHIL_CC_C2");
  //
  auto kMedCastiron = matmgr.getTGeoMedium("SHIL_CAST_IRON0");
  auto kMedCastironSh = matmgr.getTGeoMedium("SHIL_CAST_IRON2");
  //
  const Float_t kDegRad = TMath::Pi() / 180.;
  const Float_t kAngle02 = TMath::Tan(2.00 * kDegRad);
  const Float_t kAngle0071 = TMath::Tan(0.71 * kDegRad);
 
  //    FA Tungsten Tail           //
  //    Drawing ALIP2A__0049       //
  //    Drawing ALIP2A__0111       //
  //
  //    The tail as built is shorter than in drawing ALIP2A__0049.
  //    The CDD data base has to be updated !
  //
  //    Inner radius at the entrance of the flange
  Float_t rInFaWTail1 = 13.98 / 2.;
  //    Outer radius at the entrance of the flange
  Float_t rOuFaWTail1 = 52.00 / 2.;
  //    Outer radius at the end of the section inside the FA
  Float_t rOuFaWTail2 = 35.27 / 2.;
  //    Length of the Flange section inside the FA
  Float_t dzFaWTail1 = 6.00;
  //    Length of the Flange section ouside the FA
  Float_t dzFaWTail2 = 12.70;
  //    Inner radius at the end of the section inside the FA
  Float_t rInFaWTail2 = rInFaWTail1 + dzFaWTail1 * kAngle0071;
  //    Inner radius at the end of the flange
  Float_t rInFaWTail3 = rInFaWTail2 + dzFaWTail2 * kAngle0071;
  //    Outer radius at the end of the flange
  Float_t rOuFaWTail3 = rOuFaWTail2 + dzFaWTail2 * kAngle02;
  //    Outer radius of the recess for station 1
  Float_t rOuFaWTailR = 30.8 / 2.;
  //    Length of the recess
  Float_t dzFaWTailR = 36.00;
  //    Inner radiues at the end of the recess
  Float_t rInFaWTail4 = rInFaWTail3 + dzFaWTailR * kAngle0071;
  //    Outer radius at the end of the recess
  Float_t rOuFaWTail4 = rOuFaWTail3 + dzFaWTailR * kAngle02;
  //    Inner radius of the straight section
  Float_t rInFaWTailS = 22.30 / 2.;
  //    Length of the bulge
  Float_t dzFaWTailB = 13.0;
  //    Outer radius at the end of the bulge
  Float_t rOuFaWTailB = rOuFaWTail4 + dzFaWTailB * kAngle02;
  //    Outer radius at the end of the tail
  Float_t rOuFaWTailE = 31.6 / 2.;
  //    Total length of the tail
  const Float_t dzFaWTail = 70.7 - 6.;
 
  TGeoPcon* shFaWTail = new TGeoPcon(0., 360., 8);
  z = 0.;
  //    Flange section inside FA (moved to ABSO)
  //  shFaWTail->DefineSection(0, z, rInFaWTail1, rOuFaWTail1);
  //  z += dzFaWTail1;
  //  shFaWTail->DefineSection(1, z, rInFaWTail2, rOuFaWTail1);
  shFaWTail->DefineSection(0, z, rInFaWTail2, rOuFaWTail2);
  //    Flange section outside FA
  z += dzFaWTail2;
  shFaWTail->DefineSection(1, z, rInFaWTail3, rOuFaWTail3);
  shFaWTail->DefineSection(2, z, rInFaWTail3, rOuFaWTailR);
  //    Recess Station 1
  z += dzFaWTailR;
  shFaWTail->DefineSection(3, z, rInFaWTail4, rOuFaWTailR);
  shFaWTail->DefineSection(4, z, rInFaWTailS, rOuFaWTail4);
  //    Bulge
  z += dzFaWTailB;
  shFaWTail->DefineSection(5, z, rInFaWTailS, rOuFaWTailB);
  shFaWTail->DefineSection(6, z, rInFaWTailS, rOuFaWTailE);
  //    End
  z = dzFaWTail;
  shFaWTail->DefineSection(7, z, rInFaWTailS, rOuFaWTailE);
 
  TGeoVolume* voFaWTail = new TGeoVolume("YFaWTail", shFaWTail, kMedNiW);
  //
  //    Define an inner region with higher transport cuts
  TGeoPcon* shFaWTailI = new TGeoPcon(0., 360., 4);
  z = 0.;
  dr = 3.5;
  shFaWTailI->DefineSection(0, z, rInFaWTail2, rInFaWTail2 + dr);
  z += (dzFaWTail2 + dzFaWTailR);
  shFaWTailI->DefineSection(1, z, rInFaWTail4, rInFaWTail4 + dr);
  shFaWTailI->DefineSection(2, z, rInFaWTailS, rInFaWTailS + dr);
  z = dzFaWTail;
  shFaWTailI->DefineSection(3, z, rInFaWTailS, rInFaWTailS + dr);
  TGeoVolume* voFaWTailI = new TGeoVolume("YFaWTailI", shFaWTailI, kMedNiWsh);
  voFaWTail->AddNode(voFaWTailI, 1, gGeoIdentity);
 
  //                               //
  // Recess Station 1              //
  // Drawing ALIP2A__0260          //
 
  //    FA W-Ring 2                //
  //    Drawing ALIP2A__0220       //
  const Float_t kFaWring2Rinner = 15.41;
  const Float_t kFaWring2Router = 18.40;
  const Float_t kFaWring2HWidth = 3.75;
  const Float_t kFaWring2Cutoffx = 3.35;
  const Float_t kFaWring2Cutoffy = 3.35;
  TGeoTubeSeg* shFaWring2a = new TGeoTubeSeg(kFaWring2Rinner, kFaWring2Router, kFaWring2HWidth, 0., 90.);
  shFaWring2a->SetName("shFaWring2a");
  TGeoBBox* shFaWring2b = new TGeoBBox(kFaWring2Router / 2., kFaWring2Router / 2., kFaWring2HWidth);
  shFaWring2b->SetName("shFaWring2b");
  TGeoTranslation* trFaWring2b = new TGeoTranslation("trFaWring2b", kFaWring2Router / 2. + kFaWring2Cutoffx,
                                                     kFaWring2Router / 2. + kFaWring2Cutoffy, 0.);
  trFaWring2b->RegisterYourself();
  TGeoCompositeShape* shFaWring2 = new TGeoCompositeShape("shFaWring2", "(shFaWring2a)*(shFaWring2b:trFaWring2b)");
  TGeoVolume* voFaWring2 = new TGeoVolume("YFA_WRING2", shFaWring2, kMedNiW);
 
  //    FA W-Ring 3                //
  //    Drawing ALIP2A__0219       //
  const Float_t kFaWring3Rinner = 15.41;
  const Float_t kFaWring3Router = 18.40;
  const Float_t kFaWring3HWidth = 3.75;
  const Float_t kFaWring3Cutoffx = 3.35;
  const Float_t kFaWring3Cutoffy = 3.35;
  TGeoTubeSeg* shFaWring3a = new TGeoTubeSeg(kFaWring3Rinner, kFaWring3Router, kFaWring3HWidth, 0., 90.);
  shFaWring3a->SetName("shFaWring3a");
  TGeoBBox* shFaWring3b = new TGeoBBox(kFaWring3Router / 2., kFaWring3Router / 2., kFaWring3HWidth);
  shFaWring3b->SetName("shFaWring3b");
  TGeoTranslation* trFaWring3b = new TGeoTranslation("trFaWring3b", kFaWring3Router / 2. + kFaWring3Cutoffx,
                                                     kFaWring3Router / 2. + kFaWring3Cutoffy, 0.);
  trFaWring3b->RegisterYourself();
  TGeoCompositeShape* shFaWring3 = new TGeoCompositeShape("shFaWring3", "(shFaWring3a)*(shFaWring3b:trFaWring3b)");
  TGeoVolume* voFaWring3 = new TGeoVolume("YFA_WRING3", shFaWring3, kMedNiW);
 
  //    FA W-Ring 5                //
  //    Drawing ALIP2A__0221       //
  const Float_t kFaWring5Rinner = 15.41;
  const Float_t kFaWring5Router = 18.67;
  const Float_t kFaWring5HWidth = 1.08;
  TGeoVolume* voFaWring5 =
    new TGeoVolume("YFA_WRING5", new TGeoTube(kFaWring5Rinner, kFaWring5Router, kFaWring5HWidth), kMedNiW);
 
  //
  // Position the rings in the assembly
  //
  TGeoVolumeAssembly* asFaExtraShield = new TGeoVolumeAssembly("YCRE");
  // Distance between rings
  const Float_t kFaDWrings = 1.92;
  //
  dz = 0.;
 
  dz += kFaWring2HWidth;
  asFaExtraShield->AddNode(voFaWring2, 1, new TGeoCombiTrans(0., 0., dz, rot180));
  asFaExtraShield->AddNode(voFaWring2, 2, new TGeoCombiTrans(0., 0., dz, rot000));
  dz += kFaWring2HWidth;
  dz += kFaDWrings;
  dz += kFaWring3HWidth;
  asFaExtraShield->AddNode(voFaWring3, 1, new TGeoCombiTrans(0., 0., dz, rot090));
  asFaExtraShield->AddNode(voFaWring3, 2, new TGeoCombiTrans(0., 0., dz, rot270));
  dz += kFaWring3HWidth;
  dz += kFaWring5HWidth;
  asFaExtraShield->AddNode(voFaWring5, 1, new TGeoTranslation(0., 0., dz));
  dz += kFaWring5HWidth;
  dz += kFaWring3HWidth;
  asFaExtraShield->AddNode(voFaWring3, 3, new TGeoCombiTrans(0., 0., dz, rot180));
  asFaExtraShield->AddNode(voFaWring3, 4, new TGeoCombiTrans(0., 0., dz, rot000));
  dz += kFaWring3HWidth;
  dz += kFaDWrings;
  dz += kFaWring2HWidth;
  asFaExtraShield->AddNode(voFaWring2, 3, new TGeoCombiTrans(0., 0., dz, rot090));
  asFaExtraShield->AddNode(voFaWring2, 4, new TGeoCombiTrans(0., 0., dz, rot270));
  dz += kFaWring2HWidth;
 
  //                SAA1               //
 
  //          FA/SAA1  W Joint         //
  //          Drawing ALIP2A__0060     //
 
  // Length of flange FA side
  Float_t dzFaSaa1F1 = 2.8;
  // Inner radius of flange FA side
  Float_t rInFaSaa1F1 = 32.0 / 2.;
  // Outer radius of flange FA side
  Float_t rOuFaSaa1F1 = 39.5 / 2.;
  // Length of first straight section
  Float_t dzFaSaa1S1 = 18.5 - dzFaSaa1F1;
  // Inner radius of first straight section
  Float_t rInFaSaa1S1 = 22.3 / 2.;
  // Length of 45 deg transition region
  Float_t dzFaSaa1T1 = 2.2;
  // Inner radius of second straight section
  Float_t rInFaSaa1S2 = 17.9 / 2.;
  // Length of second straight section
  Float_t dzFaSaa1S2 = 10.1;
  // Length of flange SAA1 side
  //    Float_t dzFaSaa1F2  =  4.0;
  // Inner radius of flange FA side
  Float_t rInFaSaa1F2 = 25.2 / 2.;
  // Length of joint
  const Float_t dzFaSaa1 = 34.8;
  // Outer Radius at the end of the joint
  Float_t rOuFaSaa1E = 41.93 / 2.;
 
  TGeoPcon* shFaSaa1 = new TGeoPcon(0., 360., 8);
  z = 0;
  // Flange FA side
  shFaSaa1->DefineSection(0, z, rInFaSaa1F1, rOuFaSaa1F1 - 0.01);
  z += dzFaSaa1F1;
  shFaSaa1->DefineSection(1, z, rInFaSaa1F1, 40.0);
  shFaSaa1->DefineSection(2, z, rInFaSaa1S1, 40.0);
  // First straight section
  z += dzFaSaa1S1;
  shFaSaa1->DefineSection(3, z, rInFaSaa1S1, 40.0);
  // 45 deg transition region
  z += dzFaSaa1T1;
  shFaSaa1->DefineSection(4, z, rInFaSaa1S2, 40.0);
  // Second straight section
  z += dzFaSaa1S2;
  shFaSaa1->DefineSection(5, z, rInFaSaa1S2, 40.0);
  shFaSaa1->DefineSection(6, z, rInFaSaa1F2, 40.0);
  // Flange SAA1 side
  z = dzFaSaa1;
  shFaSaa1->DefineSection(7, z, rInFaSaa1F2, rOuFaSaa1E - 0.01);
 
  // Outer 2 deg line
  for (Int_t i = 1; i < 7; i++) {
    Double_t zp = shFaSaa1->GetZ(i);
    Double_t r1 = shFaSaa1->GetRmin(i);
    Double_t r2 = 39.5 / 2. + zp * TMath::Tan(2. * kDegRad) - 0.01;
    shFaSaa1->DefineSection(i, zp, r1, r2);
  }
  TGeoVolume* voFaSaa1 = new TGeoVolume("YFASAA1", shFaSaa1, kMedNiWsh);
  //
  // Outer region with lower transport cuts
  TGeoCone* shFaSaa1O =
    new TGeoCone(dzFaSaa1 / 2., rOuFaSaa1F1 - 3.5, rOuFaSaa1F1 - 0.01, rOuFaSaa1E - 3.5, rOuFaSaa1E - 0.01);
  TGeoVolume* voFaSaa1O = new TGeoVolume("YFASAA1O", shFaSaa1O, kMedNiW);
  voFaSaa1->AddNode(voFaSaa1O, 1, new TGeoTranslation(0., 0., dzFaSaa1 / 2.));
 
  //    SAA1 Steel Envelope        //
  //    Drawing ALIP2A__0039       //
 
  Float_t rOut; // Outer radius
                // Thickness of the steel envelope
  Float_t dSt = 4.;
  // 4 Section
  // z-positions
  Float_t zSaa1StEnv[5] = {111.2, 113.7, 229.3, 195.0};
  // Radii
  // 1
  Float_t rOuSaa1StEnv1 = 40.4 / 2.;
  Float_t rInSaa1StEnv1 = rOuSaa1StEnv1 - dSt - 0.05;
  // 2
  Float_t rInSaa1StEnv2 = 41.7 / 2.;
  Float_t rOuSaa1StEnv2 = rInSaa1StEnv2 + dSt / TMath::Cos(2.0 * kDegRad) - 0.05;
  // 3
  Float_t rOuSaa1StEnv3 = 57.6 / 2.;
  Float_t rInSaa1StEnv3 = rOuSaa1StEnv3 - dSt + 0.05;
  // 4
  Float_t rInSaa1StEnv4 = 63.4 / 2.;
  Float_t rOuSaa1StEnv4 = rInSaa1StEnv4 + dSt / TMath::Cos(1.6 * kDegRad) - 0.05;
  // end
  Float_t rInSaa1StEnv5 = 74.28 / 2.;
  Float_t rOuSaa1StEnv5 = rInSaa1StEnv5 + dSt / TMath::Cos(1.6 * kDegRad) - 0.05;
  // Relative starting position
  Float_t zSaa1StEnvS = 3.;
 
  TGeoPcon* shSaa1StEnv = new TGeoPcon(0., 360., 11);
  // 1st Section
  z = zSaa1StEnvS;
  shSaa1StEnv->DefineSection(0, z, rInSaa1StEnv1, rOuSaa1StEnv1);
  z += (zSaa1StEnv[0] - dSt);
  shSaa1StEnv->DefineSection(1, z, rInSaa1StEnv1, rOuSaa1StEnv1);
  // 1 - 2
  shSaa1StEnv->DefineSection(2, z, rInSaa1StEnv1, rOuSaa1StEnv2);
  z += dSt;
  shSaa1StEnv->DefineSection(3, z, rInSaa1StEnv1, rOuSaa1StEnv2);
  // 2nd Section
  shSaa1StEnv->DefineSection(4, z, rInSaa1StEnv2, rOuSaa1StEnv2);
  z += zSaa1StEnv[1];
  shSaa1StEnv->DefineSection(5, z, rInSaa1StEnv3, rOuSaa1StEnv3);
  // 3rd Section
  z += (zSaa1StEnv[2] - dSt);
  shSaa1StEnv->DefineSection(6, z, rInSaa1StEnv3, rOuSaa1StEnv3);
  // 3 - 4
  shSaa1StEnv->DefineSection(7, z, rInSaa1StEnv3, rOuSaa1StEnv4);
  z += dSt;
  shSaa1StEnv->DefineSection(8, z, rInSaa1StEnv3, rOuSaa1StEnv4);
  // 4th Section
  shSaa1StEnv->DefineSection(9, z, rInSaa1StEnv4, rOuSaa1StEnv4);
  z += zSaa1StEnv[3];
  shSaa1StEnv->DefineSection(10, z, rInSaa1StEnv5, rOuSaa1StEnv5);
  TGeoVolume* voSaa1StEnv = new TGeoVolume("YSAA1_SteelEnvelope", shSaa1StEnv, kMedSteel);
 
  //    SAA1 W-Pipe                //
  //    Drawing ALIP2A__0059       //
  //
  //    Flange FA side
  //    Length of first section
  Float_t dzSaa1WPipeF1 = 0.9;
  //    Outer radius
  Float_t rOuSaa1WPipeF1 = 24.5 / 2.;
  //    Inner Radius
  Float_t rInSaa1WPipeF1 = 22.0 / 2.;
  //    Length of second section
  Float_t dzSaa1WPipeF11 = 2.1;
  //    Inner Radius
  Float_t rInSaa1WPipeF11 = 18.5 / 2.;
  //
  //    Central tube
  //    Length
  Float_t dzSaa1WPipeC = 111.2;
  //    Inner Radius at the end
  Float_t rInSaa1WPipeC = 22.0 / 2.;
  //    Outer Radius
  Float_t rOuSaa1WPipeC = 31.9 / 2.;
  //
  //    Flange SAA2 Side
  //    Length
  Float_t dzSaa1WPipeF2 = 6.0;
  //    Outer radius
  Float_t rOuSaa1WPipeF2 = 41.56 / 2.;
 
  //
  TGeoPcon* shSaa1WPipe = new TGeoPcon(0., 360., 8);
  z = 0.;
  // Flange FA side first section
  shSaa1WPipe->DefineSection(0, z, rInSaa1WPipeF1, rOuSaa1WPipeF1);
  z += dzSaa1WPipeF1;
  shSaa1WPipe->DefineSection(1, z, rInSaa1WPipeF1, rOuSaa1WPipeF1);
  // Flange FA side second section
  shSaa1WPipe->DefineSection(2, z, rInSaa1WPipeF11, rOuSaa1WPipeF1);
  z += dzSaa1WPipeF11;
  shSaa1WPipe->DefineSection(3, z, rInSaa1WPipeF11, rOuSaa1WPipeF1);
  // Central Section
  shSaa1WPipe->DefineSection(4, z, rInSaa1WPipeF11, rOuSaa1WPipeC);
  z += dzSaa1WPipeC;
  shSaa1WPipe->DefineSection(5, z, rInSaa1WPipeC, rOuSaa1WPipeC);
  // Flange SAA2 side
  shSaa1WPipe->DefineSection(6, z, rInSaa1WPipeC, rOuSaa1WPipeF2);
  z += dzSaa1WPipeF2;
  shSaa1WPipe->DefineSection(7, z, rInSaa1WPipeC, rOuSaa1WPipeF2);
 
  TGeoVolume* voSaa1WPipe = new TGeoVolume("YSAA1_WPipe", shSaa1WPipe, kMedNiW);
  //
  // Inner region with higher transport cuts
  TGeoTube* shSaa1WPipeI = new TGeoTube(rInSaa1WPipeC, rOuSaa1WPipeC, dzSaa1WPipeC / 2.);
  TGeoVolume* voSaa1WPipeI = new TGeoVolume("YSAA1_WPipeI", shSaa1WPipeI, kMedNiWsh);
  voSaa1WPipe->AddNode(voSaa1WPipeI, 1, new TGeoTranslation(0., 0., dzSaa1WPipeF1 + dzSaa1WPipeF11 + dzSaa1WPipeC / 2));
 
  //    SAA1 Pb Components         //
  //    Drawing ALIP2A__0078       //
  //
  //    Inner angle
  Float_t tanAlpha = TMath::Tan(1.69 / 2. * kDegRad);
  Float_t tanBeta = TMath::Tan(3.20 / 2. * kDegRad);
  //
  //    1st Section 2deg opening cone
  //    Length
  Float_t dzSaa1PbComp1 = 100.23;
  //    Inner radius at entrance
  Float_t rInSaa1PbComp1 = 22.0 / 2.; // It's 21 cm diameter in the drawing. Is this a typo ??!!
                                      //    Outer radius at entrance
  Float_t rOuSaa1PbComp1 = 42.0 / 2.;
  //
  //    2nd Section: Straight Section
  //    Length
  Float_t dzSaa1PbComp2 = 236.77;
  //    Inner radius
  Float_t rInSaa1PbComp2 = rInSaa1PbComp1 + dzSaa1PbComp1 * tanAlpha;
  //    Outer radius
  Float_t rOuSaa1PbComp2 = 49.0 / 2.;
  //
  //    3rd Section: 1.6deg opening cone until bellow
  //    Length
  Float_t dzSaa1PbComp3 = 175.6;
  //    Inner radius
  Float_t rInSaa1PbComp3 = rInSaa1PbComp2 + dzSaa1PbComp2 * tanAlpha;
  //    Outer radius
  Float_t rOuSaa1PbComp3 = 62.8 / 2.;
  //
  //   4th Section: Bellow region
  Float_t dzSaa1PbComp4 = 26.4;
  //    Inner radius
  Float_t rInSaa1PbComp4 = 37.1 / 2.;
  Float_t rInSaa1PbCompB = 43.0 / 2.;
  //    Outer radius
  Float_t rOuSaa1PbComp4 = rOuSaa1PbComp3 + dzSaa1PbComp3 * tanBeta;
  //
  //   5th Section: Flange SAA2 side
  //   1st detail
  Float_t dzSaa1PbCompF1 = 4.;
  Float_t rOuSaa1PbCompF1 = 74.1 / 2.;
  //   2nd detail
  Float_t dzSaa1PbCompF2 = 3.;
  Float_t rOuSaa1PbCompF2 = 66.0 / 2.;
  Float_t rOuSaa1PbCompF3 = 58.0 / 2.;
 
  TGeoPcon* shSaa1PbComp = new TGeoPcon(0., 360., 11);
  z = 120.2;
  // 2 deg opening cone
  shSaa1PbComp->DefineSection(0, z, rInSaa1PbComp1, rOuSaa1PbComp1);
  z += dzSaa1PbComp1;
  shSaa1PbComp->DefineSection(1, z, rInSaa1PbComp2, rOuSaa1PbComp2);
  // Straight section
  z += dzSaa1PbComp2;
  shSaa1PbComp->DefineSection(2, z, rInSaa1PbComp3, rOuSaa1PbComp2);
  // 1.6 deg opening cone
  shSaa1PbComp->DefineSection(3, z, rInSaa1PbComp3, rOuSaa1PbComp3);
  z += dzSaa1PbComp3;
  shSaa1PbComp->DefineSection(4, z, rInSaa1PbComp4, rOuSaa1PbComp4);
  // Bellow region until outer flange
  shSaa1PbComp->DefineSection(5, z, rInSaa1PbCompB, rOuSaa1PbComp4);
  z += (dzSaa1PbComp4 - dzSaa1PbCompF1 - dzSaa1PbCompF2);
  shSaa1PbComp->DefineSection(6, z, rInSaa1PbCompB, rOuSaa1PbCompF1);
  shSaa1PbComp->DefineSection(7, z, rInSaa1PbCompB, rOuSaa1PbCompF2);
  // Flange first step
  z += dzSaa1PbCompF1;
  shSaa1PbComp->DefineSection(8, z, rInSaa1PbCompB, rOuSaa1PbCompF2);
  shSaa1PbComp->DefineSection(9, z, rInSaa1PbCompB, rOuSaa1PbCompF3);
  // Flange second step
  z += dzSaa1PbCompF2;
  shSaa1PbComp->DefineSection(10, z, rInSaa1PbCompB, rOuSaa1PbCompF3);
 
  TGeoVolume* voSaa1PbComp = new TGeoVolume("YSAA1_PbComp", shSaa1PbComp, kMedPb);
  //
  // Inner region with higher transport cuts
  TGeoPcon* shSaa1PbCompI = MakeShapeFromTemplate(shSaa1PbComp, 0., -3.);
  TGeoVolume* voSaa1PbCompI = new TGeoVolume("YSAA1_PbCompI", shSaa1PbCompI, kMedPbSh);
  voSaa1PbComp->AddNode(voSaa1PbCompI, 1, gGeoIdentity);
 
  //    SAA1 W-Cone                //
  //    Drawing ALIP2A__0058       //
  // Length of the Cone
  Float_t dzSaa1WCone = 52.9;
  // Inner and outer radii
  Float_t rInSaa1WCone1 = 20.4;
  Float_t rOuSaa1WCone1 = rInSaa1WCone1 + 0.97;
  Float_t rOuSaa1WCone2 = rInSaa1WCone1 + 2.80;
  // relative z-position
  Float_t zSaa1WCone = 9.3;
 
  TGeoPcon* shSaa1WCone = new TGeoPcon(0., 360., 2);
  z = zSaa1WCone;
  shSaa1WCone->DefineSection(0, z, rInSaa1WCone1, rOuSaa1WCone1);
  z += dzSaa1WCone;
  shSaa1WCone->DefineSection(1, z, rInSaa1WCone1, rOuSaa1WCone2);
  TGeoVolume* voSaa1WCone = new TGeoVolume("YSAA1_WCone", shSaa1WCone, kMedNiW);
 
  //    SAA1 Steel-Ring            //
  //    Drawing ALIP2A__0040       //
  //
  //    Length of the ring
  Float_t dzSaa1StRing = 4.;
  //    Inner and outer radius
  Float_t rInSaa1String = 33.0;
  Float_t rOuSaa1String = 41.1;
  //    Relative z-position
  Float_t zSaa1StRing = 652.2;
  TGeoPcon* shSaa1StRing = new TGeoPcon(0., 360., 2);
  z = zSaa1StRing;
  shSaa1StRing->DefineSection(0, z, rInSaa1String, rOuSaa1String);
  z += dzSaa1StRing;
  shSaa1StRing->DefineSection(1, z, rInSaa1String, rOuSaa1String);
  TGeoVolume* voSaa1StRing = new TGeoVolume("YSAA1_StRing", shSaa1StRing, kMedSteel);
 
  //    SAA1 Inner Tube            //
  //    Drawing ALIP2A__0082       //
  //
  // Length of saa2:               659.2 cm
  // Length of inner tube:         631.9 cm
  // Lenth of bellow cavern:        27.3 cm
  // Radius at entrance 18.5/2,  d = 0.3
  // Radius at exit     37.1/2,  d = 0.3
  //
  Float_t dzSaa1InnerTube = 631.9 / 2.; // Half length of the tube
  Float_t rInSaa1InnerTube = 18.2 / 2.; // Radius at entrance
  Float_t rOuSaa1InnerTube = 36.8 / 2.; // Radius at exit
  Float_t dSaa1InnerTube = 0.2;         // Thickness
  TGeoVolume* voSaa1InnerTube =
    new TGeoVolume("YSAA1_InnerTube", new TGeoCone(dzSaa1InnerTube, rInSaa1InnerTube - dSaa1InnerTube, rInSaa1InnerTube, rOuSaa1InnerTube - dSaa1InnerTube, rOuSaa1InnerTube),
                   kMedSteelSh);
 
  //    SAA1 Outer Shape           //
  //    Drawing ALIP2A__0107       //
  // Total length
  const Float_t dzSaa1 = 659.2;
  //
  TGeoPcon* shSaa1M = new TGeoPcon(0., 360., 20);
  Float_t kSec = 0.2; // security distance to avoid trivial extrusions
  Float_t rmin = rInSaa1InnerTube - dSaa1InnerTube - kSec;
  rmax = rOuSaa1InnerTube - dSaa1InnerTube - kSec;
  z = 0.;
  shSaa1M->DefineSection(0, z, rmin, rOuSaa1WPipeF1);
  z += dzSaa1WPipeF1;
  shSaa1M->DefineSection(1, z, rmin, rOuSaa1WPipeF1);
  shSaa1M->DefineSection(2, z, 0., rOuSaa1WPipeF1);
  z += dzSaa1WPipeF11;
  shSaa1M->DefineSection(3, z, 0., rOuSaa1WPipeF1);
  shSaa1M->DefineSection(4, z, 0., rOuSaa1StEnv1);
  z = zSaa1WCone;
  shSaa1M->DefineSection(5, z, 0., rOuSaa1StEnv1);
  shSaa1M->DefineSection(6, z, 0., rOuSaa1WCone1);
  z += dzSaa1WCone;
  shSaa1M->DefineSection(7, z, 0., rOuSaa1WCone2);
  shSaa1M->DefineSection(8, z, 0., rOuSaa1StEnv1);
  z = zSaa1StEnv[0] - dSt + zSaa1StEnvS;
  shSaa1M->DefineSection(9, z, 0., rOuSaa1StEnv1);
  shSaa1M->DefineSection(10, z, 0., rOuSaa1StEnv2);
  z += (zSaa1StEnv[1] + dSt);
  shSaa1M->DefineSection(11, z, 0., rOuSaa1StEnv3);
  z += (zSaa1StEnv[2] - dSt);
  shSaa1M->DefineSection(12, z, 0., rOuSaa1StEnv3);
  shSaa1M->DefineSection(13, z, 0., rOuSaa1StEnv4);
 
  z += (zSaa1StEnv[3] - dSt + dzSaa1PbCompF1 + dzSaa1PbCompF2 - dzSaa1PbComp4);
  Float_t rmaxSaa1 = shSaa1M->GetRmax(13) + (z - shSaa1M->GetZ(13)) * TMath::Tan(1.6 * kDegRad);
 
  shSaa1M->DefineSection(14, z, 0., rmaxSaa1);
  shSaa1M->DefineSection(15, z, rmax, rmaxSaa1);
  z = zSaa1StRing;
  shSaa1M->DefineSection(16, z, rmax + 0.4, rOuSaa1String);
  z += dzSaa1PbCompF1;
  shSaa1M->DefineSection(17, z, rmax + 0.4, rOuSaa1String);
  shSaa1M->DefineSection(18, z, rmax + 0.4, rOuSaa1PbCompF3);
  z += dzSaa1PbCompF2;
  shSaa1M->DefineSection(19, z, rmax + 0.4, rOuSaa1PbCompF3);
 
  //
  //    Inner 1.69deg line
  for (Int_t i = 2; i < 15; i++) {
    Double_t zp = shSaa1M->GetZ(i);
    Double_t r2 = shSaa1M->GetRmax(i);
    Double_t r1 = rmin + (zp - 0.9) * TMath::Tan(1.686 / 2. * kDegRad) - kSec;
    shSaa1M->DefineSection(i, zp, r1, r2);
  }
 
  TGeoVolume* voSaa1M = new TGeoVolume("YSAA1M", shSaa1M, kMedAir);
  voSaa1M->SetVisibility(0);
 
  //                               //
  // Recess Station 2              //
  // Drawing ALIP2A__0260          //
  //    SAA1 W-Ring 1              //
  //    Drawing ALIP2A__0217       //
  Float_t saa1Wring1Width = 5.85;
  TGeoPcon* shSaa1Wring1 = new TGeoPcon(0., 360., 2);
  shSaa1Wring1->DefineSection(0, 0.00, 20.31, 23.175);
  shSaa1Wring1->DefineSection(1, saa1Wring1Width, 20.31, 23.400);
  TGeoVolume* voSaa1Wring1 = new TGeoVolume("YSAA1_WRING1", shSaa1Wring1, kMedNiW);
 
  //    SAA1 W-Ring 2              //
  //    Drawing ALIP2A__0055       //
  Float_t saa1Wring2Rinner = 20.31;
  Float_t saa1Wring2Router = 23.40;
  Float_t saa1Wring2HWidth = 3.75;
  Float_t saa1Wring2Cutoffx = 4.9;
  Float_t saa1Wring2Cutoffy = 4.9;
  TGeoTubeSeg* shSaa1Wring2a = new TGeoTubeSeg(saa1Wring2Rinner, saa1Wring2Router, saa1Wring2HWidth, 0., 90.);
  shSaa1Wring2a->SetName("shSaa1Wring2a");
  TGeoBBox* shSaa1Wring2b = new TGeoBBox(saa1Wring2Router / 2., saa1Wring2Router / 2., saa1Wring2HWidth);
  shSaa1Wring2b->SetName("shSaa1Wring2b");
  TGeoTranslation* trSaa1Wring2b = new TGeoTranslation("trSaa1Wring2b", saa1Wring2Router / 2. + saa1Wring2Cutoffx,
                                                       saa1Wring2Router / 2. + saa1Wring2Cutoffy, 0.);
  trSaa1Wring2b->RegisterYourself();
  TGeoCompositeShape* shSaa1Wring2 =
    new TGeoCompositeShape("shSaa1Wring2", "(shSaa1Wring2a)*(shSaa1Wring2b:trSaa1Wring2b)");
  TGeoVolume* voSaa1Wring2 = new TGeoVolume("YSAA1_WRING2", shSaa1Wring2, kMedNiW);
 
  //    SAA1 W-Ring 3              //
  //    Drawing ALIP2A__0216       //
 
  Float_t saa1Wring3Rinner = 20.31;
  Float_t saa1Wring3Router = 23.40;
  Float_t saa1Wring3HWidth = 3.75;
  Float_t saa1Wring3Cutoffx = 4.50;
  Float_t saa1Wring3Cutoffy = 4.50;
  TGeoTubeSeg* shSaa1Wring3a = new TGeoTubeSeg(saa1Wring3Rinner, saa1Wring3Router, saa1Wring3HWidth, 0., 90.);
  shSaa1Wring3a->SetName("shSaa1Wring3a");
  TGeoBBox* shSaa1Wring3b = new TGeoBBox(saa1Wring3Router / 2., saa1Wring3Router / 2., saa1Wring3HWidth);
  shSaa1Wring3b->SetName("shSaa1Wring3b");
  TGeoTranslation* trSaa1Wring3b = new TGeoTranslation("trSaa1Wring3b", saa1Wring3Router / 2. + saa1Wring3Cutoffx,
                                                       saa1Wring3Router / 2. + saa1Wring3Cutoffy, 0.);
  trSaa1Wring3b->RegisterYourself();
  TGeoCompositeShape* shSaa1Wring3 =
    new TGeoCompositeShape("shSaa1Wring3", "(shSaa1Wring3a)*(shSaa1Wring3b:trSaa1Wring3b)");
  TGeoVolume* voSaa1Wring3 = new TGeoVolume("YSAA1_WRING3", shSaa1Wring3, kMedNiW);
 
  //    SAA1 W-Ring 4              //
  //    Drawing ALIP2A__0215       //
  Float_t saa1Wring4Width = 5.85;
  TGeoPcon* shSaa1Wring4 = new TGeoPcon(0., 360., 5);
  shSaa1Wring4->DefineSection(0, 0.00, 20.31, 23.40);
  shSaa1Wring4->DefineSection(1, 1.00, 20.31, 23.40);
  shSaa1Wring4->DefineSection(2, 1.00, 20.31, 24.50);
  shSaa1Wring4->DefineSection(3, 4.85, 20.31, 24.80);
  shSaa1Wring4->DefineSection(4, 5.85, 24.10, 24.80);
  TGeoVolume* voSaa1Wring4 = new TGeoVolume("YSAA1_WRING4", shSaa1Wring4, kMedNiW);
 
  //    SAA1 W-Ring 5              //
  //    Drawing ALIP2A__0218       //
  Float_t saa1Wring5Rinner = 20.31;
  Float_t saa1Wring5Router = 23.40;
  Float_t saa1Wring5HWidth = 0.85;
  TGeoVolume* voSaa1Wring5 =
    new TGeoVolume("YSAA1_WRING5", new TGeoTube(saa1Wring5Rinner, saa1Wring5Router, saa1Wring5HWidth), kMedNiW);
  //
  // Position the rings in the assembly
  //
  TGeoVolumeAssembly* asSaa1ExtraShield = new TGeoVolumeAssembly("YSAA1ExtraShield");
  // Distance between rings
  Float_t saa1DWrings = 2.3;
  //
  dz = -(saa1Wring1Width + 6. * saa1Wring2HWidth + 2. * saa1Wring3HWidth + saa1Wring4Width + 2. * saa1Wring5HWidth +
         2. * saa1DWrings) /
       2.;
  asSaa1ExtraShield->AddNode(voSaa1Wring1, 1, new TGeoTranslation(0., 0., dz));
  dz += saa1Wring1Width;
  dz += saa1Wring2HWidth;
  asSaa1ExtraShield->AddNode(voSaa1Wring2, 1, new TGeoCombiTrans(0., 0., dz, rot000));
  asSaa1ExtraShield->AddNode(voSaa1Wring2, 2, new TGeoCombiTrans(0., 0., dz, rot180));
  dz += saa1Wring2HWidth;
  dz += saa1DWrings;
  dz += saa1Wring2HWidth;
  asSaa1ExtraShield->AddNode(voSaa1Wring2, 3, new TGeoCombiTrans(0., 0., dz, rot090));
  asSaa1ExtraShield->AddNode(voSaa1Wring2, 4, new TGeoCombiTrans(0., 0., dz, rot270));
  dz += saa1Wring2HWidth;
  dz += saa1Wring5HWidth;
  asSaa1ExtraShield->AddNode(voSaa1Wring5, 1, new TGeoTranslation(0., 0., dz));
  dz += saa1Wring5HWidth;
  dz += saa1Wring2HWidth;
  asSaa1ExtraShield->AddNode(voSaa1Wring2, 5, new TGeoCombiTrans(0., 0., dz, rot000));
  asSaa1ExtraShield->AddNode(voSaa1Wring2, 6, new TGeoCombiTrans(0., 0., dz, rot180));
  dz += saa1Wring2HWidth;
  dz += saa1DWrings;
  dz += saa1Wring3HWidth;
  asSaa1ExtraShield->AddNode(voSaa1Wring3, 1, new TGeoCombiTrans(0., 0., dz, rot090));
  asSaa1ExtraShield->AddNode(voSaa1Wring3, 2, new TGeoCombiTrans(0., 0., dz, rot270));
  dz += saa1Wring3HWidth;
  asSaa1ExtraShield->AddNode(voSaa1Wring4, 1, new TGeoTranslation(0., 0., dz));
  dz += saa1Wring4Width;
  const Float_t saa1ExtraShieldL = 48;
  //
  // Assemble SAA1
  voSaa1M->AddNode(voSaa1StEnv, 1, gGeoIdentity);
  voSaa1M->AddNode(voSaa1WPipe, 1, gGeoIdentity);
  voSaa1M->AddNode(voSaa1PbComp, 1, gGeoIdentity);
  voSaa1M->AddNode(voSaa1WCone, 1, gGeoIdentity);
  voSaa1M->AddNode(voSaa1StRing, 1, gGeoIdentity);
  voSaa1M->AddNode(voSaa1InnerTube, 1, new TGeoTranslation(0., 0., dzSaa1InnerTube + 0.9));
  TGeoVolumeAssembly* voSaa1 = new TGeoVolumeAssembly("YSAA1");
  voSaa1->AddNode(voSaa1M, 1, gGeoIdentity);
 
  //          SAA1/SAA2  Pb Joint      //
  //          Drawing ALIP2A__0081     //
  //
  // Outer radius
  Float_t rOuSaa1Saa2 = 70.0 / 2.;
  // Flange SAA1 side
  Float_t dzSaa1Saa2F1 = 3.;
  Float_t rInSaa1Saa2F1 = 58.5 / 2.;
  // 1st Central Section
  Float_t dzSaa1Saa2C1 = 19.3;
  Float_t rInSaa1Saa2C1 = 42.8 / 2.;
  // Transition Region
  Float_t dzSaa1Saa2T = 3.3;
  // 1st Central Section
  Float_t dzSaa1Saa2C2 = 6.2;
  Float_t rInSaa1Saa2C2 = 36.2 / 2.;
  // Flange SAA2 side
  Float_t dzSaa1Saa2F2 = 3.1;
  Float_t rInSaa1Saa2F2 = 54.1 / 2.;
  // Total length
  const Float_t dzSaa1Saa2 = 34.9;
 
  TGeoPcon* shSaa1Saa2Pb = new TGeoPcon(0., 360., 8);
  z = 0.;
  // Flange SAA1 side
  shSaa1Saa2Pb->DefineSection(0, z, rInSaa1Saa2F1, rOuSaa1Saa2);
  z += dzSaa1Saa2F1;
  shSaa1Saa2Pb->DefineSection(1, z, rInSaa1Saa2F1, rOuSaa1Saa2);
  shSaa1Saa2Pb->DefineSection(2, z, rInSaa1Saa2C1, rOuSaa1Saa2);
  // Central region 1
  z += dzSaa1Saa2C1;
  shSaa1Saa2Pb->DefineSection(3, z, rInSaa1Saa2C1, rOuSaa1Saa2);
  // 45 deg transition
  z += dzSaa1Saa2T;
  shSaa1Saa2Pb->DefineSection(4, z, rInSaa1Saa2C2, rOuSaa1Saa2);
  z += dzSaa1Saa2C2;
  shSaa1Saa2Pb->DefineSection(5, z, rInSaa1Saa2C2, rOuSaa1Saa2);
  shSaa1Saa2Pb->DefineSection(6, z, rInSaa1Saa2F2, rOuSaa1Saa2);
  z += dzSaa1Saa2F2;
  shSaa1Saa2Pb->DefineSection(7, z, rInSaa1Saa2F2, rOuSaa1Saa2);
  TGeoVolume* voSaa1Saa2Pb = new TGeoVolume("YSAA1SAA2Pb", shSaa1Saa2Pb, kMedPb);
  //
  //    Mother volume and outer steel envelope
  Float_t rOuSaa1Saa2Steel = 36.9;
 
  TGeoPcon* shSaa1Saa2 = MakeShapeFromTemplate(shSaa1Saa2Pb, 0., rOuSaa1Saa2Steel - rOuSaa1Saa2);
  TGeoVolume* voSaa1Saa2 = new TGeoVolume("YSAA1SAA2", shSaa1Saa2, kMedSteel);
  voSaa1Saa2->AddNode(voSaa1Saa2Pb, 1, gGeoIdentity);
  //
  //    Inner region with higher transport cuts
  //
  TGeoPcon* shSaa1Saa2I = MakeShapeFromTemplate(shSaa1Saa2Pb, 0., -3.);
  TGeoVolume* voSaa1Saa2I = new TGeoVolume("YSAA1_SAA2I", shSaa1Saa2I, kMedPbSh);
  voSaa1Saa2Pb->AddNode(voSaa1Saa2I, 1, gGeoIdentity);
 
  //                SAA2               //
 
  //    SAA2 Steel Envelope        //
  //    Drawing ALIP2A__0041       //
  dSt = 4.; // Thickness of steel envelope
  // Length of the first section
  Float_t dzSaa2StEnv1 = 163.15;
  Float_t rInSaa2StEnv1 = 65.8 / 2.;
  // Length of the second section
  Float_t dzSaa2StEnv2 = 340.35 - 4.;
  Float_t rInSaa2StEnv2 = 87.2 / 2.;
  // Rel. starting position
  Float_t zSaa2StEnv = 3.;
 
  TGeoPcon* shSaa2StEnv = new TGeoPcon(0., 360., 6);
  // First Section
  z = zSaa2StEnv;
  shSaa2StEnv->DefineSection(0, z, rInSaa2StEnv1, rInSaa2StEnv1 + dSt);
  z += dzSaa2StEnv1;
  shSaa2StEnv->DefineSection(1, z, rInSaa2StEnv1, rInSaa2StEnv1 + dSt);
  // Transition region
  shSaa2StEnv->DefineSection(2, z, rInSaa2StEnv1, rInSaa2StEnv2 + dSt);
  z += dSt;
  shSaa2StEnv->DefineSection(3, z, rInSaa2StEnv1, rInSaa2StEnv2 + dSt);
  // Second section
  shSaa2StEnv->DefineSection(4, z, rInSaa2StEnv2, rInSaa2StEnv2 + dSt);
  z += dzSaa2StEnv2;
  shSaa2StEnv->DefineSection(5, z, rInSaa2StEnv2, rInSaa2StEnv2 + dSt);
 
  TGeoVolume* voSaa2StEnv = new TGeoVolume("YSAA2_SteelEnvelope", shSaa2StEnv, kMedSteel);
 
  //    SAA2 Pb Ring               //
  //    Drawing ALIP2A__0080       //
  //    Drawing ALIP2A__0111       //
  //
  // Rel. position in z
  Float_t zSaa2PbRing = 35.25;
  // Length
  Float_t dzSaa2PbRing = 65.90;
  // Inner radius
  Float_t rInSaa2PbRing = 37.00;
  // Outer radius at front
  Float_t rOuSaa2PbRingF = 42.74;
  // Outer Rradius at rear
  Float_t rOuSaa2PbRingR = 44.58;
 
  TGeoPcon* shSaa2PbRing = new TGeoPcon(0., 360., 2);
  z = zSaa2PbRing;
  shSaa2PbRing->DefineSection(0, z, rInSaa2PbRing, rOuSaa2PbRingF);
  z += dzSaa2PbRing;
  shSaa2PbRing->DefineSection(1, z, rInSaa2PbRing, rOuSaa2PbRingR);
 
  TGeoVolume* voSaa2PbRing = new TGeoVolume("YSAA2_PbRing", shSaa2PbRing, kMedPb);
 
  //    SAA2 Pb Components         //
  //    Drawing ALIP2A__0079       //
  tanAlpha = TMath::Tan(1.89 / 2. * kDegRad);
  TGeoPcon* shSaa2PbComp = new TGeoPcon(0., 360., 16);
  // Total length
  const Float_t dzSaa2PbComp = 512.;
  // Length of 1st bellow recess
  Float_t dzSaa2PbCompB1 = 24.;
  // Length of 2nd bellow recess
  Float_t dzSaa2PbCompB2 = 27.;
  // Flange on the SAA1 side Detail A
  // 1st Step
  Float_t dzSaa2PbCompA1 = 1.5;
  Float_t rInSaa2PbCompA1 = 43.0 / 2.;
  Float_t rOuSaa2PbCompA1 = 53.0 / 2.;
  // 2nd Step
  Float_t dzSaa2PbCompA2 = 1.5;
  Float_t rInSaa2PbCompA2 = 36.8 / 2.;
  Float_t rOuSaa2PbCompA2 = rOuSaa2PbCompA1;
  // Straight section
  Float_t dzSaa2PbCompA3 = 21.0;
  Float_t rInSaa2PbCompA3 = rInSaa2PbCompA2;
  Float_t rOuSaa2PbCompA3 = 65.2 / 2.;
  //
  // 1st Section (outer straight, inner 1.89/2. deg opening cone)
  // Length
  Float_t dzSaa2PbComp1 = 146.15;
  // Inner radius at the end
  Float_t rInSaa2PbComp1 = rInSaa2PbCompA3 + dzSaa2PbComp1 * tanAlpha;
  // Outer radius
  Float_t rOuSaa2PbComp1 = rOuSaa2PbCompA3;
  //
  // 2nd Section (outer straight, inner 1.89/2. deg opening cone)
  // Length
  Float_t dzSaa2PbComp2 = (dzSaa2PbComp - dzSaa2PbComp1 - dzSaa2PbCompB1 - dzSaa2PbCompB2);
  // Inner radius at the end
  Float_t rInSaa2PbComp2 = rInSaa2PbComp1 + dzSaa2PbComp2 * tanAlpha;
  // Outer radius
  Float_t rOuSaa2PbComp2 = 86.6 / 2.;
  //
  // Flange on the SAA3 side (Detail E)
  //
  // Straight Section
  // Length  dzSaa2PbCompB2 - 8.8 = 27 - 8.8 = 18.2
  Float_t dzSaa2PbCompE1 = 18.2;
  Float_t rInSaa2PbCompE1 = 52.0 / 2.;
  Float_t rOuSaa2PbCompE1 = 86.6 / 2.;
  // 45 deg transition
  Float_t dzSaa2PbCompE2 = 2.7;
  // 1st Step
  Float_t dzSaa2PbCompE3 = 0.6;
  Float_t rInSaa2PbCompE3 = 52.0 / 2. + dzSaa2PbCompE2;
  Float_t rOuSaa2PbCompE3 = 83.0 / 2.;
  // 2nd Step
  Float_t dzSaa2PbCompE4 = 4.0;
  Float_t rOuSaa2PbCompE4 = 61.6 / 2.;
  // end
  Float_t dzSaa2PbCompE5 = 1.5;
 
  //
  // Flange on SAA1 side (Detail A)
  z = 0.;
  // 1st Step
  shSaa2PbComp->DefineSection(0, z, rInSaa2PbCompA1, rOuSaa2PbCompA1);
  z += dzSaa2PbCompA1;
  shSaa2PbComp->DefineSection(1, z, rInSaa2PbCompA1, rOuSaa2PbCompA1);
  shSaa2PbComp->DefineSection(2, z, rInSaa2PbCompA2, rOuSaa2PbCompA2);
  // 2nd Step
  z += dzSaa2PbCompA2;
  shSaa2PbComp->DefineSection(3, z, rInSaa2PbCompA2, rOuSaa2PbCompA2);
  shSaa2PbComp->DefineSection(4, z, rInSaa2PbCompA3, rOuSaa2PbCompA3);
  // straight section
  z += dzSaa2PbCompA3;
  shSaa2PbComp->DefineSection(5, z, rInSaa2PbCompA3, rOuSaa2PbCompA3);
  //
  // Section 1
  z += dzSaa2PbComp1;
  shSaa2PbComp->DefineSection(6, z, rInSaa2PbComp1, rOuSaa2PbComp1);
  shSaa2PbComp->DefineSection(7, z, rInSaa2PbComp1, rOuSaa2PbComp2);
  //
  // Section 2
  z += dzSaa2PbComp2;
  shSaa2PbComp->DefineSection(8, z, rInSaa2PbComp2, rOuSaa2PbComp2);
  //
  // Flange SAA3 side (Detail E)
  z += dzSaa2PbCompE1;
  shSaa2PbComp->DefineSection(9, z, rInSaa2PbCompE1, rOuSaa2PbCompE1);
  // 45 deg transition
  z += dzSaa2PbCompE2;
  shSaa2PbComp->DefineSection(10, z, rInSaa2PbCompE3, rOuSaa2PbCompE1);
  // 1st step
  z += dzSaa2PbCompE3;
  shSaa2PbComp->DefineSection(11, z, rInSaa2PbCompE3, rOuSaa2PbCompE1);
  shSaa2PbComp->DefineSection(12, z, rInSaa2PbCompE3, rOuSaa2PbCompE3);
  // 2nd step
  z += dzSaa2PbCompE4;
  shSaa2PbComp->DefineSection(13, z, rInSaa2PbCompE3, rOuSaa2PbCompE3);
  shSaa2PbComp->DefineSection(14, z, rInSaa2PbCompE3, rOuSaa2PbCompE4);
  // end
  z += dzSaa2PbCompE5;
  shSaa2PbComp->DefineSection(15, z, rInSaa2PbCompE3, rOuSaa2PbCompE4);
 
  TGeoVolume* voSaa2PbComp = new TGeoVolume("YSAA2_PbComp", shSaa2PbComp, kMedPbSh);
 
  //    SAA2 Inner Tube            //
  //    Drawing ALIP2A__0083       //
  //
  //
  //
  // Length of saa2:               512.0 cm
  // Length of inner tube:         501.7 cm
  // Lenth of bellow recess:        10.3 cm   ( 1.5 + 8.8)
  // Radius at entrance 36.8/2,  d = 0.1
  // Radius at exit     52.0/2,  d = 0.1
  //
  const Float_t kSaa2InnerTubeL = 501.7;        // Length of the tube
  const Float_t kSaa2InnerTubeRmin = 36.6 / 2.; // Radius at entrance
  const Float_t kSaa2InnerTubeRmax = 51.8 / 2.; // Radius at exit
  const Float_t kSaa2InnerTubeD = 0.2;          // Thickness
  TGeoPcon* shSaa2InnerTube = new TGeoPcon(0., 360., 4);
  z = 0.;
  shSaa2InnerTube->DefineSection(0, z, kSaa2InnerTubeRmin - kSaa2InnerTubeD, kSaa2InnerTubeRmin);
  z += dzSaa2PbCompA2 + dzSaa2PbCompA3;
  shSaa2InnerTube->DefineSection(1, z, kSaa2InnerTubeRmin - kSaa2InnerTubeD, kSaa2InnerTubeRmin);
  z = kSaa2InnerTubeL - dzSaa2PbCompE1;
  shSaa2InnerTube->DefineSection(2, z, kSaa2InnerTubeRmax - kSaa2InnerTubeD, kSaa2InnerTubeRmax);
  z = kSaa2InnerTubeL;
  shSaa2InnerTube->DefineSection(3, z, kSaa2InnerTubeRmax - kSaa2InnerTubeD, kSaa2InnerTubeRmax);
  TGeoVolume* voSaa2InnerTube = new TGeoVolume("YSAA2_InnerTube", shSaa2InnerTube, kMedSteelSh);
 
  //    SAA2 Steel Ring            //
  //    Drawing ALIP2A__0042       //
  //  HalfWidth
  Float_t dzSaa2SteelRing = 2.;
  TGeoTube* shSaa2SteelRing = new TGeoTube(41.6, 47.6, dzSaa2SteelRing);
  TGeoVolume* voSaa2SteelRing = new TGeoVolume("YSAA2_SteelRing", shSaa2SteelRing, kMedSteel);
 
  //    SAA2 Outer Shape           //
  //    Drawing ALIP2A__0108       //
 
  TGeoPcon* shSaa2 = new TGeoPcon(0., 360., 16);
  kSec = 0.02; // security distance to avoid trivial extrusions
  rmin = kSaa2InnerTubeRmin - kSaa2InnerTubeD - kSec;
  rmax = kSaa2InnerTubeRmax - kSaa2InnerTubeD - kSec;
  // Flange SAA1 side
  z = 0.;
  shSaa2->DefineSection(0, z, rmin, rOuSaa2PbCompA1);
  z += dzSaa2PbCompA1 + dzSaa2PbCompA2;
  shSaa2->DefineSection(1, z, rmin, rOuSaa2PbCompA1);
  shSaa2->DefineSection(2, z, rmin, rInSaa2StEnv1 + dSt);
  z += dzSaa2PbCompA3;
  shSaa2->DefineSection(3, z, rmin, rInSaa2StEnv1 + dSt);
  z = zSaa2PbRing;
  shSaa2->DefineSection(4, z, 0., rInSaa2StEnv1 + dSt);
  shSaa2->DefineSection(5, z, 0., rOuSaa2PbRingF);
  z += dzSaa2PbRing;
  shSaa2->DefineSection(6, z, 0., rOuSaa2PbRingR);
  shSaa2->DefineSection(7, z, 0., rInSaa2StEnv1 + dSt);
  z = dzSaa2PbCompA1 + dzSaa2PbCompA2 + dzSaa2StEnv1;
  shSaa2->DefineSection(8, z, 0., rInSaa2StEnv1 + dSt);
  shSaa2->DefineSection(9, z, 0., rInSaa2StEnv2 + dSt);
  z = dzSaa2PbComp - dzSaa2PbCompB2;
  shSaa2->DefineSection(10, z, rmax, rInSaa2StEnv2 + dSt);
  z += dzSaa2PbCompE1;
  shSaa2->DefineSection(11, z, rmax, rInSaa2StEnv2 + dSt);
  z += dzSaa2PbCompE2;
  shSaa2->DefineSection(12, z, rInSaa2PbCompE3, rInSaa2StEnv2 + dSt);
  z += (dzSaa2PbCompE3 + dzSaa2PbCompE4);
  shSaa2->DefineSection(13, z, rInSaa2PbCompE3, rInSaa2StEnv2 + dSt);
  shSaa2->DefineSection(14, z, rInSaa2PbCompE3, rOuSaa2PbCompE4);
  z += dzSaa2PbCompE5;
  shSaa2->DefineSection(15, z, rInSaa2PbCompE3, rOuSaa2PbCompE4);
 
  TGeoVolume* voSaa2 = new TGeoVolume("YSAA2", shSaa2, kMedAir);
  voSaa2->SetVisibility(0);
  // Inner 1.89/2 deg line
  Double_t zref = dzSaa2PbCompA1 + dzSaa2PbCompA2 + dzSaa2PbCompA3;
  for (Int_t i = 4; i < 10; i++) {
    Double_t zp = shSaa2->GetZ(i);
    Double_t r2 = shSaa2->GetRmax(i);
    Double_t r1 = rmin + (zp - zref) * TMath::Tan(1.89 / 2. * kDegRad) - kSec;
    shSaa2->DefineSection(i, zp, r1, r2);
  }
 
  //
  //    Assemble SAA2
  voSaa2->AddNode(voSaa2StEnv, 1, gGeoIdentity);
  voSaa2->AddNode(voSaa2PbRing, 1, gGeoIdentity);
  voSaa2->AddNode(voSaa2PbComp, 1, gGeoIdentity);
  voSaa2->AddNode(voSaa2InnerTube, 1, new TGeoTranslation(0., 0., dzSaa2PbCompA1));
  z = (dzSaa2PbComp - dzSaa2PbCompE4 - dzSaa2PbCompE5) + dzSaa2SteelRing;
  voSaa2->AddNode(voSaa2SteelRing, 1, new TGeoTranslation(0., 0., z));
 
  //                SAA3               //
  //
  //
  //  This is a study performed by S. Maridor
  //  The SAA3 has not yet been designed !!!!!!!!
  //
  //    SAA3 Outer Shape           //
  //    Drawing ALIP2A__0288       //
  TGeoVolumeAssembly* voSaa3 = new TGeoVolumeAssembly("YSAA3");
  //    SAA3 Concrete cone         //
  //    Drawing ALIP2A__0284       //
  //    Block
  TGeoBBox* shSaa3CCBlockO = new TGeoBBox(80. / 2., 80. / 2., 80. / 2.); // modified by ernesto.calvo@pucp.edu.pe
  shSaa3CCBlockO->SetName("Saa3CCBlockO");
 
  TGeoPcon* shSaa3InnerRegion = new TGeoPcon(0., 360., 4);
  // + 10.0 cm // ecv
  shSaa3InnerRegion->DefineSection(0, -50.0, 0., 27.1);
  shSaa3InnerRegion->DefineSection(1, -13.0, 0., 27.1);
  shSaa3InnerRegion->DefineSection(2, 39.1, 0., 12.3);
  shSaa3InnerRegion->DefineSection(3, 80.0, 0., 12.3);
  shSaa3InnerRegion->SetName("Saa3InnerRegion");
 
  TGeoCompositeShape* shSaa3CCBlock = new TGeoCompositeShape("Saa3CCBlock", "Saa3CCBlockO-Saa3InnerRegion");
  TGeoVolume* voSaa3CCBlock = new TGeoVolume("YSAA3CCBlock", shSaa3CCBlock, kMedConcSh);
 
  voSaa3->AddNode(voSaa3CCBlock, 1, gGeoIdentity);
 
  //    Plate 1: 240 cm x 80 cm x  80 cm (x 2)
  TGeoVolume* voSaa3SteelPlate1 =
    new TGeoVolume("YSAA3SteelPlate1", new TGeoBBox(240. / 2., 80. / 2., 80. / 2.), kMedSteelSh);
  TGeoVolume* voSaa3SteelPlate11 =
    new TGeoVolume("YSAA3SteelPlate11", new TGeoBBox(240. / 2., 80. / 2., 10. / 2.), kMedSteel);
  voSaa3SteelPlate11->SetVisContainers(kTRUE);
  voSaa3SteelPlate1->SetVisibility(kTRUE);
  voSaa3SteelPlate1->AddNode(voSaa3SteelPlate11, 1, new TGeoTranslation(0., 0., -35.));
  voSaa3->AddNode(voSaa3SteelPlate1, 1, new TGeoTranslation(0., +80., 0.));
  voSaa3->AddNode(voSaa3SteelPlate1, 2, new TGeoTranslation(0., -80., 0.));
 
  //    Plate 2:  80 cm x 80 cm x  80 cm (x 2)
  TGeoVolume* voSaa3SteelPlate2 =
    new TGeoVolume("YSAA3SteelPlate2", new TGeoBBox(80. / 2., 80. / 2., 80. / 2.), kMedSteelSh);
  TGeoVolume* voSaa3SteelPlate21 =
    new TGeoVolume("YSAA3SteelPlate21", new TGeoBBox(80. / 2., 80. / 2., 10. / 2.), kMedSteel);
  voSaa3SteelPlate2->AddNode(voSaa3SteelPlate21, 1, new TGeoTranslation(0., 0., -35.));
  voSaa3->AddNode(voSaa3SteelPlate2, 1, new TGeoTranslation(+80, 0., 0.));
  voSaa3->AddNode(voSaa3SteelPlate2, 2, new TGeoTranslation(-80, 0., 0.));
 
  //    Muon Filter                //
  //    Drawing ALIP2A__0105       //
  // Half Length
  const Float_t dzMuonFilter = 60.;
 
  TGeoBBox* shMuonFilterO = new TGeoBBox(550. / 2., 620. / 2., dzMuonFilter);
  shMuonFilterO->SetName("FilterO");
  TGeoCombiTrans* trFilter = new TGeoCombiTrans("trFilter", 0., -dzMuonFilter * TMath::Tan(alhc * kDegrad), 0., rotlhc);
  trFilter->RegisterYourself();
  TGeoTube* shMuonFilterI = new TGeoTube(0., 48.8, dzMuonFilter + 20.);
  shMuonFilterI->SetName("FilterI");
  TGeoCompositeShape* shMuonFilter = new TGeoCompositeShape("MuonFilter", "FilterO-FilterI:trFilter");
  //
  // !!!!! Needs to be inclined
  TGeoVolume* voMuonFilter = new TGeoVolume("YMuonFilter", shMuonFilter, kMedCastiron);
 
  // Inner part with higher transport cuts
  Float_t dzMuonFilterH = 50.;
  TGeoBBox* shMuonFilterOH = new TGeoBBox(550. / 2., 620. / 2., dzMuonFilterH);
  shMuonFilterOH->SetName("FilterOH");
  TGeoTube* shMuonFilterIH = new TGeoTube(0., 50., dzMuonFilterH + 5.);
  shMuonFilterIH->SetName("FilterIH");
  TGeoCompositeShape* shMuonFilterH = new TGeoCompositeShape("MuonFilterH", "FilterOH-FilterIH:trFilter");
  TGeoVolume* voMuonFilterH = new TGeoVolume("YMuonFilterH", shMuonFilterH, kMedCastironSh);
  voMuonFilter->AddNode(voMuonFilterH, 1, gGeoIdentity);
 
  //
  TGeoVolumeAssembly* voSaa = new TGeoVolumeAssembly("YSAA");
  //
  //    Starting position of the FA Flange/Tail
  const Float_t ziFaWTail = 505.0;
  //    End of the FA Flange/Tail
  const Float_t zoFaWTail = ziFaWTail + dzFaWTail;
  //    Starting position of the FA/SAA1 Joint (2.8 cm overlap with tail)
  const Float_t ozFaSaa1 = 2.8;
  const Float_t ziFaSaa1 = zoFaWTail - ozFaSaa1;
  //    End of the FA/SAA1 Joint
  const Float_t zoFaSaa1 = ziFaSaa1 + dzFaSaa1;
  //    Starting position of SAA1 (2.0 cm overlap with joint)
  const Float_t ozSaa1 = 2.;
  const Float_t ziSaa1 = zoFaSaa1 - ozSaa1;
  //    End of SAA1
  const Float_t zoSaa1 = ziSaa1 + dzSaa1;
  //    Starting position of SAA1/SAA2 Joint (1.95 cm overlap with SAA1)
  const Float_t ziSaa1Saa2 = zoSaa1 - 1.95;
  //    End of SAA1/SAA2 Joint
  const Float_t zoSaa1Saa2 = ziSaa1Saa2 + dzSaa1Saa2;
  //    Starting position of SAA2 (3.1 cm overlap with the joint)
  const Float_t ziSaa2 = zoSaa1Saa2 - 3.1;
  //    End of SAA2
  const Float_t zoSaa2 = ziSaa2 + dzSaa2PbComp;
  //    Position of SAA3
  const Float_t zcSaa3 =
    zoSaa2 + 40.; // ecv (put the CC block 5cm closer to the Wall)  // <--cancell that!! (cancelled)
                  //    Position of the Muon Filter
  const Float_t zcFilter = 1465.9 + dzMuonFilter;
 
  voSaa->AddNode(voFaWTail, 1, new TGeoTranslation(0., 0., ziFaWTail));
  voSaa->AddNode(voFaSaa1, 1, new TGeoTranslation(0., 0., ziFaSaa1));
  voSaa->AddNode(voSaa1, 1, new TGeoTranslation(0., 0., ziSaa1));
  voSaa->AddNode(voSaa1Saa2, 1, new TGeoTranslation(0., 0., ziSaa1Saa2 - 0.1));
  voSaa->AddNode(voSaa2, 1, new TGeoTranslation(0., 0., ziSaa2));
 
  //
  //    Add Saa3 to Saa and Saa to top volume
  //
  voSaa->AddNode(voSaa3, 1, new TGeoTranslation(0., 0., zcSaa3));
 
  TGeoRotation* rotxz = new TGeoRotation("rotxz", 90., 0., 90., 90., 180., 0.);
  top->AddNode(voSaa, 1, new TGeoCombiTrans(0., 0., 0., rotxz));
 
  //
  //  Mother volume for muon stations 1+2 and shielding material placed between the quadrants
  //
  // Position of the dipole
  Float_t ziDipole = 741.;
 
  TGeoPcon* shYOUT1 = new TGeoPcon(0., 360., 22);
  Float_t eps = 1.e-2;
  // FA Tail Section
  for (Int_t iz = 1; iz < 9; iz++) {
    z = shFaWTail->GetZ(iz - 1);
    if (iz == 8) {
      z -= ozFaSaa1;
    }
    shYOUT1->DefineSection(iz - 1, z + ziFaWTail, shFaWTail->GetRmax(iz - 1) + eps, 150.);
  }
  // FA-SAA1 Joint
  z = shYOUT1->GetZ(7);
  Int_t izc = 8;
  for (Int_t iz = 9; iz < 17; iz++) {
    if (iz == 11 || iz == 15) {
      continue;
    }
    shYOUT1->DefineSection(izc, z + shFaSaa1->GetZ(iz - 9), shFaSaa1->GetRmax(iz - 9) + eps, 150.);
    izc++;
  }
  z = shYOUT1->GetZ(13) - ozSaa1;
  // SAA1  - Dipole
  for (Int_t iz = 15; iz < 22; iz++) {
    shYOUT1->DefineSection(iz - 1, z + shSaa1M->GetZ(iz - 11), shSaa1M->GetRmax(iz - 11) + eps, 150.);
  }
  // Distance between dipole and start of SAA1 2deg opening cone
  dz = ziDipole - (zSaa1StEnv[0] - dSt + zSaa1StEnvS + ziSaa1);
  rOut = rOuSaa1StEnv2 + dz * TMath::Tan(2. * kDegRad);
 
  shYOUT1->DefineSection(21, ziDipole, rOut + eps, 150.);
 
  InvertPcon(shYOUT1);
  TGeoVolume* voYOUT1 = new TGeoVolume("YOUT1", shYOUT1, kMedAirMu);
  voYOUT1->SetVisibility(0);
 
  voYOUT1->AddNode(asSaa1ExtraShield, 1,
                   new TGeoCombiTrans(0., 0., -(100.7 + 56.2 + saa1ExtraShieldL / 2. + ziFaWTail), rotxz));
  voYOUT1->AddNode(asFaExtraShield, 1,
                   new TGeoCombiTrans(0., 0., -(16.41 - 7.46 + kFaWring2HWidth + ziFaWTail), rotxz));
  top->AddNode(voYOUT1, 1, gGeoIdentity);
  //
  //  Mother volume for muon stations 4+5 and trigger stations.
  //
  Float_t zoDipole = 1249.;
 
  TGeoPcon* shYOUT21 = new TGeoPcon(0., 360., 12);
  z = zoDipole;
  shYOUT21->DefineSection(0, z, rOuSaa1String, 375.);
  //    Start of SAA1-SAA2
  z = ziSaa1Saa2;
  shYOUT21->DefineSection(1, z, rOuSaa1String, 375.);
  shYOUT21->DefineSection(2, z, rOuSaa1Saa2Steel, 375.);
  //    End of SAA1-SAA2
  z = ziSaa2;
  shYOUT21->DefineSection(3, z, rOuSaa1Saa2Steel, 375.);
  //    SAA2
  z = ziSaa2 + zSaa2PbRing;
  shYOUT21->DefineSection(4, z, rInSaa2StEnv1 + dSt, 375.);
  //    Pb Cone
  shYOUT21->DefineSection(5, z, rOuSaa2PbRingF, 375.);
  rmin = rOuSaa2PbRingF + (1380. - z) * TMath::Tan(1.6 * kDegRad);
  shYOUT21->DefineSection(6, 1380., rmin, 375.);
  z = ziSaa2 + zSaa2PbRing + dzSaa2PbRing;
  shYOUT21->DefineSection(7, z, rOuSaa2PbRingR, 375.);
  //    Straight Sections
  shYOUT21->DefineSection(8, z, rInSaa2StEnv1 + dSt, 460.);
  z = ziSaa2 + dzSaa2StEnv1;
  shYOUT21->DefineSection(9, z, rInSaa2StEnv1 + dSt, 460.);
  shYOUT21->DefineSection(10, z, rInSaa2StEnv2 + dSt, 460.);
  z += dzSaa2StEnv2;
  shYOUT21->DefineSection(11, z, rInSaa2StEnv2 + dSt, 460.);
 
  InvertPcon(shYOUT21);
  shYOUT21->SetName("shYOUT21");
 
  TGeoBBox* shYOUT22 = new TGeoBBox(460., 200., 65. - 1.5);
  shYOUT22->SetName("shYOUT22");
 
  TGeoTranslation* tYOUT22 = new TGeoTranslation(0., -310. - 200., -zcFilter);
  tYOUT22->SetName("tYOUT22");
  tYOUT22->RegisterYourself();
 
  TGeoCompositeShape* shYOUT2 = new TGeoCompositeShape("shYOUT2", "shYOUT21-shYOUT22:tYOUT22");
 
  TGeoVolume* voYOUT2 = new TGeoVolume("YOUT2", shYOUT2, kMedAirMu);
  voYOUT2->SetVisibility(1);
  voYOUT2->AddNode(voMuonFilter, 1,
                   new TGeoCombiTrans(0., dzMuonFilter * TMath::Tan(alhc * kDegrad), -zcFilter, rotxzlhc));
  top->AddNode(voYOUT2, 1, gGeoIdentity);
}
 
void Shil::createMaterials()
{
  auto& matmgr = o2::base::MaterialManager::Instance();
 
  //
  // Defines materials for the muon shield
  //
  Int_t isxfld1 = 2.;
  Float_t sxmgmx = 10.;
  o2::base::Detector::initFieldTrackingParams(isxfld1, sxmgmx);
  Int_t isxfld2 = 2.; // ((AliMagF*)TGeoGlobalMagField::Instance()->GetField())->PrecInteg();
 
  // Steel
  Float_t asteel[4] = {55.847, 51.9961, 58.6934, 28.0855};
  Float_t zsteel[4] = {26., 24., 28., 14.};
  Float_t wsteel[4] = {.715, .18, .1, .005};
  // PbW
  Float_t apbw[2] = {207.2, 183.85};
  Float_t zpbw[2] = {82., 74.};
  Float_t wpbw[2] = {.5, .5};
  // Concrete
  Float_t aconc[10] = {1., 12.01, 15.994, 22.99, 24.305, 26.98, 28.086, 39.1, 40.08, 55.85};
  Float_t zconc[10] = {1., 6., 8., 11., 12., 13., 14., 19., 20., 26.};
  Float_t wconc[10] = {.01, .001, .529107, .016, .002, .033872, .337021, .013, .044, .014};
  // Ni-Cu-W alloy
  Float_t aniwcu[3] = {58.6934, 183.84, 63.546};
  Float_t zniwcu[3] = {28., 74., 29.};
  Float_t wniwcu[3] = {0.015, 0.95, 0.035};
  //
  // Insulation powder
  //                    Si         O       Ti     Al
  Float_t ains[4] = {28.0855, 15.9994, 47.867, 26.982};
  Float_t zins[4] = {14., 8., 22., 13.};
  Float_t wins[4] = {0.3019, 0.4887, 0.1914, 0.018};
  //
  // Air
  //
  Float_t aAir[4] = {12.0107, 14.0067, 15.9994, 39.948};
  Float_t zAir[4] = {6., 7., 8., 18.};
  Float_t wAir[4] = {0.000124, 0.755267, 0.231781, 0.012827};
  Float_t dAir = 1.20479E-3;
  Float_t dAir1 = 1.20479E-10;
  //
  // Cast iron
  //
  Float_t acasti[4] = {55.847, 12.011, 28.085, 54.938};
  Float_t zcasti[4] = {26., 6., 14., 25.};
  Float_t wcasti[4] = {0.929, 0.035, 0.031, 0.005};
 
  // ****************
  //     Defines tracking media parameters.
  //     Les valeurs sont commentees pour laisser le defaut
  //     a GEANT (version 3-21, page CONS200), f.m.
  Float_t epsil, stmin, tmaxfd, deemax, stemax;
  epsil = .001;  // Tracking precision,
  stemax = -1.;  // Maximum displacement for multiple scat
  tmaxfd = -20.; // Maximum angle due to field deflection
  deemax = -.3;  // Maximum fractional energy loss, DLS
  stmin = -.8;
  // ***************
 
  //    Aluminum
  matmgr.Material("SHIL", 9, "ALU1      ", 26.98, 13., 2.7, 8.9, 37.2);
  matmgr.Material("SHIL", 29, "ALU2      ", 26.98, 13., 2.7, 8.9, 37.2);
  matmgr.Material("SHIL", 49, "ALU3      ", 26.98, 13., 2.7, 8.9, 37.2);
  matmgr.Medium("SHIL", 9, "ALU_C0          ", 9, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
  matmgr.Medium("SHIL", 29, "ALU_C1          ", 29, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
  matmgr.Medium("SHIL", 49, "ALU_C2          ", 49, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
 
  //    Iron
  matmgr.Material("SHIL", 10, "IRON1     ", 55.85, 26., 7.87, 1.76, 17.1);
  matmgr.Material("SHIL", 30, "IRON2     ", 55.85, 26., 7.87, 1.76, 17.1);
  matmgr.Material("SHIL", 50, "IRON3     ", 55.85, 26., 7.87, 1.76, 17.1);
  matmgr.Medium("SHIL", 10, "FE_C0           ", 10, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
  matmgr.Medium("SHIL", 30, "FE_C1           ", 30, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
  matmgr.Medium("SHIL", 50, "FE_C2           ", 50, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
 
  //    Copper
  matmgr.Material("SHIL", 11, "COPPER1   ", 63.55, 29., 8.96, 1.43, 15.1);
  matmgr.Material("SHIL", 31, "COPPER2   ", 63.55, 29., 8.96, 1.43, 15.1);
  matmgr.Material("SHIL", 51, "COPPER3   ", 63.55, 29., 8.96, 1.43, 15.1);
  matmgr.Medium("SHIL", 11, "Cu_C0           ", 11, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
  matmgr.Medium("SHIL", 31, "Cu_C1           ", 31, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
  matmgr.Medium("SHIL", 51, "Cu_C2           ", 51, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
 
  //    Tungsten
  matmgr.Material("SHIL", 12, "TUNGSTEN1 ", 183.85, 74., 19.3, .35, 10.3);
  matmgr.Material("SHIL", 32, "TUNGSTEN2 ", 183.85, 74., 19.3, .35, 10.3);
  matmgr.Material("SHIL", 52, "TUNGSTEN3 ", 183.85, 74., 19.3, .35, 10.3);
  matmgr.Medium("SHIL", 12, "W_C0            ", 12, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
  matmgr.Medium("SHIL", 32, "W_C1            ", 32, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
  matmgr.Medium("SHIL", 52, "W_C2            ", 52, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
 
  //    Lead
  matmgr.Material("SHIL", 13, "LEAD1     ", 207.19, 82., 11.35, .56, 18.5);
  matmgr.Material("SHIL", 33, "LEAD2     ", 207.19, 82., 11.35, .56, 18.5);
  matmgr.Material("SHIL", 53, "LEAD3     ", 207.19, 82., 11.35, .56, 18.5);
  matmgr.Medium("SHIL", 13, "PB_C0           ", 13, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
  matmgr.Medium("SHIL", 33, "PB_C1           ", 33, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
  matmgr.Medium("SHIL", 53, "PB_C2           ", 53, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
 
  //    Insulation Powder
  matmgr.Mixture("SHIL", 14, "INSULATION1", ains, zins, 0.41, 4, wins);
  matmgr.Mixture("SHIL", 34, "INSULATION2", ains, zins, 0.41, 4, wins);
  matmgr.Mixture("SHIL", 54, "INSULATION3", ains, zins, 0.41, 4, wins);
  matmgr.Medium("SHIL", 14, "INS_C0          ", 14, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
  matmgr.Medium("SHIL", 34, "INS_C1          ", 34, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
  matmgr.Medium("SHIL", 54, "INS_C2          ", 54, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
 
  //    Air
  matmgr.Mixture("SHIL", 15, "AIR1      ", aAir, zAir, dAir, 4, wAir);
  matmgr.Mixture("SHIL", 35, "AIR2      ", aAir, zAir, dAir, 4, wAir);
  matmgr.Mixture("SHIL", 55, "AIR3      ", aAir, zAir, dAir, 4, wAir);
  matmgr.Mixture("SHIL", 75, "AIR_MUON  ", aAir, zAir, dAir, 4, wAir);
  matmgr.Medium("SHIL", 15, "AIR_C0          ", 15, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
  matmgr.Medium("SHIL", 35, "AIR_C1          ", 35, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
  matmgr.Medium("SHIL", 55, "AIR_C2          ", 55, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
  matmgr.Medium("SHIL", 75, "AIR_MUON        ", 75, 0, isxfld2, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
 
  //    Vacuum
  matmgr.Mixture("SHIL", 16, "VACUUM1 ", aAir, zAir, dAir1, 4, wAir);
  matmgr.Mixture("SHIL", 36, "VACUUM2 ", aAir, zAir, dAir1, 4, wAir);
  matmgr.Mixture("SHIL", 56, "VACUUM3 ", aAir, zAir, dAir1, 4, wAir);
  matmgr.Medium("SHIL", 16, "VA_C0           ", 16, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
  matmgr.Medium("SHIL", 36, "VA_C1           ", 36, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
  matmgr.Medium("SHIL", 56, "VA_C2           ", 56, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
 
  //     Stainless Steel
  matmgr.Mixture("SHIL", 19, "STAINLESS STEEL1", asteel, zsteel, 7.88, 4, wsteel);
  matmgr.Mixture("SHIL", 39, "STAINLESS STEEL2", asteel, zsteel, 7.88, 4, wsteel);
  matmgr.Mixture("SHIL", 59, "STAINLESS STEEL3", asteel, zsteel, 7.88, 4, wsteel);
  matmgr.Medium("SHIL", 19, "ST_C0           ", 19, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
  matmgr.Medium("SHIL", 39, "ST_C1           ", 39, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
  matmgr.Medium("SHIL", 59, "ST_C3           ", 59, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
 
  //    Lead/Tungsten
  matmgr.Mixture("SHIL", 20, "LEAD/TUNGSTEN1", apbw, zpbw, 15.325, 2, wpbw);
  matmgr.Mixture("SHIL", 40, "LEAD/TUNGSTEN2", apbw, zpbw, 15.325, 2, wpbw);
  matmgr.Mixture("SHIL", 60, "LEAD/TUNGSTEN3", apbw, zpbw, 15.325, 2, wpbw);
  matmgr.Medium("SHIL", 20, "PB/W0           ", 20, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
  matmgr.Medium("SHIL", 40, "PB/W1           ", 40, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
  matmgr.Medium("SHIL", 60, "PB/W3           ", 60, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
 
  //    Ni/Tungsten
  //     Ni-W-Cu
  matmgr.Mixture("SHIL", 21, "Ni-W-Cu1", aniwcu, zniwcu, 18.78, 3, wniwcu);
  matmgr.Mixture("SHIL", 41, "Ni-W-Cu2", aniwcu, zniwcu, 18.78, 3, wniwcu);
  matmgr.Mixture("SHIL", 61, "Ni-W-Cu3", aniwcu, zniwcu, 18.78, 3, wniwcu);
  matmgr.Medium("SHIL", 21, "Ni/W0           ", 21, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
  matmgr.Medium("SHIL", 41, "Ni/W1           ", 41, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
  matmgr.Medium("SHIL", 61, "Ni/W3           ", 61, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
 
  //    Concrete
  matmgr.Mixture("SHIL", 17, "CONCRETE1", aconc, zconc, 2.35, 10, wconc);
  matmgr.Mixture("SHIL", 37, "CONCRETE2", aconc, zconc, 2.35, 10, wconc);
  matmgr.Mixture("SHIL", 57, "CONCRETE3", aconc, zconc, 2.35, 10, wconc);
  matmgr.Medium("SHIL", 17, "CC_C0           ", 17, 0, 0, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
  matmgr.Medium("SHIL", 37, "CC_C1           ", 37, 0, 0, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
  matmgr.Medium("SHIL", 57, "CC_C2           ", 57, 0, 0, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
 
  //    Cast iron
  matmgr.Mixture("SHIL", 18, "CAST IRON1", acasti, zcasti, 7.2, 4, wcasti);
  matmgr.Mixture("SHIL", 38, "CAST IRON2", acasti, zcasti, 7.2, 4, wcasti);
  matmgr.Mixture("SHIL", 58, "CAST IRON3", acasti, zcasti, 7.2, 4, wcasti);
  matmgr.Medium("SHIL", 18, "CAST_IRON0      ", 18, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
  matmgr.Medium("SHIL", 38, "CAST_IRON1      ", 38, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
  matmgr.Medium("SHIL", 58, "CAST_IRON2      ", 58, 0, isxfld1, sxmgmx, tmaxfd, stemax, deemax, epsil, stmin);
}
 
FairModule* Shil::CloneModule() const { return new Shil(*this); }
ClassImp(o2::passive::Shil);