Be able to set recombination scheme (Ruediger)
authorloizides <loizides@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 1 Oct 2013 10:55:10 +0000 (10:55 +0000)
committerloizides <loizides@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 1 Oct 2013 10:55:10 +0000 (10:55 +0000)
PWGJE/EMCALJetTasks/AliEmcalJetTask.cxx
PWGJE/EMCALJetTasks/AliEmcalJetTask.h

index 08753c2..c0ef3b5 100644 (file)
@@ -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();
 
index b4a4968..88e8d74 100644 (file)
@@ -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