Manual linearization is removed in case of pass2, i.e. saturation correction is done...

diff --git a/STEER/AliESDUtils.cxx b/STEER/AliESDUtils.cxx
index 75535f87a1c3bf2556fdc2a870e816f277578743..909d292021701315e5eaa54e01cf23b5eb15de4c 100644 (file)
--- a/STEER/AliESDUtils.cxx
+++ b/STEER/AliESDUtils.cxx
@@ -119,10 +119,15 @@ Float_t AliESDUtils::GetCorrV0(const AliESDEvent* esd, Float_t &v0CorrResc, Floa
   Float_t multChCorrResc[64];
   AliESDVZERO* esdV0 = esd->GetVZEROData();
   for(Int_t i = 0; i < 64; ++i) {
-    Double_t b = (esdV0->GetMultiplicity(i)*par1[i]-par0[i]);
-    Double_t s = (b*b-4.*par2[i]*esdV0->GetMultiplicity(i)*esdV0->GetMultiplicity(i));
-    Double_t n = (s<0) ? -b : (-b + TMath::Sqrt(s));
-    multChCorr[i] = 2.*esdV0->GetMultiplicity(i)/n*par0[i];
+    if (esdV0->TestBit(AliESDVZERO::kCorrectedForSaturation)) {
+      multChCorr[i] = esdV0->GetMultiplicity(i);
+    }
+    else {
+      Double_t b = (esdV0->GetMultiplicity(i)*par1[i]-par0[i]);
+      Double_t s = (b*b-4.*par2[i]*esdV0->GetMultiplicity(i)*esdV0->GetMultiplicity(i));
+      Double_t n = (s<0) ? -b : (-b + TMath::Sqrt(s));
+      multChCorr[i] = 2.*esdV0->GetMultiplicity(i)/n*par0[i];
+    }
     multCorr += multChCorr[i];
     multChCorrResc[i] = multChCorr[i]/par0[i]/64.;
     multCorr2 += multChCorrResc[i];