Line data Source code
1 : // $Id$
2 :
3 : //**************************************************************************
4 : //* This file is property of and copyright by the ALICE HLT Project *
5 : //* ALICE Experiment at CERN, All rights reserved. *
6 : //* *
7 : //* Primary Authors: Zhong-Bao Yin <Zhong-Bao.Yin@cern.ch> *
8 : //* Matthias Richter <Matthias.Richter@ift.uib.no> *
9 : //* for The ALICE HLT Project. *
10 : //* *
11 : //* Permission to use, copy, modify and distribute this software and its *
12 : //* documentation strictly for non-commercial purposes is hereby granted *
13 : //* without fee, provided that the above copyright notice appears in all *
14 : //* copies and that both the copyright notice and this permission notice *
15 : //* appear in the supporting documentation. The authors make no claims *
16 : //* about the suitability of this software for any purpose. It is *
17 : //* provided "as is" without express or implied warranty. *
18 : //**************************************************************************
19 :
20 : /** @file AliHLTTPCQADataMaker.cxx
21 : @author Zhongbao Yin, Matthias Richter
22 : @date 2009-05-14
23 : @brief Container for the HLT offline QA
24 : */
25 : #include "AliHLTTPCQADataMaker.h"
26 : #include "AliESDEvent.h"
27 : #include <iostream>
28 :
29 : #include "TH1F.h"
30 : #include "TH2F.h"
31 :
32 : #include "AliESDtrack.h"
33 :
34 : using namespace std;
35 :
36 : /** ROOT macro for the implementation of ROOT specific class methods */
37 6 : ClassImp(AliHLTTPCQADataMaker)
38 :
39 : AliHLTTPCQADataMaker::AliHLTTPCQADataMaker()
40 2 : : AliHLTQADataMakerBase()
41 10 : {
42 : // see header file for class documentation
43 : // or
44 : // refer to README to build package
45 : // or
46 : // visit http://web.ift.uib.no/~kjeks/doc/alice-hlt
47 4 : }
48 :
49 : AliHLTTPCQADataMaker::~AliHLTTPCQADataMaker()
50 0 : {
51 : // see header file for class documentation
52 0 : }
53 :
54 : void AliHLTTPCQADataMaker::StartOfDetectorCycle()
55 : {
56 : // see header file for class documentation
57 0 : }
58 :
59 : void AliHLTTPCQADataMaker::EndOfDetectorCycle(AliQAv1::TASKINDEX_t, TObjArray** /*list*/)
60 : {
61 : // see header file for class documentation
62 0 : for (int itc=-1;itc<GetNTrigClasses();itc++) {
63 0 : TH1* hp = GetESDsData(kPHLT,itc);
64 0 : TH1* hr = GetESDsData(kPRatio,itc);
65 0 : TH1* ho = GetESDsData(kPOffline,itc);
66 : //
67 0 : if (hp && hr && ho && hp->GetEntries()) {
68 0 : hp->Sumw2();
69 0 : hr->Add(ho);
70 0 : hr->Divide(hp);
71 0 : }
72 : //
73 0 : hp = GetESDsData(kPHLTFired,itc);
74 0 : hr = GetESDsData(kPRatioFired,itc);
75 0 : ho = GetESDsData(kPOfflineFired,itc);
76 : //
77 0 : if (hp && hr && ho && hp->GetEntries()) {
78 0 : hp->Sumw2();
79 0 : hr->Add(ho);
80 0 : hr->Divide(hp);
81 0 : }
82 : }
83 0 : }
84 :
85 : void AliHLTTPCQADataMaker::MakeRaws(AliRawReader * rawReader)
86 : {
87 : // see header file for class documentation
88 : if (!rawReader) return;
89 0 : }
90 :
91 : void AliHLTTPCQADataMaker::InitESDs(){
92 :
93 : //create ESDs histograms in ESDs subdir
94 : const Bool_t expert = kTRUE ;
95 : const Bool_t image = kTRUE ;
96 :
97 : TH2F * histESDMultiplicity =
98 0 : new TH2F("hESDMultiplicity",
99 : "Number of tracks from all events; Number of offline tracks; Number of HLT tracks",
100 : 300, 0, 300, 300, 0, 300);
101 0 : histESDMultiplicity->Sumw2();
102 0 : Add2ESDsList(histESDMultiplicity, kMultiplicity, !expert, image);
103 :
104 : TH2F * histESDMultiplicityFired =
105 0 : new TH2F("hESDMultiplicityFired",
106 : "Number of tracks from HLT triggered events; Number of offline tracks; Number of HLT tracks",
107 : 300, 0, 300, 300, 0, 300);
108 0 : histESDMultiplicityFired->Sumw2();
109 0 : Add2ESDsList(histESDMultiplicityFired, kMultiplicityFired, !expert, image);
110 :
111 : TH2F * histESDNCls =
112 0 : new TH2F("hESDNCls", "Mean number of TPC clusters from all events; Offline; HLT",
113 : 200, 0, 200, 200, 0, 200);
114 0 : histESDNCls->Sumw2();
115 0 : Add2ESDsList(histESDNCls, kNCls, !expert, image);
116 :
117 : TH2F * histESDNClsFired =
118 0 : new TH2F("hESDNClsFired", "Mean number of TPC clusters from triggered events; Offline; HLT",
119 : 200, 0, 200, 200, 0, 200);
120 0 : histESDNClsFired->Sumw2();
121 0 : Add2ESDsList(histESDNClsFired, kNClsFired, !expert, image);
122 :
123 : TH1F * histPHLT
124 0 : = new TH1F("hPHLT", "P distribution for all events from HLT; P [GeV/c]",
125 : 100, -0.5, 99.5);
126 0 : histPHLT->Sumw2();
127 0 : Add2ESDsList(histPHLT, kPHLT, !expert, image);
128 :
129 : TH1F * histPOffline
130 0 : = new TH1F("hPOffline",
131 : "P distribution for all events from offline; P [GeV/c]",
132 : 100, -0.5, 99.5);
133 0 : histPOffline->Sumw2();
134 0 : Add2ESDsList(histPOffline, kPOffline, !expert, image);
135 :
136 : TH1F * histPHLTFired
137 0 : = new TH1F("hPHLTFired",
138 : "P distribution for fired events from HLT; P [GeV/c]",
139 : 100, -0.5, 99.5);
140 0 : histPHLTFired->Sumw2();
141 0 : Add2ESDsList(histPHLTFired, kPHLTFired, !expert, image);
142 :
143 : TH1F * histPOfflineFired
144 0 : = new TH1F("hPOfflineFired",
145 : "P distribution for fired events from offline; P [GeV/c]",
146 : 100, -0.5, 99.5);
147 0 : histPOfflineFired->Sumw2();
148 0 : Add2ESDsList(histPOfflineFired, kPOfflineFired, !expert, image);
149 :
150 : TH1F * histPRatio
151 0 : = new TH1F("hPRatio",
152 : "Ratio of P distribution for all events; P [GeV/c]",
153 : 100, -0.5, 99.5);
154 0 : histPRatio->Sumw2();
155 0 : Add2ESDsList(histPRatio, kPRatio, !expert, image);
156 :
157 : TH1F * histPRatioFired
158 0 : = new TH1F("hPRatioFired",
159 : "Ratio of P distribution for fired events; P [GeV/c]",
160 : 100, -0.5, 99.5);
161 0 : histPRatioFired->Sumw2();
162 0 : Add2ESDsList(histPRatioFired, kPRatioFired, !expert, image);
163 :
164 : TH1F * histPtHLT
165 0 : = new TH1F("hPtHLT",
166 : "P_{T} distribution for all events from HLT; P_{T} [GeV/c]",
167 : 200, -100.5, 99.5);
168 0 : histPtHLT->Sumw2();
169 0 : Add2ESDsList(histPtHLT, kPtHLT, !expert, image);
170 :
171 : TH1F * histPtOffline
172 0 : = new TH1F("hPtOffline",
173 : "P_{T} distribution for all events from offline; P_{T} [GeV/c]",
174 : 200, -100.5, 99.5);
175 0 : histPtOffline->Sumw2();
176 0 : Add2ESDsList(histPtOffline, kPtOffline, !expert, image);
177 :
178 : TH1F * histPtHLTFired
179 0 : = new TH1F("hPtHLTFired",
180 : "P_{T} distribution for fired events from HLT; P_{T} [GeV/c]",
181 : 200, -100.5, 99.5);
182 0 : histPtHLTFired->Sumw2();
183 0 : Add2ESDsList(histPtHLTFired, kPtHLTFired, !expert, image);
184 :
185 : TH1F * histPtOfflineFired
186 0 : = new TH1F("hPtOfflineFired",
187 : "P_{T} distribution for fired events from offline; P_{T} [GeV/c]",
188 : 200, -100.5, 99.5);
189 0 : histPtOfflineFired->Sumw2();
190 0 : Add2ESDsList(histPtOfflineFired, kPtOfflineFired, !expert, image);
191 :
192 : TH1F * histNClsPerTrkHLT
193 0 : = new TH1F("hNClsPerTrkHLT", "Clusters per HLT track; N cluster; Counts",
194 : 200, 0, 200);
195 0 : histNClsPerTrkHLT->Sumw2();
196 0 : Add2ESDsList(histNClsPerTrkHLT, kNClsPerTrkHLT, !expert, image);
197 :
198 : TH1F * histNClsPerTrkOffline
199 0 : = new TH1F("hNClsPerTrkOffline",
200 : "Clusters per offline track; N cluster; Counts",
201 : 200, 0, 200);
202 0 : histNClsPerTrkOffline->Sumw2();
203 0 : Add2ESDsList(histNClsPerTrkOffline, kNClsPerTrkOffline, !expert, image);
204 :
205 : TH1F * histNClsPerTrkHLTFired
206 0 : = new TH1F("hNClsPerTrkHLTFired",
207 : "Clusters per HLT track from HLT fired events; N cluster; Counts",
208 : 200, 0, 200);
209 0 : histNClsPerTrkHLTFired->Sumw2();
210 0 : Add2ESDsList(histNClsPerTrkHLTFired, kNClsPerTrkHLTFired, !expert, image);
211 :
212 : TH1F * histNClsPerTrkOfflineFired
213 0 : = new TH1F("hNClsPerTrkOfflineFired",
214 : "Clusters per offline track from HLT fired events; N cluster; Counts",
215 : 200, 0, 200);
216 0 : histNClsPerTrkOfflineFired->Sumw2();
217 0 : Add2ESDsList(histNClsPerTrkOfflineFired, kNClsPerTrkOfflineFired, !expert, image);
218 :
219 : TH1F * histPhiHLT =
220 0 : new TH1F("hPhiHLT", "Phi distribution of HLT tracks; Phi; Counts",
221 : 360, 0, 360);
222 0 : histPhiHLT->Sumw2();
223 0 : Add2ESDsList(histPhiHLT, kPhiHLT, !expert, image);
224 :
225 : TH1F * histPhiOffline =
226 0 : new TH1F("hPhiOffline",
227 : "Phi distribution of offline tracks; Phi; Counts",
228 : 360, 0, 360);
229 0 : histPhiOffline->Sumw2();
230 0 : Add2ESDsList(histPhiOffline, kPhiOffline, !expert, image);
231 :
232 : TH1F * histPhiHLTFired =
233 0 : new TH1F("hPhiHLTFired", "Phi distribution of HLT tracks from HLT fired event ; Phi; Counts",
234 : 360, 0, 360);
235 0 : histPhiHLTFired->Sumw2();
236 0 : Add2ESDsList(histPhiHLTFired, kPhiHLTFired, !expert, image);
237 :
238 : TH1F * histPhiOfflineFired =
239 0 : new TH1F("hPhiOfflineFired",
240 : "Phi distribution of offline tracks from HLT fired events; Phi; Counts",
241 : 360, 0, 360);
242 0 : histPhiOfflineFired->Sumw2();
243 0 : Add2ESDsList(histPhiOfflineFired, kPhiOfflineFired, !expert, image);
244 :
245 : TH1F * histEtaHLT =
246 0 : new TH1F("hEtaHLT", "Eta distribution of HLT tracks; Eta; Counts",
247 : 200, -1, 1);
248 0 : histEtaHLT->Sumw2();
249 0 : Add2ESDsList(histEtaHLT, kEtaHLT, !expert, image);
250 :
251 : TH1F * histEtaOffline =
252 0 : new TH1F("hEtaOffline", "Eta distribution of offline tracks; Eta; Counts",
253 : 200, -1, 1);
254 0 : histEtaHLT->Sumw2();
255 0 : Add2ESDsList(histEtaOffline, kEtaOffline, !expert, image);
256 :
257 : TH1F * histEtaHLTFired =
258 0 : new TH1F("hEtaHLTFired",
259 : "Eta distribution of HLT tracks from HLT fired events; Eta; Counts",
260 : 200, -1, 1);
261 0 : histEtaHLTFired->Sumw2();
262 0 : Add2ESDsList(histEtaHLTFired, kEtaHLTFired, !expert, image);
263 :
264 : TH1F * histEtaOfflineFired =
265 0 : new TH1F("hEtaOfflineFired",
266 : "Eta distribution of offline tracks from HLT fired events; Eta; Counts",
267 : 200, -1, 1);
268 0 : histEtaHLTFired->Sumw2();
269 0 : Add2ESDsList(histEtaOfflineFired, kEtaOfflineFired, !expert, image);
270 :
271 0 : }
272 :
273 : void AliHLTTPCQADataMaker::MakeESDs(AliESDEvent * esd, AliESDEvent* hltesd)
274 : {
275 : // HLT QA on ESDs
276 0 : if (!esd || !hltesd) {
277 0 : AliError("invalid parameter: missing ESDs");
278 0 : return;
279 : }
280 :
281 : // make QA data from ESDs
282 :
283 0 : const Int_t nESDTracks = esd->GetNumberOfTracks();
284 0 : const Int_t nHLTesdTracks = hltesd->GetNumberOfTracks();
285 0 : FillESDsData(kMultiplicity,nESDTracks, nHLTesdTracks);
286 :
287 : Int_t nClsHLT = 0;
288 : Int_t nClsOffline = 0;
289 :
290 0 : for(Int_t i = 0; i < nESDTracks; i++){
291 0 : AliESDtrack * esdTrk = esd->GetTrack(i);
292 0 : FillESDsData(kPOffline,esdTrk->P());
293 0 : FillESDsData(kPtOffline,esdTrk->GetSignedPt());
294 0 : FillESDsData(kNClsPerTrkOffline,esdTrk->GetTPCNcls());
295 0 : FillESDsData(kPhiOffline,esdTrk->Phi()*TMath::RadToDeg());
296 0 : FillESDsData(kEtaOffline,esdTrk->Eta());
297 0 : nClsOffline+=esdTrk->GetTPCNcls();
298 : }
299 :
300 0 : for(Int_t i = 0; i < nHLTesdTracks; i++){
301 0 : AliESDtrack * hltEsdTrk = hltesd->GetTrack(i);
302 0 : FillESDsData(kPHLT,hltEsdTrk->P());
303 0 : FillESDsData(kPtHLT,hltEsdTrk->GetSignedPt());
304 0 : FillESDsData(kNClsPerTrkHLT,hltEsdTrk->GetTPCNcls());
305 0 : FillESDsData(kPhiHLT,hltEsdTrk->Phi()*TMath::RadToDeg());
306 0 : FillESDsData(kEtaHLT,hltEsdTrk->Eta());
307 0 : nClsHLT += hltEsdTrk->GetTPCNcls();
308 : }
309 :
310 0 : if(nESDTracks)
311 0 : nClsOffline /= nESDTracks;
312 0 : if(nHLTesdTracks)
313 0 : nClsHLT /= nHLTesdTracks;
314 :
315 0 : FillESDsData(kNCls,nClsOffline, nClsHLT);
316 :
317 :
318 0 : if(hltesd->IsHLTTriggerFired()){
319 0 : FillESDsData(kMultiplicityFired,nESDTracks, nHLTesdTracks);
320 0 : FillESDsData(kNClsFired,nClsOffline, nClsHLT);
321 :
322 0 : for(Int_t i = 0; i < nESDTracks; i++){
323 0 : AliESDtrack * esdTrk = esd->GetTrack(i);
324 0 : FillESDsData(kPOfflineFired,esdTrk->P());
325 0 : FillESDsData(kPtOfflineFired,esdTrk->GetSignedPt());
326 0 : FillESDsData(kNClsPerTrkOfflineFired,esdTrk->GetTPCNcls());
327 0 : FillESDsData(kPhiOfflineFired,esdTrk->Phi()*TMath::RadToDeg());
328 0 : FillESDsData(kEtaOfflineFired,esdTrk->Eta());
329 : }
330 :
331 0 : for(Int_t i = 0; i < nHLTesdTracks; i++){
332 0 : AliESDtrack * hltEsdTrk = hltesd->GetTrack(i);
333 0 : FillESDsData(kPHLTFired,hltEsdTrk->P());
334 0 : FillESDsData(kPtHLTFired,hltEsdTrk->GetSignedPt());
335 0 : FillESDsData(kNClsPerTrkHLTFired,hltEsdTrk->GetTPCNcls());
336 0 : FillESDsData(kPhiHLTFired,hltEsdTrk->Phi()*TMath::RadToDeg());
337 0 : FillESDsData(kEtaHLTFired,hltEsdTrk->Eta());
338 : }
339 :
340 :
341 0 : }
342 :
343 0 : }
|
