[u/mrichter/AliRoot.git] / MUON / AliMUONv2.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 AliMUONv2
// ---------------
// Inherits from AliMUONv1 but with a more detailed
// geometrical description of station 1 

#include <algorithm>
#include <string>

#include <TVector2.h>
#include <TClonesArray.h>
#include <TLorentzVector.h>
#include <TArrayI.h>
#include <Riostream.h>
#include <TSystem.h>

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

#include "AliMUONv2.h"
#include "AliMUONConstants.h"
#include "AliMUONHit.h"
#include "AliRun.h"
#include "AliMagF.h"
#include "AliConst.h" 
#include "AliMUONChamber.h"

ClassImp(AliMUONv2)

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

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

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

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

const GReal_t AliMUONv2::fgkDeltaQuadLHC=2.6;  // LHC Origin wrt Quadrant Origin
const GReal_t AliMUONv2::fgkFrameOffset=5.0;  

const char* AliMUONv2::fgkHoleName="MCHL";      
const char* AliMUONv2::fgkDaughterName="MCDB";  
const char  AliMUONv2::fgkFoamLayerSuffix='F';  // prefix for automatic volume naming
const char* AliMUONv2::fgkQuadrantMLayerName="SQM";
const char* AliMUONv2::fgkQuadrantNLayerName="SQN";
const char* AliMUONv2::fgkQuadrantFLayerName="SQF";

//______________________________________________________________________________
AliMUONv2::AliMUONv2()
  : AliMUONv1()
{
// Default Constructor
// --
   fChamberV2[0] = 0;
   fChamberV2[1] = 0;
   
   // keep secondaries
   SetIshunt(0);
 
   // set path to mapping data files
   if (! gSystem->Getenv("MINSTALL")) {    
     TString dirPath = gSystem->Getenv("ALICE_ROOT");
     dirPath += "/MUON/mapping"; 
     AliMpFiles::Instance()->SetTopPath(dirPath);
     gSystem->Setenv("MINSTALL", dirPath.Data());
     //cout << "AliMpFiles top path set to " << dirPath << endl;	  
   }
   //else
   //  cout << gSystem->Getenv("MINSTALL") << endl;	  
}
 
//______________________________________________________________________________
AliMUONv2::AliMUONv2(const char *name, const char *title)
  : AliMUONv1(name,title)
{
   fChamberV2[0] = 0;
   fChamberV2[1] = 0;
   
   // keep secondaries
   SetIshunt(0);
   
   // set path to mapping data files
   if (! gSystem->Getenv("MINSTALL")) {    
     TString dirPath = gSystem->Getenv("ALICE_ROOT");
     dirPath += "/MUON/mapping"; 
     AliMpFiles::Instance()->SetTopPath(dirPath);
     gSystem->Setenv("MINSTALL", dirPath.Data());
     //cout << "AliMpFiles top path set to " << dirPath << endl;	  
   }
   //else
   //  cout << gSystem->Getenv("MINSTALL") << endl;	    	  
}
 
//______________________________________________________________________________
AliMUONv2::AliMUONv2(const AliMUONv2& rMUON):AliMUONv1(rMUON)
{
// Dummy copy constructor
}

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

//
//  Private methods
//

//______________________________________________________________________________
void AliMUONv2::CreateHole()
{
// Create all the elements found inside a foam hole
// --
  Int_t* idtmed = fIdtmed->GetArray()-1099;
  Int_t idAir  = idtmed[1100];    // medium 1
  Int_t idCopper  = idtmed[1109]; // medium 10 = 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 AliMUONv2::CreateDaughterBoard()
{
// Create all the elements in a daughter board
// --
  Int_t* idtmed = fIdtmed->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

  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 AliMUONv2::CreateInnerLayers()
{
// Create the layer of sensitive volumes with gas
// and the copper layer.
// --

// Gas Medium
  Int_t* idtmed = fIdtmed->GetArray()-1099; 
  Int_t idArCO2  = idtmed[1108];  // medium 9 (ArCO2 80%) 
  Int_t idCopper  = idtmed[1109]; // medium 10 = 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 AliMUONv2::CreateQuadrant(Int_t chamber)
{
// create the quadrant (bending and non-bending planes)
// for the given chamber
// --

  CreateFrame(chamber);

  TSpecialMap specialMap;
  specialMap[1001] = AliMUONSt1SpecialMotif(TVector2( 0.1, 0.84), 90.);
  specialMap[1002] = AliMUONSt1SpecialMotif(TVector2( 0.5, 0.36));
  specialMap[1003] = AliMUONSt1SpecialMotif(TVector2(1.01, 0.36));
  AliMpReader reader1(kStation1, kBendingPlane);
  AliMpSector* sector1 = reader1.BuildSector();

  Bool_t reflectZ = true;
  TVector3 where = TVector3(2.5+0.1+0.56+0.001, 2.5+0.1+0.001, 0.);
  PlaceSector(sector1, specialMap, where, reflectZ, chamber);
  
  specialMap.clear();
  specialMap[4001] = AliMUONSt1SpecialMotif(TVector2(1.01,0.59),90.);
  specialMap[4002] = AliMUONSt1SpecialMotif(TVector2(1.96, 0.17));
  specialMap[4003] = AliMUONSt1SpecialMotif(TVector2(1.61,-1.18));
  specialMap[4004] = AliMUONSt1SpecialMotif(TVector2(0.2 ,-0.08));
  specialMap[4005] = AliMUONSt1SpecialMotif(TVector2(0.2 , 0.25));
  specialMap[4006] = AliMUONSt1SpecialMotif(TVector2(0.28, 0.21));
  AliMpReader reader2(kStation1, kNonBendingPlane);
  AliMpSector* sector2 = reader2.BuildSector();

  reflectZ = false;
  where = TVector3(where.X()+0.63/2.,where.Y()+0.42/2., 0.); //add a half pad shift
  PlaceSector(sector2, specialMap, where, reflectZ, chamber);
}

//______________________________________________________________________________
void AliMUONv2::CreateFoamBox(const char* name,const  TVector2& dimensions)
{
// create all the elements in the copper plane
// --

  Int_t* idtmed = fIdtmed->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

  // 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 AliMUONv2::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 AliMUONv2::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 fgkNearFarLHC=2.4;    // Near and Far TUBS Origin wrt LHC Origin

  // tracking medias
  Int_t* idtmed = fIdtmed->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
  
  
// Rotation Matrices  
      Int_t rot1, rot2, rot3;    
      
//   Rotation matrices  
     AliMatrix(rot1,  90.,  90., 90., 180.,  0., 0.); // +90 deg in x-y plane
     AliMatrix(rot2,  90.,  45., 90., 135.,  0., 0.); // +45 deg in x-y plane
Commit	Line	Data
ba030c0e	1	/**************************************************************************
	2	* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	3	* *
	4	* Author: The ALICE Off-line Project. *
	5	* Contributors are mentioned in the code where appropriate. *
	6	* *
	7	* Permission to use, copy, modify and distribute this software and its *
	8	* documentation strictly for non-commercial purposes is hereby granted *
	9	* without fee, provided that the above copyright notice appears in all *
	10	* copies and that both the copyright notice and this permission notice *
	11	* appear in the supporting documentation. The authors make no claims *
	12	* about the suitability of this software for any purpose. It is *
	13	* provided "as is" without express or implied warranty. *
	14	**************************************************************************/
	15
88cb7938	16	/* $Id$ */
ba030c0e	17
	18	// Authors: David Guez, Ivana Hrivnacova, Marion MacCormick; IPN Orsay
	19	//
	20	// Class AliMUONv2
	21	// ---------------
	22	// Inherits from AliMUONv1 but with a more detailed
	23	// geometrical description of station 1
	24
	25	#include <algorithm>
5f91c9e8	26	#include <string>
ba030c0e	27
	28	#include <TVector2.h>
	29	#include <TClonesArray.h>
	30	#include <TLorentzVector.h>
5f91c9e8	31	#include <TArrayI.h>
	32	#include <Riostream.h>
	33	#include <TSystem.h>
ba030c0e	34
5f91c9e8	35	#include "AliMpFiles.h"
	36	#include "AliMpReader.h"
	37	#include "AliMpSector.h"
	38	#include "AliMpRow.h"
	39	#include "AliMpVRowSegment.h"
	40	#include "AliMpMotifMap.h"
	41	#include "AliMpMotifPosition.h"
ba030c0e	42
	43	#include "AliMUONv2.h"
	44	#include "AliMUONConstants.h"
	45	#include "AliMUONHit.h"
	46	#include "AliRun.h"
	47	#include "AliMagF.h"
	48	#include "AliConst.h"
5f91c9e8	49	#include "AliMUONChamber.h"
ba030c0e	50
	51	ClassImp(AliMUONv2)
	52
5f91c9e8	53	// Thickness Constants
	54	const GReal_t AliMUONv2::fgkHzPadPlane=0.0148/2.; //Pad plane
	55	const GReal_t AliMUONv2::fgkHzFoam = 2.083/2.; //Foam of mechanicalplane
	56	const GReal_t AliMUONv2::fgkHzFR4 = 0.0031/2.; //FR4 of mechanical plane
	57	const GReal_t AliMUONv2::fgkHzSnPb = 0.0091/2.; //Pad/Kapton connection (66 pt)
	58	const GReal_t AliMUONv2::fgkHzKapton = 0.0122/2.; //Kapton
	59	const GReal_t AliMUONv2::fgkHzBergPlastic = 0.3062/2.;//Berg connector
ba030c0e	60	const GReal_t AliMUONv2::fgkHzBergCopper = 0.1882/2.; //Berg connector
5f91c9e8	61	const GReal_t AliMUONv2::fgkHzDaughter = 0.0156/2.; //Daughter board
	62	const GReal_t AliMUONv2::fgkHzGas = 0.2/2.; //Gas thickness
	63
	64	// Quadrant Mother volume - TUBS1 - Middle layer of model
	65	const GReal_t AliMUONv2::fgkMotherIR1 = 18.3;
	66	const GReal_t AliMUONv2::fgkMotherOR1 = 105.673;
	67	const GReal_t AliMUONv2::fgkMotherThick1 = 6.5/2;
	68	const GReal_t AliMUONv2::fgkMotherPhiL1 = 0.;
	69	const GReal_t AliMUONv2::fgkMotherPhiU1 = 90.;
	70
	71	// Quadrant Mother volume - TUBS2 - near and far layers of model
	72	const GReal_t AliMUONv2::fgkMotherIR2 = 20.7;
	73	const GReal_t AliMUONv2::fgkMotherOR2 = 100.073;
	74	const GReal_t AliMUONv2::fgkMotherThick2 = 3.0/2;
	75	const GReal_t AliMUONv2::fgkMotherPhiL2 = 0.;
	76	const GReal_t AliMUONv2::fgkMotherPhiU2 = 90.;
	77
	78	// Sensitive copper pads, foam layer, PCB and electronics model parameters
ba030c0e	79	const GReal_t AliMUONv2::fgkHxHole=1.5/2.;
	80	const GReal_t AliMUONv2::fgkHyHole=6./2.;
	81	const GReal_t AliMUONv2::fgkHxBergPlastic=0.74/2.;
	82	const GReal_t AliMUONv2::fgkHyBergPlastic=5.09/2.;
	83	const GReal_t AliMUONv2::fgkHxBergCopper=0.25/2.;
	84	const GReal_t AliMUONv2::fgkHyBergCopper=3.6/2.;
	85	const GReal_t AliMUONv2::fgkHxKapton=0.8/2.;
	86	const GReal_t AliMUONv2::fgkHyKapton=5.7/2.;
	87	const GReal_t AliMUONv2::fgkHxDaughter=2.3/2.;
	88	const GReal_t AliMUONv2::fgkHyDaughter=6.3/2.;
	89	const GReal_t AliMUONv2::fgkOffsetX=1.46;
	90	const GReal_t AliMUONv2::fgkOffsetY=0.71;
5f91c9e8	91	const GReal_t AliMUONv2::fgkDeltaFilleEtamX=1.46;
5f91c9e8	92	const GReal_t AliMUONv2::fgkDeltaFilleEtamY=0.051;
ba030c0e	93
5f91c9e8	94	const GReal_t AliMUONv2::fgkDeltaQuadLHC=2.6; // LHC Origin wrt Quadrant Origin
5f91c9e8	95	const GReal_t AliMUONv2::fgkFrameOffset=5.0;
ba030c0e	96
5f91c9e8	97	const char* AliMUONv2::fgkHoleName="MCHL";
	98	const char* AliMUONv2::fgkDaughterName="MCDB";
	99	const char AliMUONv2::fgkFoamLayerSuffix='F'; // prefix for automatic volume naming
	100	const char* AliMUONv2::fgkQuadrantMLayerName="SQM";
	101	const char* AliMUONv2::fgkQuadrantNLayerName="SQN";
	102	const char* AliMUONv2::fgkQuadrantFLayerName="SQF";
	103
	104	//______________________________________________________________________________
ba030c0e	105	AliMUONv2::AliMUONv2()
5f91c9e8	106	: AliMUONv1()
ba030c0e	107	{
ba030c0e	108	// Default Constructor
5f91c9e8	109	// --
	110	fChamberV2[0] = 0;
	111	fChamberV2[1] = 0;
	112
ba030c0e	113	// keep secondaries
ba030c0e	114	SetIshunt(0);
5f91c9e8	115
	116	// set path to mapping data files
	117	if (! gSystem->Getenv("MINSTALL")) {
	118	TString dirPath = gSystem->Getenv("ALICE_ROOT");
	119	dirPath += "/MUON/mapping";
	120	AliMpFiles::Instance()->SetTopPath(dirPath);
	121	gSystem->Setenv("MINSTALL", dirPath.Data());
	122	//cout << "AliMpFiles top path set to " << dirPath << endl;
	123	}
	124	//else
	125	// cout << gSystem->Getenv("MINSTALL") << endl;
ba030c0e	126	}
ba030c0e	127
5f91c9e8	128	//______________________________________________________________________________
ba030c0e	129	AliMUONv2::AliMUONv2(const char name, const char title)
5f91c9e8	130	: AliMUONv1(name,title)
ba030c0e	131	{
5f91c9e8	132	fChamberV2[0] = 0;
	133	fChamberV2[1] = 0;
	134
ba030c0e	135	// keep secondaries
ba030c0e	136	SetIshunt(0);
5f91c9e8	137
	138	// set path to mapping data files
	139	if (! gSystem->Getenv("MINSTALL")) {
	140	TString dirPath = gSystem->Getenv("ALICE_ROOT");
	141	dirPath += "/MUON/mapping";
	142	AliMpFiles::Instance()->SetTopPath(dirPath);
	143	gSystem->Setenv("MINSTALL", dirPath.Data());
	144	//cout << "AliMpFiles top path set to " << dirPath << endl;
	145	}
	146	//else
	147	// cout << gSystem->Getenv("MINSTALL") << endl;
ba030c0e	148	}
ba030c0e	149
5f91c9e8	150	//______________________________________________________________________________
6aaf81e6	151	AliMUONv2::AliMUONv2(const AliMUONv2& rMUON):AliMUONv1(rMUON)
ba030c0e	152	{
	153	// Dummy copy constructor
	154	}
	155
5f91c9e8	156	//______________________________________________________________________________
ba030c0e	157	AliMUONv2::~AliMUONv2()
	158	{
	159	// Destructor
ba030c0e	160	}
ba030c0e	161
5f91c9e8	162	//
	163	// Private methods
	164	//
	165
	166	//______________________________________________________________________________
	167	void AliMUONv2::CreateHole()
ba030c0e	168	{
5f91c9e8	169	// Create all the elements found inside a foam hole
ba030c0e	170	// --
ba030c0e	171	Int_t* idtmed = fIdtmed->GetArray()-1099;
5f91c9e8	172	Int_t idAir = idtmed[1100]; // medium 1
5f91c9e8	173	Int_t idCopper = idtmed[1109]; // medium 10 = copper
ba030c0e	174
5f91c9e8	175	GReal_t par[3];
	176	GReal_t posX,posY,posZ;
	177
	178	par[0] = fgkHxHole;
	179	par[1] = fgkHyHole;
	180	par[2] = fgkHzFoam;
	181	gMC->Gsvolu(fgkHoleName,"BOX",idAir,par,3);
ba030c0e	182
5f91c9e8	183	par[0] = fgkHxKapton;
	184	par[1] = fgkHyKapton;
	185	par[2] = fgkHzSnPb;
	186	gMC->Gsvolu("SNPB", "BOX", idCopper, par, 3);
	187	posX = 0.;
	188	posY = 0.;
	189	posZ = -fgkHzFoam+fgkHzSnPb;
	190	gMC->Gspos("SNPB",1,fgkHoleName, posX, posY, posZ, 0,"ONLY");
ba030c0e	191
5f91c9e8	192	par[0] = fgkHxHole;
	193	par[1] = fgkHyBergPlastic;
	194	par[2] = fgkHzKapton;
	195	gMC->Gsvolu("KAPT", "BOX", idCopper, par, 3);
	196	posX = 0.;
	197	posY = 0.;
	198	posZ = 0.;
	199	gMC->Gspos("KAPT",1,fgkHoleName, posX, posY, posZ, 0,"ONLY");
	200	}
ba030c0e	201
5f91c9e8	202	//______________________________________________________________________________
	203	void AliMUONv2::CreateDaughterBoard()
	204	{
	205	// Create all the elements in a daughter board
	206	// --
	207	Int_t* idtmed = fIdtmed->GetArray()-1099;
	208	Int_t idAir = idtmed[1100]; // medium 1
	209	Int_t idCopper = idtmed[1109]; // medium 10 = copper
	210	Int_t idPlastic =idtmed[1116]; // medium 17 = Plastic
ba030c0e	211
5f91c9e8	212	GReal_t par[3];
	213	GReal_t posX,posY,posZ;
	214
	215	par[0]=fgkHxDaughter;
	216	par[1]=fgkHyDaughter;
	217	par[2]=TotalHzDaughter();
	218	gMC->Gsvolu(fgkDaughterName,"BOX",idAir,par,3);
ba030c0e	219
5f91c9e8	220	par[0]=fgkHxBergPlastic;
	221	par[1]=fgkHyBergPlastic;
	222	par[2]=fgkHzBergPlastic;
	223	gMC->Gsvolu("BRGP","BOX",idPlastic,par,3);
	224	posX=0.;
	225	posY=0.;
	226	posZ = -TotalHzDaughter() + fgkHzBergPlastic;
	227	gMC->Gspos("BRGP",1,fgkDaughterName,posX,posY,posZ,0,"ONLY");
	228
	229	par[0]=fgkHxBergCopper;
	230	par[1]=fgkHyBergCopper;
	231	par[2]=fgkHzBergCopper;
	232	gMC->Gsvolu("BRGC","BOX",idCopper,par,3);
	233	posX=0.;
	234	posY=0.;
	235	posZ=0.;
	236	gMC->Gspos("BRGC",1,"BRGP",posX,posY,posZ,0,"ONLY");
	237
	238	par[0]=fgkHxDaughter;
	239	par[1]=fgkHyDaughter;
	240	par[2]=fgkHzDaughter;
	241	gMC->Gsvolu("DGHT","BOX",idCopper,par,3);
	242	posX=0.;
	243	posY=0.;
	244	posZ = -TotalHzDaughter() + 2.*fgkHzBergPlastic + fgkHzDaughter;
	245	gMC->Gspos("DGHT",1,fgkDaughterName,posX,posY,posZ,0,"ONLY");
	246	}
	247
	248	//______________________________________________________________________________
	249	void AliMUONv2::CreateInnerLayers()
	250	{
	251	// Create the layer of sensitive volumes with gas
	252	// and the copper layer.
	253	// --
	254
	255	// Gas Medium
	256	Int_t* idtmed = fIdtmed->GetArray()-1099;
	257	Int_t idArCO2 = idtmed[1108]; // medium 9 (ArCO2 80%)
	258	Int_t idCopper = idtmed[1109]; // medium 10 = copper
	259
	260	Float_t par[11];
	261
	262	//Make gas volume - composed of 11 trapezoids
	263	// section 1 of 11
	264	par[0] = fgkHzGas;
	265	par[1] = 0.;
	266	par[2] = 0.;
	267	par[3] = 71.33/2.;
	268	par[4] = 9.76/2.;
	269	par[5] = 48.77/2.;
	270	par[6] = 15.3;
	271	par[7] = 71.33/2.;
	272	par[8] = 9.76/2.;
	273	par[9] = 48.77/2.;
	274	par[10] = 15.3;
	275
	276	gMC->Gsvolu("SA1G", "TRAP", idArCO2, par, 11);
	277	gMC->Gsvolu("SA2G", "TRAP", idArCO2, par, 11);
ba030c0e	278
5f91c9e8	279	par[0] = fgkHzPadPlane;
	280	gMC->Gsvolu("SA1C", "TRAP", idCopper,par, 11);
	281
	282	// section 2 of 11
	283	par[0] = fgkHzGas;
	284	par[1] = 0.;
	285	par[2] = 0.;
	286	par[3] = 79.68/2.;
	287	par[4] = 10.4/2.;
	288	par[5] = 57.0/2.;
	289	par[6] = 0.;
	290	par[7] = 79.68/2.;
	291	par[8] = 10.4/2.;
	292	par[9] = 57.0/2.;
	293	par[10] = 0.;
	294	gMC->Gsvolu("SB1G", "TRAP", idArCO2, par, 11);
	295	gMC->Gsvolu("SB2G", "TRAP", idArCO2, par, 11);
	296
	297	par[0] = fgkHzPadPlane;
	298	gMC->Gsvolu("SB1C", "TRAP", idCopper,par, 11);
	299
	300	// section 3 of 11
	301	par[0] = fgkHzGas;
	302	par[1] = 0.;
	303	par[2] = 0.;
	304	par[3] = 71.33/2.;
	305	par[4] = 48.77/2.;
	306	par[5] = 9.73/2.;
	307	par[6] = -15.3;
	308	par[7] = 71.33/2.;
	309	par[8] = 48.77/2.;
	310	par[9] = 9.73/2.;
	311	par[10] = -15.3;
	312
	313	gMC->Gsvolu("SC1G", "TRAP", idArCO2, par, 11);
	314	gMC->Gsvolu("SC2G", "TRAP", idArCO2, par, 11);
	315
	316	par[0] = fgkHzPadPlane;
	317	gMC->Gsvolu("SC1C", "TRAP", idCopper,par, 11);
	318
	319	// section 4 of 11
	320	par[0] = fgkHzGas;
	321	par[1] = 0.;
	322	par[2] = 0.;
	323	par[3] = 6.00/2.;
	324	par[4] = 0.;
	325	par[5] = 1.56/2.;
	326	par[6] = 7.41;
	327	par[7] = 6.00/2.;
	328	par[8] = 0.;
	329	par[9] = 1.56/2.;
	330	par[10] = 7.41;
	331	gMC->Gsvolu("SD1G", "TRAP", idArCO2, par, 11);
	332	gMC->Gsvolu("SD2G", "TRAP", idArCO2, par, 11);
	333
	334	par[0] = fgkHzPadPlane;
	335	gMC->Gsvolu("SD1C", "TRAP", idCopper,par, 11);
	336
	337	// section 5 of 11
	338	par[0] = fgkHzGas;
	339	par[1] = 0.;
	340	par[2] = 0.;
	341	par[3] = 1.516/2.;
	342	par[4] = 0.;
343	par[5] = 0.829/2.;
344	par[6] = 15.3;
345	par[7] = 1.516/2.;
346	par[8] = 0.;
347	par[9] = 0.829/2.;
348	par[10] = 15.3;
349	gMC->Gsvolu("SE1G", "TRAP", idArCO2, par, 11);
350	gMC->Gsvolu("SE2G", "TRAP", idArCO2, par, 11);
351
352	par[0] = fgkHzPadPlane;
353	gMC->Gsvolu("SE1C", "TRAP", idCopper,par, 11);
354
355	// section 6 of 11
356	par[0] = fgkHzGas;
357	par[1] = 0.;
358	par[2] = 0.;
359	par[3] = 3.92/2.;
360	par[4] = 0.;
361	par[5] = 0.562/2.;
362	par[6] = -4.1;
363	par[7] = 3.92/2.;
364	par[8] = 0.;
365	par[9] = 0.562/2.;
366	par[10] = -4.1;
367	gMC->Gsvolu("SF1G", "TRAP", idArCO2, par, 11);
368	gMC->Gsvolu("SF2G", "TRAP", idArCO2, par, 11);
369
370	par[0] = fgkHzPadPlane;
371	gMC->Gsvolu("SF1C", "TRAP", idCopper,par, 11);
372
373	// section 7 of 11
374	par[0] = fgkHzGas;
375	par[1] = 0.;
376	par[2] = 0.;
377	par[3] = 0.941/2.;
378	par[4] = 0.562/2.;
379	par[5] = 0.;
380	par[6] = -16.6;
381	par[7] = 0.941/2.;
382	par[8] = 0.562/2.;
383	par[9] = 0.;
384	par[10] =-16.6;
385	gMC->Gsvolu("SG1G", "TRAP", idArCO2, par, 11);
386	gMC->Gsvolu("SG2G", "TRAP", idArCO2, par, 11);
387
388	par[0] = fgkHzPadPlane;
389	gMC->Gsvolu("SG1C", "TRAP", idCopper,par, 11);
390
391	// section 8 of 11
392	par[0] = fgkHzGas;
393	par[1] = 0.;
394	par[2] = 0.;
395	par[3] = 3.94/2.;
396	par[4] = 0.57/2.;
397	par[5] = 0.;
398	par[6] = 4.14;
399	par[7] = 3.94/2.;
400	par[8] = 0.57/2.;
401	par[9] = 0.;
402	par[10] = 4.14;
403	gMC->Gsvolu("SH1G", "TRAP", idArCO2, par, 11);
404	gMC->Gsvolu("SH2G", "TRAP", idArCO2, par, 11);
405
406	par[0] = fgkHzPadPlane;
407	gMC->Gsvolu("SH1C", "TRAP", idCopper,par, 11);
408
409	// section 9 of 11
410	par[0] = fgkHzGas;
411	par[1] = 0.;
412	par[2] = 0.;
413	par[3] = 0.95/2.;
414	par[4] = 0.;
415	par[5] = 0.57/2;
416	par[6] = 16.7;
417	par[7] = 0.95/2.;
418	par[8] = 0.;
419	par[9] = 0.57/2;
420	par[10] = 16.7;
421	gMC->Gsvolu("SI1G", "TRAP", idArCO2, par, 11);
422	gMC->Gsvolu("SI2G", "TRAP", idArCO2, par, 11);
423
424	par[0] = fgkHzPadPlane;
425	gMC->Gsvolu("SI1C", "TRAP", idCopper,par, 11);
426
427	// section 10 of 11
428	par[0] = fgkHzGas;
429	par[1] = 0.;
430	par[2] = 0.;
431	par[3] = 1.49/2.;
432	par[4] = 0.;
433	par[5] = 0.817/2.;
434	par[6] = -15.4;
435	par[7] = 1.49/2.;
436	par[8] = 0.;
437	par[9] = 0.817/2.;
438	par[10] = -15.4;
439	gMC->Gsvolu("SJ1G", "TRAP", idArCO2, par, 11);
440	gMC->Gsvolu("SJ2G", "TRAP", idArCO2, par, 11);
441
442	par[0] = fgkHzPadPlane;
443	gMC->Gsvolu("SJ1C", "TRAP", idCopper,par, 11);
444
445	// section 11 of 11
446	par[0] = fgkHzGas;
447	par[1] = 0.;
448	par[2] = 0.;
449	par[3] = 5.93/2.;
450	par[4] = 0.;
451	par[5] = 1.49/2.;
452	par[6] = -7.16;
453	par[7] = 5.93/2.;
454	par[8] = 0.;
455	par[9] = 1.49/2.;
456	par[10] = -7.16;
457	gMC->Gsvolu("SK1G", "TRAP", idArCO2, par, 11);
458	gMC->Gsvolu("SK2G", "TRAP", idArCO2, par, 11);
459
460	par[0] = fgkHzPadPlane;
461	gMC->Gsvolu("SK1C", "TRAP", idCopper,par, 11);
ba030c0e	462	}
5f91c9e8	463
5f91c9e8	464	//______________________________________________________________________________
ba030c0e	465	void AliMUONv2::CreateQuadrant(Int_t chamber)
	466	{
	467	// create the quadrant (bending and non-bending planes)
	468	// for the given chamber
	469	// --
5f91c9e8	470
ba030c0e	471	CreateFrame(chamber);
ba030c0e	472
5f91c9e8	473	TSpecialMap specialMap;
	474	specialMap[1001] = AliMUONSt1SpecialMotif(TVector2( 0.1, 0.84), 90.);
	475	specialMap[1002] = AliMUONSt1SpecialMotif(TVector2( 0.5, 0.36));
	476	specialMap[1003] = AliMUONSt1SpecialMotif(TVector2(1.01, 0.36));
4a147321	477	AliMpReader reader1(kStation1, kBendingPlane);
5f91c9e8	478	AliMpSector* sector1 = reader1.BuildSector();
ba030c0e	479
5f91c9e8	480	Bool_t reflectZ = true;
	481	TVector3 where = TVector3(2.5+0.1+0.56+0.001, 2.5+0.1+0.001, 0.);
	482	PlaceSector(sector1, specialMap, where, reflectZ, chamber);
ba030c0e	483
ba030c0e	484	specialMap.clear();
	485	specialMap[4001] = AliMUONSt1SpecialMotif(TVector2(1.01,0.59),90.);
	486	specialMap[4002] = AliMUONSt1SpecialMotif(TVector2(1.96, 0.17));
	487	specialMap[4003] = AliMUONSt1SpecialMotif(TVector2(1.61,-1.18));
	488	specialMap[4004] = AliMUONSt1SpecialMotif(TVector2(0.2 ,-0.08));
	489	specialMap[4005] = AliMUONSt1SpecialMotif(TVector2(0.2 , 0.25));
	490	specialMap[4006] = AliMUONSt1SpecialMotif(TVector2(0.28, 0.21));
4a147321	491	AliMpReader reader2(kStation1, kNonBendingPlane);
5f91c9e8	492	AliMpSector* sector2 = reader2.BuildSector();
ba030c0e	493
5f91c9e8	494	reflectZ = false;
	495	where = TVector3(where.X()+0.63/2.,where.Y()+0.42/2., 0.); //add a half pad shift
	496	PlaceSector(sector2, specialMap, where, reflectZ, chamber);
ba030c0e	497	}
ba030c0e	498
5f91c9e8	499	//______________________________________________________________________________
5f91c9e8	500	void AliMUONv2::CreateFoamBox(const char* name,const TVector2& dimensions)
ba030c0e	501	{
5f91c9e8	502	// create all the elements in the copper plane
5f91c9e8	503	// --
ba030c0e	504
5f91c9e8	505	Int_t* idtmed = fIdtmed->GetArray()-1099;
	506	Int_t idAir = idtmed[1100]; // medium 1
	507	Int_t idFoam = idtmed[1115]; // medium 16 = Foam
	508	Int_t idFR4 = idtmed[1114]; // medium 15 = FR4
ba030c0e	509
5f91c9e8	510	// mother volume
	511	GReal_t par[3];
	512	par[0] = dimensions.X();
	513	par[1] = dimensions.Y();
	514	par[2] = TotalHzPlane();
	515	gMC->Gsvolu(name,"BOX",idAir,par,3);
	516
	517	// foam layer
	518	GReal_t posX,posY,posZ;
	519	char eName[5];
	520	strcpy(eName,name);
	521	eName[3]=fgkFoamLayerSuffix;
	522	par[0] = dimensions.X();
	523	par[1] = dimensions.Y();
	524	par[2] = fgkHzFoam;
	525	gMC->Gsvolu(eName,"BOX",idFoam,par,3);
	526	posX=0.;
	527	posY=0.;
	528	posZ = -TotalHzPlane() + fgkHzFoam;
	529	gMC->Gspos(eName,1,name,posX,posY,posZ,0,"ONLY");
ba030c0e	530
5f91c9e8	531	// mechanical plane FR4 layer
	532	eName[3]='R';
	533	par[0] = dimensions.X();
	534	par[1] = dimensions.Y();
	535	par[2] = fgkHzFR4;
	536	gMC->Gsvolu(eName,"BOX",idFR4,par,3);
	537	posX=0.;
	538	posY=0.;
	539	posZ = -TotalHzPlane()+ 2.*fgkHzFoam + fgkHzFR4;
	540	gMC->Gspos(eName,1,name,posX,posY,posZ,0,"ONLY");
	541	}
ba030c0e	542
5f91c9e8	543	//______________________________________________________________________________
	544	void AliMUONv2::CreatePlaneSegment(const char* name,const TVector2& dimensions,
	545	Int_t nofHoles)
	546	{
	547	// Create a segment of a plane (this includes a foam layer,
	548	// holes in the foam to feed the kaptons through, kapton connectors
	549	// and the mother board.)
	550	// --
ba030c0e	551
5f91c9e8	552	CreateFoamBox(name,dimensions);
ba030c0e	553
5f91c9e8	554	char eName[5];
	555	strcpy(eName,name);
	556	eName[3]=fgkFoamLayerSuffix;
ba030c0e	557
5f91c9e8	558	for (Int_t holeNum=0;holeNum<nofHoles;holeNum++) {
	559	GReal_t posX = ((2.holeNum+1.)/nofHoles-1.)dimensions.X();
	560	GReal_t posY = 0.;
	561	GReal_t posZ = 0.;
ba030c0e	562
5f91c9e8	563	gMC->Gspos(fgkHoleName,holeNum+1,eName,posX,posY,posZ,0,"ONLY");
	564	}
	565	}
ba030c0e	566
5f91c9e8	567	//______________________________________________________________________________
	568	void AliMUONv2::CreateFrame(Int_t chamber)
	569	{
	570	// Create the non-sensitive elements of the frame for the <chamber>
	571	//
	572	//
	573	// Model and notation:
	574	//
	575	// The Quadrant volume name starts with SQ
	576	// The volume segments are numbered 00 to XX.
	577	//
	578	// OutTopFrame
	579	// (SQ02-16)
	580	// ------------
	581	// OutEdgeFrame / \|
	582	// (SQ17-24) / \| InVFrame (SQ00-01)
	583	// / \|
	584	// \| \|
	585	// OutVFrame \| _- -
	586	// (SQ25-39) \| \| InArcFrame (SQ42-45)
	587	// \| \|
	588	// -------------
	589	// InHFrame (SQ40-41)
	590	//
	591	//
	592	// 06 February 2003 - Overlapping volumes resolved.
	593	// One quarter chamber is comprised of three TUBS volumes: SQMx, SQNx, and SQFx,
	594	// where SQMx is the Quadrant Middle layer for chamber <x> ( posZ in [-3.25,3.25]),
	595	// SQNx is the Quadrant Near side layer for chamber <x> ( posZ in [-6.25,3-.25) ), and
	596	// SQFx is the Quadrant Far side layer for chamber <x> ( posZ in (3.25,6.25] ).
	597	//---
ba030c0e	598
5f91c9e8	599	const Float_t fgkNearFarLHC=2.4; // Near and Far TUBS Origin wrt LHC Origin
ba030c0e	600
ba030c0e	601	// tracking medias
ba030c0e	602	Int_t* idtmed = fIdtmed->GetArray()-1099;
5f91c9e8	603
ba030c0e	604	Int_t idAir = idtmed[1100]; // medium 1
5f91c9e8	605	Int_t idFrameEpoxy = idtmed[1115]; // medium 16 = Frame Epoxy ME730
	606	Int_t idInox = idtmed[1116]; // medium 17 Stainless Steel (18%Cr,9%Ni,Fe)
	607	Int_t idFR4 = idtmed[1110]; // medium 11 FR4
	608	Int_t idCopper = idtmed[1109]; // medium 10 Copper
	609	Int_t idAlu = idtmed[1103]; // medium 4 Aluminium
	610
	611
	612	// Rotation Matrices
	613	Int_t rot1, rot2, rot3;
	614
	615	// Rotation matrices
	616	AliMatrix(rot1, 90., 90., 90., 180., 0., 0.); // +90 deg in x-y plane
	617	AliMatrix(rot2, 90., 45., 90., 135., 0., 0.); // +45 deg in x-y plane
	618