Modified cut classes for Sigma* analysis to include cut on the TPC Nclusters. Config...
authorfbellini <fbellini@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 4 Mar 2013 13:33:21 +0000 (13:33 +0000)
committerfbellini <fbellini@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 4 Mar 2013 13:33:21 +0000 (13:33 +0000)
PWGLF/RESONANCES/AliRsnCutDaughterSigmaStar2010PP.cxx
PWGLF/RESONANCES/AliRsnCutDaughterSigmaStar2010PP.h
PWGLF/RESONANCES/AliRsnCutV0.cxx
PWGLF/RESONANCES/AliRsnCutV0.h
PWGLF/RESONANCES/macros/mini/AddAnalysisTaskSigmaStarPP.C
PWGLF/RESONANCES/macros/mini/AddAnalysisTaskSigmaStarPbPb.C
PWGLF/RESONANCES/macros/mini/ConfigSigmaStar.C
PWGLF/RESONANCES/macros/mini/ConfigSigmaStarPbPb.C

index 1860eb5..26d1741 100644 (file)
@@ -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);
index 12a1c37..b688cf8 100644 (file)
@@ -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
 };
 
 //__________________________________________________________________________________________________
index 9f5b276..efb6e2b 100644 (file)
@@ -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;
index 5ff476d..0506fb2 100644 (file)
@@ -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
index 28f1fe8..4513b7b 100644 (file)
@@ -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 --------------------------------------------------------------------------------
index cc75f91..33f42d4 100644 (file)
@@ -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 --------------------------------------------------------------------------------
index 96ee172..2524d2b 100644 (file)
@@ -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
index 0b8f428..b9f35ad 100644 (file)
@@ -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);