Small corrections to the TENDER code, which will support AliTenderInputHander(current...
[u/mrichter/AliRoot.git] / ANALYSIS / AliAnalysisTaskME.cxx
CommitLineData
b81460ad 1/**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
3 * *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
6 * *
7 * Permission to use, copy, modify and distribute this software and its *
8 * documentation strictly for non-commercial purposes is hereby granted *
9 * without fee, provided that the above copyright notice appears in all *
10 * copies and that both the copyright notice and this permission notice *
11 * appear in the supporting documentation. The authors make no claims *
12 * about the suitability of this software for any purpose. It is *
13 * provided "as is" without express or implied warranty. *
14 **************************************************************************/
15
16/* $Id$ */
17
18#include <TROOT.h>
19#include <TSystem.h>
20#include <TInterpreter.h>
21#include <TChain.h>
22#include <TFile.h>
23#include <TList.h>
24
25#include "AliAnalysisTaskME.h"
26#include "AliAnalysisManager.h"
79478008 27#include "AliAnalysisDataSlot.h"
b81460ad 28#include "AliAODEvent.h"
cd17c2ba 29#include "AliVEvent.h"
b81460ad 30#include "AliAODHandler.h"
cd17c2ba 31#include "AliMultiEventInputHandler.h"
b81460ad 32#include "AliLog.h"
33
34
df8176e9 35ClassImp(AliAnalysisTaskME)
b81460ad 36
37////////////////////////////////////////////////////////////////////////
38
39AliAnalysisTaskME::AliAnalysisTaskME():
40 AliAnalysisTask(),
41 fDebug(0),
42 fEntry(0),
43 fFreshBufferOnly(kFALSE),
44 fInputHandler(0x0),
45 fOutputAOD(0x0),
3c329509 46 fTreeA(0x0),
47 fOfflineTriggerMask(0)
b81460ad 48{
49 // Default constructor
50}
51
52AliAnalysisTaskME::AliAnalysisTaskME(const char* name):
53 AliAnalysisTask(name, "AnalysisTaskME"),
54 fDebug(0),
55 fEntry(0),
56 fFreshBufferOnly(kFALSE),
57 fInputHandler(0x0),
58 fOutputAOD(0x0),
3c329509 59 fTreeA(0x0),
60 fOfflineTriggerMask(0)
b81460ad 61{
62 // Default constructor
63 DefineInput (0, TChain::Class());
64 DefineOutput(0, TTree::Class());
65}
66
67AliAnalysisTaskME::AliAnalysisTaskME(const AliAnalysisTaskME& obj):
68 AliAnalysisTask(obj),
69 fDebug(0),
70 fEntry(0),
71 fFreshBufferOnly(kFALSE),
72 fInputHandler(0x0),
73 fOutputAOD(0x0),
3c329509 74 fTreeA(0x0),
75 fOfflineTriggerMask(0)
b81460ad 76{
77// Copy constructor
78 fDebug = obj.fDebug;
79 fEntry = obj.fEntry;
80 fInputHandler = obj.fInputHandler;
81 fOutputAOD = obj.fOutputAOD;
3c329509 82 fTreeA = obj.fTreeA;
83 fOfflineTriggerMask = obj.fOfflineTriggerMask;
b81460ad 84}
85
86
87AliAnalysisTaskME& AliAnalysisTaskME::operator=(const AliAnalysisTaskME& other)
88{
89// Assignment
90 AliAnalysisTask::operator=(other);
91 fDebug = other.fDebug;
92 fEntry = other.fEntry;
93 fFreshBufferOnly = other.fFreshBufferOnly;
94 fInputHandler = other.fInputHandler;
95 fOutputAOD = other.fOutputAOD;
96 fTreeA = other.fTreeA;
3c329509 97 fOfflineTriggerMask = other.fOfflineTriggerMask;
b81460ad 98 return *this;
99}
100
101
102void AliAnalysisTaskME::ConnectInputData(Option_t* /*option*/)
103{
104// Connect the input data
105 if (fDebug > 1) printf("AnalysisTaskME::ConnectInputData() \n");
106//
107// Multi AOD
108//
cd17c2ba 109 fInputHandler = dynamic_cast<AliMultiEventInputHandler*>
b81460ad 110 ((AliAnalysisManager::GetAnalysisManager())->GetInputEventHandler());
94d4e23c 111 if (fInputHandler == 0) {
cd17c2ba 112 AliFatal("Event Handler has to be MultiEventInputHandler !");
94d4e23c 113 } else {
114 // Check that we have an event pool
115 if (!fInputHandler->GetEventPool()) {
116 fInputHandler->SetEventPool(AliAnalysisManager::GetAnalysisManager()->GetEventPool());
117 if (!fInputHandler->GetEventPool())
cd17c2ba 118 AliFatal("MultiEventInputHandler has no EventPool connected !");
94d4e23c 119 }
120 }
b81460ad 121}
122
123void AliAnalysisTaskME::CreateOutputObjects()
124{
125// Create the output container
126//
127// Default AOD
128 if (fDebug > 1) printf("AnalysisTaskME::CreateOutPutData() \n");
129
130 AliAODHandler* handler = (AliAODHandler*)
131 ((AliAnalysisManager::GetAnalysisManager())->GetOutputEventHandler());
132
133 if (handler) {
134 fOutputAOD = handler->GetAOD();
135 fTreeA = handler->GetTree();
136 } else {
137 AliWarning("No AOD Event Handler connected.") ;
138 }
139 UserCreateOutputObjects();
140}
141
142void AliAnalysisTaskME::Exec(Option_t* option)
143{
144//
145// Exec analysis of one event
e7d9460a 146
b81460ad 147 if (fDebug > 1) AliInfo("AliAnalysisTaskME::Exec() \n");
148 if( fInputHandler )
149 fEntry = fInputHandler->GetReadEntry();
150 if ( !((Entry()-1)%100) && fDebug > 0)
151 AliInfo(Form("%s ----> Processing event # %lld", CurrentFileName(), Entry()));
15917707 152
153 AliAODHandler* outputHandler = (AliAODHandler*)
154 ((AliAnalysisManager::GetAnalysisManager())->GetOutputEventHandler());
3c329509 155//
156// Was event selected ? If no event selection mechanism, the event SHOULD be selected (AG)
157 UInt_t isSelected = AliVEvent::kAny;
9daedfd1 158 if(fInputHandler && fInputHandler->GetEventSelection()) {
3c329509 159 // Get the actual offline trigger mask for the event and AND it with the
160 // requested mask. If no mask requested select by default the event.
161 if (fOfflineTriggerMask)
162 isSelected = fOfflineTriggerMask & fInputHandler->IsEventSelected();
163 }
9daedfd1 164
3c329509 165 if (!isSelected) {
166 if (fDebug > 1) AliInfo("Event rejected \n");
167 fInputHandler->EventSkipped();
168 return;
169 }
b81460ad 170// Call the user analysis
3c329509 171
9daedfd1 172 if (fInputHandler && fInputHandler->IsBufferReady()) {
b81460ad 173 if ((fFreshBufferOnly && fInputHandler->IsFreshBuffer()) || !fFreshBufferOnly)
174 {
cd17c2ba 175 if (outputHandler) outputHandler->SetFillAOD(kTRUE);
b81460ad 176 UserExec(option);
79478008 177 // Added protection in case the derived task is not an AOD producer.
178 AliAnalysisDataSlot *out0 = GetOutputSlot(0);
179 if (out0 && out0->IsConnected()) PostData(0, fTreeA);
15917707 180 } else {
cd17c2ba 181 if (outputHandler) outputHandler->SetFillAOD(kFALSE);
b81460ad 182 }
e7d9460a 183 } else {
184 AliInfo(Form("Waiting for buffer to be ready !\n"));
b81460ad 185 }
186}
187
188const char* AliAnalysisTaskME::CurrentFileName()
189{
190// Returns the current file name
191 if(fInputHandler )
192 return fInputHandler->GetTree()->GetCurrentFile()->GetName();
193 else return "";
194}
195
c8fe2783 196void AliAnalysisTaskME::AddAODBranch(const char* cname, void* addobj, const char *fname)
b81460ad 197{
198 // Add a new branch to the aod tree
199 AliAODHandler* handler = (AliAODHandler*)
200 ((AliAnalysisManager::GetAnalysisManager())->GetOutputEventHandler());
201 if (handler) {
c8fe2783 202 handler->AddBranch(cname, addobj, fname);
b81460ad 203 }
204}
205
cd17c2ba 206AliVEvent* AliAnalysisTaskME::GetEvent(Int_t iev)
b81460ad 207{
208 // Get an event from the input handler
209 return (fInputHandler->GetEvent(iev));
210}
211