Skip to content

[PWGLF] Update datamodel of hyphe4s analysis #11811

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
maciacco merged 2 commits into from
Jun 28, 2025
Merged

[PWGLF] Update datamodel of hyphe4s analysis #11811

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
c5eded1
Update datamodel of hyphe4s analysis
wang-yuanzhe Jun 27, 2025
462f1a6
Fix processData
wang-yuanzhe Jun 27, 2025
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
15 changes: 13 additions & 2 deletions PWGLF/DataModel/LFHyperhelium4sigmaTables.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -30,6 +30,9 @@ DECLARE_SOA_COLUMN(ZPrimVtx, zPrimVtx, float); // Pri
DECLARE_SOA_COLUMN(XDecVtx, xDecVtx, float); // Decay vertex of the candidate (x direction)
DECLARE_SOA_COLUMN(YDecVtx, yDecVtx, float); // Decay vertex of the candidate (y direction)
DECLARE_SOA_COLUMN(ZDecVtx, zDecVtx, float); // Decay vertex of the candidate (z direction)
DECLARE_SOA_COLUMN(XMoth, xMoth, float); // X of the mother track at the radii of ITS layer which has the outermost update
DECLARE_SOA_COLUMN(YMoth, yMoth, float); // Y of the mother track at the radii of ITS layer which has the outermost update
DECLARE_SOA_COLUMN(ZMoth, zMoth, float); // Z of the mother track at the radii of ITS layer which has the outermost update
DECLARE_SOA_COLUMN(PxMoth, pxMoth, float); //! Px of the mother track at the decay vertex
DECLARE_SOA_COLUMN(PyMoth, pyMoth, float); //! Py of the mother track at the decay vertex
DECLARE_SOA_COLUMN(PzMoth, pzMoth, float); //! Pz of the mother track at the decay vertex
Expand All @@ -56,6 +59,9 @@ DECLARE_SOA_COLUMN(TrueZDecVtx, trueZDecVtx, float); // true z decay
DECLARE_SOA_COLUMN(GenPxMoth, genPxMoth, float); // Generated px of the mother track
DECLARE_SOA_COLUMN(GenPyMoth, genPyMoth, float); // Generated py of the mother track
DECLARE_SOA_COLUMN(GenPzMoth, genPzMoth, float); // Generated pz of the mother track
DECLARE_SOA_COLUMN(TruePxMoth, truePxMoth, float); // true px of the mother track at the decay vertex
DECLARE_SOA_COLUMN(TruePyMoth, truePyMoth, float); // true py of the mother track at the decay vertex
DECLARE_SOA_COLUMN(TruePzMoth, truePzMoth, float); // true pz of the mother track at the decay vertex
DECLARE_SOA_COLUMN(GenPxAlpha, genPxAlpha, float); // true px of the daughter alpha track
DECLARE_SOA_COLUMN(GenPyAlpha, genPyAlpha, float); // true py of the daughter alpha track
DECLARE_SOA_COLUMN(GenPzAlpha, genPzAlpha, float); // true pz of the daughter alpha track
Expand All @@ -68,7 +74,9 @@ DECLARE_SOA_TABLE(He4S2BCands, "AOD", "HE4S2BCANDS",
o2::soa::Index<>,
he4scand::XPrimVtx, he4scand::YPrimVtx, he4scand::ZPrimVtx,
he4scand::XDecVtx, he4scand::YDecVtx, he4scand::ZDecVtx,
he4scand::IsMatter, he4scand::PxMoth, he4scand::PyMoth, he4scand::PzMoth,
he4scand::IsMatter,
he4scand::XMoth, he4scand::YMoth, he4scand::ZMoth,
he4scand::PxMoth, he4scand::PyMoth, he4scand::PzMoth,
he4scand::PxAlpha, he4scand::PyAlpha, he4scand::PzAlpha,
he4scand::DcaMothPv, he4scand::DcaAlphaPv, he4scand::DcaKinkTopo,
he4scand::ItsChi2Moth, he4scand::ItsClusterSizesMoth, he4scand::ItsClusterSizesAlpha,
Expand All @@ -78,14 +86,17 @@ DECLARE_SOA_TABLE(MCHe4S2BCands, "AOD", "MCHE4S2BCANDS",
o2::soa::Index<>,
he4scand::XPrimVtx, he4scand::YPrimVtx, he4scand::ZPrimVtx,
he4scand::XDecVtx, he4scand::YDecVtx, he4scand::ZDecVtx,
he4scand::IsMatter, he4scand::PxMoth, he4scand::PyMoth, he4scand::PzMoth,
he4scand::IsMatter,
he4scand::XMoth, he4scand::YMoth, he4scand::ZMoth,
he4scand::PxMoth, he4scand::PyMoth, he4scand::PzMoth,
he4scand::PxAlpha, he4scand::PyAlpha, he4scand::PzAlpha,
he4scand::DcaMothPv, he4scand::DcaAlphaPv, he4scand::DcaKinkTopo,
he4scand::ItsChi2Moth, he4scand::ItsClusterSizesMoth, he4scand::ItsClusterSizesAlpha,
he4scand::NSigmaTPCAlpha, he4scand::NSigmaITSAlpha,
he4scand::IsSignal, he4scand::IsSignalReco, he4scand::IsCollReco, he4scand::IsSurvEvSelection,
he4scand::TrueXDecVtx, he4scand::TrueYDecVtx, he4scand::TrueZDecVtx,
he4scand::GenPxMoth, he4scand::GenPyMoth, he4scand::GenPzMoth,
he4scand::TruePxMoth, he4scand::TruePyMoth, he4scand::TruePzMoth,
he4scand::GenPxAlpha, he4scand::GenPyAlpha, he4scand::GenPzAlpha,
he4scand::IsMothReco, he4scand::RecoPtMoth, he4scand::RecoPzMoth);

Expand Down
71 changes: 54 additions & 17 deletions PWGLF/TableProducer/Nuspex/hyperhelium4sigmaRecoTask.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -68,7 +68,7 @@ enum DaughterType {

namespace
{
// constexpr std::array<float, 7> LayerRadii{2.33959f, 3.14076f, 3.91924f, 19.6213f, 24.5597f, 34.388f, 39.3329f};
constexpr std::array<float, 7> LayerRadii{2.33959f, 3.14076f, 3.91924f, 19.6213f, 24.5597f, 34.388f, 39.3329f};
constexpr int kITSLayers = 7;
constexpr int kITSInnerBarrelLayers = 3;
// constexpr int kITSOuterBarrelLayers = 4;
Expand Down Expand Up @@ -231,6 +231,7 @@ struct Hyphe4sCandidate {

std::array<float, 3> primVtx = {0.0f, 0.0f, 0.0f};
std::array<float, 3> decVtx = {0.0f, 0.0f, 0.0f};
std::array<float, 3> lastPosMoth = {0.0f, 0.0f, 0.0f}; // last position of mother track at the radii of ITS layer which has the outermost update
std::array<float, 3> momMoth = {0.0f, 0.0f, 0.0f};
std::array<float, 3> momDaug = {0.0f, 0.0f, 0.0f};

Expand All @@ -251,7 +252,8 @@ struct Hyphe4sCandidate {
bool isSurvEvSelection = false;

std::array<float, 3> trueDecVtx = {0.0f, 0.0f, 0.0f};
std::array<float, 3> gMomMoth = {0.0f, 0.0f, 0.0f};
std::array<float, 3> gMomMoth = {0.0f, 0.0f, 0.0f}; // generated mother momentum
std::array<float, 3> trueMomMoth = {0.0f, 0.0f, 0.0f}; // true mother momentum at decay vertex
std::array<float, 3> gMomDau = {0.0f, 0.0f, 0.0f};

bool isMothReco = false;
Expand Down Expand Up @@ -322,8 +324,8 @@ struct Hyperhelium4sigmaRecoTask {
registry.add<TH2>("h2TrueMotherDiffPyVsRecSVR", ";Reconstruced SV R (cm);#Delta #it{p}_{T} (GeV/#it{c});", HistType::kTH2F, {radiusAxis, diffPyAxis});
registry.add<TH2>("h2TrueMotherDiffPtVsRecSVR", ";Reconstruced SV R (cm);#Delta #it{p}_{T} (GeV/#it{c});", HistType::kTH2F, {radiusAxis, diffPtAxis});
registry.add<TH2>("h2TrueMotherDiffPzVsRecSVR", ";Reconstruced SV R (cm);#Delta #it{p}_{z} (GeV/#it{c});", HistType::kTH2F, {radiusAxis, diffPzAxis});
registry.add<TH2>("h2TrueMotherDiffTglVsRecSVR", ";Reconstruced SV R (cm);#Delta tan#lambda;", HistType::kTH2F, {radiusAxis, {200, -1.f, 1.f}});
registry.add<TH2>("h2TrueMotherDiffEtaVsRecSVR", ";Reconstruced SV R (cm);#Delta #eta;", HistType::kTH2F, {radiusAxis, {200, -1.f, 1.f}});
registry.add<TH2>("h2TrueMotherDiffTglVsRecSVR", ";Reconstruced SV R (cm);#Delta tan#lambda;", HistType::kTH2F, {radiusAxis, {200, -0.1, 0.1}});
registry.add<TH2>("h2TrueMotherDiffEtaVsRecSVR", ";Reconstruced SV R (cm);#Delta #eta;", HistType::kTH2F, {radiusAxis, {200, -0.1, 0.1}});
registry.add<TH1>("hDiffDauPx", ";#Delta p_{x} (GeV/#it{c}); ", HistType::kTH1D, {{200, -10, 10}});
registry.add<TH1>("hDiffDauPy", ";#Delta p_{y} (GeV/#it{c}); ", HistType::kTH1D, {{200, -10, 10}});
registry.add<TH1>("hDiffDauPz", ";#Delta p_{z} (GeV/#it{c}); ", HistType::kTH1D, {{200, -10, 10}});
Expand Down Expand Up @@ -394,6 +396,22 @@ struct Hyperhelium4sigmaRecoTask {
hyphe4sCand.itsClusterSizeDau = trackDau.itsClusterSizes();
hyphe4sCand.nSigmaTPCDau = trackDau.tpcNSigmaAl();
hyphe4sCand.nSigmaITSDau = itsResponse.nSigmaITS<o2::track::PID::Alpha>(trackDau);

int lastLayerMoth = 0;
for (int i = 6; i >= 0; i--) {
if (trackMoth.itsClusterMap() & (1 << i)) {
lastLayerMoth = i;
break;
}
}
auto trackparMother = getTrackParCov(trackMoth);
o2::base::Propagator::Instance()->PropagateToXBxByBz(trackparMother, LayerRadii[lastLayerMoth]);
std::array<float, 9> vecLab{0.f};
if (trackparMother.getPosDirGlo(vecLab)) {
hyphe4sCand.lastPosMoth[0] = vecLab[0];
hyphe4sCand.lastPosMoth[1] = vecLab[1];
hyphe4sCand.lastPosMoth[2] = vecLab[2];
}
}

template <typename TTrack>
Expand All @@ -407,20 +425,23 @@ struct Hyperhelium4sigmaRecoTask {
}

template <typename TMCParticle>
void fillCandidateMCInfo(Hyphe4sCandidate& hyphe4sCand, TMCParticle const& mcMothTrack, TMCParticle const& mcDauTrack)
void fillCandidateMCInfo(Hyphe4sCandidate& hyphe4sCand, TMCParticle const& mcMothTrack, TMCParticle const& mcDauTrack, TMCParticle const& mcNeutDauTrack)
{
hyphe4sCand.trueDecVtx[0] = mcDauTrack.vx();
hyphe4sCand.trueDecVtx[1] = mcDauTrack.vy();
hyphe4sCand.trueDecVtx[2] = mcDauTrack.vz();
hyphe4sCand.gMomMoth[0] = mcMothTrack.px();
hyphe4sCand.gMomMoth[1] = mcMothTrack.py();
hyphe4sCand.gMomMoth[2] = mcMothTrack.pz();
hyphe4sCand.trueMomMoth[0] = mcDauTrack.px() + mcNeutDauTrack.px();
hyphe4sCand.trueMomMoth[1] = mcDauTrack.py() + mcNeutDauTrack.py();
hyphe4sCand.trueMomMoth[2] = mcDauTrack.pz() + mcNeutDauTrack.pz();
hyphe4sCand.gMomDau[0] = mcDauTrack.px();
hyphe4sCand.gMomDau[1] = mcDauTrack.py();
hyphe4sCand.gMomDau[2] = mcDauTrack.pz();
}

void processData(CollisionsFull const& collisions, aod::KinkCands const& KinkCands, FullTracksExtIU const&)
void processData(CollisionsFull const& collisions, aod::KinkCands const& KinkCands, FullTracksExtIU const&, aod::BCs const&)
{
for (const auto& collision : collisions) {
registry.fill(HIST("hEventCounter"), 0);
Expand Down Expand Up @@ -448,14 +469,18 @@ struct Hyperhelium4sigmaRecoTask {
registry.fill(HIST("h2MassHyperhelium4sigmaPt"), kinkCand.mothSign() * kinkCand.ptMoth(), invMass);
registry.fill(HIST("h2NSigmaAlPt"), kinkCand.mothSign() * kinkCand.ptDaug(), dauTrack.tpcNSigmaAl());

auto bc = collision.bc_as<aod::BCs>();
initCCDB(bc);
auto motherTrack = kinkCand.trackMoth_as<FullTracksExtIU>();
Hyphe4sCandidate hyphe4sCand;
fillCandidate(hyphe4sCand, collision, kinkCand, motherTrack, dauTrack);

outputDataTable(
hyphe4sCand.primVtx[0], hyphe4sCand.primVtx[1], hyphe4sCand.primVtx[2],
hyphe4sCand.decVtx[0], hyphe4sCand.decVtx[1], hyphe4sCand.decVtx[2],
hyphe4sCand.isMatter, hyphe4sCand.momMoth[0], hyphe4sCand.momMoth[1], hyphe4sCand.momMoth[2],
hyphe4sCand.isMatter,
hyphe4sCand.lastPosMoth[0], hyphe4sCand.lastPosMoth[1], hyphe4sCand.lastPosMoth[2],
hyphe4sCand.momMoth[0], hyphe4sCand.momMoth[1], hyphe4sCand.momMoth[2],
hyphe4sCand.momDaug[0], hyphe4sCand.momDaug[1], hyphe4sCand.momDaug[2],
hyphe4sCand.dcaXYMothPv, hyphe4sCand.dcaXYDauPv, hyphe4sCand.dcaKinkTopo,
hyphe4sCand.chi2ITSMoth, hyphe4sCand.itsClusterSizeMoth, hyphe4sCand.itsClusterSizeDau,
Expand Down Expand Up @@ -524,13 +549,16 @@ struct Hyperhelium4sigmaRecoTask {
registry.fill(HIST("h2MassHyperhelium4sigmaPt"), kinkCand.mothSign() * kinkCand.ptMoth(), invMass);
registry.fill(HIST("h2NSigmaAlPt"), kinkCand.mothSign() * kinkCand.ptDaug(), dauTrack.tpcNSigmaAl());

auto bc = collision.bc_as<aod::BCs>();
initCCDB(bc);
Hyphe4sCandidate hyphe4sCand;
fillCandidate(hyphe4sCand, collision, kinkCand, motherTrack, dauTrack);

// qa for true signal
if (isTrueSignal) {
auto mcMotherTrack = motherTrack.mcParticle_as<aod::McParticles>();
auto mcDauTrack = dauTrack.mcParticle_as<aod::McParticles>();
auto mcNeutTrack = particlesMC.rawIteratorAt(dauIDList[2]);
float posDecVtx[3] = {kinkCand.xDecVtx() + collision.posX(), kinkCand.yDecVtx() + collision.posY(), kinkCand.zDecVtx() + collision.posZ()};
float recSVR = std::sqrt(posDecVtx[0] * posDecVtx[0] + posDecVtx[1] * posDecVtx[1]);
registry.fill(HIST("hDiffSVx"), posDecVtx[0] - mcDauTrack.vx());
Expand All @@ -551,11 +579,9 @@ struct Hyperhelium4sigmaRecoTask {
hyphe4sCand.isSignalReco = true;
hyphe4sCand.isCollReco = true;
hyphe4sCand.isSurvEvSelection = true;
fillCandidateMCInfo(hyphe4sCand, mcMotherTrack, mcDauTrack);
fillCandidateMCInfo(hyphe4sCand, mcMotherTrack, mcDauTrack, mcNeutTrack);
mcHe4sIndices.push_back(mcMotherTrack.globalIndex());

auto bc = collision.bc_as<aod::BCs>();
initCCDB(bc);
std::array<float, 2> dcaInfo;
auto mcMotherTrackPar = getTrackParFromMC(mcMotherTrack);
o2::base::Propagator::Instance()->propagateToDCABxByBz({posDecVtx[0], posDecVtx[1], posDecVtx[2]}, mcMotherTrackPar, 2.f, matCorr, &dcaInfo);
Expand Down Expand Up @@ -591,14 +617,17 @@ struct Hyperhelium4sigmaRecoTask {
outputMCTable(
hyphe4sCand.primVtx[0], hyphe4sCand.primVtx[1], hyphe4sCand.primVtx[2],
hyphe4sCand.decVtx[0], hyphe4sCand.decVtx[1], hyphe4sCand.decVtx[2],
hyphe4sCand.isMatter, hyphe4sCand.momMoth[0], hyphe4sCand.momMoth[1], hyphe4sCand.momMoth[2],
hyphe4sCand.isMatter,
hyphe4sCand.lastPosMoth[0], hyphe4sCand.lastPosMoth[1], hyphe4sCand.lastPosMoth[2],
hyphe4sCand.momMoth[0], hyphe4sCand.momMoth[1], hyphe4sCand.momMoth[2],
hyphe4sCand.momDaug[0], hyphe4sCand.momDaug[1], hyphe4sCand.momDaug[2],
hyphe4sCand.dcaXYMothPv, hyphe4sCand.dcaXYDauPv, hyphe4sCand.dcaKinkTopo,
hyphe4sCand.chi2ITSMoth, hyphe4sCand.itsClusterSizeMoth, hyphe4sCand.itsClusterSizeDau,
hyphe4sCand.nSigmaTPCDau, hyphe4sCand.nSigmaITSDau,
hyphe4sCand.isSignal, hyphe4sCand.isSignalReco, hyphe4sCand.isCollReco, hyphe4sCand.isSurvEvSelection,
hyphe4sCand.trueDecVtx[0], hyphe4sCand.trueDecVtx[1], hyphe4sCand.trueDecVtx[2],
hyphe4sCand.gMomMoth[0], hyphe4sCand.gMomMoth[1], hyphe4sCand.gMomMoth[2],
hyphe4sCand.trueMomMoth[0], hyphe4sCand.trueMomMoth[1], hyphe4sCand.trueMomMoth[2],
hyphe4sCand.gMomDau[0], hyphe4sCand.gMomDau[1], hyphe4sCand.gMomDau[2],
hyphe4sCand.isMothReco, hyphe4sCand.ptMoth, hyphe4sCand.pzMoth);
}
Expand All @@ -615,7 +644,8 @@ struct Hyperhelium4sigmaRecoTask {

Hyphe4sCandidate hyphe4sCand;
auto mcDauTrack = particlesMC.rawIteratorAt(dauIDList[0]);
fillCandidateMCInfo(hyphe4sCand, mcparticle, mcDauTrack);
auto mcNeutTrack = particlesMC.rawIteratorAt(dauIDList[2]);
fillCandidateMCInfo(hyphe4sCand, mcparticle, mcDauTrack, mcNeutTrack);

if (mcPartIndices[mcparticle.globalIndex()] != -1) {
auto mothTrack = tracks.rawIteratorAt(mcPartIndices[mcparticle.globalIndex()]);
Expand All @@ -625,14 +655,17 @@ struct Hyperhelium4sigmaRecoTask {
outputMCTable(
-1, -1, -1,
-1, -1, -1,
-1, -1, -1, -1,
-1,
-1, -1, -1,
-1, -1, -1,
-1, -1, -1,
-1, -1, -1,
-1, -1, -1,
-1, -1,
true, false, isReconstructedMCCollisions[mcparticle.mcCollisionId()], isSelectedMCCollisions[mcparticle.mcCollisionId()],
hyphe4sCand.trueDecVtx[0], hyphe4sCand.trueDecVtx[1], hyphe4sCand.trueDecVtx[2],
hyphe4sCand.gMomMoth[0], hyphe4sCand.gMomMoth[1], hyphe4sCand.gMomMoth[2],
hyphe4sCand.trueMomMoth[0], hyphe4sCand.trueMomMoth[1], hyphe4sCand.trueMomMoth[2],
hyphe4sCand.gMomDau[0], hyphe4sCand.gMomDau[1], hyphe4sCand.gMomDau[2],
hyphe4sCand.isMothReco, hyphe4sCand.ptMoth, hyphe4sCand.pzMoth);
}
Expand Down Expand Up @@ -731,6 +764,7 @@ struct Hyperhelium4sigmaQa {
hDauTPCNSigma[kDauTriton] = recoQAHist.add<TH2>("hDauTritonTPCNSigma", "", HistType::kTH2F, {rigidityAxis, nsigmaAxis});
hDauTPCNSigma[kDauProton] = recoQAHist.add<TH2>("hDauProtonTPCNSigma", "", HistType::kTH2F, {rigidityAxis, nsigmaAxis});
hDauTPCNSigma[kDauChargedPion] = recoQAHist.add<TH2>("hDauPionTPCNSigma", "", HistType::kTH2F, {rigidityAxis, nsigmaAxis});
recoQAHist.add<TH2>("hDauAlphaITSNSigmaCheck", "", HistType::kTH2F, {rigidityAxis, itsnsigmaAxis});

hRecoMotherCounter = recoQAHist.add<TH1>("hRecoMotherCounter", "", HistType::kTH1F, {{9, 0.f, 9.f}});
hRecoDauAlphaCounter = recoQAHist.add<TH1>("hRecoDauAlphaCounter", "", HistType::kTH1F, {{9, 0.f, 9.f}});
Expand All @@ -750,8 +784,8 @@ struct Hyperhelium4sigmaQa {
recoQAHist.add<TH1>("hMotherITSCls", "", HistType::kTH1F, {{8, 0.f, 8.f}});
recoQAHist.add<TH1>("hDauAlphaIsPVContributer", "", HistType::kTH1F, {{2, 0.f, 2.f}});
recoQAHist.add<TH1>("hDauAlphaITSCls", "", HistType::kTH1F, {{8, 0.f, 8.f}});
recoQAHist.add<TH2>("hDauAlphaPVsITSNSigma", "", HistType::kTH2F, {rigidityAxis, itsnsigmaAxis});
recoQAHist.add<TH2>("h2BCandDauAlphaPVsITSNSigma", "", HistType::kTH2F, {rigidityAxis, nsigmaAxis});
recoQAHist.add<TH2>("hDauAlphaITSNSigma", "", HistType::kTH2F, {rigidityAxis, itsnsigmaAxis});
recoQAHist.add<TH2>("hReco2BDauAlphaPVsITSNSigma", "", HistType::kTH2F, {rigidityAxis, itsnsigmaAxis});
recoQAHist.add<TH1>("hReco2BCandidateCount", "", HistType::kTH1F, {{4, 0.f, 4.f}});
}
}
Expand Down Expand Up @@ -950,6 +984,9 @@ struct Hyperhelium4sigmaQa {
if (track.hasTPC()) {
hDauTPCNSigma[type]->Fill(track.p() * track.sign(), tpcNSigma);
}
if (type == kDauAlpha && track.itsNCls() > kITSLayers - 2) {
recoQAHist.fill(HIST("hDauAlphaITSNSigmaCheck"), track.p() * track.sign(), itsResponse.nSigmaITS<o2::track::PID::Alpha>(track));
}
}
}
}
Expand Down Expand Up @@ -1000,7 +1037,7 @@ struct Hyperhelium4sigmaQa {

float itsNSigma = itsResponse.nSigmaITS<o2::track::PID::Alpha>(daughterTrack);
if (daughterTrack.hasITS()) {
recoQAHist.fill(HIST("hDauAlphaPVsITSNSigma"), daughterTrack.sign() * daughterTrack.p(), itsNSigma);
recoQAHist.fill(HIST("hDauAlphaITSNSigma"), daughterTrack.sign() * daughterTrack.p(), itsNSigma);
recoQAHist.fill(HIST("hDauAlphaITSCls"), daughterTrack.itsNCls());
}

Expand All @@ -1013,7 +1050,7 @@ struct Hyperhelium4sigmaQa {

if (isMoth && isDaug) {
recoQAHist.fill(HIST("hReco2BCandidateCount"), 3.5);
recoQAHist.fill(HIST("h2BCandDauAlphaPVsITSNSigma"), daughterTrack.sign() * daughterTrack.p(), itsNSigma);
recoQAHist.fill(HIST("hReco2BDauAlphaPVsITSNSigma"), daughterTrack.sign() * daughterTrack.p(), itsNSigma);
}
}
}
Expand Down
Loading