Minimum energy for jets can be now required (D. Perrino)
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 12 Mar 2008 20:47:36 +0000 (20:47 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 12 Mar 2008 20:47:36 +0000 (20:47 +0000)
JETAN/AliAnalysisTaskJets.cxx
JETAN/AliAnalysisTaskJets.h
JETAN/AliDAJetFinder.cxx

index d185bb4..4fe4e75 100644 (file)
@@ -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()
index 7ff8860..cbd9a02 100644 (file)
@@ -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
 };
  
index 38e8d1c..4e97585 100644 (file)
@@ -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 {