Minimum energy for jets can be now required (D. Perrino)

diff --git a/JETAN/AliAnalysisTaskJets.cxx b/JETAN/AliAnalysisTaskJets.cxx
index d185bb47878078bdb3f3e5e5b67c33d995745b3a..4fe4e75e0e8e59a721636520f44bb58bac650e42 100644 (file)
--- a/JETAN/AliAnalysisTaskJets.cxx
+++ b/JETAN/AliAnalysisTaskJets.cxx
@@ -29,6 +29,7 @@
 #include "AliESDEvent.h"
 #include "AliESD.h"
 #include "AliAODEvent.h"
+#include "AliAODJet.h"
 #include "AliAODHandler.h"
 #include "AliMCEventHandler.h"
 #include "AliESDInputHandler.h"
@@ -75,6 +76,9 @@ void AliAnalysisTaskJets::UserCreateOutputObjects()
     fListOfHistos = new TList();
     fHistos       = new AliJetHistos();
     fHistos->AddHistosToList(fListOfHistos);
+
+    fNewJets      = new TClonesArray("AliAODJet", 1000);
+//    AliAnalysisTaskSE::AddAODBranch("NewJets", "TClonesArray", &fNewJets);
 }
 
 void AliAnalysisTaskJets::Init()

diff --git a/JETAN/AliAnalysisTaskJets.h b/JETAN/AliAnalysisTaskJets.h
index 7ff88605044f1a8d48a3b651e9d2404f5e8c16c8..cbd9a02e9f015364b5e755b7bc80b984001c1260 100644 (file)
--- a/JETAN/AliAnalysisTaskJets.h
+++ b/JETAN/AliAnalysisTaskJets.h
@@ -10,6 +10,7 @@ class AliESDEvent;
 class TTree;
 class AliAODEvent;
 class AliJetHistos;
+class TClonesArray;
 
 
 class AliAnalysisTaskJets : public AliAnalysisTaskSE
@@ -24,17 +25,17 @@ class AliAnalysisTaskJets : public AliAnalysisTaskSE
     virtual void LocalInit() {Init();}
     virtual void UserExec(Option_t *option);
     virtual void Terminate(Option_t *option);
-    virtual void SetDebugLevel(Int_t level) {fDebug = level;}
  private:
   AliAnalysisTaskJets(const AliAnalysisTaskJets &det);
   AliAnalysisTaskJets &operator=(const AliAnalysisTaskJets &det);
     
  private:
-    Int_t         fDebug;        //  Debug flag
     AliJetFinder* fJetFinder;    //  Pointer to the jet finder 
     TTree*        fTree;         //! The input tree
     AliJetHistos* fHistos;       //  Histogram manager class
     TList*        fListOfHistos; //  Output list of histograms
+    TClonesArray* fNewJets;
+    
     ClassDef(AliAnalysisTaskJets, 2); // Analysis task for standard jet analysis
 };
  

diff --git a/JETAN/AliDAJetFinder.cxx b/JETAN/AliDAJetFinder.cxx
index 38e8d1c3b49c369f53ac4df633689e9ac4de1a83..4e9758546ae6dd7ab8f7182eaef018a668e92b74 100644 (file)
--- a/JETAN/AliDAJetFinder.cxx
+++ b/JETAN/AliDAJetFinder.cxx
@@ -426,6 +426,10 @@ void AliDAJetFinder::StoreJets(Int_t nk,Int_t *xx)
                                }
                                etDensMed/=extSurf;
                                etNoBg[iClust]=(*fY)(3,iClust)-etDensMed*surf[iClust];
+                               if (etNoBg[iClust]<((AliDAJetHeader*)fHeader)->GetEtMin()){
+                                       isJet[iClust]=kFALSE;
+                                       iClust=-1;
+                               }
                        }
                }
        } else {