@@ -350,7 +350,7 @@ struct JetHadronRecoil {
350350 registry.fill (HIST (" hJetPhi" ), jet.phi (), weight);
351351
352352 double dR = getWTAaxisDifference (jet, jetsWTA, particles);
353-
353+
354354 registry.fill (HIST (" hDeltaRPart" ), dR, weight);
355355 registry.fill (HIST (" hDeltaRpTPart" ), jet.pt (), dR, weight);
356356 if (nTT > 0 ) {
@@ -725,12 +725,13 @@ struct JetHadronRecoil {
725725 PROCESS_SWITCH (JetHadronRecoil, processRecoilJetsMCPMCDMatchedWeighted, " process MC matched with event weights (recoil jets)" , false );
726726
727727 template <typename T, typename U, typename X>
728- double getWTAaxisDifference (T const & jet, U const & jetsWTA, X const & /* tracks or particles*/ , bool isMatched = false ) {
728+ double getWTAaxisDifference (T const & jet, U const & jetsWTA, X const & /* tracks or particles*/ , bool isMatched = false )
729+ {
729730 double deltaPhi;
730731 double deltaEta;
731732 double deltaY;
732733 double dR;
733- if (wtaMethod == 0 ) {
734+ if (wtaMethod == 0 ) {
734735 // get WTA matched jet - should just be one jet matched geometrically
735736 if (isMatched) {
736737 // response - requires alternative method
@@ -741,17 +742,15 @@ struct JetHadronRecoil {
741742 break ;
742743 }
743744 }
744- }
745- else {
745+ } else {
746746 // MCP or MCD
747747 for (const auto & jetWTA : jet.template matchedJetGeo_as <std::decay_t <U>>()) {
748748 deltaPhi = RecoDecay::constrainAngle (jetWTA.phi () - jet.phi (), -o2::constants::math::PI);
749749 deltaEta = jetWTA.eta () - jet.eta ();
750750 dR = RecoDecay::sqrtSumOfSquares (deltaPhi, deltaEta);
751751 }
752752 }
753- }
754- else if (wtaMethod == 1 ) {
753+ } else if (wtaMethod == 1 ) {
755754 // recluster jet
756755 jetConstituents.clear ();
757756 for (auto & jetConstituent : jet.template tracks_as <X>()) {
@@ -764,13 +763,12 @@ struct JetHadronRecoil {
764763 deltaPhi = RecoDecay::constrainAngle (jet.phi () - jetReclustered[0 ].phi (), -o2::constants::math::PI);
765764 deltaY = jet.y () - jetReclustered[0 ].rap ();
766765 dR = RecoDecay::sqrtSumOfSquares (deltaPhi, deltaY);
767- LOG (debug) << " orig. jet n const = " << jet.tracksIds ().size () << " pt = " << jet.pt () << " eta = " << jet.eta () << " phi = " << jet.phi ();
768- LOG (debug) << " recl. jet n const = " << clusterSeq.constituents (jetReclustered[0 ]).size () << " pt = " << jetReclustered[0 ].pt () << " eta = " << jetReclustered[0 ].eta () << " phi = " << jetReclustered[0 ].phi ();
769- LOG (debug) << " distance = " << dR;
766+ LOG (debug) << " orig. jet n const = " << jet.tracksIds ().size () << " pt = " << jet.pt () << " eta = " << jet.eta () << " phi = " << jet.phi ();
767+ LOG (debug) << " recl. jet n const = " << clusterSeq.constituents (jetReclustered[0 ]).size () << " pt = " << jetReclustered[0 ].pt () << " eta = " << jetReclustered[0 ].eta () << " phi = " << jetReclustered[0 ].phi ();
768+ LOG (debug) << " distance = " << dR;
770769 }
771- return dR;
770+ return dR;
772771 }
773-
774772};
775773
776774WorkflowSpec defineDataProcessing (ConfigContext const & cfgc) { return WorkflowSpec{adaptAnalysisTask<JetHadronRecoil>(cfgc)}; }
0 commit comments