diff --git a/PWGJE/Tasks/jetCorrelationD0.cxx b/PWGJE/Tasks/jetCorrelationD0.cxx
index c04500a87b4..3f877916a32 100644
--- a/PWGJE/Tasks/jetCorrelationD0.cxx
+++ b/PWGJE/Tasks/jetCorrelationD0.cxx
@@ -212,9 +212,9 @@ struct JetCorrelationD0 {
       if (jetBase.pt() > pTHatMaxMCD * pTHat || pTHat < pTHatAbsoluteMin) { // reject events that are too hard / soft
         return;
       }
-      if (jetBase.has_matchedJetGeo()) {                                            // geometric matching
+      if (jetBase.has_matchedJetGeo()) { // geometric matching
         for (auto& jetTag : jetBase.template matchedJetGeo_as<std::decay_t<U>>()) {
-          if (jetTag.pt() > pTHatMaxMCP * pTHat) {                                  // cuts overly hard jets from jettag (mcp) sample
+          if (jetTag.pt() > pTHatMaxMCP * pTHat) { // cuts overly hard jets from jettag (mcp) sample
             continue;
           }
           registry.fill(HIST("hPtMatched"), jetBase.pt() - (rho * jetBase.area()), jetTag.pt(), weight);
@@ -278,7 +278,7 @@ struct JetCorrelationD0 {
 
     for (const auto& d0Candidate : d0Candidates) {
       const auto scores = d0Candidate.mlScores();
-      if(d0Candidate.pt() < d0PtCutMin) {
+      if (d0Candidate.pt() < d0PtCutMin) {
         return;
       }
       fillD0Histograms(d0Candidate, scores);
@@ -291,11 +291,11 @@ struct JetCorrelationD0 {
               d0Candidate.eta(),
               d0Candidate.phi());
       for (const auto& jet : jets) {
-        if(jet.pt() < jetPtCutMin) {
+        if (jet.pt() < jetPtCutMin) {
           return;
         }
         float dphi = RecoDecay::constrainAngle(jet.phi() - d0Candidate.phi());
-                if (abs(dphi - M_PI) > (M_PI /2)) {
+        if (abs(dphi - M_PI) > (M_PI / 2)) {
           return;
         }
         fillJetHistograms(jet, dphi);
@@ -318,7 +318,7 @@ struct JetCorrelationD0 {
     tableCollision(collision.posZ());
 
     for (const auto& d0MCPCandidate : d0MCPCandidates) {
-      if(d0MCPCandidate.pt() < d0PtCutMin) {
+      if (d0MCPCandidate.pt() < d0PtCutMin) {
         return;
       }
       tableD0MCParticle(tableCollision.lastIndex(),
@@ -328,11 +328,11 @@ struct JetCorrelationD0 {
                         d0MCPCandidate.phi());
 
       for (const auto& jet : jets) {
-        if(jet.pt() < jetPtCutMin) {
+        if (jet.pt() < jetPtCutMin) {
           return;
         }
         float dphi = RecoDecay::constrainAngle(jet.phi() - d0MCPCandidate.phi());
-                if (abs(dphi - M_PI) > (M_PI /2)) {
+        if (abs(dphi - M_PI) > (M_PI / 2)) {
           return;
         }
         fillJetHistograms(jet, dphi);
@@ -406,7 +406,7 @@ struct JetCorrelationD0 {
         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);
 
         // Jet matching
-        fillMatchedHistograms(mcdJets, mcpJets); // Do I need to include pthat cuts in loop rather than this function to actually do jet matching? 
+        fillMatchedHistograms(mcdJets, mcpJets); // Do I need to include pthat cuts in loop rather than this function to actually do jet matching?
 
         for (const auto& mcpJet : mcpJets) {
           float dphi = RecoDecay::constrainAngle(mcpJet.phi() - d0MCDCandidate.phi());