trend plots for DET and ESD performance tasks
[u/mrichter/AliRoot.git] / TRD / qaRec / AliTRDcheckDET.cxx
1 #include <TAxis.h>
2 #include <TCanvas.h>
3 #include <TFile.h>
4 #include <TH1F.h>
5 #include <TH1I.h>
6 #include <TF1.h>
7 #include <TGaxis.h>
8 #include <TGraph.h>
9 #include <TLegend.h>
10 #include <TMath.h>
11 #include <TMap.h>
12 #include <TObjArray.h>
13 #include <TObject.h>
14 #include <TObjString.h>
15
16 #include <TPad.h>
17 #include <TProfile.h>
18 #include <TProfile2D.h>
19 #include <TROOT.h>
20
21 #include "AliLog.h"
22 #include "AliTRDcluster.h"
23 #include "AliESDHeader.h"
24 #include "AliESDRun.h"
25 #include "AliESDtrack.h"
26 #include "AliTRDgeometry.h"
27 #include "AliTRDpadPlane.h"
28 #include "AliTRDSimParam.h"
29 #include "AliTRDseedV1.h"
30 #include "AliTRDtrackV1.h"
31 #include "AliTRDtrackerV1.h"
32 #include "AliTRDReconstructor.h"
33 #include "AliTrackReference.h"
34 #include "AliTrackPointArray.h"
35 #include "AliTracker.h"
36 #include "TTreeStream.h"
37
38 #include "info/AliTRDtrackInfo.h"
39 #include "info/AliTRDeventInfo.h"
40 #include "AliTRDcheckDET.h"
41
42 #include <cstdio>
43 #include <iostream>
44
45 ////////////////////////////////////////////////////////////////////////////
46 //                                                                        //
47 //  Reconstruction QA                                                     //
48 //                                                                        //
49 //  Task doing basic checks for tracking and detector performance         //
50 //                                                                        //
51 //  Authors:                                                              //
52 //    Anton Andronic <A.Andronic@gsi.de>                                  //
53 //    Alexandru Bercuci <A.Bercuci@gsi.de>                                //
54 //    Markus Fasel <M.Fasel@gsi.de>                                       //
55 //                                                                        //
56 ////////////////////////////////////////////////////////////////////////////
57
58 //_______________________________________________________
59 AliTRDcheckDET::AliTRDcheckDET():
60   AliTRDrecoTask("checkDET", "Basic TRD data checker")
61   ,fEventInfo(0x0)
62   ,fTriggerNames(0x0)
63   ,fReconstructor(0x0)
64   ,fGeo(0x0)
65 {
66   //
67   // Default constructor
68   //
69   DefineInput(1,AliTRDeventInfo::Class());
70   fReconstructor = new AliTRDReconstructor;
71   fReconstructor->SetRecoParam(AliTRDrecoParam::GetLowFluxParam());
72   fGeo = new AliTRDgeometry;
73   InitFunctorList();
74 }
75
76 //_______________________________________________________
77 AliTRDcheckDET::~AliTRDcheckDET(){
78   //
79   // Destructor
80   // 
81   if(fTriggerNames) delete fTriggerNames;
82   delete fReconstructor;
83   delete fGeo;
84 }
85
86 //_______________________________________________________
87 void AliTRDcheckDET::ConnectInputData(Option_t *opt){
88   //
89   // Connect the Input data with the task
90   //
91   AliTRDrecoTask::ConnectInputData(opt);
92   fEventInfo = dynamic_cast<AliTRDeventInfo *>(GetInputData(1));
93 }
94
95 //_______________________________________________________
96 void AliTRDcheckDET::CreateOutputObjects(){
97   //
98   // Create Output Objects
99   //
100   OpenFile(0,"RECREATE");
101   fContainer = Histos();
102   if(!fTriggerNames) fTriggerNames = new TMap();
103 }
104
105 //_______________________________________________________
106 void AliTRDcheckDET::Exec(Option_t *opt){
107   //
108   // Execution function
109   // Filling TRD quality histos
110   //
111   if(!HasMCdata() && fEventInfo->GetEventHeader()->GetEventType() != 7) return; // For real data we select only physical events
112   AliTRDrecoTask::Exec(opt);  
113   Int_t nTracks = 0;            // Count the number of tracks per event
114   Int_t triggermask = fEventInfo->GetEventHeader()->GetTriggerMask();
115   TString triggername =  fEventInfo->GetRunInfo()->GetFiredTriggerClasses(triggermask);
116   if(fDebugLevel > 6)printf("Trigger cluster: %d, Trigger class: %s\n", triggermask, triggername.Data());
117   dynamic_cast<TH1F *>(fContainer->UncheckedAt(kNeventsTrigger))->Fill(triggermask);
118   for(Int_t iti = 0; iti < fTracks->GetEntriesFast(); iti++){
119     if(!fTracks->UncheckedAt(iti)) continue;
120     AliTRDtrackInfo *fTrackInfo = dynamic_cast<AliTRDtrackInfo *>(fTracks->UncheckedAt(iti));
121     if(!fTrackInfo->GetTrack()) continue;
122     nTracks++;
123   }
124   if(nTracks){
125     dynamic_cast<TH1F *>(fContainer->UncheckedAt(kNeventsTriggerTracks))->Fill(triggermask);
126     dynamic_cast<TH1F *>(fContainer->UncheckedAt(kNtracksEvent))->Fill(nTracks);
127   }
128   if(triggermask <= 20 && !fTriggerNames->FindObject(Form("%d", triggermask))){
129     fTriggerNames->Add(new TObjString(Form("%d", triggermask)), new TObjString(triggername));
130     // also set the label for both histograms
131     TH1 *histo = dynamic_cast<TH1F *>(fContainer->UncheckedAt(kNeventsTriggerTracks));
132     histo->GetXaxis()->SetBinLabel(histo->FindBin(triggermask), triggername);
133     histo = dynamic_cast<TH1F *>(fContainer->UncheckedAt(kNeventsTrigger));
134     histo->GetXaxis()->SetBinLabel(histo->FindBin(triggermask), triggername);
135   }
136   PostData(0, fContainer);
137 }
138
139 //_______________________________________________________
140 void AliTRDcheckDET::Terminate(Option_t *){
141   //
142   // Terminate function
143   //
144 }
145
146 //_______________________________________________________
147 Bool_t AliTRDcheckDET::PostProcess(){
148   //
149   // Do Postprocessing (for the moment set the number of Reference histograms)
150   //
151   
152   TH1 * h = 0x0;
153   
154   // Calculate of the trigger clusters purity
155   h = dynamic_cast<TH1F *>(fContainer->FindObject("hEventsTrigger"));
156   TH1F *h1 = dynamic_cast<TH1F *>(fContainer->FindObject("hEventsTriggerTracks"));
157   h1->Divide(h);
158   Float_t purities[20], val = 0;
159   TString triggernames[20];
160   Int_t nTriggerClasses = 0;
161   for(Int_t ibin = 1; ibin <= h->GetNbinsX(); ibin++){
162     if((val = h1->GetBinContent(ibin))){
163       purities[nTriggerClasses] = val;
164       triggernames[nTriggerClasses] = h1->GetXaxis()->GetBinLabel(ibin);
165       nTriggerClasses++;
166     }
167   }
168   h = dynamic_cast<TH1F *>(fContainer->UncheckedAt(kTriggerPurity));
169   TAxis *ax = h->GetXaxis();
170   for(Int_t itrg = 0; itrg < nTriggerClasses; itrg++){
171     h->Fill(itrg, purities[itrg]);
172     ax->SetBinLabel(itrg+1, triggernames[itrg].Data());
173   }
174   ax->SetRangeUser(-0.5, nTriggerClasses+.5);
175   h->GetYaxis()->SetRangeUser(0,1);
176
177   fNRefFigures = 14;
178
179   return kTRUE;
180 }
181
182 //_______________________________________________________
183 Bool_t AliTRDcheckDET::GetRefFigure(Int_t ifig){
184   //
185   // Setting Reference Figures
186   //
187   gPad->SetLogy(0);
188   gPad->SetLogx(0);
189   TH1 *h = 0x0;
190   switch(ifig){
191   case kNclustersTrack:
192     (h=(TH1F*)fContainer->FindObject("hNcls"))->Draw("pl");
193     PutTrendValue("NClustersTrack", h->GetMean(), h->GetRMS());
194     return kTRUE;
195   case kNclustersTracklet:
196     (h =(TH1F*)fContainer->FindObject("hNclTls"))->Draw("pc");
197     PutTrendValue("NClustersTracklet", h->GetMean(), h->GetRMS());
198     return kTRUE;
199   case kNtrackletsTrack:
200     h=MakePlotNTracklets();
201     PutTrendValue("NTrackletsTracklet", h->GetMean(), h->GetRMS());
202     return kTRUE;
203   case kNtrackletsCross:
204     h = (TH1F*)fContainer->FindObject("hNtlsCross");
205     if(!MakeBarPlot(h, kRed)) break;
206     PutTrendValue("NTrackletsCross", h->GetMean(), h->GetRMS());
207     return kTRUE;
208   case kNtrackletsFindable:
209     h = (TH1F*)fContainer->FindObject("hNtlsFindable");
210     if(!MakeBarPlot(h, kGreen)) break;
211     PutTrendValue("NTrackletsFindable", h->GetMean(), h->GetRMS());
212     return kTRUE;
213   case kNtracksEvent:
214     (h = (TH1F*)fContainer->FindObject("hNtrks"))->Draw("pl");
215     PutTrendValue("NTracksEvent", h->GetMean(), h->GetRMS());
216     return kTRUE;
217   case kNtracksSector:
218     h = (TH1F*)fContainer->FindObject("hNtrksSector");
219     if(!MakeBarPlot(h, kGreen)) break;
220     PutTrendValue("NTracksSector", h->GetMean(2), h->GetRMS(2));
221     return kTRUE;
222   case kTrackStatus:
223     ((TH1I *)fContainer->FindObject("hTrackStatus"))->Draw("c");
224     return kTRUE;
225   case kTrackletStatus:
226     ((TH2S *)fContainer->FindObject("hTrackletStatus"))->Draw("colz");
227     return kTRUE;
228   case kChi2:
229     MakePlotChi2();
230     return kTRUE;
231   case kPH:
232     MakePlotPulseHeight();
233     return kTRUE;
234   case kChargeCluster:
235     (h = (TH1F*)fContainer->FindObject("hQcl"))->Draw("c");
236     gPad->SetLogy(1);
237     PutTrendValue("ChargeCluster", h->GetMaximumBin(), h->GetRMS());
238     return kTRUE;
239   case kChargeTracklet:
240     (h=(TH1F*)fContainer->FindObject("hQtrklt"))->Draw("c");
241     PutTrendValue("ChargeTracklet", h->GetMaximumBin(), h->GetRMS());
242     return kTRUE;
243   case kNeventsTrigger:
244     ((TH1F*)fContainer->FindObject("hEventsTrigger"))->Draw("");
245     return kTRUE;
246   case kNeventsTriggerTracks:
247     ((TH1F*)fContainer->FindObject("hEventsTriggerTracks"))->Draw("");
248     return kTRUE;
249   case kTriggerPurity: 
250     if(!MakeBarPlot((TH1F*)fContainer->FindObject("hTriggerPurity"), kGreen)) break;
251     break;
252   default:
253     break;
254   }
255   AliInfo(Form("Reference plot [%d] missing result", ifig));
256   return kFALSE;
257 }
258
259 //_______________________________________________________
260 TObjArray *AliTRDcheckDET::Histos(){
261   //
262   // Create QA histograms
263   //
264   if(fContainer) return fContainer;
265   
266   fContainer = new TObjArray(20);
267   //fContainer->SetOwner(kTRUE);
268
269   // Register Histograms
270   TH1 * h = NULL;
271   TAxis *axis = NULL;
272   if(!(h = (TH1F *)gROOT->FindObject("hNcls"))){
273     h = new TH1F("hNcls", "N_{clusters} / track", 181, -0.5, 180.5);
274     h->GetXaxis()->SetTitle("N_{clusters}");
275     h->GetYaxis()->SetTitle("Entries");
276   } else h->Reset();
277   fContainer->AddAt(h, kNclustersTrack);
278
279   if(!(h = (TH1F *)gROOT->FindObject("hNclTls"))){
280     h = new TH1F("hNclTls","N_{clusters} / tracklet", 51, -0.5, 50.5);
281     h->GetXaxis()->SetTitle("N_{clusters}");
282     h->GetYaxis()->SetTitle("Entries");
283   } else h->Reset();
284   fContainer->AddAt(h, kNclustersTracklet);
285
286   if(!(h = (TH1F *)gROOT->FindObject("hNtls"))){
287     h = new TH1F("hNtls", "N_{tracklets} / track", AliTRDgeometry::kNlayer, 0.5, 6.5);
288     h->GetXaxis()->SetTitle("N^{tracklet}");
289     h->GetYaxis()->SetTitle("freq. [%]");
290   } else h->Reset();
291   fContainer->AddAt(h, kNtrackletsTrack);
292
293   if(!(h = (TH1F *)gROOT->FindObject("htlsSTA"))){
294     h = new TH1F("hNtlsSTA", "N_{tracklets} / track (Stand Alone)", AliTRDgeometry::kNlayer, 0.5, 6.5);
295     h->GetXaxis()->SetTitle("N^{tracklet}");
296     h->GetYaxis()->SetTitle("freq. [%]");
297   }
298   fContainer->AddAt(h, kNtrackletsSTA);
299
300   if(!(h = (TH1F *)gROOT->FindObject("htlsBAR"))){
301     h = new TH1F("hNtlsBAR", "N_{tracklets} / track (Barrel)", AliTRDgeometry::kNlayer, 0.5, 6.5);
302     h->GetXaxis()->SetTitle("N^{tracklet}");
303     h->GetYaxis()->SetTitle("freq. [%]");
304   }
305   fContainer->AddAt(h, kNtrackletsBAR);
306
307   // 
308   if(!(h = (TH1F *)gROOT->FindObject("hNtlsCross"))){
309     h = new TH1F("hNtlsCross", "N_{tracklets}^{cross} / track", 7, -0.5, 6.5);
310     h->GetXaxis()->SetTitle("n_{row cross}");
311     h->GetYaxis()->SetTitle("freq. [%]");
312   } else h->Reset();
313   fContainer->AddAt(h, kNtrackletsCross);
314
315   if(!(h = (TH1F *)gROOT->FindObject("hNtlsFindable"))){
316     h = new TH1F("hNtlsFindable", "Found/Findable Tracklets" , 101, -0.005, 1.005);
317     h->GetXaxis()->SetTitle("r [a.u]");
318     h->GetYaxis()->SetTitle("Entries");
319   } else h->Reset();
320   fContainer->AddAt(h, kNtrackletsFindable);
321
322   if(!(h = (TH1F *)gROOT->FindObject("hNtrks"))){
323     h = new TH1F("hNtrks", "N_{tracks} / event", 100, 0, 100);
324     h->GetXaxis()->SetTitle("N_{tracks}");
325     h->GetYaxis()->SetTitle("Entries");
326   } else h->Reset();
327   fContainer->AddAt(h, kNtracksEvent);
328
329   if(!(h = (TH1F *)gROOT->FindObject("hNtrksSector"))){
330     h = new TH1F("hNtrksSector", "N_{tracks} / sector", AliTRDgeometry::kNsector, -0.5, 17.5);
331     h->GetXaxis()->SetTitle("sector");
332     h->GetYaxis()->SetTitle("freq. [%]");
333   } else h->Reset();
334   fContainer->AddAt(h, kNtracksSector);
335
336   if(!(h = (TH1I *)gROOT->FindObject("hTrackStatus"))){
337     h = new TH1I("hTrackStatus", "Track Status", 7,0,7);
338     axis = h->GetXaxis();
339     axis->SetBinLabel(axis->GetFirst() + 0, "OK");
340     axis->SetBinLabel(axis->GetFirst() + 1, "PROL");
341     axis->SetBinLabel(axis->GetFirst() + 2, "PROP");
342     axis->SetBinLabel(axis->GetFirst() + 3, "ADJ");
343     axis->SetBinLabel(axis->GetFirst() + 4, "SNP");
344     axis->SetBinLabel(axis->GetFirst() + 5, "TLIN");
345     axis->SetBinLabel(axis->GetFirst() + 6, "UP");
346   }
347   fContainer->AddAt(h, kTrackStatus);
348
349   if(!(h = (TH2S *)gROOT->FindObject("hTrackletStatus"))){
350     h = new TH2S("hTrackletStatus", "Tracklet status", 6,0,6,8,0,8);
351     axis = h->GetYaxis();
352     axis->SetBinLabel(axis->GetFirst() + 0, "OK");
353     axis->SetBinLabel(axis->GetFirst() + 1, "Geom");
354     axis->SetBinLabel(axis->GetFirst() + 2, "Boun");
355     axis->SetBinLabel(axis->GetFirst() + 3, "NoCl");
356     axis->SetBinLabel(axis->GetFirst() + 4, "NoAttach");
357     axis->SetBinLabel(axis->GetFirst() + 5, "NoClTr");
358     axis->SetBinLabel(axis->GetFirst() + 6, "NoFit");
359     axis->SetBinLabel(axis->GetFirst() + 7, "Chi2");
360   }
361   fContainer->AddAt(h, kTrackletStatus);
362
363   // <PH> histos
364   TObjArray *arr = new TObjArray(2);
365   arr->SetOwner(kTRUE);  arr->SetName("<PH>");
366   fContainer->AddAt(arr, kPH);
367   if(!(h = (TH1F *)gROOT->FindObject("hPHt"))){
368     h = new TProfile("hPHt", "<PH>", 31, -0.5, 30.5);
369     h->GetXaxis()->SetTitle("Time / 100ns");
370     h->GetYaxis()->SetTitle("<PH> [a.u]");
371   } else h->Reset();
372   arr->AddAt(h, 0);
373   if(!(h = (TH1F *)gROOT->FindObject("hPHx")))
374     h = new TProfile("hPHx", "<PH>", 31, -0.08, 4.88);
375   else h->Reset();
376   arr->AddAt(h, 1);
377
378   // Chi2 histos
379   if(!(h = (TH2S*)gROOT->FindObject("hChi2"))){
380     h = new TH2S("hChi2", "#chi^{2} per track", AliTRDgeometry::kNlayer, .5, AliTRDgeometry::kNlayer+.5, 100, 0, 50);
381     h->SetXTitle("ndf");
382     h->SetYTitle("#chi^{2}/ndf");
383     h->SetZTitle("entries");
384   } else h->Reset();
385   fContainer->AddAt(h, kChi2);
386
387   if(!(h = (TH1F *)gROOT->FindObject("hQcl"))){
388     h = new TH1F("hQcl", "Q_{cluster}", 200, 0, 1200);
389     h->GetXaxis()->SetTitle("Q_{cluster} [a.u.]");
390     h->GetYaxis()->SetTitle("Entries");
391   }else h->Reset();
392   fContainer->AddAt(h, kChargeCluster);
393
394   if(!(h = (TH1F *)gROOT->FindObject("hQtrklt"))){
395     h = new TH1F("hQtrklt", "Q_{tracklet}", 6000, 0, 6000);
396     h->GetXaxis()->SetTitle("Q_{tracklet} [a.u.]");
397     h->GetYaxis()->SetTitle("Entries");
398   }else h->Reset();
399   fContainer->AddAt(h, kChargeTracklet);
400
401
402   if(!(h = (TH1F *)gROOT->FindObject("hEventsTrigger")))
403     h = new TH1F("hEventsTrigger", "Trigger Class", 100, 0, 100);
404   else h->Reset();
405   fContainer->AddAt(h, kNeventsTrigger);
406
407   if(!(h = (TH1F *)gROOT->FindObject("hEventsTriggerTracks")))
408     h = new TH1F("hEventsTriggerTracks", "Trigger Class (Tracks)", 100, 0, 100);
409   else h->Reset();
410   fContainer->AddAt(h, kNeventsTriggerTracks);
411
412   if(!(h = (TH1F *)gROOT->FindObject("hTriggerPurity"))){
413     h = new TH1F("hTriggerPurity", "Trigger Purity", 10, -0.5, 9.5);
414     h->GetXaxis()->SetTitle("Trigger Cluster");
415     h->GetYaxis()->SetTitle("freq.");
416   } else h->Reset();
417   fContainer->AddAt(h, kTriggerPurity);
418
419   return fContainer;
420 }
421
422 /*
423 * Plotting Functions
424 */
425
426 //_______________________________________________________
427 TH1 *AliTRDcheckDET::PlotTrackStatus(const AliTRDtrackV1 *track){
428   //
429   // Plot the track status
430   //
431   if(track) fTrack = track;
432   if(!fTrack){
433     AliWarning("No Track defined.");
434     return 0x0;
435   }
436   TH1 *h = 0x0;
437   if(!(h = dynamic_cast<TH1I *>(fContainer->At(kTrackStatus)))){
438     AliWarning("No Histogram defined.");
439     return 0x0;
440   }
441   h->Fill(fTrack->GetStatusTRD());
442   return h;
443 }
444
445 //_______________________________________________________
446 TH1 *AliTRDcheckDET::PlotTrackletStatus(const AliTRDtrackV1 *track){
447   //
448   // Plot the track status
449   //
450   if(track) fTrack = track;
451   if(!fTrack){
452     AliWarning("No Track defined.");
453     return 0x0;
454   }
455   TH1 *h = 0x0;
456   if(!(h = dynamic_cast<TH2S *>(fContainer->At(kTrackletStatus)))){
457     AliWarning("No Histogram defined.");
458     return 0x0;
459   }
460   UChar_t status = 0;
461   for(Int_t il = 0; il < AliTRDgeometry::kNlayer; il++){
462     status = fTrack->GetStatusTRD(il);
463     h->Fill(il, status);
464   }
465   return h;
466 }
467
468 //_______________________________________________________
469 TH1 *AliTRDcheckDET::PlotNClustersTracklet(const AliTRDtrackV1 *track){
470   //
471   // Plot the mean number of clusters per tracklet
472   //
473   if(track) fTrack = track;
474   if(!fTrack){
475     AliWarning("No Track defined.");
476     return 0x0;
477   }
478   TH1 *h = 0x0;
479   if(!(h = dynamic_cast<TH1F *>(fContainer->At(kNclustersTracklet)))){
480     AliWarning("No Histogram defined.");
481     return 0x0;
482   }
483   AliTRDseedV1 *tracklet = 0x0;
484   for(Int_t itl = 0; itl < AliTRDgeometry::kNlayer; itl++){
485     if(!(tracklet = fTrack->GetTracklet(itl)) || !tracklet->IsOK()) continue;
486     h->Fill(tracklet->GetN2());
487   }
488   return h;
489 }
490
491 //_______________________________________________________
492 TH1 *AliTRDcheckDET::PlotNClustersTrack(const AliTRDtrackV1 *track){
493   //
494   // Plot the number of clusters in one track
495   //
496   if(track) fTrack = track;
497   if(!fTrack){
498     AliWarning("No Track defined.");
499     return 0x0;
500   }
501   TH1 *h = 0x0;
502   if(!(h = dynamic_cast<TH1F *>(fContainer->At(kNclustersTrack)))){
503     AliWarning("No Histogram defined.");
504     return 0x0;
505   }
506   
507   Int_t nclusters = 0;
508   AliTRDseedV1 *tracklet = 0x0;
509   for(Int_t itl = 0; itl < AliTRDgeometry::kNlayer; itl++){
510     if(!(tracklet = fTrack->GetTracklet(itl)) || !tracklet->IsOK()) continue;
511     nclusters += tracklet->GetN();
512     if(fDebugLevel > 2){
513       Int_t crossing = Int_t(tracklet->IsRowCross());
514       Int_t detector = tracklet->GetDetector();
515       Float_t theta = TMath::ATan(tracklet->GetZref(1));
516       Float_t phi = TMath::ATan(tracklet->GetYref(1));
517       Float_t momentum = 0.;
518       Int_t pdg = 0;
519       Int_t kinkIndex = fESD ? fESD->GetKinkIndex() : 0;
520       UShort_t TPCncls = fESD ? fESD->GetTPCncls() : 0;
521       if(fMC){
522         if(fMC->GetTrackRef()) momentum = fMC->GetTrackRef()->P();
523         pdg = fMC->GetPDG();
524       }
525       (*fDebugStream) << "NClustersTrack"
526         << "Detector="  << detector
527         << "crossing="  << crossing
528         << "momentum="  << momentum
529         << "pdg="                               << pdg
530         << "theta="                     << theta
531         << "phi="                               << phi
532         << "kinkIndex=" << kinkIndex
533         << "TPCncls="           << TPCncls
534         << "nclusters=" << nclusters
535         << "\n";
536     }
537   }
538   h->Fill(nclusters);
539   return h;
540 }
541
542
543 //_______________________________________________________
544 TH1 *AliTRDcheckDET::PlotNTrackletsTrack(const AliTRDtrackV1 *track){
545   //
546   // Plot the number of tracklets
547   //
548   if(track) fTrack = track;
549   if(!fTrack){
550     AliWarning("No Track defined.");
551     return 0x0;
552   }
553   TH1 *h = 0x0, *hMethod = 0x0;
554   if(!(h = dynamic_cast<TH1F *>(fContainer->At(kNtrackletsTrack)))){
555     AliWarning("No Histogram defined.");
556     return 0x0;
557   }
558   Int_t status = fESD->GetStatus();
559 /*  printf("in/out/refit/pid: TRD[%d|%d|%d|%d]\n", status &AliESDtrack::kTRDin ? 1 : 0, status &AliESDtrack::kTRDout ? 1 : 0, status &AliESDtrack::kTRDrefit ? 1 : 0, status &AliESDtrack::kTRDpid ? 1 : 0);*/
560   if((status & AliESDtrack::kTRDin) != 0){
561     // Full BarrelTrack
562     if(!(hMethod = dynamic_cast<TH1F *>(fContainer->At(kNtrackletsBAR))))
563       AliWarning("Method: Barrel.  Histogram not processed!");
564   } else {
565     // Stand alone Track
566     if(!(hMethod = dynamic_cast<TH1F *>(fContainer->At(kNtrackletsSTA))))
567       AliWarning("Method: StandAlone.  Histogram not processed!");
568   }
569   Int_t nTracklets = fTrack->GetNumberOfTracklets();
570   h->Fill(nTracklets);
571   hMethod->Fill(nTracklets);
572   if(fDebugLevel > 3){
573     if(nTracklets == 1){
574       // If we have one Tracklet, check in which layer this happens
575       Int_t layer = -1;
576       AliTRDseedV1 *tracklet = 0x0;
577       for(Int_t il = 0; il < AliTRDgeometry::kNlayer; il++){
578         if((tracklet = fTrack->GetTracklet(il)) && tracklet->IsOK()){layer =  il; break;}
579       }
580       (*fDebugStream) << "NTrackletsTrack"
581         << "Layer=" << layer
582         << "\n";
583     }
584   }
585   return h;
586 }
587
588
589 //_______________________________________________________
590 TH1 *AliTRDcheckDET::PlotNTrackletsRowCross(const AliTRDtrackV1 *track){
591   //
592   // Plot the number of tracklets
593   //
594   if(track) fTrack = track;
595   if(!fTrack){
596     AliWarning("No Track defined.");
597     return 0x0;
598   }
599   TH1 *h = 0x0;
600   if(!(h = dynamic_cast<TH1F *>(fContainer->At(kNtrackletsCross)))){
601     AliWarning("No Histogram defined.");
602     return 0x0;
603   }
604
605   Int_t ncross = 0;
606   AliTRDseedV1 *tracklet = 0x0;
607   for(Int_t il = 0; il < AliTRDgeometry::kNlayer; il++){
608     if(!(tracklet = fTrack->GetTracklet(il)) || !tracklet->IsOK()) continue;
609
610     if(tracklet->IsRowCross()) ncross++;
611   }
612   h->Fill(ncross);
613   return h;
614 }
615
616 //_______________________________________________________
617 TH1 *AliTRDcheckDET::PlotFindableTracklets(const AliTRDtrackV1 *track){
618   //
619   // Plots the ratio of number of tracklets vs.
620   // number of findable tracklets
621   //
622   // Findable tracklets are defined as track prolongation
623   // to layer i does not hit the dead area +- epsilon
624   //
625   // In order to check whether tracklet hist active area in Layer i, 
626   // the track is refitted and the fitted position + an uncertainty 
627   // range is compared to the chamber border (also with a different
628   // uncertainty)
629   //
630   // For the track fit two cases are distinguished:
631   // If the track is a stand alone track (defined by the status bit 
632   // encoding, then the track is fitted with the tilted Rieman model
633   // Otherwise the track is fitted with the Kalman fitter in two steps:
634   // Since the track parameters are give at the outer point, we first 
635   // fit in direction inwards. Afterwards we fit again in direction outwards
636   // to extrapolate the track to layers which are not reached by the track
637   // For the Kalman model, the radial track points have to be shifted by
638   // a distance epsilon in the direction that we want to fit
639   //
640   const Float_t epsilon = 0.01;   // dead area tolerance
641   const Float_t epsilon_R = 1;    // shift in radial direction of the anode wire position (Kalman filter only)
642   const Float_t delta_y = 0.7;    // Tolerance in the track position in y-direction
643   const Float_t delta_z = 7.0;    // Tolerance in the track position in z-direction (Padlength)
644   Double_t x_anode[AliTRDgeometry::kNlayer] = {300.2, 312.8, 325.4, 338.0, 350.6, 363.2}; // Take the default X0
645  
646   if(track) fTrack = track;
647   if(!fTrack){
648     AliWarning("No Track defined.");
649     return 0x0;
650   }
651   TH1 *h = 0x0;
652   if(!(h = dynamic_cast<TH1F *>(fContainer->At(kNtrackletsFindable)))){
653     AliWarning("No Histogram defined.");
654     return 0x0;
655   }
656   Int_t nFound = 0, nFindable = 0;
657   Int_t stack = -1;
658   Double_t ymin = 0., ymax = 0., zmin = 0., zmax = 0.;
659   Double_t y = 0., z = 0.;
660   AliTRDseedV1 *tracklet = 0x0;
661   AliTRDpadPlane *pp;  
662   for(Int_t il = 0; il < AliTRDgeometry::kNlayer; il++){
663     if((tracklet = fTrack->GetTracklet(il)) && tracklet->IsOK()){
664       tracklet->SetReconstructor(fReconstructor);
665       nFound++;
666     }
667   }
668   // 2 Different cases:
669   // 1st stand alone: here we cannot propagate, but be can do a Tilted Rieman Fit
670   // 2nd barrel track: here we propagate the track to the layers
671   AliTrackPoint points[6];
672   Float_t xyz[3];
673   memset(xyz, 0, sizeof(Float_t) * 3);
674   if(((fESD->GetStatus() & AliESDtrack::kTRDout) > 0) && !((fESD->GetStatus() & AliESDtrack::kTRDin) > 0)){
675     // stand alone track
676     for(Int_t il = 0; il < AliTRDgeometry::kNlayer; il++){
677       xyz[0] = x_anode[il];
678       points[il].SetXYZ(xyz);
679     }
680     AliTRDtrackerV1::FitRiemanTilt(const_cast<AliTRDtrackV1 *>(fTrack), 0x0, kTRUE, 6, points);
681   } else {
682     // barrel track
683     //
684     // 2 Steps:
685     // -> Kalman inwards
686     // -> Kalman outwards
687     AliTRDtrackV1 copy_track(*fTrack);  // Do Kalman on a (non-constant) copy of the track
688     AliTrackPoint points_inward[6], points_outward[6];
689     for(Int_t il = AliTRDgeometry::kNlayer; il--;){
690       // In order to avoid complications in the Kalman filter if the track points have the same radial
691       // position like the tracklets, we have to shift the radial postion of the anode wire by epsilon
692       // in the direction we want to go
693       // The track points have to be in reverse order for the Kalman Filter inwards
694       xyz[0] = x_anode[AliTRDgeometry::kNlayer - il - 1] - epsilon_R;
695       points_inward[il].SetXYZ(xyz);
696       xyz[0] = x_anode[il] + epsilon_R;
697       points_outward[il].SetXYZ(xyz);
698     }
699     /*for(Int_t ipt = 0; ipt < AliTRDgeometry::kNlayer; ipt++)
700       printf("%d. X = %f\n", ipt, points[ipt].GetX());*/
701     // Kalman inwards
702     AliTRDtrackerV1::FitKalman(&copy_track, 0x0, kFALSE, 6, points_inward);
703     memcpy(points, points_inward, sizeof(AliTrackPoint) * 6); // Preliminary store the inward results in the Array points
704     // Kalman outwards
705     AliTRDtrackerV1::FitKalman(&copy_track, 0x0, kTRUE, 6, points_inward);
706     memcpy(points, points_outward, sizeof(AliTrackPoint) * AliTRDgeometry::kNlayer);
707   }
708   for(Int_t il = 0; il < AliTRDgeometry::kNlayer; il++){
709     y = points[il].GetY();
710     z = points[il].GetZ();
711     if((stack = fGeo->GetStack(z, il)) < 0) continue; // Not findable
712     pp = fGeo->GetPadPlane(il, stack);
713     ymin = pp->GetCol0() + epsilon;
714     ymax = pp->GetColEnd() - epsilon; 
715     zmin = pp->GetRowEnd() + epsilon; 
716     zmax = pp->GetRow0() - epsilon;
717     // ignore y-crossing (material)
718     if((z + delta_z > zmin && z - delta_z < zmax) && (y + delta_y > ymin && y - delta_y < ymax)) nFindable++;
719       if(fDebugLevel > 3){
720         Double_t pos_tracklet[2] = {tracklet ? tracklet->GetYfit(0) : 0, tracklet ? tracklet->GetZfit(0) : 0};
721         Int_t hasTracklet = tracklet ? 1 : 0;
722         (*fDebugStream)   << "FindableTracklets"
723           << "layer="     << il
724           << "ytracklet=" << pos_tracklet[0]
725           << "ytrack="    << y
726           << "ztracklet=" << pos_tracklet[1]
727           << "ztrack="    << z
728           << "tracklet="  << hasTracklet
729           << "\n";
730       }
731   }
732   
733   h->Fill(nFindable > 0 ? TMath::Min(nFound/static_cast<Double_t>(nFindable), 1.) : 1);
734   if(fDebugLevel > 2) AliInfo(Form("Findable[Found]: %d[%d|%f]", nFindable, nFound, nFound/static_cast<Float_t>(nFindable > 0 ? nFindable : 1)));
735   return h;
736 }
737
738
739 //_______________________________________________________
740 TH1 *AliTRDcheckDET::PlotChi2(const AliTRDtrackV1 *track){
741   //
742   // Plot the chi2 of the track
743   //
744   if(track) fTrack = track;
745   if(!fTrack){
746     AliWarning("No Track defined.");
747     return 0x0;
748   }
749   TH1 *h = 0x0;
750   if(!(h = dynamic_cast<TH2S*>(fContainer->At(kChi2)))) {
751     AliWarning("No Histogram defined.");
752     return 0x0;
753   }
754   Int_t n = fTrack->GetNumberOfTracklets();
755   if(!n) return 0x0;
756
757   h->Fill(n, fTrack->GetChi2()/n);
758   return h;
759 }
760
761
762 //_______________________________________________________
763 TH1 *AliTRDcheckDET::PlotPHt(const AliTRDtrackV1 *track){
764   //
765   // Plot the average pulse height
766   //
767   if(track) fTrack = track;
768   if(!fTrack){
769     AliWarning("No Track defined.");
770     return 0x0;
771   }
772   TProfile *h = 0x0;
773   if(!(h = dynamic_cast<TProfile *>(((TObjArray*)(fContainer->At(kPH)))->At(0)))){
774     AliWarning("No Histogram defined.");
775     return 0x0;
776   }
777   AliTRDseedV1 *tracklet = 0x0;
778   AliTRDcluster *c = 0x0;
779   for(Int_t itl = 0; itl < AliTRDgeometry::kNlayer; itl++){
780     if(!(tracklet = fTrack->GetTracklet(itl)) || !tracklet->IsOK())continue;
781     Int_t crossing = Int_t(tracklet->IsRowCross());
782     Int_t detector = tracklet->GetDetector();
783     tracklet->ResetClusterIter();
784     while((c = tracklet->NextCluster())){
785       if(!c->IsInChamber()) continue;
786       Int_t localtime        = c->GetLocalTimeBin();
787       Double_t absolute_charge = TMath::Abs(c->GetQ());
788       h->Fill(localtime, absolute_charge);
789       if(fDebugLevel > 3){
790         Double_t distance[2];
791         GetDistanceToTracklet(distance, tracklet, c);
792         Float_t theta = TMath::ATan(tracklet->GetZref(1));
793         Float_t phi = TMath::ATan(tracklet->GetYref(1));
794         Float_t momentum = 0.;
795         Int_t pdg = 0;
796         Int_t kinkIndex = fESD ? fESD->GetKinkIndex() : 0;
797         UShort_t TPCncls = fESD ? fESD->GetTPCncls() : 0;
798         if(fMC){
799           if(fMC->GetTrackRef()) momentum = fMC->GetTrackRef()->P();
800           pdg = fMC->GetPDG();
801         }
802         (*fDebugStream) << "PHt"
803           << "Detector="        << detector
804           << "crossing="        << crossing
805           << "Timebin="         << localtime
806           << "Charge="          << absolute_charge
807           << "momentum="        << momentum
808           << "pdg="                             << pdg
809           << "theta="                   << theta
810           << "phi="                             << phi
811           << "kinkIndex="       << kinkIndex
812           << "TPCncls="         << TPCncls
813           << "dy="        << distance[0]
814           << "dz="        << distance[1]
815           << "c.="        << c
816           << "\n";
817       }
818     }
819   }
820   return h;
821 }
822
823 //_______________________________________________________
824 TH1 *AliTRDcheckDET::PlotPHx(const AliTRDtrackV1 *track){
825   //
826   // Plots the average pulse height vs the distance from the anode wire
827   // (plus const anode wire offset)
828   //
829   if(track) fTrack = track;
830   if(!fTrack){
831     AliWarning("No Track defined.");
832     return 0x0;
833   }
834   TProfile *h = 0x0;
835   if(!(h = dynamic_cast<TProfile *>(((TObjArray*)(fContainer->At(kPH)))->At(1)))){
836     AliWarning("No Histogram defined.");
837     return 0x0;
838   }
839   Float_t offset = .5*AliTRDgeometry::CamHght();
840   AliTRDseedV1 *tracklet = 0x0;
841   AliTRDcluster *c = 0x0;
842   Double_t distance = 0;
843   Double_t x, y;
844   for(Int_t itl = 0; itl < AliTRDgeometry::kNlayer; itl++){
845     if(!(tracklet = fTrack->GetTracklet(itl)) || !(tracklet->IsOK())) continue;
846     tracklet->ResetClusterIter();
847     while((c = tracklet->NextCluster())){
848       if(!c->IsInChamber()) continue;
849       x = c->GetX()-AliTRDcluster::GetXcorr(c->GetLocalTimeBin());
850       y = c->GetY()-AliTRDcluster::GetYcorr(AliTRDgeometry::GetLayer(c->GetDetector()), c->GetCenter());
851
852       distance = tracklet->GetX0() - (c->GetX() + 0.3) + offset;
853       h->Fill(distance, TMath::Abs(c->GetQ()));
854     }
855   }  
856   return h;
857 }
858
859 //_______________________________________________________
860 TH1 *AliTRDcheckDET::PlotChargeCluster(const AliTRDtrackV1 *track){
861   //
862   // Plot the cluster charge
863   //
864   if(track) fTrack = track;
865   if(!fTrack){
866     AliWarning("No Track defined.");
867     return 0x0;
868   }
869   TH1 *h = 0x0;
870   if(!(h = dynamic_cast<TH1F *>(fContainer->At(kChargeCluster)))){
871     AliWarning("No Histogram defined.");
872     return 0x0;
873   }
874   AliTRDseedV1 *tracklet = 0x0;
875   AliTRDcluster *c = 0x0;
876   for(Int_t itl = 0; itl < AliTRDgeometry::kNlayer; itl++){
877     if(!(tracklet = fTrack->GetTracklet(itl)) || !tracklet->IsOK())continue;
878     for(Int_t itime = 0; itime < AliTRDtrackerV1::GetNTimeBins(); itime++){
879       if(!(c = tracklet->GetClusters(itime))) continue;
880       h->Fill(c->GetQ());
881     }
882   }
883   return h;
884 }
885
886 //_______________________________________________________
887 TH1 *AliTRDcheckDET::PlotChargeTracklet(const AliTRDtrackV1 *track){
888   //
889   // Plot the charge deposit per chamber
890   //
891   if(track) fTrack = track;
892   if(!fTrack){
893     AliWarning("No Track defined.");
894     return 0x0;
895   }
896   TH1 *h = 0x0;
897   if(!(h = dynamic_cast<TH1F *>(fContainer->At(kChargeTracklet)))){
898     AliWarning("No Histogram defined.");
899     return 0x0;
900   }
901   AliTRDseedV1 *tracklet = 0x0;
902   AliTRDcluster *c = 0x0;
903   Double_t Qtot = 0;
904   Int_t nTracklets =fTrack->GetNumberOfTracklets();
905   for(Int_t itl = 0x0; itl < AliTRDgeometry::kNlayer; itl++){
906     if(!(tracklet = fTrack->GetTracklet(itl)) || !tracklet->IsOK()) continue;
907     Qtot = 0.;
908     for(Int_t ic = AliTRDseedV1::kNclusters; ic--;){
909       if(!(c = tracklet->GetClusters(ic))) continue;
910       Qtot += TMath::Abs(c->GetQ());
911     }
912     h->Fill(Qtot);
913     if(fDebugLevel > 3){
914       Int_t crossing = (Int_t)tracklet->IsRowCross();
915       Int_t detector = tracklet->GetDetector();
916       Float_t theta = TMath::ATan(tracklet->GetZfit(1));
917       Float_t phi = TMath::ATan(tracklet->GetYfit(1));
918       Float_t momentum = 0.;
919       Int_t pdg = 0;
920       Int_t kinkIndex = fESD ? fESD->GetKinkIndex() : 0;
921       UShort_t TPCncls = fESD ? fESD->GetTPCncls() : 0;
922       if(fMC){
923               if(fMC->GetTrackRef()) momentum = fMC->GetTrackRef()->P();
924         pdg = fMC->GetPDG();
925       }
926       (*fDebugStream) << "ChargeTracklet"
927         << "Detector="  << detector
928         << "crossing="  << crossing
929         << "momentum="  << momentum
930         << "nTracklets="<< nTracklets
931         << "pdg="                               << pdg
932         << "theta="                     << theta
933         << "phi="                               << phi
934         << "kinkIndex=" << kinkIndex
935         << "TPCncls="           << TPCncls
936         << "QT="        << Qtot
937         << "\n";
938     }
939   }
940   return h;
941 }
942
943 //_______________________________________________________
944 TH1 *AliTRDcheckDET::PlotNTracksSector(const AliTRDtrackV1 *track){
945   //
946   // Plot the number of tracks per Sector
947   //
948   if(track) fTrack = track;
949   if(!fTrack){
950     AliWarning("No Track defined.");
951     return 0x0;
952   }
953   TH1 *h = 0x0;
954   if(!(h = dynamic_cast<TH1F *>(fContainer->At(kNtracksSector)))){
955     AliWarning("No Histogram defined.");
956     return 0x0;
957   }
958
959   // TODO we should compare with
960   // sector = Int_t(track->GetAlpha() / AliTRDgeometry::GetAlpha());
961
962   AliTRDseedV1 *tracklet = 0x0;
963   Int_t sector = -1;
964   for(Int_t itl = 0; itl < AliTRDgeometry::kNlayer; itl++){
965     if(!(tracklet = fTrack->GetTracklet(itl)) || !tracklet->IsOK()) continue;
966     sector = static_cast<Int_t>(tracklet->GetDetector()/AliTRDgeometry::kNdets);
967     break;
968   }
969   h->Fill(sector);
970   return h;
971 }
972
973
974 //________________________________________________________
975 void AliTRDcheckDET::SetRecoParam(AliTRDrecoParam *r)
976 {
977
978   fReconstructor->SetRecoParam(r);
979 }
980
981 //________________________________________________________
982 void AliTRDcheckDET::GetDistanceToTracklet(Double_t *dist, AliTRDseedV1 *tracklet, AliTRDcluster *c)
983 {
984   Float_t x = c->GetX();
985   dist[0] = c->GetY() - tracklet->GetYat(x);
986   dist[1] = c->GetZ() - tracklet->GetZat(x);
987 }
988
989
990 //_______________________________________________________
991 TH1* AliTRDcheckDET::MakePlotChi2()
992 {
993 // Plot chi2/track normalized to number of degree of freedom 
994 // (tracklets) and compare with the theoretical distribution.
995 // 
996 // Alex Bercuci <A.Bercuci@gsi.de>
997
998   TH2S *h2 = (TH2S*)fContainer->At(kChi2);
999   TF1 f("fChi2", "[0]*pow(x, [1]-1)*exp(-0.5*x)", 0., 50.);
1000   TLegend *leg = new TLegend(.7,.7,.95,.95);
1001   leg->SetBorderSize(1); leg->SetHeader("Tracklets per Track");
1002   TH1D *h1 = 0x0;
1003   Bool_t kFIRST = kTRUE;
1004   for(Int_t il=1; il<=h2->GetNbinsX(); il++){
1005     h1 = h2->ProjectionY(Form("pyChi2%d", il), il, il);
1006     if(h1->Integral()<50) continue;
1007     h1->Scale(1./h1->Integral());
1008     h1->SetMarkerStyle(7);h1->SetMarkerColor(il);
1009     h1->SetLineColor(il);h1->SetLineStyle(2);
1010     f.SetParameter(1, .5*il);f.SetLineColor(il);
1011     h1->Fit(&f, "QW+", kFIRST ? "pc": "pcsame");
1012     leg->AddEntry(h1, Form("%d", il), "l");
1013     if(kFIRST){
1014       h1->GetXaxis()->SetRangeUser(0., 25.);
1015     }
1016     kFIRST = kFALSE;
1017   }
1018   leg->Draw();
1019   gPad->SetLogy();
1020   return h1;
1021 }
1022
1023
1024 //________________________________________________________
1025 TH1* AliTRDcheckDET::MakePlotNTracklets(){
1026   //
1027   // Make nice bar plot of the number of tracklets in each method
1028   //
1029   TH1F *hBAR = (TH1F *)fContainer->FindObject("hNtlsBAR");
1030   TH1F *hSTA = (TH1F *)fContainer->FindObject("hNtlsSTA");
1031   TH1F *hCON = (TH1F *)fContainer->FindObject("hNtls");
1032   TLegend *leg = new TLegend(0.6, 0.75, 0.89, 0.89);
1033
1034   Float_t scale = hCON->Integral();
1035   hCON->Scale(100./scale);
1036   hCON->SetFillColor(kRed);hCON->SetLineColor(kRed);
1037   hCON->SetBarWidth(0.2);
1038   hCON->SetBarOffset(0.6);
1039   hCON->SetStats(kFALSE);
1040   hCON->GetYaxis()->SetRangeUser(0.,40.);
1041   hCON->GetYaxis()->SetTitleOffset(1.2);
1042   hCON->Draw("bar1"); leg->AddEntry(hCON, "Total", "f");
1043
1044   hBAR->Scale(100./scale);
1045   hBAR->SetFillColor(kGreen);hBAR->SetLineColor(kGreen);
1046   hBAR->SetBarWidth(0.2);
1047   hBAR->SetBarOffset(0.2);
1048   hBAR->SetTitle("");
1049   hBAR->SetStats(kFALSE);
1050   hBAR->GetYaxis()->SetRangeUser(0.,40.);
1051   hBAR->GetYaxis()->SetTitleOffset(1.2);
1052   hBAR->Draw("bar1same"); leg->AddEntry(hBAR, "Barrel", "f");
1053
1054   hSTA->Scale(100./scale);
1055   hSTA->SetFillColor(kBlue);hSTA->SetLineColor(kBlue);
1056   hSTA->SetBarWidth(0.2);
1057   hSTA->SetBarOffset(0.4);
1058   hSTA->SetTitle("");
1059   hSTA->SetStats(kFALSE);
1060   hSTA->GetYaxis()->SetRangeUser(0.,40.);
1061   hSTA->GetYaxis()->SetTitleOffset(1.2);
1062   hSTA->Draw("bar1same"); leg->AddEntry(hSTA, "Stand Alone", "f");
1063   leg->Draw();
1064   gPad->Update();
1065   return hCON;
1066 }
1067
1068 //________________________________________________________
1069 TH1* AliTRDcheckDET::MakePlotPulseHeight(){
1070   //
1071   // Create Plot of the Pluse Height Spectrum
1072   //
1073   TH1 *h, *h1, *h2;
1074   TObjArray *arr = (TObjArray*)fContainer->FindObject("<PH>");
1075   h = (TH1F*)arr->At(0);
1076   h->SetMarkerStyle(24);
1077   h->SetMarkerColor(kBlack);
1078   h->SetLineColor(kBlack);
1079   h->Draw("e1");
1080 //   copy the second histogram in a new one with the same x-dimension as the phs with respect to time
1081   h1 = (TH1F *)arr->At(1);
1082   h2 = new TH1F("hphs1","Average PH", 31, -0.5, 30.5);
1083   for(Int_t ibin = h1->GetXaxis()->GetFirst(); ibin < h1->GetNbinsX(); ibin++) 
1084     h2->SetBinContent(ibin, h1->GetBinContent(ibin));
1085   h2->SetMarkerStyle(22);
1086   h2->SetMarkerColor(kBlue);
1087   h2->SetLineColor(kBlue);
1088   h2->Draw("e1same");
1089   gPad->Update();
1090 //   create axis according to the histogram dimensions of the original second histogram
1091   TGaxis *axis = new TGaxis(gPad->GetUxmin(),
1092                     gPad->GetUymax(),
1093                     gPad->GetUxmax(),
1094                     gPad->GetUymax(),
1095                     -0.08, 4.88, 510,"-L");
1096   axis->SetLineColor(kBlue);
1097   axis->SetLabelColor(kBlue);
1098   axis->SetTextColor(kBlue);
1099   axis->SetTitle("x_{0}-x_{c} [cm]");
1100   axis->Draw();
1101   return h1;
1102 }
1103
1104 //________________________________________________________
1105 Bool_t AliTRDcheckDET::MakeBarPlot(TH1 *histo, Int_t color){
1106   //
1107   // Draw nice bar plots
1108   //
1109   if(!histo->GetEntries()) return kFALSE;
1110   histo->Scale(100./histo->Integral());
1111   histo->SetFillColor(color);
1112   histo->SetBarOffset(.2);
1113   histo->SetBarWidth(.6);
1114   histo->Draw("bar1");
1115   return kTRUE;
1116 }