Bug fix in the copy constr of the esd friend. Noticed while running threaded online...

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

diff --git a/STEER/AliESDZDC.cxx b/STEER/AliESDZDC.cxx
index dc431b9ed5f1a676cf25efaa58d1e307e23291a6..bf0dc943f36ec23e7a9bb29188a4d79a5ab8068f 100644 (file)
--- a/STEER/AliESDZDC.cxx
+++ b/STEER/AliESDZDC.cxx
@@ -35,11 +35,27 @@ AliESDZDC::AliESDZDC() :
   fZDCN2Energy(0),
   fZDCP2Energy(0),
   fZDCEMEnergy(0),
-  fZDCParticipants(0)
+  fZDCEMEnergy1(0),
+  fZDCParticipants(0),
+  fZDCPartSideA(0),
+  fZDCPartSideC(0),
+  fImpactParameter(0),
+  fImpactParamSideA(0),
+  fImpactParamSideC(0),
+  fESDQuality(0)
 {
-  for(int i = 0;i<4;i++) fZN1TowerEnergy[i] = fZN2TowerEnergy[i] = 0.;
+  for(Int_t i=0; i<5; i++){
+    fZN1TowerEnergy[i] = fZN2TowerEnergy[i] = 0.;
+    fZP1TowerEnergy[i] = fZP2TowerEnergy[i] = 0.;
+    fZN1TowerEnergyLR[i] = fZN2TowerEnergyLR[i] = 0.;
+    fZP1TowerEnergyLR[i] = fZP2TowerEnergyLR[i] = 0.;
+  }
+  for(Int_t i=0; i<2; i++){
+    fZNACentrCoord[i] = fZNCCentrCoord[i] = 0.;
+  }
 }
 
+//______________________________________________________________________________
 AliESDZDC::AliESDZDC(const AliESDZDC& zdc) :
   TObject(zdc),
   fZDCN1Energy(zdc.fZDCN1Energy),
@@ -47,11 +63,33 @@ AliESDZDC::AliESDZDC(const AliESDZDC& zdc) :
   fZDCN2Energy(zdc.fZDCN2Energy),
   fZDCP2Energy(zdc.fZDCP2Energy),
   fZDCEMEnergy(zdc.fZDCEMEnergy),
-  fZDCParticipants(zdc.fZDCParticipants)
+  fZDCEMEnergy1(zdc.fZDCEMEnergy1),
+  fZDCParticipants(zdc.fZDCParticipants),
+  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++){
+     fZN1TowerEnergy[i] = zdc.fZN1TowerEnergy[i];
+     fZN2TowerEnergy[i] = zdc.fZN2TowerEnergy[i];
+     fZP1TowerEnergy[i] = zdc.fZP1TowerEnergy[i];
+     fZP2TowerEnergy[i] = zdc.fZP2TowerEnergy[i];
+     fZN1TowerEnergyLR[i] = zdc.fZN1TowerEnergyLR[i];
+     fZN2TowerEnergyLR[i] = zdc.fZN2TowerEnergyLR[i];
+     fZP1TowerEnergyLR[i] = zdc.fZP1TowerEnergyLR[i];
+     fZP2TowerEnergyLR[i] = zdc.fZP2TowerEnergyLR[i];
+  }
+  for(Int_t i=0; i<2; i++){
+    fZNACentrCoord[i] = zdc.fZNACentrCoord[i];
+    fZNCCentrCoord[i] = zdc.fZNCCentrCoord[i];
+  }
 }
 
+//______________________________________________________________________________
 AliESDZDC& AliESDZDC::operator=(const AliESDZDC&zdc)
 {
   // assigment operator
@@ -62,15 +100,49 @@ AliESDZDC& AliESDZDC::operator=(const AliESDZDC&zdc)
     fZDCN2Energy = zdc.fZDCN2Energy;
     fZDCP2Energy = zdc.fZDCP2Energy;
     fZDCEMEnergy = zdc.fZDCEMEnergy;
-    fZDCParticipants = zdc.fZDCParticipants;
-    for(Int_t i=0; i<4; i++){
+    fZDCEMEnergy1 = zdc.fZDCEMEnergy1;
+    for(Int_t i=0; i<5; i++){
        fZN1TowerEnergy[i] = zdc.fZN1TowerEnergy[i];
        fZN2TowerEnergy[i] = zdc.fZN2TowerEnergy[i];
+       fZP1TowerEnergy[i] = zdc.fZP1TowerEnergy[i];
+       fZP2TowerEnergy[i] = zdc.fZP2TowerEnergy[i];
+       fZN1TowerEnergyLR[i] = zdc.fZN1TowerEnergyLR[i];
+       fZN2TowerEnergyLR[i] = zdc.fZN2TowerEnergyLR[i];
+       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;
   } 
   return *this;
 }
 
+//______________________________________________________________________________
+void AliESDZDC::Copy(TObject &obj) const {
+  
+  // this overwrites the virtual TOBject::Copy()
+  // to allow run time copying without casting
+  // in AliESDEvent
+
+  if(this==&obj)return;
+  AliESDZDC *robj = dynamic_cast<AliESDZDC*>(&obj);
+  if(!robj)return; // not an AliESDZDC
+  *robj = *this;
+
+}
+
 
 //______________________________________________________________________________
 void AliESDZDC::Reset()
@@ -81,14 +153,84 @@ void AliESDZDC::Reset()
   fZDCN2Energy=0;
   fZDCP2Energy=0;
   fZDCEMEnergy=0;
+  fZDCEMEnergy1=0;
+  for(Int_t i=0; i<5; i++){
+    fZN1TowerEnergy[i] = fZN2TowerEnergy[i] = 0.;
+    fZP1TowerEnergy[i] = fZP2TowerEnergy[i] = 0.;
+    fZN1TowerEnergyLR[i] = fZN2TowerEnergyLR[i] = 0.;
+    fZP1TowerEnergyLR[i] = fZP2TowerEnergyLR[i] = 0.;
+  }
   fZDCParticipants=0;  
-  for(int i=0; i<24; i++) fZN1TowerEnergy[i] = fZN2TowerEnergy[i] = 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;
 }
 
 //______________________________________________________________________________
 void AliESDZDC::Print(const Option_t *) const
 {
+  //  Print ESD for the ZDC
+  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);
 }
 
+//______________________________________________________________________________
+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.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;
+  }
+  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;
+}
 
+//______________________________________________________________________________
+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.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;
+  }
+  //
+  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;
+}