Use kAny for SPD by default

diff --git a/PWG3/vertexingHF/AddTaskD0Mass.C b/PWG3/vertexingHF/AddTaskD0Mass.C
index a895de265bdf0120d947e16c92c3157b1b46b798..f62e8aab306959cbea73169f7a300cae860c5c1f 100644 (file)
--- a/PWG3/vertexingHF/AddTaskD0Mass.C
+++ b/PWG3/vertexingHF/AddTaskD0Mass.C
@@ -1,4 +1,4 @@
-AliAnalysisTaskSED0Mass *AddTaskD0Mass(Int_t flag=0/*0 = D0,1 = LS*/,Bool_t readMC=kTRUE,Bool_t cutOnDistr=kFALSE)
+AliAnalysisTaskSED0Mass *AddTaskD0Mass(Int_t flag=0/*0 = D0,1 = LS*/,Bool_t readMC=kFALSE,Bool_t cutOnDistr=kFALSE)
 {
   //
   // AddTask for the AliAnalysisTaskSE for D0 candidates
@@ -71,9 +71,9 @@ AliAnalysisTaskSED0Mass *AddTaskD0Mass(Int_t flag=0/*0 = D0,1 = LS*/,Bool_t read
   //default
   esdTrackCuts->SetRequireTPCRefit(kTRUE);
   esdTrackCuts->SetRequireITSRefit(kTRUE);
-  esdTrackCuts->SetMinNClustersITS(5); // default is 5
+  esdTrackCuts->SetMinNClustersITS(4); // default is 5
   esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,
-                                        AliESDtrackCuts::kBoth); //test d0 asimmetry
+                                        AliESDtrackCuts::kAny); 
  // default is kBoth, otherwise kAny
   esdTrackCuts->SetMinDCAToVertexXY(0.);
   esdTrackCuts->SetPtRange(0.3,1.e10);
@@ -203,7 +203,7 @@ AliAnalysisTaskSED0Mass *AddTaskD0Mass(Int_t flag=0/*0 = D0,1 = LS*/,Bool_t read
 
   //cout<<"\tbefore SetCuts : npt = "<<RDHFD0toKpi->GetNPtBins()<<endl;
   RDHFD0toKpi->SetCuts(nvars,nptbins,rdcutsvalmine);
-  RDHFD0toKpi->PrintAll();
+  //  RDHFD0toKpi->PrintAll();
 
   // Aanalysis task    
   TString taskname="MassAndDistrAnalysis";

diff --git a/PWG3/vertexingHF/AliAnalysisTaskSED0Mass.cxx b/PWG3/vertexingHF/AliAnalysisTaskSED0Mass.cxx
index d43e8b324587dc1b359404fbd11764a2a68b786b..656827b3411c582cfd0bfa85f09ead8173728681 100644 (file)
--- a/PWG3/vertexingHF/AliAnalysisTaskSED0Mass.cxx
+++ b/PWG3/vertexingHF/AliAnalysisTaskSED0Mass.cxx
@@ -560,6 +560,11 @@ void AliAnalysisTaskSED0Mass::UserExec(Option_t */*option*/)
     Double_t eta1=d->EtaProng(1);
    
     if (TMath::Abs(eta0) < 0.9 && TMath::Abs(eta1) < 0.9) {
+      //apply cuts on tracks
+      Int_t isSelected = fCuts->IsSelected(d,AliRDHFCuts::kTracks);
+      if (!isSelected) return;
+      if(fDebug>1) cout<<"tracks selected";
+
       FillVarHists(d,mcArray,fCuts,fDistr);
       FillMassHists(d,mcArray,fCuts,fOutputMass);
     }
@@ -602,8 +607,9 @@ void AliAnalysisTaskSED0Mass::FillVarHists(AliAODRecoDecayHF2Prong *part, TClone
   //Double_t pt = d->Pt(); //mother pt
   Bool_t isSelected=kFALSE;
 
+
   if(fCutOnDistr){
-    isSelected = cuts->IsSelected(part,AliRDHFCuts::kAll);
+    isSelected = cuts->IsSelected(part,AliRDHFCuts::kCandidate);
     if (!isSelected){
       //cout<<"Not Selected"<<endl;
       return;
@@ -794,7 +800,7 @@ void AliAnalysisTaskSED0Mass::FillMassHists(AliAODRecoDecayHF2Prong *part, TClon
 
   Double_t mPDG=TDatabasePDG::Instance()->GetParticle(421)->Mass();
 
-  Int_t isSelected=cuts->IsSelected(part,AliRDHFCuts::kAll); //selected
+  Int_t isSelected=cuts->IsSelected(part,AliRDHFCuts::kCandidate); //selected
   //cout<<"is selected = "<<isSelected<<endl;
 
   //cout<<"check cuts = "<<endl;
@@ -999,4 +1005,3 @@ void AliAnalysisTaskSED0Mass::Terminate(Option_t */*option*/)
   return;
 }
 
-