[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.                  *
 **************************************************************************/

/*
$Log$
*/

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