Skip to content

Yuanzhe dev #5

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
wang-yuanzhe merged 9 commits into from
Sep 4, 2023
Merged

Yuanzhe dev #5

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
9 commits
Select commit Hold shift + click to select a range
a2433d8
Add QA, Decay3Body datalink and Mc labels
wang-yuanzhe Aug 1, 2023
2f08658
Merge branch 'AliceO2Group:master' into yuanzhe_dev
wang-yuanzhe Aug 1, 2023
b871272
Please consider the following formatting changes
alibuild Aug 1, 2023
fcaea81
Merge pull request #4 from alibuild/alibot-cleanup-3235
wang-yuanzhe Aug 1, 2023
e9caad1
Add MC Dalitz plot
wang-yuanzhe Aug 1, 2023
0a46d00
Change getMassPDG to physics constant
wang-yuanzhe Aug 1, 2023
d7138c4
Change getMassPDG to physics constant
wang-yuanzhe Aug 1, 2023
9f94340
Merge branch 'master' of github.com:wang-yuanzhe/O2Physics into yuanz…
wang-yuanzhe Sep 4, 2023
ae92bf9
MC sigma Update
wang-yuanzhe Sep 4, 2023
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
8 changes: 4 additions & 4 deletions PWGLF/TableProducer/hypertriton3bodybuilder.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -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})"}}}},
Expand Down Expand Up @@ -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) {
Expand All @@ -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) {
Expand Down
16 changes: 8 additions & 8 deletions PWGLF/Tasks/hypertriton3bodyMCQA.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -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;

Expand Down Expand Up @@ -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}}}},
},
};

Expand All @@ -398,11 +398,11 @@ struct hypertriton3bodyLabelCheck {
registry.fill(HIST("hLabeledVtxCounter"), 1.5);
auto mcparticle = vtx.mcParticle_as<aod::McParticles>();
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());
}
}
}
Expand Down
2 changes: 1 addition & 1 deletion PWGLF/Tasks/hypertriton3bodyanalysis.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -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;

Expand Down