//*-- Author: Sahal Yacoob (LBL / UCT)
// and : Yves Schutz (SUBATECH)
// and : Jennifer Klay (LBL)
-// SHASHLYK : Aleksei Pavlinov (WSU)
+// SHASHLYK : Aleksei Pavlinov (WSU)
//
#include <assert.h>
-// --- AliRoot header files ---
+// --- Root header files ---
#include <Riostream.h>
#include <TBrowser.h>
#include <TClonesArray.h>
#include <TGeoManager.h>
#include <TGeoMatrix.h>
#include <TGeoNode.h>
+#include <TList.h>
#include <TMatrixD.h>
#include <TObjArray.h>
#include <TObjString.h>
+#include <TVector2.h>
#include <TVector3.h>
// -- ALICE Headers.
// 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";
+//
+// Usage:
+// You can create the AliEMCALGeometry object independently from anything.
+// You have to use just the correct name of geometry. If name is empty string the
+// default name of geometry will be used.
+//
+// AliEMCALGeometry* g = AliEMCALGeometry::GetInstance(name,title); // first time
+// ..
+// g = AliEMCALGeometry::GetInstance(); // after first time
+//
+// 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();
AliEMCALGeometry::AliEMCALGeometry()
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),fNTRU(0),fNTRUEta(0),fNTRUPhi(0),fTrd1Angle(0.),f2Trd1Dx2(0.),
+ fNCells(0),fNCellsInSupMod(0),fNCellsInModule(0),fNTRUEta(0),fNTRUPhi(0),
+ fNCellsInTRUEta(0), fNCellsInTRUPhi(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)
+ fShishKebabTrd1Modules(0), fNAdditionalOpts(0),
+ fILOSS(-1), fIHADR(-1)
{
// default ctor only for internal usage (singleton)
// must be kept public for root persistency purposes, but should never be called by the outside world
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),fNTRU(0),fNTRUEta(0),fNTRUPhi(0),fTrd1Angle(0.),f2Trd1Dx2(0.),
+ fNCells(0),fNCellsInSupMod(0),fNCellsInModule(0),fNTRUEta(0),fNTRUPhi(0),
+ fNCellsInTRUEta(0), fNCellsInTRUPhi(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)
+ fShishKebabTrd1Modules(0),fNAdditionalOpts(0),
+ fILOSS(-1), fIHADR(-1)
{
// ctor only for internal usage (singleton)
AliDebug(2, Form("AliEMCALGeometry(%s,%s) ", name,title));
fNCells(geom.fNCells),
fNCellsInSupMod(geom.fNCellsInSupMod),
fNCellsInModule(geom.fNCellsInModule),
- fNTRU(geom.fNTRU),
fNTRUEta(geom.fNTRUEta),
fNTRUPhi(geom.fNTRUPhi),
+ fNCellsInTRUEta(geom.fNCellsInTRUEta),
+ fNCellsInTRUPhi(geom.fNCellsInTRUPhi),
fTrd1Angle(geom.fTrd1Angle),
f2Trd1Dx2(geom.f2Trd1Dx2),
fPhiGapForSM(geom.fPhiGapForSM),
fEtaCentersOfCells(geom.fEtaCentersOfCells),
fPhiCentersOfCells(geom.fPhiCentersOfCells),
fShishKebabTrd1Modules(geom.fShishKebabTrd1Modules),
- fNAdditionalOpts(geom.fNAdditionalOpts)
+ fNAdditionalOpts(geom.fNAdditionalOpts),
+ fILOSS(geom.fILOSS), fIHADR(geom.fIHADR)
{
//copy ctor
}
// Oct 30,2006 - SHISH_TRD1_CURRENT_1X1, SHISH_TRD1_CURRENT_2X2 or SHISH_TRD1_CURRENT_3X3;
//
- fAdditionalOpts[0] = "nl="; // number of sampling layers (fNECLayers)
- fAdditionalOpts[1] = "pbTh="; // cm, Thickness of the Pb (fECPbRadThick)
- fAdditionalOpts[2] = "scTh="; // cm, Thickness of the Sc (fECScintThick)
- fAdditionalOpts[3] = "latSS="; // cm, Thickness of lateral steel strip (fLateralSteelStrip)
+ fAdditionalOpts[0] = "nl="; // number of sampling layers (fNECLayers)
+ fAdditionalOpts[1] = "pbTh="; // cm, Thickness of the Pb (fECPbRadThick)
+ fAdditionalOpts[2] = "scTh="; // cm, Thickness of the Sc (fECScintThick)
+ fAdditionalOpts[3] = "latSS="; // cm, Thickness of lateral steel strip (fLateralSteelStrip)
+ fAdditionalOpts[4] = "allILOSS="; // = 0,1,2,3,4 (4 - energy loss without fluctuation)
+ fAdditionalOpts[5] = "allIHADR="; // = 0,1,2 (0 - no hadronic interaction)
fNAdditionalOpts = sizeof(fAdditionalOpts) / sizeof(char*);
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
- for(int i=0; i<12; i++) fMatrixOfSM[i] = 0;
// geometry
if(fGeoName.Contains("SHISH")){ // Only shahslyk now
fECScintThick = fECPbRadThickness = 0.5;
}
if(fGeoName.Contains("WSUC")){ // 18-may-05 - about common structure
- fShellThickness = 30.; // should be change
+ fShellThickness = 30.; // should be change
fNPhi = fNZ = 4;
}
fPhiBoundariesOfSM[10] = fPhiBoundariesOfSM[11] - TMath::ATan2((fParSM[1]) , fIPDistance);
fPhiCentersOfSM[5] = (fPhiBoundariesOfSM[10]+fPhiBoundariesOfSM[11])/2.;
- fgInit = kTRUE;
-
//TRU parameters. These parameters values are not the final ones.
- fNTRU = 3 ;
fNTRUEta = 3 ;
fNTRUPhi = 1 ;
+ fNCellsInTRUEta = 16 ;
+ fNCellsInTRUPhi = 24 ;
+ if(fGeoName.Contains("WSUC")) fNumberOfSuperModules = 1; // Jul 12, 2007
+
+ fgInit = kTRUE;
}
void AliEMCALGeometry::PrintGeometry()
{
// Separate routine is callable from broswer; Nov 7,2006
- printf("\nInit: geometry of EMCAL named %s is as follows:\n", fGeoName.Data());
+ printf("\nInit: geometry of EMCAL named %s :\n", fGeoName.Data());
+ if(fArrayOpts) {
+ for(Int_t i=0; i<fArrayOpts->GetEntries(); i++){
+ TObjString *o = (TObjString*)fArrayOpts->At(i);
+ printf(" %i : %s \n", i, o->String().Data());
+ }
+ }
printf("Granularity: %d in eta and %d in phi\n", GetNZ(), GetNPhi()) ;
printf("Layout: phi = (%7.1f, %7.1f), eta = (%5.2f, %5.2f), IP = %7.2f -> for EMCAL envelope only\n",
GetArm1PhiMin(), GetArm1PhiMax(),GetArm1EtaMin(), GetArm1EtaMax(), GetIPDistance() );
printf(" fLongModuleSize %6.3f cm \n", fLongModuleSize);
printf(" #supermodule in phi direction %i \n", fNPhiSuperModule );
}
+ 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((iphi+1)%12 == 0) printf("\n");
}
printf("\n");
+
}
printf("\n Cells grid in phi directions : size %i\n", fCentersOfCellsPhiDir.GetSize());
void AliEMCALGeometry::CheckAdditionalOptions()
{
// Feb 06,2006
- //Additional options that
- //can be used to select
- //the specific geometry of
- //EMCAL to run
-
+ // Additional options that
+ // can be used to select
+ // the specific geometry of
+ // EMCAL to run
+ // Dec 27,2006
+ // adeed allILOSS= and allIHADR= for MIP investigation
fArrayOpts = new TObjArray;
Int_t nopt = AliEMCALHistoUtilities::ParseString(fGeoName, *fArrayOpts);
if(nopt==1) { // no aditional option(s)
} else if(addOpt.Contains("LATSS=",TString::kIgnoreCase)) {// Thickness of lateral steel strip (fLateralSteelStrip)
sscanf(addOpt.Data(),"LATSS=%f", &fLateralSteelStrip);
AliDebug(2,Form(" fLateralSteelStrip %f (new) \n", fLateralSteelStrip));
+ } else if(addOpt.Contains("ILOSS=",TString::kIgnoreCase)) {// As in Geant
+ sscanf(addOpt.Data(),"ALLILOSS=%i", &fILOSS);
+ AliDebug(2,Form(" fILOSS %i \n", fILOSS));
+ } else if(addOpt.Contains("IHADR=",TString::kIgnoreCase)) {// As in Geant
+ sscanf(addOpt.Data(),"ALLIHADR=%i", &fIHADR);
+ AliDebug(2,Form(" fIHADR %i \n", fIHADR));
}
}
}
}
}
-//____________________________________________________________________________
-void AliEMCALGeometry::FillTRU(const TClonesArray * digits, TClonesArray * ampmatrix, TClonesArray * timeRmatrix) {
-
-
-// Orders digits ampitudes list in fNTRU TRUs (384 cells) per supermodule.
-// Each TRU is a TMatrixD, and they are kept in TClonesArrays. The number of
-// TRU in phi is fNTRUPhi, and the number of TRU in eta is fNTRUEta.
-// Last 2 modules are half size in Phi, I considered that the number of TRU
-// is maintained for the last modules but decision not taken. If different,
-// then this must be changed.
-
-
- //Check data members
-
- if(fNTRUEta*fNTRUPhi != fNTRU)
- Error("FillTRU"," Wrong number of TRUS per Eta or Phi");
-
- //Initilize and declare variables
- //List of TRU matrices initialized to 0.
- Int_t nCellsPhi = fNPhi*2/fNTRUPhi;
- Int_t nCellsPhi2 = fNPhi/fNTRUPhi; //HalfSize modules
- Int_t nCellsEta = fNZ*2/fNTRUEta;
- Int_t id = -1;
- Float_t amp = -1;
- Float_t timeR = -1;
- Int_t iSupMod = -1;
- Int_t nModule = -1;
- Int_t nIphi = -1;
- Int_t nIeta = -1;
- Int_t iphi = -1;
- Int_t ieta = -1;
-
- //List of TRU matrices initialized to 0.
- for(Int_t k = 0; k < fNTRU*fNumberOfSuperModules; k++){
- TMatrixD * amptrus = new TMatrixD(nCellsPhi,nCellsEta) ;
- TMatrixD * timeRtrus = new TMatrixD(nCellsPhi,nCellsEta) ;
- for(Int_t i = 0; i < nCellsPhi; i++){
- for(Int_t j = 0; j < nCellsEta; j++){
- (*amptrus)(i,j) = 0.0;
- (*timeRtrus)(i,j) = 0.0;
- }
- }
- new((*ampmatrix)[k]) TMatrixD(*amptrus) ;
- new((*timeRmatrix)[k]) TMatrixD(*timeRtrus) ;
- }
-
- AliEMCALDigit * dig ;
-
- //Digits loop to fill TRU matrices with amplitudes.
- for(Int_t idig = 0 ; idig < digits->GetEntriesFast() ; idig++){
-
- dig = dynamic_cast<AliEMCALDigit *>(digits->At(idig)) ;
- amp = dig->GetAmp() ; // Energy of the digit (arbitrary units)
- id = dig->GetId() ; // Id label of the cell
- timeR = dig->GetTimeR() ; // Earliest time of the digit
-
- //Get eta and phi cell position in supermodule
- Bool_t bCell = GetCellIndex(id, iSupMod, nModule, nIphi, nIeta) ;
- if(!bCell)
- Error("FillTRU","Wrong cell id number") ;
-
- GetCellPhiEtaIndexInSModule(iSupMod,nModule,nIphi, nIeta,iphi,ieta);
-
- //Check to which TRU in the supermodule belongs the cell.
- //Supermodules are divided in a TRU matrix of dimension
- //(fNTRUPhi,fNTRUEta).
- //Each TRU is a cell matrix of dimension (nCellsPhi,nCellsEta)
-
- //First calculate the row and column in the supermodule
- //of the TRU to which the cell belongs.
- Int_t col = ieta/nCellsEta;
- Int_t row = iphi/nCellsPhi;
- if(iSupMod > 9)
- row = iphi/nCellsPhi2;
- //Calculate label number of the TRU
- Int_t itru = row + col*fNTRUPhi + iSupMod*fNTRU ;
-
- //Fill TRU matrix with cell values
- TMatrixD * amptrus = dynamic_cast<TMatrixD *>(ampmatrix->At(itru)) ;
- TMatrixD * timeRtrus = dynamic_cast<TMatrixD *>(timeRmatrix->At(itru)) ;
-
- //Calculate row and column of the cell inside the TRU with number itru
- Int_t irow = iphi - row * nCellsPhi;
- if(iSupMod > 9)
- irow = iphi - row * nCellsPhi2;
- Int_t icol = ieta - col * nCellsEta;
-
- (*amptrus)(irow,icol) = amp ;
- (*timeRtrus)(irow,icol) = timeR ;
-
- }
-}
-
//______________________________________________________________________
void AliEMCALGeometry::GetCellPhiEtaIndexInSModuleFromTRUIndex(const Int_t itru, const Int_t iphitru, const Int_t ietatru, Int_t &iphiSM, Int_t &ietaSM) const
{
Int_t col = itru/ fNTRUPhi ;
Int_t row = itru - col*fNTRUPhi ;
- //Calculate the (eta,phi) index in SM
- Int_t nCellsPhi = fNPhi*2/fNTRUPhi;
- Int_t nCellsEta = fNZ*2/fNTRUEta;
-
- iphiSM = nCellsPhi*row + iphitru ;
- ietaSM = nCellsEta*col + ietatru ;
+ iphiSM = fNCellsInTRUPhi*row + iphitru ;
+ ietaSM = fNCellsInTRUEta*col + ietatru ;
}
//______________________________________________________________________
AliEMCALGeometry * rv = 0;
if ( fgGeom == 0 ) {
- if ( strcmp(name,"") == 0 ) rv = 0;
- else {
- fgGeom = new AliEMCALGeometry(name, title);
- if ( fgInit ) rv = (AliEMCALGeometry * ) fgGeom;
- else {
- rv = 0;
- delete fgGeom;
- fgGeom = 0;
- } // end if fgInit
- } // end if strcmp(name,"")
+ if ( strcmp(name,"") == 0 ) { // get default geometry
+ fgGeom = new AliEMCALGeometry(fgDefaultGeometryName, title);
+ } else {
+ fgGeom = new AliEMCALGeometry(name, title);
+ } // end if strcmp(name,"")
+ if ( fgInit ) rv = (AliEMCALGeometry * ) fgGeom;
+ else {
+ rv = 0;
+ delete fgGeom;
+ fgGeom = 0;
+ } // end if fgInit
}else{
if ( strcmp(fgGeom->GetName(), name) != 0) {
printf("\ncurrent geometry is %s : ", fgGeom->GetName());
return 0;
Double_t phi = TMath::ATan2(y,x) * 180./TMath::Pi();
+ if (phi < 0) phi += 360; // phi should go from 0 to 360 in this case
if (phi > fArm1PhiMin && phi < fArm1PhiMax)
return 1;
}
// Shift index taking into account the difference between standard SM
// and SM of half size in phi direction
- const Int_t phiIndexShift = fCentersOfCellsPhiDir.GetSize()/4; // Nov 22, 2006; was 6 for cas 2X2
+ const Int_t kphiIndexShift = fCentersOfCellsPhiDir.GetSize()/4; // Nov 22, 2006; was 6 for cas 2X2
static Int_t nSupMod, nModule, nIphi, nIeta, iphi, ieta;
if(!CheckAbsCellId(absId)) return kFALSE;
if(nSupMod<10) {
yr = fCentersOfCellsPhiDir.At(iphi);
} else {
- yr = fCentersOfCellsPhiDir.At(iphi + phiIndexShift);
+ yr = fCentersOfCellsPhiDir.At(iphi + kphiIndexShift);
}
AliDebug(1,Form("absId %i nSupMod %i iphi %i ieta %i xr %f yr %f zr %f ",absId,nSupMod,iphi,ieta,xr,yr,zr));
// 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
+ // Look to see what the relative
+ // position inside a given cell is
+ // for a recpoint.
+ // In:
+ // absId - cell is as in Geant, 0<= absId < fNCells;
+ // e - cluster energy
+ // OUT:
+ // xr,yr,zr - x,y,z coordinates of cell with absId inside SM
+
+ // Shift index taking into account the difference between standard SM
+ // and SM of half size in phi direction
+ const Int_t kphiIndexShift = fCentersOfCellsPhiDir.GetSize()/4; // Nov 22, 2006; was 6 for cas 2X2
+ static Int_t nSupMod, nModule, nIphi, nIeta, iphi, ieta;
+ static Int_t iphim, ietam;
+ static AliEMCALShishKebabTrd1Module *mod = 0;
+ static TVector2 v;
+ if(!CheckAbsCellId(absId)) return kFALSE;
+
+ GetCellIndex(absId, nSupMod, nModule, nIphi, nIeta);
+ GetModulePhiEtaIndexInSModule(nSupMod, nModule, iphim, ietam);
+ GetCellPhiEtaIndexInSModule(nSupMod,nModule,nIphi,nIeta, iphi, ieta);
+
+ mod = GetShishKebabModule(ietam);
+ mod->GetPositionAtCenterCellLine(nIeta, distEff, v);
+ xr = v.Y() - fParSM[0];
+ zr = v.X() - fParSM[2];
+
+ if(nSupMod<10) {
+ yr = fCentersOfCellsPhiDir.At(iphi);
+ } else {
+ yr = fCentersOfCellsPhiDir.At(iphi + kphiIndexShift);
+ }
+ AliDebug(1,Form("absId %i nSupMod %i iphi %i ieta %i xr %f yr %f zr %f ",absId,nSupMod,iphi,ieta,xr,yr,zr));
+
+ 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.
+ // Look to see what the relative
+ // position inside a given cell is
+ // for a recpoint.
+ // In:
+ // absId - cell is as in Geant, 0<= absId < fNCells;
+ // maxAbsId - abs id of cell with highest energy
+ // e - cluster energy
+ // OUT:
+ // xr,yr,zr - x,y,z coordinates of cell with absId inside SM
+
+ // Shift index taking into account the difference between standard SM
+ // and SM of half size in phi direction
+ const Int_t kphiIndexShift = fCentersOfCellsPhiDir.GetSize()/4; // Nov 22, 2006; was 6 for cas 2X2
+ static Int_t nSupMod, nModule, nIphi, nIeta, iphi, ieta;
+ static Int_t iphim, ietam;
+ static AliEMCALShishKebabTrd1Module *mod = 0;
+ static TVector2 v;
+
+ static Int_t nSupModM, nModuleM, nIphiM, nIetaM, iphiM, ietaM;
+ static Int_t iphimM, ietamM, maxAbsIdCopy=-1;
+ static AliEMCALShishKebabTrd1Module *modM = 0;
+ static Double_t distCorr;
+
+ if(!CheckAbsCellId(absId)) return kFALSE;
+
+ GetCellIndex(absId, nSupMod, nModule, nIphi, nIeta);
+ GetModulePhiEtaIndexInSModule(nSupMod, nModule, iphim, ietam);
+ GetCellPhiEtaIndexInSModule(nSupMod,nModule,nIphi,nIeta, iphi, ieta);
+ mod = GetShishKebabModule(ietam);
+
+ if(absId != maxAbsId) {
+ distCorr = 0.;
+ if(maxAbsIdCopy != maxAbsId) {
+ GetCellIndex(maxAbsId, nSupModM, nModuleM, nIphiM, nIetaM);
+ GetModulePhiEtaIndexInSModule(nSupModM, nModuleM, iphimM, ietamM);
+ GetCellPhiEtaIndexInSModule(nSupModM,nModuleM,nIphiM,nIetaM, iphiM, ietaM);
+ modM = GetShishKebabModule(ietamM); // do I need this ?
+ maxAbsIdCopy = maxAbsId;
+ }
+ if(ietamM !=0) {
+ distCorr = GetEtaModuleSize()*(ietam-ietamM)/TMath::Tan(modM->GetTheta()); // Stay here
+ //printf(" distCorr %f | dist %f | ietam %i -> etamM %i\n", distCorr, dist, ietam, ietamM);
+ }
+ // distEff += distCorr;
+ }
+ // Bad resolution in this case, strong bias vs phi
+ // distEff = 0.0;
+ mod->GetPositionAtCenterCellLine(nIeta, distEff, v); // Stay here
+ xr = v.Y() - fParSM[0];
+ zr = v.X() - fParSM[2];
+
+ if(nSupMod<10) {
+ yr = fCentersOfCellsPhiDir.At(iphi);
+ } else {
+ yr = fCentersOfCellsPhiDir.At(iphi + kphiIndexShift);
+ }
+ AliDebug(1,Form("absId %i nSupMod %i iphi %i ieta %i xr %f yr %f zr %f ",absId,nSupMod,iphi,ieta,xr,yr,zr));
+
+ return kTRUE;
+}
+
void AliEMCALGeometry::CreateListOfTrd1Modules()
{
// Generate the list of Trd1 modules
//
AliDebug(2,Form(" Cells grid in phi directions : size %i\n", fCentersOfCellsPhiDir.GetSize()));
Int_t ind=0; // this is phi index
- Int_t iphi=0, ieta=0, nModule=0, iphiTemp;
+ Int_t ieta=0, nModule=0, iphiTemp;
Double_t xr, zr, theta, phi, eta, r, x,y;
TVector3 vglob;
- Double_t ytCenterModule, ytCenterCell;
+ Double_t ytCenterModule=0.0, ytCenterCell=0.0;
fCentersOfCellsPhiDir.Set(fNPhi*fNPHIdiv);
fPhiCentersOfCells.Set(fNPhi*fNPHIdiv);
- Double_t R0 = GetIPDistance() + GetLongModuleSize()/2.;
+ Double_t r0 = GetIPDistance() + GetLongModuleSize()/2.;
for(Int_t it=0; it<fNPhi; it++) { // cycle on modules
ytCenterModule = -fParSM[1] + fPhiModuleSize*(2*it+1)/2; // center of module
for(Int_t ic=0; ic<fNPHIdiv; ic++) { // cycle on cells in module
// Define grid on phi direction
// Grid is not the same for different eta bin;
// Effect is small but is still here
- phi = TMath::ATan2(ytCenterCell, R0);
+ phi = TMath::ATan2(ytCenterCell, r0);
fPhiCentersOfCells.AddAt(phi, ind);
AliDebug(2,Form(" ind %2.2i : y %8.3f ", ind, fCentersOfCellsPhiDir.At(ind)));
}
-void AliEMCALGeometry::GetTransformationForSM()
-{
- //Uses the geometry manager to
- //load the transformation matrix
- //for the supermodules
-
- static Bool_t transInit=kFALSE;
- if(transInit) return;
-
- int i=0;
- if(gGeoManager == 0) {
- Info("CreateTransformationForSM() "," Load geometry : TGeoManager::Import()");
- assert(0);
- }
- TGeoNode *tn = gGeoManager->GetTopNode();
- TGeoNode *node=0, *xen1 = 0;
- for(i=0; i<tn->GetNdaughters(); i++) {
- node = tn->GetDaughter(i);
- TString ns(node->GetName());
- if(ns.Contains(GetNameOfEMCALEnvelope())) {
- xen1 = node;
- break;
- }
- }
- if(!xen1) {
- Info("CreateTransformationForSM() "," geometry has not EMCAL envelope with name %s",
- GetNameOfEMCALEnvelope());
- assert(0);
- }
- printf(" i %i : EMCAL Envelope is %s : #SM %i \n", i, xen1->GetName(), xen1->GetNdaughters());
- for(i=0; i<xen1->GetNdaughters(); i++) {
- TGeoNodeMatrix *sm = (TGeoNodeMatrix*)xen1->GetDaughter(i);
- fMatrixOfSM[i] = sm->GetMatrix();
- //Compiler doesn't like this syntax...
- // printf(" %i : matrix %x \n", i, fMatrixOfSM[i]);
- }
- transInit = kTRUE;
-}
-
void AliEMCALGeometry::GetGlobal(const Double_t *loc, Double_t *glob, int ind) const
{
// Figure out the global numbering
// of a given supermodule from the
- // local numbering
- // Alice numbering - Jun 03,2006
- // if(fMatrixOfSM[0] == 0) GetTransformationForSM();
+ // local numbering and the transformation
+ // matrix stored by the geometry manager (allows for misaligned
+ // geometry)
if(ind>=0 && ind < GetNumberOfSuperModules()) {
- fMatrixOfSM[ind]->LocalToMaster(loc, glob);
+ TString volpath = "ALIC_1/XEN1_1/SMOD_";
+ volpath += ind+1;
+
+ if(GetKey110DEG() && ind>=10) {
+ volpath = "ALIC_1/XEN1_1/SM10_";
+ volpath += ind-10+1;
+ }
+
+ if(!gGeoManager->cd(volpath.Data()))
+ AliFatal(Form("AliEMCALGeometry::GeoManager cannot find path %s!",volpath.Data()));
+
+ TGeoHMatrix* m = gGeoManager->GetCurrentMatrix();
+ if(m) {
+ m->LocalToMaster(loc, glob);
+ } else {
+ AliFatal("Geo matrixes are not loaded \n") ;
+ }
}
}
glob[0]=glob[1]=glob[2]=0.0; // bad case
if(RelPosCellInSModule(absId, loc)) {
GetCellIndex(absId, nSupMod, nModule, nIphi, nIeta);
- fMatrixOfSM[nSupMod]->LocalToMaster(loc, glob);
+
+ TString volpath = "ALIC_1/XEN1_1/SMOD_";
+ volpath += (nSupMod+1);
+
+ if(GetKey110DEG() && nSupMod>=10) {
+ volpath = "ALIC_1/XEN1_1/SM10_";
+ volpath += (nSupMod-10+1);
+ }
+ if(!gGeoManager->cd(volpath.Data()))
+ AliFatal(Form("GeoManager cannot find path %s!",volpath.Data()));
+
+ TGeoHMatrix* m = gGeoManager->GetCurrentMatrix();
+ if(m) {
+ m->LocalToMaster(loc, glob);
+ } else {
+ AliFatal("Geo matrixes are not loaded \n") ;
+ }
}
}
+//___________________________________________________________________
void AliEMCALGeometry::GetGlobal(Int_t absId , TVector3 &vglob) const
{
// Alice numbering scheme - Jun 03, 2006
}
-void AliEMCALGeometry::GetGlobal(const AliRecPoint *rp, TVector3 &vglob) const
+//____________________________________________________________________________
+void AliEMCALGeometry::GetGlobal(const AliRecPoint* /*rp*/, TVector3& /* vglob */) const
+{
+ AliFatal(Form("Please use GetGlobalEMCAL(recPoint,gpos) instead of GetGlobal!"));
+}
+
+//_________________________________________________________________________________
+void AliEMCALGeometry::GetGlobalEMCAL(const AliEMCALRecPoint *rp, TVector3 &vglob) const
{
// Figure out the global numbering
// of a given supermodule from the
static TVector3 vloc;
static Int_t nSupMod, nModule, nIphi, nIeta;
- AliRecPoint *rpTmp = (AliRecPoint*)rp; // const_cast ??
- if(!rpTmp) return;
- AliEMCALRecPoint *rpEmc = (AliEMCALRecPoint*)rpTmp;
+ const AliEMCALRecPoint *rpTmp = rp;
+ const AliEMCALRecPoint *rpEmc = rpTmp;
GetCellIndex(rpEmc->GetAbsId(0), nSupMod, nModule, nIphi, nIeta);
rpTmp->GetLocalPosition(vloc);
return kFALSE;
}
-AliEMCALShishKebabTrd1Module* AliEMCALGeometry::GetShishKebabModule(Int_t neta)
+AliEMCALShishKebabTrd1Module* AliEMCALGeometry::GetShishKebabModule(Int_t neta) const
{
//This method was too long to be
//included in the header file - the
return trd1;
}
-void AliEMCALGeometry::Browse(TBrowser* b)
+void AliEMCALGeometry::Browse(TBrowser* b) const
{
+ //Browse the modules
if(fShishKebabTrd1Modules) b->Add(fShishKebabTrd1Modules);
}
Bool_t AliEMCALGeometry::IsFolder() const
{
+ //Check if fShishKebabTrd1Modules is in folder
if(fShishKebabTrd1Modules) return kTRUE;
else return kFALSE;
}
+
+Double_t AliEMCALGeometry::GetPhiCenterOfSM(Int_t nsupmod) const
+{
+ //returns center of supermodule in phi
+ static int i = nsupmod/2;
+ return fPhiCentersOfSM[i];
+
+}