Changeset d244bc9 in git for modules/FastJetFinder.cc

Timestamp:

Dec 4, 2014, 3:25:06 PM (10 years ago)

Author:

Michele <michele.selvaggi@…>

Branches:

ImprovedOutputFile, Timing, dual_readout, llp, master

Children:

51f7063

Parents:

6fb1a5d

Message:

added exception in FastJetFinder when winner-take-all recombination scheme is used that prevents jets from having mass = 0

File:

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

modules/FastJetFinder.cc

@@ -223 +223 @@
   Candidate *candidate, *constituent;
   TLorentzVector momentum;
+  
+  TLorentzVector constmomentum;
+  
   Double_t deta, dphi, detaMax, dphiMax;
   Double_t time, weightTime, avTime;
@@ -286 +289 @@
   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();
@@ -297 +304 @@
     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());