//_________________________________________________________________________
// Geometry class for EMCAL : singleton
// EMCAL consists of layers of scintillator and lead
+// with scintillator fiber arranged as "shish-kebab" skewers
// Places the the Barrel Geometry of The EMCAL at Midrapidity
// between 80 and 180(or 190) degrees of Phi and
// -0.7 to 0.7 in eta
-// Number of Modules and Layers may be controlled by
-// the name of the instance defined
+//
// EMCAL geometry tree:
// EMCAL -> superModule -> module -> tower(cell)
// Indexes
// absId -> nSupMod -> nModule -> (nIphi,nIeta)
//
+// Name choices:
+// EMCAL_PDC06 (geometry used for PDC06 simulations, kept for backward compatibility)
+// = equivalent to SHISH_77_TRD1_2X2_FINAL_110DEG in old notation
+// EMCAL_COMPLETE (geometry for expected complete detector)
+// = equivalent to SHISH_77_TRD1_2X2_FINAL_110DEG scTh=0.176 pbTh=0.144
+// in old notation
+// EMCAL_WSUC (Wayne State test stand)
+// = no definite equivalent in old notation, was only used by
+// Aleksei, but kept for testing purposes
+//
+// etc.
+//
+//
+//
//*-- Author: Sahal Yacoob (LBL / UCT)
// and : Yves Schutz (SUBATECH)
// and : Jennifer Klay (LBL)
-// SHASHLYK : Aleksei Pavlinov (WSU)
+// and : Aleksei Pavlinov (WSU)
//
-#include <assert.h>
+#include <cassert>
// --- Root header files ---
#include <Riostream.h>
// these initialisations are needed for a singleton
AliEMCALGeometry *AliEMCALGeometry::fgGeom = 0;
Bool_t AliEMCALGeometry::fgInit = kFALSE;
-Char_t* AliEMCALGeometry::fgDefaultGeometryName = "SHISH_77_TRD1_2X2_FINAL_110DEG scTh=0.176 pbTh=0.144";
+Char_t* AliEMCALGeometry::fgDefaultGeometryName = "EMCAL_COMPLETE";
//
// Usage:
// You can create the AliEMCALGeometry object independently from anything.
// MC: If you work with MC data you have to get geometry the next way:
// == =============================
// AliRunLoader *rl = AliRunLoader::GetRunLoader();
-// AliEMCALGeometry *geom = dynamic_cast<AliEMCAL*>(rl->GetAliRun()->GetDetector("EMCAL"))->GetGeometry();
-// TGeoManager::Import("geometry.root");
+// AliEMCALGeometry *geom = dynamic_cast<AliEMCAL*>(rl->GetAliRun()->GetDetector("EMCAL"))->GetGeometry();
+// TGeoManager::Import("geometry.root");
AliEMCALGeometry::AliEMCALGeometry()
: AliGeometry(),
- fGeoName(0),fArrayOpts(0),fAlFrontThick(0.),fECPbRadThickness(0.),fECScintThick(0.),
+ fGeoName(0),fArrayOpts(0),fNAdditionalOpts(0),fECPbRadThickness(0.),fECScintThick(0.),
fNECLayers(0),fArm1PhiMin(0.),fArm1PhiMax(0.),fArm1EtaMin(0.),fArm1EtaMax(0.),fIPDistance(0.),
- fShellThickness(0.),fZLength(0.),fGap2Active(0.),fNZ(0),fNPhi(0),fSampling(0.),fNumberOfSuperModules(0),
- fSteelFrontThick(0.),fFrontSteelStrip(0.),fLateralSteelStrip(0.),fPassiveScintThick(0.),fPhiModuleSize(0.),
- fEtaModuleSize(0.),fPhiTileSize(0.),fEtaTileSize(0.),fLongModuleSize(0.),fNPhiSuperModule(0),fNPHIdiv(0),fNETAdiv(0),
- fNCells(0),fNCellsInSupMod(0),fNCellsInModule(0),fNTRUEta(0),fNTRUPhi(0),
- fNCellsInTRUEta(0), fNCellsInTRUPhi(0), fTrd1Angle(0.),f2Trd1Dx2(0.),
+ fShellThickness(0.),fZLength(0.),fNZ(0),fNPhi(0),fSampling(0.),fNumberOfSuperModules(0),
+ fFrontSteelStrip(0.),fLateralSteelStrip(0.),fPassiveScintThick(0.),fPhiModuleSize(0.),
+ fEtaModuleSize(0.),fPhiTileSize(0.),fEtaTileSize(0.),fLongModuleSize(0.),fNPhiSuperModule(0),
+ fNPHIdiv(0),fNETAdiv(0), fNCells(0),fNCellsInSupMod(0),fNCellsInModule(0),
+ // Trigger staff
+ fNTRUEta(0),
+ fNTRUPhi(0),
+ fNModulesInTRUEta(0),
+ fNModulesInTRUPhi(0),
+ fNEtaSubOfTRU(0),
+ //
+ fTrd1Angle(0.),f2Trd1Dx2(0.),
fPhiGapForSM(0.),fKey110DEG(0),fPhiBoundariesOfSM(0), fPhiCentersOfSM(0),fEtaMaxOfTRD1(0),
- fTrd2AngleY(0.),f2Trd2Dy2(0.),fEmptySpace(0.),fTubsR(0.),fTubsTurnAngle(0.),fCentersOfCellsEtaDir(0),
- fCentersOfCellsXDir(0),fCentersOfCellsPhiDir(0),fEtaCentersOfCells(0),fPhiCentersOfCells(0),
- fShishKebabTrd1Modules(0), fNAdditionalOpts(0),
- fILOSS(-1), fIHADR(-1)
+ fCentersOfCellsEtaDir(0), fCentersOfCellsXDir(0),fCentersOfCellsPhiDir(0),
+ fEtaCentersOfCells(0),fPhiCentersOfCells(0),fShishKebabTrd1Modules(0),
+ fILOSS(-1), fIHADR(-1),
+ //obsolete member data
+ fAlFrontThick(0.), fGap2Active(0.), fSteelFrontThick(0.), fTrd2AngleY(0.),
+ f2Trd2Dy2(0.), fEmptySpace(0.), fTubsR(0.), fTubsTurnAngle(0.)
{
// default ctor only for internal usage (singleton)
- // must be kept public for root persistency purposes, but should never be called by the outside world
- // CreateListOfTrd1Modules();
+ // must be kept public for root persistency purposes,
+ // but should never be called by the outside world
+
AliDebug(2, "AliEMCALGeometry : default ctor ");
}
//______________________________________________________________________
AliEMCALGeometry::AliEMCALGeometry(const Text_t* name, const Text_t* title)
: AliGeometry(name, title),
- fGeoName(0),fArrayOpts(0),fAlFrontThick(0.),fECPbRadThickness(0.),fECScintThick(0.),
+ fGeoName(0),fArrayOpts(0),fNAdditionalOpts(0),fECPbRadThickness(0.),fECScintThick(0.),
fNECLayers(0),fArm1PhiMin(0.),fArm1PhiMax(0.),fArm1EtaMin(0.),fArm1EtaMax(0.),fIPDistance(0.),
- fShellThickness(0.),fZLength(0.),fGap2Active(0.),fNZ(0),fNPhi(0),fSampling(0.),fNumberOfSuperModules(0),
- fSteelFrontThick(0.),fFrontSteelStrip(0.),fLateralSteelStrip(0.),fPassiveScintThick(0.),fPhiModuleSize(0.),
- fEtaModuleSize(0.),fPhiTileSize(0.),fEtaTileSize(0.),fLongModuleSize(0.),fNPhiSuperModule(0),fNPHIdiv(0),fNETAdiv(0),
- fNCells(0),fNCellsInSupMod(0),fNCellsInModule(0),fNTRUEta(0),fNTRUPhi(0),
- fNCellsInTRUEta(0), fNCellsInTRUPhi(0), fTrd1Angle(0.),f2Trd1Dx2(0.),
+ fShellThickness(0.),fZLength(0.),fNZ(0),fNPhi(0),fSampling(0.),fNumberOfSuperModules(0),
+ fFrontSteelStrip(0.),fLateralSteelStrip(0.),fPassiveScintThick(0.),fPhiModuleSize(0.),
+ fEtaModuleSize(0.),fPhiTileSize(0.),fEtaTileSize(0.),fLongModuleSize(0.),fNPhiSuperModule(0),
+ fNPHIdiv(0),fNETAdiv(0), fNCells(0),fNCellsInSupMod(0),fNCellsInModule(0),
+ // Trigger staff
+ fNTRUEta(0),
+ fNTRUPhi(0),
+ fNModulesInTRUEta(0),
+ fNModulesInTRUPhi(0),
+ fNEtaSubOfTRU(0),
+ //
+ fTrd1Angle(0.),f2Trd1Dx2(0.),
fPhiGapForSM(0.),fKey110DEG(0),fPhiBoundariesOfSM(0), fPhiCentersOfSM(0), fEtaMaxOfTRD1(0),
- fTrd2AngleY(0.),f2Trd2Dy2(0.),fEmptySpace(0.),fTubsR(0.),fTubsTurnAngle(0.),fCentersOfCellsEtaDir(0),
- fCentersOfCellsXDir(0),fCentersOfCellsPhiDir(0),fEtaCentersOfCells(0),fPhiCentersOfCells(0),
- fShishKebabTrd1Modules(0),fNAdditionalOpts(0),
- fILOSS(-1), fIHADR(-1)
+ fCentersOfCellsEtaDir(0),fCentersOfCellsXDir(0),fCentersOfCellsPhiDir(0),
+ fEtaCentersOfCells(0),fPhiCentersOfCells(0),fShishKebabTrd1Modules(0),
+ fILOSS(-1), fIHADR(-1),
+ //obsolete member data
+ fAlFrontThick(0.), fGap2Active(0.), fSteelFrontThick(0.), fTrd2AngleY(0.),
+ f2Trd2Dy2(0.), fEmptySpace(0.), fTubsR(0.), fTubsTurnAngle(0.)
{
// ctor only for internal usage (singleton)
AliDebug(2, Form("AliEMCALGeometry(%s,%s) ", name,title));
: AliGeometry(geom),
fGeoName(geom.fGeoName),
fArrayOpts(geom.fArrayOpts),
- fAlFrontThick(geom.fAlFrontThick),
+ fNAdditionalOpts(geom.fNAdditionalOpts),
fECPbRadThickness(geom.fECPbRadThickness),
fECScintThick(geom.fECScintThick),
fNECLayers(geom.fNECLayers),
fIPDistance(geom.fIPDistance),
fShellThickness(geom.fShellThickness),
fZLength(geom.fZLength),
- fGap2Active(geom.fGap2Active),
fNZ(geom.fNZ),
fNPhi(geom.fNPhi),
fSampling(geom.fSampling),
fNumberOfSuperModules(geom.fNumberOfSuperModules),
- fSteelFrontThick(geom.fSteelFrontThick),
fFrontSteelStrip(geom.fFrontSteelStrip),
fLateralSteelStrip(geom.fLateralSteelStrip),
fPassiveScintThick(geom.fPassiveScintThick),
fNCells(geom.fNCells),
fNCellsInSupMod(geom.fNCellsInSupMod),
fNCellsInModule(geom.fNCellsInModule),
+ // Trigger staff
fNTRUEta(geom.fNTRUEta),
fNTRUPhi(geom.fNTRUPhi),
- fNCellsInTRUEta(geom.fNCellsInTRUEta),
- fNCellsInTRUPhi(geom.fNCellsInTRUPhi),
+ fNModulesInTRUEta(geom.fNModulesInTRUEta),
+ fNModulesInTRUPhi(geom.fNModulesInTRUPhi),
+ fNEtaSubOfTRU(geom.fNEtaSubOfTRU),
+ //
fTrd1Angle(geom.fTrd1Angle),
f2Trd1Dx2(geom.f2Trd1Dx2),
fPhiGapForSM(geom.fPhiGapForSM),
fPhiBoundariesOfSM(geom.fPhiBoundariesOfSM),
fPhiCentersOfSM(geom.fPhiCentersOfSM),
fEtaMaxOfTRD1(geom.fEtaMaxOfTRD1),
- fTrd2AngleY(geom.fTrd2AngleY),
- f2Trd2Dy2(geom.f2Trd2Dy2),
- fEmptySpace(geom.fEmptySpace),
- fTubsR(geom.fTubsR),
- fTubsTurnAngle(geom.fTubsTurnAngle),
fCentersOfCellsEtaDir(geom.fCentersOfCellsEtaDir),
fCentersOfCellsXDir(geom.fCentersOfCellsXDir),
fCentersOfCellsPhiDir(geom.fCentersOfCellsPhiDir),
fEtaCentersOfCells(geom.fEtaCentersOfCells),
fPhiCentersOfCells(geom.fPhiCentersOfCells),
fShishKebabTrd1Modules(geom.fShishKebabTrd1Modules),
- fNAdditionalOpts(geom.fNAdditionalOpts),
- fILOSS(geom.fILOSS), fIHADR(geom.fIHADR)
+ fILOSS(geom.fILOSS), fIHADR(geom.fIHADR),
+ //obsolete member data
+ fAlFrontThick(geom.fAlFrontThick),
+ fGap2Active(geom.fGap2Active),
+ fSteelFrontThick(geom.fSteelFrontThick),
+ fTrd2AngleY(geom.fTrd2AngleY),
+ f2Trd2Dy2(geom.f2Trd2Dy2),
+ fEmptySpace(geom.fEmptySpace),
+ fTubsR(geom.fTubsR),
+ fTubsTurnAngle(geom.fTubsTurnAngle)
{
//copy ctor
}
AliEMCALGeometry::~AliEMCALGeometry(void){
// dtor
}
+
//______________________________________________________________________
void AliEMCALGeometry::Init(void){
- // Initializes the EMCAL parameters
- // naming convention : GUV_WX_N_ gives the composition of a tower
- // WX inform about the composition of the EM calorimeter section:
- // thickness in mm of Pb radiator (W) and of scintillator (X), and number of scintillator layers (N)
- // New geometry: EMCAL_55_25
- // 24-aug-04 for shish-kebab
- // SHISH_25 or SHISH_62
- // 11-oct-05 - correction for pre final design
- // Feb 06,2006 - decrease the weight of EMCAL
//
- // Oct 30,2006 - SHISH_TRD1_CURRENT_1X1, SHISH_TRD1_CURRENT_2X2 or SHISH_TRD1_CURRENT_3X3;
+ // Initializes the EMCAL parameters based on the name
+ // Only Shashlyk geometry is available, but various combinations of
+ // layers and number of supermodules can be selected with additional
+ // options or geometry name
//
fAdditionalOpts[0] = "nl="; // number of sampling layers (fNECLayers)
fNAdditionalOpts = sizeof(fAdditionalOpts) / sizeof(char*);
+ // geometry
fgInit = kFALSE; // Assume failed until proven otherwise.
fGeoName = GetName();
fGeoName.ToUpper();
- fKey110DEG = 0;
- if(fGeoName.Contains("110DEG") || fGeoName.Contains("CURRENT")) fKey110DEG = 1; // for GetAbsCellId
- fShishKebabTrd1Modules = 0;
- fTrd2AngleY = f2Trd2Dy2 = fEmptySpace = fTubsR = fTubsTurnAngle = 0;
-
- fNZ = 114; // granularity along Z (eta)
- fNPhi = 168; // granularity in phi (azimuth)
- fArm1PhiMin = 80.0; // degrees, Starting EMCAL Phi position
- fArm1PhiMax = 190.0; // degrees, Ending EMCAL Phi position
- fArm1EtaMin = -0.7; // pseudorapidity, Starting EMCAL Eta position
- fArm1EtaMax = +0.7; // pseudorapidity, Ending EMCAL Eta position
- fIPDistance = 454.0; // cm, Radial distance to inner surface of EMCAL
- fPhiGapForSM = 0.; // cm, only for final TRD1 geometry
- // geometry
- if(fGeoName.Contains("SHISH")){ // Only shahslyk now
- // 7-sep-05; integration issue
- fArm1PhiMin = 80.0; // 60 -> 80
- fArm1PhiMax = 180.0; // 180 -> 190
-
- fNumberOfSuperModules = 10; // 12 = 6 * 2 (6 in phi, 2 in Z);
- fSteelFrontThick = 2.54; // 9-sep-04
- fIPDistance = 460.0;
- fFrontSteelStrip = fPassiveScintThick = 0.0; // 13-may-05
- fLateralSteelStrip = 0.025; // before MAY 2005
- fPhiModuleSize = fEtaModuleSize = 11.4;
- fPhiTileSize = fEtaTileSize = 5.52; // (11.4-5.52*2)/2. = 0.18 cm (wall thickness)
- fNPhi = 14;
- fNZ = 30;
- fAlFrontThick = fGap2Active = 0;
- fNPHIdiv = fNETAdiv = 2;
-
- fNECLayers = 62;
- fECScintThick = fECPbRadThickness = 0.2;
- fSampling = 1.; // 30-aug-04 - should be calculated
- if(fGeoName.Contains("TWIST")) { // all about EMCAL module
- fNZ = 27; // 16-sep-04
- } else if(fGeoName.Contains("TRD")) {
- fIPDistance = 428.0; // 11-may-05
- fSteelFrontThick = 0.0; // 3.17 -> 0.0; 28-mar-05 : no stell plate
- fNPhi = 12;
- fSampling = 12.327;
- fPhiModuleSize = fEtaModuleSize = 12.26;
- fNZ = 26; // 11-oct-04
- fTrd1Angle = 1.3; // in degree
-// 18-nov-04; 1./0.08112=12.327
-// http://pdsfweb01.nersc.gov/~pavlinov/ALICE/SHISHKEBAB/RES/linearityAndResolutionForTRD1.html
- if(fGeoName.Contains("TRD1")) { // 30-jan-05
- // for final design
- fPhiGapForSM = 2.; // cm, only for final TRD1 geometry
- if(fGeoName.Contains("MAY05") || fGeoName.Contains("WSUC") || fGeoName.Contains("FINAL") || fGeoName.Contains("CURRENT")){
- fNumberOfSuperModules = 12; // 20-may-05
- if(fGeoName.Contains("WSUC")) fNumberOfSuperModules = 1; // 27-may-05
- fNECLayers = 77; // (13-may-05 from V.Petrov)
- fPhiModuleSize = 12.5; // 20-may-05 - rectangular shape
- fEtaModuleSize = 11.9;
- fECScintThick = fECPbRadThickness = 0.16;// (13-may-05 from V.Petrov)
- fFrontSteelStrip = 0.025;// 0.025cm = 0.25mm (13-may-05 from V.Petrov)
- fLateralSteelStrip = 0.01; // 0.01cm = 0.1mm (13-may-05 from V.Petrov) - was 0.025
- fPassiveScintThick = 0.8; // 0.8cm = 8mm (13-may-05 from V.Petrov)
- fNZ = 24;
- fTrd1Angle = 1.5; // 1.3 or 1.5
-
- if(fGeoName.Contains("FINAL") || fGeoName.Contains("CURRENT")) { // 9-sep-05
- fNumberOfSuperModules = 10;
- if(GetKey110DEG()) {
- fNumberOfSuperModules = 12;// last two modules have size 10 degree in phi (180<phi<190)
- fArm1PhiMax = 200.0; // for XEN1 and turn angle of super modules
- }
- if(fGeoName.Contains("FINAL")) {
- fPhiModuleSize = 12.26 - fPhiGapForSM / Float_t(fNPhi); // first assumption
- } else if(fGeoName.Contains("CURRENT")) {
- fECScintThick = 0.176; // 10% of weight reduction
- fECPbRadThickness = 0.144; //
- fLateralSteelStrip = 0.015; // 0.015cm = 0.15mm (Oct 30, from Fred)
- fPhiModuleSize = 12.00;
- fPhiGapForSM = (12.26 - fPhiModuleSize)*fNPhi; // have to check
- }
- fEtaModuleSize = fPhiModuleSize;
- if(fGeoName.Contains("HUGE")) fNECLayers *= 3; // 28-oct-05 for analysing leakage
- }
- }
- } else if(fGeoName.Contains("TRD2")) { // 30-jan-05
- fSteelFrontThick = 0.0; // 11-mar-05
- fIPDistance+= fSteelFrontThick; // 1-feb-05 - compensate absence of steel plate
- fTrd1Angle = 1.64; // 1.3->1.64
- fTrd2AngleY = fTrd1Angle; // symmetric case now
- fEmptySpace = 0.2; // 2 mm
- fTubsR = fIPDistance; // 31-jan-05 - as for Fred case
-
- fPhiModuleSize = fTubsR*2.*TMath::Tan(fTrd2AngleY*TMath::DegToRad()/2.);
- fPhiModuleSize -= fEmptySpace/2.; // 11-mar-05
- fEtaModuleSize = fPhiModuleSize; // 20-may-05
- fTubsTurnAngle = 3.;
- }
- fNPHIdiv = fNETAdiv = 2; // 13-oct-04 - division again
- if(fGeoName.Contains("3X3")) { // 23-nov-04
- fNPHIdiv = fNETAdiv = 3;
- } else if(fGeoName.Contains("4X4")) {
- fNPHIdiv = fNETAdiv = 4;
- } else if(fGeoName.Contains("1X1")) {
- fNPHIdiv = fNETAdiv = 1;
- }
- }
- if(fGeoName.Contains("25")){
- fNECLayers = 25;
- fECScintThick = fECPbRadThickness = 0.5;
- }
- if(fGeoName.Contains("WSUC")){ // 18-may-05 - about common structure
- fShellThickness = 30.; // should be change
- fNPhi = fNZ = 4;
+ //Convert old geometry names to new ones
+ if(fGeoName.Contains("SHISH_77_TRD1_2X2_FINAL_110DEG")) {
+ if(fGeoName.Contains("PBTH=0.144") && fGeoName.Contains("SCTH=0.176")) {
+ fGeoName = "EMCAL_COMPLETE";
+ } else {
+ fGeoName = "EMCAL_PDC06";
}
+ }
+ if(fGeoName.Contains("WSUC")) fGeoName = "EMCAL_WSUC";
- CheckAdditionalOptions();
- DefineSamplingFraction();
+ //check that we have a valid geometry name
+ if(!(fGeoName.Contains("EMCAL_PDC06") || fGeoName.Contains("EMCAL_COMPLETE") || fGeoName.Contains("EMCAL_WSUC"))) {
+ Fatal("Init", "%s is an undefined geometry!", fGeoName.Data()) ;
+ }
- fPhiTileSize = fPhiModuleSize/double(fNPHIdiv) - fLateralSteelStrip; // 13-may-05
- fEtaTileSize = fEtaModuleSize/double(fNETAdiv) - fLateralSteelStrip; // 13-may-05
+ // Option to know whether we have the "half" supermodule(s) or not
+ fKey110DEG = 0;
+ if(fGeoName.Contains("COMPLETE") || fGeoName.Contains("PDC06")) fKey110DEG = 1; // for GetAbsCellId
+ fShishKebabTrd1Modules = 0;
- // constant for transition absid <--> indexes
- fNCellsInModule = fNPHIdiv*fNETAdiv;
- fNCellsInSupMod = fNCellsInModule*fNPhi*fNZ;
- fNCells = fNCellsInSupMod*fNumberOfSuperModules;
- if(GetKey110DEG()) fNCells -= fNCellsInSupMod;
+ // JLK 13-Apr-2008
+ //default parameters are those of EMCAL_COMPLETE geometry
+ //all others render variations from these at the end of
+ //geometry-name specific options
+
+ fNumberOfSuperModules = 12; // 12 = 6 * 2 (6 in phi, 2 in Z)
+ fNPhi = 12; // module granularity in phi within smod (azimuth)
+ fNZ = 24; // module granularity along Z within smod (eta)
+ fNPHIdiv = fNETAdiv = 2; // tower granularity within module
+ fArm1PhiMin = 80.0; // degrees, Starting EMCAL Phi position
+ fArm1PhiMax = 200.0; // degrees, Ending EMCAL Phi position
+ fArm1EtaMin = -0.7; // pseudorapidity, Starting EMCAL Eta position
+ fArm1EtaMax = +0.7; // pseudorapidity, Ending EMCAL Eta position
+ fIPDistance = 428.0; // cm, radial distance to front face from nominal vertex point
+ fPhiGapForSM = 2.; // cm, only for final TRD1 geometry
+ fFrontSteelStrip = 0.025; // 0.025cm = 0.25mm (13-may-05 from V.Petrov)
+ fPassiveScintThick = 0.8; // 0.8cm = 8mm (13-may-05 from V.Petrov)
+ fLateralSteelStrip = 0.01; // 0.01cm = 0.1mm (13-may-05 from V.Petrov) - was 0.025
+ fTrd1Angle = 1.5; // in degrees
+
+ fSampling = 1.; // should be calculated with call to DefineSamplingFraction()
+ fNECLayers = 77; // (13-may-05 from V.Petrov) - can be changed with additional options
+ fECScintThick = 0.176; // scintillator layer thickness
+ fECPbRadThickness = 0.144; // lead layer thickness
+
+ fPhiModuleSize = 12.26 - fPhiGapForSM / Float_t(fNPhi); // first assumption
+ fEtaModuleSize = fPhiModuleSize;
+
+ //needs to be called for each geometry and before setting geometry
+ //parameters which can depend on the outcome
+ CheckAdditionalOptions();
+
+ //modifications to the above for PDC06 geometry
+ if(fGeoName.Contains("PDC06")){ // 18-may-05 - about common structure
+ fECScintThick = fECPbRadThickness = 0.16;// (13-may-05 from V.Petrov)
+ CheckAdditionalOptions();
+ }
- fLongModuleSize = fNECLayers*(fECScintThick + fECPbRadThickness);
- if(fGeoName.Contains("MAY05")) fLongModuleSize += (fFrontSteelStrip + fPassiveScintThick);
+ //modifications to the above for WSUC geometry
+ if(fGeoName.Contains("WSUC")){ // 18-may-05 - about common structure
+ fPhiModuleSize = 12.5; // 20-may-05 - rectangular shape
+ fEtaModuleSize = 11.9;
+ fECScintThick = fECPbRadThickness = 0.16;// (13-may-05 from V.Petrov)
+ fNumberOfSuperModules = 1; // 27-may-05
+ fShellThickness = 30.; // should be change
+ fNPhi = fNZ = 4;
+ CheckAdditionalOptions();
+ }
- // 30-sep-04
- if(fGeoName.Contains("TRD")) {
- f2Trd1Dx2 = fEtaModuleSize + 2.*fLongModuleSize*TMath::Tan(fTrd1Angle*TMath::DegToRad()/2.);
- if(fGeoName.Contains("TRD2")) { // 27-jan-05
- f2Trd2Dy2 = fPhiModuleSize + 2.*fLongModuleSize*TMath::Tan(fTrd2AngleY*TMath::DegToRad()/2.);
- }
- }
- } else Fatal("Init", "%s is an undefined geometry!", fGeoName.Data()) ;
+ // constant for transition absid <--> indexes
+ fNCellsInModule = fNPHIdiv*fNETAdiv;
+ fNCellsInSupMod = fNCellsInModule*fNPhi*fNZ;
+ fNCells = fNCellsInSupMod*fNumberOfSuperModules;
+ if(GetKey110DEG()) fNCells -= fNCellsInSupMod;
fNPhiSuperModule = fNumberOfSuperModules/2;
if(fNPhiSuperModule<1) fNPhiSuperModule = 1;
+
+ fPhiTileSize = fPhiModuleSize/double(fNPHIdiv) - fLateralSteelStrip; // 13-may-05
+ fEtaTileSize = fEtaModuleSize/double(fNETAdiv) - fLateralSteelStrip; // 13-may-05
- fShellThickness = fAlFrontThick + fGap2Active + fNECLayers*GetECScintThick()+(fNECLayers-1)*GetECPbRadThick();
- if(fGeoName.Contains("SHISH")) {
- fShellThickness = fSteelFrontThick + fLongModuleSize;
- if(fGeoName.Contains("TWIST")) { // 13-sep-04
- fShellThickness = TMath::Sqrt(fLongModuleSize*fLongModuleSize + fPhiModuleSize*fEtaModuleSize);
- fShellThickness += fSteelFrontThick;
- } else if(fGeoName.Contains("TRD")) { // 1-oct-04
- fShellThickness = TMath::Sqrt(fLongModuleSize*fLongModuleSize + f2Trd1Dx2*f2Trd1Dx2);
- fShellThickness += fSteelFrontThick;
- // Local coordinates
- fParSM[0] = GetShellThickness()/2.;
- fParSM[1] = GetPhiModuleSize() * GetNPhi()/2.;
- fParSM[2] = 350./2.;
- }
- }
-
+ fLongModuleSize = fNECLayers*(fECScintThick + fECPbRadThickness);
+ f2Trd1Dx2 = fEtaModuleSize + 2.*fLongModuleSize*TMath::Tan(fTrd1Angle*TMath::DegToRad()/2.);
+ if(!fGeoName.Contains("WSUC")) fShellThickness = TMath::Sqrt(fLongModuleSize*fLongModuleSize + f2Trd1Dx2*f2Trd1Dx2);
+
fZLength = 2.*ZFromEtaR(fIPDistance+fShellThickness,fArm1EtaMax); // Z coverage
fEnvelop[0] = fIPDistance; // mother volume inner radius
fEnvelop[1] = fIPDistance + fShellThickness; // mother volume outer r.
fEnvelop[2] = 1.00001*fZLength; // add some padding for mother volume.
- fNumberOfSuperModules = 12;
+ // Local coordinates
+ fParSM[0] = GetShellThickness()/2.;
+ fParSM[1] = GetPhiModuleSize() * GetNPhi()/2.;
+ fParSM[2] = 350./2.;
// SM phi boundaries - (0,1),(2,3) .. (10,11) - has the same boundaries; Nov 7, 2006
fPhiBoundariesOfSM.Set(fNumberOfSuperModules);
fPhiCentersOfSM.Set(fNumberOfSuperModules/2);
fPhiBoundariesOfSM[0] = TMath::PiOver2() - TMath::ATan2(fParSM[1] , fIPDistance); // 1th and 2th modules)
- fPhiBoundariesOfSM[1] = TMath::PiOver2() + TMath::ATan2(fParSM[1] , fIPDistance);
fPhiCentersOfSM[0] = TMath::PiOver2();
- for(int i=1; i<=4; i++) { // from 2th ro 9th
- fPhiBoundariesOfSM[2*i] = fPhiBoundariesOfSM[0] + 20.*TMath::DegToRad()*i;
- fPhiBoundariesOfSM[2*i+1] = fPhiBoundariesOfSM[1] + 20.*TMath::DegToRad()*i;
- fPhiCentersOfSM[i] = fPhiCentersOfSM[0] + 20.*TMath::DegToRad()*i;
+ if(fNumberOfSuperModules > 1)
+ fPhiBoundariesOfSM[1] = TMath::PiOver2() + TMath::ATan2(fParSM[1] , fIPDistance);
+ if(fNumberOfSuperModules > 2) {
+ for(int i=1; i<=4; i++) { // from 2th ro 9th
+ fPhiBoundariesOfSM[2*i] = fPhiBoundariesOfSM[0] + 20.*TMath::DegToRad()*i;
+ fPhiBoundariesOfSM[2*i+1] = fPhiBoundariesOfSM[1] + 20.*TMath::DegToRad()*i;
+ fPhiCentersOfSM[i] = fPhiCentersOfSM[0] + 20.*TMath::DegToRad()*i;
+ }
+ }
+ if(fNumberOfSuperModules > 10) {
+ fPhiBoundariesOfSM[11] = 190.*TMath::DegToRad();
+ fPhiBoundariesOfSM[10] = fPhiBoundariesOfSM[11] - TMath::ATan2((fParSM[1]) , fIPDistance);
+ fPhiCentersOfSM[5] = (fPhiBoundariesOfSM[10]+fPhiBoundariesOfSM[11])/2.;
}
- fPhiBoundariesOfSM[11] = 190.*TMath::DegToRad();
- fPhiBoundariesOfSM[10] = fPhiBoundariesOfSM[11] - TMath::ATan2((fParSM[1]) , fIPDistance);
- fPhiCentersOfSM[5] = (fPhiBoundariesOfSM[10]+fPhiBoundariesOfSM[11])/2.;
-
- //TRU parameters. These parameters values are not the final ones.
- fNTRUEta = 3 ;
- fNTRUPhi = 1 ;
- fNCellsInTRUEta = 16 ;
- fNCellsInTRUPhi = 24 ;
- if(fGeoName.Contains("WSUC")) fNumberOfSuperModules = 1; // Jul 12, 2007
+ //called after setting of scintillator and lead layer parameters
+ DefineSamplingFraction();
+
+ // TRU parameters - Apr 29,08 by PAI.
+ // These parameters values was updated at Nov 05, 2007
+ // As is on Olivier BOURRION (LPSC) ppt preasentation
+ // at ALICE trigger meeting at 13th-14th March
+ fNTRUEta = 1; // was 3
+ fNTRUPhi = 3; // was 1
+ fNModulesInTRUEta = 24; // was 8
+ fNModulesInTRUPhi = 4; // was 12
+ // Jet trigger
+ // 3*6*10 + 2*6*2 = 204 -> matrix (nphi(17), neta(12))
+ fNEtaSubOfTRU = 6;
fgInit = kTRUE;
}
printf( " ECAL : %d x (%f cm Pb, %f cm Sc) \n",
GetNECLayers(), GetECPbRadThick(), GetECScintThick() ) ;
printf(" fSampling %5.2f \n", fSampling );
- if(fGeoName.Contains("SHISH")){
- printf(" fIPDistance %6.3f cm \n", fIPDistance);
- if(fSteelFrontThick>0.)
- printf(" fSteelFrontThick %6.3f cm \n", fSteelFrontThick);
- printf(" fNPhi %i | fNZ %i \n", fNPhi, fNZ);
- printf(" fNCellsInModule %i : fNCellsInSupMod %i : fNCells %i\n",fNCellsInModule, fNCellsInSupMod, fNCells);
- if(fGeoName.Contains("MAY05")){
- printf(" fFrontSteelStrip %6.4f cm (thickness of front steel strip)\n",
- fFrontSteelStrip);
- printf(" fLateralSteelStrip %6.4f cm (thickness of lateral steel strip)\n",
- fLateralSteelStrip);
- printf(" fPassiveScintThick %6.4f cm (thickness of front passive Sc tile)\n",
- fPassiveScintThick);
+ printf(" fIPDistance %6.3f cm \n", fIPDistance);
+ printf(" fNPhi %i | fNZ %i \n", fNPhi, fNZ);
+ printf(" fNCellsInModule %i : fNCellsInSupMod %i : fNCells %i\n",fNCellsInModule, fNCellsInSupMod, fNCells);
+ printf(" X:Y module size %6.3f , %6.3f cm \n", fPhiModuleSize, fEtaModuleSize);
+ printf(" X:Y tile size %6.3f , %6.3f cm \n", fPhiTileSize, fEtaTileSize);
+ printf(" #of sampling layers %i(fNECLayers) \n", fNECLayers);
+ printf(" fLongModuleSize %6.3f cm \n", fLongModuleSize);
+ printf(" #supermodule in phi direction %i \n", fNPhiSuperModule );
+ printf(" fILOSS %i : fIHADR %i \n", fILOSS, fIHADR);
+ printf(" fTrd1Angle %7.4f\n", fTrd1Angle);
+ printf(" f2Trd1Dx2 %7.4f\n", f2Trd1Dx2);
+ printf("SM dimensions(TRD1) : dx %7.2f dy %7.2f dz %7.2f (SMOD, BOX)\n",
+ fParSM[0],fParSM[1],fParSM[2]);
+ printf(" fPhiGapForSM %7.4f cm (%7.4f <- phi size in degree)\n",
+ fPhiGapForSM, TMath::ATan2(fPhiGapForSM,fIPDistance)*TMath::RadToDeg());
+ if(GetKey110DEG()) printf(" Last two modules have size 10 degree in phi (180<phi<190)\n");
+ printf(" phi SM boundaries \n");
+ for(int i=0; i<fPhiBoundariesOfSM.GetSize()/2.; i++) {
+ printf(" %i : %7.5f(%7.2f) -> %7.5f(%7.2f) : center %7.5f(%7.2f) \n", i,
+ fPhiBoundariesOfSM[2*i], fPhiBoundariesOfSM[2*i]*TMath::RadToDeg(),
+ fPhiBoundariesOfSM[2*i+1], fPhiBoundariesOfSM[2*i+1]*TMath::RadToDeg(),
+ fPhiCentersOfSM[i], fPhiCentersOfSM[i]*TMath::RadToDeg());
+ }
+ printf(" fShishKebabTrd1Modules has %i modules : max eta %5.4f \n",
+ fShishKebabTrd1Modules->GetSize(),fEtaMaxOfTRD1);
+
+ printf("\n Cells grid in eta directions : size %i\n", fCentersOfCellsEtaDir.GetSize());
+ for(Int_t i=0; i<fCentersOfCellsEtaDir.GetSize(); i++) {
+ printf(" ind %2.2i : z %8.3f : x %8.3f \n", i,
+ fCentersOfCellsEtaDir.At(i),fCentersOfCellsXDir.At(i));
+ int ind=0; // Nov 21,2006
+ for(Int_t iphi=0; iphi<fCentersOfCellsPhiDir.GetSize(); iphi++) {
+ ind = iphi*fCentersOfCellsEtaDir.GetSize() + i;
+ printf("%6.4f ", fEtaCentersOfCells[ind]);
+ if((iphi+1)%12 == 0) printf("\n");
}
- printf(" X:Y module size %6.3f , %6.3f cm \n", fPhiModuleSize, fEtaModuleSize);
- printf(" X:Y tile size %6.3f , %6.3f cm \n", fPhiTileSize, fEtaTileSize);
- printf(" #of sampling layers %i(fNECLayers) \n", fNECLayers);
- printf(" fLongModuleSize %6.3f cm \n", fLongModuleSize);
- printf(" #supermodule in phi direction %i \n", fNPhiSuperModule );
+ printf("\n");
+
}
- printf(" fILOSS %i : fIHADR %i \n", fILOSS, fIHADR);
- if(fGeoName.Contains("TRD")) {
- printf(" fTrd1Angle %7.4f\n", fTrd1Angle);
- printf(" f2Trd1Dx2 %7.4f\n", f2Trd1Dx2);
- if(fGeoName.Contains("TRD2")) {
- printf(" fTrd2AngleY %7.4f\n", fTrd2AngleY);
- printf(" f2Trd2Dy2 %7.4f\n", f2Trd2Dy2);
- printf(" fTubsR %7.2f cm\n", fTubsR);
- printf(" fTubsTurnAngle %7.4f\n", fTubsTurnAngle);
- printf(" fEmptySpace %7.4f cm\n", fEmptySpace);
- } else if(fGeoName.Contains("TRD1")){
- printf("SM dimensions(TRD1) : dx %7.2f dy %7.2f dz %7.2f (SMOD, BOX)\n",
- fParSM[0],fParSM[1],fParSM[2]);
- printf(" fPhiGapForSM %7.4f cm (%7.4f <- phi size in degree)\n",
- fPhiGapForSM, TMath::ATan2(fPhiGapForSM,fIPDistance)*TMath::RadToDeg());
- if(GetKey110DEG()) printf(" Last two modules have size 10 degree in phi (180<phi<190)\n");
- printf(" phi SM boundaries \n");
- for(int i=0; i<fPhiBoundariesOfSM.GetSize()/2.; i++) {
- printf(" %i : %7.5f(%7.2f) -> %7.5f(%7.2f) : center %7.5f(%7.2f) \n", i,
- fPhiBoundariesOfSM[2*i], fPhiBoundariesOfSM[2*i]*TMath::RadToDeg(),
- fPhiBoundariesOfSM[2*i+1], fPhiBoundariesOfSM[2*i+1]*TMath::RadToDeg(),
- fPhiCentersOfSM[i], fPhiCentersOfSM[i]*TMath::RadToDeg());
- }
- printf(" fShishKebabTrd1Modules has %i modules : max eta %5.4f \n",
- fShishKebabTrd1Modules->GetSize(),fEtaMaxOfTRD1);
-
- printf("\n Cells grid in eta directions : size %i\n", fCentersOfCellsEtaDir.GetSize());
- for(Int_t i=0; i<fCentersOfCellsEtaDir.GetSize(); i++) {
- printf(" ind %2.2i : z %8.3f : x %8.3f \n", i,
- fCentersOfCellsEtaDir.At(i),fCentersOfCellsXDir.At(i));
- int ind=0; // Nov 21,2006
- for(Int_t iphi=0; iphi<fCentersOfCellsPhiDir.GetSize(); iphi++) {
- ind = iphi*fCentersOfCellsEtaDir.GetSize() + i;
- printf("%6.4f ", fEtaCentersOfCells[ind]);
- if((iphi+1)%12 == 0) printf("\n");
- }
- printf("\n");
-
- }
- printf("\n Cells grid in phi directions : size %i\n", fCentersOfCellsPhiDir.GetSize());
- for(Int_t i=0; i<fCentersOfCellsPhiDir.GetSize(); i++) {
- double phi=fPhiCentersOfCells.At(i);
- printf(" ind %2.2i : y %8.3f : phi %7.5f(%6.2f) \n", i, fCentersOfCellsPhiDir.At(i),
- phi, phi*TMath::RadToDeg());
- }
- }
+ printf("\n Cells grid in phi directions : size %i\n", fCentersOfCellsPhiDir.GetSize());
+ for(Int_t i=0; i<fCentersOfCellsPhiDir.GetSize(); i++) {
+ double phi=fPhiCentersOfCells.At(i);
+ printf(" ind %2.2i : y %8.3f : phi %7.5f(%6.2f) \n", i, fCentersOfCellsPhiDir.At(i),
+ phi, phi*TMath::RadToDeg());
}
+
}
+//______________________________________________________________________
void AliEMCALGeometry::PrintCellIndexes(Int_t absId, int pri, char *tit)
{
// Service methods
} else if(fNECLayers == 61) { // 20% layer reduction
fSampling = 12.80;
} else if(fNECLayers == 77) {
- if (fECScintThick>0.175 && fECScintThick<0.177) { // 10% Pb thicknes reduction
+ if (fECScintThick>0.159 && fECScintThick<0.161) { // original sampling fraction, equal layers
+ fSampling = 12.327; // fECScintThick = fECPbRadThickness = 0.160;
+ } else if (fECScintThick>0.175 && fECScintThick<0.177) { // 10% Pb thicknes reduction
fSampling = 10.5; // fECScintThick = 0.176, fECPbRadThickness=0.144;
} else if(fECScintThick>0.191 && fECScintThick<0.193) { // 20% Pb thicknes reduction
fSampling = 8.93; // fECScintThick = 0.192, fECPbRadThickness=0.128;
}
+
}
}
//______________________________________________________________________
-void AliEMCALGeometry::GetCellPhiEtaIndexInSModuleFromTRUIndex(const Int_t itru, const Int_t iphitru, const Int_t ietatru, Int_t &iphiSM, Int_t &ietaSM) const
+void AliEMCALGeometry::GetModulePhiEtaIndexInSModuleFromTRUIndex(const Int_t itru, const Int_t iphitru, const Int_t ietatru, Int_t &iphiSM, Int_t &ietaSM) const
{
- // This method transforms the (eta,phi) index of cells in a
+ // This method transforms the (eta,phi) index of module in a
// TRU matrix into Super Module (eta,phi) index.
// Calculate in which row and column where the TRU are
// ordered in the SM
- Int_t col = itru/ fNTRUPhi ;
+ Int_t col = itru/ fNTRUPhi ; // indexes of TRU in SM
Int_t row = itru - col*fNTRUPhi ;
- iphiSM = fNCellsInTRUPhi*row + iphitru ;
- ietaSM = fNCellsInTRUEta*col + ietatru ;
+ iphiSM = fNModulesInTRUPhi*row + iphitru ;
+ ietaSM = fNModulesInTRUEta*col + ietatru ;
+ //printf(" GetModulePhiEtaIndexInSModuleFromTRUIndex : itru %2i iphitru %2i ietatru %2i iphiSM %2i ietaSM %2i \n",
+ // itru, iphitru, ietatru, iphiSM, ietaSM);
}
//______________________________________________________________________
return rv;
}
+//_____________________________________________________________________________
Bool_t AliEMCALGeometry::IsInEMCAL(Double_t x, Double_t y, Double_t z) const {
// Checks whether point is inside the EMCal volume, used in AliEMCALv*.cxx
//
}
return 0;
}
-// ==
//
// == Shish-kebab cases ==
//
+//________________________________________________________________________________________________
Int_t AliEMCALGeometry::GetAbsCellId(Int_t nSupMod, Int_t nModule, Int_t nIphi, Int_t nIeta) const
{
// 27-aug-04;
return id;
}
+//________________________________________________________________________________________________
Bool_t AliEMCALGeometry::CheckAbsCellId(Int_t absId) const
{
// May 31, 2006; only trd1 now
else return kTRUE;
}
+//________________________________________________________________________________________________
Bool_t AliEMCALGeometry::GetCellIndex(Int_t absId,Int_t &nSupMod,Int_t &nModule,Int_t &nIphi,Int_t &nIeta) const
{
// 21-sep-04; 19-oct-05;
return kTRUE;
}
+//________________________________________________________________________________________________
void AliEMCALGeometry::GetModulePhiEtaIndexInSModule(Int_t nSupMod, Int_t nModule, int &iphim, int &ietam) const
{
// added nSupMod; - 19-oct-05 !
iphim = nModule%nphi;
}
+//________________________________________________________________________________________________
void AliEMCALGeometry::GetCellPhiEtaIndexInSModule(Int_t nSupMod, Int_t nModule, Int_t nIphi, Int_t nIeta,
int &iphi, int &ieta) const
{
nSupMod, nModule, nIphi, nIeta, ieta, iphi));
}
+//________________________________________________________________________________________________
Int_t AliEMCALGeometry::GetSuperModuleNumber(Int_t absId) const
{
// Return the number of the supermodule given the absolute
return nSupMod;
}
+//________________________________________________________________________________________________
void AliEMCALGeometry::GetModuleIndexesFromCellIndexesInSModule(Int_t nSupMod, Int_t iphi, Int_t ieta,
Int_t &iphim, Int_t &ietam, Int_t &nModule) const
{
nModule = ietam * nphi + iphim;
}
+//________________________________________________________________________________________________
Int_t AliEMCALGeometry::GetAbsCellIdFromCellIndexes(Int_t nSupMod, Int_t iphi, Int_t ieta) const
{
// Transition from super module number(nSupMod) and cell indexes (ieta,iphi) to absId
// Methods for AliEMCALRecPoint - Feb 19, 2006
+//________________________________________________________________________________________________
Bool_t AliEMCALGeometry::RelPosCellInSModule(Int_t absId, Double_t &xr, Double_t &yr, Double_t &zr) const
{
// Look to see what the relative
return kTRUE;
}
+//________________________________________________________________________________________________
Bool_t AliEMCALGeometry::RelPosCellInSModule(Int_t absId, Double_t loc[3]) const
{
// Alice numbering scheme - Jun 03, 2006
return kFALSE;
}
+//________________________________________________________________________________________________
Bool_t AliEMCALGeometry::RelPosCellInSModule(Int_t absId, TVector3 &vloc) const
{
static Double_t loc[3];
// Alice numbering scheme - Jun 03, 2006
}
+//________________________________________________________________________________________________
Bool_t AliEMCALGeometry::RelPosCellInSModule(Int_t absId, Double_t distEff, Double_t &xr, Double_t &yr, Double_t &zr) const
{
// Jul 30, 2007 - taking into account position of shower max
return kTRUE;
}
+//________________________________________________________________________________________________
Bool_t AliEMCALGeometry::RelPosCellInSModule(Int_t absId, Int_t maxAbsId, Double_t distEff, Double_t &xr, Double_t &yr, Double_t &zr) const
{
// Jul 31, 2007 - taking into account position of shower max and apply coor2.
return kTRUE;
}
+//________________________________________________________________________________________________
void AliEMCALGeometry::CreateListOfTrd1Modules()
{
// Generate the list of Trd1 modules
}
+//________________________________________________________________________________________________
void AliEMCALGeometry::GetGlobal(const Double_t *loc, Double_t *glob, int ind) const
{
// Figure out the global numbering
}
}
+//________________________________________________________________________________________________
void AliEMCALGeometry::GetGlobal(const TVector3 &vloc, TVector3 &vglob, int ind) const
{
//Figure out the global numbering
vglob.SetXYZ(tglob[0], tglob[1], tglob[2]);
}
+//________________________________________________________________________________________________
void AliEMCALGeometry::GetGlobal(Int_t absId , double glob[3]) const
{
// Alice numbering scheme - Jun 03, 2006
GetGlobal(vloc, vglob, nSupMod);
}
+//________________________________________________________________________________________________
void AliEMCALGeometry::EtaPhiFromIndex(Int_t absId,Double_t &eta,Double_t &phi) const
{
// Nov 16, 2006- float to double
phi = vglob.Phi();
}
+//________________________________________________________________________________________________
void AliEMCALGeometry::EtaPhiFromIndex(Int_t absId,Float_t &eta,Float_t &phi) const
{
// Nov 16,2006 - should be discard in future
phi = float(vglob.Phi());
}
+//________________________________________________________________________________________________
Bool_t AliEMCALGeometry::GetPhiBoundariesOfSM(Int_t nSupMod, Double_t &phiMin, Double_t &phiMax) const
{
// 0<= nSupMod <=11; phi in rad
return kTRUE;
}
+//________________________________________________________________________________________________
Bool_t AliEMCALGeometry::GetPhiBoundariesOfSMGap(Int_t nPhiSec, Double_t &phiMin, Double_t &phiMax) const
{
// 0<= nPhiSec <=4; phi in rad
return kTRUE;
}
+//________________________________________________________________________________________________
Bool_t AliEMCALGeometry::SuperModuleNumberFromEtaPhi(Double_t eta, Double_t phi, Int_t &nSupMod) const
{
// Return false if phi belongs a phi cracks between SM
return kFALSE;
}
+//________________________________________________________________________________________________
Bool_t AliEMCALGeometry::GetAbsCellIdFromEtaPhi(Double_t eta, Double_t phi, Int_t &absId) const
{
// Nov 17,2006
return kFALSE;
}
+//________________________________________________________________________________________________
AliEMCALShishKebabTrd1Module* AliEMCALGeometry::GetShishKebabModule(Int_t neta) const
{
//This method was too long to be
return trd1;
}
-void AliEMCALGeometry::Browse(TBrowser* b) const
+//________________________________________________________________________________________________
+Int_t AliEMCALGeometry::GetAbsTRUNumberFromNumberInSm(const Int_t row, const Int_t col, const Int_t sm)
+{ // Nov 6, 2007
+ Int_t itru = row + col*GetNModulesInTRUPhi() + sm*GetNTRU();
+ // printf(" GetAbsTRUNumberFromNumberInSm : row %2i col %2i sm %2i -> itru %2i\n", row, col, sm, itru);
+ return itru;
+}
+
+//________________________________________________________________________________________________
+void AliEMCALGeometry::Browse(TBrowser* b)
{
//Browse the modules
if(fShishKebabTrd1Modules) b->Add(fShishKebabTrd1Modules);
}
+//________________________________________________________________________________________________
Bool_t AliEMCALGeometry::IsFolder() const
{
//Check if fShishKebabTrd1Modules is in folder
else return kFALSE;
}
+//________________________________________________________________________________________________
Double_t AliEMCALGeometry::GetPhiCenterOfSM(Int_t nsupmod) const
{
//returns center of supermodule in phi