From: kleinb <kleinb@f7af4fe6-9843-0410-8265-dc069ae4e863>
Date: Wed, 22 Sep 2010 08:39:07 +0000 (+0000)
Subject: correcting jet area for fastjet, need definition of ghost area
X-Git-Url: http://git.uio.no/git/?a=commitdiff_plain;h=4dd6159dc7c5f4064f8e21dc0baef1af6e8b027b;p=u%2Fmrichter%2FAliRoot.git

correcting jet area for fastjet, need definition of ghost area
---

diff --git a/JETAN/AliAnalysisTaskJetCluster.cxx b/JETAN/AliAnalysisTaskJetCluster.cxx
index 58674864e79..36cf719c03c 100644
--- a/JETAN/AliAnalysisTaskJetCluster.cxx
+++ b/JETAN/AliAnalysisTaskJetCluster.cxx
@@ -628,9 +628,20 @@ void AliAnalysisTaskJetCluster::UserExec(Option_t */*option*/)
   }
   
   // run fast jet
+  // employ setters for these...
+
+  double ghostEtamax = 0.9;
+  double ghostArea   = 0.01;
+  int    activeAreaRepeats = 1;
+  // now create the object that holds info about ghosts                         
+  fastjet::GhostedAreaSpec ghostSpec(ghostEtamax, activeAreaRepeats, ghostArea)\
+;
+  fastjet::AreaType areaType =   fastjet::active_area;
+  fastjet::AreaDefinition areaDef = fastjet::AreaDefinition(areaType,ghostSpec);
+  
   fastjet::JetDefinition jetDef(fAlgorithm, fRparam, fRecombScheme, fStrategy);
   vector <fastjet::PseudoJet> inclusiveJets, sortedJets;
-  fastjet::ClusterSequence clustSeq(inputParticlesRec, jetDef);
+  fastjet::ClusterSequenceArea clustSeq(inputParticlesRec, jetDef,areaDef);
   
   inclusiveJets = clustSeq.inclusive_jets();
   sortedJets    = sorted_by_pt(inclusiveJets);
@@ -703,6 +714,7 @@ void AliAnalysisTaskJetCluster::UserExec(Option_t */*option*/)
      if(tmpPt>fJetOutputMinPt&&jarray){
        aodOutJet =  new ((*jarray)[nAodOutJets++]) AliAODJet(tmpRec);
        Double_t area=clustSeq.area(sortedJets[j]);
+       
        aodOutJet->SetEffArea(area,0);
      }