@@ -212,9 +212,9 @@ struct JetCorrelationD0 {
212212 if (jetBase.pt () > pTHatMaxMCD * pTHat || pTHat < pTHatAbsoluteMin) { // reject events that are too hard / soft
213213 return ;
214214 }
215- if (jetBase.has_matchedJetGeo ()) { // geometric matching
215+ if (jetBase.has_matchedJetGeo ()) { // geometric matching
216216 for (auto & jetTag : jetBase.template matchedJetGeo_as <std::decay_t <U>>()) {
217- if (jetTag.pt () > pTHatMaxMCP * pTHat) { // cuts overly hard jets from jettag (mcp) sample
217+ if (jetTag.pt () > pTHatMaxMCP * pTHat) { // cuts overly hard jets from jettag (mcp) sample
218218 continue ;
219219 }
220220 registry.fill (HIST (" hPtMatched" ), jetBase.pt () - (rho * jetBase.area ()), jetTag.pt (), weight);
@@ -278,7 +278,7 @@ struct JetCorrelationD0 {
278278
279279 for (const auto & d0Candidate : d0Candidates) {
280280 const auto scores = d0Candidate.mlScores ();
281- if (d0Candidate.pt () < d0PtCutMin) {
281+ if (d0Candidate.pt () < d0PtCutMin) {
282282 return ;
283283 }
284284 fillD0Histograms (d0Candidate, scores);
@@ -291,11 +291,11 @@ struct JetCorrelationD0 {
291291 d0Candidate.eta (),
292292 d0Candidate.phi ());
293293 for (const auto & jet : jets) {
294- if (jet.pt () < jetPtCutMin) {
294+ if (jet.pt () < jetPtCutMin) {
295295 return ;
296296 }
297297 float dphi = RecoDecay::constrainAngle (jet.phi () - d0Candidate.phi ());
298- if (abs (dphi - M_PI) > (M_PI /2 )) {
298+ if (abs (dphi - M_PI) > (M_PI / 2 )) {
299299 return ;
300300 }
301301 fillJetHistograms (jet, dphi);
@@ -314,11 +314,11 @@ struct JetCorrelationD0 {
314314 aod::CandidatesD0MCP const & d0MCPCandidates,
315315 soa::Filtered<soa::Join<aod::ChargedMCParticleLevelJets, aod::ChargedMCParticleLevelJetConstituents>> const & jets)
316316 {
317- // applyCollisionSelections(collision, eventSelectionBits);
317+ // applyCollisionSelections(collision, eventSelectionBits);
318318 tableCollision (collision.posZ ());
319319
320320 for (const auto & d0MCPCandidate : d0MCPCandidates) {
321- if (d0MCPCandidate.pt () < d0PtCutMin) {
321+ if (d0MCPCandidate.pt () < d0PtCutMin) {
322322 return ;
323323 }
324324 tableD0MCParticle (tableCollision.lastIndex (),
@@ -328,11 +328,11 @@ struct JetCorrelationD0 {
328328 d0MCPCandidate.phi ());
329329
330330 for (const auto & jet : jets) {
331- if (jet.pt () < jetPtCutMin) {
331+ if (jet.pt () < jetPtCutMin) {
332332 return ;
333333 }
334334 float dphi = RecoDecay::constrainAngle (jet.phi () - d0MCPCandidate.phi ());
335- if (abs (dphi - M_PI) > (M_PI /2 )) {
335+ if (abs (dphi - M_PI) > (M_PI / 2 )) {
336336 return ;
337337 }
338338 fillJetHistograms (jet, dphi);
@@ -356,7 +356,7 @@ struct JetCorrelationD0 {
356356 soa::Filtered<soa::Join<aod::ChargedMCParticleLevelJets, aod::ChargedMCParticleLevelJetConstituents, aod::ChargedMCDetectorLevelJetsMatchedToChargedMCParticleLevelJets>> const & mcpJets,
357357 soa::Filtered<soa::Join<aod::ChargedMCDetectorLevelJets, aod::ChargedMCDetectorLevelJetConstituents, aod::ChargedMCParticleLevelJetsMatchedToChargedMCDetectorLevelJets>> const & mcdJets)
358358 {
359- // applyCollisionSelections(collision, eventSelectionBits);
359+ // applyCollisionSelections(collision, eventSelectionBits);
360360
361361 const auto CollIdx = collision.mcCollisionId ();
362362
@@ -406,7 +406,7 @@ struct JetCorrelationD0 {
406406 LOGF (info, " Collision ID %i, D0 pt %.2f, D0 eta %.2f, D0 phi %.2f, MCP origin %hhd, Reflection %i" , CollIdx, particleMother.pt (), particleMother.eta (), particleMother.phi (), particleMother.originMcGen (), reflection);
407407
408408 // Jet matching
409- fillMatchedHistograms (mcdJets, mcpJets); // Do I need to include pthat cuts in loop rather than this function to actually do jet matching?
409+ fillMatchedHistograms (mcdJets, mcpJets); // Do I need to include pthat cuts in loop rather than this function to actually do jet matching?
410410
411411 for (const auto & mcpJet : mcpJets) {
412412 float dphi = RecoDecay::constrainAngle (mcpJet.phi () - d0MCDCandidate.phi ());
0 commit comments