Skip to content

Commit 67f5ed3

Browse files
fmazzascFrancesco Mazzaschi
andauthored
Add QA histogram for non-prompt nuclei (#11437)
Co-authored-by: Francesco Mazzaschi <fmazzasc@alipap1.cern.ch>
1 parent 9e330c5 commit 67f5ed3

File tree

1 file changed

+3
-0
lines changed

1 file changed

+3
-0
lines changed

PWGLF/TableProducer/Nuspex/nucleiSpectra.cxx

Lines changed: 3 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -516,6 +516,7 @@ struct nucleiSpectra {
516516
spectra.add("hTpcSignalData", "Specific energy loss", HistType::kTH2F, {{600, -6., 6., "#it{p} (GeV/#it{c})"}, {1400, 0, 1400, "d#it{E} / d#it{X} (a. u.)"}});
517517
spectra.add("hTpcSignalDataSelected", "Specific energy loss for selected particles", HistType::kTH2F, {{600, -6., 6., "#it{p} (GeV/#it{c})"}, {1400, 0, 1400, "d#it{E} / d#it{X} (a. u.)"}});
518518
spectra.add("hTofSignalData", "TOF beta", HistType::kTH2F, {{500, 0., 5., "#it{p} (GeV/#it{c})"}, {750, 0, 1.5, "TOF #beta"}});
519+
spectra.add("hMCDCAxyVsDecRadius", "DCA xy vs DecayRadius", HistType::kTH2F, {{20, 0., 0.2, "Decay radius"}, {150, -0.05, 0.05, "DCA_{xy} (cm)"}});
519520
for (int iC{0}; iC < 2; ++iC) {
520521
nuclei::hGloTOFtracks[iC] = spectra.add<TH2>(fmt::format("hTPCTOFtracks{}", nuclei::matter[iC]).data(), fmt::format("Global vs TOF matched {} tracks in a collision", nuclei::chargeLabelNames[iC]).data(), HistType::kTH2D, {{300, -0.5, 300.5, "Number of global tracks"}, {300, -0.5, 300.5, "Number of TOF matched tracks"}});
521522

@@ -991,6 +992,8 @@ struct nucleiSpectra {
991992
if (std::find(nuclei::hfMothCodes.begin(), nuclei::hfMothCodes.end(), std::abs(motherparticle.pdgCode())) != nuclei::hfMothCodes.end()) {
992993
c.flags |= kIsSecondaryFromWeakDecay;
993994
MotherpdgCode = motherparticle.pdgCode();
995+
float decRadius = std::hypot(particle.vx() - motherparticle.vx(), particle.vy() - motherparticle.vy());
996+
spectra.fill(HIST("hMCDCAxyVsDecRadius"), decRadius, c.DCAxy);
994997
break;
995998
}
996999
}

0 commit comments

Comments
 (0)