Be able to set recombination scheme (Ruediger)

diff --git a/PWGJE/EMCALJetTasks/AliEmcalJetTask.cxx b/PWGJE/EMCALJetTasks/AliEmcalJetTask.cxx
index 08753c23b4ff9f6ca9ba38cc8448abc4befdf60a..c0ef3b551d7595a1000a6f1afcdad84c8063a8b6 100644 (file)
--- a/PWGJE/EMCALJetTasks/AliEmcalJetTask.cxx
+++ b/PWGJE/EMCALJetTasks/AliEmcalJetTask.cxx
@@ -56,6 +56,7 @@ AliEmcalJetTask::AliEmcalJetTask() :
   fJetEtaMax(+1),
   fGhostArea(0.005),
   fMinMCLabel(0),
+  fRecombScheme(-1),
   fTrackEfficiency(1.),
   fIsInit(0),
   fIsPSelSet(0),
@@ -94,6 +95,7 @@ AliEmcalJetTask::AliEmcalJetTask(const char *name) :
   fJetEtaMax(+1),
   fGhostArea(0.005),
   fMinMCLabel(0),
+  fRecombScheme(-1),
   fTrackEfficiency(1.),
   fIsInit(0),
   fIsPSelSet(0),
@@ -178,7 +180,9 @@ void AliEmcalJetTask::FindJets()
   fjw.SetAreaType(fastjet::active_area_explicit_ghosts);
   fjw.SetGhostArea(fGhostArea);
   fjw.SetR(fRadius);
-  fjw.SetAlgorithm(jalgo);  
+  fjw.SetAlgorithm(jalgo);
+  if(fRecombScheme>0)
+    fjw.SetRecombScheme(static_cast<fastjet::RecombinationScheme>(fRecombScheme)); 
   fjw.SetMaxRap(fEtaMax);
   fjw.Clear();
 

diff --git a/PWGJE/EMCALJetTasks/AliEmcalJetTask.h b/PWGJE/EMCALJetTasks/AliEmcalJetTask.h
index b4a4968e46da8bef6c77c33ba6a2afc73a640442..88e8d74b21838cb0c0dfd1f1de84006960509995 100644 (file)
--- a/PWGJE/EMCALJetTasks/AliEmcalJetTask.h
+++ b/PWGJE/EMCALJetTasks/AliEmcalJetTask.h
@@ -62,6 +62,7 @@ class AliEmcalJetTask : public AliAnalysisTaskSE {
   void                   SetJetPhiRange(Double_t pmi, Double_t pma) {fJetPhiMin = pmi; fJetPhiMax = pma; }
   void                   SetGhostArea(Double_t gharea)    { fGhostArea      = gharea;  }
   void                   SetMinMCLabel(Int_t s)           { fMinMCLabel     = s     ;  }
+  void                   SetRecombSheme(Int_t scheme)     { fRecombScheme   = scheme;  }
   void                   SelectCollisionCandidates(UInt_t offlineTriggerMask = AliVEvent::kMB)
   {
     if(!fIsPSelSet)
@@ -105,6 +106,7 @@ class AliEmcalJetTask : public AliAnalysisTaskSE {
   Double_t               fJetEtaMax;              // maximum eta to keep jet in output
   Double_t               fGhostArea;              // ghost area
   Int_t                  fMinMCLabel;             // minimum MC label value for the tracks/clusters being considered MC particles
+  Int_t                  fRecombScheme;           // recombination scheme used by fastjet
   Double_t               fTrackEfficiency;        // artificial tracking inefficiency (0...1)
   Bool_t                 fIsInit;                 //!=true if already initialized
   Bool_t                 fIsPSelSet;              //!=true if physics selection was set