Skip to content

Commit 0247fa2

Browse files
romainschotterromainschotter
authored
[PWGLF] Add missing rapidity cut on secondary Lambda (#17660)
1 parent 66f232d commit 0247fa2

1 file changed

Lines changed: 47 additions & 41 deletions

File tree

‎PWGLF/Tasks/Strangeness/derivedlambdakzeroanalysis.cxx‎

Lines changed: 47 additions & 41 deletions
Original file line numberDiff line numberDiff line change
@@ -2985,7 +2985,7 @@ struct derivedlambdakzeroanalysis {
29852985
histos.fill(HIST("GeneralQA/h2dArmenterosAll"), v0.alpha(), v0.qtarm());
29862986

29872987
float ptmc = RecoDecay::sqrtSumOfSquares(v0MC.pxPosMC() + v0MC.pxNegMC(), v0MC.pyPosMC() + v0MC.pyNegMC());
2988-
float ymc = 1e-3;
2988+
float ymc = 1e3;
29892989
if (v0MC.pdgCode() == PDG_t::kK0Short)
29902990
ymc = RecoDecay::y(std::array{v0MC.pxPosMC() + v0MC.pxNegMC(), v0MC.pyPosMC() + v0MC.pyNegMC(), v0MC.pzPosMC() + v0MC.pzNegMC()}, o2::constants::physics::MassKaonNeutral);
29912991
else if (std::abs(v0MC.pdgCode()) == PDG_t::kLambda0)
@@ -3242,23 +3242,26 @@ struct derivedlambdakzeroanalysis {
32423242

32433243
if (doSecondaryV0s && std::abs(cascMC.pdgCodeV0()) == kLambda0) {
32443244
float v0PtMc = std::hypot(cascMC.pxPosMC() + cascMC.pxNegMC(), cascMC.pyPosMC() + cascMC.pyNegMC());
3245-
if (cascMC.pdgCodeV0() == kLambda0) {
3246-
histos.fill(HIST("h2dGenSecLambdaVsMultMC_RecoedEvt"), mcCollision.multMCNParticlesEta05(), v0PtMc);
3247-
}
3248-
if (cascMC.pdgCodeV0() == kLambda0Bar) {
3249-
histos.fill(HIST("h2dGenSecAntiLambdaVsMultMC_RecoedEvt"), mcCollision.multMCNParticlesEta05(), v0PtMc);
3250-
}
3251-
if (cascMC.pdgCode() == PDG_t::kXiMinus && cascMC.pdgCodeV0() == kLambda0) {
3252-
histos.fill(HIST("h2dGenSecLambdaFromXiVsMultMC_RecoedEvt"), mcCollision.multMCNParticlesEta05(), v0PtMc);
3253-
}
3254-
if (cascMC.pdgCode() == PDG_t::kXiPlusBar && cascMC.pdgCodeV0() == kLambda0Bar) {
3255-
histos.fill(HIST("h2dGenSecAntiLambdaFromXiVsMultMC_RecoedEvt"), mcCollision.multMCNParticlesEta05(), v0PtMc);
3256-
}
3257-
if (cascMC.pdgCode() == PDG_t::kOmegaMinus && cascMC.pdgCodeV0() == kLambda0) {
3258-
histos.fill(HIST("h2dGenSecLambdaFromOmegaVsMultMC_RecoedEvt"), mcCollision.multMCNParticlesEta05(), v0PtMc);
3259-
}
3260-
if (cascMC.pdgCode() == PDG_t::kOmegaPlusBar && cascMC.pdgCodeV0() == kLambda0Bar) {
3261-
histos.fill(HIST("h2dGenSecAntiLambdaFromOmegaVsMultMC_RecoedEvt"), mcCollision.multMCNParticlesEta05(), v0PtMc);
3245+
float v0RapMc = RecoDecay::y(std::array{cascMC.pxPosMC() + cascMC.pxNegMC(), cascMC.pyPosMC() + cascMC.pyNegMC(), cascMC.pzPosMC() + cascMC.pzNegMC()}, o2::constants::physics::MassLambda);
3246+
if (v0Selections.rapidityMinCut < v0RapMc && v0RapMc < v0Selections.rapidityMaxCut) {
3247+
if (cascMC.pdgCodeV0() == kLambda0) {
3248+
histos.fill(HIST("h2dGenSecLambdaVsMultMC_RecoedEvt"), mcCollision.multMCNParticlesEta05(), v0PtMc);
3249+
}
3250+
if (cascMC.pdgCodeV0() == kLambda0Bar) {
3251+
histos.fill(HIST("h2dGenSecAntiLambdaVsMultMC_RecoedEvt"), mcCollision.multMCNParticlesEta05(), v0PtMc);
3252+
}
3253+
if (cascMC.pdgCode() == PDG_t::kXiMinus && cascMC.pdgCodeV0() == kLambda0) {
3254+
histos.fill(HIST("h2dGenSecLambdaFromXiVsMultMC_RecoedEvt"), mcCollision.multMCNParticlesEta05(), v0PtMc);
3255+
}
3256+
if (cascMC.pdgCode() == PDG_t::kXiPlusBar && cascMC.pdgCodeV0() == kLambda0Bar) {
3257+
histos.fill(HIST("h2dGenSecAntiLambdaFromXiVsMultMC_RecoedEvt"), mcCollision.multMCNParticlesEta05(), v0PtMc);
3258+
}
3259+
if (cascMC.pdgCode() == PDG_t::kOmegaMinus && cascMC.pdgCodeV0() == kLambda0) {
3260+
histos.fill(HIST("h2dGenSecLambdaFromOmegaVsMultMC_RecoedEvt"), mcCollision.multMCNParticlesEta05(), v0PtMc);
3261+
}
3262+
if (cascMC.pdgCode() == PDG_t::kOmegaPlusBar && cascMC.pdgCodeV0() == kLambda0Bar) {
3263+
histos.fill(HIST("h2dGenSecAntiLambdaFromOmegaVsMultMC_RecoedEvt"), mcCollision.multMCNParticlesEta05(), v0PtMc);
3264+
}
32623265
}
32633266
}
32643267
}
@@ -3314,29 +3317,32 @@ struct derivedlambdakzeroanalysis {
33143317

33153318
if (doSecondaryV0s && std::abs(cascMC.pdgCodeV0()) == kLambda0) {
33163319
float v0PtMc = std::hypot(cascMC.pxPosMC() + cascMC.pxNegMC(), cascMC.pyPosMC() + cascMC.pyNegMC());
3317-
if (cascMC.pdgCodeV0() == kLambda0) {
3318-
histos.fill(HIST("h2dGenSecLambda"), centrality, v0PtMc);
3319-
histos.fill(HIST("h2dGenSecLambdaVsMultMC"), mcCollision.multMCNParticlesEta05(), v0PtMc);
3320-
}
3321-
if (cascMC.pdgCodeV0() == kLambda0Bar) {
3322-
histos.fill(HIST("h2dGenSecAntiLambda"), centrality, v0PtMc);
3323-
histos.fill(HIST("h2dGenSecAntiLambdaVsMultMC"), mcCollision.multMCNParticlesEta05(), v0PtMc);
3324-
}
3325-
if (cascMC.pdgCode() == PDG_t::kXiMinus && cascMC.pdgCodeV0() == kLambda0) {
3326-
histos.fill(HIST("h2dGenSecLambdaFromXi"), centrality, v0PtMc);
3327-
histos.fill(HIST("h2dGenSecLambdaFromXiVsMultMC"), mcCollision.multMCNParticlesEta05(), v0PtMc);
3328-
}
3329-
if (cascMC.pdgCode() == PDG_t::kXiPlusBar && cascMC.pdgCodeV0() == kLambda0Bar) {
3330-
histos.fill(HIST("h2dGenSecAntiLambdaFromXi"), centrality, v0PtMc);
3331-
histos.fill(HIST("h2dGenSecAntiLambdaFromXiVsMultMC"), mcCollision.multMCNParticlesEta05(), v0PtMc);
3332-
}
3333-
if (cascMC.pdgCode() == PDG_t::kOmegaMinus && cascMC.pdgCodeV0() == kLambda0) {
3334-
histos.fill(HIST("h2dGenSecLambdaFromOmega"), centrality, v0PtMc);
3335-
histos.fill(HIST("h2dGenSecLambdaFromOmegaVsMultMC"), mcCollision.multMCNParticlesEta05(), v0PtMc);
3336-
}
3337-
if (cascMC.pdgCode() == PDG_t::kOmegaPlusBar && cascMC.pdgCodeV0() == kLambda0Bar) {
3338-
histos.fill(HIST("h2dGenSecAntiLambdaFromOmega"), centrality, v0PtMc);
3339-
histos.fill(HIST("h2dGenSecAntiLambdaFromOmegaVsMultMC"), mcCollision.multMCNParticlesEta05(), v0PtMc);
3320+
float v0RapMc = RecoDecay::y(std::array{cascMC.pxPosMC() + cascMC.pxNegMC(), cascMC.pyPosMC() + cascMC.pyNegMC(), cascMC.pzPosMC() + cascMC.pzNegMC()}, o2::constants::physics::MassLambda);
3321+
if (v0Selections.rapidityMinCut < v0RapMc && v0RapMc < v0Selections.rapidityMaxCut) {
3322+
if (cascMC.pdgCodeV0() == kLambda0) {
3323+
histos.fill(HIST("h2dGenSecLambda"), centrality, v0PtMc);
3324+
histos.fill(HIST("h2dGenSecLambdaVsMultMC"), mcCollision.multMCNParticlesEta05(), v0PtMc);
3325+
}
3326+
if (cascMC.pdgCodeV0() == kLambda0Bar) {
3327+
histos.fill(HIST("h2dGenSecAntiLambda"), centrality, v0PtMc);
3328+
histos.fill(HIST("h2dGenSecAntiLambdaVsMultMC"), mcCollision.multMCNParticlesEta05(), v0PtMc);
3329+
}
3330+
if (cascMC.pdgCode() == PDG_t::kXiMinus && cascMC.pdgCodeV0() == kLambda0) {
3331+
histos.fill(HIST("h2dGenSecLambdaFromXi"), centrality, v0PtMc);
3332+
histos.fill(HIST("h2dGenSecLambdaFromXiVsMultMC"), mcCollision.multMCNParticlesEta05(), v0PtMc);
3333+
}
3334+
if (cascMC.pdgCode() == PDG_t::kXiPlusBar && cascMC.pdgCodeV0() == kLambda0Bar) {
3335+
histos.fill(HIST("h2dGenSecAntiLambdaFromXi"), centrality, v0PtMc);
3336+
histos.fill(HIST("h2dGenSecAntiLambdaFromXiVsMultMC"), mcCollision.multMCNParticlesEta05(), v0PtMc);
3337+
}
3338+
if (cascMC.pdgCode() == PDG_t::kOmegaMinus && cascMC.pdgCodeV0() == kLambda0) {
3339+
histos.fill(HIST("h2dGenSecLambdaFromOmega"), centrality, v0PtMc);
3340+
histos.fill(HIST("h2dGenSecLambdaFromOmegaVsMultMC"), mcCollision.multMCNParticlesEta05(), v0PtMc);
3341+
}
3342+
if (cascMC.pdgCode() == PDG_t::kOmegaPlusBar && cascMC.pdgCodeV0() == kLambda0Bar) {
3343+
histos.fill(HIST("h2dGenSecAntiLambdaFromOmega"), centrality, v0PtMc);
3344+
histos.fill(HIST("h2dGenSecAntiLambdaFromOmegaVsMultMC"), mcCollision.multMCNParticlesEta05(), v0PtMc);
3345+
}
33403346
}
33413347
}
33423348
}

0 commit comments

Comments
 (0)