Swap D0/Dobar hypothesis for reflections (Carlos)

diff --git a/PWGHF/vertexingHF/charmFlow/AliAnalysisTaskFlowD2H.cxx b/PWGHF/vertexingHF/charmFlow/AliAnalysisTaskFlowD2H.cxx
index ab03bdf9fdee7742b864eb73bd66bd8212083f62..0bf84b266bb3a7c24583c2f24327ca54b5b6be61 100644 (file)
--- a/PWGHF/vertexingHF/charmFlow/AliAnalysisTaskFlowD2H.cxx
+++ b/PWGHF/vertexingHF/charmFlow/AliAnalysisTaskFlowD2H.cxx
@@ -63,7 +63,7 @@ ClassImp(AliAnalysisTaskFlowD2H)
 AliAnalysisTaskFlowD2H::AliAnalysisTaskFlowD2H() :
   AliAnalysisTaskSE(), fTPCEvent(NULL), fVZEEvent(NULL), 
   fCutsTPC(NULL), fCutsVZE(NULL), fNoPOIs(NULL), fCutsPOI(NULL),
-  fSource(0), fDebugV2(kFALSE), fMassBins(0), fMinMass(0.),
+  fSource(0), fDebugV2(kFALSE), fSwap(kFALSE), fMassBins(0), fMinMass(0.),
   fMaxMass(0.), fPtBinWidth(0), fHList(NULL), fEvent(NULL), 
   fCC(NULL), fRFPMTPC(NULL), fRFPPhiTPC(NULL), fCandidates(NULL)
 {
@@ -78,7 +78,7 @@ AliAnalysisTaskFlowD2H::AliAnalysisTaskFlowD2H(const char *name,
                                               Int_t specie) :
   AliAnalysisTaskSE(name), fTPCEvent(NULL), fVZEEvent(NULL), 
   fCutsTPC(cutsTPC), fCutsVZE(cutsVZE), fNoPOIs(NULL), fCutsPOI(cutsPOIs),
-  fSource(specie), fDebugV2(kFALSE), fMassBins(0), fMinMass(0.),
+  fSource(specie), fDebugV2(kFALSE), fSwap(kFALSE), fMassBins(0), fMinMass(0.),
   fMaxMass(0.), fPtBinWidth(0), fHList(NULL), fEvent(NULL),
   fCC(NULL), fRFPMTPC(NULL), fRFPPhiTPC(NULL), fCandidates(NULL)
 {
@@ -234,6 +234,7 @@ void AliAnalysisTaskFlowD2H::UserExec(Option_t *)
 void AliAnalysisTaskFlowD2H::FillD0toKpi(const AliAODEvent *theAOD)
 {
   // Fill D0->Kpi histos
+  fSwap = kTRUE;
   TList *listHF = (TList*) theAOD->GetList();
   if(!listHF) return;
   TClonesArray *listDzero = (TClonesArray*) listHF->FindObject("D0toKpi");
@@ -254,6 +255,10 @@ void AliAnalysisTaskFlowD2H::FillD0toKpi(const AliAODEvent *theAOD)
     Int_t topCut = fCutsD0toKpi->IsSelected( d0cand, AliRDHFCuts::kAll, NULL );
     if(!topCut) continue;
     Double_t massD0=topCut>1?d0cand->InvMassD0bar():d0cand->InvMassD0();
+    if((fSwap)&&(topCut==3)) {
+      massD0=d0cand->InvMassD0();
+    }
+
     // TO HANDLE AUTOCORRELATIONS
     nIDs[0] = ( (AliAODTrack*)d0cand->GetDaughter(0) )->GetID();
     nIDs[1] = ( (AliAODTrack*)d0cand->GetDaughter(1) )->GetID();

diff --git a/PWGHF/vertexingHF/charmFlow/AliAnalysisTaskFlowD2H.h b/PWGHF/vertexingHF/charmFlow/AliAnalysisTaskFlowD2H.h
index d9b77de82f6634f63ef6d7f4c5998d70ab4b85b2..53761a75dfddf4d7e1b8f6201cc247abda870583 100644 (file)
--- a/PWGHF/vertexingHF/charmFlow/AliAnalysisTaskFlowD2H.h
+++ b/PWGHF/vertexingHF/charmFlow/AliAnalysisTaskFlowD2H.h
@@ -35,6 +35,7 @@ class AliAnalysisTaskFlowD2H : public AliAnalysisTaskSE {
                            AliRDHFCuts *cutsPOIs, 
                            Int_t specie );
     void SetDebug() {fDebugV2 = true;}
+    void SetSwapAsumption() {fSwap = true;}
     virtual ~AliAnalysisTaskFlowD2H();
     virtual void UserCreateOutputObjects();
     virtual void UserExec(Option_t *);
@@ -65,6 +66,7 @@ class AliAnalysisTaskFlowD2H : public AliAnalysisTaskSE {
     AliRDHFCuts *fCutsPOI; // cuts for POIs
     Int_t  fSource; // AliRDHFCuts::ESele
     Bool_t fDebugV2; // fully talkative task
+    Bool_t fSwap;   // swap assumption (for neutral)
 
     Int_t fMassBins; // configures mass bins for the analysis
     Double_t fMinMass; // configures mass range for the analysis
@@ -79,7 +81,7 @@ class AliAnalysisTaskFlowD2H : public AliAnalysisTaskSE {
 
     TObjArray *fCandidates; // Array of selected candidates
 
-  ClassDef(AliAnalysisTaskFlowD2H, 3);
+  ClassDef(AliAnalysisTaskFlowD2H, 4);
 };
 
 #endif