diff --git a/PWGLF/TableProducer/hypertriton3bodybuilder.cxx b/PWGLF/TableProducer/hypertriton3bodybuilder.cxx index 04c6370c477..c2c1651e222 100644 --- a/PWGLF/TableProducer/hypertriton3bodybuilder.cxx +++ b/PWGLF/TableProducer/hypertriton3bodybuilder.cxx @@ -305,8 +305,8 @@ struct hypertriton3bodyLabelBuilder { {"hLabelCounter", "hLabelCounter", {HistType::kTH1F, {{2, 0.0f, 2.0f}}}}, {"hHypertritonCounter", "hHypertritonCounter", {HistType::kTH1F, {{4, 0.0f, 4.0f}}}}, {"hPIDCounter", "hPIDCounter", {HistType::kTH1F, {{6, 0.0f, 6.0f}}}}, - {"hHypertriton", "hHypertriton", {HistType::kTH1F, {{100, 0.0f, 10.0f}}}}, - {"hAntiHypertriton", "hAntiHypertriton", {HistType::kTH1F, {{100, 0.0f, 10.0f}}}}, + {"hHypertritonMCPt", "hHypertritonMCPt", {HistType::kTH1F, {{100, 0.0f, 10.0f}}}}, + {"hAntiHypertritonMCPt", "hAntiHypertritonMCPt", {HistType::kTH1F, {{100, 0.0f, 10.0f}}}}, {"hHypertritonMass", "hHypertritonMass", {HistType::kTH1F, {{40, 2.95f, 3.05f}}}}, {"hAntiHypertritonMass", "hAntiHypertritonMass", {HistType::kTH1F, {{40, 2.95f, 3.05f}}}}, {"h3dTotalTrueHypertriton", "h3dTotalTrueHypertriton", {HistType::kTH3F, {{50, 0, 50, "ct(cm)"}, {200, 0.0f, 10.0f, "#it{p}_{T} (GeV/c)"}, {40, 2.95f, 3.05f, "Inv. Mass (GeV/c^{2})"}}}}, @@ -394,7 +394,7 @@ struct hypertriton3bodyLabelBuilder { if (lPDG == 1010010030 && lMCTrack0.pdgCode() == 2212 && lMCTrack1.pdgCode() == -211 && lMCTrack2.pdgCode() == 1000010020) { double hypertritonMCMass = RecoDecay::m(array{array{lMCTrack0.px(), lMCTrack0.py(), lMCTrack0.pz()}, array{lMCTrack1.px(), lMCTrack1.py(), lMCTrack1.pz()}, array{lMCTrack2.px(), lMCTrack2.py(), lMCTrack2.pz()}}, array{o2::constants::physics::MassProton, o2::constants::physics::MassPionCharged, o2::constants::physics::MassDeuteron}); registry.fill(HIST("hHypertritonCounter"), 0.5); - registry.fill(HIST("hHypertriton"), lPt); + registry.fill(HIST("hHypertritonMCPt"), lPt); registry.fill(HIST("hHypertritonMass"), hypertritonMCMass); registry.fill(HIST("h3dTotalTrueHypertriton"), MClifetime, lPt, hypertritonMCMass); if (TMath::Abs(lTrack0.tpcNSigmaPr()) > TpcPidNsigmaCut) { @@ -413,7 +413,7 @@ struct hypertriton3bodyLabelBuilder { if (lPDG == -1010010030 && lMCTrack0.pdgCode() == 211 && lMCTrack1.pdgCode() == -2212 && lMCTrack2.pdgCode() == -1000010020) { double antiHypertritonMCMass = RecoDecay::m(array{array{lMCTrack0.px(), lMCTrack0.py(), lMCTrack0.pz()}, array{lMCTrack1.px(), lMCTrack1.py(), lMCTrack1.pz()}, array{lMCTrack2.px(), lMCTrack2.py(), lMCTrack2.pz()}}, array{o2::constants::physics::MassPionCharged, o2::constants::physics::MassProton, o2::constants::physics::MassDeuteron}); registry.fill(HIST("hHypertritonCounter"), 2.5); - registry.fill(HIST("hAntiHypertriton"), lPt); + registry.fill(HIST("hAntiHypertritonMCPt"), lPt); registry.fill(HIST("hAntiHypertritonMass"), antiHypertritonMCMass); registry.fill(HIST("h3dTotalTrueHypertriton"), MClifetime, lPt, antiHypertritonMCMass); if (TMath::Abs(lTrack0.tpcNSigmaPi()) > TpcPidNsigmaCut) { diff --git a/PWGLF/Tasks/hypertriton3bodyMCQA.cxx b/PWGLF/Tasks/hypertriton3bodyMCQA.cxx index 05afc8b3cb5..e9e0c40557b 100644 --- a/PWGLF/Tasks/hypertriton3bodyMCQA.cxx +++ b/PWGLF/Tasks/hypertriton3bodyMCQA.cxx @@ -254,7 +254,7 @@ struct hypertriton3bodyAnalysisMc { } // 3sigma region for Dalitz plot - double mcsigma = 0.0018; // from local MC + double mcsigma = 0.0015; // from local MC double lowerlimit = o2::constants::physics::MassHyperTriton - 3 * mcsigma; double upperlimit = o2::constants::physics::MassHyperTriton + 3 * mcsigma; @@ -372,9 +372,9 @@ struct hypertriton3bodyLabelCheck { "registry", { {"hLabeledVtxCounter", "hLabeledVtxCounter", {HistType::kTH1F, {{2, 0.0f, 2.0f}}}}, - {"hMassTruesH3L", "hMassTruesH3L", {HistType::kTH1F, {{100, 2.95f, 3.05f}}}}, - {"hMassTruesH3LMatter", "hMassTruesH3LMatter", {HistType::kTH1F, {{100, 2.95f, 3.05f}}}}, - {"hMassTruesH3LAntiMatter", "hMassTruesH3LAntiMatter", {HistType::kTH1F, {{100, 2.95f, 3.05f}}}}, + {"hMassTrueH3L", "hMassTrueH3L", {HistType::kTH1F, {{100, 2.95f, 3.05f}}}}, + {"hMassTrueH3LMatter", "hMassTrueH3LMatter", {HistType::kTH1F, {{100, 2.95f, 3.05f}}}}, + {"hMassTrueH3LAntiMatter", "hMassTrueH3LAntiMatter", {HistType::kTH1F, {{100, 2.95f, 3.05f}}}}, }, }; @@ -398,11 +398,11 @@ struct hypertriton3bodyLabelCheck { registry.fill(HIST("hLabeledVtxCounter"), 1.5); auto mcparticle = vtx.mcParticle_as(); if (mcparticle.pdgCode() > 0) { - registry.fill(HIST("hMassTruesH3L"), vtx.mHypertriton()); - registry.fill(HIST("hMassTruesH3LMatter"), vtx.mHypertriton()); + registry.fill(HIST("hMassTrueH3L"), vtx.mHypertriton()); + registry.fill(HIST("hMassTrueH3LMatter"), vtx.mHypertriton()); } else { - registry.fill(HIST("hMassTruesH3L"), vtx.mAntiHypertriton()); - registry.fill(HIST("hMassTruesH3LAntiMatter"), vtx.mAntiHypertriton()); + registry.fill(HIST("hMassTrueH3L"), vtx.mAntiHypertriton()); + registry.fill(HIST("hMassTrueH3LAntiMatter"), vtx.mAntiHypertriton()); } } } diff --git a/PWGLF/Tasks/hypertriton3bodyanalysis.cxx b/PWGLF/Tasks/hypertriton3bodyanalysis.cxx index ca2b304b33d..472bd06ca46 100644 --- a/PWGLF/Tasks/hypertriton3bodyanalysis.cxx +++ b/PWGLF/Tasks/hypertriton3bodyanalysis.cxx @@ -199,7 +199,7 @@ struct hypertriton3bodyAnalysis { registry.fill(HIST("hSelectedCandidatesCounter"), 5.5); // 3sigma region for Dalitz plot - double mcsigma = 0.0018; // from local MC + double mcsigma = 0.0015; // from local MC double lowerlimit = o2::constants::physics::MassHyperTriton - 3 * mcsigma; double upperlimit = o2::constants::physics::MassHyperTriton + 3 * mcsigma;