Changeset ce4feac in git

Timestamp:

Oct 6, 2015, 10:18:59 AM (9 years ago)

Author:

Michele Selvaggi <michele.selvaggi@…>

Branches:

ImprovedOutputFile, Timing, dual_readout, llp, master

Children:

b3c0535

Parents:

4f00e0b

Message:

fill NNeutrals and NCharged in FastJetFinder

File:

• modules/FastJetFinder.cc (modified) (5 diffs)

modules/FastJetFinder.cc

@@ -282 +282 @@
   Double_t deta, dphi, detaMax, dphiMax;
   Double_t time, timeWeight;
-  Int_t number;
+  Int_t number, ncharged, nneutrals;
   Int_t charge; 
   Double_t rho = 0.0;
@@ -340 +340 @@
   detaMax = 0.0;
   dphiMax = 0.0;
+  
   for(itOutputList = outputList.begin(); itOutputList != outputList.end(); ++itOutputList)
   {
@@ -356 +357 @@
 
     charge = 0;
+
+    ncharged = 0;
+    nneutrals = 0;
 
     inputList.clear();
@@ -370 +374 @@
       if(dphi > dphiMax) dphiMax = dphi;
 
+      if(constituent->Charge == 0) nneutrals++;
+      else ncharged++;
+
       time += TMath::Sqrt(constituent->Momentum.E())*(constituent->Position.T());
       timeWeight += TMath::Sqrt(constituent->Momentum.E());
@@ -385 +392 @@
     candidate->DeltaPhi = dphiMax;
     candidate->Charge = charge; 
+    candidate->NNeutrals = nneutrals;
+    candidate->NCharged = ncharged;
+    
     //------------------------------------
     // Trimming