STEER/ESD/AliESDZDC.h

   1 // -*- mode: C++ -*-
   2 #ifndef ALIESDZDC_H
   3 #define ALIESDZDC_H
   4
   5 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   6  * See cxx source for full Copyright notice                               */
   7
   8 //-------------------------------------------------------------------------
   9 //                      Implementation of   Class AliESDZDC
  10 //   This is a class that summarizes the ZDC data for the ESD
  11 //   Origin: Christian Klein-Boesing, CERN, Christian.Klein-Boesing@cern.ch
  12 //   *** 15/10/2009 Scaler added to AliESDZDC class ***
  13 //   *** Scaler format:  32 floats from VME scaler  ***
  14 //-------------------------------------------------------------------------
  15
  16 #include <AliVZDC.h>
  17
  18 class AliESDZDC: public AliVZDC {
  19 public:
  20
  21   AliESDZDC();
  22   AliESDZDC(const AliESDZDC& zdc);
  23   AliESDZDC& operator=(const AliESDZDC& zdc);
  24
  25
  26   // Getters
  27
  28   virtual Short_t  GetZDCParticipants() const {return fZDCParticipants;}
  29   virtual Short_t  GetZDCPartSideA() const {return fZDCPartSideA;}
  30   virtual Short_t  GetZDCPartSideC() const {return fZDCPartSideC;}
  31   virtual Double_t GetImpactParameter() const {return fImpactParameter;}
  32   virtual Double_t GetImpactParamSideA() const {return fImpactParamSideA;}
  33   virtual Double_t GetImpactParamSideC() const {return fImpactParamSideC;}
  34
  35   virtual Double_t GetZNCEnergy() const {return (Double_t) fZDCN1Energy;}
  36   virtual Double_t GetZPCEnergy() const {return (Double_t) fZDCP1Energy;}
  37   virtual Double_t GetZNAEnergy() const {return (Double_t) fZDCN2Energy;}
  38   virtual Double_t GetZPAEnergy() const {return (Double_t) fZDCP2Energy;}
  39   virtual Double_t GetZEM1Energy() const {return (Double_t) fZDCEMEnergy;}
  40   virtual Double_t GetZEM2Energy() const {return (Double_t) fZDCEMEnergy1;}
  41
  42   virtual const Double_t *GetZNCTowerEnergy() const {return fZN1TowerEnergy;}
  43   virtual const Double_t *GetZNATowerEnergy() const {return fZN2TowerEnergy;}
  44   virtual const Double_t *GetZPCTowerEnergy() const {return fZP1TowerEnergy;}
  45   virtual const Double_t *GetZPATowerEnergy() const {return fZP2TowerEnergy;}
  46   virtual const Double_t *GetZNCTowerEnergyLR() const {return fZN1TowerEnergyLR;}
  47   virtual const Double_t *GetZNATowerEnergyLR() const {return fZN2TowerEnergyLR;}
  48   const Double_t *GetZPCTowerEnergyLR() const {return fZP1TowerEnergyLR;}
  49   const Double_t *GetZPATowerEnergyLR() const {return fZP2TowerEnergyLR;}
  50
  51   UInt_t   GetESDQuality()  const {return fESDQuality;}
  52   Double_t GetZDCN1Energy() const {return fZDCN1Energy;}
  53   Double_t GetZDCP1Energy() const {return fZDCP1Energy;}
  54   Double_t GetZDCN2Energy() const {return fZDCN2Energy;}
  55   Double_t GetZDCP2Energy() const {return fZDCP2Energy;}
  56   Double_t GetZDCEMEnergy(Int_t i) const
  57            {if(i==0){return fZDCEMEnergy;} else if(i==1){return fZDCEMEnergy1;}
  58            return 0;}
  59
  60   const Double_t * GetZN1TowerEnergy() const {return fZN1TowerEnergy;}
  61   const Double_t * GetZN2TowerEnergy() const {return fZN2TowerEnergy;}
  62   const Double_t * GetZP1TowerEnergy() const {return fZP1TowerEnergy;}
  63   const Double_t * GetZP2TowerEnergy() const {return fZP2TowerEnergy;}
  64   const Double_t * GetZN1TowerEnergyLR() const {return fZN1TowerEnergyLR;}
  65   const Double_t * GetZN2TowerEnergyLR() const {return fZN2TowerEnergyLR;}
  66   const Double_t * GetZP1TowerEnergyLR() const {return fZP1TowerEnergyLR;}
  67   const Double_t * GetZP2TowerEnergyLR() const {return fZP2TowerEnergyLR;}
  68   //
  69   virtual Bool_t GetZNCentroidInPbPb(Float_t beamEne, Double_t centrZNC[2], Double_t centrZNA[2]);
  70   virtual Bool_t GetZNCentroidInpp(Double_t centrZNC[2], Double_t centrZNA[2]);
  71   //
  72
  73   UInt_t GetZDCScaler(Int_t i)  const {return fVMEScaler[i];}
  74   const UInt_t* GetZDCScaler()  const {return fVMEScaler;}
  75
  76   Int_t GetZDCTDCData(Int_t i, Int_t j) const {return fZDCTDCData[i][j];}
  77   Float_t GetZDCTDCCorrected(Int_t i, Int_t j) const
  78   {
  79     if(AliESDZDC::kCorrectedTDCFilled && (fZDCTDCData[i][j]!=0)) return fZDCTDCCorrected[i][j];
  80            else return 0.;
  81   }
  82
  83   Float_t GetZNTDCSum(Int_t ihit) const
  84   {
  85     if(ihit<4 && AliESDZDC::kCorrectedTDCFilled && (fZDCTDCData[10][ihit]!=0) && (fZDCTDCData[12][ihit]!=0)) return (Float_t) (fZDCTDCCorrected[10][ihit]+fZDCTDCCorrected[12][ihit]);
  86     else return 0.;
  87   }
  88   Float_t GetZNTDCDiff(Int_t ihit) const
  89   {
  90     if(ihit<4 && AliESDZDC::kCorrectedTDCFilled &&  (fZDCTDCData[10][ihit]!=0) && (fZDCTDCData[12][ihit]!=0)) return (Float_t) (fZDCTDCCorrected[12][ihit]-fZDCTDCCorrected[10][ihit]);
  91     else return 0.;
  92   }
  93
  94   virtual Float_t GetZDCTimeSum() const
  95           {if(AliESDZDC::kCorrectedTDCFilled && (fZDCTDCData[10][0]!=0) && (fZDCTDCData[12][0]!=0)) return (Float_t) (fZDCTDCCorrected[10][0]+fZDCTDCCorrected[12][0]);
  96            else return 0.;}
  97   virtual Float_t GetZDCTimeDiff() const
  98           {if(AliESDZDC::kCorrectedTDCFilled &&  (fZDCTDCData[10][0]!=0) && (fZDCTDCData[12][0]!=0)) return (Float_t) (fZDCTDCCorrected[12][0]-fZDCTDCCorrected[10][0]);
  99            else return 0.;}
 100
 101   void  SetZDC(Double_t n1Energy, Double_t p1Energy,
 102                Double_t emEnergy0, Double_t emEnergy1,
 103                Double_t n2Energy, Double_t p2Energy,
 104                Short_t participants, Short_t nPartA, Short_t nPartC,
 105                Double_t b, Double_t bA, Double_t bC, UInt_t recoFlag)
 106         {fZDCN1Energy=n1Energy; fZDCP1Energy=p1Energy;
 107          fZDCEMEnergy=emEnergy0; fZDCEMEnergy1=emEnergy1;
 108          fZDCN2Energy=n2Energy; fZDCP2Energy=p2Energy;
 109          fZDCParticipants=participants; fZDCPartSideA=nPartA; fZDCPartSideC=nPartC;
 110          fImpactParameter=b; fImpactParamSideA=bA, fImpactParamSideC=bC,
 111          fESDQuality=recoFlag;}
 112   //
 113   void  SetZN1TowerEnergy(const Float_t tow1[5])
 114           {for(Int_t i=0; i<5; i++) fZN1TowerEnergy[i] = tow1[i];}
 115   void  SetZN2TowerEnergy(const Float_t tow2[5])
 116           {for(Int_t i=0; i<5; i++) fZN2TowerEnergy[i] = tow2[i];}
 117   void  SetZP1TowerEnergy(const Float_t tow1[5])
 118           {for(Int_t i=0; i<5; i++) fZP1TowerEnergy[i] = tow1[i];}
 119   void  SetZP2TowerEnergy(const Float_t tow2[5])
 120           {for(Int_t i=0; i<5; i++) fZP2TowerEnergy[i] = tow2[i];}
 121   void  SetZN1TowerEnergyLR(const Float_t tow1[5])
 122           {for(Int_t i=0; i<5; i++) fZN1TowerEnergyLR[i] = tow1[i];}
 123   void  SetZN2TowerEnergyLR(const Float_t tow2[5])
 124           {for(Int_t i=0; i<5; i++) fZN2TowerEnergyLR[i] = tow2[i];}
 125   void  SetZP1TowerEnergyLR(const Float_t tow1[5])
 126           {for(Int_t i=0; i<5; i++) fZP1TowerEnergyLR[i] = tow1[i];}
 127   void  SetZP2TowerEnergyLR(const Float_t tow2[5])
 128           {for(Int_t i=0; i<5; i++) fZP2TowerEnergyLR[i] = tow2[i];}
 129   void  SetZNACentroid(const Float_t centrCoord[2])
 130            {for(Int_t i=0; i<2; i++) fZNACentrCoord[i] = centrCoord[i];}
 131   void  SetZNCCentroid(const Float_t centrCoord[2])
 132            {for(Int_t i=0; i<2; i++) fZNCCentrCoord[i] = centrCoord[i];}
 133
 134   void SetZDCScaler(const UInt_t count[32])
 135        {for(Int_t k=0; k<32; k++) fVMEScaler[k] = count[k];}
 136
 137   void SetZDCTDCData(const Int_t values[32][4])
 138        {for(Int_t k=0; k<32; k++)
 139            for(Int_t j=0; j<4; j++) fZDCTDCData[k][j] = values[k][j];}
 140
 141   void SetZDCTDCCorrected(const Float_t values[32][4])
 142        {for(Int_t k=0; k<32; k++)
 143            for(Int_t j=0; j<4; j++) fZDCTDCCorrected[k][j] = values[k][j];}
 144
 145   void    Reset();
 146   void    Print(const Option_t *opt=0) const;
 147 private:
 148   virtual void Copy(TObject &obj) const;
 149 private:
 150
 151   Double32_t   fZDCN1Energy;  // reconstructed energy in the neutron ZDC
 152   Double32_t   fZDCP1Energy;  // reconstructed energy in the proton ZDC
 153   Double32_t   fZDCN2Energy;  // reconstructed energy in the neutron ZDC
 154   Double32_t   fZDCP2Energy;  // reconstructed energy in the proton ZDC
 155   Double32_t   fZDCEMEnergy;  // signal in the electromagnetic ZDCs
 156   Double32_t   fZDCEMEnergy1; // second EM signal,cannot change fZDCEMEnergy to array (not backward compatible)
 157   Double32_t   fZN1TowerEnergy[5];// reco E in 5 ZN1 sectors - high gain chain
 158   Double32_t   fZN2TowerEnergy[5];// reco E in 5 ZN2 sectors - high gain chain
 159   Double32_t   fZP1TowerEnergy[5];// reco E in 5 ZP1 sectors - high gain chain
 160   Double32_t   fZP2TowerEnergy[5];// reco E in 5 ZP2 sectors - high gain chain
 161   Double32_t   fZN1TowerEnergyLR[5];// reco E in 5 ZN1 sectors - low gain chain
 162   Double32_t   fZN2TowerEnergyLR[5];// reco E in 5 ZN2 sectors - low gain chain
 163   Double32_t   fZP1TowerEnergyLR[5];// reco E in 5 ZP1 sectors - low gain chain
 164   Double32_t   fZP2TowerEnergyLR[5];// reco E in 5 ZP2 sectors - low gain chain
 165   Short_t      fZDCParticipants;    // number of participants estimated by the ZDC (ONLY in A-A)
 166   Short_t      fZDCPartSideA;     // number of participants estimated by the ZDC (ONLY in A-A)
 167   Short_t      fZDCPartSideC;     // number of participants estimated by the ZDC (ONLY in A-A)
 168   Double32_t   fImpactParameter;  // impact parameter estimated by the ZDC (ONLY in A-A)
 169   Double32_t   fImpactParamSideA; // impact parameter estimated by the ZDC (ONLY in A-A)
 170   Double32_t   fImpactParamSideC; // impact parameter estimated by the ZDC (ONLY in A-A)
 171   Double32_t   fZNACentrCoord[2]; // Coordinates of the centroid over ZNC
 172   Double32_t   fZNCCentrCoord[2]; // Coordinates of the centroid over ZNA
 173   UInt_t       fESDQuality;       // flags from reconstruction
 174   UInt_t       fVMEScaler[32];    // counts from VME scaler
 175   Int_t        fZDCTDCData[32][4];     // ZDC TDC data
 176   Float_t      fZDCTDCCorrected[32][4];// ZDC TDC data in ns corrected 4 phase shift
 177
 178   ClassDef(AliESDZDC,17)
 179 };
 180
 181 #endif
 182