#include <assert.h>
-// --- AliRoot header files ---
+// --- Root header files ---
#include <Riostream.h>
#include <TBrowser.h>
#include <TClonesArray.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";
+Char_t* AliEMCALGeometry::fgDefaultGeometryName = "SHISH_77_TRD1_2X2_FINAL_110DEG scTh=0.176 pbTh=0.144";
//
// Usage:
// You can create the AliEMCALGeometry object independently from anything.
// == =============================
// AliRunLoader *rl = AliRunLoader::GetRunLoader();
// AliEMCALGeometry *geom = dynamic_cast<AliEMCAL*>(rl->GetAliRun()->GetDetector("EMCAL"))->GetGeometry();
-
+// TGeoManager::Import("geometry.root");
AliEMCALGeometry::AliEMCALGeometry()
: AliGeometry(),
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
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
// 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()) ;
fNCellsInTRUEta = 16 ;
fNCellsInTRUPhi = 24 ;
- // Define TGeoMatrix of SM - Jan 19, 2007 (just fro TRD1)
- if(fGeoName.Contains("TRD1")) { // copy code from AliEMCALv0::CreateSmod()
- int nphism = GetNumberOfSuperModules()/2;
- double dphi = (GetArm1PhiMax() - GetArm1PhiMin())/nphism;
- double rpos = (GetEnvelop(0) + GetEnvelop(1))/2.;
- double phi, phiRad, xpos, ypos, zpos;
- for(int i=0; i<nphism; i++){
- phi = GetArm1PhiMin() + dphi*(2*i+1)/2.; // phi= 90, 110, 130, 150, 170, 190
- phiRad = phi*TMath::Pi()/180.;
- xpos = rpos * TMath::Cos(phiRad);
- ypos = rpos * TMath::Sin(phiRad);
- zpos = fParSM[2];
- if(i==5) {
- xpos += (fParSM[1]/2. * TMath::Sin(phiRad));
- ypos -= (fParSM[1]/2. * TMath::Cos(phiRad));
- }
- // pozitive z
- int ind = 2*i;
- TGeoRotation *geoRot0 = new TGeoRotation("geoRot0", 90.0, phi, 90.0, 90.0+phi, 0.0, 0.0);
- fMatrixOfSM[ind] = new TGeoCombiTrans(Form("EmcalSM%2.2i",ind),
- xpos,ypos, zpos, geoRot0);
- // negaive z
- ind++;
- double phiy = 90. + phi + 180.;
- if(phiy>=360.) phiy -= 360.;
- TGeoRotation *geoRot1 = new TGeoRotation("geoRot1", 90.0, phi, 90.0, phiy, 180.0, 0.0);
- fMatrixOfSM[ind] = new TGeoCombiTrans(Form("EmcalSM%2.2i",ind),
- xpos,ypos,-zpos, geoRot1);
- } // for
- }
-
if(fGeoName.Contains("WSUC")) fNumberOfSuperModules = 1; // Jul 12, 2007
fgInit = kTRUE;
- AliInfo(" is ended");
}
void AliEMCALGeometry::PrintGeometry()
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")){
+ 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",
printf("\n");
}
- printf(" Matrix transformation\n");
- for(Int_t i=0; i<12; i++) {
- TGeoMatrix *m = fMatrixOfSM[i];
- if(m==0) continue;
- const double *xyz = m->GetTranslation();
- printf(" %2.2i %s %s x %7.2f y %7.2f z %7.2f\n",
- i, m->GetName(), m->ClassName(), xyz[0],xyz[1],xyz[2]);
- }
printf("\n Cells grid in phi directions : size %i\n", fCentersOfCellsPhiDir.GetSize());
for(Int_t i=0; i<fCentersOfCellsPhiDir.GetSize(); i++) {
}
}
}
- cout<<endl;
}
void AliEMCALGeometry::PrintCellIndexes(Int_t absId, int pri, char *tit)
}else{
if ( strcmp(fgGeom->GetName(), name) != 0) {
printf("\ncurrent geometry is %s : ", fgGeom->GetName());
- printf(" you cannot call %s ", name);
+ printf(" you cannot call %s ",name);
}else{
rv = (AliEMCALGeometry *) fgGeom;
} // end
// 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));
// 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;
static Int_t iphim, ietam;
static AliEMCALShishKebabTrd1Module *mod = 0;
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));
// 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;
static Int_t iphim, ietam;
static AliEMCALShishKebabTrd1Module *mod = 0;
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));
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
- // Unused after 19 Jan, 2007 - keep for compatibility;
-
- return;
- 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") ;
+ }
}
}
static Int_t nSupMod, nModule, nIphi, nIeta;
static double loc[3];
+ if (!gGeoManager || !gGeoManager->IsClosed()) {
+ AliError("Can't get the global coordinates! gGeoManager doesn't exist or it is still open!");
+ return;
+ }
+
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);
void AliEMCALGeometry::Browse(TBrowser* b)
{
+ //Browse the modules
if(fShishKebabTrd1Modules) b->Add(fShishKebabTrd1Modules);
- for(int i=0; i<fNumberOfSuperModules; i++) {
- if(fMatrixOfSM[i]) b->Add(fMatrixOfSM[i]);
- }
}
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
+ int i = nsupmod/2;
+ return fPhiCentersOfSM[i];
+
+}