clean up the EMCAL name scheme : super module -> module -> tower (or cell)
authorpavlinov <pavlinov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 19 Dec 2006 02:34:13 +0000 (02:34 +0000)
committerpavlinov <pavlinov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 19 Dec 2006 02:34:13 +0000 (02:34 +0000)
EMCAL/AliEMCAL.cxx
EMCAL/AliEMCALClusterizerv1.cxx
EMCAL/AliEMCALDigitizer.cxx
EMCAL/AliEMCALGeometry.cxx
EMCAL/AliEMCALGeometry.h
EMCAL/AliEMCALGeometryOfflineTrd1.cxx
EMCAL/AliEMCALGeometryOfflineTrd1.h
EMCAL/AliEMCALRecPoint.cxx
EMCAL/AliEMCALv2.cxx

index 1a43c31..7916780 100644 (file)
@@ -17,6 +17,9 @@
 /* History of cvs commits:
  *
  * $Log$
+ * Revision 1.47  2006/12/05 17:12:03  gustavo
+ * Updated AliEMCAL::Digits2Raw, reads first provisional RCU mapping files to make Raw data with new AliCaloAltroMapping and AliCaloRawStream
+ *
  *
  */
 //_________________________________________________________________________
@@ -309,9 +312,9 @@ void AliEMCAL::Digits2Raw()
     Int_t nIeta = 0;
     Int_t iphi = 0;
     Int_t ieta = 0;
-    Int_t nTower = 0;
-    geom->GetCellIndex(digit->GetId(), nSM, nTower, nIphi, nIeta);
-    geom->GetCellPhiEtaIndexInSModule(nSM, nTower, nIphi, nIeta,iphi, ieta) ;
+    Int_t nModule = 0;
+    geom->GetCellIndex(digit->GetId(), nSM, nModule, nIphi, nIeta);
+    geom->GetCellPhiEtaIndexInSModule(nSM, nModule, nIphi, nIeta,iphi, ieta) ;
     
     //Check which is the RCU of the cell.
     Int_t iRCU = -111;
index 7591979..a298ffe 100644 (file)
@@ -186,19 +186,19 @@ Float_t  AliEMCALClusterizerv1::Calibrate(Int_t amp, Int_t AbsId)
       AliFatal("Did not get geometry from EMCALLoader") ;
     
     Int_t iSupMod = -1;
-    Int_t nTower  = -1;
+    Int_t nModule  = -1;
     Int_t nIphi   = -1;
     Int_t nIeta   = -1;
     Int_t iphi    = -1;
     Int_t ieta    = -1;
     
-    Bool_t bCell = fGeom->GetCellIndex(AbsId, iSupMod, nTower, nIphi, nIeta) ;
+    Bool_t bCell = fGeom->GetCellIndex(AbsId, iSupMod, nModule, nIphi, nIeta) ;
     if(!bCell) {
       fGeom->PrintGeometry();
       Error("Calibrate()"," Wrong cell id number : %i", AbsId);
       assert(0);
     }
-    fGeom->GetCellPhiEtaIndexInSModule(iSupMod,nTower,nIphi, nIeta,iphi,ieta);
+    fGeom->GetCellPhiEtaIndexInSModule(iSupMod,nModule,nIphi, nIeta,iphi,ieta);
 
     fADCchannelECA  = fCalibData->GetADCchannel (iSupMod,ieta,iphi);
     fADCpedestalECA = fCalibData->GetADCpedestal(iSupMod,ieta,iphi);
@@ -438,17 +438,17 @@ Int_t AliEMCALClusterizerv1::AreNeighbours(AliEMCALDigit * d1, AliEMCALDigit * d
   //                                      which is compared to a digit (d2)  not yet in a cluster  
 
   static Int_t rv; 
-  static Int_t nSupMod1=0, nTower1=0, nIphi1=0, nIeta1=0, iphi1=0, ieta1=0;
-  static Int_t nSupMod2=0, nTower2=0, nIphi2=0, nIeta2=0, iphi2=0, ieta2=0;
+  static Int_t nSupMod1=0, nModule1=0, nIphi1=0, nIeta1=0, iphi1=0, ieta1=0;
+  static Int_t nSupMod2=0, nModule2=0, nIphi2=0, nIeta2=0, iphi2=0, ieta2=0;
   static Int_t rowdiff, coldiff;
   rv = 0 ; 
 
-  fGeom->GetCellIndex(d1->GetId(), nSupMod1,nTower1,nIphi1,nIeta1);
-  fGeom->GetCellIndex(d2->GetId(), nSupMod2,nTower2,nIphi2,nIeta2);
+  fGeom->GetCellIndex(d1->GetId(), nSupMod1,nModule1,nIphi1,nIeta1);
+  fGeom->GetCellIndex(d2->GetId(), nSupMod2,nModule2,nIphi2,nIeta2);
   if(nSupMod1 != nSupMod2) return 2; // different SM
 
-  fGeom->GetCellPhiEtaIndexInSModule(nSupMod1,nTower1,nIphi1,nIeta1, iphi1,ieta1);
-  fGeom->GetCellPhiEtaIndexInSModule(nSupMod2,nTower2,nIphi2,nIeta2, iphi2,ieta2);
+  fGeom->GetCellPhiEtaIndexInSModule(nSupMod1,nModule1,nIphi1,nIeta1, iphi1,ieta1);
+  fGeom->GetCellPhiEtaIndexInSModule(nSupMod2,nModule2,nIphi2,nIeta2, iphi2,ieta2);
 
   rowdiff = TMath::Abs(iphi1 - iphi2);  
   coldiff = TMath::Abs(ieta1 - ieta2) ;  
@@ -482,21 +482,21 @@ Int_t AliEMCALClusterizerv1::AreInGroup(AliEMCALDigit * d1, AliEMCALDigit * d2)
   //find another in the list that will match?  How about my TowerGroup search?
 
   static Int_t rv;
-  static Int_t nSupMod1=0, nTower1=0, nIphi1=0, nIeta1=0, iphi1=0, ieta1=0;
-  static Int_t nSupMod2=0, nTower2=0, nIphi2=0, nIeta2=0, iphi2=0, ieta2=0;
+  static Int_t nSupMod1=0, nModule1=0, nIphi1=0, nIeta1=0, iphi1=0, ieta1=0;
+  static Int_t nSupMod2=0, nModule2=0, nIphi2=0, nIeta2=0, iphi2=0, ieta2=0;
   static Int_t towerGroup1 = -1, towerGroup2 = -1;
   rv = 0 ;
 
-  fGeom->GetCellIndex(d1->GetId(), nSupMod1,nTower1,nIphi1,nIeta1);
-  fGeom->GetCellIndex(d2->GetId(), nSupMod2,nTower2,nIphi2,nIeta2);
+  fGeom->GetCellIndex(d1->GetId(), nSupMod1,nModule1,nIphi1,nIeta1);
+  fGeom->GetCellIndex(d2->GetId(), nSupMod2,nModule2,nIphi2,nIeta2);
   if(nSupMod1 != nSupMod2) return 2; // different SM
 
-  static Int_t nTowerInSM = fGeom->GetNCellsInSupMod()/fGeom->GetNCellsInTower();
+  static Int_t nTowerInSM = fGeom->GetNCellsInSupMod()/fGeom->GetNCellsInModule();
 
   //figure out which tower grouping each digit belongs to
   for(int it = 0; it < nTowerInSM/fNTowerInGroup; it++) {
-    if(nTower1 <= nTowerInSM - it*fNTowerInGroup) towerGroup1 = it;
-    if(nTower2 <= nTowerInSM - it*fNTowerInGroup) towerGroup2 = it;
+    if(nModule1 <= nTowerInSM - it*fNTowerInGroup) towerGroup1 = it;
+    if(nModule2 <= nTowerInSM - it*fNTowerInGroup) towerGroup2 = it;
   }
   if(towerGroup1 != towerGroup2) return 3; //different Towergroup
 
index 90ba49f..275170f 100644 (file)
@@ -464,17 +464,17 @@ Int_t AliEMCALDigitizer::DigitizeEnergy(Float_t energy, Int_t AbsId)
     AliFatal("Did not get geometry from EMCALLoader");
 
   Int_t iSupMod = -1;
-  Int_t nTower  = -1;
+  Int_t nModule  = -1;
   Int_t nIphi   = -1;
   Int_t nIeta   = -1;
   Int_t iphi    = -1;
   Int_t ieta    = -1;
   Int_t channel = -999; 
 
-  Bool_t bCell = geom->GetCellIndex(AbsId, iSupMod, nTower, nIphi, nIeta) ;
+  Bool_t bCell = geom->GetCellIndex(AbsId, iSupMod, nModule, nIphi, nIeta) ;
   if(!bCell)
     Error("DigitizeEnergy","Wrong cell id number : AbsId %i ", AbsId) ;
-  geom->GetCellPhiEtaIndexInSModule(iSupMod,nTower,nIphi, nIeta,iphi,ieta);
+  geom->GetCellPhiEtaIndexInSModule(iSupMod,nModule,nIphi, nIeta,iphi,ieta);
   
   if(emcalLoader->CalibData()) {
     fADCpedestalEC = emcalLoader->CalibData()
index cfc169c..0585295 100644 (file)
@@ -24,9 +24,9 @@
 // Number of Modules and Layers may be controlled by 
 // the name of the instance defined               
 //     EMCAL geometry tree:
-//     EMCAL -> superModule -> module -> cell
+//     EMCAL -> superModule -> module -> tower(cell)
 //     Indexes
-//     absId -> nSupMod     -> nTower -> (nIphi,nIeta)
+//     absId -> nSupMod     -> nModule -> (nIphi,nIeta)
 //
 //*-- Author: Sahal Yacoob (LBL / UCT)
 //     and  : Yves Schutz (SUBATECH)
@@ -73,7 +73,7 @@ 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),fNCellsInTower(0),fNTRU(0),fNTRUEta(0),fNTRUPhi(0),fTrd1Angle(0.),f2Trd1Dx2(0.),
+    fNCells(0),fNCellsInSupMod(0),fNCellsInModule(0),fNTRU(0),fNTRUEta(0),fNTRUPhi(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),
@@ -92,7 +92,7 @@ AliEMCALGeometry::AliEMCALGeometry(const Text_t* name, const Text_t* title)
     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),fNCellsInTower(0),fNTRU(0),fNTRUEta(0),fNTRUPhi(0),fTrd1Angle(0.),f2Trd1Dx2(0.),
+    fNCells(0),fNCellsInSupMod(0),fNCellsInModule(0),fNTRU(0),fNTRUEta(0),fNTRUPhi(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),
@@ -145,7 +145,7 @@ AliEMCALGeometry::AliEMCALGeometry(const AliEMCALGeometry& geom)
     fNETAdiv(geom.fNETAdiv),
     fNCells(geom.fNCells),
     fNCellsInSupMod(geom.fNCellsInSupMod),
-    fNCellsInTower(geom.fNCellsInTower),
+    fNCellsInModule(geom.fNCellsInModule),
     fNTRU(geom.fNTRU),
     fNTRUEta(geom.fNTRUEta),
     fNTRUPhi(geom.fNTRUPhi),
@@ -322,8 +322,8 @@ void AliEMCALGeometry::Init(void){
     fEtaTileSize = fEtaModuleSize/double(fNETAdiv) - fLateralSteelStrip; // 13-may-05 
 
     // constant for transition absid <--> indexes
-    fNCellsInTower  = fNPHIdiv*fNETAdiv;
-    fNCellsInSupMod = fNCellsInTower*fNPhi*fNZ;
+    fNCellsInModule  = fNPHIdiv*fNETAdiv;
+    fNCellsInSupMod = fNCellsInModule*fNPhi*fNZ;
     fNCells         = fNCellsInSupMod*fNumberOfSuperModules;
     if(GetKey110DEG()) fNCells -= fNCellsInSupMod;
 
@@ -405,7 +405,7 @@ void AliEMCALGeometry::PrintGeometry()
     if(fSteelFrontThick>0.) 
     printf(" fSteelFrontThick  %6.3f cm \n", fSteelFrontThick);
     printf(" fNPhi %i   |  fNZ %i \n", fNPhi, fNZ);
-    printf(" fNCellsInTower %i : fNCellsInSupMod %i : fNCells %i\n",fNCellsInTower, fNCellsInSupMod, fNCells);
+    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);
@@ -471,14 +471,14 @@ void AliEMCALGeometry::PrintGeometry()
 void AliEMCALGeometry::PrintCellIndexes(Int_t absId, int pri, char *tit)
 {
   // Service methods
-  Int_t nSupMod, nTower, nIphi, nIeta;
+  Int_t nSupMod, nModule, nIphi, nIeta;
   Int_t iphi, ieta;
   TVector3 vg;
 
-  GetCellIndex(absId,  nSupMod, nTower, nIphi, nIeta);
-  printf(" %s | absId : %i -> nSupMod %i nTower %i nIphi %i nIeta %i \n", tit, absId,  nSupMod, nTower, nIphi, nIeta);
+  GetCellIndex(absId,  nSupMod, nModule, nIphi, nIeta);
+  printf(" %s | absId : %i -> nSupMod %i nModule %i nIphi %i nIeta %i \n", tit, absId,  nSupMod, nModule, nIphi, nIeta);
   if(pri>0) {
-    GetCellPhiEtaIndexInSModule(nSupMod,nTower,nIphi,nIeta, iphi,ieta);
+    GetCellPhiEtaIndexInSModule(nSupMod,nModule,nIphi,nIeta, iphi,ieta);
     printf(" local SM index : iphi %i : ieta %i \n", iphi,ieta);
     GetGlobal(absId, vg);
     printf(" vglob : mag %7.2f : perp %7.2f : z %7.2f : eta %6.4f : phi %6.4f(%6.2f) \n", 
@@ -582,7 +582,7 @@ void AliEMCALGeometry::FillTRU(const TClonesArray * digits, TClonesArray * ampma
   Float_t amp   = -1;
   Float_t timeR = -1;
   Int_t iSupMod = -1;
-  Int_t nTower  = -1;
+  Int_t nModule  = -1;
   Int_t nIphi   = -1;
   Int_t nIeta   = -1;
   Int_t iphi    = -1;
@@ -613,11 +613,11 @@ void AliEMCALGeometry::FillTRU(const TClonesArray * digits, TClonesArray * ampma
     timeR  = dig->GetTimeR() ; // Earliest time of the digit
    
     //Get eta and phi cell position in supermodule
-    Bool_t bCell = GetCellIndex(id, iSupMod, nTower, nIphi, nIeta) ;
+    Bool_t bCell = GetCellIndex(id, iSupMod, nModule, nIphi, nIeta) ;
     if(!bCell)
       Error("FillTRU","Wrong cell id number") ;
     
-    GetCellPhiEtaIndexInSModule(iSupMod,nTower,nIphi, nIeta,iphi,ieta);
+    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 
@@ -738,14 +738,14 @@ Bool_t AliEMCALGeometry::IsInEMCAL(Double_t x, Double_t y, Double_t z) const {
 //
 // == Shish-kebab cases ==
 //
-Int_t AliEMCALGeometry::GetAbsCellId(Int_t nSupMod, Int_t nTower, Int_t nIphi, Int_t nIeta) const
+Int_t AliEMCALGeometry::GetAbsCellId(Int_t nSupMod, Int_t nModule, Int_t nIphi, Int_t nIeta) const
 { 
   // 27-aug-04; 
   // corr. 21-sep-04; 
   //       13-oct-05; 110 degree case
   // May 31, 2006; ALICE numbering scheme:
   // 0 <= nSupMod < fNumberOfSuperModules
-  // 0 <= nTower  < fNPHI * fNZ ( fNPHI * fNZ/2 for fKey110DEG=1)
+  // 0 <= nModule  < fNPHI * fNZ ( fNPHI * fNZ/2 for fKey110DEG=1)
   // 0 <= nIphi   < fNPHIdiv
   // 0 <= nIeta   < fNETAdiv
   // 0 <= absid   < fNCells
@@ -755,7 +755,7 @@ Int_t AliEMCALGeometry::GetAbsCellId(Int_t nSupMod, Int_t nTower, Int_t nIphi, I
   } else {
     id  = fNCellsInSupMod*nSupMod;
   }
-  id += fNCellsInTower *nTower;
+  id += fNCellsInModule *nModule;
   id += fNPHIdiv *nIphi;
   id += nIeta;
   if(id<0 || id >= fNCells) {
@@ -763,7 +763,7 @@ Int_t AliEMCALGeometry::GetAbsCellId(Int_t nSupMod, Int_t nTower, Int_t nIphi, I
 //     printf("    id      %6i(will be force to -1)\n", id);
 //     printf("    fNCells %6i\n", fNCells);
 //     printf("    nSupMod %6i\n", nSupMod);
-//     printf("    nTower  %6i\n", nTower);
+//     printf("    nModule  %6i\n", nModule);
 //     printf("    nIphi   %6i\n", nIphi);
 //     printf("    nIeta   %6i\n", nIeta);
     id = -TMath::Abs(id); // if negative something wrong
@@ -778,7 +778,7 @@ Bool_t  AliEMCALGeometry::CheckAbsCellId(Int_t absId) const
   else                            return kTRUE;
 }
 
-Bool_t AliEMCALGeometry::GetCellIndex(Int_t absId,Int_t &nSupMod,Int_t &nTower,Int_t &nIphi,Int_t &nIeta) const
+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;
   // May 31, 2006; ALICE numbering scheme:
@@ -787,7 +787,7 @@ Bool_t AliEMCALGeometry::GetCellIndex(Int_t absId,Int_t &nSupMod,Int_t &nTower,I
   // absId   - cell is as in Geant,     0<= absId   < fNCells;
   // Out:
   // nSupMod - super module(SM) number, 0<= nSupMod < fNumberOfSuperModules;
-  // nTower  - module number in SM,     0<= nTower  < fNCellsInSupMod/fNCellsInSupMod or(/2) for tow last SM (10th and 11th);
+  // nModule  - module number in SM,     0<= nModule  < fNCellsInSupMod/fNCellsInSupMod or(/2) for tow last SM (10th and 11th);
   // nIphi   - cell number in phi driection inside module; 0<= nIphi < fNPHIdiv; 
   // nIeta   - cell number in eta driection inside module; 0<= nIeta < fNETAdiv; 
   // 
@@ -803,15 +803,15 @@ Bool_t AliEMCALGeometry::GetCellIndex(Int_t absId,Int_t &nSupMod,Int_t &nTower,I
     tmp     = absId % fNCellsInSupMod;
   }
 
-  nTower  = tmp / fNCellsInTower;
-  tmp     = tmp % fNCellsInTower;
+  nModule  = tmp / fNCellsInModule;
+  tmp     = tmp % fNCellsInModule;
   nIphi   = tmp / fNPHIdiv;
   nIeta   = tmp % fNPHIdiv;
 
   return kTRUE;
 }
 
-void AliEMCALGeometry::GetModulePhiEtaIndexInSModule(Int_t nSupMod, Int_t nTower,  int &iphim, int &ietam) const
+void AliEMCALGeometry::GetModulePhiEtaIndexInSModule(Int_t nSupMod, Int_t nModule,  int &iphim, int &ietam) const
 { 
   // added nSupMod; - 19-oct-05 !
   // Alice numbering scheme        - Jun 01,2006 
@@ -823,11 +823,11 @@ void AliEMCALGeometry::GetModulePhiEtaIndexInSModule(Int_t nSupMod, Int_t nTower
   if(fKey110DEG == 1 && nSupMod>=10) nphi = fNPhi/2;
   else                               nphi = fNPhi;
 
-  ietam = nTower/nphi;
-  iphim = nTower%nphi;
+  ietam = nModule/nphi;
+  iphim = nModule%nphi;
 }
 
-void AliEMCALGeometry::GetCellPhiEtaIndexInSModule(Int_t nSupMod, Int_t nTower, Int_t nIphi, Int_t nIeta, 
+void AliEMCALGeometry::GetCellPhiEtaIndexInSModule(Int_t nSupMod, Int_t nModule, Int_t nIphi, Int_t nIeta, 
 int &iphi, int &ieta) const
 { 
   // 
@@ -835,7 +835,7 @@ int &iphi, int &ieta) const
   // Alice numbering scheme  - Jun 01,2006 
   // IN:
   // nSupMod - super module(SM) number, 0<= nSupMod < fNumberOfSuperModules;
-  // nTower  - module number in SM,     0<= nTower  < fNCellsInSupMod/fNCellsInSupMod or(/2) for tow last SM (10th and 11th);
+  // nModule  - module number in SM,     0<= nModule  < fNCellsInSupMod/fNCellsInSupMod or(/2) for tow last SM (10th and 11th);
   // nIphi   - cell number in phi driection inside module; 0<= nIphi < fNPHIdiv; 
   // nIeta   - cell number in eta driection inside module; 0<= nIeta < fNETAdiv; 
   // 
@@ -846,14 +846,14 @@ int &iphi, int &ieta) const
   //
   static Int_t iphim, ietam;
 
-  GetModulePhiEtaIndexInSModule(nSupMod,nTower, iphim, ietam); 
+  GetModulePhiEtaIndexInSModule(nSupMod,nModule, iphim, ietam); 
   //  ieta  = ietam*fNETAdiv + (1-nIeta); // x(module) = -z(SM) 
   ieta  = ietam*fNETAdiv + (fNETAdiv - 1 - nIeta); // x(module) = -z(SM) 
   iphi  = iphim*fNPHIdiv + nIphi;     // y(module) =  y(SM) 
 
   if(iphi<0 || ieta<0)
-  AliDebug(1,Form(" nSupMod %i nTower %i nIphi %i nIeta %i => ieta %i iphi %i\n", 
-  nSupMod, nTower, nIphi, nIeta, ieta, iphi));
+  AliDebug(1,Form(" nSupMod %i nModule %i nIphi %i nIeta %i => ieta %i iphi %i\n", 
+  nSupMod, nModule, nIphi, nIeta, ieta, iphi));
 }
 
 Int_t  AliEMCALGeometry::GetSuperModuleNumber(Int_t absId)  const
@@ -861,36 +861,36 @@ Int_t  AliEMCALGeometry::GetSuperModuleNumber(Int_t absId)  const
   // Return the number of the  supermodule given the absolute
   // ALICE numbering id
 
-  static Int_t nSupMod, nTower, nIphi, nIeta;
-  GetCellIndex(absId, nSupMod, nTower, nIphi, nIeta);
+  static Int_t nSupMod, nModule, nIphi, nIeta;
+  GetCellIndex(absId, nSupMod, nModule, nIphi, nIeta);
   return nSupMod;
 } 
 
 void  AliEMCALGeometry::GetModuleIndexesFromCellIndexesInSModule(Int_t nSupMod, Int_t iphi, Int_t ieta, 
-                       Int_t &iphim, Int_t &ietam, Int_t &nTower) const
+                       Int_t &iphim, Int_t &ietam, Int_t &nModule) const
 {
-  // Transition from cell indexes (ieta,iphi) to module indexes (ietam,iphim, nTower)
+  // Transition from cell indexes (ieta,iphi) to module indexes (ietam,iphim, nModule)
   static Int_t nphi;
   nphi  = GetNumberOfModuleInPhiDirection(nSupMod);  
 
   ietam  = ieta/fNETAdiv;
   iphim  = iphi/fNPHIdiv;
-  nTower = ietam * nphi + iphim; 
+  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
-  static Int_t ietam, iphim, nTower;
+  static Int_t ietam, iphim, nModule;
   static Int_t nIeta, nIphi; // cell indexes in module
 
-  GetModuleIndexesFromCellIndexesInSModule(nSupMod, iphi, ieta, ietam, iphim, nTower);
+  GetModuleIndexesFromCellIndexesInSModule(nSupMod, iphi, ieta, ietam, iphim, nModule);
 
   nIeta = ieta%fNETAdiv;
   nIeta = fNETAdiv - 1 - nIeta;
   nIphi = iphi%fNPHIdiv;
 
-  return GetAbsCellId(nSupMod, nTower, nIphi, nIeta);
+  return GetAbsCellId(nSupMod, nModule, nIphi, nIeta);
 }
 
 
@@ -909,11 +909,11 @@ Bool_t AliEMCALGeometry::RelPosCellInSModule(Int_t absId, Double_t &xr, Double_t
   // 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
-  static Int_t nSupMod, nTower, nIphi, nIeta, iphi, ieta;
+  static Int_t nSupMod, nModule, nIphi, nIeta, iphi, ieta;
   if(!CheckAbsCellId(absId)) return kFALSE;
 
-  GetCellIndex(absId, nSupMod, nTower, nIphi, nIeta);
-  GetCellPhiEtaIndexInSModule(nSupMod,nTower,nIphi,nIeta, iphi, ieta); 
+  GetCellIndex(absId, nSupMod, nModule, nIphi, nIeta);
+  GetCellPhiEtaIndexInSModule(nSupMod,nModule,nIphi,nIeta, iphi, ieta); 
  
   xr = fCentersOfCellsXDir.At(ieta);
   zr = fCentersOfCellsEtaDir.At(ieta);
@@ -991,7 +991,7 @@ void AliEMCALGeometry::CreateListOfTrd1Modules()
   // 
   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, nTower=0, iphiTemp;
+  Int_t iphi=0, ieta=0, nModule=0, iphiTemp;
   Double_t xr, zr, theta, phi, eta, r, x,y;
   TVector3 vglob;
   Double_t ytCenterModule, ytCenterCell;
@@ -1028,17 +1028,17 @@ void AliEMCALGeometry::CreateListOfTrd1Modules()
   AliDebug(2,Form(" Cells grid in eta directions : size %i\n", fCentersOfCellsEtaDir.GetSize()));
   for(Int_t it=0; it<fNZ; it++) {
     AliEMCALShishKebabTrd1Module *trd1 = GetShishKebabModule(it);
-    nTower = fNPhi*it;
+    nModule = fNPhi*it;
     for(Int_t ic=0; ic<fNETAdiv; ic++) {
       if(fNPHIdiv==2) {
         trd1->GetCenterOfCellInLocalCoordinateofSM(ic, xr, zr);      // case of 2X2
-        GetCellPhiEtaIndexInSModule(0, nTower, 0, ic, iphiTemp, ieta); 
+        GetCellPhiEtaIndexInSModule(0, nModule, 0, ic, iphiTemp, ieta); 
       } if(fNPHIdiv==3) {
         trd1->GetCenterOfCellInLocalCoordinateofSM_3X3(ic, xr, zr);  // case of 3X3
-        GetCellPhiEtaIndexInSModule(0, nTower, 0, ic, iphiTemp, ieta); 
+        GetCellPhiEtaIndexInSModule(0, nModule, 0, ic, iphiTemp, ieta); 
       } if(fNPHIdiv==1) {
         trd1->GetCenterOfCellInLocalCoordinateofSM_1X1(xr, zr);      // case of 1X1
-        GetCellPhiEtaIndexInSModule(0, nTower, 0, ic, iphiTemp, ieta); 
+        GetCellPhiEtaIndexInSModule(0, nModule, 0, ic, iphiTemp, ieta); 
       }
       fCentersOfCellsXDir.AddAt(float(xr) - fParSM[0],ieta);
       fCentersOfCellsEtaDir.AddAt(float(zr) - fParSM[2],ieta);
index db79b9e..c77227e 100644 (file)
@@ -142,17 +142,31 @@ public:
   TGeoMatrix *GetTransformationForSM(int i) {
   if(i>=0 && i<GetNumberOfSuperModules()) return fMatrixOfSM[i]; 
                                         else return 0;}
-  // abs id <-> indexes; Shish-kebab case, only TRD1 now.
-  // EMCAL -> Super Module -> module -> tower(or cell) - logic tree of EMCAL
   // May 31, 2006; ALICE numbering scheme: 
   // see ALICE-INT-2003-038: ALICE Coordinate System and Software Numbering Convention
   // All indexes are stared from zero now.
-  Int_t   GetAbsCellId(Int_t nSupMod, Int_t nTower, Int_t nIphi, Int_t nIeta) const;
+  // 
+  // abs id <-> indexes; Shish-kebab case, only TRD1 now.
+  // EMCAL -> Super Module -> module -> tower(or cell) - logic tree of EMCAL
+  // 
+  //**  Usual name of variable - Dec 18,2006 **
+  //  nSupMod - index of super module (SM)
+  //  nModule - index of module in SM
+  //  nIphi   - phi index of tower(cell) in module
+  //  nIeta   - eta index of tower(cell) in module
+  //  
+  //  Inside SM
+  //  iphim   - phi index of module in SM  
+  //  ietam   - eta index of module in SM  
+  //
+  //  iphi    - phi index of tower(cell) in SM  
+  //  ieta    - eta index of tower(cell) in SM  
+  Int_t   GetAbsCellId(Int_t nSupMod, Int_t nModule, Int_t nIphi, Int_t nIeta) const;
   Bool_t  CheckAbsCellId(Int_t absId) const;
-  Bool_t  GetCellIndex(Int_t absId, Int_t &nSupMod, Int_t &nTower, Int_t &nIphi, Int_t &nIeta) const;
+  Bool_t  GetCellIndex(Int_t absId, Int_t &nSupMod, Int_t &nModule, Int_t &nIphi, Int_t &nIeta) const;
   // Local coordinate of Super Module 
-  void    GetModulePhiEtaIndexInSModule(Int_t nSupMod, Int_t nTower, Int_t &iphim, Int_t &ietam) const;
-  void    GetCellPhiEtaIndexInSModule(Int_t nSupMod, Int_t nTower, Int_t nIphi, Int_t nIeta,
+  void    GetModulePhiEtaIndexInSModule(Int_t nSupMod, Int_t nModule, Int_t &iphim, Int_t &ietam) const;
+  void    GetCellPhiEtaIndexInSModule(Int_t nSupMod, Int_t nModule, Int_t nIphi, Int_t nIeta,
                                       Int_t &iphi, Int_t &ieta) const ;
   Int_t   GetSuperModuleNumber(Int_t absId)  const;
   Int_t   GetNumberOfModuleInPhiDirection(Int_t nSupMod)  const
@@ -163,7 +177,7 @@ public:
   }
   // From cell indexes to abs cell id
   void    GetModuleIndexesFromCellIndexesInSModule(Int_t nSupMod, Int_t iphi, Int_t ieta, 
-         Int_t &iphim, Int_t &ietam, Int_t &nTower) const;
+         Int_t &iphim, Int_t &ietam, Int_t &nModule) const;
   Int_t   GetAbsCellIdFromCellIndexes(Int_t nSupMod, Int_t iphi, Int_t ieta) const;
   // Methods for AliEMCALRecPoint - Feb 19, 2006
   Bool_t   RelPosCellInSModule(Int_t absId, Double_t &xr, Double_t &yr, Double_t &zr) const;
@@ -187,7 +201,7 @@ public:
   void SetSampling(Float_t samp) { fSampling = samp; printf("SetSampling: Sampling factor set to %f", fSampling) ; }
 
   Int_t GetNCellsInSupMod() const {return fNCellsInSupMod;}
-  Int_t GetNCellsInTower()  const {return fNCellsInTower; }
+  Int_t GetNCellsInModule()  const {return fNCellsInModule; }
   Int_t GetKey110DEG()      const {return fKey110DEG;}
 
   AliEMCALGeometry(); // default ctor only for internal usage (singleton)
@@ -244,7 +258,7 @@ private:
   //
   Int_t   fNCells;                       // number of cells in calo
   Int_t   fNCellsInSupMod;               // number cell in super module
-  Int_t   fNCellsInTower;                // number cell in tower(or module)
+  Int_t   fNCellsInModule;               // number cell in module)
   //TRU parameters
   Int_t   fNTRU ;                        //! Number of TRUs per module
   Int_t   fNTRUEta ;                     //! Number of cell rows per Z in one TRU
index 203a803..b64f125 100644 (file)
@@ -134,10 +134,10 @@ void AliEMCALGeometryOfflineTrd1::Init()
   // Fill fXYZofCells
   fXYZofCells = new TObjArray(fGeometry->GetNCells());
   fXYZofCells->SetName("CellsInGC"); 
-  Int_t nSupMod, nTower, nIphi, nIeta, iphi, ieta;
+  Int_t nSupMod, nModule, nIphi, nIeta, iphi, ieta;
   for(Int_t absId=1; absId<=fGeometry->GetNCells(); absId++){
-    if(fGeometry->GetCellIndex(absId, nSupMod,nTower,nIphi,nIeta)){
-      fGeometry->GetCellPhiEtaIndexInSModule(nSupMod,nTower,nIphi,nIeta, iphi,ieta);
+    if(fGeometry->GetCellIndex(absId, nSupMod,nModule,nIphi,nIeta)){
+      fGeometry->GetCellPhiEtaIndexInSModule(nSupMod,nModule,nIphi,nIeta, iphi,ieta);
       TVector3 *v = new TVector3;
       v->SetX(fSMPositionEta[ieta-1].Y()); 
       v->SetZ(fSMPositionEta[ieta-1].X()); 
@@ -149,13 +149,13 @@ void AliEMCALGeometryOfflineTrd1::Init()
 }
 
 //___________________________________________________________________________
-TVector3& AliEMCALGeometryOfflineTrd1::PosInSuperModule(int nSupMod, Int_t nTower, Int_t nIphi, Int_t nIeta)
+TVector3& AliEMCALGeometryOfflineTrd1::PosInSuperModule(int nSupMod, Int_t nModule, Int_t nIphi, Int_t nIeta)
 { 
   //return location of position within supermodule
   // 10-nov-04
   static Int_t iphi, ieta;
   static TVector3 v;
-  fGeometry->GetCellPhiEtaIndexInSModule(nSupMod, nTower,nIphi,nIeta, iphi,ieta);
+  fGeometry->GetCellPhiEtaIndexInSModule(nSupMod, nModule,nIphi,nIeta, iphi,ieta);
 
   // x-radius; y-phi; eta-z;
   v.SetXYZ(fSMPositionEta[ieta].Y(), fSMPositionPhi[iphi], fSMPositionEta[ieta].X());
@@ -177,13 +177,13 @@ double &lphi, double &leta)
 }
 
 //___________________________________________________________________________
-void AliEMCALGeometryOfflineTrd1::PositionInSuperModule(int nSupMod, Int_t nTower, Int_t nIphi, Int_t nIeta,
+void AliEMCALGeometryOfflineTrd1::PositionInSuperModule(int nSupMod, Int_t nModule, Int_t nIphi, Int_t nIeta,
 double &lphi, double &leta)
 {
   //return location of position within supermodule
 
   static Int_t iphi,ieta;
-  fGeometry->GetCellPhiEtaIndexInSModule(nSupMod, nTower,nIphi,nIeta,iphi,ieta);
+  fGeometry->GetCellPhiEtaIndexInSModule(nSupMod, nModule,nIphi,nIeta,iphi,ieta);
   PositionInSuperModule(iphi,ieta, lphi,leta);
 }
 
@@ -226,12 +226,12 @@ void AliEMCALGeometryOfflineTrd1::PrintSuperModule()
 void AliEMCALGeometryOfflineTrd1::PrintCell(Int_t absId)
 {
   //utility method for printing cell info
-  Int_t nSupMod, nTower, nIphi, nIeta, iphi, ieta;
-  if(fGeometry->GetCellIndex(absId, nSupMod,nTower,nIphi,nIeta)){
-     fGeometry->GetCellPhiEtaIndexInSModule(nSupMod,nTower,nIphi,nIeta, iphi,ieta);
+  Int_t nSupMod, nModule, nIphi, nIeta, iphi, ieta;
+  if(fGeometry->GetCellIndex(absId, nSupMod,nModule,nIphi,nIeta)){
+     fGeometry->GetCellPhiEtaIndexInSModule(nSupMod,nModule,nIphi,nIeta, iphi,ieta);
      TVector3 *v = CellPosition(absId);
      printf("(%5i) X %8.3f Y %8.3f Z %8.3f | #sup.Mod %2i #tower %3i nIphi %1i nIeta %1i | iphi %2i ieta %2i\n",
-           absId, v->X(),v->Y(),v->Z(), nSupMod,nTower,nIphi,nIeta, iphi,ieta);
+           absId, v->X(),v->Y(),v->Z(), nSupMod,nModule,nIphi,nIeta, iphi,ieta);
   } else {
     Warning("PrintCell","Wrong abs id %i\n",absId); 
   }
index df42d83..4241933 100644 (file)
@@ -33,12 +33,12 @@ class AliEMCALGeometryOfflineTrd1 : public TNamed {
   virtual ~AliEMCALGeometryOfflineTrd1() { /* nothing */ };
   static   AliEMCALGeometryOfflineTrd1* GetInstance();
   // positon in SuperModule
-  TVector3&  PosInSuperModule(int nSupMod, int nTower, int nIphi, int nIeta); 
+  TVector3&  PosInSuperModule(int nSupMod, int nModule, int nIphi, int nIeta); 
 
  public:
   // One Super Module
   void PositionInSuperModule(int iphi, int ieta, double &lphi, double &leta);
-  void PositionInSuperModule(int nSupMod, int nTower, int nIphi, int nIeta, double &lphi, double &leta);
+  void PositionInSuperModule(int nSupMod, int nModule, int nIphi, int nIeta, double &lphi, double &leta);
   // Position towers(cells)
   TVector3* CellPosition(int absId); // from 0 to fGeometry->GetNCells()
   // Global System
index 2ceb877..c87cc39 100644 (file)
@@ -220,18 +220,18 @@ Bool_t AliEMCALRecPoint::AreNeighbours(AliEMCALDigit * digit1, AliEMCALDigit * d
   // A neighbour is defined as being two digits which share a corner
   
   static Bool_t areNeighbours = kFALSE ;
-  static Int_t nSupMod=0, nTower=0, nIphi=0, nIeta=0;
-  static int nSupMod1=0, nTower1=0, nIphi1=0, nIeta1=0;
+  static Int_t nSupMod=0, nModule=0, nIphi=0, nIeta=0;
+  static int nSupMod1=0, nModule1=0, nIphi1=0, nIeta1=0;
   static Int_t relid1[2] , relid2[2] ; // ieta, iphi
   static Int_t rowdiff=0, coldiff=0;
 
   areNeighbours = kFALSE ;
 
-  fGeomPtr->GetCellIndex(digit1->GetId(), nSupMod,nTower,nIphi,nIeta);
-  fGeomPtr->GetCellPhiEtaIndexInSModule(nSupMod,nTower,nIphi,nIeta, relid1[0],relid1[1]);
+  fGeomPtr->GetCellIndex(digit1->GetId(), nSupMod,nModule,nIphi,nIeta);
+  fGeomPtr->GetCellPhiEtaIndexInSModule(nSupMod,nModule,nIphi,nIeta, relid1[0],relid1[1]);
 
-  fGeomPtr->GetCellIndex(digit2->GetId(), nSupMod1,nTower1,nIphi1,nIeta1);
-  fGeomPtr->GetCellPhiEtaIndexInSModule(nSupMod1,nTower1,nIphi1,nIeta1, relid2[0],relid2[1]);
+  fGeomPtr->GetCellIndex(digit2->GetId(), nSupMod1,nModule1,nIphi1,nIeta1);
+  fGeomPtr->GetCellPhiEtaIndexInSModule(nSupMod1,nModule1,nIphi1,nIeta1, relid2[0],relid2[1]);
   
   rowdiff = TMath::Abs( relid1[0] - relid2[0] ) ;  
   coldiff = TMath::Abs( relid1[1] - relid2[1] ) ;  
@@ -414,15 +414,15 @@ void  AliEMCALRecPoint::EvalDispersion(Float_t logWeight, TClonesArray * digits)
   
   // Calculates the dispersion in cell units 
   Double_t etai, phii, etaMean=0.0, phiMean=0.0; 
-  int nSupMod=0, nTower=0, nIphi=0, nIeta=0;
+  int nSupMod=0, nModule=0, nIphi=0, nIeta=0;
   int iphi=0, ieta=0;
   // Calculate mean values
   for(iDigit=0; iDigit < fMulDigit; iDigit++) {
     digit = (AliEMCALDigit *) digits->At(fDigitsList[iDigit])  ;
 
     if (fAmp>0 && fEnergyList[iDigit]>0) {
-      fGeomPtr->GetCellIndex(digit->GetId(), nSupMod,nTower,nIphi,nIeta);
-      fGeomPtr->GetCellPhiEtaIndexInSModule(nSupMod,nTower,nIphi,nIeta, iphi,ieta);
+      fGeomPtr->GetCellIndex(digit->GetId(), nSupMod,nModule,nIphi,nIeta);
+      fGeomPtr->GetCellPhiEtaIndexInSModule(nSupMod,nModule,nIphi,nIeta, iphi,ieta);
       etai=(Double_t)ieta;
       phii=(Double_t)iphi;
       w = TMath::Max(0.,logWeight+TMath::Log(fEnergyList[iDigit]/fAmp ) ) ;
@@ -444,8 +444,8 @@ void  AliEMCALRecPoint::EvalDispersion(Float_t logWeight, TClonesArray * digits)
     digit = (AliEMCALDigit *) digits->At(fDigitsList[iDigit])  ;
 
     if (fAmp>0 && fEnergyList[iDigit]>0) {
-      fGeomPtr->GetCellIndex(digit->GetId(), nSupMod,nTower,nIphi,nIeta);
-      fGeomPtr->GetCellPhiEtaIndexInSModule(nSupMod,nTower,nIphi,nIeta, iphi,ieta);
+      fGeomPtr->GetCellIndex(digit->GetId(), nSupMod,nModule,nIphi,nIeta);
+      fGeomPtr->GetCellPhiEtaIndexInSModule(nSupMod,nModule,nIphi,nIeta, iphi,ieta);
       etai=(Double_t)ieta;
       phii=(Double_t)iphi;
       w = TMath::Max(0.,logWeight+TMath::Log(fEnergyList[iDigit]/fAmp ) ) ;
@@ -574,7 +574,7 @@ void  AliEMCALRecPoint::EvalElipsAxis(Float_t logWeight,TClonesArray * digits)
   AliEMCALDigit * digit = 0;
 
   Double_t etai , phii, w; 
-  int nSupMod=0, nTower=0, nIphi=0, nIeta=0;
+  int nSupMod=0, nModule=0, nIphi=0, nIeta=0;
   int iphi=0, ieta=0;
   for(Int_t iDigit=0; iDigit<fMulDigit; iDigit++) {
     digit = (AliEMCALDigit *) digits->At(fDigitsList[iDigit])  ;
@@ -583,10 +583,10 @@ void  AliEMCALRecPoint::EvalElipsAxis(Float_t logWeight,TClonesArray * digits)
     // Nov 15,2006 - use cell numbers as coordinates
     // Copied for shish-kebab geometry, ieta,iphi is cast as double as eta,phi
     // We can use the eta,phi(or coordinates) of cell
-      nSupMod = nTower = nIphi = nIeta = iphi = ieta = 0;
+      nSupMod = nModule = nIphi = nIeta = iphi = ieta = 0;
 
-      fGeomPtr->GetCellIndex(digit->GetId(), nSupMod,nTower,nIphi,nIeta);
-      fGeomPtr->GetCellPhiEtaIndexInSModule(nSupMod,nTower,nIphi,nIeta, iphi,ieta);
+      fGeomPtr->GetCellIndex(digit->GetId(), nSupMod,nModule,nIphi,nIeta);
+      fGeomPtr->GetCellPhiEtaIndexInSModule(nSupMod,nModule,nIphi,nIeta, iphi,ieta);
       etai=(Double_t)ieta;
       phii=(Double_t)iphi;
     } else { // 
index 6038a3f..fb93fe3 100644 (file)
@@ -496,7 +496,7 @@ void AliEMCALv2::DrawAlicWithHits(int mode)
   }
 
   TClonesArray *hits = Hits();
-  Int_t nhits = hits->GetEntries(), absId, nSupMod, nTower, nIphi, nIeta, iphi, ieta;
+  Int_t nhits = hits->GetEntries(), absId, nSupMod, nModule, nIphi, nIeta, iphi, ieta;
   AliEMCALHit *hit = 0;
   Double_t de, des=0.;
   for(Int_t i=0; i<nhits; i++) {
@@ -504,11 +504,11 @@ void AliEMCALv2::DrawAlicWithHits(int mode)
     absId = hit->GetId();
     de    = hit->GetEnergy();
     des += de;
-    if(fGeometry->GetCellIndex(absId, nSupMod, nTower, nIphi, nIeta)){
+    if(fGeometry->GetCellIndex(absId, nSupMod, nModule, nIphi, nIeta)){
       //      printf(" de %f abs id %i smod %i tower %i | cell iphi %i : ieta %i\n",
-      // de, absId, nSupMod, nTower, nIphi, nIeta); 
+      // de, absId, nSupMod, nModule, nIphi, nIeta); 
       if(nSupMod==3) {
-        fGeometry->GetCellPhiEtaIndexInSModule(nSupMod,nTower,nIphi,nIeta, iphi,ieta);
+        fGeometry->GetCellPhiEtaIndexInSModule(nSupMod,nModule,nIphi,nIeta, iphi,ieta);
         // printf(" iphi %i : ieta %i\n", iphi,ieta);
         h2->Fill(double(ieta),double(iphi), de);
       }
@@ -540,18 +540,18 @@ void AliEMCALv2::TestIndexTransition(int pri, int idmax)
     return; 
   }
 
-  Int_t nSupMod, nTower, nIphi, nIeta, idNew, nGood=0;
+  Int_t nSupMod, nModule, nIphi, nIeta, idNew, nGood=0;
   if(idmax==0) idmax = fGeometry->GetNCells();
   for(Int_t id=1; id<=idmax; id++) {
-    if(!fGeometry->GetCellIndex(id, nSupMod, nTower, nIphi, nIeta)){
+    if(!fGeometry->GetCellIndex(id, nSupMod, nModule, nIphi, nIeta)){
       printf(" Wrong abs ID %i : #cells %i\n", id, fGeometry->GetNCells());
       break;  
     }
-    idNew = fGeometry->GetAbsCellId(nSupMod, nTower, nIphi, nIeta);
+    idNew = fGeometry->GetAbsCellId(nSupMod, nModule, nIphi, nIeta);
     if(id != idNew || pri>0) {
       printf(" ID %i : %i <- new id\n", id, idNew);
       printf(" nSupMod   %i ",  nSupMod);
-      printf(" nTower    %i ", nTower);
+      printf(" nModule    %i ", nModule);
       printf(" nIphi     %i ", nIphi);
       printf(" nIeta     %i \n", nIeta);