@@ -530,16 +530,17 @@ struct JetHadronRecoil {
530530
531531 float dphip = RecoDecay::constrainAngle (jetTag.phi () - phiTTPart);
532532 dRp = getWTAaxisDifference (jetTag, mcpjetsWTA, particles, true );
533-
534- registry.fill (HIST (" hPtMatched" ), jetBase.pt () - (rho * jetBase.area ()), jetTag.pt (), weight);
535- registry.fill (HIST (" hPhiMatched" ), dphi, dphip, weight);
536- registry.fill (HIST (" hPtResolution" ), jetTag.pt (), (jetTag.pt () - (jetBase.pt () - (rho * jetBase.area ()))) / jetTag.pt (), weight);
537- registry.fill (HIST (" hPhiResolution" ), dphip, dphip - dphi, weight);
538- registry.fill (HIST (" hDeltaRMatched" ), dR, dRp, weight);
539- registry.fill (HIST (" hDeltaRResolution" ), jetTag.pt (), dRp - dR, weight);
540- registry.fill (HIST (" hFullMatching" ), jetBase.pt () - (rho * jetBase.area ()), jetTag.pt (), dphi, dphip, dR, dRp, weight);
541- registry.fill (HIST (" hPtMatched1d" ), jetTag.pt (), weight);
542- registry.fill (HIST (" hDeltaRMatched1d" ), dRp, weight);
533+ if ((std::abs (dphi - o2::constants::math::PI) < 0.6 ) || (std::abs (dphip - o2::constants::math::PI) < 0.6 )) {
534+ registry.fill (HIST (" hPtMatched" ), jetBase.pt () - (rho * jetBase.area ()), jetTag.pt (), weight);
535+ registry.fill (HIST (" hPhiMatched" ), dphi, dphip, weight);
536+ registry.fill (HIST (" hPtResolution" ), jetTag.pt (), (jetTag.pt () - (jetBase.pt () - (rho * jetBase.area ()))) / jetTag.pt (), weight);
537+ registry.fill (HIST (" hPhiResolution" ), dphip, dphip - dphi, weight);
538+ registry.fill (HIST (" hDeltaRMatched" ), dR, dRp, weight);
539+ registry.fill (HIST (" hDeltaRResolution" ), jetTag.pt (), dRp - dR, weight);
540+ registry.fill (HIST (" hFullMatching" ), jetBase.pt () - (rho * jetBase.area ()), jetTag.pt (), dphi, dphip, dR, dRp, weight);
541+ registry.fill (HIST (" hPtMatched1d" ), jetTag.pt (), weight);
542+ registry.fill (HIST (" hDeltaRMatched1d" ), dRp, weight);
543+ }
543544 }
544545 }
545546 }
0 commit comments