From 3bb3cc92f3f7d8243949d17282cfcd7e83564b9f Mon Sep 17 00:00:00 2001
From: perrotta <perrotta@cern.ch>
Date: Sat, 21 Nov 2020 18:38:17 +0100
Subject: [PATCH] Update PFClusterFromHGCalMultiCluster.cc

Check if the multicluster is empty and move to the next one before  having to compute the trackster PID, as it was implemented in the 11_1 "backport" PR #31907 (it avoids computing the probabilities when not needed)
@rovere @felicepantaleo
---
 .../plugins/PFClusterFromHGCalMultiCluster.cc         | 11 ++++-------
 1 file changed, 4 insertions(+), 7 deletions(-)

diff --git a/RecoParticleFlow/PFClusterProducer/plugins/PFClusterFromHGCalMultiCluster.cc b/RecoParticleFlow/PFClusterProducer/plugins/PFClusterFromHGCalMultiCluster.cc
index dd4352de9719e..4d78d73b01e75 100644
--- a/RecoParticleFlow/PFClusterProducer/plugins/PFClusterFromHGCalMultiCluster.cc
+++ b/RecoParticleFlow/PFClusterProducer/plugins/PFClusterFromHGCalMultiCluster.cc
@@ -26,26 +26,23 @@ void PFClusterFromHGCalMultiCluster::buildClusters(const edm::Handle<reco::PFRec
   }
 
   int iMultiClus = -1;
-
   for (const auto& mcl : hgcalMultiClusters) {
     iMultiClus++;
-
+    // Skip empty multiclusters
+    if (mcl.hitsAndFractions().empty()) {
+      continue;
+    }
     // Filter using trackster PID
     if (filterByTracksterPID_) {
       float probTotal = 0.0f;
-
       for (int cat : filter_on_categories_) {
         probTotal += tracksters[iMultiClus].id_probabilities(cat);
       }
-
       if (probTotal < pid_threshold_) {
         continue;
       }
     }
 
-    if (mcl.hitsAndFractions().empty()) {
-      continue;
-    }
     DetId seed;
     double energy = 0.0, highest_energy = 0.0;
     output.emplace_back();