Changeset 1317 in svn for trunk/modules/JetPileUpSubtractor.cc

Timestamp:

Nov 4, 2013, 1:14:33 PM (11 years ago)

Author:

Pavel Demin

Message:

fix loop in JetPileUpSubtractor and add rho eta ranges for Isolation

File:

• trunk/modules/JetPileUpSubtractor.cc (modified) (3 diffs)

trunk/modules/JetPileUpSubtractor.cc

@@ -83 +83 @@
 void JetPileUpSubtractor::Process()
 {
-  Candidate *candidate;
+  Candidate *candidate, *object;
   TLorentzVector momentum, area;
   Double_t eta = 0.0;
   Double_t rho = 0.0;
-
-  if(!fRhoInputArray) return;
 
   // loop over all input candidates
@@ -100 +98 @@
     // find rho
     rho = 0.0;
-    while((candidate = static_cast<Candidate*>(fItRhoInputArray->Next())))
+    if(fRhoInputArray)
     {
-      if(eta >= candidate->Edges[0] && eta < candidate->Edges[1])
+      fItRhoInputArray->Reset();
+      while((object = static_cast<Candidate*>(fItRhoInputArray->Next())))
       {
-        rho = candidate->Momentum.Pt();
+        if(eta >= object->Edges[0] && eta < object->Edges[1])
+        {
+          rho = object->Momentum.Pt();
+        }
       }
-    }  
+    }
 
     // apply pile-up correction
@@ -112 +114 @@
 
     momentum -= rho * area;
-  
+
     if(momentum.Pt() <= fJetPTMin) continue;