Skip to content

Commit 4def2ea

Browse files
joonsukbaejoonsukbae
committed
PWGJE: updates for background subtraction methods
1 parent 8551e88 commit 4def2ea

File tree

3 files changed

+133
-52
lines changed

3 files changed

+133
-52
lines changed

‎PWGJE/DataModel/JetSubtraction.h‎

Lines changed: 20 additions & 8 deletions
Original file line numberDiff line numberDiff line change
@@ -93,54 +93,66 @@ DECLARE_SOA_TABLE(BkgChargedMcRhos, "AOD", "BkgCMcRho",
9393
bkgrho::RhoM,
9494
o2::soa::Marker<1>);
9595

96-
DECLARE_SOA_TABLE(BkgD0Rhos, "AOD", "BkgD0Rho",
96+
DECLARE_SOA_TABLE(BkgChargedPerpRhos, "AOD", "BkgCPerpRho",
9797
o2::soa::Index<>,
9898
bkgrho::Rho,
9999
bkgrho::RhoM,
100100
o2::soa::Marker<2>);
101101

102-
DECLARE_SOA_TABLE(BkgD0McRhos, "AOD", "BkgD0McRho",
102+
DECLARE_SOA_TABLE(BkgChargedMcPerpRhos, "AOD", "BkgCMcPerpRho",
103103
o2::soa::Index<>,
104104
bkgrho::Rho,
105105
bkgrho::RhoM,
106106
o2::soa::Marker<3>);
107107

108-
DECLARE_SOA_TABLE(BkgLcRhos, "AOD", "BkgLcRho",
108+
DECLARE_SOA_TABLE(BkgD0Rhos, "AOD", "BkgD0Rho",
109109
o2::soa::Index<>,
110110
bkgrho::Rho,
111111
bkgrho::RhoM,
112112
o2::soa::Marker<4>);
113113

114-
DECLARE_SOA_TABLE(BkgLcMcRhos, "AOD", "BkgLcMcRho",
114+
DECLARE_SOA_TABLE(BkgD0McRhos, "AOD", "BkgD0McRho",
115115
o2::soa::Index<>,
116116
bkgrho::Rho,
117117
bkgrho::RhoM,
118118
o2::soa::Marker<5>);
119119

120-
DECLARE_SOA_TABLE(BkgBplusRhos, "AOD", "BkgRho",
120+
DECLARE_SOA_TABLE(BkgLcRhos, "AOD", "BkgLcRho",
121121
o2::soa::Index<>,
122122
bkgrho::Rho,
123123
bkgrho::RhoM,
124124
o2::soa::Marker<6>);
125125

126-
DECLARE_SOA_TABLE(BkgBplusMcRhos, "AOD", "BkgBPMcRho",
126+
DECLARE_SOA_TABLE(BkgLcMcRhos, "AOD", "BkgLcMcRho",
127127
o2::soa::Index<>,
128128
bkgrho::Rho,
129129
bkgrho::RhoM,
130130
o2::soa::Marker<7>);
131131

132-
DECLARE_SOA_TABLE(BkgDielectronRhos, "AOD", "BkgDIELRho",
132+
DECLARE_SOA_TABLE(BkgBplusRhos, "AOD", "BkgRho",
133133
o2::soa::Index<>,
134134
bkgrho::Rho,
135135
bkgrho::RhoM,
136136
o2::soa::Marker<8>);
137137

138-
DECLARE_SOA_TABLE(BkgDielectronMcRhos, "AOD", "BkgDIELMcRho",
138+
DECLARE_SOA_TABLE(BkgBplusMcRhos, "AOD", "BkgBPMcRho",
139139
o2::soa::Index<>,
140140
bkgrho::Rho,
141141
bkgrho::RhoM,
142142
o2::soa::Marker<9>);
143143

144+
DECLARE_SOA_TABLE(BkgDielectronRhos, "AOD", "BkgDIELRho",
145+
o2::soa::Index<>,
146+
bkgrho::Rho,
147+
bkgrho::RhoM,
148+
o2::soa::Marker<10>);
149+
150+
DECLARE_SOA_TABLE(BkgDielectronMcRhos, "AOD", "BkgDIELMcRho",
151+
o2::soa::Index<>,
152+
bkgrho::Rho,
153+
bkgrho::RhoM,
154+
o2::soa::Marker<11>);
155+
144156
DECLARE_SOA_TABLE(JTrackSubs, "AOD", "JTrackSubs",
145157
o2::soa::Index<>,
146158
bkgcharged::JCollisionId,

‎PWGJE/TableProducer/rhoEstimator.cxx‎

Lines changed: 22 additions & 18 deletions
Original file line numberDiff line numberDiff line change
@@ -32,6 +32,8 @@ using namespace o2::framework::expressions;
3232
struct RhoEstimatorTask {
3333
Produces<aod::BkgChargedRhos> rhoChargedTable;
3434
Produces<aod::BkgChargedMcRhos> rhoChargedMcTable;
35+
Produces<aod::BkgChargedPerpRhos> rhoChargedPerpTable;
36+
Produces<aod::BkgChargedMcPerpRhos> rhoChargedMcPerpTable;
3537
Produces<aod::BkgD0Rhos> rhoD0Table;
3638
Produces<aod::BkgD0McRhos> rhoD0McTable;
3739
Produces<aod::BkgLcRhos> rhoLcTable;
@@ -88,7 +90,9 @@ struct RhoEstimatorTask {
8890
inputParticles.clear();
8991
jetfindingutilities::analyseTracks<soa::Filtered<aod::JetTracks>, soa::Filtered<aod::JetTracks>::iterator>(inputParticles, tracks, trackSelection, trackingEfficiency);
9092
auto [rho, rhoM] = bkgSub.estimateRhoAreaMedian(inputParticles, doSparse);
93+
auto [rhoPerpCone, rhoMPerpCone] = bkgSub.estimateRhoPerpCone(inputParticles, inputParticles);
9194
rhoChargedTable(rho, rhoM);
95+
rhoChargedPerpTable(rhoPerpCone, rhoMPerpCone);
9296
}
9397
PROCESS_SWITCH(RhoEstimatorTask, processChargedCollisions, "Fill rho tables for collisions using charged tracks", true);
9498

@@ -153,31 +157,31 @@ struct RhoEstimatorTask {
153157
}
154158
PROCESS_SWITCH(RhoEstimatorTask, processLcMcCollisions, "Fill rho tables for collisions with Lc MCP candidates", false);
155159

156-
void processBplusCollisions(aod::JetCollision const&, soa::Filtered<aod::JetTracks> const& tracks, aod::CandidatesBplusData const& candidates)
157-
{
160+
void processBplusCollisions(aod::JetCollision const&, soa::Filtered<aod::JetTracks> const& tracks, aod::CandidatesBplusData const& candidates)
161+
{
162+
inputParticles.clear();
163+
for (auto& candidate : candidates) {
158164
inputParticles.clear();
159-
for (auto& candidate : candidates) {
160-
inputParticles.clear();
161-
jetfindingutilities::analyseTracks(inputParticles, tracks, trackSelection, trackingEfficiency, std::optional{candidate});
165+
jetfindingutilities::analyseTracks(inputParticles, tracks, trackSelection, trackingEfficiency, std::optional{candidate});
162166

163-
auto [rho, rhoM] = bkgSub.estimateRhoAreaMedian(inputParticles, doSparse);
164-
rhoBplusTable(rho, rhoM);
165-
}
167+
auto [rho, rhoM] = bkgSub.estimateRhoAreaMedian(inputParticles, doSparse);
168+
rhoBplusTable(rho, rhoM);
166169
}
167-
PROCESS_SWITCH(RhoEstimatorTask, processBplusCollisions, "Fill rho tables for collisions with Bplus candidates", false);
170+
}
171+
PROCESS_SWITCH(RhoEstimatorTask, processBplusCollisions, "Fill rho tables for collisions with Bplus candidates", false);
168172

169-
void processBplusMcCollisions(aod::JetMcCollision const&, soa::Filtered<aod::JetParticles> const& particles, aod::CandidatesBplusMCP const& candidates)
170-
{
173+
void processBplusMcCollisions(aod::JetMcCollision const&, soa::Filtered<aod::JetParticles> const& particles, aod::CandidatesBplusMCP const& candidates)
174+
{
175+
inputParticles.clear();
176+
for (auto& candidate : candidates) {
171177
inputParticles.clear();
172-
for (auto& candidate : candidates) {
173-
inputParticles.clear();
174-
jetfindingutilities::analyseParticles<true>(inputParticles, particleSelection, 1, particles, pdgDatabase, std::optional{candidate});
178+
jetfindingutilities::analyseParticles<true>(inputParticles, particleSelection, 1, particles, pdgDatabase, std::optional{candidate});
175179

176-
auto [rho, rhoM] = bkgSub.estimateRhoAreaMedian(inputParticles, doSparse);
177-
rhoBplusMcTable(rho, rhoM);
178-
}
180+
auto [rho, rhoM] = bkgSub.estimateRhoAreaMedian(inputParticles, doSparse);
181+
rhoBplusMcTable(rho, rhoM);
179182
}
180-
PROCESS_SWITCH(RhoEstimatorTask, processBplusMcCollisions, "Fill rho tables for collisions with Bplus MCP candidates", false);
183+
}
184+
PROCESS_SWITCH(RhoEstimatorTask, processBplusMcCollisions, "Fill rho tables for collisions with Bplus MCP candidates", false);
181185

182186
void processDielectronCollisions(aod::JetCollision const&, soa::Filtered<aod::JetTracks> const& tracks, aod::CandidatesDielectronData const& candidates)
183187
{

‎PWGJE/Tasks/jetFinderQA.cxx‎

Lines changed: 91 additions & 26 deletions
Original file line numberDiff line numberDiff line change
@@ -156,8 +156,10 @@ struct JetFinderQATask {
156156
registry.add("h3_centrality_occupancy_jet_pt", "centrality; occupancy; #it{p}_{T,jet} (GeV/#it{c})", {HistType::kTH3F, {{120, -10.0, 110.0}, {60, 0, 30000}, jetPtAxis}});
157157
}
158158

159-
if (doprocessJetsRhoAreaSubData || doprocessJetsRhoAreaSubMCD) {
159+
if (doprocessJetsRhoAreaSubData || doprocessJetsRhoAreaSubMCD || doprocessJetsRhoPerpSubData || doprocessJetsRhoPerpSubMCD) {
160160

161+
registry.add("h2_leadingjet_pt_rho", "#it{p}_{T,leading jet} vs #rho (GeV/#it{c});#it{p}_{T,leading jet} (GeV/#it{c});#rho (GeV/#it{c})", {HistType::kTH2F, {jetPtAxis, {1000, 0.0, 10.0}}});
162+
registry.add("h2_leadingjet_pt_rhoM", "#it{p}_{T,leading jet} vs #rho_{M} (GeV/#it{c});#it{p}_{T,jet} (GeV/#it{c});#rho_{M} (GeV/#it{c})", {HistType::kTH2F, {jetPtAxis, {1000, 0.0, 10.0}}});
161163
registry.add("h_jet_pt_rhoareasubtracted", "jet pT;#it{p}_{T,jet} (GeV/#it{c});entries", {HistType::kTH1F, {jetPtAxisRhoAreaSub}});
162164
registry.add("h_jet_eta_rhoareasubtracted", "jet #eta;#eta_{jet};entries", {HistType::kTH1F, {jetEtaAxis}});
163165
registry.add("h_jet_phi_rhoareasubtracted", "jet #varphi;#varphi_{jet};entries", {HistType::kTH1F, {{160, -1.0, 7.}}});
@@ -200,19 +202,19 @@ struct JetFinderQATask {
200202
}
201203

202204
if (doprocessRho) {
203-
registry.add("h2_centrality_ntracks", "; centrality; N_{tracks};", {HistType::kTH2F, {{1100, 0., 110.0}, {10000, 0.0, 10000.0}}});
205+
registry.add("h2_centrality_ntracks", "; centrality; N_{tracks};", {HistType::kTH2F, {{1200, -10.0, 110.0}, {10000, 0.0, 10000.0}}});
204206
registry.add("h2_ntracks_rho", "; N_{tracks}; #it{rho} (GeV/area);", {HistType::kTH2F, {{10000, 0.0, 10000.0}, {400, 0.0, 400.0}}});
205207
registry.add("h2_ntracks_rhom", "; N_{tracks}; #it{rho}_{m} (GeV/area);", {HistType::kTH2F, {{10000, 0.0, 10000.0}, {100, 0.0, 100.0}}});
206-
registry.add("h2_centrality_rho", "; centrality; #it{rho} (GeV/area);", {HistType::kTH2F, {{1100, 0., 110.}, {400, 0., 400.0}}});
207-
registry.add("h2_centrality_rhom", ";centrality; #it{rho}_{m} (GeV/area)", {HistType::kTH2F, {{1100, 0., 110.}, {100, 0., 100.0}}});
208+
registry.add("h2_centrality_rho", "; centrality; #it{rho} (GeV/area);", {HistType::kTH2F, {{1200, -10.0, 110.0}, {400, 0., 400.0}}});
209+
registry.add("h2_centrality_rhom", ";centrality; #it{rho}_{m} (GeV/area)", {HistType::kTH2F, {{1200, -10.0, 110.0}, {100, 0., 100.0}}});
208210
}
209211

210212
if (doprocessRandomConeData || doprocessRandomConeMCD) {
211-
registry.add("h2_centrality_rhorandomcone", "; centrality; #it{p}_{T,random cone} - #it{area, random cone} * #it{rho} (GeV/c);", {HistType::kTH2F, {{1100, 0., 110.}, {800, -400.0, 400.0}}});
212-
registry.add("h2_centrality_rhorandomconerandomtrackdirection", "; centrality; #it{p}_{T,random cone} - #it{area, random cone} * #it{rho} (GeV/c);", {HistType::kTH2F, {{1100, 0., 110.}, {800, -400.0, 400.0}}});
213-
registry.add("h2_centrality_rhorandomconewithoutleadingjet", "; centrality; #it{p}_{T,random cone} - #it{area, random cone} * #it{rho} (GeV/c);", {HistType::kTH2F, {{1100, 0., 110.}, {800, -400.0, 400.0}}});
214-
registry.add("h2_centrality_rhorandomconerandomtrackdirectionwithoutoneleadingjets", "; centrality; #it{p}_{T,random cone} - #it{area, random cone} * #it{rho} (GeV/c);", {HistType::kTH2F, {{1100, 0., 110.}, {800, -400.0, 400.0}}});
215-
registry.add("h2_centrality_rhorandomconerandomtrackdirectionwithouttwoleadingjets", "; centrality; #it{p}_{T,random cone} - #it{area, random cone} * #it{rho} (GeV/c);", {HistType::kTH2F, {{1100, 0., 110.}, {800, -400.0, 400.0}}});
213+
registry.add("h2_centrality_rhorandomcone", "; centrality; #it{p}_{T,random cone} - #it{area, random cone} * #it{rho} (GeV/c);", {HistType::kTH2F, {{1200, -10.0, 110.0}, {800, -400.0, 400.0}}});
214+
registry.add("h2_centrality_rhorandomconerandomtrackdirection", "; centrality; #it{p}_{T,random cone} - #it{area, random cone} * #it{rho} (GeV/c);", {HistType::kTH2F, {{1200, -10.0, 110.0}, {800, -400.0, 400.0}}});
215+
registry.add("h2_centrality_rhorandomconewithoutleadingjet", "; centrality; #it{p}_{T,random cone} - #it{area, random cone} * #it{rho} (GeV/c);", {HistType::kTH2F, {{1200, -10.0, 110.0}, {800, -400.0, 400.0}}});
216+
registry.add("h2_centrality_rhorandomconerandomtrackdirectionwithoutoneleadingjets", "; centrality; #it{p}_{T,random cone} - #it{area, random cone} * #it{rho} (GeV/c);", {HistType::kTH2F, {{1200, -10.0, 110.0}, {800, -400.0, 400.0}}});
217+
registry.add("h2_centrality_rhorandomconerandomtrackdirectionwithouttwoleadingjets", "; centrality; #it{p}_{T,random cone} - #it{area, random cone} * #it{rho} (GeV/c);", {HistType::kTH2F, {{1200, -10.0, 110.0}, {800, -400.0, 400.0}}});
216218
}
217219

218220
if (doprocessJetsMCP || doprocessJetsMCPWeighted) {
@@ -664,29 +666,30 @@ struct JetFinderQATask {
664666
}
665667
}
666668
}
667-
}
668-
669-
registry.fill(HIST("h2_centrality_rhorandomconewithoutleadingjet"), collision.centrality(), randomConePt - M_PI * randomConeR * randomConeR * collision.rho());
670669

671-
// randomised eta,phi for tracks, to assess part of fluctuations coming from statistically independently emitted particles, removing tracks from 2 leading jets
672-
double randomConePtWithoutOneLeadJet = 0;
673-
double randomConePtWithoutTwoLeadJet = 0;
674-
for (auto const& track : tracks) {
675-
if (jetderiveddatautilities::selectTrack(track, trackSelection)) {
676-
float dPhi = RecoDecay::constrainAngle(randomNumber.Uniform(0.0, 2 * M_PI) - randomConePhi, static_cast<float>(-M_PI)); // ignores actual phi of track
677-
float dEta = randomNumber.Uniform(trackEtaMin, trackEtaMax) - randomConeEta; // ignores actual eta of track
678-
if (TMath::Sqrt(dEta * dEta + dPhi * dPhi) < randomConeR) {
679-
if (!trackIsInJet(track, jets.iteratorAt(0))) {
680-
randomConePtWithoutOneLeadJet += track.pt();
681-
if (!trackIsInJet(track, jets.iteratorAt(1))) {
682-
randomConePtWithoutTwoLeadJet += track.pt();
670+
registry.fill(HIST("h2_centrality_rhorandomconewithoutleadingjet"), collision.centrality(), randomConePt - M_PI * randomConeR * randomConeR * collision.rho());
671+
672+
// randomised eta,phi for tracks, to assess part of fluctuations coming from statistically independently emitted particles, removing tracks from 2 leading jets
673+
double randomConePtWithoutOneLeadJet = 0;
674+
double randomConePtWithoutTwoLeadJet = 0;
675+
for (auto const& track : tracks) {
676+
if (jetderiveddatautilities::selectTrack(track, trackSelection)) {
677+
float dPhi = RecoDecay::constrainAngle(randomNumber.Uniform(0.0, 2 * M_PI) - randomConePhi, static_cast<float>(-M_PI)); // ignores actual phi of track
678+
float dEta = randomNumber.Uniform(trackEtaMin, trackEtaMax) - randomConeEta; // ignores actual eta of track
679+
if (TMath::Sqrt(dEta * dEta + dPhi * dPhi) < randomConeR) {
680+
if (!trackIsInJet(track, jets.iteratorAt(0))) {
681+
randomConePtWithoutOneLeadJet += track.pt();
682+
if (jets.size() > 1 && !trackIsInJet(track, jets.iteratorAt(1))) {
683+
randomConePtWithoutTwoLeadJet += track.pt();
684+
}
683685
}
684686
}
685687
}
686688
}
689+
690+
registry.fill(HIST("h2_centrality_rhorandomconerandomtrackdirectionwithoutoneleadingjets"), collision.centrality(), randomConePtWithoutOneLeadJet - M_PI * randomConeR * randomConeR * collision.rho());
691+
registry.fill(HIST("h2_centrality_rhorandomconerandomtrackdirectionwithouttwoleadingjets"), collision.centrality(), randomConePtWithoutTwoLeadJet - M_PI * randomConeR * randomConeR * collision.rho());
687692
}
688-
registry.fill(HIST("h2_centrality_rhorandomconerandomtrackdirectionwithoutoneleadingjets"), collision.centrality(), randomConePtWithoutOneLeadJet - M_PI * randomConeR * randomConeR * collision.rho());
689-
registry.fill(HIST("h2_centrality_rhorandomconerandomtrackdirectionwithouttwoleadingjets"), collision.centrality(), randomConePtWithoutTwoLeadJet - M_PI * randomConeR * randomConeR * collision.rho());
690693
}
691694

692695
void processJetsData(soa::Filtered<aod::JetCollisions>::iterator const& collision, soa::Join<aod::ChargedJets, aod::ChargedJetConstituents> const& jets, aod::JetTracks const&)
@@ -713,13 +716,19 @@ struct JetFinderQATask {
713716
if (collision.trackOccupancyInTimeRange() < trackOccupancyInTimeRangeMin || trackOccupancyInTimeRangeMax < collision.trackOccupancyInTimeRange()) {
714717
return;
715718
}
719+
bool leadingJetFilled = false;
716720
for (auto jet : jets) {
717721
if (!jetfindingutilities::isInEtaAcceptance(jet, jetEtaMin, jetEtaMax, trackEtaMin, trackEtaMax)) {
718722
continue;
719723
}
720724
if (!isAcceptedJet<aod::JetTracks>(jet)) {
721725
continue;
722726
}
727+
if (!leadingJetFilled) {
728+
registry.fill(HIST("h2_leadingjet_pt_rho"), jet.pt(), collision.rho());
729+
registry.fill(HIST("h2_leadingjet_pt_rhoM"), jet.pt(), collision.rhoM());
730+
leadingJetFilled = true;
731+
}
723732
fillRhoAreaSubtractedHistograms(jet, collision.centrality(), collision.trackOccupancyInTimeRange(), collision.rho());
724733
}
725734
}
@@ -732,18 +741,74 @@ struct JetFinderQATask {
732741
if (collision.trackOccupancyInTimeRange() < trackOccupancyInTimeRangeMin || trackOccupancyInTimeRangeMax < collision.trackOccupancyInTimeRange()) {
733742
return;
734743
}
744+
bool leadingJetFilled = false;
735745
for (auto jet : jets) {
736746
if (!jetfindingutilities::isInEtaAcceptance(jet, jetEtaMin, jetEtaMax, trackEtaMin, trackEtaMax)) {
737747
continue;
738748
}
739749
if (!isAcceptedJet<aod::JetTracks>(jet)) {
740750
continue;
741751
}
752+
if (!leadingJetFilled) {
753+
registry.fill(HIST("h2_leadingjet_pt_rho"), jet.pt(), collision.rho());
754+
registry.fill(HIST("h2_leadingjet_pt_rhoM"), jet.pt(), collision.rhoM());
755+
leadingJetFilled = true;
756+
}
742757
fillRhoAreaSubtractedHistograms(jet, collision.centrality(), collision.trackOccupancyInTimeRange(), collision.rho());
743758
}
744759
}
745760
PROCESS_SWITCH(JetFinderQATask, processJetsRhoAreaSubMCD, "jet finder QA for rho-area subtracted mcd jets", false);
746761

762+
void processJetsRhoPerpSubData(soa::Filtered<soa::Join<aod::JetCollisions, aod::BkgChargedPerpRhos>>::iterator const& collision,
763+
soa::Join<aod::ChargedJets, aod::ChargedJetConstituents> const& jets,
764+
aod::JetTracks const&)
765+
{
766+
if (collision.trackOccupancyInTimeRange() < trackOccupancyInTimeRangeMin || trackOccupancyInTimeRangeMax < collision.trackOccupancyInTimeRange()) {
767+
return;
768+
}
769+
bool leadingJetFilled = false;
770+
for (auto jet : jets) {
771+
if (!jetfindingutilities::isInEtaAcceptance(jet, jetEtaMin, jetEtaMax, trackEtaMin, trackEtaMax)) {
772+
continue;
773+
}
774+
if (!isAcceptedJet<aod::JetTracks>(jet)) {
775+
continue;
776+
}
777+
if (!leadingJetFilled) {
778+
registry.fill(HIST("h2_leadingjet_pt_rho"), jet.pt(), collision.rho());
779+
registry.fill(HIST("h2_leadingjet_pt_rhoM"), jet.pt(), collision.rhoM());
780+
leadingJetFilled = true;
781+
}
782+
fillRhoAreaSubtractedHistograms(jet, collision.centrality(), collision.trackOccupancyInTimeRange(), collision.rho());
783+
}
784+
}
785+
PROCESS_SWITCH(JetFinderQATask, processJetsRhoPerpSubData, "jet finder QA for rho-perpendicular-cone subtracted jets", false);
786+
787+
void processJetsRhoPerpSubMCD(soa::Filtered<soa::Join<aod::JetCollisions, aod::BkgChargedPerpRhos>>::iterator const& collision,
788+
soa::Join<aod::ChargedMCDetectorLevelJets, aod::ChargedMCDetectorLevelJetConstituents> const& jets,
789+
aod::JetTracks const&)
790+
{
791+
if (collision.trackOccupancyInTimeRange() < trackOccupancyInTimeRangeMin || trackOccupancyInTimeRangeMax < collision.trackOccupancyInTimeRange()) {
792+
return;
793+
}
794+
bool leadingJetFilled = false;
795+
for (auto jet : jets) {
796+
if (!jetfindingutilities::isInEtaAcceptance(jet, jetEtaMin, jetEtaMax, trackEtaMin, trackEtaMax)) {
797+
continue;
798+
}
799+
if (!isAcceptedJet<aod::JetTracks>(jet)) {
800+
continue;
801+
}
802+
if (!leadingJetFilled) {
803+
registry.fill(HIST("h2_leadingjet_pt_rho"), jet.pt(), collision.rho());
804+
registry.fill(HIST("h2_leadingjet_pt_rhoM"), jet.pt(), collision.rhoM());
805+
leadingJetFilled = true;
806+
}
807+
fillRhoAreaSubtractedHistograms(jet, collision.centrality(), collision.trackOccupancyInTimeRange(), collision.rho());
808+
}
809+
}
810+
PROCESS_SWITCH(JetFinderQATask, processJetsRhoPerpSubMCD, "jet finder QA for rho-perpendicular-cone subtracted mcd jets", false);
811+
747812
void processEvtWiseConstSubJetsData(soa::Filtered<aod::JetCollisions>::iterator const& collision, soa::Join<aod::ChargedEventWiseSubtractedJets, aod::ChargedEventWiseSubtractedJetConstituents> const& jets, aod::JetTracksSub const&)
748813
{
749814
if (collision.trackOccupancyInTimeRange() < trackOccupancyInTimeRangeMin || trackOccupancyInTimeRangeMax < collision.trackOccupancyInTimeRange()) {

0 commit comments

Comments
 (0)