Skip to content

Commit

Permalink
PWGJE:adding track histograms to QA task
Browse files Browse the repository at this point in the history
  • Loading branch information
Nima Zardoshti committed Oct 7, 2023
1 parent b92e7a3 commit 02e1a18
Showing 1 changed file with 22 additions and 0 deletions.
22 changes: 22 additions & 0 deletions PWGJE/Tasks/jetfinderQA.cxx
Original file line number Diff line number Diff line change
Expand Up @@ -83,6 +83,9 @@ struct JetFinderQATask {
registry.add("h3_jet_r_jet_pt_track_eta", "#it{R}_{jet};#it{p}_{T,jet} (GeV/#it{c});#eta_{jet tracks}", {HistType::kTH3F, {{jetRadiiBins, ""}, {200, 0., 200.}, {100, -1.0, 1.0}}});
registry.add("h3_jet_r_jet_pt_track_phi", "#it{R}_{jet};#it{p}_{T,jet} (GeV/#it{c});#varphi_{jet tracks}", {HistType::kTH3F, {{jetRadiiBins, ""}, {200, 0., 200.}, {80, -1.0, 7.}}});
registry.add("h3_jet_r_jet_pt_leadingtrack_pt", "#it{R}_{jet};#it{p}_{T,jet} (GeV/#it{c}); #it{p}_{T,leading track} (GeV/#it{c})", {HistType::kTH3F, {{jetRadiiBins, ""}, {200, 0.0, 200.0}, {200, 0.0, 200.0}}});
registry.add("h_track_pt", "track pT;#it{p}_{T,track} (GeV/#it{c});entries", {HistType::kTH1F, {{200, 0., 200.}}});
registry.add("h_track_eta", "track #eta;#eta_{track};entries", {HistType::kTH1F, {{100, -1.0, 1.0}}});
registry.add("h_track_phi", "track #varphi;#varphi_{track};entries", {HistType::kTH1F, {{80, -1.0, 7.}}});
}

if (doprocessJetsMCP || doprocessJetsMCPWeighted) {
Expand Down Expand Up @@ -235,6 +238,9 @@ struct JetFinderQATask {
if (!selectTrack(track)) {
continue;
}
registry.fill(HIST("h_track_pt"), track.pt());
registry.fill(HIST("h_track_eta"), track.eta());
registry.fill(HIST("h_track_phi"), track.phi());
if (track.pt() > leadingTrackpT)
leadingTrackpT = track.pt();
}
Expand All @@ -245,12 +251,28 @@ struct JetFinderQATask {
void processJetsMCD(soa::Join<aod::ChargedMCDetectorLevelJets, aod::ChargedMCDetectorLevelJetConstituents>::iterator const& jet, JetTracks const& tracks)
{
fillHistograms(jet);
for (auto const& track : tracks) {
if (!selectTrack(track)) {
continue;
}
registry.fill(HIST("h_track_pt"), track.pt());
registry.fill(HIST("h_track_eta"), track.eta());
registry.fill(HIST("h_track_phi"), track.phi());
}
}
PROCESS_SWITCH(JetFinderQATask, processJetsMCD, "jet finder QA mcd", false);

void processJetsMCDWeighted(soa::Join<aod::ChargedMCDetectorLevelJets, aod::ChargedMCDetectorLevelJetConstituents, aod::ChargedMCDetectorLevelJetEventWeights>::iterator const& jet, JetTracks const& tracks)
{
fillHistograms(jet, jet.eventWeight());
for (auto const& track : tracks) {
if (!selectTrack(track)) {
continue;
}
registry.fill(HIST("h_track_pt"), track.pt(), jet.eventWeight());
registry.fill(HIST("h_track_eta"), track.eta(), jet.eventWeight());
registry.fill(HIST("h_track_phi"), track.phi(), jet.eventWeight());
}
}
PROCESS_SWITCH(JetFinderQATask, processJetsMCDWeighted, "jet finder QA mcd with weighted events", false);

Expand Down

0 comments on commit 02e1a18

Please sign in to comment.