promoted 2 methods to public needed by AMORE

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

diff --git a/STEER/AliESDZDC.cxx b/STEER/AliESDZDC.cxx
index d5cb31b61fb9621379c73db5ce63f9ddcd718452..f66d79d2a12ba65419e9263c9724f303ad2ca6f3 100644 (file)
--- a/STEER/AliESDZDC.cxx
+++ b/STEER/AliESDZDC.cxx
@@ -27,8 +27,6 @@
 
 ClassImp(AliESDZDC)
 
-//Float_t funCParCentr(const Int_t n) {return 1.8936-0.7126/(n + 0.7179);}
-
 //______________________________________________________________________________
 AliESDZDC::AliESDZDC() :
   TObject(),
@@ -39,7 +37,12 @@ AliESDZDC::AliESDZDC() :
   fZDCEMEnergy(0),
   fZDCEMEnergy1(0),
   fZDCParticipants(0),
-  fZDCParticipants2(0)
+  fZDCPartSideA(0),
+  fZDCPartSideC(0),
+  fImpactParameter(0),
+  fImpactParamSideA(0),
+  fImpactParamSideC(0),
+  fESDQuality(0)
 {
   for(Int_t i=0; i<5; i++){
     fZN1TowerEnergy[i] = fZN2TowerEnergy[i] = 0.;
@@ -50,6 +53,7 @@ AliESDZDC::AliESDZDC() :
   for(Int_t i=0; i<2; i++){
     fZNACentrCoord[i] = fZNCCentrCoord[i] = 0.;
   }
+  for(Int_t i=0; i<32; i++) fVMEScaler[i]=0;
 }
 
 //______________________________________________________________________________
@@ -62,7 +66,12 @@ AliESDZDC::AliESDZDC(const AliESDZDC& zdc) :
   fZDCEMEnergy(zdc.fZDCEMEnergy),
   fZDCEMEnergy1(zdc.fZDCEMEnergy1),
   fZDCParticipants(zdc.fZDCParticipants),
-  fZDCParticipants2(zdc.fZDCParticipants2)
+  fZDCPartSideA(zdc.fZDCPartSideA),
+  fZDCPartSideC(zdc.fZDCPartSideC),
+  fImpactParameter(zdc.fImpactParameter),
+  fImpactParamSideA(zdc.fImpactParamSideA),
+  fImpactParamSideC(zdc.fImpactParamSideC),
+  fESDQuality(zdc.fESDQuality)
 {
   // copy constructor
   for(Int_t i=0; i<5; i++){
@@ -79,6 +88,7 @@ AliESDZDC::AliESDZDC(const AliESDZDC& zdc) :
     fZNACentrCoord[i] = zdc.fZNACentrCoord[i];
     fZNCCentrCoord[i] = zdc.fZNCCentrCoord[i];
   }
+  for(Int_t i=0; i<32; i++) fVMEScaler[i] = zdc.fVMEScaler[i];
 }
 
 //______________________________________________________________________________
@@ -91,8 +101,6 @@ AliESDZDC& AliESDZDC::operator=(const AliESDZDC&zdc)
     fZDCP1Energy = zdc.fZDCP1Energy;
     fZDCN2Energy = zdc.fZDCN2Energy;
     fZDCP2Energy = zdc.fZDCP2Energy;
-    fZDCParticipants = zdc.fZDCParticipants;
-    fZDCParticipants2 = zdc.fZDCParticipants2;
     fZDCEMEnergy = zdc.fZDCEMEnergy;
     fZDCEMEnergy1 = zdc.fZDCEMEnergy1;
     for(Int_t i=0; i<5; i++){
@@ -105,10 +113,21 @@ AliESDZDC& AliESDZDC::operator=(const AliESDZDC&zdc)
        fZP1TowerEnergyLR[i] = zdc.fZP1TowerEnergyLR[i];
        fZP2TowerEnergyLR[i] = zdc.fZP2TowerEnergyLR[i];
     }
+    //
+    fZDCParticipants = zdc.fZDCParticipants;
+    fZDCPartSideA = zdc.fZDCPartSideA;
+    fZDCPartSideC = zdc.fZDCPartSideC;
+    fImpactParameter = zdc.fImpactParameter;
+    fImpactParamSideA = zdc.fImpactParamSideA;
+    fImpactParamSideC = zdc.fImpactParamSideC;
+    //
     for(Int_t i=0; i<2; i++){
          fZNACentrCoord[i] = zdc.fZNACentrCoord[i];
          fZNCCentrCoord[i] = zdc.fZNCCentrCoord[i];
     }
+    //
+    fESDQuality = zdc.fESDQuality;
+    for(Int_t i=0; i<32; i++) fVMEScaler[i] = zdc.fVMEScaler[i];
   } 
   return *this;
 }
@@ -136,8 +155,6 @@ void AliESDZDC::Reset()
   fZDCP1Energy=0;
   fZDCN2Energy=0;
   fZDCP2Energy=0;
-  fZDCParticipants=0;  
-  fZDCParticipants2=0;  
   fZDCEMEnergy=0;
   fZDCEMEnergy1=0;
   for(Int_t i=0; i<5; i++){
@@ -146,18 +163,31 @@ void AliESDZDC::Reset()
     fZN1TowerEnergyLR[i] = fZN2TowerEnergyLR[i] = 0.;
     fZP1TowerEnergyLR[i] = fZP2TowerEnergyLR[i] = 0.;
   }
+  fZDCParticipants=0;  
+  fZDCPartSideA=0;  
+  fZDCPartSideC=0;  
+  fImpactParameter=0;
+  fImpactParamSideA=0;
+  fImpactParamSideC=0;
   for(Int_t i=0; i<2; i++){
        fZNACentrCoord[i] = fZNCCentrCoord[i] = 0.;
   }
+  fESDQuality=0;
+  for(Int_t i=0; i<32; i++) fVMEScaler[i] = 0;
 }
 
 //______________________________________________________________________________
 void AliESDZDC::Print(const Option_t *) const
 {
   //  Print ESD for the ZDC
-  printf("\n \t ZN1Energy = %f TeV, ZP1Energy = %f TeV, ZN2Energy = %f TeV,"
-  " ZP2Energy = %f, Npart side C = %d, Npart side A = %d\n",
-  fZDCN1Energy,fZDCP1Energy,fZDCN2Energy,fZDCP2Energy,fZDCParticipants,fZDCParticipants2);
+  printf("\n \t E_{ZNC} = %f TeV, E_{ZPC} = %f TeV, E_{ZNA} = %f TeV, E_{ZPA} = %f TeV,"
+  " E_{ZEM} = %f GeV, Npart = %d, b = %1.2f fm\n",
+  fZDCN1Energy/1000.,fZDCP1Energy/1000.,fZDCN2Energy/1000.,fZDCP2Energy/1000.,
+  fZDCEMEnergy+fZDCEMEnergy1, fZDCParticipants,fImpactParameter);
+  //
+  printf(" ### fVMEScaler: \n");
+  for(Int_t i=0; i<32; i++) printf("\t datum %d: %d \n",i,fVMEScaler[i]);
+  printf("\n");
 }
 
 //______________________________________________________________________________
@@ -165,19 +195,23 @@ Double32_t * AliESDZDC::GetZNCCentroid()
 {
   // Provide coordinates of centroid over ZN (side C) front face
   Float_t x[4] = {-1.75, 1.75, -1.75, 1.75};
-  Float_t y[4] = {-1.75, -1.72, 1.75, 1.75};
-  Float_t NumX=0., NumY=0., Den=0.;
-  Float_t alpha=0.395, w;
+  Float_t y[4] = {-1.75, -1.75, 1.75, 1.75};
+  Float_t numX=0., numY=0., den=0.;
+  Float_t c, alpha=0.395, w;
+  //
   for(Int_t i=0; i<4; i++){
     if(fZN1TowerEnergy[i+1]<0.) fZN1TowerEnergy[i+1]=0.;
     w = TMath::Power(fZN1TowerEnergy[i+1], alpha);
-    NumX += x[i]*w;
-    NumY += y[i]*w;
-    Den += w;
+    numX += x[i]*w;
+    numY += y[i]*w;
+    den += w;
   }
-  if(Den!=0){
-    fZNCCentrCoord[0] = 1.8665*NumX/Den;
-    fZNCCentrCoord[1] = 1.8665*NumY/Den;
+  if(den!=0){
+    // ATTENTION! Needs to be changed if E_beam(A-A) != 2.76 A TeV !!!!
+    Float_t nSpecn = fZDCN1Energy/2760.;
+    c = 1.89358-0.71262/(nSpecn+0.71789);
+    fZNCCentrCoord[0] = c*numX/den;
+    fZNCCentrCoord[1] = c*numY/den;
   }
   return fZNCCentrCoord;
 }
@@ -187,20 +221,24 @@ Double32_t * AliESDZDC::GetZNACentroid()
 {
   // Provide coordinates of centroid over ZN (side A) front face
   Float_t x[4] = {-1.75, 1.75, -1.75, 1.75};
-  Float_t y[4] = {-1.75, -1.72, 1.75, 1.75};
-  Float_t NumX=0., NumY=0., Den=0.;
-  Float_t alpha=0.395, w;
+  Float_t y[4] = {-1.75, -1.75, 1.75, 1.75};
+  Float_t numX=0., numY=0., den=0.;
+  Float_t c, alpha=0.395, w;
+
   for(Int_t i=0; i<4; i++){
     if(fZN2TowerEnergy[i+1]<0.) fZN2TowerEnergy[i+1]=0.;
     w = TMath::Power(fZN2TowerEnergy[i+1], alpha);
-    NumX += x[i]*w;
-    NumY += y[i]*w;
-    Den += w;
+    numX += x[i]*w;
+    numY += y[i]*w;
+    den += w;
   }
   //
-  if(Den!=0){
-    fZNACentrCoord[0] = 1.8665*NumX/Den;
-    fZNACentrCoord[1] = 1.8665*NumY/Den;
+  if(den!=0){
+    // ATTENTION! Needs to be changed if E_beam(A-A) != 2.76 A TeV !!!!
+    Float_t nSpecn = fZDCN2Energy/2760.;
+    c = 1.89358-0.71262/(nSpecn+0.71789);
+    fZNACentrCoord[0] = c*numX/den;
+    fZNACentrCoord[1] = c*numY/den;
   }
   return fZNACentrCoord;
 }