New Calorimeter Calibration directory, now only PHOS
[u/mrichter/AliRoot.git] / PWG4 / CaloCalib / macros / anaTask.C
1 void anaTask()
2 {
3   
4   AliCDBManager::Instance()->SetDefaultStorage("local://$ALICE_ROOT/OCDB"); 
5   AliCDBManager::Instance()->SetSpecificStorage("PHOS/*","local://./");
6   
7   // load analysis framework
8   gSystem->Load("libANALYSIS");
9   gSystem->Load("libANALYSISalice"); //AliAnalysisTaskSE
10   gSystem->Load("libPWG4CaloCalib");
11   
12   gROOT->LoadMacro("$ALICE_ROOT/PWG0/CreateESDChain.C");
13   TChain* chain = CreateESDChain("list.txt", 15);
14   
15   // for includes use either global setting in $HOME/.rootrc
16   // ACLiC.IncludePaths: -I$ALICE_ROOT/include
17   // or in each macro
18   // gSystem->AddIncludePath("$ALICE_ROOT/include");
19   
20   // Create the analysis manager
21   AliAnalysisManager *mgr = new AliAnalysisManager("PHOSPi0Calib");
22   
23   //Input event handler
24   AliVEventHandler* esdH = new AliESDInputHandler();
25   mgr->SetInputEventHandler(esdH);
26   
27   //Output event handler
28   AliAODHandler* aodoutHandler   = new AliAODHandler();
29   aodoutHandler->SetOutputFileName("aod.root");
30   mgr->SetOutputEventHandler(aodoutHandler);
31   
32   // ESD filter task
33   gROOT->LoadMacro("$ALICE_ROOT/PWG4/macros/AddTaskESDfilter.C");
34   AliAnalysisTaskESDfilter *esdfilter = AddTaskESDfilter(kFALSE);
35   
36   // Calibration task 
37   //
38   AliAnalysisTaskPHOSPi0CalibSelection *task = new AliAnalysisTaskPHOSPi0CalibSelection("PHOSPi0CalibSelection");
39   task->SetClusterMinEnergy(0.4); 
40   mgr->AddTask(task);
41   
42   // Create containers for input/output
43   AliAnalysisDataContainer *cinput   = mgr->GetCommonInputContainer();
44   AliAnalysisDataContainer *coutput2 = mgr->CreateContainer("histos", TList::Class(),  AliAnalysisManager::kOutputContainer, "PHOShistos.root");
45   
46   // Connect input/output
47   mgr->ConnectInput   (task, 0, cinput);
48   mgr->ConnectOutput(task, 1, coutput2);
49   
50   // Enable debug printouts
51   //mgr->SetDebugLevel(10);
52   
53   if (!mgr->InitAnalysis())
54     return;
55   
56   mgr->PrintStatus();
57   
58   mgr->StartAnalysis("local", chain);
59 }