Skip to content

[PWGDQ] remove the kPhyiscalPrimay restriction for muon from various sources in MCSignals #11316

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
iarsene merged 3 commits into from
May 26, 2025
Merged

[PWGDQ] remove the kPhyiscalPrimay restriction for muon from various sources in MCSignals #11316

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
9061aa5
Remove the kPhysicalPrimary restriction for muFromXX signals
MaolinZH May 22, 2025
35f3034
Add two cuts for type0 and type3 selections combined with muonQuality…
MaolinZH May 22, 2025
b0e35b0
Change the tempory histograms subgroup tpcpid_nsigel into tpcpidvspt
MaolinZH May 22, 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
11 changes: 11 additions & 0 deletions PWGDQ/Core/CutsLibrary.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -3188,6 +3188,17 @@ AnalysisCompositeCut* o2::aod::dqcuts::GetCompositeCut(const char* cutName)
return cut;
}

if (!nameStr.compare("muonQualityCutsMUONStandalone")) {
cut->AddCut(GetAnalysisCut("matchedMchMid"));
cut->AddCut(GetAnalysisCut("muonQualityCuts"));
return cut;
}

if (!nameStr.compare("muonQualityCutsGlobal")) {
cut->AddCut(GetAnalysisCut("matchedGlobal"));
cut->AddCut(GetAnalysisCut("muonQualityCuts"));
return cut;
}
// -----------------------------------------------------------
// Pair cuts
if (!nameStr.compare("pairNoCut")) {
Expand Down
9 changes: 5 additions & 4 deletions PWGDQ/Core/HistogramsLibrary.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -565,10 +565,11 @@ void o2::aod::dqhistograms::DefineHistograms(HistogramManager* hm, const char* h
hm->AddHistogram(histClass, "TPCnSigKa_Corr_pIN", "TPC n-#sigma(K) Corr. vs pIN", false, 100, 0.0, 10.0, VarManager::kPin, 100, -5.0, 5.0, VarManager::kTPCnSigmaKa_Corr);
hm->AddHistogram(histClass, "TPCnSigPr_Corr_pIN", "TPC n-#sigma(p) Corr. vs pIN", false, 100, 0.0, 10.0, VarManager::kPin, 100, -5.0, 5.0, VarManager::kTPCnSigmaPr_Corr);
}
if (subGroupStr.Contains("tpcpid_nsigel")) {
hm->AddHistogram(histClass, "TPCnSigEl_TPCnSigPi", "TPC n-#sigma(e) . vs TPC n-#sigma(#pi)", false, 100, -5.0, 5.0, VarManager::kTPCnSigmaEl, 100, -5.0, 5.0, VarManager::kTPCnSigmaPi);
hm->AddHistogram(histClass, "TPCnSigEl_TPCnSigKa", "TPC n-#sigma(e) . vs TPC n-#sigma(K)", false, 100, -5.0, 5.0, VarManager::kTPCnSigmaEl, 100, -5.0, 5.0, VarManager::kTPCnSigmaKa);
hm->AddHistogram(histClass, "TPCnSigEl_TPCnSigPr", "TPC n-#sigma(e) . vs TPC n-#sigma(p)", false, 100, -5.0, 5.0, VarManager::kTPCnSigmaEl, 100, -5.0, 5.0, VarManager::kTPCnSigmaPr);
if (subGroupStr.Contains("tpcpidvspt")) {
hm->AddHistogram(histClass, "TPCnSigEl_Pt", "TPC n-#sigma(e). vs Pt", false, 200, 0.0, 20.0, VarManager::kPt, 100, -5.0, 5.0, VarManager::kTPCnSigmaEl);
hm->AddHistogram(histClass, "TPCnSigPi_Pt", "TPC n-#sigma(#pi). vs Pt", false, 200, 0.0, 20.0, VarManager::kPt, 100, -5.0, 5.0, VarManager::kTPCnSigmaPi);
hm->AddHistogram(histClass, "TPCnSigKa_Pt", "TPC n-#sigma(K). vs Pt", false, 200, 0.0, 20.0, VarManager::kPt, 100, -5.0, 5.0, VarManager::kTPCnSigmaKa);
hm->AddHistogram(histClass, "TPCnSigPr_Pt", "TPC n-#sigma(p). vs Pt", false, 200, 0.0, 20.0, VarManager::kPt, 100, -5.0, 5.0, VarManager::kTPCnSigmaPr);
Comment on lines -568 to +572

@iarsene iarsene May 22, 2025

Copy link
Copy Markdown
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Are these your changes or its because of not rebasing before pushing your changes ?

@MaolinZH MaolinZH May 23, 2025

Copy link
Copy Markdown
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Hi Ionut, this subgroup was added by me for some tempory tests. And since I assume no one else is using them I change them for some other tests.

}
}
}
Expand Down
12 changes: 1 addition & 11 deletions PWGDQ/Core/MCSignalLibrary.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -610,19 +610,16 @@ MCSignal* o2::aod::dqmcsignals::GetMCSignal(const char* name)
}
if (!nameStr.compare("eFromPromptLM")) {
MCProng prong(2, {11, 101}, {true, true}, {false, false}, {0, 0}, {0, 0}, {false, false}, false, {502, 402}, {true, true});
prong.SetSourceBit(0, MCProng::kPhysicalPrimary);
signal = new MCSignal(name, "Electrons from light mesons without B/D in decay history", {prong}, {-1});
return signal;
}
if (!nameStr.compare("eFromHbtoLM")) {
MCProng prong(2, {11, 101}, {true, true}, {false, false}, {0, 0}, {0, 0}, {false, false}, false, {502}, {false});
prong.SetSourceBit(0, MCProng::kPhysicalPrimary);
signal = new MCSignal(name, "Electrons from light mesons with B hadron in decay history", {prong}, {-1});
return signal;
}
if (!nameStr.compare("eFromHctoLM")) {
MCProng prong(2, {11, 101, 402}, {true, true, true}, {false, false, false}, {0, 0, 0}, {0, 0, 0}, {false, false, false}, false, {502}, {true});
prong.SetSourceBit(0, MCProng::kPhysicalPrimary);
signal = new MCSignal(name, "Electrons from light mesons from D hadron decays and no B in decay history", {prong}, {-1});
return signal;
}
Expand Down Expand Up @@ -655,44 +652,37 @@ MCSignal* o2::aod::dqmcsignals::GetMCSignal(const char* name)
}
if (!nameStr.compare("muFromHb")) {
MCProng prong(2, {13, 502}, {true, true}, {false, false}, {0, 0}, {0, 0}, {false, false});
prong.SetSourceBit(0, MCProng::kPhysicalPrimary);
signal = new MCSignal(name, "muons from b->mu", {prong}, {-1});
return signal;
}
if (!nameStr.compare("muFromPromptHc")) {
MCProng prong(2, {13, 402}, {true, true}, {false, false}, {0, 0}, {0, 0}, {false, false}, false, {502}, {true});
prong.SetSourceBit(0, MCProng::kPhysicalPrimary);
signal = new MCSignal(name, "muons from c->mu, without beauty in decay history", {prong}, {-1});
return signal;
}
if (!nameStr.compare("muFromHbtoHc")) {
MCProng prong(3, {13, 402, 502}, {true, true, true}, {false, false, false}, {0, 0, 0}, {0, 0, 0}, {false, false, false});
prong.SetSourceBit(0, MCProng::kPhysicalPrimary);
signal = new MCSignal(name, "muons from b->c->mu", {prong}, {-1});
return signal;
}
if (!nameStr.compare("secondaryMuon")) {
MCProng prong(1);
prong.SetPDGcode(0, 13, true);
MCProng prong(1, {13}, {true}, {false}, {0}, {0}, {false});
prong.SetSourceBit(0, MCProng::kProducedInTransport);
signal = new MCSignal(name, "muons produced during transport in detector", {prong}, {-1});
return signal;
}
if (!nameStr.compare("muFromPromptLM")) {
MCProng prong(2, {13, 101}, {true, true}, {false, false}, {0, 0}, {0, 0}, {false, false}, false, {502, 402}, {true, true});
prong.SetSourceBit(0, MCProng::kPhysicalPrimary);
signal = new MCSignal(name, "muons from light mesons without B/D in decay history", {prong}, {-1});
return signal;
}
if (!nameStr.compare("muFromHbtoLM")) {
MCProng prong(2, {13, 101}, {true, true}, {false, false}, {0, 0}, {0, 0}, {false, false}, false, {502}, {false});
prong.SetSourceBit(0, MCProng::kPhysicalPrimary);
signal = new MCSignal(name, "muons from light mesons with B hadron in decay history", {prong}, {-1});
return signal;
}
if (!nameStr.compare("muFromHctoLM")) {
MCProng prong(2, {13, 101, 402}, {true, true, true}, {false, false, false}, {0, 0, 0}, {0, 0, 0}, {false, false, false}, false, {502}, {true});
prong.SetSourceBit(0, MCProng::kPhysicalPrimary);
signal = new MCSignal(name, "muons from light mesons from D hadron decays and no B in decay history", {prong}, {-1});
return signal;
}
Expand Down