Bug #88017: Adding the equalization factors to AOD/VEvent.

[u/mrichter/AliRoot.git] / STEER / AOD / AliAODEvent.cxx

diff --git a/STEER/AOD/AliAODEvent.cxx b/STEER/AOD/AliAODEvent.cxx
index 6f9b7c071fb318ce4bb4cd73f94fe744287773ed..40e553f679ea988cc4b4a619470ab9f1b1022fbf 100644 (file)
--- a/STEER/AOD/AliAODEvent.cxx
+++ b/STEER/AOD/AliAODEvent.cxx
@@ -819,3 +819,20 @@ Bool_t AliAODEvent::IsPileupFromSPDInMultBins() const {
     else return IsPileupFromSPD(5,0.8);
 }
 
+Float_t AliAODEvent::GetVZEROEqMultiplicity(Int_t i) const
+{
+  // Get VZERO Multiplicity for channel i
+  // Themethod uses the equalization factors
+  // stored in the ESD-run object in order to
+  // get equal multiplicities within a VZERO rins (1/8 of VZERO)
+  if (!fAODVZERO || !fHeader) return -1;
+
+  Int_t ring = i/8;
+  Float_t factorSum = 0;
+  for(Int_t j = 8*ring; j < (8*ring+8); ++j) {
+    factorSum += fHeader->GetVZEROEqFactors(j);
+  }
+  Float_t factor = fHeader->GetVZEROEqFactors(i)*8./factorSum;
+
+  return (fAODVZERO->GetMultiplicity(i)/factor);
+}