Merge remote-tracking branch 'origin/master' into mergingFlat

[u/mrichter/AliRoot.git] / PWGLF / RESONANCES / AliRsnCutDaughterSigmaStar2010PP.cxx

diff --git a/PWGLF/RESONANCES/AliRsnCutDaughterSigmaStar2010PP.cxx b/PWGLF/RESONANCES/AliRsnCutDaughterSigmaStar2010PP.cxx
index 23a60a48910363ef349c212aeae18b51eeb3eda2..9e985df0b8d562e394f792b1dfa95550ae4506e4 100644 (file)
--- a/PWGLF/RESONANCES/AliRsnCutDaughterSigmaStar2010PP.cxx
+++ b/PWGLF/RESONANCES/AliRsnCutDaughterSigmaStar2010PP.cxx
@@ -19,7 +19,10 @@ AliRsnCutDaughterSigmaStar2010PP::AliRsnCutDaughterSigmaStar2010PP(const char *n
    AliRsnCut(name, AliRsnTarget::kDaughter),
    fPID(pid),
    fCutQuality(Form("%sQuality", name)),
-   fPIDCut(3.0)
+   fPIDCut(3.0),
+   fMinTPCcluster(70),
+   //fDCARptFormula("0.0182+0.0350/pt^1.01")
+   fDCARmax(0.05)
 {
 //
 // Constructor
@@ -28,13 +31,13 @@ AliRsnCutDaughterSigmaStar2010PP::AliRsnCutDaughterSigmaStar2010PP(const char *n
 
    fCutQuality.SetPtRange(0.15, 1E+20);
    fCutQuality.SetEtaRange(-0.8, 0.8);
-   fCutQuality.SetDCARmax(0.05);
-   //fCutQuality.SetDCARPtFormula("0.0182+0.0350/pt^1.01");
+   fCutQuality.SetDCARmax(fDCARmax);
+   //fCutQuality.SetDCARPtFormula(fDCARptFormula);
    fCutQuality.SetDCAZmax(2.0);
    fCutQuality.SetSPDminNClusters(1);
    fCutQuality.SetITSminNClusters(0);
    fCutQuality.SetITSmaxChi2(1E+20);
-   fCutQuality.SetTPCminNClusters(70);
+   fCutQuality.SetTPCminNClusters(fMinTPCcluster);
    fCutQuality.SetTPCmaxChi2(4.0);
    fCutQuality.SetRejectKinkDaughters();
    fCutQuality.SetAODTestFilterBit(5);
@@ -57,14 +60,22 @@ Bool_t AliRsnCutDaughterSigmaStar2010PP::IsSelected(TObject *obj)
       return kFALSE;
    }
 
+   AliDebugClass(2, "Checking status");
+
+
    // check flags
    if ((track->GetStatus() & AliESDtrack::kTPCin   ) == 0) return kFALSE;
    if ((track->GetStatus() & AliESDtrack::kTPCrefit) == 0) return kFALSE;
    if ((track->GetStatus() & AliESDtrack::kITSrefit) == 0) return kFALSE;
 
+   AliDebugClass(2, "Checking quality cuts");
+
    // quality
    if (!fCutQuality.IsSelected(obj)) return kFALSE;
 
+   AliDebugClass(2, "...passed");
+
+
    // check initialization of PID object
    AliPIDResponse *pid = fEvent->GetPIDResponse();
    if (!pid) {
@@ -94,22 +105,7 @@ Bool_t AliRsnCutDaughterSigmaStar2010PP::IsSelected(TObject *obj)
    //For the moment TOF is not used - PID ONLY WITH TPC - 3 sigmas in the whole range
 
    // TPC:
-   // below 350 MeV: 5sigma
-   // between 350 and 500 MeV: 3sigma
-   // pions above 500 MeV: 2sigma
-   // kaons between 500 and 700 MeV: 2sigma
-   // kaons above 700 MeV: rejected
-   /*if (pTPC <= 0.35)
-      maxTPC = 5.0;
-   else if (pTPC <= 0.5)
-      maxTPC = 3.0;
-   else if (pTPC > 0.5 && fPID == AliPID::kPion)
-      maxTPC = 2.0;
-   else if (pTPC > 0.5 && pTPC <= 0.7 && fPID == AliPID::kKaon)
-      maxTPC = 2.0;
-   else
-      return kFALSE;*/
-   //maxTPC = 3.0;
+   
    maxTPC = fPIDCut;
    return (nsTPC <= maxTPC);