]>
Commit | Line | Data |
---|---|---|
b47f90d3 | 1 | /* |
78abca19 | 2 | Example usage: |
3 | ||
4 | 0. Load neccessary libraries | |
5 | ||
b47f90d3 | 6 | gSystem->Load("libANALYSIS"); |
7 | gSystem->Load("libANALYSISalice"); | |
8 | gSystem->Load("libTPCcalib"); | |
78abca19 | 9 | gSystem->Load("libXrdClient.so"); |
10 | gSystem->Load("libNetx.so"); | |
11 | TGrid::Connect("alien://",0,0,"t"); | |
12 | gSystem->Load("$ROOTSYS/lib/libXrdClient.so"); | |
13 | // | |
14 | // 1. Make list of the files | |
15 | // | |
16 | .L $ALICE_ROOT/TPC/macros/testTPC/AlienToolkit.cxx+ | |
17 | gSystem->Load("libXrdClient.so"); | |
18 | gSystem->Load("libNetx.so"); | |
19 | AlienToolkit toolkit; | |
20 | char *path = "/alice/cern.ch/user/a/amatyja/alice/data/" | |
21 | toolkit.MakeCollection(path,"32129*Krypton.root"); | |
22 | toolkit.MakeCollection(path,"32231*Krypton.root"); | |
23 | toolkit.MakeCollection(path,"32284*Krypton.root"); | |
24 | toolkit.PrintPFN(); > list.txt | |
25 | ||
26 | // | |
27 | // 2. Initialization of proof | |
28 | // | |
29 | TProofMgr * proofmgr = TProof::Mgr("lxgrid5.gsi.de"); | |
30 | TProof * proof = proofmgr->CreateSession(); | |
31 | proof->SetParameter("PROOF_MaxSlavesPerNode", (Long_t)1000); | |
32 | .L /u/miranov/macros/ProofEnableAliRoot.C | |
33 | ProofEnableAliRoot("/usr/local/grid/AliRoot/HEAD0108"); | |
34 | gProof->Exec("gSystem->Load(\"libANALYSIS.so\")",kTRUE); | |
35 | gProof->Exec("gSystem->Load(\"libSTAT.so\")",kTRUE); | |
36 | gProof->Exec("gSystem->Load(\"libTPCcalib.so\")",kTRUE); | |
37 | ||
b47f90d3 | 38 | // |
78abca19 | 39 | // 3 Run analysis on PROOF |
b47f90d3 | 40 | // |
78abca19 | 41 | // |
42 | // Create chain of input files | |
43 | // | |
44 | gSystem->AddIncludePath("-I$ALICE_ROOT/TPC/macros"); | |
45 | gROOT->LoadMacro("$ALICE_ROOT/TPC/macros/AliXRDPROOFtoolkit.cxx+"); | |
46 | .L $ALICE_ROOT/TPC/macros/RunAliTPCCalibKrTask.C | |
47 | RunAliTPCCalibKrTask(kTRUE); | |
48 | ||
b47f90d3 | 49 | |
78abca19 | 50 | //.. |
51 | // | |
52 | // | |
53 | ||
b47f90d3 | 54 | //Example usage |
55 | TFile f("KrHisto.root"); | |
56 | AliTPCCalibKr *kr = f.Get("AliTPCCalibKr"); | |
57 | ||
58 | kr->ProjectHisto(kr->GetHistoKr(71),"aaa",30,36,30,40)->Draw() | |
59 | // | |
60 | // | |
61 | // | |
62 | MakeTree(); | |
63 | //default cuts | |
64 | TCut cutKr("cutKr","entries.fElements<5000&&fitNormChi2.fElements<3&&fitNormChi2.fElements>0.2&&abs(fitRMS.fElements/fitMean.fElements-0.06)<0.025"); | |
65 | ||
cb6e6903 | 66 | TObjArray * array = AliTPCCalibViewerGUI::ShowGUI("kryptonTree.root"); |
67 | AliTPCCalibViewerGUI * viewer = (AliTPCCalibViewerGUI*)array->At(0); | |
68 | TTree * tree = viewer->GetViewer()->GetTree(); | |
69 | ||
70 | tree->SetAlias("cutAll","abs(fitNormChi2.fElements-2.)<1.8&&entries.fElements/entries_Median.fElements<4&&entries.fElements/entries_Median.fElements>0.4&&fitRMS.fElements/fitMean.fElements<0.09&&fitRMS.fElements/fitMean.fElements>0.02") | |
71 | ||
78abca19 | 72 | */ |
b47f90d3 | 73 | |
78abca19 | 74 | TChain * chain = 0; |
b47f90d3 | 75 | |
78abca19 | 76 | void RunAliTPCCalibKrTask(Bool_t bProof = kFALSE) |
ffbf9f24 | 77 | { |
ffbf9f24 | 78 | // |
78abca19 | 79 | AliXRDPROOFtoolkit tool; |
80 | chain = tool.MakeChain("list.txt","Kr","",20000,0); | |
81 | chain->Lookup(); | |
82 | chain->SetBranchStatus("Cl.fCluster",kFALSE); | |
83 | ||
ffbf9f24 | 84 | // |
85 | // Create the analysis manager | |
86 | // | |
87 | AliAnalysisManager *mgr = new AliAnalysisManager("testAnalysis"); | |
88 | ||
89 | // Calibration component | |
90 | AliTPCCalibKr *calibObj = new AliTPCCalibKr; | |
91 | //calibObj->SetASide(kFALSE); | |
92 | ||
93 | // Add task | |
94 | AliTPCCalibKrTask *task = new AliTPCCalibKrTask; | |
95 | task->SetInputChain(chain); | |
96 | task->SetTPCCalibKr(calibObj); | |
97 | mgr->AddTask(task); | |
98 | ||
99 | // Attach input | |
100 | cInput = mgr->CreateContainer("cInput", TChain::Class(), AliAnalysisManager::kInputContainer); | |
101 | mgr->ConnectInput(task, 0, cInput); | |
102 | ||
103 | // Attach output | |
104 | cOutput = mgr->CreateContainer("cOutput", TList::Class(), AliAnalysisManager::kOutputContainer,"outHistFile.root"); | |
105 | mgr->ConnectOutput(task, 0, cOutput); | |
78abca19 | 106 | // |
107 | cOutput->SetSpecialOutput(kTRUE); | |
108 | cOutput->SetFileName("CalibObjectFile.root"); | |
109 | // | |
ffbf9f24 | 110 | // Run analysis |
78abca19 | 111 | if (!mgr->InitAnalysis()) return; |
ffbf9f24 | 112 | mgr->PrintStatus(); |
78abca19 | 113 | mgr->SetDebugLevel(1); |
114 | ||
115 | if(bProof) { | |
116 | mgr->StartAnalysis("proof", chain); | |
117 | } | |
ffbf9f24 | 118 | else mgr->StartAnalysis("local", chain); |
119 | } | |
b47f90d3 | 120 | |
121 | ||
122 | ||
123 | void MakeTree(){ | |
124 | ||
125 | TFile fpad("calibKr.root"); | |
126 | AliTPCPreprocessorOnline * preprocesor = new AliTPCPreprocessorOnline; | |
127 | preprocesor->AddComponent(spectrMean->Clone()); | |
128 | preprocesor->AddComponent(spectrRMS->Clone()); | |
129 | preprocesor->AddComponent(fitMean->Clone()); | |
130 | preprocesor->AddComponent(fitRMS->Clone()); | |
131 | preprocesor->AddComponent(fitNormChi2->Clone()); | |
132 | preprocesor->AddComponent(entries->Clone()); | |
133 | preprocesor->DumpToFile("kryptonTree.root"); | |
134 | ||
135 | } |