esd_tracks.C

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

diff --git a/STEER/AliESDZDC.cxx b/STEER/AliESDZDC.cxx
index b1a8fa51068bcb6a6daa2f6726053b4bb236b530..d5cb31b61fb9621379c73db5ce63f9ddcd718452 100644 (file)
--- a/STEER/AliESDZDC.cxx
+++ b/STEER/AliESDZDC.cxx
@@ -27,6 +27,8 @@
 
 ClassImp(AliESDZDC)
 
+//Float_t funCParCentr(const Int_t n) {return 1.8936-0.7126/(n + 0.7179);}
+
 //______________________________________________________________________________
 AliESDZDC::AliESDZDC() :
   TObject(),
@@ -36,16 +38,21 @@ AliESDZDC::AliESDZDC() :
   fZDCP2Energy(0),
   fZDCEMEnergy(0),
   fZDCEMEnergy1(0),
-  fZDCParticipants(0)
+  fZDCParticipants(0),
+  fZDCParticipants2(0)
 {
-  for(int i=0; i<5; i++){
+  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),
@@ -54,10 +61,11 @@ AliESDZDC::AliESDZDC(const AliESDZDC& zdc) :
   fZDCP2Energy(zdc.fZDCP2Energy),
   fZDCEMEnergy(zdc.fZDCEMEnergy),
   fZDCEMEnergy1(zdc.fZDCEMEnergy1),
-  fZDCParticipants(zdc.fZDCParticipants)
+  fZDCParticipants(zdc.fZDCParticipants),
+  fZDCParticipants2(zdc.fZDCParticipants2)
 {
   // copy constructor
-  for(int i=0; i<5; i++){
+  for(Int_t i=0; i<5; i++){
      fZN1TowerEnergy[i] = zdc.fZN1TowerEnergy[i];
      fZN2TowerEnergy[i] = zdc.fZN2TowerEnergy[i];
      fZP1TowerEnergy[i] = zdc.fZP1TowerEnergy[i];
@@ -67,8 +75,13 @@ AliESDZDC::AliESDZDC(const AliESDZDC& zdc) :
      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
@@ -79,6 +92,7 @@ AliESDZDC& AliESDZDC::operator=(const AliESDZDC&zdc)
     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++){
@@ -91,10 +105,28 @@ AliESDZDC& AliESDZDC::operator=(const AliESDZDC&zdc)
        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];
+    }
   } 
   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()
@@ -105,14 +137,18 @@ void AliESDZDC::Reset()
   fZDCN2Energy=0;
   fZDCP2Energy=0;
   fZDCParticipants=0;  
+  fZDCParticipants2=0;  
   fZDCEMEnergy=0;
   fZDCEMEnergy1=0;
-  for(int i=0; i<4; i++){
+  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.;
+  }
 }
 
 //______________________________________________________________________________
@@ -120,8 +156,51 @@ 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, Nparticipants = %d\n",
-  fZDCN1Energy,fZDCP1Energy,fZDCN2Energy,fZDCP2Energy,fZDCParticipants);
+  " ZP2Energy = %f, Npart side C = %d, Npart side A = %d\n",
+  fZDCN1Energy,fZDCP1Energy,fZDCN2Energy,fZDCP2Energy,fZDCParticipants,fZDCParticipants2);
 }
 
+//______________________________________________________________________________
+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;
+  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){
+    fZNCCentrCoord[0] = 1.8665*NumX/Den;
+    fZNCCentrCoord[1] = 1.8665*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.72, 1.75, 1.75};
+  Float_t NumX=0., NumY=0., Den=0.;
+  Float_t 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){
+    fZNACentrCoord[0] = 1.8665*NumX/Den;
+    fZNACentrCoord[1] = 1.8665*NumY/Den;
+  }
+  return fZNACentrCoord;
+}