]> git.uio.no Git - u/mrichter/AliRoot.git/blame - STEER/AliGlobalQADataMaker.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Adding a few more histogramms.
[u/mrichter/AliRoot.git] / STEER / AliGlobalQADataMaker.cxx
CommitLineData
c0af1be2 1/*
2 The class for calculating the global (not detector specific) quality assurance.
3 It reuses the following TLists from its base class
4 AliQADataMaker::fRecPointsQAList (for keeping the track residuals)
5 AliQADataMaker::fESDsQAList (for keeping global ESD QA data)
6*/
7
542c6c71 8#include <TPDGCode.h>
c0af1be2 9#include <TH1F.h>
2937b8ff 10#include <TH2F.h>
c0af1be2 11
fb774831 12#include "AliQAChecker.h"
c0af1be2 13#include "AliGlobalQADataMaker.h"
14#include "AliGeomManager.h"
d7259b2e 15#include "AliESDEvent.h"
542c6c71 16#include "AliESDv0.h"
fb774831 17#include "AliRawReader.h"
2937b8ff 18#include "AliESDVZERO.h"
19#include "AliMultiplicity.h"
c0af1be2 20
21ClassImp(AliGlobalQADataMaker)
22
fb774831 23//____________________________________________________________________________
4e25ac79 24void AliGlobalQADataMaker::EndOfDetectorCycle(AliQAv1::TASKINDEX_t task, TObjArray ** list)
fb774831 25{
26 //Detector specific actions at end of cycle
27 // do the QA checking
4e25ac79 28 AliQAChecker::Instance()->Run(AliQAv1::kGLOBAL, task, list) ;
fb774831 29}
30
31//____________________________________________________________________________
32void AliGlobalQADataMaker::InitRaws()
33{
34 // create Raws histograms in Raws subdir
35}
36
37//____________________________________________________________________________
c0af1be2 38void AliGlobalQADataMaker::InitRecPoints() {
39 //------------------------------------------------------
d7259b2e 40 // This function books the histograms of *track*residuals*
c0af1be2 41 // as a part of global QA
42 //------------------------------------------------------
a6e0ebfe 43 const Char_t *name[]={
dd1b15c4 44 "hGlobalSPD1ResidualsY","SPD1ResidualsZ",
45 "hGlobalSPD2ResidualsY","SPD2ResidualsZ",
46 "hGlobalSDD1ResidualsY","SDD1ResidualsZ",
47 "hGlobalSDD2ResidualsY","SDD2ResidualsZ",
48 "hGlobalSSD1ResidualsY","SSD1ResidualsZ",
49 "hGlobalSSD2ResidualsY","SSD2ResidualsZ",
50
51 "hGlobalTPC1ResidualsY","TPC1ResidualsZ",
52 "hGlobalTPC2ResidualsY","TPC2ResidualsZ",
53
54 "hGlobalTRD1ResidualsY","TRD1ResidualsZ",
55 "hGlobalTRD2ResidualsY","TRD2ResidualsZ",
56 "hGlobalTRD3ResidualsY","TRD3ResidualsZ",
57 "hGlobalTRD4ResidualsY","TRD4ResidualsZ",
58 "hGlobalTRD5ResidualsY","TRD5ResidualsZ",
59 "hGlobalTRD6ResidualsY","TRD6ResidualsZ",
60
61 "hGlobalTOFResidualsY","TOFResidualsZ",
62
63 "hGlobalPHOS1ResidualsY","PHOS1ResidualsZ",
64 "hGlobalPHOS2ResidualsY","PHOS2ResidualsZ",
65
66 "hGlobalHMPIDResidualsY","HMPIDResidualsZ",
67
68 "hGlobalMUONResidualsY","MUONResidualsZ",
69
70 "hGlobalEMCALResidualsY","EMCALResidualsZ"
71 };
72 const Char_t *title[]={
c0af1be2 73 "SPD1 residuals Y","SPD1 residuals Z",
74 "SPD2 residuals Y","SPD2 residuals Z",
75 "SDD1 residuals Y","SDD1 residuals Z",
76 "SDD2 residuals Y","SDD2 residuals Z",
77 "SSD1 residuals Y","SSD1 residuals Z",
78 "SSD2 residuals Y","SSD2 residuals Z",
dd1b15c4 79
c0af1be2 80 "TPC1 residuals Y","TPC1 residuals Z",
81 "TPC2 residuals Y","TPC2 residuals Z",
dd1b15c4 82
c0af1be2 83 "TRD1 residuals Y","TRD1 residuals Z",
84 "TRD2 residuals Y","TRD2 residuals Z",
85 "TRD3 residuals Y","TRD3 residuals Z",
86 "TRD4 residuals Y","TRD4 residuals Z",
87 "TRD5 residuals Y","TRD5 residuals Z",
88 "TRD6 residuals Y","TRD6 residuals Z",
dd1b15c4 89
c0af1be2 90 "TOF residuals Y","TOF residuals Z",
dd1b15c4 91
c0af1be2 92 "PHOS1 residuals Y","PHOS1 residuals Z",
93 "PHOS2 residuals Y","PHOS2 residuals Z",
dd1b15c4 94
c0af1be2 95 "HMPID residuals Y","HMPID residuals Z",
dd1b15c4 96
c0af1be2 97 "MUON residuals Y","MUON residuals Z",
dd1b15c4 98
c0af1be2 99 "EMCAL residuals Y","EMCAL residuals Z"
100 };
dd1b15c4 101
c0af1be2 102 for (Int_t m=1; m<AliGeomManager::kLastLayer; m++) {
103 Int_t i=2*m-2;
dd1b15c4 104 TH1F *h=new TH1F(name[i],title[i],100,-5.,5.);
c0af1be2 105 Add2RecPointsList(h,i);
dd1b15c4 106 h=new TH1F(name[i+1],title[i+1],100,-5.,5.);
c0af1be2 107 Add2RecPointsList(h,i+1);
108 }
a935c6c0 109
49ab451d 110 Add2RecPointsList(
dd1b15c4 111 new TH1F("hGlobalSSD1AbsoluteResidualsYNegZ",
112 "SSD1 Absolute Residuals Y Neg Z",100,-2.,2.),40);
a935c6c0 113 Add2RecPointsList(
dd1b15c4 114 new TH1F("hGlobalSSD1AbsoluteResidualsZNegZ",
115 "SSD1 Absolute Residuals Z Neg Z",100,-2.,2.),41);
49ab451d 116 Add2RecPointsList(
dd1b15c4 117 new TH1F("hGlobalSSD1AbsoluteResidualsYPosZ",
118 "SSD1 Absolute Residuals Y Pos Z",100,-2.,2.),42);
a935c6c0 119 Add2RecPointsList(
dd1b15c4 120 new TH1F("hGlobalSSD1AbsoluteResidualsZPosZ",
121 "SSD1 Absolute Residuals Z Pos Z",100,-2.,2.),43);
49ab451d 122
123
124 Add2RecPointsList(
dd1b15c4 125 new TH1F("hGlobalSSD2AbsoluteResidualsYNegZ",
126 "SSD2 Absolute Residuals Y Neg Z",100,-3.,3.),44);
a935c6c0 127 Add2RecPointsList(
dd1b15c4 128 new TH1F("hGlobalSSD2AbsoluteResidualsZNegZ",
129 "SSD2 Absolute Residuals Z Neg Z",100,-3.,3.),45);
49ab451d 130 Add2RecPointsList(
dd1b15c4 131 new TH1F("hGlobalSSD2AbsoluteResidualsYPosZ",
132 "SSD2 Absolute Residuals Y Pos Z",100,-3.,3.),46);
a935c6c0 133 Add2RecPointsList(
dd1b15c4 134 new TH1F("hGlobalSSD2AbsoluteResidualsZPosZ",
135 "SSD2Absolute Residuals Z Pos Z",100,-3.,3.),47);
a935c6c0 136
c0af1be2 137}
d7259b2e 138
fb774831 139//____________________________________________________________________________
d7259b2e 140void AliGlobalQADataMaker::InitESDs() {
141 //------------------------------------------------------
142 // This function books the ESD QA histograms
143 // as a part of global QA
144 //------------------------------------------------------
542c6c71 145
f728e025 146 const Bool_t expert = kTRUE ;
147 const Bool_t image = kTRUE ;
148 {// Cluster related QA
dd1b15c4 149 const Char_t *name[]={
150 "hGlobalFractionAssignedClustersITS",
151 "hGlobalFractionAssignedClustersTPC",
2937b8ff 152 "hGlobalFractionAssignedClustersTRD",
153 "hGlobalClustersPerITSModule"
dd1b15c4 154 };
155 const Char_t *title[]={
156 "Fraction of the assigned clusters in ITS",
157 "Fraction of the assigned clusters in TPC",
2937b8ff 158 "Fraction of the assigned clusters in TRD",
159 "Number of clusters per an ITS module"
dd1b15c4 160 };
f728e025 161 Add2ESDsList(new TH1F(name[0],title[0],100,0.,2.),kClr0, !expert, image);
162 Add2ESDsList(new TH1F(name[1],title[1],100,0.,2.),kClr1, !expert, image);
163 Add2ESDsList(new TH1F(name[2],title[2],100,0.,2.),kClr2, !expert, image);
2937b8ff 164 Add2ESDsList(new TH1F(name[3],title[3],2240,0.,2240.),kClr3, !expert, image);
542c6c71 165 }
166
167 {// Track related QA
dd1b15c4 168 const Char_t *name[]={
169 "hGlobalTrackAzimuthe", // kTrk0
170 "hGlobalTrackEta", // kTrk1
171 "hGlobalTPCTrackpT", // kTrk2
172 "hGlobalTPCITSMatchedpT", // kTrk3
173 "hGlobalTPCTOFMatchedpT", // kTrk4
174 "hGlobalTPCITSMatchingProbability", // kTrk5
8336b059 175 "hGlobalTPCTOFMatchingProbability", // kTrk6
176 "hGlobalTPCsideAposDCA", // kTrk7
177 "hGlobalTPCsideAnegDCA", // kTrk8
178 "hGlobalTPCsideCposDCA", // kTrk9
179 "hGlobalTPCsideCnegDCA" // kTrk10
542c6c71 180 };
dd1b15c4 181 const Char_t *title[]={
182 "Track azimuthal distribution (rad)", // kTrk0
183 "Track pseudo-rapidity distribution", // kTrk1
184 "TPC: track momentum distribution (GeV)", // kTrk2
185 "TPC-ITS matched: track momentum distribution (GeV)", // kTrk3
186 "TPC-TOF matched: track momentum distribution (GeV)", // kTrk4
187 "TPC-ITS track-matching probability", // kTrk5
8336b059 188 "TPC-TOF track-matching probability", // kTrk6
189 "TPC side A: DCA for the positive tracks (mm)", // kTrk7
190 "TPC side A: DCA for the negative tracks (mm)", // kTrk8
191 "TPC side C: DCA for the positive tracks (mm)", // kTrk9
192 "TPC side C: DCA for the negative tracks (mm)" // kTrk10
dd1b15c4 193 };
f728e025 194 Add2ESDsList(new TH1F(name[0],title[0],100, 0.,TMath::TwoPi()),kTrk0, !expert, image);
195 Add2ESDsList(new TH1F(name[1],title[1],100,-2.00,2.00),kTrk1, !expert, image);
196 Add2ESDsList(new TH1F(name[2],title[2],50, 0.20,5.00),kTrk2, !expert, image);
197 Add2ESDsList(new TH1F(name[3],title[3],50, 0.20,5.00),kTrk3, !expert, image);
198 Add2ESDsList(new TH1F(name[4],title[4],50, 0.20,5.00),kTrk4, !expert, image);
199 Add2ESDsList(new TH1F(name[5],title[5],50, 0.20,5.00),kTrk5, !expert, image);
200 Add2ESDsList(new TH1F(name[6],title[6],50, 0.20,5.00),kTrk6, !expert, image);
201 Add2ESDsList(new TH1F(name[7],title[7],50, -25.0,25.0),kTrk7, !expert, image);
202 Add2ESDsList(new TH1F(name[8],title[8],50, -25.0,25.0),kTrk8, !expert, image);
203 Add2ESDsList(new TH1F(name[9],title[9],50, -25.0,25.0),kTrk9, !expert, image);
204 Add2ESDsList(new TH1F(name[10],title[10],50, -25.0,25.0),kTrk10, !expert, image);
542c6c71 205 }
206
207 {// V0 related QA
dd1b15c4 208 const Char_t *name[]={
209 "hGlobalPromptK0sMass",
210 "hGlobalOfflineK0sMass",
211 "hGlobalPromptLambda0Lambda0BarMass",
212 "hGlobalOfflineLambda0Lambda0BarMass"
213 };
214 const Char_t *title[]={
215 "On-the-fly K0s mass (GeV)",
216 "Offline K0s mass (GeV)",
217 "On-the-fly Lambda0 + Lambda0Bar mass (GeV)",
218 "Offline Lambda0 + Lambda0Bar mass (GeV)"
219 };
f728e025 220 Add2ESDsList(new TH1F(name[0],title[0],50, 0.4477,0.5477),kK0on, !expert, image);
221 Add2ESDsList(new TH1F(name[1],title[1],50, 0.4477,0.5477),kK0off, !expert, image);
222 Add2ESDsList(new TH1F(name[2],title[2],50, 1.0657,1.1657),kL0on, !expert, image);
223 Add2ESDsList(new TH1F(name[3],title[3],50, 1.0657,1.1657),kL0off, !expert, image);
542c6c71 224 }
225
226 {// PID related QA
a6e0ebfe 227 const Char_t *name[]={
c9622234 228 "hGlobalITSdEdx",
229 "hGlobalTPCdEdx",
2937b8ff 230 "hGlobalTOFTrackingvsMeasured",
231 "hGlobalTPCdEdxvsMomentum"
232 };
dd1b15c4 233 const Char_t *title[]={
234 "ITS: dEdx (ADC) for particles with momentum 0.4 - 0.5 (GeV)",
235 "TPC: dEdx (ADC) for particles with momentum 0.4 - 0.5 (GeV)",
2937b8ff 236 "TOF: tracking - measured (ps)",
237 "TPC: dEdx (A.U.) vs momentum (GeV)"
238 };
f728e025 239 Add2ESDsList(new TH1F(name[0],title[0],50,0.00,200.),kPid0, !expert, image);
240 Add2ESDsList(new TH1F(name[1],title[1],50,0.00,100.),kPid1, !expert, image);
241 Add2ESDsList(new TH1F(name[2],title[2],50,-3500.,3500.),kPid2, !expert, image);
2937b8ff 242 Add2ESDsList(new TH2F(name[3],title[3],1500,0.05,15.,700,0.,700.),kPid3,!expert,image);
243 }
244 {// Multiplicity related QA
245 const Char_t *name[]={
246 "hGlobalV0AvsITS",
247 "hGlobalV0CvsITS"
248 };
249 const Char_t *title[]={
250 "Multiplicity: V0A vs ITS",
251 "Multiplicity: V0C vs ITS"
252 };
253 TH2F *h0=new TH2F(name[0],title[0],40,0.,40., 32,0.,32.);
254 Add2ESDsList(h0,kMlt0, !expert, image);
255 TH2F *h1=new TH2F(name[1],title[1],40,0.,40., 32,0.,32.);
256 Add2ESDsList(h1,kMlt1, !expert, image);
542c6c71 257 }
258
d7259b2e 259}
260
fb774831 261//____________________________________________________________________________
262void AliGlobalQADataMaker::MakeRaws(AliRawReader* rawReader)
263{
264 //Fill prepared histograms with Raw digit properties
265 rawReader->Reset() ;
266
267}
268
269//____________________________________________________________________________
542c6c71 270void AliGlobalQADataMaker::MakeESDs(AliESDEvent * event) {
d7259b2e 271 //-----------------------------------------------------------
272 // This function fills the ESD QA histograms
273 // as a part of global QA
274 //-----------------------------------------------------------
eca4fa66 275 // Check id histograms already created for this Event Specie
276 if ( ! GetESDsData(kClr0) )
277 InitESDs() ;
278
c9622234 279
542c6c71 280 const AliESDEvent *esd=event;
281
d7259b2e 282 Int_t ntrk=esd->GetNumberOfTracks() ;
283 for (Int_t i=0; i<ntrk; i++) {
542c6c71 284 const AliESDtrack *track=esd->GetTrack(i);
d7259b2e 285
542c6c71 286 // Cluster related QA
d7259b2e 287 if (track->IsOn(AliESDtrack::kITSrefit)) {
288 Int_t n=track->GetITSclusters(0);
542c6c71 289 GetESDsData(kClr0)->Fill(Float_t(n)/6.); //6 is the number of ITS layers
d7259b2e 290 }
291
2937b8ff 292 for (Int_t i=0; i<6; i++) {
293 Int_t idet, sts;
294 Float_t xloc,zloc;
295 if (!track->GetITSModuleIndexInfo(i,idet,sts,xloc,zloc)) continue;
296 if (i>=2) idet+=240;
297 if (i>=4) idet+=260;
298 if ((sts==1)||(sts==2)||(sts==4)) GetESDsData(kClr3)->Fill(idet);
299 }
300
d7259b2e 301 if (track->IsOn(AliESDtrack::kTPCrefit)) {
302 Int_t n =track->GetTPCNcls();
303 Int_t nf=track->GetTPCNclsF(); // number of crossed TPC pad rows
1aaf4118 304 if (nf>0) {
305 Double_t val = n*1.0/nf;
306 GetESDsData(kClr1)->Fill(val);
307 }
d7259b2e 308 }
309
310 if (track->IsOn(AliESDtrack::kTRDrefit)) {
311 Int_t n=track->GetTRDclusters(0);
e3bd5382 312 GetESDsData(kClr2)->Fill(Float_t(n)/(6*24));//(6*24) is the number of TRD time bins
542c6c71 313 }
314
315 Double_t p=track->GetP();
316
317 // Track related QA
318 if (track->IsOn(AliESDtrack::kTPCrefit)) {
319 Float_t dz[2];
320 track->GetDZ(0.,0.,0.,esd->GetMagneticField(),dz);
321 if ((TMath::Abs(dz[0])<3.) && (TMath::Abs(dz[1])<3.)) { // beam pipe
322 Double_t phi=track->Phi();
323 GetESDsData(kTrk0)->Fill(phi);
44f14df0 324 Double_t y=track->Eta();
542c6c71 325 GetESDsData(kTrk1)->Fill(y);
326
cecff4fc 327 if (TMath::Abs(y)<0.9) {
328 GetESDsData(kTrk2)->Fill(p);
329 if (track->IsOn(AliESDtrack::kITSrefit)) GetESDsData(kTrk3)->Fill(p);
4d9ba05c 330 //if (track->IsOn(AliESDtrack::kTOFout)) GetESDsData(kTrk4)->Fill(p);
331 if (track->GetTOFsignal()>0) GetESDsData(kTrk4)->Fill(p);
cecff4fc 332 }
542c6c71 333 }
334 }
8336b059 335 const AliExternalTrackParam *tpcTrack=track->GetTPCInnerParam();
336 const AliExternalTrackParam *innTrack=track->GetInnerParam();
337 if (tpcTrack)
338 if (innTrack) {
339 const AliESDVertex *vtx=esd->GetPrimaryVertex();
340 Double_t xv=vtx->GetXv();
341 Double_t yv=vtx->GetYv();
342 Double_t zv=vtx->GetZv();
343 Float_t dz[2];
344 tpcTrack->GetDZ(xv,yv,zv,esd->GetMagneticField(),dz);
345 dz[0]*=10.; // in mm
346 if (innTrack->GetZ() > 0)
347 if (innTrack->GetTgl()> 0) { // TPC side A
348 if (tpcTrack->GetSign() > 0) GetESDsData(kTrk7)->Fill(dz[0]);
349 else GetESDsData(kTrk8)->Fill(dz[0]);
350 }
351 if (innTrack->GetZ() < 0)
352 if (innTrack->GetTgl()< 0) { // TPC side C
353 if (tpcTrack->GetSign() > 0) GetESDsData(kTrk9)->Fill(dz[0]);
354 else GetESDsData(kTrk10)->Fill(dz[0]);
355 }
356 }
542c6c71 357
358 // PID related QA
359 if ((p>0.4) && (p<0.5)) {
360 if (track->IsOn(AliESDtrack::kITSpid)) {
361 Double_t dedx=track->GetITSsignal();
362 GetESDsData(kPid0)->Fill(dedx);
363 }
364 if (track->IsOn(AliESDtrack::kTPCpid)) {
365 Double_t dedx=track->GetTPCsignal();
366 GetESDsData(kPid1)->Fill(dedx);
367 }
d7259b2e 368 }
542c6c71 369 if (p>1.0) {
370 if (track->IsOn(AliESDtrack::kTOFpid)) {
371 Double_t times[10];
372 track->GetIntegratedTimes(times);
373 Double_t tof=track->GetTOFsignal();
374 GetESDsData(kPid2)->Fill(times[2]-tof);
375 }
376 }
2937b8ff 377 const AliExternalTrackParam *par=track->GetInnerParam();
378 if (par) {
379 Double_t pp=par->GetP();
380 Double_t dedx=track->GetTPCsignal();
381 TH2F *h = dynamic_cast<TH2F*>(GetESDsData(kPid3));
382 h->Fill(pp,dedx);
383 }
384
542c6c71 385 }
d7259b2e 386
2937b8ff 387 // Multiplicity related QA
388 AliESDVZERO *mltV0 =esd->GetVZEROData();
389 const AliMultiplicity *mltITS=esd->GetMultiplicity();
390 if (mltV0)
391 if (mltITS) {
392 Short_t nv0a=mltV0->GetNbPMV0A();
393 Short_t nv0c=mltV0->GetNbPMV0C();
394 Int_t nits=mltITS->GetNumberOfTracklets();
395 TH2F *h0=dynamic_cast<TH2F*>(GetESDsData(kMlt0));
396 h0->Fill(nits,nv0a);
397 TH2F *h1=dynamic_cast<TH2F*>(GetESDsData(kMlt1));
398 h1->Fill(nits,nv0c);
399 }
400
401
cecff4fc 402 TH1 *tpc=GetESDsData(kTrk2); tpc->Sumw2();
403 TH1 *its=GetESDsData(kTrk3); its->Sumw2();
404 TH1 *tof=GetESDsData(kTrk4); tof->Sumw2();
542c6c71 405 GetESDsData(kTrk5)->Divide(its,tpc,1,1.,"b");
406 GetESDsData(kTrk6)->Divide(tof,tpc,1,1.,"b");
407
408 // V0 related QA
409 Int_t nV0=esd->GetNumberOfV0s();
410 for (Int_t i=0; i<nV0; i++) {
411 Double_t mass;
412 AliESDv0 v0(*esd->GetV0(i));
413
414 v0.ChangeMassHypothesis(kK0Short);
415 mass=v0.GetEffMass();
b72280b3 416 if (v0.GetOnFlyStatus())
417 GetESDsData(kK0on)->Fill(mass);
418 else
419 GetESDsData(kK0off)->Fill(mass);
542c6c71 420
421 v0.ChangeMassHypothesis(kLambda0);
422 mass=v0.GetEffMass();
b72280b3 423 if (v0.GetOnFlyStatus())
424 GetESDsData(kL0on)->Fill(mass);
425 else
426 GetESDsData(kL0off)->Fill(mass);
542c6c71 427
428 v0.ChangeMassHypothesis(kLambda0Bar);
429 mass=v0.GetEffMass();
b72280b3 430 if (v0.GetOnFlyStatus())
431 GetESDsData(kL0on)->Fill(mass);
432 else
433 GetESDsData(kL0off)->Fill(mass);
d7259b2e 434 }
542c6c71 435
f92b626a 436}