1 void runCentrality(const char * type = "a", const char *mode="grid")
3 // Load common libraries
4 gSystem->Load("libCore.so");
5 gSystem->Load("libTree.so");
6 gSystem->Load("libGeom.so");
7 gSystem->Load("libVMC.so");
8 gSystem->Load("libPhysics.so");
9 gSystem->Load("libSTEERBase");
10 gSystem->Load("libESD");
11 gSystem->Load("libAOD");
12 gSystem->Load("libANALYSIS");
13 gSystem->Load("libANALYSISalice");
14 // Use AliRoot includes to compile our task
15 gROOT->ProcessLine(".include $ALICE_ROOT/include");
18 TString filenameStr = Form("%s",type);
19 filenameStr = TString("LHC10g2")+filenameStr;
20 const char * filename = filenameStr.Data();
22 // Create and configure the alien handler plugin
23 gROOT->LoadMacro("CreateAlienHandlerAOD.C");
24 AliAnalysisGrid *alienHandler = CreateAlienHandlerAOD(filename);
25 if(!alienHandler) return;
27 // Create the analysis manager
28 AliAnalysisManager *mgr = new AliAnalysisManager("AnalysisManager");
30 // Connect plug-in to the analysis manager
31 mgr->SetGridHandler(alienHandler);
34 gROOT->LoadMacro("AliAnalysisTaskCentrality.cxx++g");
35 AliAnalysisTaskCentrality *task = new AliAnalysisTaskCentrality("CentralityTask");
36 // Writing (or not) output tree
37 //task->SetTreeFilling(writeTree);
41 AliMCEventHandler* mcHandler = new AliMCEventHandler();
42 mgr->SetMCtruthEventHandler(mcHandler);
44 AliESDInputHandler* esdH = new AliESDInputHandler();
45 mgr->SetInputEventHandler(esdH);
48 //gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskPhysicsSelection.C");
49 //AliPhysicsSelectionTask *physSelTask = AddTaskPhysicsSelection(kFALSE);
50 // Selecting collision candidates
51 //task->SelectCollisionCandidates();
53 // Create containers for input/output
54 mgr->ConnectInput(task, 0, mgr->GetCommonInputContainer());
56 AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("coutput1",TList::Class(),
57 AliAnalysisManager::kOutputContainer,
59 mgr->ConnectOutput(task, 1, coutput1);
61 // AliAnalysisDataContainer *coutput2 = mgr->CreateContainer("coutput2",TTree::Class(),
62 // AliAnalysisManager::kOutputContainer,
64 // // coutput2->SetSpecialOutput();
65 // mgr->ConnectOutput(task, 2, coutput2);
67 // Enable debug printouts
68 mgr->SetDebugLevel(2);
70 if(!mgr->InitAnalysis()) return;
73 // Start analysis in grid.
74 mgr->StartAnalysis(mode);