Skip to content

Commit ecdcd98

Browse files
scannitoscannito
authored
[PWGLF] Fix logic in particle id (#12035)
1 parent 0e9ab89 commit ecdcd98

File tree

1 file changed

+27
-57
lines changed

1 file changed

+27
-57
lines changed

‎PWGLF/Tasks/Strangeness/phik0shortanalysis.cxx‎

Lines changed: 27 additions & 57 deletions
Original file line numberDiff line numberDiff line change
@@ -2843,28 +2843,16 @@ struct Phik0shortanalysis {
28432843
continue;
28442844

28452845
// Phi selection
2846-
if (mcParticle.pdgCode() != o2::constants::physics::Pdg::kPhi)
2847-
continue;
2848-
if (mcParticle.pt() < minPhiPt)
2849-
continue;
2850-
2851-
mcPhiHist.fill(HIST("h3PhiMCGenRecoCheckNewProc"), genmultiplicity, mcParticle.pt(), mcParticle.y());
2846+
if (mcParticle.pdgCode() == o2::constants::physics::Pdg::kPhi && mcParticle.pt() >= minPhiPt)
2847+
mcPhiHist.fill(HIST("h3PhiMCGenRecoCheckNewProc"), genmultiplicity, mcParticle.pt(), mcParticle.y());
28522848

28532849
// K0S selection
2854-
if (mcParticle.pdgCode() != PDG_t::kK0Short)
2855-
continue;
2856-
if (!mcParticle.isPhysicalPrimary() || mcParticle.pt() < v0Configs.v0SettingMinPt)
2857-
continue;
2858-
2859-
mcK0SHist.fill(HIST("h3K0SMCGenRecoCheckNewProc"), genmultiplicity, mcParticle.pt(), mcParticle.y());
2850+
if (mcParticle.pdgCode() == PDG_t::kK0Short && mcParticle.isPhysicalPrimary() && mcParticle.pt() >= v0Configs.v0SettingMinPt)
2851+
mcK0SHist.fill(HIST("h3K0SMCGenRecoCheckNewProc"), genmultiplicity, mcParticle.pt(), mcParticle.y());
28602852

28612853
// Pion selection
2862-
if (std::abs(mcParticle.pdgCode()) != PDG_t::kPiPlus)
2863-
continue;
2864-
if (!mcParticle.isPhysicalPrimary() || mcParticle.pt() < trackConfigs.cMinPionPtcut)
2865-
continue;
2866-
2867-
mcPionHist.fill(HIST("h3PiMCGenRecoCheckNewProc"), genmultiplicity, mcParticle.pt(), mcParticle.y());
2854+
if (std::abs(mcParticle.pdgCode()) == PDG_t::kPiPlus && mcParticle.isPhysicalPrimary() && mcParticle.pt() >= trackConfigs.cMinPionPtcut)
2855+
mcPionHist.fill(HIST("h3PiMCGenRecoCheckNewProc"), genmultiplicity, mcParticle.pt(), mcParticle.y());
28682856
}
28692857
}
28702858

@@ -2891,25 +2879,16 @@ struct Phik0shortanalysis {
28912879
continue;
28922880

28932881
// Phi selection
2894-
if (mcParticle.pdgCode() != o2::constants::physics::Pdg::kPhi)
2895-
continue;
2896-
if (mcParticle.pt() < minPhiPt)
2897-
continue;
2898-
mcPhiHist.fill(HIST("h3PhiMCGenRecoNewProc"), genmultiplicity, mcParticle.pt(), mcParticle.y());
2882+
if (mcParticle.pdgCode() == o2::constants::physics::Pdg::kPhi && mcParticle.pt() >= minPhiPt)
2883+
mcPhiHist.fill(HIST("h3PhiMCGenRecoNewProc"), genmultiplicity, mcParticle.pt(), mcParticle.y());
28992884

29002885
// K0S selection
2901-
if (mcParticle.pdgCode() != PDG_t::kK0Short)
2902-
continue;
2903-
if (!mcParticle.isPhysicalPrimary() || mcParticle.pt() < v0Configs.v0SettingMinPt)
2904-
continue;
2905-
mcK0SHist.fill(HIST("h3K0SMCGenRecoNewProc"), genmultiplicity, mcParticle.pt(), mcParticle.y());
2886+
if (mcParticle.pdgCode() == PDG_t::kK0Short && mcParticle.isPhysicalPrimary() && mcParticle.pt() >= v0Configs.v0SettingMinPt)
2887+
mcK0SHist.fill(HIST("h3K0SMCGenRecoNewProc"), genmultiplicity, mcParticle.pt(), mcParticle.y());
29062888

29072889
// Pion selection
2908-
if (std::abs(mcParticle.pdgCode()) != PDG_t::kPiPlus)
2909-
continue;
2910-
if (!mcParticle.isPhysicalPrimary() || mcParticle.pt() < trackConfigs.cMinPionPtcut)
2911-
continue;
2912-
mcPionHist.fill(HIST("h3PiMCGenRecoNewProc"), genmultiplicity, mcParticle.pt(), mcParticle.y());
2890+
if (std::abs(mcParticle.pdgCode()) == PDG_t::kPiPlus && mcParticle.isPhysicalPrimary() && mcParticle.pt() >= trackConfigs.cMinPionPtcut)
2891+
mcPionHist.fill(HIST("h3PiMCGenRecoNewProc"), genmultiplicity, mcParticle.pt(), mcParticle.y());
29132892
}
29142893

29152894
numberAssocColl++;
@@ -2929,34 +2908,25 @@ struct Phik0shortanalysis {
29292908
continue;
29302909

29312910
// Phi selection
2932-
if (mcParticle.pdgCode() != o2::constants::physics::Pdg::kPhi)
2933-
continue;
2934-
if (mcParticle.pt() < minPhiPt)
2935-
continue;
2936-
2937-
mcPhiHist.fill(HIST("h3PhiMCGenNewProc"), genmultiplicity, mcParticle.pt(), mcParticle.y());
2938-
if (numberAssocColl > 0)
2939-
mcPhiHist.fill(HIST("h3PhiMCGenAssocRecoNewProc"), genmultiplicity, mcParticle.pt(), mcParticle.y());
2911+
if (mcParticle.pdgCode() == o2::constants::physics::Pdg::kPhi && mcParticle.pt() >= minPhiPt) {
2912+
mcPhiHist.fill(HIST("h3PhiMCGenNewProc"), genmultiplicity, mcParticle.pt(), mcParticle.y());
2913+
if (numberAssocColl > 0)
2914+
mcPhiHist.fill(HIST("h3PhiMCGenAssocRecoNewProc"), genmultiplicity, mcParticle.pt(), mcParticle.y());
2915+
}
29402916

29412917
// K0S selection
2942-
if (mcParticle.pdgCode() != PDG_t::kK0Short)
2943-
continue;
2944-
if (!mcParticle.isPhysicalPrimary() || mcParticle.pt() < v0Configs.v0SettingMinPt)
2945-
continue;
2946-
2947-
mcK0SHist.fill(HIST("h3K0SMCGenNewProc"), genmultiplicity, mcParticle.pt(), mcParticle.y());
2948-
if (numberAssocColl > 0)
2949-
mcK0SHist.fill(HIST("h3K0SMCGenAssocRecoNewProc"), genmultiplicity, mcParticle.pt(), mcParticle.y());
2918+
if (mcParticle.pdgCode() == PDG_t::kK0Short && mcParticle.isPhysicalPrimary() && mcParticle.pt() >= v0Configs.v0SettingMinPt) {
2919+
mcK0SHist.fill(HIST("h3K0SMCGenNewProc"), genmultiplicity, mcParticle.pt(), mcParticle.y());
2920+
if (numberAssocColl > 0)
2921+
mcK0SHist.fill(HIST("h3K0SMCGenAssocRecoNewProc"), genmultiplicity, mcParticle.pt(), mcParticle.y());
2922+
}
29502923

29512924
// Pion selection
2952-
if (std::abs(mcParticle.pdgCode()) != PDG_t::kPiPlus)
2953-
continue;
2954-
if (!mcParticle.isPhysicalPrimary() || mcParticle.pt() < trackConfigs.cMinPionPtcut)
2955-
continue;
2956-
2957-
mcPionHist.fill(HIST("h3PiMCGenNewProc"), genmultiplicity, mcParticle.pt(), mcParticle.y());
2958-
if (numberAssocColl > 0)
2959-
mcPionHist.fill(HIST("h3PiMCGenAssocRecoNewProc"), genmultiplicity, mcParticle.pt(), mcParticle.y());
2925+
if (std::abs(mcParticle.pdgCode()) == PDG_t::kPiPlus && mcParticle.isPhysicalPrimary() && mcParticle.pt() >= trackConfigs.cMinPionPtcut) {
2926+
mcPionHist.fill(HIST("h3PiMCGenNewProc"), genmultiplicity, mcParticle.pt(), mcParticle.y());
2927+
if (numberAssocColl > 0)
2928+
mcPionHist.fill(HIST("h3PiMCGenAssocRecoNewProc"), genmultiplicity, mcParticle.pt(), mcParticle.y());
2929+
}
29602930
}
29612931
}
29622932

0 commit comments

Comments
 (0)