Changes in modules/FastJetFinder.cc [ba75867:f319c1d] in git

File:

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

modules/FastJetFinder.cc

@@ -66 +66 @@
 #include "fastjet/contribs/Nsubjettiness/ExtraRecombiners.hh"
 
+#include "fastjet/contribs/ValenciaPlugin/ValenciaPlugin.hh"
+
 #include "fastjet/tools/Filter.hh"
 #include "fastjet/tools/Pruner.hh"
@@ -79 +81 @@
 FastJetFinder::FastJetFinder() :
   fPlugin(0), fRecomb(0), fAxesDef(0), fMeasureDef(0), fNjettinessPlugin(0), 
-  fDefinition(0), fAreaDefinition(0), fItInputArray(0)
+  fDefinition(0), fAreaDefinition(0), fItInputArray(0), fValenciaPlugin(0)
 {
 
@@ -118 +120 @@
   fJetPTMin = GetDouble("JetPTMin", 10.0);
 
+ 
   //-- N(sub)jettiness parameters --
 
@@ -125 +128 @@
   fRcutOff = GetDouble("RcutOff", 0.8); // used only if Njettiness is used as jet clustering algo (case 8)
   fN = GetInt("N", 2);                  // used only if Njettiness is used as jet clustering algo (case 8)
-     
+
+  //-- Exclusive clustering for e+e- collisions --
+  
+  fNJets = GetInt("NJets",2);
+  fExclusiveClustering = GetBool("ExclusiveClustering", false);
+
+  //-- Valencia Linear Collider algorithm
+  fGamma = GetDouble("Gamma", 1.0);
+  //fBeta parameter see above
+  
   fMeasureDef = new NormalizedMeasure(fBeta, fParameterR);
    
@@ -236 +248 @@
       fDefinition = new JetDefinition(fNjettinessPlugin);
       break;
+  case 9:
+      fValenciaPlugin = new ValenciaPlugin(fParameterR, fBeta, fGamma);
+      fDefinition = new JetDefinition(fValenciaPlugin);
+      break;
+
   }
 
@@ -294 +311 @@
   if(fAxesDef) delete fAxesDef;
   if(fMeasureDef) delete fMeasureDef;
+  if(fValenciaPlugin) delete static_cast<JetDefinition::Plugin*>(fValenciaPlugin);
+
 }
 
@@ -357 +376 @@
 
   outputList.clear();
-  outputList = sorted_by_pt(sequence->inclusive_jets(fJetPTMin));
-
+
+  if(fExclusiveClustering)
+    {
+      outputList = sorted_by_pt(sequence->exclusive_jets( fNJets ));
+    }
+  else
+    {
+      outputList = sorted_by_pt(sequence->inclusive_jets(fJetPTMin));
+    }
 
   // loop over all jets and export them
@@ -494 +520 @@
       candidate->NSubJetsSoftDropped = softdrop_jet.pieces().size();
 
-      candidate->SoftDroppedJet = candidate->SoftDroppedP4[0];
-
       for (size_t i = 0; i < subjets.size()  and i < 4; i++)
       {
             if(subjets.at(i).pt() < 0) continue ; 
             candidate->SoftDroppedP4[i+1].SetPtEtaPhiM(subjets.at(i).pt(), subjets.at(i).eta(), subjets.at(i).phi(), subjets.at(i).m());
-            if(i==0) candidate->SoftDroppedSubJet1 = candidate->SoftDroppedP4[i+1];
-            if(i==1) candidate->SoftDroppedSubJet2 = candidate->SoftDroppedP4[i+1];
       }
     }