3 #include "AliPhysicsSelectionTask.h"
10 #include <AliESDEvent.h>
11 #include <AliHeader.h>
13 #include "AliPhysicsSelection.h"
14 #include "AliAnalysisManager.h"
15 #include "AliInputEventHandler.h"
17 //#include "AliBackgroundSelection.h"
19 ClassImp(AliPhysicsSelectionTask)
21 AliPhysicsSelectionTask::AliPhysicsSelectionTask() :
22 AliAnalysisTaskSE("AliPhysicsSelectionTask"),
28 // Default event handler
32 AliPhysicsSelectionTask::AliPhysicsSelectionTask(const char* opt) :
33 AliAnalysisTaskSE("AliPhysicsSelectionTask"),
36 fPhysicsSelection(new AliPhysicsSelection())
39 // Constructor. Initialization of pointers
41 AliInputEventHandler* handler = dynamic_cast<AliInputEventHandler*> (AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler());
43 handler->SetEventSelection(fPhysicsSelection);
44 AliInfo("Physics Event Selection enabled.");
46 AliError("No input event handler connected to analysis manager. No Physics Event Selection.");
48 // Define input and output slots here
49 DefineOutput(1, TList::Class());
51 AliLog::SetClassDebugLevel("AliPhysicsSelectionTask", AliLog::kWarning);
54 AliPhysicsSelectionTask::~AliPhysicsSelectionTask()
60 // histograms are in the output list and deleted when the output
61 // list is deleted by the TSelector dtor
69 void AliPhysicsSelectionTask::UserCreateOutputObjects()
71 // create result objects and add to output list
73 Printf("AliPhysicsSelectionTask::CreateOutputObjects");
78 if (!fPhysicsSelection)
79 fPhysicsSelection = new AliPhysicsSelection;
81 fOutput->Add(fPhysicsSelection);
84 void AliPhysicsSelectionTask::UserExec(Option_t*)
88 // AliPhysicsSelection::IsCollisionCandidate is called from the event handler
89 // post the data here anyway!
93 void AliPhysicsSelectionTask::Terminate(Option_t *)
95 // The Terminate() function is the last function to be called during
96 // a query. It always runs on the client, it can be used to present
97 // the results graphically or save the results to file.
99 fOutput = dynamic_cast<TList*> (GetOutputData(1));
101 Printf("ERROR: fOutput not available");
105 fPhysicsSelection = dynamic_cast<AliPhysicsSelection*> (fOutput->FindObject("AliPhysicsSelection"));
108 TFile* fout = new TFile("event_stat.root", "RECREATE");
110 if (fPhysicsSelection)
112 fPhysicsSelection->Print();
113 fPhysicsSelection->SaveHistograms();
119 Printf("Writting result to event_stat.root");