Skip to content

Commit 7eb9af2

Browse files
author
Chiara De Martin
committed
fix selection based on ctau of cascades
1 parent c5313c8 commit 7eb9af2

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)