From 94e5e7bcae101c71e433b46129f1546766d3ed59 Mon Sep 17 00:00:00 2001 From: Prottay Das Date: Thu, 14 Aug 2025 19:17:13 +0200 Subject: [PATCH 1/4] removed the cuts on findable clusters --- PWGLF/Tasks/Strangeness/lambdapolsp.cxx | 24 +++++++++++++----------- 1 file changed, 13 insertions(+), 11 deletions(-) diff --git a/PWGLF/Tasks/Strangeness/lambdapolsp.cxx b/PWGLF/Tasks/Strangeness/lambdapolsp.cxx index 492f71e8d7d..b13147da53f 100644 --- a/PWGLF/Tasks/Strangeness/lambdapolsp.cxx +++ b/PWGLF/Tasks/Strangeness/lambdapolsp.cxx @@ -142,8 +142,8 @@ struct lambdapolsp { Configurable ConfDaughEta{"ConfDaughEta", 0.8f, "V0 Daugh sel: max eta"}; Configurable cfgDaughPrPt{"cfgDaughPrPt", 0.4, "minimum daughter proton pt"}; Configurable cfgDaughPiPt{"cfgDaughPiPt", 0.2, "minimum daughter pion pt"}; - Configurable rcrfc{"rcrfc", 0.8f, "Ratio of CR to FC"}; - Configurable ConfDaughTPCnclsMin{"ConfDaughTPCnclsMin", 50.f, "V0 Daugh sel: Min. nCls TPC"}; + // Configurable rcrfc{"rcrfc", 0.8f, "Ratio of CR to FC"}; + // Configurable ConfDaughTPCnclsMin{"ConfDaughTPCnclsMin", 50.f, "V0 Daugh sel: Min. nCls TPC"}; Configurable ConfDaughPIDCuts{"ConfDaughPIDCuts", 3, "PID selections for Lambda daughters"}; Configurable usesubdet{"usesubdet", false, "use subdet"}; Configurable useAccCorr{"useAccCorr", false, "use acceptance correction"}; @@ -483,12 +483,13 @@ struct lambdapolsp { /*if (TMath::Abs(eta) > ConfDaughEta) { return false; }*/ + /* if (tpcNClsF < ConfDaughTPCnclsMin) { return false; } if (track.tpcCrossedRowsOverFindableCls() < rcrfc) { return false; - } + }*/ if (pid == 0 && TMath::Abs(track.tpcNSigmaPr()) > ConfDaughPIDCuts) { return false; @@ -541,13 +542,14 @@ struct lambdapolsp { if (posTrackExtra.tpcNClsCrossedRows() < cfgTPCcluster || negTrackExtra.tpcNClsCrossedRows() < cfgTPCcluster) { return false; } + /* if (posTrackExtra.tpcNClsFound() < ConfDaughTPCnclsMin || negTrackExtra.tpcNClsFound() < ConfDaughTPCnclsMin) { return false; } if (posTrackExtra.tpcCrossedRowsOverFindableCls() < rcrfc || negTrackExtra.tpcCrossedRowsOverFindableCls() < rcrfc) { return false; } - + */ // check TPC PID if (pid == 0 && ((std::abs(posTrackExtra.tpcNSigmaPr()) > ConfDaughPIDCuts) || (std::abs(negTrackExtra.tpcNSigmaPi()) > ConfDaughPIDCuts))) { return false; @@ -586,13 +588,13 @@ struct lambdapolsp { // check TPC tracking properties if (posTrackExtra.tpcNClsCrossedRows() < cfgTPCcluster || negTrackExtra.tpcNClsCrossedRows() < cfgTPCcluster) { return false; - } - if (posTrackExtra.tpcNClsFound() < ConfDaughTPCnclsMin || negTrackExtra.tpcNClsFound() < ConfDaughTPCnclsMin) { - return false; - } - if (posTrackExtra.tpcCrossedRowsOverFindableCls() < rcrfc || negTrackExtra.tpcCrossedRowsOverFindableCls() < rcrfc) { - return false; - } + } /* + if (posTrackExtra.tpcNClsFound() < ConfDaughTPCnclsMin || negTrackExtra.tpcNClsFound() < ConfDaughTPCnclsMin) { + return false; + } + if (posTrackExtra.tpcCrossedRowsOverFindableCls() < rcrfc || negTrackExtra.tpcCrossedRowsOverFindableCls() < rcrfc) { + return false; + }*/ // check TPC PID if (((std::abs(posTrackExtra.tpcNSigmaPi()) > ConfDaughPIDCuts) || (std::abs(negTrackExtra.tpcNSigmaPi()) > ConfDaughPIDCuts))) { return false; From 4543138a796a090709b17fba5c378bde622d551a Mon Sep 17 00:00:00 2001 From: Prottay Das Date: Thu, 14 Aug 2025 19:32:46 +0200 Subject: [PATCH 2/4] updated --- PWGLF/Tasks/Strangeness/lambdapolsp.cxx | 6 ++++-- 1 file changed, 4 insertions(+), 2 deletions(-) diff --git a/PWGLF/Tasks/Strangeness/lambdapolsp.cxx b/PWGLF/Tasks/Strangeness/lambdapolsp.cxx index b13147da53f..afca5b2c784 100644 --- a/PWGLF/Tasks/Strangeness/lambdapolsp.cxx +++ b/PWGLF/Tasks/Strangeness/lambdapolsp.cxx @@ -588,13 +588,15 @@ struct lambdapolsp { // check TPC tracking properties if (posTrackExtra.tpcNClsCrossedRows() < cfgTPCcluster || negTrackExtra.tpcNClsCrossedRows() < cfgTPCcluster) { return false; - } /* + } + /* if (posTrackExtra.tpcNClsFound() < ConfDaughTPCnclsMin || negTrackExtra.tpcNClsFound() < ConfDaughTPCnclsMin) { return false; } if (posTrackExtra.tpcCrossedRowsOverFindableCls() < rcrfc || negTrackExtra.tpcCrossedRowsOverFindableCls() < rcrfc) { return false; - }*/ + } + */ // check TPC PID if (((std::abs(posTrackExtra.tpcNSigmaPi()) > ConfDaughPIDCuts) || (std::abs(negTrackExtra.tpcNSigmaPi()) > ConfDaughPIDCuts))) { return false; From 6a7a792ef379c34c020699da480665b7572b5198 Mon Sep 17 00:00:00 2001 From: Prottay Das Date: Sat, 16 Aug 2025 11:37:00 +0200 Subject: [PATCH 3/4] added different centrality estimators --- PWGLF/Tasks/Strangeness/lambdapolsp.cxx | 53 +++++++++++++++++-------- 1 file changed, 36 insertions(+), 17 deletions(-) diff --git a/PWGLF/Tasks/Strangeness/lambdapolsp.cxx b/PWGLF/Tasks/Strangeness/lambdapolsp.cxx index afca5b2c784..b9466d04e97 100644 --- a/PWGLF/Tasks/Strangeness/lambdapolsp.cxx +++ b/PWGLF/Tasks/Strangeness/lambdapolsp.cxx @@ -99,6 +99,7 @@ struct lambdapolsp { Configurable cfgMaxOccupancy{"cfgMaxOccupancy", 1000, "maximum occupancy of tracks in neighbouring collisions in a given time range"}; Configurable cfgMinOccupancy{"cfgMinOccupancy", 0, "maximum occupancy of tracks in neighbouring collisions in a given time range"}; Configurable sys{"sys", 1, "flag to select systematic source"}; + Configurable centestim{"centestim", 0, "flag to select centrality estimator"}; Configurable dosystematic{"dosystematic", false, "flag to perform systematic study"}; Configurable needetaaxis{"needetaaxis", false, "flag to use last axis"}; struct : ConfigurableGroup { @@ -142,8 +143,8 @@ struct lambdapolsp { Configurable ConfDaughEta{"ConfDaughEta", 0.8f, "V0 Daugh sel: max eta"}; Configurable cfgDaughPrPt{"cfgDaughPrPt", 0.4, "minimum daughter proton pt"}; Configurable cfgDaughPiPt{"cfgDaughPiPt", 0.2, "minimum daughter pion pt"}; - // Configurable rcrfc{"rcrfc", 0.8f, "Ratio of CR to FC"}; - // Configurable ConfDaughTPCnclsMin{"ConfDaughTPCnclsMin", 50.f, "V0 Daugh sel: Min. nCls TPC"}; + Configurable rcrfc{"rcrfc", 0.8f, "Ratio of CR to FC"}; + Configurable ConfDaughTPCnclsMin{"ConfDaughTPCnclsMin", 50.f, "V0 Daugh sel: Min. nCls TPC"}; Configurable ConfDaughPIDCuts{"ConfDaughPIDCuts", 3, "PID selections for Lambda daughters"}; Configurable usesubdet{"usesubdet", false, "use subdet"}; Configurable useAccCorr{"useAccCorr", false, "use acceptance correction"}; @@ -483,13 +484,13 @@ struct lambdapolsp { /*if (TMath::Abs(eta) > ConfDaughEta) { return false; }*/ - /* + if (tpcNClsF < ConfDaughTPCnclsMin) { return false; } if (track.tpcCrossedRowsOverFindableCls() < rcrfc) { return false; - }*/ + } if (pid == 0 && TMath::Abs(track.tpcNSigmaPr()) > ConfDaughPIDCuts) { return false; @@ -542,14 +543,14 @@ struct lambdapolsp { if (posTrackExtra.tpcNClsCrossedRows() < cfgTPCcluster || negTrackExtra.tpcNClsCrossedRows() < cfgTPCcluster) { return false; } - /* + if (posTrackExtra.tpcNClsFound() < ConfDaughTPCnclsMin || negTrackExtra.tpcNClsFound() < ConfDaughTPCnclsMin) { return false; } if (posTrackExtra.tpcCrossedRowsOverFindableCls() < rcrfc || negTrackExtra.tpcCrossedRowsOverFindableCls() < rcrfc) { return false; } - */ + // check TPC PID if (pid == 0 && ((std::abs(posTrackExtra.tpcNSigmaPr()) > ConfDaughPIDCuts) || (std::abs(negTrackExtra.tpcNSigmaPi()) > ConfDaughPIDCuts))) { return false; @@ -589,14 +590,14 @@ struct lambdapolsp { if (posTrackExtra.tpcNClsCrossedRows() < cfgTPCcluster || negTrackExtra.tpcNClsCrossedRows() < cfgTPCcluster) { return false; } - /* - if (posTrackExtra.tpcNClsFound() < ConfDaughTPCnclsMin || negTrackExtra.tpcNClsFound() < ConfDaughTPCnclsMin) { - return false; - } - if (posTrackExtra.tpcCrossedRowsOverFindableCls() < rcrfc || negTrackExtra.tpcCrossedRowsOverFindableCls() < rcrfc) { - return false; - } - */ + + if (posTrackExtra.tpcNClsFound() < ConfDaughTPCnclsMin || negTrackExtra.tpcNClsFound() < ConfDaughTPCnclsMin) { + return false; + } + if (posTrackExtra.tpcCrossedRowsOverFindableCls() < rcrfc || negTrackExtra.tpcCrossedRowsOverFindableCls() < rcrfc) { + return false; + } + // check TPC PID if (((std::abs(posTrackExtra.tpcNSigmaPi()) > ConfDaughPIDCuts) || (std::abs(negTrackExtra.tpcNSigmaPi()) > ConfDaughPIDCuts))) { return false; @@ -802,7 +803,16 @@ struct lambdapolsp { if (!collision.sel8()) { return; } - auto centrality = collision.centFT0C(); + double centrality; + if (centestim == 0) + centrality = collision.centFT0C(); + else if (centestim == 1) + centrality = collision.centFT0M(); + else if (centestim == 2) + centrality = collision.centFT0A(); + else if (centestim == 3) + centrality = collision.centFV0A(); + // histos.fill(HIST("hCentrality0"), centrality); if (!collision.triggereventsp()) { return; @@ -1203,7 +1213,17 @@ struct lambdapolsp { if (!collision.sel8()) { return; } - auto centrality = collision.centFT0C(); + double centrality; + if (centestim == 0) + centrality = collision.centFT0C(); + else if (centestim == 1) + centrality = collision.centFT0M(); + else if (centestim == 2) + centrality = collision.centFT0A(); + else if (centestim == 3) + centrality = collision.centFV0A(); + + // auto centrality = collision.centFT0C(); if (!collision.triggereventsp()) { // provided by StraZDCSP return; } @@ -1621,7 +1641,6 @@ struct lambdapolsp { if (additionalEvSel4 && !collision2.selection_bit(o2::aod::evsel::kIsGoodITSLayersAll)) { continue; } - auto centrality = collision1.centFT0C(); auto qxZDCA = collision1.qxZDCA(); auto qxZDCC = collision1.qxZDCC(); From 9cac2d0560abc20020c2197a1f21761f75be7caf Mon Sep 17 00:00:00 2001 From: Prottay Das Date: Sat, 16 Aug 2025 19:04:02 +0200 Subject: [PATCH 4/4] initialised cent --- PWGLF/Tasks/Strangeness/lambdapolsp.cxx | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/PWGLF/Tasks/Strangeness/lambdapolsp.cxx b/PWGLF/Tasks/Strangeness/lambdapolsp.cxx index b9466d04e97..c9af2944e34 100644 --- a/PWGLF/Tasks/Strangeness/lambdapolsp.cxx +++ b/PWGLF/Tasks/Strangeness/lambdapolsp.cxx @@ -803,7 +803,7 @@ struct lambdapolsp { if (!collision.sel8()) { return; } - double centrality; + double centrality = -999.; if (centestim == 0) centrality = collision.centFT0C(); else if (centestim == 1) @@ -1213,7 +1213,7 @@ struct lambdapolsp { if (!collision.sel8()) { return; } - double centrality; + double centrality = -999.; if (centestim == 0) centrality = collision.centFT0C(); else if (centestim == 1)