Skip to content

Navigation Menu

Sign in
Appearance settings

Search code, repositories, users, issues, pull requests...

Provide feedback

We read every piece of feedback, and take your input very seriously.

Saved searches

Use saved searches to filter your results more quickly
Appearance settings

PWGHF: Little fix of the J/ψ → e+ e- code. #5265

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
jgrosseo merged 1 commit into from
Jan 24, 2021
Merged

PWGHF: Little fix of the J/ψ → e+ e- code. #5265

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
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
20 changes: 11 additions & 9 deletions 20 Analysis/Tasks/PWGHF/HFCandidateCreator2Prong.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -159,14 +159,15 @@ struct HFCandidateCreator2ProngMC {

// D0(bar) → π± K∓
//Printf("Checking D0(bar) → π± K∓");
if (RecoDecay::getMatchedMCRec(particlesMC, std::move(arrayDaughters), 421, array{+kPiPlus, -kKPlus}, true, &sign) > -1) {
if (RecoDecay::getMatchedMCRec(particlesMC, arrayDaughters, 421, array{+kPiPlus, -kKPlus}, true, &sign) > -1) {
result = sign * D0ToPiK;
}
// Jpsi → e+e-

// J/ψ → e+ e-
if (result == N2ProngDecays) {
//Printf("Checking Jpsi → e+e-");
if (RecoDecay::getMatchedMCRec(particlesMC, std::move(arrayDaughters), 443, array{+kElectron, -kElectron}, true, &sign) > -1) {
result = sign * JpsiToEE;
//Printf("Checking J/ψ → e+ e-");
if (RecoDecay::getMatchedMCRec(particlesMC, std::move(arrayDaughters), 443, array{+kElectron, -kElectron}, true) > -1) {
result = JpsiToEE;
}
}

Expand All @@ -183,11 +184,12 @@ struct HFCandidateCreator2ProngMC {
if (RecoDecay::isMatchedMCGen(particlesMC, particle, 421, array{+kPiPlus, -kKPlus}, true, &sign)) {
result = sign * D0ToPiK;
}
// Jpsi → e+e-

// J/ψ → e+ e-
if (result == N2ProngDecays) {
//Printf("Checking Jpsi → e+e-");
if (RecoDecay::isMatchedMCGen(particlesMC, particle, 443, array{+kElectron, -kElectron}, true, &sign)) {
result = sign * JpsiToEE;
//Printf("Checking J/ψ → e+ e-");
if (RecoDecay::isMatchedMCGen(particlesMC, particle, 443, array{+kElectron, -kElectron}, true)) {
result = JpsiToEE;
}
}

Expand Down
6 changes: 1 addition & 5 deletions 6 Analysis/Tasks/PWGHF/HFJpsiToEECandidateSelector.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -134,11 +134,9 @@ struct HFJpsiToEECandidateSelector {
return true;
}

/// Check if track is compatible with given TPC Nsigma cut for a given flavour hypothesis
/// Check if track is compatible with given TPC Nsigma cut for the electron hypothesis
/// \param track is the track
/// \param nPDG is the flavour hypothesis PDG number
/// \param nSigmaCut is the nsigma threshold to test against
/// \note nPDG=11 electron
/// \return true if track satisfies TPC PID hypothesis for given Nsigma cut
template <typename T>
bool selectionPIDTPC(const T& track, int nSigmaCut)
Expand All @@ -148,8 +146,6 @@ struct HFJpsiToEECandidateSelector {

/// PID selection on daughter track
/// \param track is the daughter track
/// \param nPDG is the PDG code of the flavour hypothesis
/// \note nPDG=11 electron
/// \return 1 if successful PID match, 0 if successful PID rejection, -1 if no PID info
template <typename T>
int selectionPID(const T& track)
Expand Down
Morty Proxy This is a proxified and sanitized view of the page, visit original site.