Improving MUONrectracks (Christophe)

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

#ifdef ST1_WITH_STL
  #include <vector>
#endif

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

#include <TVector2.h>
#include <TVector3.h>
#include <TGeoMatrix.h>
#include <TClonesArray.h>
#include <Riostream.h>
#include <TSystem.h>
#include <TVirtualMC.h>

#include "AliMpFiles.h"
#include "AliMpSectorReader.h"
#include "AliMpSector.h"
#include "AliMpRow.h"
#include "AliMpVRowSegment.h"
#include "AliMpMotifMap.h"
#include "AliMpMotifPosition.h"

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

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

ClassImp(AliMUONSt1GeometryBuilderV2)

// Thickness Constants
const GReal_t AliMUONSt1GeometryBuilderV2::fgkHzPadPlane=0.0148/2.;     //Pad plane
const GReal_t AliMUONSt1GeometryBuilderV2::fgkHzFoam = 2.083/2.;        //Foam of mechanicalplane
const GReal_t AliMUONSt1GeometryBuilderV2::fgkHzFR4 = 0.0031/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.46;
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="MCHL";      
const char* AliMUONSt1GeometryBuilderV2::fgkDaughterName="MCDB";  
const char  AliMUONSt1GeometryBuilderV2::fgkFoamLayerSuffix='F';  // prefix for automatic volume naming
const char* AliMUONSt1GeometryBuilderV2::fgkQuadrantEnvelopeName="SE";
const char* AliMUONSt1GeometryBuilderV2::fgkQuadrantMLayerName="SQM";
const char* AliMUONSt1GeometryBuilderV2::fgkQuadrantNLayerName="SQN";
const char* AliMUONSt1GeometryBuilderV2::fgkQuadrantFLayerName="SQF";
const Int_t AliMUONSt1GeometryBuilderV2::fgkDaughterCopyNoOffset=1000;

//______________________________________________________________________________
AliMUONSt1GeometryBuilderV2::AliMUONSt1GeometryBuilderV2(AliMUON* muon)
  : AliMUONVGeometryBuilder(0, 1),
    fMUON(muon)
{
   // set path to mapping data files
   if (! gSystem->Getenv("MINSTALL")) {    
     TString dirPath = gSystem->Getenv("ALICE_ROOT");
     dirPath += "/MUON/mapping"; 
     AliMpFiles::SetTopPath(dirPath);
     gSystem->Setenv("MINSTALL", dirPath.Data());
     //cout << "AliMpFiles top path set to " << dirPath << endl;	  
   }
   //else
   //  cout << gSystem->Getenv("MINSTALL") << endl;	    	  
}
 
//______________________________________________________________________________
AliMUONSt1GeometryBuilderV2::AliMUONSt1GeometryBuilderV2()
  : AliMUONVGeometryBuilder(),
    fMUON(0)
{
// Default Constructor
// --
}
 
//______________________________________________________________________________
AliMUONSt1GeometryBuilderV2::AliMUONSt1GeometryBuilderV2(const AliMUONSt1GeometryBuilderV2& rhs)
  : AliMUONVGeometryBuilder(rhs)
{
// Dummy copy constructor

 AliFatal("Copy constructor is not implemented.");
}

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


//______________________________________________________________________________
AliMUONSt1GeometryBuilderV2& 
AliMUONSt1GeometryBuilderV2::operator = (const AliMUONSt1GeometryBuilderV2& rhs) 
{
  // check assignement to self
  if (this == &rhs) return *this;

  AliFatal("Assignment operator is not implemented.");
    
  return *this;  
}

//
//  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("KAPT", "BOX", idCopper, par, 3);
  posX = 0.;
  posY = 0.;
  posZ = 0.;
  gMC->Gspos("KAPT",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("BRGP","BOX",idPlastic,par,3);
  posX=0.;
  posY=0.;
  posZ = -TotalHzDaughter() + fgkHzBergPlastic;
  gMC->Gspos("BRGP",1,fgkDaughterName,posX,posY,posZ,0,"ONLY");

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

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

//______________________________________________________________________________
void AliMUONSt1GeometryBuilderV2::CreateInnerLayers()
{
// Create the layer of sensitive volumes with gas
// and the copper layer.
// --

// Gas Medium
  Int_t* idtmed = fMUON->GetIdtmed()->GetArray()-1099; 
  //Int_t idArCO2  = idtmed[1108];  // medium 9 (ArCO2 80%) 
  //Int_t idCopper  = idtmed[1109]; // medium 10 = copper
  Int_t idArCO2   = idtmed[1124]; // medium 25 (ArCO2 80%) 
  Int_t idCopper  = idtmed[1121]; // medium 22 = copper

  Float_t par[11];

//Make gas volume - composed of 11 trapezoids
// section 1 of 11
    par[0] = fgkHzGas;
    par[1] = 0.;
    par[2] = 0.;
    par[3] = 71.33/2.;
    par[4] = 9.76/2.;
    par[5] = 48.77/2.;
    par[6] = 15.3;
    par[7] = 71.33/2.;
    par[8] = 9.76/2.;
    par[9] = 48.77/2.;
    par[10] = 15.3;        

  gMC->Gsvolu("SA1G", "TRAP", idArCO2, par, 11);
  gMC->Gsvolu("SA2G", "TRAP", idArCO2, par, 11);
  
  par[0] = fgkHzPadPlane;
  gMC->Gsvolu("SA1C", "TRAP", idCopper,par, 11);

// section 2 of 11  
    par[0] = fgkHzGas;
    par[1] = 0.;
    par[2] = 0.;
    par[3] = 79.68/2.;
    par[4] = 10.4/2.;
    par[5] = 57.0/2.;
    par[6] = 0.;  
    par[7] = 79.68/2.; 
    par[8] = 10.4/2.;
    par[9] = 57.0/2.;
    par[10] = 0.;  
  gMC->Gsvolu("SB1G", "TRAP", idArCO2, par, 11);  
  gMC->Gsvolu("SB2G", "TRAP", idArCO2, par, 11);

  par[0] = fgkHzPadPlane;
  gMC->Gsvolu("SB1C", "TRAP", idCopper,par, 11);

// section 3 of 11
    par[0] = fgkHzGas;
    par[1] = 0.;
    par[2] = 0.;
    par[3] = 71.33/2.;
    par[4] = 48.77/2.;
    par[5] = 9.73/2.;
    par[6] = -15.3;
    par[7] = 71.33/2.;
    par[8] = 48.77/2.;
    par[9] = 9.73/2.;
    par[10] = -15.3;   
 
  gMC->Gsvolu("SC1G", "TRAP", idArCO2, par, 11);  
  gMC->Gsvolu("SC2G", "TRAP", idArCO2, par, 11);

  par[0] = fgkHzPadPlane;
  gMC->Gsvolu("SC1C", "TRAP", idCopper,par, 11);

// section 4 of 11
    par[0] = fgkHzGas;
    par[1] = 0.;
    par[2] = 0.;
    par[3] = 6.00/2.;
    par[4] = 0.;
    par[5] = 1.56/2.;
    par[6] = 7.41; 
    par[7] = 6.00/2.; 
    par[8] = 0.;
    par[9] = 1.56/2.;
    par[10] = 7.41;    
  gMC->Gsvolu("SD1G", "TRAP", idArCO2, par, 11);  
  gMC->Gsvolu("SD2G", "TRAP", idArCO2, par, 11);

  par[0] = fgkHzPadPlane;
  gMC->Gsvolu("SD1C", "TRAP", idCopper,par, 11);

// section 5 of 11  
    par[0] = fgkHzGas;
    par[1] = 0.;
    par[2] = 0.;
    par[3] = 1.516/2.;
    par[4] = 0.;
    par[5] = 0.829/2.;
    par[6] = 15.3;
    par[7] = 1.516/2.;
    par[8] = 0.;
    par[9] = 0.829/2.;
    par[10] = 15.3;   
  gMC->Gsvolu("SE1G", "TRAP", idArCO2, par, 11);  
  gMC->Gsvolu("SE2G", "TRAP", idArCO2, par, 11);

  par[0] = fgkHzPadPlane;
  gMC->Gsvolu("SE1C", "TRAP", idCopper,par, 11);

// section 6 of 11
    par[0] = fgkHzGas;
    par[1] = 0.;
    par[2] = 0.;
    par[3] = 3.92/2.;
    par[4] = 0.;
    par[5] = 0.562/2.;
    par[6] = -4.1;
    par[7] = 3.92/2.;
    par[8] = 0.;
    par[9] = 0.562/2.;
    par[10] = -4.1;   
  gMC->Gsvolu("SF1G", "TRAP", idArCO2, par, 11);  
  gMC->Gsvolu("SF2G", "TRAP", idArCO2, par, 11);
    
  par[0] = fgkHzPadPlane;
  gMC->Gsvolu("SF1C", "TRAP", idCopper,par, 11);

// section 7 of 11
    par[0] = fgkHzGas;
    par[1] = 0.;
    par[2] = 0.;
    par[3] = 0.941/2.;
    par[4] = 0.562/2.;
    par[5] = 0.;
    par[6] = -16.6; 
    par[7] = 0.941/2.;
    par[8] = 0.562/2.;
    par[9] = 0.;
    par[10] =-16.6;    
  gMC->Gsvolu("SG1G", "TRAP", idArCO2, par, 11);  
  gMC->Gsvolu("SG2G", "TRAP", idArCO2, par, 11);

  par[0] = fgkHzPadPlane;
  gMC->Gsvolu("SG1C", "TRAP", idCopper,par, 11);

// section 8 of 11
    par[0] = fgkHzGas;
    par[1] = 0.;
    par[2] = 0.;
    par[3] = 3.94/2.;
    par[4] = 0.57/2.;
    par[5] = 0.;
    par[6] = 4.14; 
    par[7] = 3.94/2.; 
    par[8] = 0.57/2.;
    par[9] = 0.;
    par[10] = 4.14;    
  gMC->Gsvolu("SH1G", "TRAP", idArCO2, par, 11);  
  gMC->Gsvolu("SH2G", "TRAP", idArCO2, par, 11);

  par[0] = fgkHzPadPlane;
  gMC->Gsvolu("SH1C", "TRAP", idCopper,par, 11);

// section 9 of 11  
    par[0] = fgkHzGas;
    par[1] = 0.;
    par[2] = 0.;
    par[3] = 0.95/2.;
    par[4] = 0.;
    par[5] = 0.57/2;
    par[6] = 16.7;
    par[7] = 0.95/2.;
    par[8] = 0.;
    par[9] = 0.57/2;
    par[10] = 16.7;   
  gMC->Gsvolu("SI1G", "TRAP", idArCO2, par, 11);  
  gMC->Gsvolu("SI2G", "TRAP", idArCO2, par, 11);

  par[0] = fgkHzPadPlane;
  gMC->Gsvolu("SI1C", "TRAP", idCopper,par, 11);

// section 10 of 11
    par[0] = fgkHzGas;
    par[1] = 0.;
    par[2] = 0.;
    par[3] = 1.49/2.;
    par[4] = 0.;
    par[5] = 0.817/2.;
    par[6] = -15.4;
    par[7] = 1.49/2.;
    par[8] = 0.;
    par[9] = 0.817/2.;
    par[10] = -15.4;   
  gMC->Gsvolu("SJ1G", "TRAP", idArCO2, par, 11);  
  gMC->Gsvolu("SJ2G", "TRAP", idArCO2, par, 11);
    
  par[0] = fgkHzPadPlane;
  gMC->Gsvolu("SJ1C", "TRAP", idCopper,par, 11);

// section 11 of 11
    par[0] = fgkHzGas;
    par[1] = 0.;
    par[2] = 0.;
    par[3] = 5.93/2.;
    par[4] = 0.;
    par[5] = 1.49/2.;
    par[6] = -7.16; 
    par[7] = 5.93/2.;
    par[8] = 0.;
    par[9] = 1.49/2.;
    par[10] = -7.16;    
  gMC->Gsvolu("SK1G", "TRAP", idArCO2, par, 11);  
  gMC->Gsvolu("SK2G", "TRAP", idArCO2, par, 11);

  par[0] = fgkHzPadPlane;
  gMC->Gsvolu("SK1C", "TRAP", idCopper,par, 11);
}

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

  CreateFrame(chamber);

#ifdef ST1_WITH_STL
  SpecialMap specialMap;
  specialMap[76] = AliMUONSt1SpecialMotif(TVector2( 0.1, 0.84), 90.);
  specialMap[75] = AliMUONSt1SpecialMotif(TVector2( 0.5, 0.36));
  specialMap[47] = AliMUONSt1SpecialMotif(TVector2(1.01, 0.36));
#endif
  
#ifdef ST1_WITH_ROOT
  SpecialMap specialMap;
  specialMap.Add(76, (Long_t) new AliMUONSt1SpecialMotif(TVector2( 0.1, 0.84), 90.));
  specialMap.Add(75, (Long_t) new AliMUONSt1SpecialMotif(TVector2( 0.5, 0.36)));
  specialMap.Add(47, (Long_t) new AliMUONSt1SpecialMotif(TVector2(1.01, 0.36)));
#endif

  AliMpSectorReader reader1(kStation1, kBendingPlane);
  AliMpSector* sector1 = reader1.BuildSector();

  //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(sector1, specialMap, where, reflectZ, chamber);
  
#ifdef ST1_WITH_STL
  specialMap.clear();
  specialMap[76] = AliMUONSt1SpecialMotif(TVector2(1.01,0.59),90.);
  specialMap[75] = AliMUONSt1SpecialMotif(TVector2(1.96, 0.17));
  specialMap[47] = AliMUONSt1SpecialMotif(TVector2(2.18,-0.98));
  specialMap[20] = AliMUONSt1SpecialMotif(TVector2(0.2 ,-0.08));
  specialMap[46] = AliMUONSt1SpecialMotif(TVector2(0.2 , 0.25));
  specialMap[74] = AliMUONSt1SpecialMotif(TVector2(0.28, 0.21));
      // Fix (7) - overlap of SQ42 with MCHL (after moving the whole sector
      // in the true position)   
      // Was: specialMap[47] = AliMUONSt1SpecialMotif(TVector2(1.61,-1.18));
#endif

#ifdef ST1_WITH_ROOT
  specialMap.Delete();
  specialMap.Add(76,(Long_t) new AliMUONSt1SpecialMotif(TVector2(1.01,0.59),90.));
  specialMap.Add(75,(Long_t) new AliMUONSt1SpecialMotif(TVector2(1.96, 0.17)));
  specialMap.Add(47,(Long_t) new AliMUONSt1SpecialMotif(TVector2(2.18,-0.98)));
  specialMap.Add(20,(Long_t) new AliMUONSt1SpecialMotif(TVector2(0.2 ,-0.08)));
  specialMap.Add(46,(Long_t) new AliMUONSt1SpecialMotif(TVector2(0.2 , 0.25)));
  specialMap.Add(74,(Long_t) new AliMUONSt1SpecialMotif(TVector2(0.28, 0.21)));  
      // Fix (7) - overlap of SQ42 with MCHL (after moving the whole sector
      // in the true position)   
      // Was: specialMap.Add(47,(Long_t) new AliMUONSt1SpecialMotif(TVector2(1.61,-1.18)));
#endif

  AliMpSectorReader reader2(kStation1, kNonBendingPlane);
  AliMpSector* sector2 = reader2.BuildSector();
  
  //reflectZ = false;
  reflectZ = true;
  TVector2 offset = sector2->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(sector2, specialMap, where, reflectZ, chamber);

#ifdef ST1_WITH_ROOT
  specialMap.Delete();
#endif
}

//______________________________________________________________________________
void AliMUONSt1GeometryBuilderV2::CreateFoamBox(const char* name,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(name,"BOX",idAir,par,3);
  
  // foam layer
  GReal_t posX,posY,posZ;
  char eName[5];
  strcpy(eName,name);
  eName[3]=fgkFoamLayerSuffix;
  par[0] = dimensions.X();
  par[1] = dimensions.Y();
  par[2] = fgkHzFoam;
  gMC->Gsvolu(eName,"BOX",idFoam,par,3);
  posX=0.;
  posY=0.;
  posZ = -TotalHzPlane() + fgkHzFoam;
  gMC->Gspos(eName,1,name,posX,posY,posZ,0,"ONLY");

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

//______________________________________________________________________________
void AliMUONSt1GeometryBuilderV2::CreatePlaneSegment(const char* name,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(name,dimensions);

  char eName[5];
  strcpy(eName,name);
  eName[3]=fgkFoamLayerSuffix;
  
  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,eName,posX,posY,posZ,0,"ONLY");
  }
}

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

  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
  
  
// Rotation Matrices  
      Int_t rot1, rot2, rot3;    
      
//   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