AliCentrality for ESD and AOD analysis

[u/mrichter/AliRoot.git] / STEER / AliESDZDC.cxx

diff --git a/STEER/AliESDZDC.cxx b/STEER/AliESDZDC.cxx
index 02f40e0fdcc07a138b99ef35985da2cd2310c43f..c570602ffd88b6ba6d5d962102cb165b57a43836 100644 (file)
--- a/STEER/AliESDZDC.cxx
+++ b/STEER/AliESDZDC.cxx
@@ -55,7 +55,10 @@ AliESDZDC::AliESDZDC() :
   }
   for(Int_t i=0; i<32; i++){
     fVMEScaler[i]=0;
-    fZDCTDCData[i]=0;
+    for(Int_t y=0; y<4; y++){
+      fZDCTDCData[i][y]=0.;
+      fZDCTDCCorrected[i][y]=0.;
+    }
   }
 }
 
@@ -93,7 +96,10 @@ AliESDZDC::AliESDZDC(const AliESDZDC& zdc) :
   }
   for(Int_t i=0; i<32; i++){
     fVMEScaler[i] = zdc.fVMEScaler[i];
-    fZDCTDCData[i] = zdc.fZDCTDCData[i];
+    for(Int_t y=0; y<4; y++){
+       fZDCTDCData[i][y] = zdc.fZDCTDCData[i][y];
+       fZDCTDCCorrected[i][y] = zdc.fZDCTDCCorrected[i][y];
+    }
   }
 }
 
@@ -135,7 +141,10 @@ AliESDZDC& AliESDZDC::operator=(const AliESDZDC&zdc)
     fESDQuality = zdc.fESDQuality;
     for(Int_t i=0; i<32; i++){
       fVMEScaler[i] = zdc.fVMEScaler[i];
-      fZDCTDCData[i] = zdc.fZDCTDCData[i];
+      for(Int_t y=0; y<4; y++){ 
+         fZDCTDCData[i][y] = zdc.fZDCTDCData[i][y];
+         fZDCTDCCorrected[i][y] = zdc.fZDCTDCCorrected[i][y];
+      }
     }
   } 
   return *this;
@@ -184,7 +193,10 @@ void AliESDZDC::Reset()
   fESDQuality=0;
   for(Int_t i=0; i<32; i++){
      fVMEScaler[i] = 0;
-     fZDCTDCData[i] = 0;
+     for(Int_t y=0; y<4; y++){
+        fZDCTDCData[i][y] = 0.;
+        fZDCTDCCorrected[i][y] = 0.;
+     }
   }
 }
 
@@ -202,7 +214,10 @@ void AliESDZDC::Print(const Option_t *) const
   printf("\n");
   //
   printf(" ### TDCData (!=0): \n");
-  for(Int_t i=0; i<32; i++) if(fZDCTDCData[i]!=0) printf("\t %d \n",fZDCTDCData[i]);
+  for(Int_t i=0; i<32; i++){
+    for(Int_t j=0; j<4; j++)
+      if(fZDCTDCCorrected[i][j]!=0) printf("\t %1.0f \n",fZDCTDCCorrected[i][j]);
+  }
   printf("\n");
 }
 
@@ -230,7 +245,7 @@ Bool_t AliESDZDC::GetZNCentroidInPbPb(Float_t beamEne, Double_t centrZNC[2], Dou
       denZNC += wZNC;
     }
     if(fZN2TowerEnergy[i+1]>0.) {
-      wZNA = TMath::Power(fZN1TowerEnergy[i+1], alpha);
+      wZNA = TMath::Power(fZN2TowerEnergy[i+1], alpha);
       numXZNA += x[i]*wZNA;
       numYZNA += y[i]*wZNA;
       denZNA += wZNA;
@@ -247,10 +262,10 @@ Bool_t AliESDZDC::GetZNCentroidInPbPb(Float_t beamEne, Double_t centrZNC[2], Dou
     fZNCCentrCoord[0] = fZNCCentrCoord[1] = 999.;
   }
   if(denZNA!=0){
-    Float_t nSpecnA = fZDCN1Energy/beamEne;
+    Float_t nSpecnA = fZDCN2Energy/beamEne;
     cZNA = 1.89358-0.71262/(nSpecnA+0.71789);
-    fZNCCentrCoord[0] = cZNA*numXZNA/denZNA;
-    fZNCCentrCoord[1] = cZNA*numYZNA/denZNA;
+    fZNACentrCoord[0] = cZNA*numXZNA/denZNA;
+    fZNACentrCoord[1] = cZNA*numYZNA/denZNA;
   } 
   else{
     fZNACentrCoord[0] = fZNACentrCoord[1] = 999.;