From 0caa916b85ea54b165b0284e603cd10593ef172d Mon Sep 17 00:00:00 2001
From: cvetan <cvetan@f7af4fe6-9843-0410-8265-dc069ae4e863>
Date: Thu, 16 Dec 2010 13:10:56 +0000
Subject: [PATCH 1/1] Restoring previous fixes that were lost during one of the
 last commits in this class. Some minor fixes + corrected V0 as a static
 method so that it can be used in the monitoring until pass2 of the HI data.

---
 ANALYSIS/AliCentralitySelectionTask.cxx | 16 ++++------------
 ANALYSIS/AliCentralitySelectionTask.h   |  2 +-
 2 files changed, 5 insertions(+), 13 deletions(-)

diff --git a/ANALYSIS/AliCentralitySelectionTask.cxx b/ANALYSIS/AliCentralitySelectionTask.cxx
index 5b06b7e9cc3..cd590253d23 100644
--- a/ANALYSIS/AliCentralitySelectionTask.cxx
+++ b/ANALYSIS/AliCentralitySelectionTask.cxx
@@ -357,8 +357,7 @@ void AliCentralitySelectionTask::UserExec(Option_t */*option*/)
     multV0C=esdV0->GetMTotV0C();
 
     float v0CorrR;
-    v0Corr = GetCorrV0(esd,v0CorrR,fRunNo);
-    v0Corr = (Short_t)v0Corr;
+    v0Corr = (Short_t)GetCorrV0(esd,v0CorrR);
     v0CorrResc = (Short_t)v0CorrR;
 
     // ***** Vertex Info
@@ -575,7 +574,7 @@ Int_t AliCentralitySelectionTask::SetupRun(AliESDEvent* esd)
   return 0;
 }
 //________________________________________________________________________
-Float_t AliCentralitySelectionTask::GetCorrV0(const AliESDEvent* esd, float &v0CorrResc, int run) 
+Float_t AliCentralitySelectionTask::GetCorrV0(const AliESDEvent* esd, float &v0CorrResc) 
 {
   // correct V0 non-linearity, prepare a version rescaled to SPD2 corr
   Double_t *par0;
@@ -647,7 +646,7 @@ Float_t AliCentralitySelectionTask::GetCorrV0(const AliESDEvent* esd, float &v0C
 			       3.92e-09 , 1.18e-09 , 2.26e-09 , 1.57e-09 , 2.02e-09 , 2.71e-09 , 2.99e-09 , 3.04e-09 }; 
   
   
-  if (run==137161) {
+  if (esd->GetRunNumber() <= 137165) {
     par0=par0_137161;
     par1=par1_137161;
     par2=par2_137161;
@@ -664,14 +663,7 @@ Float_t AliCentralitySelectionTask::GetCorrV0(const AliESDEvent* esd, float &v0C
   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;
-    if (s<0) {
-      printf("FPE %d %.2f %.2f %.2e\n",i,esdV0->GetMultiplicity(i),b,(b*b-4.*par2[i]*esdV0->GetMultiplicity(i)*esdV0->GetMultiplicity(i)));
-      n = -b;
-    }
-    else {
-      n = (-b + TMath::Sqrt(s));
-    }
+    Double_t n = (s<0) ? -b : (-b + TMath::Sqrt(s));
     multChCorr[i] = 2.*esdV0->GetMultiplicity(i)/n*par0[i];
     multCorr += multChCorr[i];
     multCorr2 += (multChCorr[i]/par0[i]/64.);
diff --git a/ANALYSIS/AliCentralitySelectionTask.h b/ANALYSIS/AliCentralitySelectionTask.h
index e62af0878a8..b63b84ab335 100644
--- a/ANALYSIS/AliCentralitySelectionTask.h
+++ b/ANALYSIS/AliCentralitySelectionTask.h
@@ -42,7 +42,7 @@ class AliCentralitySelectionTask : public AliAnalysisTaskSE {
   void ReadCentralityHistos(TString filename);
   void ReadCentralityHistos2(TString filename);
   
-  Float_t GetCorrV0(const AliESDEvent* esd, float &v0CorrResc, int run);
+  static Float_t GetCorrV0(const AliESDEvent* esd, float &v0CorrResc);
   Float_t GetCorrSPD2(Float_t spd2raw,Float_t zv);
 
  private:
-- 
2.39.3