#include "AliAODInputHandler.h"
#include "AliMCEventHandler.h"
#include "AliInputEventHandler.h"
+#include "AliMultiInputEventHandler.h"
#include "AliESDInputHandler.h"
#include "AliMCEvent.h"
#include "AliStack.h"
fTreeA(0x0),
fCurrentRunNumber(-1),
fHistosQA(0x0),
- fOfflineTriggerMask(0)
+ fOfflineTriggerMask(0),
+ fMultiInputHandler(0),
+ fMCEventHandler(0)
{
// Default constructor
}
fTreeA(0x0),
fCurrentRunNumber(-1),
fHistosQA(0x0),
- fOfflineTriggerMask(0)
+ fOfflineTriggerMask(0),
+ fMultiInputHandler(0),
+ fMCEventHandler(0)
{
// Default constructor
DefineInput (0, TChain::Class());
fTreeA(0x0),
fCurrentRunNumber(-1),
fHistosQA(0x0),
- fOfflineTriggerMask(0)
+ fOfflineTriggerMask(0),
+ fMultiInputHandler(obj.fMultiInputHandler),
+ fMCEventHandler(obj.fMCEventHandler)
{
// Copy constructor
fDebug = obj.fDebug;
fCurrentRunNumber = other.fCurrentRunNumber;
fHistosQA = other.fHistosQA;
fOfflineTriggerMask = other.fOfflineTriggerMask;
+ fMultiInputHandler = other.fMultiInputHandler;
+ fMCEventHandler = other.fMCEventHandler;
return *this;
}
{
// Connect the input data
if (fDebug > 1) printf("AnalysisTaskSE::ConnectInputData() \n");
-//
-// ESD
-//
- fInputHandler = (AliInputEventHandler*)
- ((AliAnalysisManager::GetAnalysisManager())->GetInputEventHandler());
-//
-// Monte Carlo
-//
- AliMCEventHandler* mcH = 0;
- mcH = (AliMCEventHandler*) ((AliAnalysisManager::GetAnalysisManager())->GetMCtruthEventHandler());
- if (mcH) {
- fMCEvent = mcH->MCEvent();
- }
-
+
+ // Connect input handlers (multi input handler is handled)
+ ConnectMultiHandler();
+
if (fInputHandler) {
if ((fInputHandler->GetTree())->GetBranch("ESDfriend."))
fESDfriend = ((AliESDInputHandler*)fInputHandler)->GetESDfriend();
AliError("No Input Event Handler connected") ;
return ;
}
+ // Disconnect multi handler
+ DisconnectMultiHandler();
}
void AliAnalysisTaskSE::CreateOutputObjects()
{
//
// Exec analysis of one event
+
+ ConnectMultiHandler();
+
if ( fDebug >= 10)
printf("Task is active %5d\n", IsActive());
if (handler && aodH) {
fMCEvent = aodH->MCEvent();
Bool_t merging = aodH->GetMergeEvents();
-
- if (!(handler->IsStandard()) && !(handler->AODIsReplicated())) {
+ AliAODEvent* aod = dynamic_cast<AliAODEvent*>(InputEvent());
+
+ if (aod && !(handler->IsStandard()) && !(handler->AODIsReplicated())) {
if ((handler->NeedsHeaderReplication()) && (fgAODHeader))
{
// copy the contents by assigment
- *fgAODHeader = *(dynamic_cast<AliAODHeader*>(InputEvent()->GetHeader()));
+ *fgAODHeader = *(aod->GetHeader());
}
if ((handler->NeedsTracksBranchReplication() || merging) && (fgAODTracks))
{
- TClonesArray* tracks = (dynamic_cast<AliAODEvent*>(InputEvent()))->GetTracks();
+ TClonesArray* tracks = aod->GetTracks();
new (fgAODTracks) TClonesArray(*tracks);
}
if ((handler->NeedsVerticesBranchReplication() || merging) && (fgAODVertices))
{
- TClonesArray* vertices = (dynamic_cast<AliAODEvent*>(InputEvent()))->GetVertices();
+ TClonesArray* vertices = aod->GetVertices();
new (fgAODVertices) TClonesArray(*vertices);
}
if ((handler->NeedsV0sBranchReplication()) && (fgAODV0s))
{
- TClonesArray* v0s = (dynamic_cast<AliAODEvent*>(InputEvent()))->GetV0s();
+ TClonesArray* v0s = aod->GetV0s();
new (fgAODV0s) TClonesArray(*v0s);
}
if ((handler->NeedsTrackletsBranchReplication()) && (fgAODTracklets))
{
- *fgAODTracklets = *(dynamic_cast<AliAODEvent*>(InputEvent()))->GetTracklets();
+ *fgAODTracklets = *aod->GetTracklets();
}
if ((handler->NeedsPMDClustersBranchReplication()) && (fgAODPMDClusters))
{
- TClonesArray* pmdClusters = (dynamic_cast<AliAODEvent*>(InputEvent()))->GetPmdClusters();
+ TClonesArray* pmdClusters = aod->GetPmdClusters();
new (fgAODPMDClusters) TClonesArray(*pmdClusters);
}
if ((handler->NeedsJetsBranchReplication() || merging) && (fgAODJets))
{
- TClonesArray* jets = (dynamic_cast<AliAODEvent*>(InputEvent()))->GetJets();
+ TClonesArray* jets = aod->GetJets();
new (fgAODJets) TClonesArray(*jets);
}
if ((handler->NeedsFMDClustersBranchReplication()) && (fgAODFMDClusters))
{
- TClonesArray* fmdClusters = (dynamic_cast<AliAODEvent*>(InputEvent()))->GetFmdClusters();
+ TClonesArray* fmdClusters = aod->GetFmdClusters();
new (fgAODFMDClusters) TClonesArray(*fmdClusters);
}
if ((handler->NeedsCaloClustersBranchReplication() || merging) && (fgAODCaloClusters))
{
- TClonesArray* caloClusters = (dynamic_cast<AliAODEvent*>(InputEvent()))->GetCaloClusters();
+ TClonesArray* caloClusters = aod->GetCaloClusters();
new (fgAODCaloClusters) TClonesArray(*caloClusters);
}
if ((handler->NeedsMCParticlesBranchReplication() || merging) && (fgAODMCParticles))
{
- TClonesArray* mcParticles = (TClonesArray*) ((dynamic_cast<AliAODEvent*>(InputEvent()))->FindListObject("mcparticles"));
+ TClonesArray* mcParticles = (TClonesArray*) (aod->FindListObject("mcparticles"));
new (fgAODMCParticles) TClonesArray(*mcParticles);
}
}
// Call the user analysis
- AliMCEventHandler* mcH = 0;
- mcH = (AliMCEventHandler*) ((AliAnalysisManager::GetAnalysisManager())->GetMCtruthEventHandler());
-
- if (!mcH) {
+ if (!fMCEventHandler) {
if (isSelected)
UserExec(option);
} else {
- if (isSelected && (mcH->InitOk()))
+ if (isSelected && (fMCEventHandler->InitOk()))
UserExec(option);
}
// Added protection in case the derived task is not an AOD producer.
AliAnalysisDataSlot *out0 = GetOutputSlot(0);
if (out0 && out0->IsConnected()) PostData(0, fTreeA);
+
+ DisconnectMultiHandler();
}
const char* AliAnalysisTaskSE::CurrentFileName()
TObject *obj;
while ((obj=next())) mgr->LoadBranch(obj->GetName());
}
+
+
+//_________________________________________________________________________________________________
+void AliAnalysisTaskSE::ConnectMultiHandler()
+{
+ //
+ // Connect MultiHandler
+ //
+ fInputHandler = (AliInputEventHandler *)((AliAnalysisManager::GetAnalysisManager())->GetInputEventHandler());
+ fMultiInputHandler = dynamic_cast<AliMultiInputEventHandler *>(fInputHandler);
+ if (fMultiInputHandler) {
+ fInputHandler = dynamic_cast<AliInputEventHandler *>(fMultiInputHandler->GetFirstInputEventHandler());
+ fMCEventHandler = dynamic_cast<AliMCEventHandler *>(fMultiInputHandler->GetFirstMCEventHandler());
+ } else {
+ fMCEventHandler = dynamic_cast<AliMCEventHandler *>((AliAnalysisManager::GetAnalysisManager())->GetMCtruthEventHandler());
+ }
+ if (fMCEventHandler) fMCEvent = fMCEventHandler->MCEvent();
+}
+
+//_________________________________________________________________________________________________
+void AliAnalysisTaskSE::DisconnectMultiHandler()
+{
+ //
+ // Disconnect MultiHandler
+ //
+ if (fMultiInputHandler) fInputHandler = fMultiInputHandler;
+}