]> git.uio.no Git - u/mrichter/AliRoot.git/blame - TRD/qaRec/AliTRDtrackingResolution.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
new task introduced by Alex Wilk (calculation of reference data for NN
[u/mrichter/AliRoot.git] / TRD / qaRec / AliTRDtrackingResolution.cxx
CommitLineData
77203477 1/**************************************************************************
d2381af5 2* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
3* *
4* Author: The ALICE Off-line Project. *
5* Contributors are mentioned in the code where appropriate. *
6* *
7* Permission to use, copy, modify and distribute this software and its *
8* documentation strictly for non-commercialf purposes is hereby granted *
9* without fee, provided that the above copyright notice appears in all *
10* copies and that both the copyright notice and this permission notice *
11* appear in the supporting documentation. The authors make no claims *
12* about the suitability of this software for any purpose. It is *
13* provided "as is" without express or implied warranty. *
14**************************************************************************/
77203477 15
16/* $Id: AliTRDtrackingResolution.cxx 27496 2008-07-22 08:35:45Z cblume $ */
17
18////////////////////////////////////////////////////////////////////////////
19// //
20// Reconstruction QA //
21// //
22// Authors: //
23// Markus Fasel <M.Fasel@gsi.de> //
24// //
25////////////////////////////////////////////////////////////////////////////
26
aaf47b30 27#include <cstring>
28
29
77203477 30#include <TObjArray.h>
7102d1b1 31#include <TH2.h>
32#include <TH1.h>
33#include <TF1.h>
77203477 34#include <TProfile.h>
7102d1b1 35#include <TGraphErrors.h>
77203477 36#include <TMath.h>
aaf47b30 37#include "TTreeStream.h"
38#include "TGeoManager.h"
77203477 39
40#include "AliAnalysisManager.h"
77203477 41#include "AliTrackReference.h"
aaf47b30 42#include "AliTrackPointArray.h"
43#include "AliCDBManager.h"
44
45#include "AliTRDcluster.h"
46#include "AliTRDseedV1.h"
47#include "AliTRDtrackV1.h"
48#include "AliTRDtrackerV1.h"
49#include "AliTRDReconstructor.h"
50#include "AliTRDrecoParam.h"
77203477 51
52#include "AliTRDtrackInfo/AliTRDtrackInfo.h"
53#include "AliTRDtrackingResolution.h"
54
77203477 55ClassImp(AliTRDtrackingResolution)
56
57//________________________________________________________
3d86166d 58AliTRDtrackingResolution::AliTRDtrackingResolution()
59 :AliTRDrecoTask("Resolution", "Tracking Resolution")
aaf47b30 60 ,fReconstructor(0x0)
77203477 61{
aaf47b30 62 fReconstructor = new AliTRDReconstructor();
63 fReconstructor->SetRecoParam(AliTRDrecoParam::GetLowFluxParam());
77203477 64}
65
ed383798 66//________________________________________________________
67AliTRDtrackingResolution::~AliTRDtrackingResolution()
68{
ed383798 69 delete fReconstructor;
2c0cf367 70 if(gGeoManager) delete gGeoManager;
ed383798 71}
72
77203477 73
74//________________________________________________________
75void AliTRDtrackingResolution::CreateOutputObjects()
76{
77 // spatial resolution
77203477 78 OpenFile(0, "RECREATE");
39779ce6 79
3d86166d 80 fContainer = Histos();
874acced 81
82 // cluster to tracklet residuals [2]
765bd0ab 83 fContainer->AddAt(new TH2I("fYClRes", "Clusters Residuals", 21, -21., 21., 100, -.5, .5), kClusterYResidual);
874acced 84 // tracklet to Riemann fit residuals [2]
765bd0ab 85 fContainer->AddAt(new TH2I("fYTrkltRRes", "Tracklet Riemann Residuals", 21, -21., 21., 100, -.5, .5), kTrackletRiemanYResidual);
86 fContainer->AddAt(new TH2I("fAngleTrkltRRes", "Tracklet Riemann Angular Residuals", 21, -21., 21., 100, -.5, .5), kTrackletRiemanAngleResidual);
87 fContainer->AddAt(new TH2I("fYTrkltKRes", "Tracklet Kalman Residuals", 21, -21., 21., 100, -.5, .5), kTrackletKalmanYResidual);
88 fContainer->AddAt(new TH2I("fAngleTrkltKRes", "Tracklet Kalman Angular Residuals", 21, -21., 21., 100, -.5, .5), kTrackletKalmanAngleResidual);
39779ce6 89
90 // Resolution histos
4b8f8a35 91 if(HasMCdata()){
d2381af5 92 // tracklet resolution [0]
765bd0ab 93 fContainer->AddAt(new TH2I("fY", "Tracklet Resolution", 21, -21., 21., 100, -.5, .5), kTrackletYResolution);
d2381af5 94 // tracklet angular resolution [1]
765bd0ab 95 fContainer->AddAt(new TH2I("fPhi", "Tracklet Angular Resolution", 21, -21., 21., 100, -10., 10.), kTrackletAngleResolution);
874acced 96
97 // Riemann track resolution [y, z, angular]
765bd0ab 98 fContainer->AddAt(new TH2I("fYRT", "Track Riemann Y Resolution", 21, -21., 21., 100, -.5, .5), kTrackRYResolution);
99 fContainer->AddAt(new TH2I("fZRT", "Track Riemann Z Resolution", 21, -21., 21., 100, -.5, .5), kTrackRZResolution);
100 fContainer->AddAt(new TH2I("fPhiRT", "Track Riemann Angular Resolution", 21, -21., 21., 100, -10., 10.), kTrackRAngleResolution);
874acced 101
102 // Kalman track resolution [y, z, angular]
3d86166d 103 fContainer->AddAt(new TH2I("fYKT", "", 21, -21., 21., 100, -.5, .5), kTrackKYResolution);
104 fContainer->AddAt(new TH2I("fZKT", "", 21, -21., 21., 100, -.5, .5), kTrackKZResolution);
105 fContainer->AddAt(new TH2I("fPhiKT", "", 21, -21., 21., 100, -10., 10.), kTrackKAngleResolution);
d2381af5 106 }
d85cd79c 107
108 // CREATE GRAPHS for DISPLAY
109
110 // define iterator over graphs
111 Int_t jgraph = (Int_t)kGraphStart;
112 TH2I *h2 = (TH2I *)(fContainer->At(kClusterYResidual));
113 // clusters tracklet residuals (mean-phi)
765bd0ab 114 TH1 *h = new TH1I("h", "", 100, -40., 40.);
115 h->GetXaxis()->SetTitle("#Phi [deg]");
116 h->GetYaxis()->SetTitle("Clusters Residuals : #sigma/#mu [mm]");
117 h->GetYaxis()->SetRangeUser(-.05, 1.);
118 fContainer->AddAt(h, jgraph++);
119
d85cd79c 120 TGraphErrors *g = new TGraphErrors(h2->GetNbinsX());
121 g->SetLineColor(kGreen);
122 g->SetMarkerStyle(22);
123 g->SetMarkerColor(kGreen);
124 g->SetNameTitle("clm", "Residuals Clusters-Tracklet Mean");
125 fContainer->AddAt(g, jgraph++);
126
127 // clusters tracklet residuals (sigma-phi)
128 g = new TGraphErrors(h2->GetNbinsX());
129 g->SetLineColor(kRed);
130 g->SetMarkerStyle(23);
131 g->SetMarkerColor(kRed);
132 g->SetNameTitle("cls", "Residuals Clusters-Tracklet Sigma");
133 fContainer->AddAt(g, jgraph++);
134
135 if(HasMCdata()){
136 // tracklet y resolution
137 h2 = (TH2I*)fContainer->At(kTrackletYResolution);
765bd0ab 138 h = new TH1I("h", "", 100, -40., 40.);
139 h->GetXaxis()->SetTitle("#Phi [deg]");
140 h->GetYaxis()->SetTitle("Tracklet Resolution : #sigma/#mu [mm]");
141 h->GetYaxis()->SetRangeUser(-.05, 1.);
142 fContainer->AddAt(h, jgraph++);
143
d85cd79c 144 g = new TGraphErrors(h2->GetNbinsX());
145 g->SetLineColor(kGreen);
146 g->SetMarkerStyle(22);
147 g->SetMarkerColor(kGreen);
148 g->SetNameTitle("trkltym", "Resolution Tracklet Y Mean");
149 fContainer->AddAt(g, jgraph++);
150 g = new TGraphErrors(h2->GetNbinsX());
151 g->SetLineColor(kRed);
152 g->SetMarkerStyle(22);
153 g->SetMarkerColor(kRed);
154 g->SetNameTitle("trkltys", "Resolution Tracklet Y Sigma");
155 fContainer->AddAt(g, jgraph++);
156
157 // tracklet phi resolution
158 h2 = (TH2I*)fContainer->At(kTrackletAngleResolution);
765bd0ab 159 h = new TH1I("h", "", 100, -40., 40.);
160 h->GetXaxis()->SetTitle("#Phi [deg]");
161 h->GetYaxis()->SetTitle("Tracklet Angular Resolution : #sigma/#mu [deg]");
28efdace 162 h->GetYaxis()->SetRangeUser(-.05, .2);
765bd0ab 163 fContainer->AddAt(h, jgraph++);
164
d85cd79c 165 g = new TGraphErrors(h2->GetNbinsX());
166 g->SetLineColor(kGreen);
167 g->SetMarkerStyle(22);
168 g->SetMarkerColor(kGreen);
169 g->SetNameTitle("trkltam", "Resolution Tracklet Y Mean");
170 fContainer->AddAt(g, jgraph++);
171 g = new TGraphErrors(h2->GetNbinsX());
172 g->SetLineColor(kRed);
173 g->SetMarkerStyle(22);
174 g->SetMarkerColor(kRed);
175 g->SetNameTitle("trkltas", "Angle Resolution Tracklet Sigma");
176 fContainer->AddAt(g, jgraph++);
177 }
77203477 178}
179
180//________________________________________________________
7102d1b1 181void AliTRDtrackingResolution::Exec(Option_t *)
182{
77203477 183 // spatial Resolution: res = pos_{Tracklet}(x = x_{Anode wire}) - pos_{TrackRef}(x = x_{Anode wire})
184 // angular Resolution: res = Tracklet angle - TrackRef Angle
7102d1b1 185
39779ce6 186 Int_t nTrackInfos = fTracks->GetEntriesFast();
874acced 187 if(fDebugLevel>=2) printf("Number of Histograms: %d\n", Histos()->GetEntries());
aaf47b30 188
765bd0ab 189 Double_t p, dy, dphi, dymc, dzmc, dphimc;
190 Float_t fP[kNLayers], fY[kNLayers], fZ[kNLayers], fPhi[kNLayers], fTheta[kNLayers]; // phi/theta angle per layer
191 Bool_t fMCMap[kNLayers], fLayerMap[kNLayers]; // layer map
192
39779ce6 193 AliTrackPoint tr[kNLayers], tk[kNLayers];
765bd0ab 194 AliExternalTrackParam *fOp = 0x0;
aaf47b30 195 AliTRDtrackV1 *fTrack = 0x0;
77203477 196 AliTRDtrackInfo *fInfo = 0x0;
197 if(fDebugLevel>=2) printf("Number of TrackInfos: %d\n", nTrackInfos);
198 for(Int_t iTI = 0; iTI < nTrackInfos; iTI++){
77203477 199 // check if ESD and MC-Information are available
39779ce6 200 if(!(fInfo = dynamic_cast<AliTRDtrackInfo *>(fTracks->UncheckedAt(iTI)))) continue;
aaf47b30 201 if(!(fTrack = fInfo->GetTRDtrack())) continue;
765bd0ab 202 if(!(fOp = fInfo->GetOuterParam())) continue;
203
7102d1b1 204 if(fDebugLevel>=3) printf("\tDoing track[%d] NTrackRefs[%d]\n", iTI, fInfo->GetNTrackRefs());
205
765bd0ab 206 p = fOp->P();
39779ce6 207 Int_t npts = 0;
765bd0ab 208 memset(fP, 0, kNLayers*sizeof(Float_t));
209 memset(fY, 0, kNLayers*sizeof(Float_t));
210 memset(fZ, 0, kNLayers*sizeof(Float_t));
211 memset(fPhi, 0, kNLayers*sizeof(Float_t));
212 memset(fTheta, 0, kNLayers*sizeof(Float_t));
213 memset(fLayerMap, 0, kNLayers*sizeof(Bool_t));
214 memset(fMCMap, 0, kNLayers*sizeof(Bool_t));
aaf47b30 215 AliTRDseedV1 *fTracklet = 0x0;
77203477 216 for(Int_t iplane = 0; iplane < kNLayers; iplane++){
aaf47b30 217 if(!(fTracklet = fTrack->GetTracklet(iplane))) continue;
218 if(!fTracklet->IsOK()) continue;
219
39779ce6 220 // Book point arrays
765bd0ab 221 fLayerMap[iplane] = kTRUE;
39779ce6 222 tr[npts].SetXYZ(fTracklet->GetX0(), 0., 0.);
223 tk[npts].SetXYZ(fTracklet->GetX0(), fTracklet->GetYfit(0), fTracklet->GetZfit(0));
224 npts++;
225
7102d1b1 226 if(fDebugLevel>=4) printf("\t\tLy[%d] X0[%6.3f] Ncl[%d]\n", iplane, fTracklet->GetX0(), fTracklet->GetN());
227
765bd0ab 228 // define reference values
229 fP[iplane] = p;
230 fY[iplane] = fTracklet->GetYref(0);
231 fZ[iplane] = fTracklet->GetZref(0);
232 fPhi[iplane] = TMath::ATan(fTracklet->GetYref(1));
233 fTheta[iplane] = TMath::ATan(fTracklet->GetZref(1));
234
7102d1b1 235
39779ce6 236 // RESOLUTION (compare to MC)
4b8f8a35 237 if(HasMCdata()){
d2381af5 238 if(fInfo->GetNTrackRefs() >= 2){
765bd0ab 239 Double_t pmc, ymc, zmc, phiMC, thetaMC;
240 if(Resolution(fTracklet, fInfo, pmc, ymc, zmc, phiMC, thetaMC)){
241 fMCMap[iplane] = kTRUE;
242 fP[iplane] = pmc;
243 fY[iplane] = ymc;
244 fZ[iplane] = zmc;
245 fPhi[iplane] = phiMC;
246 fTheta[iplane] = thetaMC;
247 }
d2381af5 248 }
249 }
765bd0ab 250 Float_t phi = fPhi[iplane]*TMath::RadToDeg();
251 Float_t theta = fTheta[iplane]*TMath::RadToDeg();
aaf47b30 252
39779ce6 253 // Do clusters residuals
d2381af5 254 if(!fTracklet->Fit(kFALSE)) continue;
aaf47b30 255 AliTRDcluster *c = 0x0;
256 for(Int_t ic=AliTRDseed::knTimebins-1; ic>=0; ic--){
257 if(!(c = fTracklet->GetClusters(ic))) continue;
39779ce6 258
259 dy = fTracklet->GetYat(c->GetX()) - c->GetY();
765bd0ab 260 ((TH2I*)fContainer->At(kClusterYResidual))->Fill(phi, dy);
3d86166d 261 if(fDebugLevel>=2){
39779ce6 262 Float_t q = c->GetQ();
765bd0ab 263 (*fDebugStream) << "ResidualClusters"
264 << "ly=" << iplane
265 << "phi=" << phi
266 << "tht=" << theta
267 << "q=" << q
268 << "dy=" << dy
39779ce6 269 << "\n";
270 }
aaf47b30 271 }
272 }
39779ce6 273
d2381af5 274
39779ce6 275 // this protection we might drop TODO
276 if(fTrack->GetNumberOfTracklets() < 6) continue;
277
278 AliTRDtrackerV1::FitRiemanTilt(fTrack, 0x0, kTRUE, npts, tr);
765bd0ab 279 Int_t iref = 0;
280 for(Int_t ip=0; ip<kNLayers; ip++){
281 if(!fLayerMap[ip]) continue;
282 fTracklet = fTrack->GetTracklet(ip);
283 // recalculate fit based on the new tilt correction
284 fTracklet->Fit();
874acced 285
765bd0ab 286 dy = fTracklet->GetYfit(0) - tr[iref].GetY();
287 ((TH2I*)fContainer->At(kTrackletRiemanYResidual))->Fill(fPhi[ip]*TMath::RadToDeg(), dy);
288
289 dphi = fTracklet->GetYfit(1)- fTracklet->GetYref(1);
290 ((TH2I*)fContainer->At(kTrackletRiemanAngleResidual))->Fill(fPhi[ip]*TMath::RadToDeg(), dphi);
291
292 if(HasMCdata()){
293 dymc = fY[ip] - tr[iref].GetY();
294 ((TH2I*)fContainer->At(kTrackRYResolution))->Fill(fPhi[ip]*TMath::RadToDeg(), dymc);
295
296 dzmc = fZ[ip] - tr[iref].GetZ();
297 ((TH2I*)fContainer->At(kTrackRZResolution))->Fill(fPhi[ip]*TMath::RadToDeg(), dzmc);
298
299 dphimc = fPhi[ip] - fTracklet->GetYfit(1);
300 ((TH2I*)fContainer->At(kTrackRAngleResolution))->Fill(fPhi[ip]*TMath::RadToDeg(), dphimc);
301 }
874acced 302
765bd0ab 303 iref++;
874acced 304
3d86166d 305 if(fDebugLevel>=2){
765bd0ab 306 (*fDebugStream) << "RiemannTrack"
307 << "ly=" << ip
308 << "mc=" << fMCMap[ip]
309 << "p=" << fP[ip]
310 << "phi=" << fPhi[ip]
311 << "tht=" << fTheta[ip]
312 << "dy=" << dy
313 << "dphi=" << dphi
314 << "dymc=" << dymc
315 << "dzmc=" << dzmc
316 << "dphimc="<< dphimc
aaf47b30 317 << "\n";
318 }
319 }
320
2c0cf367 321// if(!gGeoManager) TGeoManager::Import("geometry.root");
aaf47b30 322// AliTRDtrackerV1::FitKalman(fTrack, 0x0, kFALSE, nc, tr);
323// for(Int_t ip=0; ip<nc; ip++){
324// dy = cl[ip].GetY() - tr[ip].GetY();
3d86166d 325// ((TH2I*)fContainer->At(kTrackletKalmanYResidual))->Fill(phi*TMath::RadToDeg(), dy);
aaf47b30 326// dz = cl[ip].GetZ() - tr[ip].GetZ();
327// if(fDebugLevel>=1){
765bd0ab 328// (*fDebugStream) << "KalmanTrack"
aaf47b30 329// << "dy=" << dy
330// << "dz=" << dz
331// /* << "phi=" << phi
332// << "theta=" << theta
333// << "dphi=" << dphi*/
334// << "\n";
335// }
336// }
39779ce6 337
338
77203477 339 }
3d86166d 340 PostData(0, fContainer);
77203477 341}
342
d85cd79c 343//________________________________________________________
28efdace 344void AliTRDtrackingResolution::GetRefFigure(Int_t ifig, Int_t &first, Int_t &last, Option_t */*opt*/)
d85cd79c 345{
346 switch(ifig){
347 case 0:
765bd0ab 348 first = (Int_t)kGraphStart; last = first+3;
d85cd79c 349 break;
350 case 1:
765bd0ab 351 first = (Int_t)kGraphStart+3; last = first+3;
d85cd79c 352 break;
353 case 2:
765bd0ab 354 first = (Int_t)kGraphStart+6; last = first+3;
d85cd79c 355 break;
356 default:
357 first = (Int_t)kGraphStart; last = first;
358 break;
359 }
360}
361
39779ce6 362
363//________________________________________________________
765bd0ab 364Bool_t AliTRDtrackingResolution::Resolution(AliTRDseedV1 *tracklet, AliTRDtrackInfo *fInfo, Double_t &p, Double_t &ymc, Double_t &zmc, Double_t &phi, Double_t &theta)
39779ce6 365{
366
367 AliTrackReference *fTrackRefs[2] = {0x0, 0x0}, *tempTrackRef = 0x0;
368
369 // check for 2 track ref where the radial position has a distance less than 3.7mm
370 Int_t nFound = 0;
371 for(Int_t itr = 0; itr < fInfo->GetNTrackRefs(); itr++){
372 if(!(tempTrackRef = fInfo->GetTrackRef(itr))) continue;
373 if(fDebugLevel>=5) printf("TrackRef %d: x = %f\n", itr, tempTrackRef->LocalX());
374 if(TMath::Abs(tracklet->GetX0() - tempTrackRef->LocalX()) > 3.7) continue;
375 fTrackRefs[nFound++] = tempTrackRef;
376 if(nFound == 2) break;
377 }
378 if(nFound < 2){
379 if(fDebugLevel>=4) printf("\t\tFound track crossing [%d] refX[%6.3f]\n", nFound, nFound>0 ? fTrackRefs[0]->LocalX() : 0.);
380 return kFALSE;
381 }
382 // We found 2 track refs for the tracklet, get y and z at the anode wire by a linear approximation
383
384
385 // RESOLUTION
386 Double_t dx = fTrackRefs[1]->LocalX() - fTrackRefs[0]->LocalX();
387 if(dx <= 0.){
388 if(fDebugLevel>=4) printf("\t\ttrack ref in the wrong order refX0[%6.3f] refX1[%6.3f]\n", fTrackRefs[0]->LocalX(), fTrackRefs[1]->LocalX());
389 return kFALSE;
390 }
391 Double_t dydx = (fTrackRefs[1]->LocalY() - fTrackRefs[0]->LocalY()) / dx;
392 Double_t dzdx = (fTrackRefs[1]->Z() - fTrackRefs[0]->Z()) / dx;
393 Double_t dx0 = fTrackRefs[1]->LocalX() - tracklet->GetX0();
765bd0ab 394 ymc = fTrackRefs[1]->LocalY() - dydx*dx0;
395 zmc = fTrackRefs[1]->Z() - dzdx*dx0;
39779ce6 396
397 // recalculate tracklet based on the MC info
84b3ecda 398 tracklet->SetZref(0, zmc);
765bd0ab 399 tracklet->SetZref(1, -dzdx); // TODO
39779ce6 400 tracklet->Fit();
401 Double_t dy = tracklet->GetYfit(0) - ymc;
402 Double_t dz = tracklet->GetZfit(0) - zmc;
403
404 //res_y *= 100; // in mm
765bd0ab 405 p = fTrackRefs[0]->P();
39779ce6 406
407 phi = TMath::ATan(dydx);
765bd0ab 408 theta = TMath::ATan(dzdx);
39779ce6 409 Double_t dphi = TMath::ATan(tracklet->GetYfit(1)) - phi;
410 if(fDebugLevel>=4) printf("\t\tdx[%6.4f] dy[%6.4f] dz[%6.4f] dphi[%6.4f] \n", dx, dy, dz, dphi);
411
412 // Fill Histograms
413 if(TMath::Abs(dx-3.7)<1.E-3){
3d86166d 414 ((TH2I*)fContainer->At(kTrackletYResolution))->Fill(phi*TMath::RadToDeg(), dy);
415 ((TH2I*)fContainer->At(kTrackletAngleResolution))->Fill(phi*TMath::RadToDeg(), dphi*TMath::RadToDeg());
39779ce6 416 }
417 // Fill Debug Tree
3d86166d 418 if(fDebugLevel>=2){
39779ce6 419 Int_t iplane = tracklet->GetPlane();
765bd0ab 420 (*fDebugStream) << "ResolutionTrklt"
421 << "ly=" << iplane
422 << "p=" << p
423 << "phi=" << phi
424 << "tht=" << theta
425 << "ymc=" << ymc
426 << "zmc=" << zmc
39779ce6 427 << "dx=" << dx
428 << "dy=" << dy
429 << "dz=" << dz
39779ce6 430 << "dphi=" << dphi
84b3ecda 431 << "tracklet.="<<tracklet
39779ce6 432 << "\n";
433 }
434
435 return kTRUE;
436}
437
77203477 438//________________________________________________________
d85cd79c 439Bool_t AliTRDtrackingResolution::PostProcess()
7102d1b1 440{
d85cd79c 441 //fContainer = dynamic_cast<TObjArray*>(GetOutputData(0));
765bd0ab 442 fNRefFigures = 0;
d85cd79c 443 if (!fContainer) {
444 Printf("ERROR: list not available");
445 return kFALSE;
3d86166d 446 }
7102d1b1 447
d2381af5 448 TH2I *h2 = 0x0;
449 TH1D *h = 0x0;
d85cd79c 450 TGraphErrors *gm = 0x0, *gs = 0x0;
d2381af5 451 TF1 f("f1", "gaus", -.5, .5);
874acced 452 // define iterator over graphs
453 Int_t jgraph = (Int_t)kGraphStart;
454
455 //PROCESS RESIDUAL DISTRIBUTIONS
456
457 // Clusters residuals
3d86166d 458 h2 = (TH2I *)(fContainer->At(kClusterYResidual));
765bd0ab 459 jgraph++; //skip the frame histo
d85cd79c 460 gm = (TGraphErrors*)fContainer->At(jgraph++);
461 gs = (TGraphErrors*)fContainer->At(jgraph++);
874acced 462 for(Int_t ibin = 1; ibin <= h2->GetNbinsX(); ibin++){
463 Double_t phi = h2->GetXaxis()->GetBinCenter(ibin);
464 Double_t dphi = h2->GetXaxis()->GetBinWidth(ibin)/2;
465 h = h2->ProjectionY("py", ibin, ibin);
466 h->Fit(&f, "QN", "", -0.5, 0.5);
765bd0ab 467 gm->SetPoint(ibin - 1, phi, 10.*f.GetParameter(1));
468 gm->SetPointError(ibin - 1, dphi, 10.*f.GetParError(1));
469 gs->SetPoint(ibin - 1, phi, 10.*f.GetParameter(2));
470 gs->SetPointError(ibin - 1, dphi, 10.*f.GetParError(2));
874acced 471 }
765bd0ab 472 fNRefFigures++;
d2381af5 473
874acced 474
475 //PROCESS RESOLUTION DISTRIBUTIONS
476 if(HasMCdata()){
477 // tracklet y resolution
3d86166d 478 h2 = (TH2I*)fContainer->At(kTrackletYResolution);
765bd0ab 479 jgraph++; //skip the frame histo
d85cd79c 480 gm = (TGraphErrors*)fContainer->At(jgraph++);
481 gs = (TGraphErrors*)fContainer->At(jgraph++);
d2381af5 482 for(Int_t iphi=1; iphi<=h2->GetNbinsX(); iphi++){
483 Double_t phi = h2->GetXaxis()->GetBinCenter(iphi);
484 f.SetParameter(1, 0.);f.SetParameter(2, 2.e-2);
485 h = h2->ProjectionY("py", iphi, iphi);
486 h->Fit(&f, "QN", "", -.5, .5);
487 Int_t jphi = iphi -1;
765bd0ab 488 gm->SetPoint(jphi, phi, 10.*f.GetParameter(1));
489 gm->SetPointError(jphi, 0., 10.*f.GetParError(1));
490 gs->SetPoint(jphi, phi, 10.*f.GetParameter(2));
491 gs->SetPointError(jphi, 0., 10.*f.GetParError(2));
d2381af5 492 }
765bd0ab 493 fNRefFigures++;
d2381af5 494
874acced 495 // tracklet phi resolution
3d86166d 496 h2 = (TH2I*)fContainer->At(kTrackletAngleResolution);
765bd0ab 497 jgraph++; //skip the frame histo
d85cd79c 498 gm = (TGraphErrors*)fContainer->At(jgraph++);
499 gs = (TGraphErrors*)fContainer->At(jgraph++);
d2381af5 500 for(Int_t iphi=1; iphi<=h2->GetNbinsX(); iphi++){
501 Double_t phi = h2->GetXaxis()->GetBinCenter(iphi);
502 f.SetParameter(1, 0.);f.SetParameter(2, 2.e-2);
503 h = h2->ProjectionY("py", iphi, iphi);
504 h->Fit(&f, "QN", "", -.5, .5);
505 Int_t jphi = iphi -1;
506 gm->SetPoint(jphi, phi, f.GetParameter(1));
507 gm->SetPointError(jphi, 0., f.GetParError(1));
508 gs->SetPoint(jphi, phi, f.GetParameter(2));
509 gs->SetPointError(jphi, 0., f.GetParError(2));
510 }
765bd0ab 511 fNRefFigures++;
d2381af5 512 }
765bd0ab 513
d85cd79c 514 return kTRUE;
515}
516
517
518//________________________________________________________
519void AliTRDtrackingResolution::Terminate(Option_t *)
520{
521 if(fDebugStream){
522 delete fDebugStream;
523 fDebugStream = 0x0;
524 fDebugLevel = 0;
525 }
526 if(HasPostProcess()) PostProcess();
874acced 527}
d2381af5 528
874acced 529//________________________________________________________
3d86166d 530TObjArray* AliTRDtrackingResolution::Histos()
874acced 531{
765bd0ab 532 if(!fContainer) fContainer = new TObjArray(25);
3d86166d 533 return fContainer;
77203477 534}
535
aaf47b30 536
537//________________________________________________________
538void AliTRDtrackingResolution::SetRecoParam(AliTRDrecoParam *r)
539{
3d86166d 540
aaf47b30 541 fReconstructor->SetRecoParam(r);
542}