updated

diff --git a/PWGHF/hfe/AliAnalysisTaskHFECal.cxx b/PWGHF/hfe/AliAnalysisTaskHFECal.cxx
index ea4fae2e29ead12dcbbb01fd6d218cb49eb58e86..82b79c90b47a2683cbcee1d96d1d7881cdb328bf 100644 (file)
--- a/PWGHF/hfe/AliAnalysisTaskHFECal.cxx
+++ b/PWGHF/hfe/AliAnalysisTaskHFECal.cxx
@@ -648,7 +648,13 @@ void AliAnalysisTaskHFECal::UserExec(Option_t*)
     eta = track->Eta();
     dEdx = track->GetTPCsignal();
     fTPCnSigma = fPID->GetPIDResponse() ? fPID->GetPIDResponse()->NumberOfSigmasTPC(track, AliPID::kElectron) : 1000;
-    
+    if(mcLabel==-1) // nsigma mean correction
+      {
+       double mean_corr = NsigCorr(cent);
+       //printf("correction %f\n",fTPCnSigma);
+       fTPCnSigma -= mean_corr;
+      }    
+
         double ncells = -1.0;
         double m20 = -1.0;
         double m02 = -1.0;
@@ -1701,7 +1707,16 @@ void AliAnalysisTaskHFECal::FindTriggerClusters()
 }
 
 
-
+double NsigCorr(int cent)
+{
+ double shift = 0.0;
+ if(cent>=20 && cent<30)shift = 0.156;
+ if(cent>=30 && cent<40)shift = 0.316;
+ if(cent>=40 && cent<50)shift = 0.336;
+ if(cent>=50 && cent<70)shift = 0.440;
+ if(cent>=70 && cent<90)shift = 0.534;
+ return shift;
+}
 
 
 

diff --git a/PWGHF/hfe/AliAnalysisTaskHFECal.h b/PWGHF/hfe/AliAnalysisTaskHFECal.h
index 977b1b7ccfcaab55a6ec0c333ef3d9ef18d69b10..74e97f00b1a5c393a1325b0fd105892b110fa5f8 100644 (file)
--- a/PWGHF/hfe/AliAnalysisTaskHFECal.h
+++ b/PWGHF/hfe/AliAnalysisTaskHFECal.h
@@ -59,6 +59,7 @@ class AliAnalysisTaskHFECal : public AliAnalysisTaskSE {
   double GetMCweight(double mcPi0pT);
   double GetMCweightEta(double mcEtapT);
   void FindTriggerClusters();
+  double NsigCorr(int cent);
  private:
   
   Bool_t ProcessCutStep(Int_t cutStep, AliVParticle *track);