Modified cut classes for Sigma* analysis to include cut on the TPC Nclusters. Config...

diff --git a/PWGLF/RESONANCES/AliRsnCutDaughterSigmaStar2010PP.cxx b/PWGLF/RESONANCES/AliRsnCutDaughterSigmaStar2010PP.cxx
index 1860eb5c4315855c5ff9612b44d95104aeb6b347..26d1741ab28f1b6970781b3988d66835ba998c9d 100644 (file)
--- a/PWGLF/RESONANCES/AliRsnCutDaughterSigmaStar2010PP.cxx
+++ b/PWGLF/RESONANCES/AliRsnCutDaughterSigmaStar2010PP.cxx
@@ -19,7 +19,8 @@ AliRsnCutDaughterSigmaStar2010PP::AliRsnCutDaughterSigmaStar2010PP(const char *n
    AliRsnCut(name, AliRsnTarget::kDaughter),
    fPID(pid),
    fCutQuality(Form("%sQuality", name)),
-   fPIDCut(3.0)
+   fPIDCut(3.0),
+   fMinTPCcluster(70)
 {
 //
 // Constructor
@@ -34,7 +35,7 @@ AliRsnCutDaughterSigmaStar2010PP::AliRsnCutDaughterSigmaStar2010PP(const char *n
    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);

diff --git a/PWGLF/RESONANCES/AliRsnCutDaughterSigmaStar2010PP.h b/PWGLF/RESONANCES/AliRsnCutDaughterSigmaStar2010PP.h
index 12a1c379823feac06516167954b216353a51eb10..b688cf8b344c1a66fb62bb503ca583858988dc20 100644 (file)
--- a/PWGLF/RESONANCES/AliRsnCutDaughterSigmaStar2010PP.h
+++ b/PWGLF/RESONANCES/AliRsnCutDaughterSigmaStar2010PP.h
@@ -25,7 +25,8 @@ public:
    virtual Bool_t         IsSelected(TObject *obj);
 
    void           SetPIDCut(Double_t value)                  {fPIDCut = value;}
-
+   void           SetMinTPCcluster(Int_t value)              {fMinTPCcluster = value;}
+   
 private:
 
    AliPID::EParticleType fPID;              // PID for track
@@ -36,7 +37,7 @@ private:
 protected:
 
    Double_t         fPIDCut;          // nsigmas for pions
-
+   Int_t            fMinTPCcluster;   // min allowed TPC cluster
 };
 
 //__________________________________________________________________________________________________

diff --git a/PWGLF/RESONANCES/AliRsnCutV0.cxx b/PWGLF/RESONANCES/AliRsnCutV0.cxx
index 9f5b276f3d94c5f090fb3c17fbdd89209d31df7b..efb6e2b355db894965497edf8784a5d2fb737cc5 100644 (file)
--- a/PWGLF/RESONANCES/AliRsnCutV0.cxx
+++ b/PWGLF/RESONANCES/AliRsnCutV0.cxx
@@ -39,6 +39,7 @@ AliRsnCutV0::AliRsnCutV0(const char *name, Int_t hypothesis, AliPID::EParticleTy
    fMaxDCAVertex(0.3),
    fMinCosPointAngle(0.95),
    fMaxDaughtersDCA(0.5),
+   fMinTPCcluster(70),
    fMaxRapidity(0.8),
    fPID(pid),
    fPID2(pid2),
@@ -66,6 +67,7 @@ AliRsnCutV0::AliRsnCutV0(const AliRsnCutV0 &copy) :
    fMaxDCAVertex(copy.fMaxDCAVertex),
    fMinCosPointAngle(copy.fMinCosPointAngle),
    fMaxDaughtersDCA(copy.fMaxDaughtersDCA),
+   fMinTPCcluster(copy.fMinTPCcluster),
    fMaxRapidity(copy.fMaxRapidity),
    fPID(copy.fPID),
    fPID2(copy.fPID2),
@@ -89,7 +91,7 @@ AliRsnCutV0::AliRsnCutV0(const AliRsnCutV0 &copy) :
    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);
@@ -111,6 +113,7 @@ AliRsnCutV0 &AliRsnCutV0::operator=(const AliRsnCutV0 &copy)
    fMaxDCAVertex = copy.fMaxDCAVertex;
    fMinCosPointAngle = copy.fMinCosPointAngle;
    fMaxDaughtersDCA = copy.fMaxDaughtersDCA;
+   fMinTPCcluster = copy.fMinTPCcluster;
    fMaxRapidity = copy.fMaxRapidity;
    fPID = copy.fPID;
    fPID2 = copy.fPID2;

diff --git a/PWGLF/RESONANCES/AliRsnCutV0.h b/PWGLF/RESONANCES/AliRsnCutV0.h
index 5ff476dfdf2536cb817845ea9086606de174840f..0506fb2a129a484fd0a6356b0a9531e0f498bb07 100644 (file)
--- a/PWGLF/RESONANCES/AliRsnCutV0.h
+++ b/PWGLF/RESONANCES/AliRsnCutV0.h
@@ -30,6 +30,7 @@ public:
    void           SetMaxDCAVertex(Double_t value)          {fMaxDCAVertex = value;}
    void           SetMinCosPointingAngle(Double_t value)   {fMinCosPointAngle = value;}
    void           SetMaxDaughtersDCA(Double_t value)       {fMaxDaughtersDCA = value;}
+   void           SetMinTPCcluster(Int_t value)            {fMinTPCcluster = value;}
    void           SetMaxRapidity(Double_t value)           {fMaxRapidity = value;}
 
    void           SetPIDCut1(Double_t value)               {fPIDCut1 = value;}
@@ -54,6 +55,7 @@ protected:
    Double_t         fMaxDCAVertex;     // max allowed DCA from primary vertex
    Double_t         fMinCosPointAngle; // min allowed cosine of pointing angle
    Double_t         fMaxDaughtersDCA;  // max allowed DCA between the two daughers
+   Int_t            fMinTPCcluster;    // min allowed TOC cluster
    Double_t         fMaxRapidity;      // max allowed V0 rapidity
 
    AliPID::EParticleType fPID;         // PID for track

diff --git a/PWGLF/RESONANCES/macros/mini/AddAnalysisTaskSigmaStarPP.C b/PWGLF/RESONANCES/macros/mini/AddAnalysisTaskSigmaStarPP.C
index 28f1fe88ee93242110fc8e56499278a09f18e4cb..4513b7be1c5db02cc57b9ab5251ad0977c800ad0 100644 (file)
--- a/PWGLF/RESONANCES/macros/mini/AddAnalysisTaskSigmaStarPP.C
+++ b/PWGLF/RESONANCES/macros/mini/AddAnalysisTaskSigmaStarPP.C
@@ -27,6 +27,7 @@ AliRsnMiniAnalysisTask * AddAnalysisTaskSigmaStarPP
    Float_t     lambdaDCA = 0.03,
    Float_t     lambdaCosPoinAn = 0.99,
    Float_t     lambdaDaughDCA = 0.5,
+   Int_t       NTPCcluster = 70,
    Int_t       nmix = 5,
    Float_t     maxDiffVzMix = 1.0,
    Float_t     maxDiffMultMix = 10.0,
@@ -126,7 +127,7 @@ AliRsnMiniAnalysisTask * AddAnalysisTaskSigmaStarPP
        Printf("========================== MC analysis - PID cuts not used");
    } else 
      Printf("========================== DATA analysis - PID cuts used");
-   if (!ConfigSigmaStar(task, isPP, isMC, piPIDCut, pPIDCut, aodFilterBit, piDCAcut, massTol, lambdaDCA, lambdaCosPoinAn, lambdaDaughDCA "", cutsPair)) return 0x0;
+   if (!ConfigSigmaStar(task, isPP, isMC, piPIDCut, pPIDCut, aodFilterBit, piDCAcut, massTol, lambdaDCA, lambdaCosPoinAn, lambdaDaughDCA, NTPCcluster, "", cutsPair)) return 0x0;
    
    //
    // -- CONTAINERS --------------------------------------------------------------------------------

diff --git a/PWGLF/RESONANCES/macros/mini/AddAnalysisTaskSigmaStarPbPb.C b/PWGLF/RESONANCES/macros/mini/AddAnalysisTaskSigmaStarPbPb.C
index cc75f913ffcae0a4fdee47bdb572a277f326df48..33f42d4355146352cd63d9ad119156c0b70e6490 100644 (file)
--- a/PWGLF/RESONANCES/macros/mini/AddAnalysisTaskSigmaStarPbPb.C
+++ b/PWGLF/RESONANCES/macros/mini/AddAnalysisTaskSigmaStarPbPb.C
@@ -27,6 +27,7 @@ AliRsnMiniAnalysisTask * AddAnalysisTaskSigmaStarPbPb
    Float_t     lambdaDCA = 0.03,
    Float_t     lambdaCosPoinAn = 0.99,
    Float_t     lambdaDaughDCA = 0.5,
+   Int_t       NTPCcluster = 70,
    Int_t       nmix = 5,
    Float_t     maxDiffVzMix = 1.0,
    Float_t     maxDiffMultMix = 10.0,
@@ -128,7 +129,7 @@ AliRsnMiniAnalysisTask * AddAnalysisTaskSigmaStarPbPb
    } else 
      Printf("========================== DATA analysis - PID cuts used");
      
-   if (!ConfigSigmaStarPbPb(task, isPP, isMC, piPIDCut, pPIDCut, aodFilterBit, piDCAcut, massTol, lambdaDCA, lambdaCosPoinAn, lambdaDaughDCA, "", cutsPair)) return 0x0;
+   if (!ConfigSigmaStarPbPb(task, isPP, isMC, piPIDCut, pPIDCut, aodFilterBit, piDCAcut, massTol, lambdaDCA, lambdaCosPoinAn, lambdaDaughDCA, NTPCcluster, "", cutsPair)) return 0x0;
    
    //
    // -- CONTAINERS --------------------------------------------------------------------------------

diff --git a/PWGLF/RESONANCES/macros/mini/ConfigSigmaStar.C b/PWGLF/RESONANCES/macros/mini/ConfigSigmaStar.C
index 96ee17232305e1fd3185f11cca30ea1bd9411bf2..2524d2bc4c799c08a194dd59a535206a34525176 100644 (file)
--- a/PWGLF/RESONANCES/macros/mini/ConfigSigmaStar.C
+++ b/PWGLF/RESONANCES/macros/mini/ConfigSigmaStar.C
@@ -21,6 +21,7 @@ Bool_t ConfigSigmaStar
    Float_t                 lambdaDCA,
    Float_t                 lambdaCosPoinAn,
    Float_t                 lambdaDaughDCA,
+   Int_t                   NTPCcluster,
    const char             *suffix,
    AliRsnCutSet           *cutsPair
 )
@@ -52,7 +53,7 @@ Bool_t ConfigSigmaStar
    esdTrackCuts->SetAcceptKinkDaughters(0); // 0 = kFalse
    esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD, AliESDtrackCuts::kAny);
    esdTrackCuts->SetMaxChi2PerClusterTPC(4);
-   esdTrackCuts->SetMinNClustersTPC(70);
+   esdTrackCuts->SetMinNClustersTPC(NTPCcluster);
    esdTrackCuts->SetRequireTPCRefit();
    
    // cut lambda

diff --git a/PWGLF/RESONANCES/macros/mini/ConfigSigmaStarPbPb.C b/PWGLF/RESONANCES/macros/mini/ConfigSigmaStarPbPb.C
index 0b8f42823c292be6b4e5a242a06767c9ea4926a8..b9f35adc765ded5e4fc83f7daae4dd9d1b30bf6f 100644 (file)
--- a/PWGLF/RESONANCES/macros/mini/ConfigSigmaStarPbPb.C
+++ b/PWGLF/RESONANCES/macros/mini/ConfigSigmaStarPbPb.C
@@ -21,6 +21,7 @@ Bool_t ConfigSigmaStarPbPb
    Float_t                 lambdaDCA,
    Float_t                 lambdaCosPoinAn,
    Float_t                 lambdaDaughDCA,
+   Int_t                   NTPCcluster,
    const char             *suffix,
    AliRsnCutSet           *cutsPair
 )
@@ -35,6 +36,7 @@ Bool_t ConfigSigmaStarPbPb
    // integrated pion cut
    AliRsnCutDaughterSigmaStar2010PP *cutPi = new AliRsnCutDaughterSigmaStar2010PP("cutPionForSigmaStar", AliPID::kPion);
    cutPi->SetPIDCut(piPIDCut);
+   cutPi->SetMinTPCcluster(NTPCcluster);
    AliRsnCutTrackQuality *cutQuality = (AliRsnCutTrackQuality*) cutPi->CutQuality();
    cutQuality->SetAODTestFilterBit(aodFilterBit);
    cutQuality->SetDCARmax(piDCAcut);            
@@ -52,7 +54,7 @@ Bool_t ConfigSigmaStarPbPb
    esdTrackCuts->SetAcceptKinkDaughters(0); // 0 = kFalse
    esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD, AliESDtrackCuts::kAny);
    esdTrackCuts->SetMaxChi2PerClusterTPC(4);
-   esdTrackCuts->SetMinNClustersTPC(70);
+   esdTrackCuts->SetMinNClustersTPC(NTPCcluster);
    esdTrackCuts->SetRequireTPCRefit();
    
    // cut lambda
@@ -62,6 +64,7 @@ Bool_t ConfigSigmaStarPbPb
    cutLambda->SetMaxDCAVertex(lambdaDCA);
    cutLambda->SetMinCosPointingAngle(lambdaCosPoinAn);
    cutLambda->SetMaxDaughtersDCA(lambdaDaughDCA);
+   cutLambda->SetMinTPCcluster(NTPCcluster);
    cutLambda->SetMaxRapidity(0.8);
    cutLambda->SetAODTestFilterBit(aodFilterBit);
    cutLambda->SetPIDCut1(pPIDCut);
@@ -83,6 +86,7 @@ Bool_t ConfigSigmaStarPbPb
    cutAntiLambda->SetMaxDCAVertex(lambdaDCA);
    cutAntiLambda->SetMinCosPointingAngle(lambdaCosPoinAn);
    cutAntiLambda->SetMaxDaughtersDCA(lambdaDaughDCA);
+   cutAntiLambda->SetMinTPCcluster(NTPCcluster);
    cutAntiLambda->SetMaxRapidity(0.8);
    cutAntiLambda->SetAODTestFilterBit(aodFilterBit);
    cutAntiLambda->SetPIDCut1(pPIDCut);