update of Jenn and Marco

[u/mrichter/AliRoot.git] / EMCAL / AliEMCALRecPoint.cxx

diff --git a/EMCAL/AliEMCALRecPoint.cxx b/EMCAL/AliEMCALRecPoint.cxx
index 40469d822d7099a0871cb98ee5a50b9a6e148c62..db3d007b98ed5a4e99805ed37b41c3c2b29faa67 100644 (file)
--- a/EMCAL/AliEMCALRecPoint.cxx
+++ b/EMCAL/AliEMCALRecPoint.cxx
@@ -34,7 +34,6 @@
 #include "AliEMCALGeometry.h"
 #include "AliEMCALDigit.h"
 #include "AliEMCALRecPoint.h"
-#include "AliEMCALGetter.h"
 
 ClassImp(AliEMCALRecPoint)
 
@@ -131,7 +130,8 @@ Bool_t AliEMCALRecPoint::AreNeighbours(AliEMCALDigit * digit1, AliEMCALDigit * d
   
   Bool_t areNeighbours = kFALSE ;
   
-  AliEMCALGeometry * geom =  (AliEMCALGetter::Instance())->EMCALGeometry();
+  //AliEMCALGeometry * geom =  (AliEMCALGetter::Instance())->EMCALGeometry();
+  AliEMCALGeometry * geom =  AliEMCALGeometry::GetInstance();
 
   Int_t relid1[2] ; 
     //copied for shish-kebab geometry, ieta,iphi is cast as float as eta,phi conversion
@@ -139,7 +139,7 @@ Bool_t AliEMCALRecPoint::AreNeighbours(AliEMCALDigit * digit1, AliEMCALDigit * d
     int nSupMod=0, nTower=0, nIphi=0, nIeta=0;
     int iphi=0, ieta=0;
        geom->GetCellIndex(digit1->GetId(), nSupMod,nTower,nIphi,nIeta);
-       geom->GetCellPhiEtaIndexInSModule(nTower,nIphi,nIeta, iphi,ieta);
+       geom->GetCellPhiEtaIndexInSModule(nSupMod,nTower,nIphi,nIeta, iphi,ieta);
        relid1[0]=ieta;
        relid1[1]=iphi;
 //  geom->AbsToRelNumbering(digit1->GetId(), relid1) ; 
@@ -150,7 +150,7 @@ Bool_t AliEMCALRecPoint::AreNeighbours(AliEMCALDigit * digit1, AliEMCALDigit * d
     int nSupMod1=0, nTower1=0, nIphi1=0, nIeta1=0;
     int iphi1=0, ieta1=0;
        geom->GetCellIndex(digit2->GetId(), nSupMod1,nTower1,nIphi1,nIeta1);
-       geom->GetCellPhiEtaIndexInSModule(nTower1,nIphi1,nIeta1, iphi1,ieta1);
+       geom->GetCellPhiEtaIndexInSModule(nSupMod1,nTower1,nIphi1,nIeta1, iphi1,ieta1);
        relid2[0]=ieta1;
        relid2[1]=iphi1;
 //  geom->AbsToRelNumbering(digit2->GetId(), relid2) ; 
@@ -334,8 +334,9 @@ void  AliEMCALRecPoint::EvalDispersion(Float_t logWeight, TClonesArray * digits)
 
   AliEMCALDigit * digit ;
  
-  AliEMCALGeometry * geom = (AliEMCALGetter::Instance())->EMCALGeometry();
-  
+  //AliEMCALGeometry * geom = (AliEMCALGetter::Instance())->EMCALGeometry();
+  AliEMCALGeometry * geom =  AliEMCALGeometry::GetInstance();
+
   // Calculates the centre of gravity in the local EMCAL-module coordinates   
   Int_t iDigit;
 
@@ -362,7 +363,7 @@ void  AliEMCALRecPoint::EvalDispersion(Float_t logWeight, TClonesArray * digits)
     int nSupMod=0, nTower=0, nIphi=0, nIeta=0;
     int iphi=0, ieta=0;
        geom->GetCellIndex(digit->GetId(), nSupMod,nTower,nIphi,nIeta);
-       geom->GetCellPhiEtaIndexInSModule(nTower,nIphi,nIeta, iphi,ieta);
+       geom->GetCellPhiEtaIndexInSModule(nSupMod,nTower,nIphi,nIeta, iphi,ieta);
        etai=(Float_t)ieta;
        phii=(Float_t)iphi;
        //        printf("%f,%d,%d \n", fAmp, ieta, iphi) ;
@@ -401,7 +402,8 @@ void AliEMCALRecPoint::EvalLocalPosition(Float_t logWeight, TClonesArray * digit
   //  Int_t relid[3] ;
   
   AliEMCALDigit * digit ;
-  AliEMCALGeometry * geom  =  (AliEMCALGetter::Instance())->EMCALGeometry();
+  //AliEMCALGeometry * geom  =  (AliEMCALGetter::Instance())->EMCALGeometry();
+  AliEMCALGeometry * geom =  AliEMCALGeometry::GetInstance();
   Int_t iDigit;
   Float_t cluEta = 0;
   Float_t cluPhi = 0;
@@ -417,7 +419,7 @@ void AliEMCALRecPoint::EvalLocalPosition(Float_t logWeight, TClonesArray * digit
     int nSupMod=0, nTower=0, nIphi=0, nIeta=0;
     int iphi=0, ieta=0;
        geom->GetCellIndex(digit->GetId(), nSupMod,nTower,nIphi,nIeta);
-       geom->GetCellPhiEtaIndexInSModule(nTower,nIphi,nIeta, iphi,ieta);
+       geom->GetCellPhiEtaIndexInSModule(nSupMod, nTower,nIphi,nIeta, iphi,ieta); //19-oct-05
        etai=(Float_t)ieta;
        phii=(Float_t)iphi;
 //Sub    geom->EtaPhiFromIndex(digit->GetId(), etai, phii);
@@ -455,7 +457,9 @@ void AliEMCALRecPoint::EvalCoreEnergy(Float_t logWeight, TClonesArray * digits)
 
   AliEMCALDigit * digit ;
   const Float_t kDeg2Rad = TMath::DegToRad() ;
-  AliEMCALGeometry * geom = (AliEMCALGetter::Instance())->EMCALGeometry();    
+  //AliEMCALGeometry * geom = (AliEMCALGetter::Instance())->EMCALGeometry();    
+  AliEMCALGeometry * geom =  AliEMCALGeometry::GetInstance();
+
   Int_t iDigit;
 
   if (!fLocPos.X() || !fLocPos.Y() ) {
@@ -490,7 +494,8 @@ void  AliEMCALRecPoint::EvalElipsAxis(Float_t logWeight,TClonesArray * digits)
   const Float_t kDeg2Rad = TMath::DegToRad();
   AliEMCALDigit * digit ;
 
-  AliEMCALGeometry * geom = (AliEMCALGetter::Instance())->EMCALGeometry();
+  //AliEMCALGeometry * geom = (AliEMCALGetter::Instance())->EMCALGeometry();
+  AliEMCALGeometry * geom =  AliEMCALGeometry::GetInstance();
   TString gn(geom->GetName());
 
   Int_t iDigit;
@@ -505,7 +510,7 @@ void  AliEMCALRecPoint::EvalElipsAxis(Float_t logWeight,TClonesArray * digits)
       int nSupMod=0, nTower=0, nIphi=0, nIeta=0;
       int iphi=0, ieta=0;
       geom->GetCellIndex(digit->GetId(), nSupMod,nTower,nIphi,nIeta);
-      geom->GetCellPhiEtaIndexInSModule(nTower,nIphi,nIeta, iphi,ieta);
+      geom->GetCellPhiEtaIndexInSModule(nSupMod,nTower,nIphi,nIeta, iphi,ieta);
       etai=(Float_t)ieta;
       phii=(Float_t)iphi;
     } else {
@@ -676,7 +681,8 @@ void AliEMCALRecPoint::GetGlobalPosition(TVector3 & gpos) const
   // returns the position of the cluster in the global reference system of ALICE
   // These are now the Cartesian X, Y and Z
 
-  AliEMCALGeometry * geom = (AliEMCALGetter::Instance())->EMCALGeometry();
+  //AliEMCALGeometry * geom = (AliEMCALGetter::Instance())->EMCALGeometry();
+  AliEMCALGeometry * geom =  AliEMCALGeometry::GetInstance();
   Int_t absid = geom->TowerIndexFromEtaPhi(fLocPos.X(), TMath::RadToDeg()*fLocPos.Y());
   geom->XYZFromIndex(absid, gpos);
 }