Update

[u/mrichter/AliRoot.git] / PWGHF / hfe / AliAnalysisTaskHFEFlow.cxx

diff --git a/PWGHF/hfe/AliAnalysisTaskHFEFlow.cxx b/PWGHF/hfe/AliAnalysisTaskHFEFlow.cxx
index a64a5fbd0bb6b5e03f291832bb2d24fcb2ca5e9f..00648903a65f99c3fe6eb01f826c37ad52ee7410 100644 (file)
--- a/PWGHF/hfe/AliAnalysisTaskHFEFlow.cxx
+++ b/PWGHF/hfe/AliAnalysisTaskHFEFlow.cxx
@@ -570,9 +570,9 @@ void AliAnalysisTaskHFEFlow::UserCreateOutputObjects()
   Double_t binLimEta[nBinsEta+1];\r
   for(Int_t i=0; i<=nBinsEta; i++) binLimEta[i]=(Double_t)minEta + (maxEta-minEta)/nBinsEta*(Double_t)i ;\r
 \r
-  Int_t nBinsStep = 6;\r
+  Int_t nBinsStep = 7;\r
   Double_t minStep = 0.;\r
-  Double_t maxStep = 6.;\r
+  Double_t maxStep = 7.;\r
   Double_t binLimStep[nBinsStep+1];\r
   for(Int_t i=0; i<=nBinsStep; i++) binLimStep[i]=(Double_t)minStep + (maxStep-minStep)/nBinsStep*(Double_t)i ;\r
 \r
@@ -1134,7 +1134,9 @@ void AliAnalysisTaskHFEFlow::UserExec(Option_t */*option*/)
   else {\r
     \r
     eventPlaneV0 = TVector2::Phi_0_2pi(vEPa->GetEventplane("V0", fInputEvent,2));\r
+    //printf("eventPlaneV0 %f\n",eventPlaneV0);\r
     if(eventPlaneV0 > TMath::Pi()) eventPlaneV0 = eventPlaneV0 - TMath::Pi();\r
+    //printf("eventPlaneV0 %f\n",eventPlaneV0);\r
     eventPlaneV0A = TVector2::Phi_0_2pi(vEPa->GetEventplane("V0A", fInputEvent,2));\r
     if(eventPlaneV0A > TMath::Pi()) eventPlaneV0A = eventPlaneV0A - TMath::Pi();\r
     eventPlaneV0C = TVector2::Phi_0_2pi(vEPa->GetEventplane("V0C", fInputEvent,2));\r
@@ -1433,31 +1435,31 @@ void AliAnalysisTaskHFEFlow::UserExec(Option_t */*option*/)
        } // Quick and dirty fix to reject both kink mothers and daughters\r
       }\r
             \r
-      valuetrackingcuts[1] = 0; \r
+      valuetrackingcuts[1] = 1; \r
        if(fDebugLevel > 3) fTrackingCuts->Fill(&valuetrackingcuts[0]);    \r
       // RecPrim\r
       if(!fHFECuts->CheckParticleCuts(AliHFEcuts::kStepRecPrim + AliHFEcuts::kNcutStepsMCTrack, (TObject *)track)) continue;\r
-      valuetrackingcuts[1] = 1; \r
+      valuetrackingcuts[1] = 2; \r
        if(fDebugLevel > 3) fTrackingCuts->Fill(&valuetrackingcuts[0]);    \r
 \r
       // HFEcuts: ITS layers cuts\r
        if(!fHFECuts->CheckParticleCuts(AliHFEcuts::kStepHFEcutsITS + AliHFEcuts::kNcutStepsMCTrack, (TObject *)track)) continue;\r
-       valuetrackingcuts[1] = 2; \r
+       valuetrackingcuts[1] = 3; \r
        if(fDebugLevel > 3) fTrackingCuts->Fill(&valuetrackingcuts[0]);     \r
 \r
       // HFE cuts: TOF PID and mismatch flag\r
       if(!fHFECuts->CheckParticleCuts(AliHFEcuts::kStepHFEcutsTOF + AliHFEcuts::kNcutStepsMCTrack, (TObject *)track)) continue;\r
-      valuetrackingcuts[1] = 3; \r
+      valuetrackingcuts[1] = 4; \r
        if(fDebugLevel > 3) fTrackingCuts->Fill(&valuetrackingcuts[0]);    \r
       \r
       // HFE cuts: TPC PID cleanup\r
       if(!fHFECuts->CheckParticleCuts(AliHFEcuts::kStepHFEcutsTPC + AliHFEcuts::kNcutStepsMCTrack, (TObject *)track)) continue;\r
-      valuetrackingcuts[1] = 4; \r
+      valuetrackingcuts[1] = 5; \r
        if(fDebugLevel > 3) fTrackingCuts->Fill(&valuetrackingcuts[0]);    \r
       \r
       // HFEcuts: Nb of tracklets TRD0\r
       if(!fHFECuts->CheckParticleCuts(AliHFEcuts::kStepHFEcutsTRD + AliHFEcuts::kNcutStepsMCTrack, (TObject *)track)) continue;\r
-      valuetrackingcuts[1] = 5; \r
+      valuetrackingcuts[1] = 6; \r
        if(fDebugLevel > 3) fTrackingCuts->Fill(&valuetrackingcuts[0]);    \r
       \r
     }\r