Skip to content

Please consider the following formatting changes to #12622 #8

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
KotliarovAr merged 1 commit into from
Aug 18, 2025
Merged

Please consider the following formatting changes to #12622 #8

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
48 changes: 25 additions & 23 deletions PWGJE/Tasks/recoilJets.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -508,8 +508,9 @@ struct RecoilJets {
}

for (const auto& jet : jets) {
//skip jets which have a constituent with pT above specified cut
if(isJetWithHighPtConstituent(jet, tracks)) continue;
// skip jets which have a constituent with pT above specified cut
if (isJetWithHighPtConstituent(jet, tracks))
continue;

spectra.fill(HIST("hJetPtEtaPhiRhoArea"), jet.pt(), jet.eta(), jet.phi(),
collision.rho(), jet.area(), weight);
Expand Down Expand Up @@ -680,7 +681,7 @@ struct RecoilJets {
for (const auto& jetBase : jetsBase) {
bool bIsBaseJetRecoil =
get<1>(isRecoilJet(jetBase, phiTTSig)) && bIsThereTTSig;
dataForUnfolding(jetBase, jetsTag, bIsBaseJetRecoil, tracks, weight);
dataForUnfolding(jetBase, jetsTag, bIsBaseJetRecoil, tracks, weight);
}
}

Expand Down Expand Up @@ -911,30 +912,30 @@ struct RecoilJets {
template <typename Jet, typename Tracks>
bool isJetWithHighPtConstituent(Jet const& jet, Tracks const&)
{
bool bIsJetWithHighPtConstituent = false;
for (const auto& jetConstituent : jet.template tracks_as<Tracks>()) {
if (jetConstituent.pt() > maxJetConstituentPt) {
bIsJetWithHighPtConstituent = true;
break;
}
}
return bIsJetWithHighPtConstituent;
bool bIsJetWithHighPtConstituent = false;
for (const auto& jetConstituent : jet.template tracks_as<Tracks>()) {
if (jetConstituent.pt() > maxJetConstituentPt) {
bIsJetWithHighPtConstituent = true;
break;
}
}
return bIsJetWithHighPtConstituent;
}

template <typename PartJet, typename DetJet, typename TracksTable>
template <typename PartJet, typename DetJet, typename TracksTable>
void dataForUnfolding(PartJet const& partJet, DetJet const& detJets,
bool bIsBaseJetRecoil, TracksTable const& tracks, float weight = 1.)
bool bIsBaseJetRecoil, TracksTable const& tracks, float weight = 1.)
{

bool bIsThereMatchedJet = partJet.has_matchedJetGeo();
bool bIsThereMatchedJet = partJet.has_matchedJetGeo();

if (bIsThereMatchedJet) {
const auto& jetsMatched =
partJet.template matchedJetGeo_as<std::decay_t<DetJet>>();
if (bIsThereMatchedJet) {
const auto& jetsMatched =
partJet.template matchedJetGeo_as<std::decay_t<DetJet>>();

for (const auto& jetMatched : jetsMatched) {

//skip matches where detector level jets have a constituent with pT above specified cut
// skip matches where detector level jets have a constituent with pT above specified cut
bool skipMatchedDetJet = isJetWithHighPtConstituent(jetMatched, tracks);

if (skipMatchedDetJet) {
Expand All @@ -948,22 +949,22 @@ struct RecoilJets {
spectra.fill(HIST("hJetPt_DetLevel_vs_PartLevel"), jetMatched.pt(),
partJet.pt(), weight);
spectra.fill(HIST("hJetPt_resolution"),
(partJet.pt() - jetMatched.pt()) / partJet.pt(),
(partJet.pt() - jetMatched.pt()) / partJet.pt(),
partJet.pt(), weight);
spectra.fill(HIST("hJetPhi_resolution"),
partJet.phi() - jetMatched.phi(), partJet.pt(), weight);

if (bIsBaseJetRecoil) {
spectra.fill(HIST("hJetPt_DetLevel_vs_PartLevel_RecoilJets"),
jetMatched.pt(), partJet.pt(), weight);
spectra.fill(HIST("hJetPt_resolution_RecoilJets"),
(partJet.pt() - jetMatched.pt()) / partJet.pt(),
(partJet.pt() - jetMatched.pt()) / partJet.pt(),
partJet.pt(), weight);
spectra.fill(HIST("hJetPhi_resolution_RecoilJets"),
partJet.phi() - jetMatched.phi(), partJet.pt(), weight);
}
}
}
}
} else {
// Miss jets
spectra.fill(HIST("hMissedJets_pT"), partJet.pt(), weight);
Expand All @@ -973,7 +974,8 @@ struct RecoilJets {

// Fake jets
for (const auto& detJet : detJets) {
if(isJetWithHighPtConstituent(detJet, tracks)) continue;
if (isJetWithHighPtConstituent(detJet, tracks))
continue;

bIsThereMatchedJet = detJet.has_matchedJetGeo();
if (!bIsThereMatchedJet) {
Expand Down
Loading