]>
Commit | Line | Data |
---|---|---|
c88234ad | 1 | void runAOD(Bool_t isMC = 0) |
2 | { | |
3 | ||
4 | //gEnv->SetValue("XSec.GSI.DelegProxy", "2"); | |
5 | //TProof::Open("pverstee@alice-caf.cern.ch"); | |
6 | ||
7 | gSystem->Load("libTree.so"); | |
8 | gSystem->Load("libGeom.so"); | |
9 | gSystem->Load("libVMC.so"); | |
10 | gSystem->Load("libPhysics.so"); | |
11 | gSystem->Load("libSTEERBase.so"); | |
12 | gSystem->Load("libESD.so"); | |
13 | gSystem->Load("libAOD.so"); | |
14 | gSystem->Load("libANALYSIS.so"); | |
15 | gSystem->Load("libANALYSISalice.so"); | |
16 | gSystem->Load("libANALYSIS"); | |
17 | gSystem->Load("libANALYSISalice"); | |
18 | gROOT->LoadMacro("$ALICE_ROOT/PWG0/CreateESDChain.C"); | |
19 | ||
20 | ifstream flist; | |
21 | flist.open("filelist.txt"); | |
22 | TChain * chain = new TChain("aodTree"); | |
23 | TString line; | |
24 | while (line.ReadLine(flist)) | |
25 | { | |
26 | gSystem->ExpandPathName(line); | |
27 | chain->AddFile(line.Data()); | |
28 | cout << "Adding file " << line.Data() << endl; | |
29 | } | |
30 | // gProof->UploadPackage("AF-v4-19-04-AN"); | |
31 | // gProof->EnablePackage("AF-v4-19-04-AN"); | |
32 | gSystem->AddIncludePath("-I$ALICE_ROOT/include"); | |
33 | gStyle->SetPalette(1); | |
34 | ||
35 | ||
36 | ||
37 | AliAnalysisManager *mgr = new AliAnalysisManager("testAnalysis"); | |
38 | AliAODInputHandler* aodH = new AliAODInputHandler(); | |
39 | mgr->SetInputEventHandler(aodH); | |
40 | ||
41 | gROOT->LoadMacro("AliSpectraAODTrackCuts.cxx+g"); | |
42 | gROOT->LoadMacro("AliSpectraAODVertexCuts.cxx+g"); | |
43 | gROOT->LoadMacro("AliSpectraAODHistoManager.cxx+g"); | |
44 | gROOT->LoadMacro("AliAnalysisTaskSpectraAOD.cxx+g"); | |
45 | ||
46 | gROOT->LoadMacro("$(ALICE_ROOT)/ANALYSIS/macros/AddTaskPIDResponse.C"); | |
47 | AliAnalysisTask * taskPID = AddTaskPIDResponse(); | |
48 | mgr->AddTask(taskPID); | |
49 | ||
50 | ||
51 | AliAnalysisTaskSpectraAOD *task = new AliAnalysisTaskSpectraAOD("TaskAODExercise"); | |
52 | task->SetIsMC(isMC); | |
53 | mgr->AddTask(task); | |
54 | ||
55 | // Set the cuts | |
56 | AliSpectraAODVertexCuts * vcuts = new AliSpectraAODVertexCuts("Vertex Cuts"); | |
57 | AliSpectraAODTrackCuts * tcuts = new AliSpectraAODTrackCuts("Tracks Cuts"); | |
58 | tcuts->SetTrackType(6); | |
59 | tcuts->SetEta(1.); | |
60 | task->SetVertexCuts(vcuts); | |
61 | task->SetTrackCuts(tcuts); | |
62 | ||
63 | ||
64 | AliAnalysisDataContainer *cinput = mgr->GetCommonInputContainer(); | |
65 | AliAnalysisDataContainer *coutputpt1 = mgr->CreateContainer("chistpt", AliSpectraAODHistoManager::Class(), AliAnalysisManager::kOutputContainer, "Pt.AOD.1.root"); | |
66 | AliAnalysisDataContainer *coutputpt2 = mgr->CreateContainer("cvcutpt", AliSpectraAODVertexCuts::Class(), AliAnalysisManager::kOutputContainer, "Pt.AOD.1.root"); | |
67 | AliAnalysisDataContainer *coutputpt3 = mgr->CreateContainer("ctcutpt", AliSpectraAODTrackCuts::Class(), AliAnalysisManager::kOutputContainer, "Pt.AOD.1.root"); | |
68 | ||
69 | mgr->ConnectInput(task, 0, cinput); | |
70 | mgr->ConnectOutput(task, 1, coutputpt1); | |
71 | mgr->ConnectOutput(task, 2, coutputpt2); | |
72 | mgr->ConnectOutput(task, 3, coutputpt3); | |
73 | mgr->SetDebugLevel(2); | |
74 | ||
75 | if (!mgr->InitAnalysis()) return; | |
76 | mgr->PrintStatus(); | |
77 | mgr->StartAnalysis("local", chain, 100); | |
78 | } |