#include "AliESDEvent.h"
#include "AliAODEvent.h"
#include "AliAODHandler.h"
+#include "AliMCEventHandler.h"
+#include "AliStack.h"
ClassImp(AliAnalysisTaskJets)
//
// Default AOD
if (fDebug > 1) printf("AnalysisTaskJets::CreateOutPutData() \n");
- AliAODHandler* handler = (AliAODHandler*) ((AliAnalysisManager::GetAnalysisManager())->GetEventHandler());
+ AliAODHandler* handler = (AliAODHandler*) ((AliAnalysisManager::GetAnalysisManager())->GetOutputEventHandler());
+
fAOD = handler->GetAOD();
fTreeA = handler->GetTree();
fJetFinder->ConnectAOD(fAOD);
// Histogram
OpenFile(1);
fHisto = new TH1F("fHisto", "Jet Et", 100, 0., 100.);
+
}
void AliAnalysisTaskJets::Init()
{
// Execute analysis for current event
//
+ AliMCEventHandler* mctruth = (AliMCEventHandler*)
+ ((AliAnalysisManager::GetAnalysisManager())->GetMCtruthEventHandler());
+ AliStack* stack = mctruth->Stack();
+ printf("AliAnalysisTaskJets: Number of tracks %5d\n", stack->GetNtrack());
+
Long64_t ientry = fChain->GetReadEntry();
if (fDebug > 1) printf("Analysing event # %5d\n", (Int_t) ientry);
fJetFinder->ProcessEvent(ientry);
// Create the outputs containers
OpenFile(0) ;
- AliAODHandler* handler = (AliAODHandler*) ((AliAnalysisManager::GetAnalysisManager())->GetEventHandler());
+ AliAODHandler* handler = (AliAODHandler*) ((AliAnalysisManager::GetAnalysisManager())->GetOutputEventHandler());
fTreeA = handler->GetTree() ;
fAOD = handler->GetAOD();
fAODPhotons = fAOD->GetClusters() ;