1 void qa(Int_t runNumber) {
5 gEnv->SetValue("XSec.GSI.DelegProxy","2");
7 TProof::Mgr("aliprod@alicecaf")->SetROOTVersion("v5-24-00b-caf");
9 TProof::Open("aliprod@alicecaf");
12 gProof->UploadPackage("/afs/cern.ch/alice/caf/sw/ALICE/PARs/v4-17-Release.rec/AF-v4-17-rec.par");
13 gProof->EnablePackage("AF-v4-17-rec.par");
15 // Enable analysis libs
16 gSystem->Load("libANALYSIS.so");
17 gSystem->Load("libANALYSISalice.so");
18 gProof->Exec("gSystem->Load(\"libANALYSIS.so\");",kTRUE);
19 gProof->Exec("gSystem->Load(\"libANALYSISalice.so\");",kTRUE);
21 gProof->Load(Form("%s/PWG1/cosmic/AliAnalysisTaskCosmic.cxx++g",
22 gSystem->Getenv("ALICE_ROOT")));
24 //____________________________________________//
25 // Make the analysis manager
26 AliAnalysisManager *mgr = new AliAnalysisManager("AnalysisQAManager");
27 AliESDInputHandler* esdH = new AliESDInputHandler();
28 mgr->SetInputEventHandler(esdH);
29 mgr->SetDebugLevel(10);
31 //____________________________________________//
33 AliAnalysisTaskCosmic *task1 = new AliAnalysisTaskCosmic("TaskCosmic");
36 // Create containers for input/output
37 AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("chist1",TList::Class(),AliAnalysisManager::kOutputContainer,
38 Form("run%d.root",runNumber));
40 //____________________________________________//
41 mgr->ConnectInput (task1, 0, mgr->GetCommonInputContainer());
42 mgr->ConnectOutput (task1, 1, coutput1);
44 if (!mgr->InitAnalysis()) return;
46 mgr->StartAnalysis("proof",
47 Form("/ALIREC/aliprod/run%d",runNumber));
55 void plot(Int_t runNumber)
57 TFile* f = new TFile(Form("run%d.root",runNumber), "read");
60 TCanvas* c1 = new TCanvas("c1", "pt,eta,phi", 10, 10, 1100, 800);
63 c1->GetPad(1)->SetLogy();
65 fhPtP = (TH1F*) chist1->FindObject("fhPtPC");
66 fhPtP->SetLineColor(2);
68 fhPtN = (TH1F*) chist1->FindObject("fhPtNC");
69 fhPtN->SetLineColor(4);
73 fhPhiP = (TH1F*) chist1->FindObject("fhPhiPC");
74 fhPhiP->SetLineColor(2);
76 fhPhiN = (TH1F*) chist1->FindObject("fhPhiNC");
77 fhPhiN->SetLineColor(4);
81 fhThetaP = (TH1F*) chist1->FindObject("fhThetaPC");
82 fhThetaP->SetLineColor(2);
84 fhThetaN = (TH1F*) chist1->FindObject("fhThetaNC");
85 fhThetaN->SetLineColor(4);
86 fhThetaN->Draw("same");
90 TCanvas* c2 = new TCanvas("c2", "#Delta#phi", 10, 10, 1100, 800);
94 c2->GetPad(1)->SetLogy();
95 fhDPhiP = (TH1F*) chist1->FindObject("fhDPhiPC");
96 fhDPhiN = (TH1F*) chist1->FindObject("fhDPhiNC");
97 fhDPhiP->SetLineColor(2);
98 fhDPhiN->SetLineColor(4);
99 fhDPhiP->SetXTitle("#Delta#phi [rad]");
101 fhDPhiN->Draw("same");
105 c2->GetPad(2)->SetLogy();
107 fhDPhiZP = (TH1F*) chist1->FindObject("fhDPhiPZ");
108 fhDPhiZN = (TH1F*) chist1->FindObject("fhDPhiNZ");
109 fhDPhiBA = (TH1F*) chist1->FindObject("fhDPhi_Bad");
110 fhDPhiZP->SetLineColor(2);
111 fhDPhiZN->SetLineColor(4);
112 fhDPhiBA->SetLineColor(6);
114 fhDPhiZP->SetXTitle("#Delta#phi [rad]");
116 fhDPhiZN->Draw("same");
117 // fhDPhiBA->Draw("same");
123 c2->GetPad(3)->SetLogy();
124 fhDThetaP = (TH1F*) chist1->FindObject("fhDThetaPC");
125 fhDThetaN = (TH1F*) chist1->FindObject("fhDThetaNC");
126 fhDThetaP->SetLineColor(2);
127 fhDThetaN->SetLineColor(4);
130 fhDThetaP->SetXTitle("#Delta#theta [rad]");
132 fhDThetaN->Draw("same");
136 c2->GetPad(4)->SetLogy();
138 fhDThetaZP = (TH1F*) chist1->FindObject("fhDThetaPZ");
139 fhDThetaZN = (TH1F*) chist1->FindObject("fhDThetaNZ");
140 fhDThetaBA = (TH1F*) chist1->FindObject("fhDTheta_Bad");
141 fhDThetaZP->SetLineColor(2);
142 fhDThetaZN->SetLineColor(4);
143 fhDThetaBA->SetLineColor(6);
144 fhDThetaZP->SetXTitle("#Delta#theta [rad]");
147 fhDThetaZN->Draw("same");
148 // fhDThetaBA->Draw("same");
152 TCanvas* c3 = new TCanvas("c3", "#Delta p_{T}", 10, 10, 1100, 800);
156 c3->GetPad(1)->SetLogy();
157 fhDPtP = (TH1F*) chist1->FindObject("fhDPtPC");
158 fhDPtN = (TH1F*) chist1->FindObject("fhDPtNC");
159 fhDPtP->SetLineColor(2);
160 fhDPtN->SetLineColor(4);
163 fhDPtN->Draw("same");
167 c3->GetPad(2)->SetLogy();
169 fhDPtZP = (TH1F*) chist1->FindObject("fhDPtPZ");
170 fhDPtZN = (TH1F*) chist1->FindObject("fhDPtNZ");
171 fhDPtZP->SetLineColor(2);
172 fhDPtZN->SetLineColor(4);
174 fhDPtZN->Draw("same");
178 fpDPtP = (TH1F*) chist1->FindObject("fpDPtPC");
179 fpDPtN = (TH1F*) chist1->FindObject("fpDPtNC");
180 fpDPtP->SetLineColor(2);
181 fpDPtN->SetLineColor(4);
183 fpDPtN->Draw("same");
187 fpDPtZP = (TH1F*) chist1->FindObject("fpDPtPZ");
188 fpDPtZN = (TH1F*) chist1->FindObject("fpDPtNZ");
189 fpDPtZP->SetLineColor(2);
190 fpDPtZN->SetLineColor(4);
192 fpDPtZN->Draw("same");
196 TCanvas* c4 = new TCanvas("c4", "#Delta Z", 10, 10, 1100, 800);
200 c4->GetPad(1)->SetLogy();
201 fhDZP = (TH1F*) chist1->FindObject("fhDZPC");
202 fhDZN = (TH1F*) chist1->FindObject("fhDZNC");
203 fhDZP->SetLineColor(2);
204 fhDZN->SetLineColor(4);
211 c4->GetPad(2)->SetLogy();
213 fhDZZP = (TH1F*) chist1->FindObject("fhDZPZ");
214 fhDZZN = (TH1F*) chist1->FindObject("fhDZNZ");
215 fhDZBA = (TH1F*) chist1->FindObject("fhDZ_Bad");
217 fhDZZP->SetLineColor(2);
218 fhDZZN->SetLineColor(4);
219 fhDZBA->SetLineColor(6);
222 fhDZZN->Draw("same");
223 fhDZBA->Draw("same");
228 c4->GetPad(3)->SetLogy();
229 fhDXP = (TH1F*) chist1->FindObject("fhDXPC");
230 fhDXN = (TH1F*) chist1->FindObject("fhDXNC");
231 fhDXP->SetLineColor(2);
232 fhDXN->SetLineColor(4);
239 c4->GetPad(4)->SetLogy();
241 fhDXZP = (TH1F*) chist1->FindObject("fhDXPZ");
242 fhDXZN = (TH1F*) chist1->FindObject("fhDXNZ");
243 fhDXBA = (TH1F*) chist1->FindObject("fhDX_Bad");
245 fhDXZP->SetLineColor(2);
246 fhDXZN->SetLineColor(4);
247 fhDXBA->SetLineColor(6);
250 fhDXZN->Draw("same");
251 fhDXBA->Draw("same");
254 TCanvas* c4a = new TCanvas("c4a", "#Delta Y", 10, 10, 1100, 800);
258 c4a->GetPad(1)->SetLogy();
259 fhDYP = (TH1F*) chist1->FindObject("fhDYPC");
260 fhDYN = (TH1F*) chist1->FindObject("fhDYNC");
261 fhDYP->SetLineColor(2);
262 fhDYN->SetLineColor(4);
269 c4a->GetPad(2)->SetLogy();
271 fhDYZP = (TH1F*) chist1->FindObject("fhDYPZ");
272 fhDYZN = (TH1F*) chist1->FindObject("fhDYNZ");
273 fhDYBA = (TH1F*) chist1->FindObject("fhDY_Bad");
275 fhDYZP->SetLineColor(2);
276 fhDYZN->SetLineColor(4);
277 fhDYBA->SetLineColor(6);
280 fhDYZN->Draw("same");
281 fhDYBA->Draw("same");
284 TCanvas* c5 = new TCanvas("c5", "#Delta p_{T} (n-sigma)", 10, 10, 1100, 800);
288 fpDPtSP = (TH1F*) chist1->FindObject("fpDPtSPC");
289 fpDPtSN = (TH1F*) chist1->FindObject("fpDPtSNC");
290 fpDPtSP->SetLineColor(2);
291 fpDPtSN->SetLineColor(4);
293 fpDPtSN->Draw("same");
298 fpDPtSPZ = (TH1F*) chist1->FindObject("fpDPtSPZ");
299 fpDPtSNZ = (TH1F*) chist1->FindObject("fpDPtSNZ");
300 fpDPtSPZ->SetLineColor(2);
301 fpDPtSNZ->SetLineColor(4);
303 fpDPtSNZ->Draw("same");
305 TCanvas* c6 = new TCanvas("c6", "Dz vs z", 10, 10, 1100, 800);
306 fhDZvsZ = (TH2F*) chist1->FindObject("fhDZvsZ");
307 fhDZvsZ->SetXTitle("z_{in} * sign(z_{in}) * sign(z_{out}) [cm]");
308 fhDZvsZ->SetYTitle("#DeltaZ [cm]");
310 gStyle->SetPalette(1);
312 fhDZvsZ->Draw("colz");