Adding the possibility to filter for events containing a jet above threshold set...

diff --git a/JETAN/AliAnalysisTaskJets.cxx b/JETAN/AliAnalysisTaskJets.cxx
index 8bb2fb5af6eba8f35e2ae6097a469671b1725dde..662027616bce7a0645a408878aa369f76a65cd31 100644 (file)
--- a/JETAN/AliAnalysisTaskJets.cxx
+++ b/JETAN/AliAnalysisTaskJets.cxx
@@ -56,7 +56,8 @@ AliAnalysisTaskJets::AliAnalysisTaskJets():
   fListOfHistos(0x0),
   fChain(0x0),
   fOpt(0),
-  fReadAODFromOutput(0)
+  fReadAODFromOutput(0),
+  fFilterPt(0.)
 {
   // Default constructor
 }
@@ -72,7 +73,8 @@ AliAnalysisTaskJets::AliAnalysisTaskJets(const char* name):
   fListOfHistos(0x0),
   fChain(0x0),
   fOpt(0),
-  fReadAODFromOutput(0)
+  fReadAODFromOutput(0),
+  fFilterPt(0.)
 {
   // Default constructor
   DefineOutput(1, TList::Class());
@@ -89,7 +91,8 @@ AliAnalysisTaskJets::AliAnalysisTaskJets(const char* name, TChain* chain):
   fListOfHistos(0x0),
   fChain(chain),
   fOpt(0),
-  fReadAODFromOutput(0)
+  fReadAODFromOutput(0),
+  fFilterPt(0.)
 {
   // Default constructor
   DefineOutput(1, TList::Class());
@@ -219,6 +222,9 @@ void AliAnalysisTaskJets::UserExec(Option_t */*option*/)
   TClonesArray* jarray = 0;
   AliAODJetEventBackground*  evBkg = 0;
 
+  // only need this once
+  static AliAODHandler *aodH = dynamic_cast<AliAODHandler*>(AliAnalysisManager::GetAnalysisManager()->GetOutputEventHandler());
+
   if(fNonStdBranch.Length()==0) {
     jarray = AODEvent()->GetJets();
     evBkg = (AliAODJetEventBackground*)(AODEvent()->FindListObject(AliAODJetEventBackground::StdBranchName()));
@@ -249,6 +255,17 @@ void AliAnalysisTaskJets::UserExec(Option_t */*option*/)
   // Fill control histos
   if(jarray)fHistos->FillHistos(jarray);
 
+
+  if(aodH&&fFilterPt>0){
+    if(jarray->GetEntries()>0){
+      AliAODJet *jet = (AliAODJet*)jarray->At(0);
+      if(jet->Pt()>fFilterPt){
+       //      aodH->EnableFillAODforEvent();
+       aodH->SetFillAOD(kTRUE);
+      }
+    }
+  }
+
   // Post the data
   PostData(1, fListOfHistos);
   return;

diff --git a/JETAN/AliAnalysisTaskJets.h b/JETAN/AliAnalysisTaskJets.h
index d8f2e1161f03cf4cab9c618b0f80b3552e77b62e..6c9e0e2b0ffb95c9d7493b7f9932c3645a7d1833 100644 (file)
--- a/JETAN/AliAnalysisTaskJets.h
+++ b/JETAN/AliAnalysisTaskJets.h
@@ -36,7 +36,7 @@ class AliAnalysisTaskJets : public AliAnalysisTaskSE
     virtual void SetNonStdOutputFile(const char *c){fNonStdFile = c;}
     virtual void Terminate(Option_t *option);
     virtual void ReadAODFromOutput() {fReadAODFromOutput = kTRUE;}
-
+    virtual void SetFilterPt(Float_t f){fFilterPt = f;}
     
  private:
   AliAnalysisTaskJets(const AliAnalysisTaskJets &det);
@@ -53,8 +53,8 @@ class AliAnalysisTaskJets : public AliAnalysisTaskSE
   TChain*       fChain;             //  Chain 
   Int_t         fOpt;               //  Detector configuration used
   Bool_t        fReadAODFromOutput; //  Force reading of the AOD from the output
-  
-  ClassDef(AliAnalysisTaskJets, 4); // Analysis task for standard jet analysis
+  Float_t       fFilterPt;          // use this as a switch for writing the AOD, minium p_T of leading jet   
+  ClassDef(AliAnalysisTaskJets, 5); // Analysis task for standard jet analysis
 };
  
 #endif

diff --git a/PWG4/macros/AddTaskJets.C b/PWG4/macros/AddTaskJets.C
index d0ed23fd630accd9d1d223ea2113d988f8ae64fd..a4e60be02037bf952bc9c6b451dacf74b9cb41f4 100644 (file)
--- a/PWG4/macros/AddTaskJets.C
+++ b/PWG4/macros/AddTaskJets.C
@@ -160,8 +160,14 @@ AliAnalysisTaskJets *AddTaskJets(Char_t *jr, Char_t *jf, Float_t radius,UInt_t f
    jetana->SetJetFinder(jetFinder);\r
    jetana->SetConfigFile("");\r
    jetana->SetDebugLevel(2);\r
+   if(TMath::Abs((radius-0.4))< 0.02&&c_jf.Contains("fastjet")){\r
+     jetana->SetFilterPt(10.);\r
+   }\r
+\r
+\r
    mgr->AddTask(jetana);\r
 \r
+\r
    // Create ONLY the output containers for the data produced by the task.\r
    // Get and connect other common input/output containers via the manager as below\r
    //==============================================================================\r