Skip to content

Commit ea67991

Browse files
basiachBarbara Chytla
andauthored
[PWGCF] FemtoUniverse - Adding a pT vs invM vs mult histogram to Cascade task (#11221)
Co-authored-by: Barbara Chytla <barbara.chytla@cern.ch>
1 parent 4d2d5b1 commit ea67991

File tree

1 file changed

+5
-0
lines changed

1 file changed

+5
-0
lines changed

PWGCF/FemtoUniverse/Tasks/femtoUniversePairTaskTrackCascadeExtended.cxx

Lines changed: 5 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -211,6 +211,7 @@ struct femtoUniversePairTaskTrackCascadeExtended {
211211

212212
void init(InitContext const&)
213213
{
214+
std::vector<double> multBinning = {0.0, 5.0, 10.0, 20.0, 30.0f, 40.0, 50.0, 60.0f, 70.0, 80.0, 100.0, 200.0, 99999.0};
214215
// Axes
215216
AxisSpec aXiMassAxis = {200, 1.28f, 1.36f, "#it{M}_{inv} [GeV/#it{c}^{2}]"};
216217
AxisSpec ptAxis = {100, 0.0f, 10.0f, "#it{p}_{T} (GeV/#it{c})"};
@@ -220,6 +221,7 @@ struct femtoUniversePairTaskTrackCascadeExtended {
220221
AxisSpec aCPAAxis = {1000, 0.95f, 1.0f, "#it{cos #theta_{p}}"};
221222
AxisSpec tranRadAxis = {1000, 0.0f, 100.0f, "#it{r}_{xy} (cm)"};
222223
AxisSpec aDCAToPVAxis = {1000, -10.0f, 10.0f, "DCA to PV (cm)"};
224+
AxisSpec multAxis = {multBinning, "Multiplicity"};
223225

224226
// Histograms
225227
rXiQA.add("hMassXi", "hMassXi", {HistType::kTH1F, {aXiMassAxis}});
@@ -238,6 +240,7 @@ struct femtoUniversePairTaskTrackCascadeExtended {
238240
rXiQA.add("hDcaBachtoPV", "hDcaBachtoPV", {HistType::kTH1F, {aDCAToPVAxis}});
239241
rXiQA.add("hDcaV0toPV", "hDcaV0toPV", {HistType::kTH1F, {aDCAToPVAxis}});
240242
rXiQA.add("hInvMpT", "hInvMpT", kTH2F, {{ptAxis}, {aXiMassAxis}});
243+
rXiQA.add("hInvMpTmult", "hInvMpTmult", kTH3F, {{ptAxis}, {aXiMassAxis}, {multAxis}});
241244

242245
eventHisto.init(&qaRegistry);
243246
qaRegistry.add("Tracks_pos/nSigmaTPC", "; #it{p} (GeV/#it{c}); n#sigma_{TPC}", kTH2F, {{100, 0, 10}, {200, -4.975, 5.025}});
@@ -377,6 +380,8 @@ struct femtoUniversePairTaskTrackCascadeExtended {
377380
posChildHistos.fillQA<false, true>(posChild);
378381
negChildHistos.fillQA<false, true>(negChild);
379382
bachHistos.fillQABase<false, true>(bachelor, HIST("hBachelor"));
383+
384+
rXiQA.fill(HIST("hInvMpTmult"), part.pt(), part.mLambda(), multCol);
380385
}
381386

382387
for (const auto& part : groupPartsOne) {

0 commit comments

Comments
 (0)