Changes in modules/FastJetFinder.cc [7c0fcd5:01f457a] in git

File:

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

modules/FastJetFinder.cc

@@ -241 +241 @@
   Candidate *candidate, *constituent;
   TLorentzVector momentum;
+  
+  TLorentzVector constmomentum;
+  
   Double_t deta, dphi, detaMax, dphiMax;
   Double_t time, weightTime, avTime;
@@ -304 +307 @@
   for(itOutputList = outputList.begin(); itOutputList != outputList.end(); ++itOutputList)
   {
-    momentum.SetPxPyPzE(itOutputList->px(), itOutputList->py(), itOutputList->pz(), itOutputList->E());
+    jet = *itOutputList;
+    if(fJetAlgorithm == 7) jet = join(jet.constituents());
+    
+    momentum.SetPxPyPzE(jet.px(), jet.py(), jet.pz(), jet.E());
+    
     area.reset(0.0, 0.0, 0.0, 0.0);
     if(fAreaDefinition) area = itOutputList->area_4vector();
@@ -315 +322 @@
     inputList.clear();
     inputList = sequence->constituents(*itOutputList);
+
+    constmomentum.SetPxPyPzE(0.0,0.0,0.0,0.0);;
 
     for(itInputList = inputList.begin(); itInputList != inputList.end(); ++itInputList)
     {
       constituent = static_cast<Candidate*>(fInputArray->At(itInputList->user_index()));
+
+      constmomentum += constituent->Momentum;
 
       deta = TMath::Abs(momentum.Eta() - constituent->Momentum.Eta());