#include "AliAODInputHandler.h"
#include "AliMCEventHandler.h"
#include "AliInputEventHandler.h"
+#include "AliMultiInputEventHandler.h"
#include "AliESDInputHandler.h"
#include "AliMCEvent.h"
#include "AliStack.h"
#include "AliLog.h"
+#include "AliAODDimuon.h"
ClassImp(AliAnalysisTaskSE)
AliAODTracklets* AliAnalysisTaskSE::fgAODTracklets = NULL;
AliAODCaloCells* AliAnalysisTaskSE::fgAODEmcalCells = NULL;
AliAODCaloCells* AliAnalysisTaskSE::fgAODPhosCells = NULL;
+TClonesArray* AliAnalysisTaskSE::fgAODDimuons = NULL;
AliAnalysisTaskSE::AliAnalysisTaskSE():
AliAnalysisTask(),
fTreeA(0x0),
fCurrentRunNumber(-1),
fHistosQA(0x0),
- fSelectCollisions(0)
+ fOfflineTriggerMask(0),
+ fMultiInputHandler(0),
+ fMCEventHandler(0)
{
// Default constructor
}
fTreeA(0x0),
fCurrentRunNumber(-1),
fHistosQA(0x0),
- fSelectCollisions(0)
+ fOfflineTriggerMask(0),
+ fMultiInputHandler(0),
+ fMCEventHandler(0)
{
// Default constructor
DefineInput (0, TChain::Class());
fTreeA(0x0),
fCurrentRunNumber(-1),
fHistosQA(0x0),
- fSelectCollisions(0)
+ fOfflineTriggerMask(0),
+ fMultiInputHandler(obj.fMultiInputHandler),
+ fMCEventHandler(obj.fMCEventHandler)
{
// Copy constructor
fDebug = obj.fDebug;
fTreeA = other.fTreeA;
fCurrentRunNumber = other.fCurrentRunNumber;
fHistosQA = other.fHistosQA;
- fSelectCollisions = other.fSelectCollisions;
+ 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()
// Default AOD
if (fDebug > 1) printf("AnalysisTaskSE::CreateOutPutData() \n");
- AliAODHandler* handler = (AliAODHandler*)
+ AliAODHandler* handler = dynamic_cast<AliAODHandler*>
((AliAnalysisManager::GetAnalysisManager())->GetOutputEventHandler());
Bool_t merging = kFALSE;
fgAODMCParticles->SetName("mcparticles");
handler->AddBranch("TClonesArray", &fgAODMCParticles);
}
+ if ((handler->NeedsDimuonsBranchReplication() || merging) && !(fgAODDimuons))
+ {
+ if (fDebug > 1) AliInfo("Replicating dimuon branch\n");
+ fgAODDimuons = new TClonesArray("AliAODDimuon",0);
+ fgAODDimuons->SetName("dimuons");
+ handler->AddBranch("TClonesArray", &fgAODDimuons);
+ }
+
// cache the pointerd in the AODEvent
fOutputAOD->GetStdContent();
}
- } else {
- AliWarning("No AOD Event Handler connected.") ;
}
UserCreateOutputObjects();
}
{
//
// Exec analysis of one event
+
+ ConnectMultiHandler();
+
+ if ( fDebug >= 10)
+ printf("Task is active %5d\n", IsActive());
+
if (fDebug > 1) AliInfo("AliAnalysisTaskSE::Exec() \n");
//
- AliAODHandler* handler = (AliAODHandler*)
+ AliAODHandler* handler = dynamic_cast<AliAODHandler*>
((AliAnalysisManager::GetAnalysisManager())->GetOutputEventHandler());
+
AliAODInputHandler* aodH = dynamic_cast<AliAODInputHandler*>(fInputHandler);
//
-// Was event selected ?
- Bool_t isSelected = kTRUE;
- if( fInputHandler && fInputHandler->GetEventSelection() && fSelectCollisions) {
- isSelected = fInputHandler->IsEventSelected();
+// Was event selected ? If no event selection mechanism, the event SHOULD be selected (AG)
+ UInt_t isSelected = AliVEvent::kAny;
+ if( fInputHandler && fInputHandler->GetEventSelection()) {
+ // Get the actual offline trigger mask for the event and AND it with the
+ // requested mask. If no mask requested select by default the event.
+ if (fOfflineTriggerMask)
+ isSelected = fOfflineTriggerMask & fInputHandler->IsEventSelected();
}
+// Functionality below moved in the filter tasks (AG)
+// if (handler) handler->SetFillAOD(isSelected);
- if (handler) handler->SetFillAOD(isSelected);
-
- if( fInputHandler ) fEntry = fInputHandler->GetReadEntry();
-
+ if( fInputHandler ) {
+ fEntry = fInputHandler->GetReadEntry();
+ fESDfriend = ((AliESDInputHandler*)fInputHandler)->GetESDfriend();
+ }
+
// Notify the change of run number
- if (InputEvent()->GetRunNumber() != fCurrentRunNumber) {
+ if (InputEvent() && (InputEvent()->GetRunNumber() != fCurrentRunNumber)) {
fCurrentRunNumber = InputEvent()->GetRunNumber();
NotifyRun();
}
TClonesArray* mcParticles = (TClonesArray*) ((dynamic_cast<AliAODEvent*>(InputEvent()))->FindListObject("mcparticles"));
new (fgAODMCParticles) TClonesArray(*mcParticles);
}
-
+
+ if ((handler->NeedsDimuonsBranchReplication() || merging) && (fgAODDimuons))
+ {
+ fgAODDimuons->Clear();
+ TClonesArray& dimuons = *fgAODDimuons;
+ TClonesArray& tracksnew = *fgAODTracks;
+
+ Int_t nMuonTrack[100];
+ for(Int_t imuon = 0; imuon < 100; imuon++) nMuonTrack[imuon] = 0;
+ Int_t nMuons=0;
+ for(Int_t ii=0; ii < fgAODTracks->GetEntries(); ii++){
+ AliAODTrack *track = (AliAODTrack*) fgAODTracks->At(ii);
+ if(track->IsMuonTrack()) {
+ nMuonTrack[nMuons]= ii;
+ nMuons++;
+ }
+ }
+ Int_t jDimuons=0;
+ if(nMuons >= 2){
+ for(Int_t i = 0; i < nMuons; i++){
+ Int_t index0 = nMuonTrack[i];
+ for(Int_t j = i+1; j < nMuons; j++){
+ Int_t index1 = nMuonTrack[j];
+ tracksnew.At(index0)->ResetBit(kIsReferenced);
+ tracksnew.At(index0)->SetUniqueID(0);
+ tracksnew.At(index1)->ResetBit(kIsReferenced);
+ tracksnew.At(index1)->SetUniqueID(0);
+ new(dimuons[jDimuons++]) AliAODDimuon(tracksnew.At(index0),tracksnew.At(index1));
+ }
+ }
+ }
+ }
+
// Additional merging if needed
if (merging) {
// mcParticles
for (Int_t i = 0; i < ntr; i++) {
AliAODTrack* track = (AliAODTrack*) tracks->At(i);
AliAODTrack* newtrack = new((*fgAODTracks)[nc++]) AliAODTrack(*track);
-
newtrack->SetLabel(newtrack->GetLabel() + nc0);
}
}
// Call the user analysis
- if (!fSelectCollisions || isSelected)
- UserExec(option);
+ if (!fMCEventHandler) {
+ if (isSelected)
+ UserExec(option);
+ } else {
+ 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()
void AliAnalysisTaskSE::AddAODBranch(const char* cname, void* addobj, const char *fname)
{
// Add a new branch to the aod tree
- AliAODHandler* handler = (AliAODHandler*)
+ AliAODHandler* handler = dynamic_cast<AliAODHandler*>
((AliAnalysisManager::GetAnalysisManager())->GetOutputEventHandler());
if (handler) {
handler->AddBranch(cname, addobj, fname);
{
// Check if the output AOD handler is configured for standard or delta AOD.
// Users should first check that AODEvent() returns non-null.
- AliAODHandler* handler = (AliAODHandler*)
+ AliAODHandler* handler = dynamic_cast<AliAODHandler*>
((AliAnalysisManager::GetAnalysisManager())->GetOutputEventHandler());
if (!handler) {
Error("IsStandardAOD", "No AOD handler. Please use AODEvent() to check this first");
return (UserNotify());
}
+const AliEventTag *AliAnalysisTaskSE::EventTag() const
+{
+// Returns tag for the current event, if any. The return value should always be checked by the user.
+ if (!fInputHandler) {
+ Error("EventTag", "Input handler not yet available. Call this in UserExec");
+ return NULL;
+ }
+ return fInputHandler->GetEventTag();
+}
+
+void AliAnalysisTaskSE::LoadBranches() const
+{
+// Load all branches declared in fBranchNames data member of the parent class.
+// Should be called in UserExec.
+ if (!fInputHandler) {
+ Error("LoadBranches", "Input handler not available yet. Call this in UserExec");
+ return;
+ }
+ AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
+ if (mgr->GetAutoBranchLoading()) return;
+ TString taskbranches;
+ GetBranches(fInputHandler->GetDataType(), taskbranches);
+ if (taskbranches.IsNull()) return;
+ TObjArray *arr = taskbranches.Tokenize(",");
+ TIter next(arr);
+ 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;
+}
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff