Updated macros and DCA cut for Sigma* analysis (M.Venaruzzo)

diff --git a/PWGLF/RESONANCES/AliRsnCutDaughterSigmaStar2010PP.cxx b/PWGLF/RESONANCES/AliRsnCutDaughterSigmaStar2010PP.cxx
index b6c3c45f733a8a43ca2de1550176b6d9666d2c8a..9e985df0b8d562e394f792b1dfa95550ae4506e4 100644 (file)
--- a/PWGLF/RESONANCES/AliRsnCutDaughterSigmaStar2010PP.cxx
+++ b/PWGLF/RESONANCES/AliRsnCutDaughterSigmaStar2010PP.cxx
@@ -21,7 +21,8 @@ AliRsnCutDaughterSigmaStar2010PP::AliRsnCutDaughterSigmaStar2010PP(const char *n
    fCutQuality(Form("%sQuality", name)),
    fPIDCut(3.0),
    fMinTPCcluster(70),
-   fDCARptFormula("")
+   //fDCARptFormula("0.0182+0.0350/pt^1.01")
+   fDCARmax(0.05)
 {
 //
 // Constructor
@@ -30,7 +31,8 @@ AliRsnCutDaughterSigmaStar2010PP::AliRsnCutDaughterSigmaStar2010PP(const char *n
 
    fCutQuality.SetPtRange(0.15, 1E+20);
    fCutQuality.SetEtaRange(-0.8, 0.8);
-   fCutQuality.SetDCARPtFormula(fDCARptFormula);
+   fCutQuality.SetDCARmax(fDCARmax);
+   //fCutQuality.SetDCARPtFormula(fDCARptFormula);
    fCutQuality.SetDCAZmax(2.0);
    fCutQuality.SetSPDminNClusters(1);
    fCutQuality.SetITSminNClusters(0);

diff --git a/PWGLF/RESONANCES/AliRsnCutDaughterSigmaStar2010PP.h b/PWGLF/RESONANCES/AliRsnCutDaughterSigmaStar2010PP.h
index 9666e94b2398340a6a2606d7d87e4a60b82be9e9..8d86c5a4b00a1fbd57664f8d7a8532ac6a2e064d 100644 (file)
--- a/PWGLF/RESONANCES/AliRsnCutDaughterSigmaStar2010PP.h
+++ b/PWGLF/RESONANCES/AliRsnCutDaughterSigmaStar2010PP.h
@@ -27,7 +27,8 @@ public:
 
    void           SetPIDCut(Double_t value)                    {fPIDCut = value;}
    void           SetMinTPCcluster(Int_t value)                {fMinTPCcluster = value;}
-   void           SetDCARPtFormula(const char *formula)        {fDCARptFormula = formula;}
+   //void           SetDCARPtFormula(const char *formula)        {fDCARptFormula = formula;}
+   void           SetDCARmax(Double_t value)                   {fDCARmax = value;}
     
 private:
 
@@ -40,7 +41,8 @@ protected:
 
    Double_t         fPIDCut;                    // nsigmas for pions
    Int_t            fMinTPCcluster;             // min allowed TPC cluster
-   TString          fDCARptFormula;             // min DCAR pt dependent formula
+   Double_t         fDCARmax;                    // max DCA for tracks
+   //TString          fDCARptFormula;             // min DCAR pt dependent formula
 };
 
 //__________________________________________________________________________________________________

diff --git a/PWGLF/RESONANCES/AliRsnCutV0.cxx b/PWGLF/RESONANCES/AliRsnCutV0.cxx
index f652902a74053a9667516b607e1f8e81308ecda8..897adb44f4420ae00c0bf7888b41985fb548df32 100644 (file)
--- a/PWGLF/RESONANCES/AliRsnCutV0.cxx
+++ b/PWGLF/RESONANCES/AliRsnCutV0.cxx
@@ -40,7 +40,6 @@ AliRsnCutV0::AliRsnCutV0(const char *name, Int_t hypothesis, AliPID::EParticleTy
    fMaxDaughtersDCA(0.5),
    fMinTPCcluster(70),
    fMaxRapidity(0.8),
-   fDCARptFormula(""),
    fPID(pid),
    fPID2(pid2),
    fPIDCutProton(0),
@@ -68,7 +67,6 @@ AliRsnCutV0::AliRsnCutV0(const AliRsnCutV0 &copy) :
    fMaxDaughtersDCA(copy.fMaxDaughtersDCA),
    fMinTPCcluster(copy.fMinTPCcluster),
    fMaxRapidity(copy.fMaxRapidity),
-   fDCARptFormula(copy.fDCARptFormula),   
    fPID(copy.fPID),
    fPID2(copy.fPID2),
    fPIDCutProton(copy.fPIDCutProton),
@@ -84,8 +82,6 @@ AliRsnCutV0::AliRsnCutV0(const AliRsnCutV0 &copy) :
 //
    fCutQuality.SetPtRange(0.15, 1E+20);
    fCutQuality.SetEtaRange(-0.8, 0.8);
-   fCutQuality.SetDCARPtFormula(fDCARptFormula);
-   fCutQuality.SetDCAZmax(2.0);
    fCutQuality.SetSPDminNClusters(1);
    fCutQuality.SetITSminNClusters(0);
    fCutQuality.SetITSmaxChi2(1E+20);
@@ -113,7 +109,6 @@ AliRsnCutV0 &AliRsnCutV0::operator=(const AliRsnCutV0 &copy)
    fMaxDaughtersDCA = copy.fMaxDaughtersDCA;
    fMinTPCcluster = copy.fMinTPCcluster;
    fMaxRapidity = copy.fMaxRapidity;
-   fDCARptFormula = copy.fDCARptFormula;
    fCutQuality = copy.fCutQuality;
    fPID = copy.fPID;
    fPID2 = copy.fPID2;

diff --git a/PWGLF/RESONANCES/AliRsnCutV0.h b/PWGLF/RESONANCES/AliRsnCutV0.h
index b0ebbbfeb6ef7fb752f7c9a8dd206ee31c34555d..01f1283b49d03ba713f8e0a349fcfc900ef8b400 100644 (file)
--- a/PWGLF/RESONANCES/AliRsnCutV0.h
+++ b/PWGLF/RESONANCES/AliRsnCutV0.h
@@ -31,7 +31,6 @@ public:
    void           SetMaxDaughtersDCA(Double_t value)       {fMaxDaughtersDCA = value;}
    void           SetMinTPCcluster(Int_t value)            {fMinTPCcluster = value;}
    void           SetMaxRapidity(Double_t value)           {fMaxRapidity = value;}
-   void           SetDCARPtFormula(const char *formula)    {fDCARptFormula = formula;}
    
    void           SetPIDCutProton(Double_t value)          {fPIDCutProton = value;}
    void           SetPIDCutPion(Double_t value)            {fPIDCutPion = value;}
@@ -58,12 +57,11 @@ protected:
    Double_t         fMaxDaughtersDCA;  // max allowed DCA between the two daughers
    Int_t            fMinTPCcluster;    // min allowed TOC cluster
    Double_t         fMaxRapidity;      // max allowed V0 rapidity
-   TString          fDCARptFormula;    // min DCAR pt dependent formula
    
    AliPID::EParticleType fPID;         // PID for track
    AliPID::EParticleType fPID2;        // PID for track
 
-   Double_t         fPIDCutProton;          // nsigmas for protons
+   Double_t         fPIDCutProton;        // nsigmas for protons
    Double_t         fPIDCutPion;          // nsigmas for pions
    
    AliESDtrackCuts *fESDtrackCuts;     // quality cuts for v0 daughters

diff --git a/PWGLF/RESONANCES/macros/mini/AddAnalysisTaskSigmaStarPP.C b/PWGLF/RESONANCES/macros/mini/AddAnalysisTaskSigmaStarPP.C
index c7da47abce7f68d241de1a8d5acdc7bd909efcfa..781fff3b36d2905b2b0f4cfc182824c3d43696e8 100644 (file)
--- a/PWGLF/RESONANCES/macros/mini/AddAnalysisTaskSigmaStarPP.C
+++ b/PWGLF/RESONANCES/macros/mini/AddAnalysisTaskSigmaStarPP.C
@@ -49,7 +49,7 @@ AliRsnMiniAnalysisTask * AddAnalysisTaskSigmaStarPP
    } 
 
    // create the task and configure 
-   TString taskName = Form("SigmaStar%s%s_%.1f_%d_%.1f_%.1f_%.1f_%.2f_%.1f_%.2f_%.1f", (isPP? "pp" : "PbPb"), (isMC ? "MC" : "Data"),cutV,NTPCcluster,piPIDCut,pPIDCut,trackDCAcut,massTol,lambdaDCA,lambdaCosPoinAn,lambdaDaughDCA);
+   TString taskName = Form("SigmaStar%s%s_%.1f_%d_%.1f_%.1f_%.2f_%.2f_%.1f_%.2f_%.1f", (isPP? "pp" : "PbPb"), (isMC ? "MC" : "Data"),cutV,NTPCcluster,piPIDCut,pPIDCut,trackDCAcut,massTol,lambdaDCA,lambdaCosPoinAn,lambdaDaughDCA);
    AliRsnMiniAnalysisTask *task = new AliRsnMiniAnalysisTask(taskName.Data(), isMC);
    if (!isMC && !isPP){
      Printf(Form("========== SETTING USE CENTRALITY PATCH AOD049 : %s", (aodN==49)? "yes" : "no"));
@@ -137,7 +137,7 @@ AliRsnMiniAnalysisTask * AddAnalysisTaskSigmaStarPP
    //  outputFileName += ":Rsn";
    Printf("AddAnalysisTaskSigmaStarPP - Set OutputFileName : \n %s\n", outputFileName.Data() );
    
-   AliAnalysisDataContainer *output = mgr->CreateContainer(Form("RsnOut_%s_%.1f_%d_%.1f_%.1f_%.1f_%.2f_%.1f_%.2f_%.1f",outNameSuffix.Data(),cutV,NTPCcluster,piPIDCut,pPIDCut,trackDCAcut,massTol,lambdaDCA,lambdaCosPoinAn,lambdaDaughDCA), 
+   AliAnalysisDataContainer *output = mgr->CreateContainer(Form("RsnOut_%s_%.1f_%d_%.1f_%.1f_%.2f_%.2f_%.1f_%.2f_%.1f",outNameSuffix.Data(),cutV,NTPCcluster,piPIDCut,pPIDCut,trackDCAcut,massTol,lambdaDCA,lambdaCosPoinAn,lambdaDaughDCA), 
                                                           TList::Class(), 
                                                           AliAnalysisManager::kOutputContainer, 
                                                           outputFileName);

diff --git a/PWGLF/RESONANCES/macros/mini/AddAnalysisTaskSigmaStarPbPb.C b/PWGLF/RESONANCES/macros/mini/AddAnalysisTaskSigmaStarPbPb.C
index 5e63dfe655bc099f0f396db68b186be9d112c9c3..517c54d76a2a0bd57468d16a96ef3752ac7d9eb6 100644 (file)
--- a/PWGLF/RESONANCES/macros/mini/AddAnalysisTaskSigmaStarPbPb.C
+++ b/PWGLF/RESONANCES/macros/mini/AddAnalysisTaskSigmaStarPbPb.C
@@ -23,7 +23,7 @@ AliRsnMiniAnalysisTask * AddAnalysisTaskSigmaStarPbPb
    Int_t       aodFilterBit = 5,
    Int_t       piPIDCut = 3.0,
    Int_t       pPIDCut = 3.0,
-   Float_t     piDCAcut = 0.05,
+   Float_t     trackDCAcut = 7.0,
    Float_t     massTol = 0.01,
    Float_t     lambdaDCA = 0.3,
    Float_t     lambdaCosPoinAn = 0.99,
@@ -49,7 +49,7 @@ AliRsnMiniAnalysisTask * AddAnalysisTaskSigmaStarPbPb
    } 
    
    // create the task and configure 
-   TString taskName = Form("SigmaStar%s%s_%.1f_%d", (isPP? "pp" : "PbPb"), (isMC ? "MC" : "Data"),cutV,NTPCcluster);
+   TString taskName = Form("SigmaStar%s%s_%.1f_%d_%.1f_%.1f_%.2f_%.2f_%.1f_%.2f_%.1f", (isPP? "pp" : "PbPb"), (isMC ? "MC" : "Data"),cutV,NTPCcluster,piPIDCut,pPIDCut,trackDCAcut,massTol,lambdaDCA,lambdaCosPoinAn,lambdaDaughDCA);
    AliRsnMiniAnalysisTask *task = new AliRsnMiniAnalysisTask(taskName.Data(), isMC);
    if (!isMC && !isPP){
      Printf(Form("========== SETTING USE CENTRALITY PATCH AOD049 : %s", (aodN==49)? "yes" : "no"));
@@ -139,7 +139,7 @@ AliRsnMiniAnalysisTask * AddAnalysisTaskSigmaStarPbPb
    //  outputFileName += ":Rsn";
    Printf("AddAnalysisTaskSigmaStarPbPb - Set OutputFileName : \n %s\n", outputFileName.Data() );
    
-   AliAnalysisDataContainer *output = mgr->CreateContainer(Form("RsnOut_%s_%.1f_%d",outNameSuffix.Data(), cutV, NTPCcluster), 
+   AliAnalysisDataContainer *output = mgr->CreateContainer(Form("RsnOut_%s_%.1f_%d_%.1f_%.1f_%.2f_%.2f_%.1f_%.2f_%.1f",outNameSuffix.Data(),cutV,NTPCcluster,piPIDCut,pPIDCut,trackDCAcut,massTol,lambdaDCA,lambdaCosPoinAn,lambdaDaughDCA), 
                                                           TList::Class(), 
                                                           AliAnalysisManager::kOutputContainer, 
                                                           outputFileName);

diff --git a/PWGLF/RESONANCES/macros/mini/ConfigSigmaStar.C b/PWGLF/RESONANCES/macros/mini/ConfigSigmaStar.C
index 95b362ad5622f35fd0b28ac57810fd64219f2f3c..1b24810a821a59bf157b8e2b9aa2d34dc4e6ba12 100644 (file)
--- a/PWGLF/RESONANCES/macros/mini/ConfigSigmaStar.C
+++ b/PWGLF/RESONANCES/macros/mini/ConfigSigmaStar.C
@@ -33,17 +33,18 @@ Bool_t ConfigSigmaStar
    // -- Define track cuts -------------------------------------------------------------------------
    //
 
-   TString s = ""; s+=trackDCAcut; s+="*(0.0026+0.0050/pt^1.01)";
-
-   const char *formula = s;
+   //TString s = ""; s+=trackDCAcut; s+="*(0.0026+0.0050/pt^1.01)";
 
+   //const char *formula = s;
+   
    // 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->SetDCARPtFormula(formula);               
+   //cutQuality->SetDCARPtFormula(formula);             
+   cutQuality->SetDCARmax(trackDCAcut);                 
     
    // cut set
    AliRsnCutSet *cutSetPi = new AliRsnCutSet("setPionForSigmaStar", AliRsnTarget::kDaughter);
@@ -56,7 +57,6 @@ Bool_t ConfigSigmaStar
    AliESDtrackCuts *esdTrackCuts = new AliESDtrackCuts("qualityDaughterLambda");
    
    esdTrackCuts->SetAcceptKinkDaughters(0); // 0 = kFalse
-   esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD, AliESDtrackCuts::kAny);
    esdTrackCuts->SetMaxChi2PerClusterTPC(4);
    esdTrackCuts->SetMinNClustersTPC(NTPCcluster);
    esdTrackCuts->SetRequireTPCRefit();
@@ -70,7 +70,6 @@ Bool_t ConfigSigmaStar
    cutLambda->SetMaxDaughtersDCA(lambdaDaughDCA);
    cutLambda->SetMinTPCcluster(NTPCcluster);
    cutLambda->SetMaxRapidity(0.8);
-   cutLambda->SetDCARPtFormula(formula);
    cutLambda->SetAODTestFilterBit(aodFilterBit);
    cutLambda->SetPIDCutProton(pPIDCut);
    cutLambda->SetPIDCutPion(piPIDCut);
@@ -92,7 +91,6 @@ Bool_t ConfigSigmaStar
    cutAntiLambda->SetMaxDaughtersDCA(lambdaDaughDCA);
    cutAntiLambda->SetMinTPCcluster(NTPCcluster);
    cutAntiLambda->SetMaxRapidity(0.8);
-   cutAntiLambda->SetDCARPtFormula(formula);
    cutAntiLambda->SetAODTestFilterBit(aodFilterBit);
    cutAntiLambda->SetPIDCutProton(pPIDCut);
    cutAntiLambda->SetPIDCutPion(piPIDCut);

diff --git a/PWGLF/RESONANCES/macros/mini/ConfigSigmaStarPbPb.C b/PWGLF/RESONANCES/macros/mini/ConfigSigmaStarPbPb.C
index 531510027adf6b310eabfe30b7d21101ad75519e..9ab3bd43e6a0c86d541b8f6265018fef3ef8ad3c 100644 (file)
--- a/PWGLF/RESONANCES/macros/mini/ConfigSigmaStarPbPb.C
+++ b/PWGLF/RESONANCES/macros/mini/ConfigSigmaStarPbPb.C
@@ -33,9 +33,9 @@ Bool_t ConfigSigmaStarPbPb
    // -- Define track cuts -------------------------------------------------------------------------
    //
 
-   TString s = ""; s+=trackDCAcut; s+="*(0.0026+0.0050/pt^1.01)";
+   //TString s = ""; s+=trackDCAcut; s+="*(0.0026+0.0050/pt^1.01)";
 
-   const char *formula = s;
+   //const char *formula = s;
 
    // integrated pion cut
    AliRsnCutDaughterSigmaStar2010PP *cutPi = new AliRsnCutDaughterSigmaStar2010PP("cutPionForSigmaStar", AliPID::kPion);
@@ -43,7 +43,8 @@ Bool_t ConfigSigmaStarPbPb
    cutPi->SetMinTPCcluster(NTPCcluster);
    AliRsnCutTrackQuality *cutQuality = (AliRsnCutTrackQuality*) cutPi->CutQuality();
    cutQuality->SetAODTestFilterBit(aodFilterBit);
-   cutQuality->SetDCARPtFormula(formula);               
+   //cutQuality->SetDCARPtFormula(formula);
+   cutQuality->SetDCARmax(trackDCAcut);                 
     
    // cut set
    AliRsnCutSet *cutSetPi = new AliRsnCutSet("setPionForSigmaStar", AliRsnTarget::kDaughter);
@@ -56,7 +57,6 @@ Bool_t ConfigSigmaStarPbPb
    AliESDtrackCuts *esdTrackCuts = new AliESDtrackCuts("qualityDaughterLambda");
    
    esdTrackCuts->SetAcceptKinkDaughters(0); // 0 = kFalse
-   esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD, AliESDtrackCuts::kAny);
    esdTrackCuts->SetMaxChi2PerClusterTPC(4);
    esdTrackCuts->SetMinNClustersTPC(NTPCcluster);
    esdTrackCuts->SetRequireTPCRefit();
@@ -70,7 +70,6 @@ Bool_t ConfigSigmaStarPbPb
    cutLambda->SetMaxDaughtersDCA(lambdaDaughDCA);
    cutLambda->SetMinTPCcluster(NTPCcluster);
    cutLambda->SetMaxRapidity(0.8);
-   cutLambda->SetDCARPtFormula(formula);
    cutLambda->SetAODTestFilterBit(aodFilterBit);
    cutLambda->SetPIDCutProton(pPIDCut);
    cutLambda->SetPIDCutPion(piPIDCut);
@@ -92,7 +91,6 @@ Bool_t ConfigSigmaStarPbPb
    cutAntiLambda->SetMaxDaughtersDCA(lambdaDaughDCA);
    cutAntiLambda->SetMinTPCcluster(NTPCcluster);
    cutAntiLambda->SetMaxRapidity(0.8);
-   cutAntiLambda->SetDCARPtFormula(formula);
    cutAntiLambda->SetAODTestFilterBit(aodFilterBit);
    cutAntiLambda->SetPIDCutProton(pPIDCut);
    cutAntiLambda->SetPIDCutPion(piPIDCut);