Changes in modules/TrackPileUpSubtractor.cc [341014c:781edf86] in git

File:

• modules/TrackPileUpSubtractor.cc (modified) (9 diffs)

modules/TrackPileUpSubtractor.cc

@@ -17 +17 @@
  */
 
+
 /** \class TrackPileUpSubtractor
  *
@@ -31 +32 @@
 #include "classes/DelphesFormula.h"
 
+#include "ExRootAnalysis/ExRootResult.h"
+#include "ExRootAnalysis/ExRootFilter.h"
 #include "ExRootAnalysis/ExRootClassifier.h"
-#include "ExRootAnalysis/ExRootFilter.h"
-#include "ExRootAnalysis/ExRootResult.h"
 
+#include "TMath.h"
+#include "TString.h"
+#include "TFormula.h"
+#include "TRandom3.h"
+#include "TObjArray.h"
 #include "TDatabasePDG.h"
-#include "TFormula.h"
 #include "TLorentzVector.h"
-#include "TMath.h"
-#include "TObjArray.h"
-#include "TRandom3.h"
-#include "TString.h"
 
 #include <algorithm>
+#include <stdexcept>
 #include <iostream>
 #include <sstream>
-#include <stdexcept>
 
 using namespace std;
@@ -53 +54 @@
 
 TrackPileUpSubtractor::TrackPileUpSubtractor() :
-  fFormula(0)
+fFormula(0) 
 {
   fFormula = new DelphesFormula;
@@ -87 +88 @@
 
   size = param.GetSize();
-  for(i = 0; i < size / 2; ++i)
+  for(i = 0; i < size/2; ++i)
   {
-    array = ImportArray(param[i * 2].GetString());
+    array = ImportArray(param[i*2].GetString());
     iterator = array->MakeIterator();
 
-    fInputMap[iterator] = ExportArray(param[i * 2 + 1].GetString());
+    fInputMap[iterator] = ExportArray(param[i*2 + 1].GetString());
   }
 }
@@ -100 +101 @@
 void TrackPileUpSubtractor::Finish()
 {
-  map<TIterator *, TObjArray *>::iterator itInputMap;
+  map< TIterator *, TObjArray * >::iterator itInputMap;
   TIterator *iterator;
 
@@ -118 +119 @@
 {
   Candidate *candidate, *particle;
-  map<TIterator *, TObjArray *>::iterator itInputMap;
+  map< TIterator *, TObjArray * >::iterator itInputMap;
   TIterator *iterator;
   TObjArray *array;
-  Double_t z, zvtx = 0;
+  Double_t z, zvtx=0;
   Double_t pt, eta, phi, e;
+
 
   // find z position of primary vertex
 
   fItVertexInputArray->Reset();
-  while((candidate = static_cast<Candidate *>(fItVertexInputArray->Next())))
+  while((candidate = static_cast<Candidate*>(fItVertexInputArray->Next())))
   {
     if(!candidate->IsPU)
@@ -144 +146 @@
     // loop over all candidates
     iterator->Reset();
-    while((candidate = static_cast<Candidate *>(iterator->Next())))
+    while((candidate = static_cast<Candidate*>(iterator->Next())))
     {
-      particle = static_cast<Candidate *>(candidate->GetCandidates()->At(0));
+      particle = static_cast<Candidate*>(candidate->GetCandidates()->At(0));
       const TLorentzVector &candidateMomentum = particle->Momentum;
 
@@ -153 +155 @@
       phi = candidateMomentum.Phi();
       e = candidateMomentum.E();
-
+      
       z = particle->Position.Z();
 
@@ -159 +161 @@
       // assume perfect pile-up subtraction for tracks outside fZVertexResolution
 
-      if(candidate->Charge != 0 && candidate->IsPU && TMath::Abs(z - zvtx) > fFormula->Eval(pt, eta, phi, e) * 1.0e3)
+      if(candidate->Charge !=0 && candidate->IsPU && TMath::Abs(z-zvtx) > fFormula->Eval(pt, eta, phi, e)* 1.0e3)
       {
         candidate->IsRecoPU = 1;