New default values for baselines (F.Prino)

[u/mrichter/AliRoot.git] / MUON / AliMUONSt1GeometryBuilderV2.cxx

/**************************************************************************
 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 *                                                                        *
 * Author: The ALICE Off-line Project.                                    *
 * Contributors are mentioned in the code where appropriate.              *
 *                                                                        *
 * Permission to use, copy, modify and distribute this software and its   *
 * documentation strictly for non-commercial purposes is hereby granted   *
 * without fee, provided that the above copyright notice appears in all   *
 * copies and that both the copyright notice and this permission notice   *
 * appear in the supporting documentation. The authors make no claims     *
 * about the suitability of this software for any purpose. It is          *
 * provided "as is" without express or implied warranty.                  *
 **************************************************************************/

// $Id$
//
//-----------------------------------------------------------------------------
// Class AliMUONSt1GeometryBuilderV2
// ---------------------------------
// MUON Station1 detailed geometry construction class.
// (Originally defined in AliMUONv2.cxx - now removed.)
// Included in AliRoot 2004/01/23
// Authors: David Guez, Ivana Hrivnacova, Marion MacCormick; IPN Orsay
//-----------------------------------------------------------------------------

#include "AliMUONSt1GeometryBuilderV2.h"
#include "AliMUONSt1SpecialMotif.h"
#include "AliMUON.h"
#include "AliMUONConstants.h"
#include "AliMUONGeometryModule.h"
#include "AliMUONGeometryEnvelopeStore.h"

#include "AliMpSegmentation.h"
#include "AliMpDEManager.h"
#include "AliMpContainers.h"
#include "AliMpConstants.h"
#include "AliMpCDB.h"
#include "AliMpSectorSegmentation.h"
#include "AliMpSector.h"
#include "AliMpRow.h"
#include "AliMpVRowSegment.h"
#include "AliMpMotifMap.h"
#include "AliMpMotifPosition.h"
#include "AliMpPlaneType.h"

#include "AliRun.h"
#include "AliMagF.h"
#include "AliLog.h"

#include <TVector2.h>
#include <TVector3.h>
#include <TGeoMatrix.h>
#include <TClonesArray.h>
#include <Riostream.h>
#include <TSystem.h>
#include <TVirtualMC.h>
#include <TGeoManager.h>
#include <TGeoVolume.h>
#include <TGeoTube.h>
#include <TGeoXtru.h>
#include <TGeoCompositeShape.h>

#ifdef WITH_STL
  #include <vector>
#endif

#ifdef WITH_ROOT
  #include "TArrayI.h"
#endif

/// \cond CLASSIMP
ClassImp(AliMUONSt1GeometryBuilderV2)
/// \endcond

// Thickness Constants
const GReal_t AliMUONSt1GeometryBuilderV2::fgkHzPadPlane=0.0148/2.;     //Pad plane
const GReal_t AliMUONSt1GeometryBuilderV2::fgkHzFoam = 2.503/2.;        //Foam of mechanicalplane
const GReal_t AliMUONSt1GeometryBuilderV2::fgkHzFR4 = 0.062/2.;         //FR4 of mechanical plane
const GReal_t AliMUONSt1GeometryBuilderV2::fgkHzSnPb = 0.0091/2.;       //Pad/Kapton connection (66 pt)
const GReal_t AliMUONSt1GeometryBuilderV2::fgkHzKapton = 0.0122/2.;     //Kapton
const GReal_t AliMUONSt1GeometryBuilderV2::fgkHzBergPlastic = 0.3062/2.;//Berg connector
const GReal_t AliMUONSt1GeometryBuilderV2::fgkHzBergCopper = 0.1882/2.; //Berg connector
const GReal_t AliMUONSt1GeometryBuilderV2::fgkHzDaughter = 0.0156/2.;   //Daughter board
const GReal_t AliMUONSt1GeometryBuilderV2::fgkHzGas = 0.42/2.;          //Gas thickness

// Quadrant Mother volume - TUBS1 - Middle layer of model
const GReal_t AliMUONSt1GeometryBuilderV2::fgkMotherIR1 = 18.3;
const GReal_t AliMUONSt1GeometryBuilderV2::fgkMotherOR1 = 105.673;   
const GReal_t AliMUONSt1GeometryBuilderV2::fgkMotherThick1 = 6.5/2;  
const GReal_t AliMUONSt1GeometryBuilderV2::fgkMotherPhiL1 = 0.; 
const GReal_t AliMUONSt1GeometryBuilderV2::fgkMotherPhiU1 = 90.;

// Quadrant Mother volume - TUBS2 - near and far layers of model
const GReal_t AliMUONSt1GeometryBuilderV2::fgkMotherIR2 = 20.7;   
const GReal_t AliMUONSt1GeometryBuilderV2::fgkMotherOR2 = 100.073;   
const GReal_t AliMUONSt1GeometryBuilderV2::fgkMotherThick2 = 3.0/2; 
const GReal_t AliMUONSt1GeometryBuilderV2::fgkMotherPhiL2 = 0.; 
const GReal_t AliMUONSt1GeometryBuilderV2::fgkMotherPhiU2 = 90.;

// Sensitive copper pads, foam layer, PCB and electronics model parameters
const GReal_t AliMUONSt1GeometryBuilderV2::fgkHxHole=1.5/2.;
const GReal_t AliMUONSt1GeometryBuilderV2::fgkHyHole=6./2.;
const GReal_t AliMUONSt1GeometryBuilderV2::fgkHxBergPlastic=0.74/2.;
const GReal_t AliMUONSt1GeometryBuilderV2::fgkHyBergPlastic=5.09/2.;
const GReal_t AliMUONSt1GeometryBuilderV2::fgkHxBergCopper=0.25/2.;
const GReal_t AliMUONSt1GeometryBuilderV2::fgkHyBergCopper=3.6/2.;
const GReal_t AliMUONSt1GeometryBuilderV2::fgkHxKapton=0.8/2.;
const GReal_t AliMUONSt1GeometryBuilderV2::fgkHyKapton=5.7/2.;
const GReal_t AliMUONSt1GeometryBuilderV2::fgkHxDaughter=2.3/2.;
const GReal_t AliMUONSt1GeometryBuilderV2::fgkHyDaughter=6.3/2.;
const GReal_t AliMUONSt1GeometryBuilderV2::fgkOffsetX=1.46;
const GReal_t AliMUONSt1GeometryBuilderV2::fgkOffsetY=0.71;
const GReal_t AliMUONSt1GeometryBuilderV2::fgkDeltaFilleEtamX=1.00;
const GReal_t AliMUONSt1GeometryBuilderV2::fgkDeltaFilleEtamY=0.051;

const GReal_t AliMUONSt1GeometryBuilderV2::fgkDeltaQuadLHC=2.6;  // LHC Origin wrt Quadrant Origin
const GReal_t AliMUONSt1GeometryBuilderV2::fgkFrameOffset=5.2;
              // Fix (1) of overlap SQN* layers with SQM* ones (was 5.0)
	      
// Pad planes offsets
const GReal_t AliMUONSt1GeometryBuilderV2::fgkPadXOffsetBP =  0.50 - 0.63/2; // = 0.185
const GReal_t AliMUONSt1GeometryBuilderV2::fgkPadYOffsetBP = -0.31 - 0.42/2; // =-0.52

const char* AliMUONSt1GeometryBuilderV2::fgkHoleName="SCHL";      
const char* AliMUONSt1GeometryBuilderV2::fgkDaughterName="SCDB";  
const char* AliMUONSt1GeometryBuilderV2::fgkQuadrantEnvelopeName="SE";
const char* AliMUONSt1GeometryBuilderV2::fgkQuadrantMLayerName="SQM";
const char* AliMUONSt1GeometryBuilderV2::fgkQuadrantNLayerName="SQN";
const char* AliMUONSt1GeometryBuilderV2::fgkQuadrantFLayerName="SQF";
const char* AliMUONSt1GeometryBuilderV2::fgkQuadrantMFLayerName="SQMF";
const Int_t AliMUONSt1GeometryBuilderV2::fgkFoamBoxNameOffset=200; 
const Int_t AliMUONSt1GeometryBuilderV2::fgkFR4BoxNameOffset=400; 
const Int_t AliMUONSt1GeometryBuilderV2::fgkDaughterCopyNoOffset=1000;

//______________________________________________________________________________
AliMUONSt1GeometryBuilderV2::AliMUONSt1GeometryBuilderV2(AliMUON* muon)
  : AliMUONVGeometryBuilder(0, 2),
    fMUON(muon)
{
/// Standard constructor
}
 
//______________________________________________________________________________
AliMUONSt1GeometryBuilderV2::AliMUONSt1GeometryBuilderV2()
  : AliMUONVGeometryBuilder(),
    fMUON(0)
{
/// Default Constructor
}

//______________________________________________________________________________
AliMUONSt1GeometryBuilderV2::~AliMUONSt1GeometryBuilderV2()
{
/// Destructor
}


//
//  Private methods
//

//______________________________________________________________________________
TString 
AliMUONSt1GeometryBuilderV2::QuadrantEnvelopeName(Int_t chamber, Int_t quadrant) const
{ 
/// Generate unique envelope name from chamber Id and quadrant number

  return Form("%s%d", Form("%s%d",fgkQuadrantEnvelopeName,chamber), quadrant); 
}

//______________________________________________________________________________
void AliMUONSt1GeometryBuilderV2::CreateHole()
{
/// Create all the elements found inside a foam hole

  Int_t* idtmed = fMUON->GetIdtmed()->GetArray()-1099;
  Int_t idAir  = idtmed[1100];      // medium 1
  //Int_t idCopper  = idtmed[1109]; // medium 10 = copper 
  Int_t idCopper  = idtmed[1121]; // medium 22 = copper 

  GReal_t par[3];
  GReal_t posX,posY,posZ;
  
  par[0] = fgkHxHole;
  par[1] = fgkHyHole;
  par[2] = fgkHzFoam;
  gMC->Gsvolu(fgkHoleName,"BOX",idAir,par,3);

  par[0] = fgkHxKapton;
  par[1] = fgkHyKapton;
  par[2] = fgkHzSnPb;
  gMC->Gsvolu("SNPB", "BOX", idCopper, par, 3);
  posX = 0.;
  posY = 0.;
  posZ = -fgkHzFoam+fgkHzSnPb;
  gMC->Gspos("SNPB",1,fgkHoleName, posX, posY, posZ, 0,"ONLY");

  par[0] = fgkHxHole;
  par[1] = fgkHyBergPlastic;
  par[2] = fgkHzKapton;
  gMC->Gsvolu("SKPT", "BOX", idCopper, par, 3);
  posX = 0.;
  posY = 0.;
  posZ = 0.;
  gMC->Gspos("SKPT",1,fgkHoleName, posX, posY, posZ, 0,"ONLY");
}

//______________________________________________________________________________
void AliMUONSt1GeometryBuilderV2::CreateDaughterBoard()
{
/// Create all the elements in a daughter board

  Int_t* idtmed = fMUON->GetIdtmed()->GetArray()-1099;
  Int_t idAir  = idtmed[1100]; // medium 1
  //Int_t idCopper  = idtmed[1109]; // medium 10 = copper
  //Int_t idPlastic  =idtmed[1116]; // medium 17 = Plastic
  Int_t idCopper  = idtmed[1121]; // medium 22 = copper
  Int_t idPlastic  =idtmed[1127]; // medium 28 = Plastic

  GReal_t par[3];
  GReal_t posX,posY,posZ;

  par[0]=fgkHxDaughter;
  par[1]=fgkHyDaughter;
  par[2]=TotalHzDaughter();
  gMC->Gsvolu(fgkDaughterName,"BOX",idAir,par,3);
  
  par[0]=fgkHxBergPlastic;
  par[1]=fgkHyBergPlastic;
  par[2]=fgkHzBergPlastic;
  gMC->Gsvolu("SBGP","BOX",idPlastic,par,3);
  posX=0.;
  posY=0.;
  posZ = -TotalHzDaughter() + fgkHzBergPlastic;
  gMC->Gspos("SBGP",1,fgkDaughterName,posX,posY,posZ,0,"ONLY");

  par[0]=fgkHxBergCopper;
  par[1]=fgkHyBergCopper;
  par[2]=fgkHzBergCopper;
  gMC->Gsvolu("SBGC","BOX",idCopper,par,3);
  posX=0.;
  posY=0.;
  posZ=0.;
  gMC->Gspos("SBGC",1,"SBGP",posX,posY,posZ,0,"ONLY");

  par[0]=fgkHxDaughter;
  par[1]=fgkHyDaughter;
  par[2]=fgkHzDaughter;
  gMC->Gsvolu("SDGH","BOX",idCopper,par,3);
  posX=0.;
  posY=0.;
  posZ = -TotalHzDaughter() + 2.*fgkHzBergPlastic + fgkHzDaughter;
  gMC->Gspos("SDGH",1,fgkDaughterName,posX,posY,posZ,0,"ONLY");
}

//______________________________________________________________________________
void AliMUONSt1GeometryBuilderV2::CreateInnerLayers()
{
/// Create the layer of sensitive volumes with gas
/// and the copper layer.
/// The shape of the sensitive area is defined as an extruded
/// solid substracted with tube (to get inner circular shape). 

  TGeoMedium* kMedArCO2  = gGeoManager->GetMedium("MUON_ARG_CO2");
  TGeoMedium* kMedCopper = gGeoManager->GetMedium("MUON_COPPER_II");

  Double_t rmin = 0.0;
  Double_t rmax = fgkMotherIR1;
  Double_t hz   = fgkHzPadPlane + fgkHzGas;
  new TGeoTube("cutTube",rmin, rmax, hz); 

  Double_t maxXY = 89.0; 
  Double_t xy1   = 77.33;
  Double_t xy2   = 48.77;
  Double_t dxy1  = maxXY - xy1;
    
  Int_t nz = 2;
  Int_t nv = 6;
  Double_t vx[6] = {  0.0,   0.0,   xy2, maxXY, maxXY, dxy1 };
  Double_t vy[6] = { dxy1, maxXY, maxXY,   xy2,   0.0,  0.0 };

  TGeoXtru* xtruS1 = new TGeoXtru(nz);
  xtruS1->SetName("xtruS1");
  xtruS1->DefinePolygon(nv, vx, vy);
  xtruS1->DefineSection(0, -fgkHzGas,  0.0, 0.0, 1.0); 
  xtruS1->DefineSection(1,  fgkHzGas,  0.0, 0.0, 1.0); 
  TGeoCompositeShape* layerS1 = new TGeoCompositeShape("layerS1", "xtruS1-cutTube");
  new TGeoVolume("SA1G", layerS1, kMedArCO2 );
  
  TGeoXtru* xtruS2 = new TGeoXtru(nz);
  xtruS2->SetName("xtruS2");
  xtruS2->DefinePolygon(nv, vx, vy);
  xtruS2->DefineSection(0, -fgkHzGas,  0.0, 0.0, 1.0); 
  xtruS2->DefineSection(1,  fgkHzGas,  0.0, 0.0, 1.0); 
  TGeoCompositeShape* layerS2 = new TGeoCompositeShape("layerS2", "xtruS2-cutTube");
  new TGeoVolume("SA2G", layerS2, kMedArCO2 );

  TGeoXtru* xtruS3 = new TGeoXtru(nz);
  xtruS3->SetName("xtruS3");
  xtruS3->DefinePolygon(nv, vx, vy);
  xtruS3->DefineSection(0, -fgkHzPadPlane,  0.0, 0.0, 1.0); 
  xtruS3->DefineSection(1,  fgkHzPadPlane,  0.0, 0.0, 1.0); 
  TGeoCompositeShape* layerS3 = new TGeoCompositeShape("layerS3", "xtruS3-cutTube");
  new TGeoVolume("SA1C", layerS3, kMedCopper );
}  
  

//______________________________________________________________________________
void AliMUONSt1GeometryBuilderV2::CreateSpacer0()
{
/// The spacer volumes are defined according to the input prepared by Nicole Willis
/// without any modifications
///                                                                       <pre>
/// No.    Type  Material Center (mm)            Dimensions (mm) (half lengths)
///  5     BOX   EPOXY    408.2  430.4 522.41    5.75  1.5   25.5
///  5P    BOX   EPOXY    408.2  445.4 522.41    5.75  1.5   25.5
///  6     BOX   EPOXY    408.2  437.9 519.76    5.75  15.0   1.0
///  6P    BOX   EPOXY    408.2  437.9 525.06    5.75  15.0   1.0
///  7     CYL   INOX     408.2  437.9 522.41    r=3.0  hz=20.63
///                                                                      </pre>

  // tracking medias
  Int_t* idtmed = fMUON->GetIdtmed()->GetArray()-1099;
  Int_t idFrameEpoxy = idtmed[1123]; // medium 24 = Frame Epoxy ME730  // was 20 not 16
  Int_t idInox = idtmed[1128];       // medium 29 Stainless Steel (18%Cr,9%Ni,Fe) // was 21 not 17

  GReal_t par[3];
  par[0] = 0.575;
  par[1] = 0.150;
  par[2] = 2.550;
  gMC->Gsvolu("Spacer05","BOX",idFrameEpoxy,par,3);

  par[0] = 0.575;
  par[1] = 1.500;
  par[2] = 0.100;
  gMC->Gsvolu("Spacer06","BOX",idFrameEpoxy,par,3);

  par[0] = 0.000;
  par[1] = 0.300;
  par[2] = 2.063;
  gMC->Gsvolu("Spacer07","TUBE",idInox,par,3);
}  


//______________________________________________________________________________
void AliMUONSt1GeometryBuilderV2::CreateSpacer()
{
/// The spacer volumes are defined according to the input prepared by Nicole Willis
/// with modifications needed to fit into existing geometry.
///                                                                       <pre>
/// No.    Type  Material Center (mm)            Dimensions (mm) (half lengths)
///  5     BOX   EPOXY    408.2  430.4 522.41    5.75  1.5   25.5
///  5P    BOX   EPOXY    408.2  445.4 522.41    5.75  1.5   25.5
///  6     BOX   EPOXY    408.2  437.9 519.76    5.75  15.0   1.0
///  6P    BOX   EPOXY    408.2  437.9 525.06    5.75  15.0   1.0
///  7     CYL   INOX     408.2  437.9 522.41    r=3.0  hz=20.63
///                                                                      </pre>
/// To fit in existing volumes the volumes 5 and 7 are represented by 2 volumes
/// with half size in z (5A, &A); the dimensions of the volume 5A were also modified
/// to avoid overlaps (x made smaller, y larger to abotain the identical volume)   

  // tracking medias
  Int_t* idtmed = fMUON->GetIdtmed()->GetArray()-1099;
  Int_t idFrameEpoxy = idtmed[1123]; // medium 24 = Frame Epoxy ME730  // was 20 not 16
  Int_t idInox = idtmed[1128];       // medium 29 Stainless Steel (18%Cr,9%Ni,Fe) // was 21 not 17

  //GReal_t par[3];
  //par[0] = 0.575;
  //par[1] = 0.150;
  //par[2] = 2.550;
  //gMC->Gsvolu("Spacer5","BOX",idFrameEpoxy,par,3);

  GReal_t par[3];
  par[0] = 0.510;
  par[1] = 0.170;
  par[2] = 1.1515;
  gMC->Gsvolu("Spacer5A","BOX",idFrameEpoxy,par,3);

  par[0] = 0.510;
  par[1] = 1.500;
  par[2] = 0.100;
  gMC->Gsvolu("Spacer6","BOX",idFrameEpoxy,par,3);

  //par[0] = 0.000;
  //par[1] = 0.300;
  //par[2] = 2.063;
  //gMC->Gsvolu("Spacer7","TUBE",idInox,par,3);

  par[0] = 0.000;
  par[1] = 0.300;
  par[2] = 1.0315;
  gMC->Gsvolu("Spacer7A","TUBE",idInox,par,3);
}  

//______________________________________________________________________________
void AliMUONSt1GeometryBuilderV2::CreateQuadrant(Int_t chamber)
{
/// Create the quadrant (bending and non-bending planes)
/// for the given chamber

  // CreateQuadrantLayersAsVolumes(chamber);
  CreateQuadrantLayersAsAssemblies(chamber);

  CreateFrame(chamber);
  
  TExMap specialMap;
  specialMap.Add(76, (Long_t) new AliMUONSt1SpecialMotif(TVector2( 0.1, 0.72), 90.));
  specialMap.Add(75, (Long_t) new AliMUONSt1SpecialMotif(TVector2( 0.7, 0.36)));
  specialMap.Add(47, (Long_t) new AliMUONSt1SpecialMotif(TVector2(1.01, 0.36)));

  // Load mapping from OCDB
  if ( ! AliMpSegmentation::Instance() ) {
    AliFatal("Mapping has to be loaded first !");
  }
       
  const AliMpSectorSegmentation* kSegmentation1 
    = dynamic_cast<const AliMpSectorSegmentation*>(
        AliMpSegmentation::Instance()
          ->GetMpSegmentation(100, AliMpDEManager::GetCathod(100, AliMp::kBendingPlane))); 
  if ( ! kSegmentation1 ) {
    AliFatal("Could not access sector segmentation !");
  }

  const AliMpSector* kSector1 = kSegmentation1->GetSector();

  //Bool_t reflectZ = true;
  Bool_t reflectZ = false;
  //TVector3 where = TVector3(2.5+0.1+0.56+0.001, 2.5+0.1+0.001, 0.);
  TVector3 where = TVector3(fgkDeltaQuadLHC + fgkPadXOffsetBP, 
                            fgkDeltaQuadLHC + fgkPadYOffsetBP, 0.);
  PlaceSector(kSector1, specialMap, where, reflectZ, chamber);
  
  Int_t nb = AliMpConstants::ManuMask(AliMp::kNonBendingPlane);
  TExMapIter it(&specialMap);
  Long_t key;
  Long_t value;
  while ( it.Next(key,value) == kTRUE ) { 
    delete reinterpret_cast<AliMUONSt1SpecialMotif*>(value);
  }
  specialMap.Delete();
  specialMap.Add(76 | nb,(Long_t) new AliMUONSt1SpecialMotif(TVector2(1.01,0.51),90.));
  specialMap.Add(75 | nb,(Long_t) new AliMUONSt1SpecialMotif(TVector2(2.20,-0.08)));
  specialMap.Add(47 | nb,(Long_t) new AliMUONSt1SpecialMotif(TVector2(2.40,-1.11)));
  specialMap.Add(20 | nb,(Long_t) new AliMUONSt1SpecialMotif(TVector2(0.2 ,-0.08)));
  specialMap.Add(46 | nb,(Long_t) new AliMUONSt1SpecialMotif(TVector2(0.92 , 0.17)));
  specialMap.Add(74 | nb,(Long_t) new AliMUONSt1SpecialMotif(TVector2(0.405, -0.10)));  
      // Fix (7) - overlap of SQ42 with MCHL (after moving the whole sector
      // in the true position)   

  const AliMpSectorSegmentation* kSegmentation2 
    = dynamic_cast<const AliMpSectorSegmentation*>(
        AliMpSegmentation::Instance()
          ->GetMpSegmentation(100, AliMpDEManager::GetCathod(100, AliMp::kNonBendingPlane))); 
  if ( ! kSegmentation2 ) {
    AliFatal("Could not access sector segmentation !");
  }
    
  const AliMpSector* kSector2 = kSegmentation2->GetSector();

  //reflectZ = false;
  reflectZ = true;
  TVector2 offset = kSector2->Position();
  where = TVector3(where.X()+offset.X(), where.Y()+offset.Y(), 0.); 
      // Add the half-pad shift of the non-bending plane wrt bending plane
      // (The shift is defined in the mapping as sector offset)
      // Fix (4) - was TVector3(where.X()+0.63/2, ... - now it is -0.63/2
  PlaceSector(kSector2, specialMap, where, reflectZ, chamber);

  it.Reset();
  while ( it.Next(key,value) == kTRUE ) {
    delete reinterpret_cast<AliMUONSt1SpecialMotif*>(value);
  }
  specialMap.Delete();
}

//______________________________________________________________________________
void AliMUONSt1GeometryBuilderV2::CreateFoamBox(
                                        Int_t segNumber,
                                        const  TVector2& dimensions)
{
/// Create all the elements in the copper plane

  Int_t* idtmed = fMUON->GetIdtmed()->GetArray()-1099;
  Int_t idAir  = idtmed[1100]; // medium 1
  //Int_t idFoam = idtmed[1115]; // medium 16 = Foam
  //Int_t idFR4  = idtmed[1114]; // medium 15 = FR4
  Int_t idFoam = idtmed[1125]; // medium 26 = Foam
  Int_t idFR4  = idtmed[1122]; // medium 23 = FR4

  // mother volume
  GReal_t par[3];
  par[0] = dimensions.X();
  par[1] = dimensions.Y();
  par[2] = TotalHzPlane();
  gMC->Gsvolu(PlaneSegmentName(segNumber).Data(),"BOX",idAir,par,3);
  
  // foam layer
  par[0] = dimensions.X();
  par[1] = dimensions.Y();
  par[2] = fgkHzFoam;
  gMC->Gsvolu(FoamBoxName(segNumber).Data(),"BOX",idFoam,par,3);
  GReal_t posX,posY,posZ;
  posX=0.;
  posY=0.;
  posZ = -TotalHzPlane() + fgkHzFoam;
  gMC->Gspos(FoamBoxName(segNumber).Data(),1, 
             PlaneSegmentName(segNumber).Data(),posX,posY,posZ,0,"ONLY");

  // mechanical plane FR4 layer
  par[0] = dimensions.X();
  par[1] = dimensions.Y();
  par[2] = fgkHzFR4;
  gMC->Gsvolu(FR4BoxName(segNumber).Data(),"BOX",idFR4,par,3);
  posX=0.;
  posY=0.;
  posZ = -TotalHzPlane()+ 2.*fgkHzFoam + fgkHzFR4;
  gMC->Gspos(FR4BoxName(segNumber).Data(),1,
             PlaneSegmentName(segNumber).Data(),posX,posY,posZ,0,"ONLY");
}

//______________________________________________________________________________
void AliMUONSt1GeometryBuilderV2::CreatePlaneSegment(Int_t segNumber,
                                    const  TVector2& dimensions,
      	      	      	      	    Int_t nofHoles)
{
/// Create a segment of a plane (this includes a foam layer, 
/// holes in the foam to feed the kaptons through, kapton connectors
/// and the mother board.)
  
  CreateFoamBox(segNumber,dimensions);
  
  // Place spacer in the concrete plane segments:
  // S225 (in S025), S267 (in S067) in chamber1 and S309 (in S109). S351(in S151) 
  // in chamber2
  // The segments were found as those which caused overlaps when we placed
  // the spacer in global coordinates via PlaceSpacer0 
  //
  //    <posXYZ   X_Y_Z=" 12.6000;   0.75000;   0.0000"> <volume name="Spacer5A"/>
  //    <posXYZ   X_Y_Z=" 12.6000;  -0.75000;   0.0000"> <volume name="Spacer5A"/>
  //    <posXYZ   X_Y_Z=" 12.6000;   0.0000;    1.1515"> <volume name="Spacer6"/>
  //    <posXYZ   X_Y_Z=" 12.6000;   0.0000;    0.0000"> <volume name="Spacer7A"/>

  if ( FoamBoxName(segNumber) == "S225" || 
       FoamBoxName(segNumber) == "S267" ||
       FoamBoxName(segNumber) == "S309" ||
       FoamBoxName(segNumber) == "S351"    )  
  {
    GReal_t posX =  12.6;
    GReal_t posY =  0.75;
    GReal_t posZ = -0.1;
    if ( FoamBoxName(segNumber) == "S267" || 
         FoamBoxName(segNumber) == "S351" ) posY += fgkPadYOffsetBP;    
    gMC->Gspos("Spacer5A", 1, FoamBoxName(segNumber).Data(), posX, posY, posZ,0, "ONLY");

    posY = -0.75;
    if ( FoamBoxName(segNumber) == "S267" || 
         FoamBoxName(segNumber) == "S351" ) posY += fgkPadYOffsetBP;    
    gMC->Gspos("Spacer5A", 2, FoamBoxName(segNumber).Data(), posX, posY, posZ,0, "ONLY");

    posY = 0.0;
    posZ = 1.1515;
    if ( FoamBoxName(segNumber) == "S267" || 
         FoamBoxName(segNumber) == "S351" ) posY += fgkPadYOffsetBP;    
    gMC->Gspos("Spacer6",  1, FoamBoxName(segNumber).Data(), posX, posY, posZ,0, "ONLY");

    posY = 0.0;
    posZ = 0.0;
    if ( FoamBoxName(segNumber) == "S267" || 
         FoamBoxName(segNumber) == "S351" ) posY += fgkPadYOffsetBP;    
    gMC->Gspos("Spacer7A", 1, FoamBoxName(segNumber).Data(), posX, posY, posZ,0, "ONLY");
  }  

  for (Int_t holeNum=0;holeNum<nofHoles;holeNum++) {
    GReal_t posX = ((2.*holeNum+1.)/nofHoles-1.)*dimensions.X();
    GReal_t posY = 0.;
    GReal_t posZ = 0.;
 
    gMC->Gspos(fgkHoleName,holeNum+1,
               FoamBoxName(segNumber).Data(),posX,posY,posZ,0,"ONLY");
  }
}

//______________________________________________________________________________
void AliMUONSt1GeometryBuilderV2::CreateQuadrantLayersAsVolumes(Int_t chamber)
{
/// Create the three main layers as real volumes.
/// Not used anymore.

  // tracking medias
  Int_t* idtmed = fMUON->GetIdtmed()->GetArray()-1099;
  Int_t idAir  = idtmed[1100];       // medium 1

  Float_t par[11];
  Float_t posX,posY,posZ;

// Quadrant volume TUBS1, positioned at the end
  par[0] = fgkMotherIR1;
  par[1] = fgkMotherOR1; 
  par[2] = fgkMotherThick1;  
  par[3] = fgkMotherPhiL1; 
  par[4] = fgkMotherPhiU1;
  gMC->Gsvolu(QuadrantMLayerName(chamber),"TUBS",idAir,par,5);
  // gMC->Gsvolu(QuadrantMFLayerName(chamber),"TUBS",idAir,par,5);

// Replace the volume shape with a composite shape
// with substracted overlap with beam shield (YMOT)

  if ( gMC->IsRootGeometrySupported() ) { 

    // Get shape
    TGeoVolume* mlayer 
      = gGeoManager->FindVolumeFast(QuadrantMLayerName(chamber));
    if ( !mlayer ) {
      AliErrorStream() 
         << "Quadrant volume " << QuadrantMLayerName(chamber) << " not found" 
	 << endl;
    }
    else {
      TGeoShape* quadrant = mlayer->GetShape();
      quadrant->SetName("quadrant");	 

      // Beam shield recess
      par[0] = 0;
      par[1] = 15.4; 
      par[2] = fgkMotherThick1;  
      new TGeoTube("shield_tube", par[0], par[1], par[2]);
  
      // Displacement
      posX = 2.6;
      posY = 2.6;
      posZ = 0;
      TGeoTranslation* displacement 
        = new TGeoTranslation("TR", posX, posY, posZ);
      displacement->RegisterYourself();

      // Composite shape
      TGeoShape* composite
      = new TGeoCompositeShape("composite", "quadrant-shield_tube:TR"); 
      
      // Reset shape to volume      
      mlayer->SetShape(composite);
    }

    TGeoVolume* malayer 
      = gGeoManager->FindVolumeFast(QuadrantMFLayerName(chamber));
    if ( !malayer ) {
      AliErrorStream() 
         << "Quadrant volume " << QuadrantMFLayerName(chamber) << " not found" 
	 << endl;
    }
    else {
      TGeoShape* quadrant = malayer->GetShape();
      quadrant->SetName("quadrant");	 

      // Beam shield recess
      par[0] = 0;
      par[1] = 15.4; 
      par[2] = fgkMotherThick1;  
      new TGeoTube("shield_tube", par[0], par[1], par[2]);
  
      // Displacement
      posX = 2.6;
      posY = 2.6;
      posZ = 0;
      TGeoTranslation* displacement 
        = new TGeoTranslation("TR", posX, posY, posZ);
      displacement->RegisterYourself();

      // Composite shape
      TGeoShape* composite
      = new TGeoCompositeShape("composite", "quadrant-shield_tube:TR"); 
      
      // Reset shape to volume      
      malayer->SetShape(composite);
    }
  }  

// Quadrant volume TUBS2, positioned at the end
  par[0] = fgkMotherIR2;
  par[1] = fgkMotherOR2; 
  par[2] = fgkMotherThick2;  
  par[3] = fgkMotherPhiL2; 
  par[4] = fgkMotherPhiU2;

  gMC->Gsvolu(QuadrantNLayerName(chamber),"TUBS",idAir,par,5); 
  gMC->Gsvolu(QuadrantFLayerName(chamber),"TUBS",idAir,par,5); 
}  

//______________________________________________________________________________
void AliMUONSt1GeometryBuilderV2::CreateQuadrantLayersAsAssemblies(Int_t chamber)
{
/// Create the three main layers as assemblies

  gGeoManager->MakeVolumeAssembly(QuadrantMLayerName(chamber).Data()); 
  gGeoManager->MakeVolumeAssembly(QuadrantMFLayerName(chamber).Data()); 
  gGeoManager->MakeVolumeAssembly(QuadrantNLayerName(chamber).Data()); 
  gGeoManager->MakeVolumeAssembly(QuadrantFLayerName(chamber).Data()); 
}  

//______________________________________________________________________________
void AliMUONSt1GeometryBuilderV2::CreateFrame(Int_t chamber)
{
/// Create the non-sensitive elements of the frame for the \a chamber
///
/// Model and notation:                                                     \n
///                                                                         \n
/// The Quadrant volume name starts with SQ                                 \n
/// The volume segments are numbered 00 to XX                               \n
///                                                                         \n
///                              OutTopFrame                                \n
///                               (SQ02-16)                                 \n 
///                              ------------                               \n
///             OutEdgeFrame   /              |                             \n
///             (SQ17-24)     /               |  InVFrame (SQ00-01)         \n 
///                          /                |                             \n
///                          |                |                             \n 
///               OutVFrame  |            _- -                              \n
///               (SQ25-39)  |           |   InArcFrame (SQ42-45)           \n
///                          |           |                                  \n 
///                          -------------                                  \n 
///                        InHFrame (SQ40-41)                               \n 
///                                                                         \n                         
///                                                                         \n
/// 06 February 2003 - Overlapping volumes resolved.                        \n
/// One quarter chamber is comprised of three TUBS volumes: SQMx, SQNx, and SQFx,
/// where SQMx is the Quadrant Middle layer for chamber \a chamber ( posZ in [-3.25,3.25]),
/// SQNx is the Quadrant Near side layer for chamber \a chamber ( posZ in [-6.25,3-.25) ), and
/// SQFx is the Quadrant Far side layer for chamber \a chamber ( posZ in (3.25,6.25] ).

  // TString quadrantMLayerName = QuadrantMLayerName(chamber);

  TString quadrantMLayerName = QuadrantMFLayerName(chamber);
  TString quadrantNLayerName = QuadrantNLayerName(chamber);
  TString quadrantFLayerName = QuadrantFLayerName(chamber);

  const Float_t kNearFarLHC=2.4;    // Near and Far TUBS Origin wrt LHC Origin

  // tracking medias
  Int_t* idtmed = fMUON->GetIdtmed()->GetArray()-1099;
  
  //Int_t idAir  = idtmed[1100];       // medium 1
  //Int_t idFrameEpoxy = idtmed[1115]; // medium 16 = Frame Epoxy ME730
  //Int_t idInox = idtmed[1116];       // medium 17 Stainless Steel (18%Cr,9%Ni,Fe)
  //Int_t idFR4 = idtmed[1110];        // medium 11 FR4
  //Int_t idCopper = idtmed[1109];     // medium 10 Copper
  //Int_t idAlu = idtmed[1103];        // medium 4 Aluminium
  Int_t idFrameEpoxy = idtmed[1123]; // medium 24 = Frame Epoxy ME730  // was 20 not 16
  Int_t idInox = idtmed[1128];       // medium 29 Stainless Steel (18%Cr,9%Ni,Fe) // was 21 not 17
  Int_t idFR4 = idtmed[1122];        // medium 23 FR4  // was 15 not 11
  Int_t idCopper = idtmed[1121];     // medium 22 Copper
  Int_t idAlu = idtmed[1120];        // medium 21 Aluminium
  
  
    TGeoMedium* kMedEpoxy = gGeoManager->GetMedium("MUON_FrameCH$");
    TGeoMedium* kMedInox  = gGeoManager->GetMedium("MUON_Kapton");
    TGeoMedium* kMedAlu   = gGeoManager->GetMedium("MUON_ALU_II$");


// Rotation Matrices  
      Int_t rot1, rot2, rot3, rot4;    
      
//   Rotation matrices  
     fMUON->AliMatrix(rot1,  90.,  90., 90., 180.,  0., 0.); // +90 deg in x-y plane
     fMUON->AliMatrix(rot2,  90.,  45., 90., 135.,  0., 0.); // +45 deg in x-y plane