Skip to content

Commit a7a3b02

Browse files
ariedel-cernariedel-cern
committed
Feat: add lambda mass to cascade extras
1 parent 1d671dd commit a7a3b02

3 files changed

Lines changed: 39 additions & 23 deletions

File tree

‎PWGCF/Femto/Core/cascadeBuilder.h‎

Lines changed: 4 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -255,10 +255,10 @@ class CascadeSelection : public baseselection::BaseSelection<float, o2::analysis
255255

256256
this->addSelection(kCascadeCpaMin, cascadeSelectionNames.at(kCascadeCpaMin), config.cascadeCpaMin.value, limits::kLowerLimit, true, true, false);
257257
this->addSelection(kCascadeTransRadMin, cascadeSelectionNames.at(kCascadeTransRadMin), config.cascadeTransRadMin.value, limits::kLowerLimit, true, true, false);
258-
this->addSelection(kCascadeDcaDaughMax, cascadeSelectionNames.at(kCascadeDcaDaughMax), config.cascadeDcaDauMax.value, limits::kAbsUpperLimit, true, true, false);
258+
this->addSelection(kCascadeDcaDaughMax, cascadeSelectionNames.at(kCascadeDcaDaughMax), config.cascadeDcaDauMax.value, limits::kUpperLimit, true, true, false);
259259
this->addSelection(kLambdaCpaMin, cascadeSelectionNames.at(kLambdaCpaMin), config.lambdaCpaMin.value, limits::kLowerLimit, true, true, false);
260260
this->addSelection(kLambdaTransRadMin, cascadeSelectionNames.at(kLambdaTransRadMin), config.lambdaTransRadMin.value, limits::kLowerLimit, true, true, false);
261-
this->addSelection(kLambdaDcaDauMax, cascadeSelectionNames.at(kLambdaDcaDauMax), config.lambdaDcaDauMax.value, limits::kAbsUpperLimit, true, true, false);
261+
this->addSelection(kLambdaDcaDauMax, cascadeSelectionNames.at(kLambdaDcaDauMax), config.lambdaDcaDauMax.value, limits::kUpperLimit, true, true, false);
262262
this->addSelection(kLambdaDcaToPvMin, cascadeSelectionNames.at(kLambdaDcaToPvMin), config.lambdaDcaToPvMin.value, limits::kLowerLimit, true, true, false);
263263
this->addSelection(kDauAbsEtaMax, cascadeSelectionNames.at(kDauAbsEtaMax), config.dauAbsEtaMax.value, limits::kAbsUpperLimit, true, true, false);
264264
this->addSelection(kDauDcaMin, cascadeSelectionNames.at(kDauDcaMin), config.dauDcaMin.value, limits::kAbsLowerLimit, true, true, false);
@@ -592,6 +592,7 @@ class CascadeBuilder
592592
cascade.casccosPA(col.posX(), col.posY(), col.posZ()),
593593
cascade.dcacascdaughters(),
594594
cascade.cascradius(),
595+
cascade.mLambda(),
595596
cascade.v0cosPA(col.posX(), col.posY(), col.posZ()),
596597
cascade.dcaV0daughters(),
597598
cascade.v0radius(),
@@ -618,6 +619,7 @@ class CascadeBuilder
618619
cascade.casccosPA(col.posX(), col.posY(), col.posZ()),
619620
cascade.dcacascdaughters(),
620621
cascade.cascradius(),
622+
cascade.mLambda(),
621623
cascade.v0cosPA(col.posX(), col.posY(), col.posZ()),
622624
cascade.dcaV0daughters(),
623625
cascade.v0radius(),

‎PWGCF/Femto/Core/cascadeHistManager.h‎

Lines changed: 32 additions & 21 deletions
Original file line numberDiff line numberDiff line change
@@ -53,6 +53,7 @@ enum CascadeHist {
5353
kCosPa,
5454
kDecayDauDca,
5555
kTransRadius,
56+
kLambdaMass,
5657
kLambdaCosPa,
5758
kLambdaDauDca,
5859
kLambdaTransRadius,
@@ -64,6 +65,7 @@ enum CascadeHist {
6465
kPtVsCosPa,
6566
kPtVsMassXi,
6667
kPtVsMassOmega,
68+
kPtVsMassLambda,
6769
kMassXiVsMassOmega,
6870
// mc
6971
kOrigin,
@@ -115,7 +117,8 @@ struct ConfCascadeQaBinning : o2::framework::ConfigurableGroup {
115117
o2::framework::ConfigurableAxis transRadius{"transRadius", {{100, 0, 100}}, "Transverse radius"};
116118
o2::framework::ConfigurableAxis massXi{"massXi", {{400, 1.2f, 1.6f}}, "mass for antiparticle hypothesis"};
117119
o2::framework::ConfigurableAxis massOmega{"massOmega", {{400, 1.4f, 1.8f}}, "mass for antiparticle hypothesis"};
118-
o2::framework::ConfigurableAxis lambdaCosPa{"lambdaCosPa", {{100, 0.9, 1}}, "Cosine of poiting angle of daughter lambda"};
120+
o2::framework::ConfigurableAxis lambdaMass{"lambdaMass", {{200, 1., 1.2}}, "Mass of daughter lambda"};
121+
o2::framework::ConfigurableAxis lambdaCosPa{"lambdaCosPa", {{100, 0.9, 1}}, "Cosine of pointing angle of daughter lambda"};
119122
o2::framework::ConfigurableAxis lambdaDauDca{"lambdaDauDca", {{150, 0, 1.0}}, "DCA of lambda daughters at lambda decay vertex"};
120123
o2::framework::ConfigurableAxis lambdaTransRadius{"lambdaTransRadius", {{100, 0, 100}}, "DCA of lambda daughters at lambda decay vertex"};
121124
o2::framework::ConfigurableAxis lambdaDcaToPv{"lambdaDcaToPv", {{100, 0, 200}}, "DCA of lambda daughter from primary vertex"};
@@ -141,16 +144,18 @@ constexpr std::array<histmanager::HistInfo<CascadeHist>, kCascadeHistLast> HistT
141144
{kCosPa, o2::framework::HistType::kTH1F, "hCosPa", "Cosine of pointing angle; cos(#alpha); Entries"},
142145
{kDecayDauDca, o2::framework::HistType::kTH1F, "hDauDca", "Daughter DCA at decay vertex ; DCA_{Decay vertex} (cm); Entries"},
143146
{kTransRadius, o2::framework::HistType::kTH1F, "hTransRadius", "Transverse radius ; r_{xy} (cm); Entries"},
144-
{kLambdaCosPa, o2::framework::HistType::kTH1F, "hLambdaCosPa", "Cosine of poiting angle of daughter lambda ; cos_{#Lambda}(#alpha); Entries"},
145-
{kLambdaDauDca, o2::framework::HistType::kTH1F, "hLambdaDauDca", "Daughter DCA at #Lambda decay vertex ; DCA_{#Lambda decay vertex} (cm); Entries"},
146-
{kLambdaTransRadius, o2::framework::HistType::kTH1F, "hLambdaTransRadius", "Transverse radius of daughter #Lambda ; r_{xy} (cm); Entries"},
147-
{kLambdaDcaToPv, o2::framework::HistType::kTH1F, "hLambdaDcaToPv", "DCA to primary vertex of daughter #Lambda ; DCA (cm); Entries"},
147+
{kLambdaMass, o2::framework::HistType::kTH1F, "hLambdaMass", "Invariant mass of daughter lambda ; m_{#Lambda dau} (GeV/#it{c}^{2}); Entries"},
148+
{kLambdaCosPa, o2::framework::HistType::kTH1F, "hLambdaCosPa", "Cosine of poiting angle of daughter lambda ; cos_{#Lambda dau}(#alpha); Entries"},
149+
{kLambdaDauDca, o2::framework::HistType::kTH1F, "hLambdaDauDca", "Daughter DCA at #Lambda decay vertex ; DCA_{#Lambda dau decay vertex} (cm); Entries"},
150+
{kLambdaTransRadius, o2::framework::HistType::kTH1F, "hLambdaTransRadius", "Transverse radius of daughter #Lambda ; r_{xy,#Lambda dau} (cm); Entries"},
151+
{kLambdaDcaToPv, o2::framework::HistType::kTH1F, "hLambdaDcaToPv", "DCA to primary vertex of daughter #Lambda ; DCA_{#Lambda} (cm); Entries"},
148152
{kPtVsEta, o2::framework::HistType::kTH2F, "hPtVsEta", "p_{T} vs #eta; p_{T} (GeV/#it{c}) ; #eta"},
149153
{kPtVsPhi, o2::framework::HistType::kTH2F, "hPtVsPhi", "p_{T} vs #varphi; p_{T} (GeV/#it{c}) ; #varphi"},
150154
{kPhiVsEta, o2::framework::HistType::kTH2F, "hPhiVsEta", "#varphi vs #eta; #varphi ; #eta"},
151155
{kPtVsCosPa, o2::framework::HistType::kTH2F, "hPtVsCosPa", "Cosine of poiting angle vs p_{T}; cos(#alpha); p_{T} (GeV/#it{c})"},
152156
{kPtVsMassXi, o2::framework::HistType::kTH2F, "hPtVsMassXi", "p_{T} vs mass #Xi; p_{T} (GeV/#it{c}); m_{#Lambda#pi} (GeV/#it{c}^{2})"},
153157
{kPtVsMassOmega, o2::framework::HistType::kTH2F, "hPtVsMassOmega", "p_{T} vs mass #Omega; p_{T} (GeV/#it{c}); m_{#LambdaK} (GeV/#it{c}^{2})"},
158+
{kPtVsMassLambda, o2::framework::HistType::kTH2F, "hPtVsMassLambda", "p_{T} vs mass daughter #Lambda; p_{T} (GeV/#it{c}); m_{#Lambda dau} (GeV/#it{c}^{2})"},
154159
{kMassXiVsMassOmega, o2::framework::HistType::kTH2F, "hMassXiVsMassOmega", "mass #Xi vs mass #Omega; m_{#Lambda#pi} (GeV/#it{c}^{2}); m_{#LambdaK} (GeV/#it{c}^{2})"},
155160
{kOrigin, o2::framework::HistType::kTH1F, "hOrigin", "Status Codes (=Origin); Status Code; Entries"},
156161
{kPdg, o2::framework::HistType::kTH1F, "hPdg", "PDG Codes of reconstructed v0; PDG Code; Entries"},
@@ -189,22 +194,24 @@ constexpr std::array<histmanager::HistInfo<CascadeHist>, kCascadeHistLast> HistT
189194
{kPdgPartonicMother, {(conf).pdgCodes}},
190195

191196
// NOLINTNEXTLINE(cppcoreguidelines-macro-usage)
192-
#define CASCADE_HIST_QA_MAP(confAnalysis, confQa) \
193-
{kCosPa, {(confQa).cosPa}}, \
194-
{kDecayDauDca, {(confQa).dauDcaAtDecay}}, \
195-
{kTransRadius, {(confQa).transRadius}}, \
196-
{kLambdaCosPa, {(confQa).lambdaCosPa}}, \
197-
{kLambdaDauDca, {(confQa).lambdaDauDca}}, \
198-
{kLambdaTransRadius, {(confQa).lambdaTransRadius}}, \
199-
{kLambdaDcaToPv, {(confQa).lambdaDcaToPv}}, \
200-
{kPtVsEta, {(confAnalysis).pt, (confAnalysis).eta}}, \
201-
{kPtVsPhi, {(confAnalysis).pt, (confAnalysis).phi}}, \
202-
{kPhiVsEta, {(confAnalysis).phi, (confAnalysis).eta}}, \
203-
{kPtVsCosPa, {(confAnalysis).pt, (confQa).cosPa}}, \
204-
{kMassXi, {(confQa).massXi}}, \
205-
{kMassOmega, {(confQa).massOmega}}, \
206-
{kPtVsMassXi, {(confAnalysis).pt, (confQa).massXi}}, \
207-
{kPtVsMassOmega, {(confAnalysis).pt, (confQa).massOmega}}, \
197+
#define CASCADE_HIST_QA_MAP(confAnalysis, confQa) \
198+
{kCosPa, {(confQa).cosPa}}, \
199+
{kDecayDauDca, {(confQa).dauDcaAtDecay}}, \
200+
{kTransRadius, {(confQa).transRadius}}, \
201+
{kLambdaMass, {(confQa).lambdaMass}}, \
202+
{kLambdaCosPa, {(confQa).lambdaCosPa}}, \
203+
{kLambdaDauDca, {(confQa).lambdaDauDca}}, \
204+
{kLambdaTransRadius, {(confQa).lambdaTransRadius}}, \
205+
{kLambdaDcaToPv, {(confQa).lambdaDcaToPv}}, \
206+
{kPtVsEta, {(confAnalysis).pt, (confAnalysis).eta}}, \
207+
{kPtVsPhi, {(confAnalysis).pt, (confAnalysis).phi}}, \
208+
{kPhiVsEta, {(confAnalysis).phi, (confAnalysis).eta}}, \
209+
{kPtVsCosPa, {(confAnalysis).pt, (confQa).cosPa}}, \
210+
{kMassXi, {(confQa).massXi}}, \
211+
{kMassOmega, {(confQa).massOmega}}, \
212+
{kPtVsMassXi, {(confAnalysis).pt, (confQa).massXi}}, \
213+
{kPtVsMassOmega, {(confAnalysis).pt, (confQa).massOmega}}, \
214+
{kPtVsMassLambda, {(confAnalysis).pt, (confQa).lambdaMass}}, \
208215
{kMassXiVsMassOmega, {(confQa).massXi, (confQa).massOmega}},
209216

210217
// NOLINTNEXTLINE(cppcoreguidelines-macro-usage)
@@ -482,6 +489,7 @@ class CascadeHistManager
482489
mHistogramRegistry->add(qaDir + getHistNameV2(kCosPa, HistTable), getHistDesc(kCosPa, HistTable), getHistType(kCosPa, HistTable), {cascadeSpecs.at(kCosPa)});
483490
mHistogramRegistry->add(qaDir + getHistNameV2(kDecayDauDca, HistTable), getHistDesc(kDecayDauDca, HistTable), getHistType(kDecayDauDca, HistTable), {cascadeSpecs.at(kDecayDauDca)});
484491
mHistogramRegistry->add(qaDir + getHistNameV2(kTransRadius, HistTable), getHistDesc(kTransRadius, HistTable), getHistType(kTransRadius, HistTable), {cascadeSpecs.at(kTransRadius)});
492+
mHistogramRegistry->add(qaDir + getHistNameV2(kLambdaMass, HistTable), getHistDesc(kLambdaMass, HistTable), getHistType(kLambdaMass, HistTable), {cascadeSpecs.at(kLambdaMass)});
485493
mHistogramRegistry->add(qaDir + getHistNameV2(kLambdaCosPa, HistTable), getHistDesc(kLambdaCosPa, HistTable), getHistType(kLambdaCosPa, HistTable), {cascadeSpecs.at(kLambdaCosPa)});
486494
mHistogramRegistry->add(qaDir + getHistNameV2(kLambdaDauDca, HistTable), getHistDesc(kLambdaDauDca, HistTable), getHistType(kLambdaDauDca, HistTable), {cascadeSpecs.at(kLambdaDauDca)});
487495
mHistogramRegistry->add(qaDir + getHistNameV2(kLambdaTransRadius, HistTable), getHistDesc(kLambdaTransRadius, HistTable), getHistType(kLambdaTransRadius, HistTable), {cascadeSpecs.at(kLambdaTransRadius)});
@@ -497,6 +505,7 @@ class CascadeHistManager
497505

498506
mHistogramRegistry->add(qaDir + getHistNameV2(kPtVsMassXi, HistTable), getHistDesc(kPtVsMassXi, HistTable), getHistType(kPtVsMassXi, HistTable), {cascadeSpecs.at(kPtVsMassXi)});
499507
mHistogramRegistry->add(qaDir + getHistNameV2(kPtVsMassOmega, HistTable), getHistDesc(kPtVsMassOmega, HistTable), getHistType(kPtVsMassOmega, HistTable), {cascadeSpecs.at(kPtVsMassOmega)});
508+
mHistogramRegistry->add(qaDir + getHistNameV2(kPtVsMassLambda, HistTable), getHistDesc(kPtVsMassLambda, HistTable), getHistType(kPtVsMassLambda, HistTable), {cascadeSpecs.at(kPtVsMassLambda)});
500509
mHistogramRegistry->add(qaDir + getHistNameV2(kMassXiVsMassOmega, HistTable), getHistDesc(kMassXiVsMassOmega, HistTable), getHistType(kMassXiVsMassOmega, HistTable), {cascadeSpecs.at(kMassXiVsMassOmega)});
501510
}
502511
}
@@ -569,6 +578,7 @@ class CascadeHistManager
569578
mHistogramRegistry->fill(HIST(cascadePrefix) + HIST(QaDir) + HIST(getHistName(kCosPa, HistTable)), cascadeCandidate.cascadeCosPa());
570579
mHistogramRegistry->fill(HIST(cascadePrefix) + HIST(QaDir) + HIST(getHistName(kDecayDauDca, HistTable)), cascadeCandidate.cascadeDauDca());
571580
mHistogramRegistry->fill(HIST(cascadePrefix) + HIST(QaDir) + HIST(getHistName(kTransRadius, HistTable)), cascadeCandidate.cascadeTransRadius());
581+
mHistogramRegistry->fill(HIST(cascadePrefix) + HIST(QaDir) + HIST(getHistName(kLambdaMass, HistTable)), cascadeCandidate.lambdaMass());
572582
mHistogramRegistry->fill(HIST(cascadePrefix) + HIST(QaDir) + HIST(getHistName(kLambdaCosPa, HistTable)), cascadeCandidate.lambdaCosPa());
573583
mHistogramRegistry->fill(HIST(cascadePrefix) + HIST(QaDir) + HIST(getHistName(kLambdaDauDca, HistTable)), cascadeCandidate.lambdaDauDca());
574584
mHistogramRegistry->fill(HIST(cascadePrefix) + HIST(QaDir) + HIST(getHistName(kLambdaTransRadius, HistTable)), cascadeCandidate.lambdaTransRadius());
@@ -583,6 +593,7 @@ class CascadeHistManager
583593
mHistogramRegistry->fill(HIST(cascadePrefix) + HIST(QaDir) + HIST(getHistName(kPtVsCosPa, HistTable)), cascadeCandidate.pt(), cascadeCandidate.cascadeCosPa());
584594
mHistogramRegistry->fill(HIST(cascadePrefix) + HIST(QaDir) + HIST(getHistName(kPtVsMassXi, HistTable)), cascadeCandidate.pt(), massXi);
585595
mHistogramRegistry->fill(HIST(cascadePrefix) + HIST(QaDir) + HIST(getHistName(kPtVsMassOmega, HistTable)), cascadeCandidate.pt(), massOmega);
596+
mHistogramRegistry->fill(HIST(cascadePrefix) + HIST(QaDir) + HIST(getHistName(kPtVsMassLambda, HistTable)), cascadeCandidate.pt(), cascadeCandidate.lambdaMass());
586597
mHistogramRegistry->fill(HIST(cascadePrefix) + HIST(QaDir) + HIST(getHistName(kMassXiVsMassOmega, HistTable)), massXi, massOmega);
587598
}
588599
}

‎PWGCF/Femto/DataModel/FemtoTables.h‎

Lines changed: 3 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -625,6 +625,7 @@ DECLARE_SOA_COLUMN(MassOmega, massOmega, float); //! Mass of o
625625
DECLARE_SOA_COLUMN(CascadeCosPa, cascadeCosPa, float); //! cosine of the poiting angle at decay vertex
626626
DECLARE_SOA_COLUMN(CascadeDauDca, cascadeDauDca, float); //! Lambda daughter DCA at decay vertex
627627
DECLARE_SOA_COLUMN(CascadeTransRadius, cascadeTransRadius, float); //! Lambda transvers radius
628+
DECLARE_SOA_COLUMN(LambdaMass, lambdaMass, float); //! Lambda daughter mass
628629
DECLARE_SOA_COLUMN(LambdaCosPa, lambdaCosPa, float); //! cosine of the poiting angle at decay vertex
629630
DECLARE_SOA_COLUMN(LambdaDauDca, lambdaDauDca, float); //! Lambda daughter DCA at decay vertex
630631
DECLARE_SOA_COLUMN(LambdaTransRadius, lambdaTransRadius, float); //! Lambda transvers radius
@@ -666,6 +667,7 @@ DECLARE_SOA_TABLE_STAGED_VERSIONED(FXiExtras_001, "FXIEXTRA", 1, //! xi extra in
666667
femtocascades::CascadeCosPa,
667668
femtocascades::CascadeDauDca,
668669
femtocascades::CascadeTransRadius,
670+
femtocascades::LambdaMass,
669671
femtocascades::LambdaCosPa,
670672
femtocascades::LambdaDauDca,
671673
femtocascades::LambdaTransRadius,
@@ -702,6 +704,7 @@ DECLARE_SOA_TABLE_STAGED_VERSIONED(FOmegaExtras_001, "FOMEGAEXTRA", 1, //! omega
702704
femtocascades::CascadeCosPa,
703705
femtocascades::CascadeDauDca,
704706
femtocascades::CascadeTransRadius,
707+
femtocascades::LambdaMass,
705708
femtocascades::LambdaCosPa,
706709
femtocascades::LambdaDauDca,
707710
femtocascades::LambdaTransRadius,

0 commit comments

Comments
 (0)