Merge branch 'master' of https://ekryshen@git.cern.ch/reps/AliRoot

[u/mrichter/AliRoot.git] / PWGLF / totEt / AliAnalysisHadEt.cxx

diff --git a/PWGLF/totEt/AliAnalysisHadEt.cxx b/PWGLF/totEt/AliAnalysisHadEt.cxx
index e0cc8ef2ea0d2da48f537122d857f5c1cd6f71e6..89230e45922a28ff4215cc5c60db8d539a7db011 100644 (file)
--- a/PWGLF/totEt/AliAnalysisHadEt.cxx
+++ b/PWGLF/totEt/AliAnalysisHadEt.cxx
@@ -24,6 +24,12 @@
 #include "TParticle.h"
 #include "Rtypes.h"
 #include "AliPDG.h"
+#include "AliPIDResponse.h"
+#include "AliTPCPIDResponse.h" 
+#include "AliInputEventHandler.h"
+#include "AliAnalysisManager.h"
+#include "AliLog.h"
+#include "AliCentrality.h"
 
 using namespace std;
 
@@ -56,6 +62,7 @@ Int_t AliAnalysisHadEt::fgNumOfResBins = 80;
 
 
 AliAnalysisHadEt::AliAnalysisHadEt() : AliAnalysisEtCommon()
+                                    ,fPIDResponse(0)
                                     ,fSumEt(0)
                                     ,fSumEtAcc(0)
                                     ,fTotEt(0)
@@ -78,6 +85,7 @@ AliAnalysisHadEt::~AliAnalysisHadEt()
   delete fEsdtrackCutsITSTPC;
   delete fEsdtrackCutsITS;
   delete fEsdtrackCutsTPC;
+  delete fPIDResponse;
 }
 
 Int_t AliAnalysisHadEt::AnalyseEvent(AliVEvent *event)
@@ -91,6 +99,7 @@ Int_t AliAnalysisHadEt::AnalyseEvent(AliVEvent *event)
 void AliAnalysisHadEt::Init()
 {// clear variables, set up cuts and PDG info
   AliAnalysisEtCommon::Init();
+  
 }
 
 void AliAnalysisHadEt::ResetEventValues()
@@ -309,3 +318,22 @@ Float_t AliAnalysisHadEt::TrueP(float pTrec) const {
   if(pTrec>1.0) return pTrec;
   return pTrec/(1-599.334*pTrec+7285.15*pTrec*pTrec)+pTrec;
 }
+
+Int_t AliAnalysisHadEt::GetCentralityBin(Int_t numberofbins,AliCentrality *centrality){
+  Int_t centralitybin = -1;
+  if(numberofbins<21) centralitybin= centrality->GetCentralityClass10(fCentralityMethod);
+  else{
+    if(numberofbins<41) centralitybin= centrality->GetCentralityClass5(fCentralityMethod);
+    else{
+      Float_t centpercent = centrality->GetCentralityPercentile(fCentralityMethod);
+      centralitybin= centrality->GetCentralityClass5(fCentralityMethod);
+      Float_t centralitybinfloat = (centpercent/2.5);
+      if(centralitybin>=0){
+       centralitybin =(Int_t) (centralitybinfloat);
+       //cout<<" centbin "<<centralitybin<<" centrality "<<centpercent<<" "<<centpercent/2.5 <<endl;
+      }
+    }
+  }
+  //cout<<" centrality bin "<<centralitybin<<endl;
+  return centralitybin;
+}