#include "AliAnalysisTaskSE.h"
#include "AliAnalysisManager.h"
+#include "AliAnalysisCuts.h"
#include "AliAnalysisDataSlot.h"
+
#include "AliESDEvent.h"
+#include "AliESDfriend.h"
#include "AliESD.h"
#include "AliAODEvent.h"
#include "AliAODHeader.h"
#include "AliAODInputHandler.h"
#include "AliMCEventHandler.h"
#include "AliInputEventHandler.h"
+#include "AliESDInputHandler.h"
#include "AliMCEvent.h"
#include "AliStack.h"
#include "AliLog.h"
+#include "AliAODDimuon.h"
ClassImp(AliAnalysisTaskSE)
////////////////////////////////////////////////////////////////////////
-AliAODHeader* AliAnalysisTaskSE::fgAODHeader = NULL;
-TClonesArray* AliAnalysisTaskSE::fgAODTracks = NULL;
-TClonesArray* AliAnalysisTaskSE::fgAODVertices = NULL;
-TClonesArray* AliAnalysisTaskSE::fgAODV0s = NULL;
-TClonesArray* AliAnalysisTaskSE::fgAODPMDClusters = NULL;
-TClonesArray* AliAnalysisTaskSE::fgAODJets = NULL;
-TClonesArray* AliAnalysisTaskSE::fgAODFMDClusters = NULL;
-TClonesArray* AliAnalysisTaskSE::fgAODCaloClusters = NULL;
-TClonesArray* AliAnalysisTaskSE::fgAODMCParticles = NULL;
-AliAODTracklets* AliAnalysisTaskSE::fgAODTracklets = NULL;
-AliAODCaloCells* AliAnalysisTaskSE::fgAODEmcalCells = NULL;
-AliAODCaloCells* AliAnalysisTaskSE::fgAODPhosCells = NULL;
-
+AliAODHeader* AliAnalysisTaskSE::fgAODHeader = NULL;
+TClonesArray* AliAnalysisTaskSE::fgAODTracks = NULL;
+TClonesArray* AliAnalysisTaskSE::fgAODVertices = NULL;
+TClonesArray* AliAnalysisTaskSE::fgAODV0s = NULL;
+TClonesArray* AliAnalysisTaskSE::fgAODPMDClusters = NULL;
+TClonesArray* AliAnalysisTaskSE::fgAODJets = NULL;
+TClonesArray* AliAnalysisTaskSE::fgAODFMDClusters = NULL;
+TClonesArray* AliAnalysisTaskSE::fgAODCaloClusters = NULL;
+TClonesArray* AliAnalysisTaskSE::fgAODMCParticles = NULL;
+AliAODTracklets* AliAnalysisTaskSE::fgAODTracklets = NULL;
+AliAODCaloCells* AliAnalysisTaskSE::fgAODEmcalCells = NULL;
+AliAODCaloCells* AliAnalysisTaskSE::fgAODPhosCells = NULL;
+TClonesArray* AliAnalysisTaskSE::fgAODDimuons = NULL;
AliAnalysisTaskSE::AliAnalysisTaskSE():
AliAnalysisTask(),
fDebug(0),
fEntry(0),
fInputEvent(0x0),
+ fESDfriend(0x0),
fInputHandler(0x0),
fOutputAOD(0x0),
fMCEvent(0x0),
- fTreeA(0x0)
+ fTreeA(0x0),
+ fCurrentRunNumber(-1),
+ fHistosQA(0x0),
+ fOfflineTriggerMask(0)
{
// Default constructor
}
fDebug(0),
fEntry(0),
fInputEvent(0x0),
+ fESDfriend(0x0),
fInputHandler(0x0),
fOutputAOD(0x0),
fMCEvent(0x0),
- fTreeA(0x0)
+ fTreeA(0x0),
+ fCurrentRunNumber(-1),
+ fHistosQA(0x0),
+ fOfflineTriggerMask(0)
{
// Default constructor
DefineInput (0, TChain::Class());
fDebug(0),
fEntry(0),
fInputEvent(0x0),
+ fESDfriend(0x0),
fInputHandler(0x0),
fOutputAOD(0x0),
fMCEvent(0x0),
- fTreeA(0x0)
+ fTreeA(0x0),
+ fCurrentRunNumber(-1),
+ fHistosQA(0x0),
+ fOfflineTriggerMask(0)
{
// Copy constructor
- fDebug = obj.fDebug;
- fEntry = obj.fEntry;
- fInputEvent = obj.fInputEvent;
- fInputHandler = obj.fInputHandler;
- fOutputAOD = obj.fOutputAOD;
- fMCEvent = obj.fMCEvent;
- fTreeA = obj.fTreeA;
+ fDebug = obj.fDebug;
+ fEntry = obj.fEntry;
+ fInputEvent = obj.fInputEvent;
+ fESDfriend = obj.fESDfriend;
+ fInputHandler = obj.fInputHandler;
+ fOutputAOD = obj.fOutputAOD;
+ fMCEvent = obj.fMCEvent;
+ fTreeA = obj.fTreeA;
+ fCurrentRunNumber = obj.fCurrentRunNumber;
+ fHistosQA = obj.fHistosQA;
+
}
{
// Assignment
AliAnalysisTask::operator=(other);
- fDebug = other.fDebug;
- fEntry = other.fEntry;
- fInputEvent = other.fInputEvent;
- fInputHandler = other.fInputHandler;
- fOutputAOD = other.fOutputAOD;
- fMCEvent = other.fMCEvent;
- fTreeA = other.fTreeA;
+ fDebug = other.fDebug;
+ fEntry = other.fEntry;
+ fInputEvent = other.fInputEvent;
+ fESDfriend = other.fESDfriend;
+ fInputHandler = other.fInputHandler;
+ fOutputAOD = other.fOutputAOD;
+ fMCEvent = other.fMCEvent;
+ fTreeA = other.fTreeA;
+ fCurrentRunNumber = other.fCurrentRunNumber;
+ fHistosQA = other.fHistosQA;
+ fOfflineTriggerMask = other.fOfflineTriggerMask;
return *this;
}
}
if (fInputHandler) {
- fInputEvent = fInputHandler->GetEvent();
+ if ((fInputHandler->GetTree())->GetBranch("ESDfriend."))
+ fESDfriend = ((AliESDInputHandler*)fInputHandler)->GetESDfriend();
+
+ fInputEvent = fInputHandler->GetEvent();
} else if( fMCEvent ) {
AliWarning("No Input Event Handler connected, only MC Truth Event Handler") ;
} else {
// Default AOD
if (fDebug > 1) printf("AnalysisTaskSE::CreateOutPutData() \n");
- AliAODHandler* handler = (AliAODHandler*)
+ AliAODHandler* handler = dynamic_cast<AliAODHandler*>
((AliAnalysisManager::GetAnalysisManager())->GetOutputEventHandler());
Bool_t merging = kFALSE;
if (aodIH->GetMergeEvents()) merging = kTRUE;
}
+
// Check if AOD replication has been required
if (handler) {
fOutputAOD = handler->GetAOD();
fTreeA = handler->GetTree();
- if (!(handler->IsStandard())) {
+ if (fOutputAOD && !(handler->IsStandard())) {
if ((handler->NeedsHeaderReplication()) && !(fgAODHeader))
{
if (fDebug > 1) AliInfo("Replicating header");
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();
}
{
//
// Exec analysis of one event
+ if ( fDebug >= 10)
+ printf("Task is active %5d\n", IsActive());
+
if (fDebug > 1) AliInfo("AliAnalysisTaskSE::Exec() \n");
+//
+ AliAODHandler* handler = dynamic_cast<AliAODHandler*>
+ ((AliAnalysisManager::GetAnalysisManager())->GetOutputEventHandler());
+
+ AliAODInputHandler* aodH = dynamic_cast<AliAODInputHandler*>(fInputHandler);
+//
+// Was event selected ?
+ UInt_t isSelected = 0;
+ if( fInputHandler && fInputHandler->GetEventSelection() && fOfflineTriggerMask) {
+ isSelected = fInputHandler->IsEventSelected() & fOfflineTriggerMask;
+ }
+
+ if (handler) handler->SetFillAOD(isSelected);
if( fInputHandler ) {
- fEntry = fInputHandler->GetReadEntry();
+ fEntry = fInputHandler->GetReadEntry();
+ fESDfriend = ((AliESDInputHandler*)fInputHandler)->GetESDfriend();
}
+
+// Notify the change of run number
+ if (InputEvent() && (InputEvent()->GetRunNumber() != fCurrentRunNumber)) {
+ fCurrentRunNumber = InputEvent()->GetRunNumber();
+ NotifyRun();
+ }
else if( fMCEvent )
fEntry = fMCEvent->Header()->GetEvent();
if ( !((Entry()-1)%100) && fDebug > 0)
AliInfo(Form("%s ----> Processing event # %lld", CurrentFileName(), Entry()));
- AliAODHandler* handler = (AliAODHandler*)
- ((AliAnalysisManager::GetAnalysisManager())->GetOutputEventHandler());
- AliAODInputHandler* aodH = dynamic_cast<AliAODInputHandler*>(fInputHandler);
+
+
if (handler && aodH) {
fMCEvent = aodH->MCEvent();
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[10];
+ for(Int_t imuon = 0; imuon < 10; 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
- UserExec(option);
- // Added protection in case the derived task is not an AOD producer.
+ AliMCEventHandler* mcH = 0;
+ mcH = (AliMCEventHandler*) ((AliAnalysisManager::GetAnalysisManager())->GetMCtruthEventHandler());
+
+ if (!mcH) {
+ if (!fOfflineTriggerMask || isSelected)
+ UserExec(option);
+ } else {
+ if ((!fOfflineTriggerMask || isSelected) && (mcH->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);
}
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 handler->IsStandard();
}
+
+Bool_t AliAnalysisTaskSE::Notify()
+{
+ return (UserNotify());
+}
+
+