Update

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

diff --git a/PWGHF/hfe/AliAnalysisTaskHFEFlow.h b/PWGHF/hfe/AliAnalysisTaskHFEFlow.h
index b9b1b6a714c53fd69034aed5148e1f2d6311624d..61f67f441275066897e102651a43463c11ee875e 100644 (file)
--- a/PWGHF/hfe/AliAnalysisTaskHFEFlow.h
+++ b/PWGHF/hfe/AliAnalysisTaskHFEFlow.h
@@ -126,6 +126,7 @@ public:
   void  SetMaxopening3D(Double_t maxOpening3D) { fMaxopening3D = maxOpening3D; };\r
   void  SetMaxopeningtheta(Double_t maxOpeningtheta) { fMaxopeningtheta = maxOpeningtheta; };\r
   void  SetMaxopeningphi(Double_t maxOpeningphi) { fMaxopeningphi = maxOpeningphi; };\r
+  void  SetAlgorithmMA(Bool_t algorithmMA) { fAlgorithmMA = algorithmMA; };\r
 \r
   Int_t    LookAtNonHFE(Int_t iTrack1, AliVTrack *track1, AliVEvent *fESD, AliMCEvent *mcEvent,Int_t binct,Double_t deltaphi,Int_t source,Int_t indexmother);\r
   \r