Skip to content

Commit 9df1729

Browse files
ChiaraDeMartin95Chiara De Martin
andauthored
[PWGLF] Fix selection based on ctau of cascades (#17249)
Co-authored-by: Chiara De Martin <chdemart@alicecerno2.cern.ch>
1 parent 1da3606 commit 9df1729

1 file changed

Lines changed: 5 additions & 3 deletions

File tree

PWGLF/Tasks/QC/v0cascadesqa.cxx

Lines changed: 5 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -126,6 +126,8 @@ struct v0cascadesQA {
126126
RCTFlagsChecker rctFlagsChecker{rctConfigurations.cfgRCTLabel.value};
127127

128128
static constexpr float DefaultLifetimeCuts[1][2] = {{30., 20.}};
129+
const float ctauxiPDG = 4.91; // from PDG
130+
const float ctauomegaPDG = 2.461; // from PDG
129131

130132
struct : ConfigurableGroup {
131133
std::string prefix = "v0Selections"; // JSON group name
@@ -205,7 +207,7 @@ struct v0cascadesQA {
205207
Configurable<float> dcabachtopv{"dcabachtopv", .05, "min DCA Bachelor To PV (cm)"};
206208
Configurable<float> cascradius{"cascradius", 0.5, "minimum Cascade radius (cm)"};
207209
Configurable<float> cascradiusMax{"cascradiusMax", 1E5, "maximum Cascade radius (cm)"};
208-
Configurable<float> cascProperLifeTime{"cascProperLifeTime", 3, "maximum lifetime (ctau)"};
210+
Configurable<float> cascProperLifeTime{"cascProperLifeTime", 3, "maximum lifetime (meanctau)"};
209211

210212
// invariant mass selection
211213
Configurable<float> v0MassWindow{"v0MassWindow", 0.008, "#Lambda mass (GeV/#it{c}^{2})"};
@@ -877,9 +879,9 @@ struct v0cascadesQA {
877879
//
878880
// proper lifetime
879881
float distOverTotMom = std::sqrt(std::pow(casc.x() - collision.posX(), 2) + std::pow(casc.y() - collision.posY(), 2) + std::pow(casc.z() - collision.posZ(), 2)) / (casc.p() + 1E-10);
880-
if ((cascType == kXiM || cascType == kXiP) && distOverTotMom * o2::constants::physics::MassXiMinus > cascSelections.cascProperLifeTime)
882+
if ((cascType == kXiM || cascType == kXiP) && distOverTotMom * o2::constants::physics::MassXiMinus > (cascSelections.cascProperLifeTime * ctauxiPDG))
881883
return false;
882-
if ((cascType == kOmegaM || cascType == kOmegaP) && distOverTotMom * o2::constants::physics::MassOmegaMinus > cascSelections.cascProperLifeTime)
884+
if ((cascType == kOmegaM || cascType == kOmegaP) && distOverTotMom * o2::constants::physics::MassOmegaMinus > (cascSelections.cascProperLifeTime * ctauomegaPDG))
883885
return false;
884886

885887
// rapidity

0 commit comments

Comments
 (0)