Skip to content

Please consider the following formatting changes to #14880 #32

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
SuJeong-Ji merged 1 commit into from
Feb 10, 2026
Merged

Please consider the following formatting changes to #14880 #32

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
101 changes: 44 additions & 57 deletions PWGLF/Tasks/Resonances/chk892pp.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -191,7 +191,7 @@ struct Chk892pp {

/// PID Selections, pion
struct : ConfigurableGroup {
Configurable<bool> cfgTPConly{"cfgTPConly", true, "Use only TPC for PID"}; // bool
Configurable<bool> cfgTPConly{"cfgTPConly", true, "Use only TPC for PID"}; // bool
Configurable<float> cfgMaxTPCnSigmaPion{"cfgMaxTPCnSigmaPion", 5.0, "TPC nSigma cut for Pion"}; // TPC
Configurable<float> cfgMaxTOFnSigmaPion{"cfgMaxTOFnSigmaPion", 5.0, "TOF nSigma cut for Pion"}; // TOF
Configurable<float> cfgNsigmaCutCombinedPion{"cfgNsigmaCutCombinedPion", -999, "Combined nSigma cut for Pion"}; // Combined
Expand Down Expand Up @@ -546,8 +546,8 @@ struct Chk892pp {
if (std::abs(track.dcaXY()) > TrackCuts.cfgMaxbDCArToPVcut)
return false;
}
if (TrackCuts.cfgpTdepDCAzCut) {
// Tuned on the LHC22f anchored MC LHC23d1d on primary pions. 7 Sigmas of the resolution
if (TrackCuts.cfgpTdepDCAzCut) {
// Tuned on the LHC22f anchored MC LHC23d1d on primary pions. 7 Sigmas of the resolution
if (std::abs(track.dcaZ()) > (0.004 + (0.013 / track.pt())))
return false;
} else {
Expand Down Expand Up @@ -914,7 +914,6 @@ struct Chk892pp {
return true;
} // matchRecoToTruthKstar


void effKstarProcessGen(MCTrueTrackCandidates const& mcparts)
{
for (const auto& part : mcparts) {
Expand All @@ -925,44 +924,37 @@ struct Chk892pp {
if (std::abs(part.y()) > KstarCuts.cfgKstarMaxRap)
continue;

const int pionWanted = (part.pdgCode() > 0) ? +kPiPlus : -kPiPlus;
bool hasRightPion = false;
bool hasK0sToPipi = false;

for (const auto& d1 : part.template daughters_as<MCTrueTrackCandidates>()) {
const int pdg1 = d1.pdgCode();
if (pdg1 == pionWanted) {
hasRightPion = true;
} else if (std::abs(pdg1) == kPDGK0) {
for (const auto& d2 : d1.template daughters_as<MCTrueTrackCandidates>()) {
if (std::abs(d2.pdgCode()) == kPDGK0s) {
bool seenPip = false, seenPim = false;
for (const auto& d3 : d2.template daughters_as<MCTrueTrackCandidates>()) {
if (d3.pdgCode() == +kPiPlus)
seenPip = true;
else if (d3.pdgCode() == -kPiPlus)
seenPim = true;
}
if (seenPip && seenPim) {
hasK0sToPipi = true;
break;
}
}
}
}
if (hasRightPion && hasK0sToPipi)
break;
}

const int pionWanted = (part.pdgCode() > 0) ? +kPiPlus : -kPiPlus;
bool hasRightPion = false;
bool hasK0sToPipi = false;

for (const auto& d1 : part.template daughters_as<MCTrueTrackCandidates>()) {
const int pdg1 = d1.pdgCode();
if (pdg1 == pionWanted)
{
hasRightPion = true;
} else if (std::abs(pdg1) == kPDGK0)
{
for (const auto& d2 : d1.template daughters_as<MCTrueTrackCandidates>())
{
if (std::abs(d2.pdgCode()) == kPDGK0s)
{
bool seenPip = false, seenPim = false;
for (const auto& d3 : d2.template daughters_as<MCTrueTrackCandidates>())
{
if (d3.pdgCode() == +kPiPlus)
seenPip = true;
else if (d3.pdgCode() == -kPiPlus)
seenPim = true;
}
if (seenPip && seenPim)
{
hasK0sToPipi = true;
break;
}
}
}
}
if (hasRightPion && hasK0sToPipi) break;
}

if (!(hasRightPion && hasK0sToPipi))
continue;

if (!(hasRightPion && hasK0sToPipi))
continue;

const auto mcid = part.mcCollisionId();
if (allowedMcIds.count(mcid) == 0)
Expand All @@ -973,15 +965,14 @@ struct Chk892pp {
continue;

const float lCentrality = iter->second;

histos.fill(HIST("EffKstar/genKstar"), part.pt(), lCentrality);

if (part.vt() == 0)
{
histos.fill(HIST("EffKstar/genKstar_pri"), part.pt(), lCentrality);
}
if (part.vt() == 0) {
histos.fill(HIST("EffKstar/genKstar_pri"), part.pt(), lCentrality);
}
}
} //effKstarProcessGen
} // effKstarProcessGen

template <typename V0RangeT, typename TrkRangeT>
void effKstarProcessReco(V0RangeT const& v0s, TrkRangeT const& tracks)
Expand Down Expand Up @@ -1068,11 +1059,9 @@ struct Chk892pp {
const float lCentrality = iter->second;

histos.fill(HIST("Correction/sigLoss_num"), part.pt(), lCentrality);
if (part.vt() == 0)
{
histos.fill(HIST("Correction/sigLoss_num_pri"), part.pt(), lCentrality);
}

if (part.vt() == 0) {
histos.fill(HIST("Correction/sigLoss_num_pri"), part.pt(), lCentrality);
}
}
} // fillSigLossNum

Expand All @@ -1097,11 +1086,9 @@ struct Chk892pp {
const float lCentrality = iter->second;

histos.fill(HIST("Correction/sigLoss_den"), part.pt(), lCentrality);
if (part.vt() == 0)
{
histos.fill(HIST("Correction/sigLoss_den_pri"), part.pt(), lCentrality);
}

if (part.vt() == 0) {
histos.fill(HIST("Correction/sigLoss_den_pri"), part.pt(), lCentrality);
}
}
} // fillSigLossDen

Expand Down