]>
Commit | Line | Data |
---|---|---|
81ae6b8d | 1 | { |
2 | gSystem->Load("libPhysics.so"); | |
3 | // Connecting to the PROOF cluster | |
4 | TProof::Open("elopez@lxb6046"); | |
5 | ||
6 | // Clear packages if changing ROOT version on CAF or local | |
7 | //gProof->ClearPackages(); | |
8 | ||
9 | // Enable the needed packages | |
10 | gProof->UploadPackage("STEERBase"); | |
11 | gProof->EnablePackage("STEERBase"); | |
12 | gProof->UploadPackage("ESD"); | |
13 | gProof->EnablePackage("ESD"); | |
14 | gProof->UploadPackage("AOD"); | |
15 | gProof->EnablePackage("AOD"); | |
16 | gProof->UploadPackage("ANALYSIS"); | |
17 | gProof->EnablePackage("ANALYSIS"); | |
18 | gProof->UploadPackage("ANALYSISalice"); | |
19 | gProof->EnablePackage("ANALYSISalice"); | |
20 | gProof->UploadPackage("JETAN"); | |
21 | gProof->EnablePackage("JETAN"); | |
22 | ||
23 | // Create the analysis manager | |
24 | mgr = new AliAnalysisManager("Analysis UE test"); | |
25 | ||
26 | // Create, add task | |
27 | gProof->Load("AliKineTrackCuts.cxx+"); | |
28 | AliKineTrackCuts* trackCuts = new AliKineTrackCuts("AliKineTrackCuts", "Eta"); | |
29 | trackCuts->SetEtaRange(-1., 1.); | |
30 | ||
31 | AliAnalysisFilter* trackFilter = new AliAnalysisFilter("trackFilter"); | |
32 | trackFilter->AddCuts(trackCuts); | |
33 | ||
34 | gProof->Load("AliAnalysisTaskKineFilter.cxx+"); | |
35 | AliAnalysisTaskKineFilter *kinefilter = new AliAnalysisTaskKineFilter("Kine Filter"); | |
36 | kinefilter->SetTrackFilter(trackFilter); | |
37 | mgr->AddTask(kinefilter); | |
38 | ||
39 | // Create chain of input files | |
40 | gROOT->LoadMacro("CreateESDChain.C"); | |
41 | chain = CreateChain( "TE", "KINE82XX_30K.txt", 200); | |
42 | ||
43 | ||
44 | ||
45 | /////////////////////////////////////////////////////////////////////////////////// | |
46 | // Create the analysis manager | |
47 | // | |
48 | // Input | |
49 | // MC Truth | |
50 | AliMCEventHandler* mcHandler = new AliMCEventHandler(); | |
51 | mcHandler->SetReadTR(kFALSE); | |
52 | ||
53 | // Output | |
54 | AliAODHandler* aodHandler = new AliAODHandler(); | |
55 | aodHandler->SetOutputFileName( "aodKine.root" ); | |
56 | ||
57 | AliAnalysisManager *mgr = new AliAnalysisManager("Jet Manager", "Jet Manager"); | |
58 | // mgr->SetInputEventHandler (inpHandler); | |
59 | mgr->SetMCtruthEventHandler(mcHandler); | |
60 | mgr->SetOutputEventHandler (aodHandler); | |
61 | mgr->SetDebugLevel(10); | |
62 | ||
63 | /////////////////////////////////////////////////////////////////////////////////// | |
64 | ||
65 | // | |
66 | // Set of cuts | |
67 | // | |
68 | AliKineTrackCuts* trackCuts = new AliKineTrackCuts("AliKineTrackCuts", "Eta"); | |
69 | trackCuts->SetEtaRange(-1., 1.); | |
70 | // trackCuts->SetPtRange(0.5); | |
71 | ||
72 | AliAnalysisFilter* trackFilter = new AliAnalysisFilter("trackFilter"); | |
73 | trackFilter->AddCuts(trackCuts); | |
74 | ||
75 | ||
76 | AliAnalysisTaskKineFilter* kinefilter = new AliAnalysisTaskKineFilter("Kine Filter"); | |
77 | kinefilter->SetTrackFilter(trackFilter); | |
78 | kinefilter->SetDebugLevel(10); | |
79 | mgr->AddTask(kinefilter); | |
80 | ||
81 | // Analysis Task for Jet | |
82 | // AliAnalysisTaskJets need ConfigJetAnalysis.C macro !!!! | |
83 | AliAnalysisTaskJets *jetana = new AliAnalysisTaskJets("JetAnalysis"); | |
84 | jetana->SetDebugLevel(10); | |
85 | mgr->AddTask(jetana); | |
86 | ||
87 | // | |
88 | // Create containers for input/output | |
89 | AliAnalysisDataContainer *cinput1 = mgr->CreateContainer("cchain",TChain::Class(), | |
90 | AliAnalysisManager::kInputContainer); | |
91 | ||
92 | AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("tree", TTree::Class(), | |
93 | AliAnalysisManager::kOutputContainer, "default"); | |
94 | AliAnalysisDataContainer *coutput2 = mgr->CreateContainer("histos", TList::Class(), | |
95 | AliAnalysisManager::kOutputContainer, "histos.root"); | |
96 | ||
97 | mgr->ConnectInput (kinefilter, 0, cinput1 ); | |
98 | mgr->ConnectOutput (kinefilter, 0, coutput1 ); | |
99 | ||
100 | mgr->ConnectInput (jetana, 0, cinput1 ); | |
101 | mgr->ConnectOutput (jetana, 0, coutput1 ); | |
102 | mgr->ConnectOutput (jetana, 1, coutput2 ); | |
103 | ||
104 | ||
105 | // | |
106 | // Run the analysis | |
107 | // | |
108 | if( mgr->InitAnalysis() ) { | |
109 | mgr->PrintStatus(); | |
110 | mgr->StartAnalysis("proof", chain ); | |
111 | } | |
112 | } |