]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - JETAN/AliAnalysisTaskJets.cxx
bug fix for outer sectors, extra interface
[u/mrichter/AliRoot.git] / JETAN / AliAnalysisTaskJets.cxx
index 936b91e38c78ffe5543e1fd1ea5ce7ddfd4e8474..d92731c9b3116cc9dc39ccdada52301bededf48b 100644 (file)
@@ -35,6 +35,7 @@
 #include "AliAODEvent.h"
 #include "AliAODJetEventBackground.h"
 #include "AliAODHandler.h"
+#include "AliAODExtension.h"
 #include "AliMCEventHandler.h"
 #include "AliESDInputHandler.h"
 #include "AliMCEvent.h"
@@ -181,6 +182,7 @@ void AliAnalysisTaskJets::UserCreateOutputObjects()
 
   // Histograms
   fListOfHistos = new TList();
+  fListOfHistos->SetOwner();
   fHistos       = new AliJetHistos();
   fHistos->AddHistosToList(fListOfHistos);
   
@@ -202,6 +204,10 @@ void AliAnalysisTaskJets::UserCreateOutputObjects()
   
   if(!fAODExtension)OutputTree()->GetUserInfo()->Add(fH);
   else fAODExtension->GetTree()->GetUserInfo()->Add(fH);
+
+  // post
+  PostData(1, fListOfHistos);
+
 }
 
 //----------------------------------------------------------------
@@ -218,6 +224,12 @@ void AliAnalysisTaskJets::Init()
   AliJetReaderHeader *header = (AliJetReaderHeader*)fJetFinder->GetReader()->GetReaderHeader();
   fOpt = header->GetDetector();
 
+// AODB path for Jetan Analysis... if not alread defined, use the standard one
+if ((!(((AliJetReader*) fJetFinder->GetReader())->GetJetanOADBPath()).Length())) {
+((AliJetReader*) fJetFinder->GetReader())->SetJetanOADBPath(AliAnalysisManager::GetOADBPath());
+}
+else Info( "Init"," OADBPath for Jetan was already defined (config file?) as %s",(char*)((fJetFinder->GetReader()->GetJetanOADBPath()).Data()));
+       
   // Initialise Jet Analysis
   if(fOpt == 0) fJetFinder->Init();
   else fJetFinder->InitTask(fChain); // V2
@@ -271,7 +283,7 @@ void AliAnalysisTaskJets::UserExec(Option_t */*option*/)
   if(jarray)fHistos->FillHistos(jarray);
 
 
-  if(aodH&&fFilterPt>0){
+  if(jarray&&aodH&&fFilterPt>0){
     if(jarray->GetEntries()>0){
       AliAODJet *jet = (AliAODJet*)jarray->At(0);
       if(jet->Pt()>fFilterPt){