Changes in / [b9e210d:723b77d] in git

Location:

modules

Files:

• Isolation.cc (modified) (2 diffs)
• TrackPileUpSubtractor.cc (modified) (4 diffs)
• TrackPileUpSubtractor.h (modified) (1 diff)

modules/Isolation.cc

@@ -206 +206 @@
       {
         sumAllParticles += isolationMomentum.Pt();
-        if(isolation->Charge != 0)
-        {
-          if(isolation->IsRecoPU)
-          {
-            sumChargedPU += isolationMomentum.Pt();
-          }
-          else
-          {
-            sumCharged += isolationMomentum.Pt();
-          }
+        if(isolation->Charge !=0)
+        {
+          sumCharged += isolationMomentum.Pt();
+          if(isolation->IsRecoPU != 0) sumChargedPU += isolationMomentum.Pt();
         }
         else
@@ -240 +234 @@
 
     // correct sum for pile-up contamination
-    sumDBeta = sumCharged + TMath::Max(sumNeutral - 0.5*sumChargedPU, 0.0);
-    sumRhoCorr = sumCharged + TMath::Max(sumNeutral - TMath::Max(rho, 0.0)*fDeltaRMax*fDeltaRMax*TMath::Pi(), 0.0);
+    sumDBeta = sumCharged + TMath::Max(sumNeutral-0.5*sumChargedPU, 0.0);
+    sumRhoCorr = sumCharged + TMath::Max(sumNeutral-TMath::Max(rho, 0.0)*fDeltaRMax*fDeltaRMax*TMath::Pi(), 0.0);
     ratioDBeta = sumDBeta/candidateMomentum.Pt();
     ratioRhoCorr = sumRhoCorr/candidateMomentum.Pt();

modules/TrackPileUpSubtractor.cc

@@ -67 +67 @@
 void TrackPileUpSubtractor::Init()
 {
-  // import input array
+// import input array
 
   fVertexInputArray = ImportArray(GetString("VertexInputArray", "PileUpMerger/vertices"));
   fItVertexInputArray = fVertexInputArray->MakeIterator();
-
+  
   fZVertexResolution  = GetDouble("ZVertexResolution", 0.005)*1.0E3;
 
   fPTMin = GetDouble("PTMin", 0.);
-
   // import arrays with output from other modules
-
+   
   ExRootConfParam param = GetParam("InputArray");
   Long_t i, size;
@@ -120 +119 @@
   Double_t z, zvtx=0;
 
-
+  
   // find z position of primary vertex
-
+  
   fItVertexInputArray->Reset();
   while((candidate = static_cast<Candidate*>(fItVertexInputArray->Next())))
@@ -128 +127 @@
     if(!candidate->IsPU)
     {
-      zvtx = candidate->Position.Z();
-      // break;
+    zvtx = candidate->Position.Z();
+    // break;
     }
   }
@@ -145 +144 @@
       particle = static_cast<Candidate*>(candidate->GetCandidates()->At(0));
       z = particle->Position.Z();
-
+      
       // apply pile-up subtraction
       // assume perfect pile-up subtraction for tracks outside fZVertexResolution
-
-      if(candidate->IsPU && TMath::Abs(z-zvtx) > fZVertexResolution)
+      
+      if(candidate->IsPU && TMath::Abs(z-zvtx) > fZVertexResolution) candidate->IsRecoPU = 1;
+      else 
       {
-        candidate->IsRecoPU = 1;
-      }
-      else
-      {
-        candidate->IsRecoPU = 0;
-        if(candidate->Momentum.Pt() > fPTMin) array->Add(candidate);
+         candidate->IsRecoPU = 0;
+         if( candidate->Momentum.Pt() > fPTMin) array->Add(candidate);
       }
     }

modules/TrackPileUpSubtractor.h

@@ -50 +50 @@
   Double_t fZVertexResolution;
 
-  Double_t fPTMin;
+  Double_t fPTMin; 
 
   std::map< TIterator *, TObjArray * > fInputMap; //!