Changes in modules/Isolation.cc [65dc4e7:341014c] in git

File:

• modules/Isolation.cc (modified) (10 diffs)

modules/Isolation.cc

@@ -17 +17 @@
  */
 
-
 /** \class Isolation
  *
@@ -35 +34 @@
 #include "classes/DelphesFormula.h"
 
+#include "ExRootAnalysis/ExRootClassifier.h"
+#include "ExRootAnalysis/ExRootFilter.h"
 #include "ExRootAnalysis/ExRootResult.h"
-#include "ExRootAnalysis/ExRootFilter.h"
-#include "ExRootAnalysis/ExRootClassifier.h"
-
+
+#include "TDatabasePDG.h"
+#include "TFormula.h"
+#include "TLorentzVector.h"
 #include "TMath.h"
+#include "TObjArray.h"
+#include "TRandom3.h"
 #include "TString.h"
-#include "TFormula.h"
-#include "TRandom3.h"
-#include "TObjArray.h"
-#include "TDatabasePDG.h"
-#include "TLorentzVector.h"
 
 #include <algorithm>
-#include <stdexcept>
 #include <iostream>
 #include <sstream>
+#include <stdexcept>
 
 using namespace std;
@@ -59 +58 @@
 {
 public:
-
   IsolationClassifier() {}
 
@@ -71 +69 @@
 Int_t IsolationClassifier::GetCategory(TObject *object)
 {
-  Candidate *track = static_cast<Candidate*>(object);
+  Candidate *track = static_cast<Candidate *>(object);
   const TLorentzVector &momentum = track->Momentum;
 
@@ -171 +169 @@
   // loop over all input jets
   fItCandidateInputArray->Reset();
-  while((candidate = static_cast<Candidate*>(fItCandidateInputArray->Next())))
+  while((candidate = static_cast<Candidate *>(fItCandidateInputArray->Next())))
   {
     const TLorentzVector &candidateMomentum = candidate->Momentum;
@@ -181 +179 @@
     {
       fItRhoInputArray->Reset();
-      while((object = static_cast<Candidate*>(fItRhoInputArray->Next())))
+      while((object = static_cast<Candidate *>(fItRhoInputArray->Next())))
       {
         if(eta >= object->Edges[0] && eta < object->Edges[1])
@@ -198 +196 @@
 
     itIsolationArray.Reset();
-    while((isolation = static_cast<Candidate*>(itIsolationArray.Next())))
+    while((isolation = static_cast<Candidate *>(itIsolationArray.Next())))
     {
       const TLorentzVector &isolationMomentum = isolation->Momentum;
@@ -204 +202 @@
       if(fUseMiniCone)
       {
-         pass = candidateMomentum.DeltaR(isolationMomentum) <= fDeltaRMax &&
-         candidateMomentum.DeltaR(isolationMomentum) > fDeltaRMin;
+        pass = candidateMomentum.DeltaR(isolationMomentum) <= fDeltaRMax && candidateMomentum.DeltaR(isolationMomentum) > fDeltaRMin;
       }
       else
       {
-         pass = candidateMomentum.DeltaR(isolationMomentum) <= fDeltaRMax &&
-         candidate->GetUniqueID() != isolation->GetUniqueID();
+        pass = candidateMomentum.DeltaR(isolationMomentum) <= fDeltaRMax && candidate->GetUniqueID() != isolation->GetUniqueID();
       }
 
@@ -233 +229 @@
         }
       }
-
     }
 
-   // find rho
+    // find rho
     rho = 0.0;
     if(fRhoInputArray)
     {
       fItRhoInputArray->Reset();
-      while((object = static_cast<Candidate*>(fItRhoInputArray->Next())))
+      while((object = static_cast<Candidate *>(fItRhoInputArray->Next())))
       {
         if(eta >= object->Edges[0] && eta < object->Edges[1])
@@ -250 +245 @@
     }
 
-
-
     // correct sum for pile-up contamination
-    sumDBeta = sumChargedNoPU + TMath::Max(sumNeutral - 0.5*sumChargedPU, 0.0);
-    sumRhoCorr = sumChargedNoPU + TMath::Max(sumNeutral - TMath::Max(rho, 0.0)*fDeltaRMax*fDeltaRMax*TMath::Pi(), 0.0);
-    ratioDBeta = sumDBeta/candidateMomentum.Pt();
-    ratioRhoCorr = sumRhoCorr/candidateMomentum.Pt();
+    sumDBeta = sumChargedNoPU + TMath::Max(sumNeutral - 0.5 * sumChargedPU, 0.0);
+    sumRhoCorr = sumChargedNoPU + TMath::Max(sumNeutral - TMath::Max(rho, 0.0) * fDeltaRMax * fDeltaRMax * TMath::Pi(), 0.0);
+    ratioDBeta = sumDBeta / candidateMomentum.Pt();
+    ratioRhoCorr = sumRhoCorr / candidateMomentum.Pt();
 
     candidate->IsolationVar = ratioDBeta;