added midpoint (from Jiri)

diff --git a/PWGJE/EMCALJetTasks/AliFJWrapper.h b/PWGJE/EMCALJetTasks/AliFJWrapper.h
index e93fc88eb1fda447aff8342ee08c47ba96daa87c..1da98a02518787cdee1e4e19939a4947343017cf 100644 (file)
--- a/PWGJE/EMCALJetTasks/AliFJWrapper.h
+++ b/PWGJE/EMCALJetTasks/AliFJWrapper.h
@@ -390,6 +390,15 @@ Int_t AliFJWrapper::Run()
                                       0,    //search of stable cones - zero = until no more
                                       1.0); // this should be seed effectively for proto jets
       fJetDef = new fastjet::JetDefinition(fPlugin);
+    } else if (fPluginAlgor == 1) {
+      // CDF cone
+      // NOTE: hardcoded split parameter
+      Double_t overlap_threshold = 0.75; // NOTE: this actually splits a lot: thr/min(pt1,pt2)
+      fPlugin = new fj::CDFMidPointPlugin(fR, 
+                                      overlap_threshold,
+                                      1.0,    //search of stable cones - zero = until no more
+                                      1.0); // this should be seed effectively for proto jets
+      fJetDef = new fastjet::JetDefinition(fPlugin);
     } else {
       AliError("[e] Unrecognized plugin number!");
     }

diff --git a/PWGJE/EMCALJetTasks/FJ_includes.h b/PWGJE/EMCALJetTasks/FJ_includes.h
index d94ad9e3d2796280e822d6d88d1391996c276187..573c8bf0999c403e5e4afd6fe28b5c4cd9eb60ec 100644 (file)
--- a/PWGJE/EMCALJetTasks/FJ_includes.h
+++ b/PWGJE/EMCALJetTasks/FJ_includes.h
@@ -14,6 +14,7 @@
 #include <fastjet/ClusterSequenceArea.hh>
 #include <fastjet/AreaDefinition.hh>
 #include <fastjet/SISConePlugin.hh>
+#include <fastjet/CDFMidPointPlugin.hh>
 #ifdef FASTJET_VERSION
 #include <fastjet/tools/JetMedianBackgroundEstimator.hh>
 #endif