Adjust needed statistics for Vd calibration and normalization with TPC signal for...
[u/mrichter/AliRoot.git] / TRD / qaGui / AliTRDqaGuiClusters.cxx
1 /**************************************************************************
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-commercial 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  **************************************************************************/
15
16 /* $Id: AliTRDqaGuiClusters.cxx 23871 2008-02-12 11:48:20Z hristov $ */
17
18 //////////////////////////////////////////////////////////////////////////////////
19 //
20 // This class is a Graphical User Interface for the Quality Monitorig 
21 // of clusters on the full detector level.
22 // It displays histograms created by the AliTRDQADataMakerRec 
23 // run during the reconstruction 
24 //
25 // S. Radomski 
26 // Uni-Heidelberg
27 // Feb. 2008
28 // 
29 //////////////////////////////////////////////////////////////////////////////////
30
31 #include "AliTRDqaGuiClusters.h"
32
33 #include "TH1D.h"
34 #include "TH2D.h"
35 #include "TFile.h"
36 #include "TCanvas.h"
37 #include "TRootEmbeddedCanvas.h"
38
39 ClassImp(AliTRDqaGuiClusters)
40
41 const Int_t AliTRDqaGuiClusters::fgkLogList[4] = {0,1,0,1};
42
43 //////////////////////////////////////////////////////////////////////////////////
44 AliTRDqaGuiClusters::AliTRDqaGuiClusters() 
45   : TGCompositeFrame() {
46   //
47   // Default constructor
48   //
49
50   for (Int_t i = 0; i < 4; i++) {
51     fNameList[i]   = 0x0;
52     fHistList[i]   = 0x0;
53     fCanvasList[i] = 0x0;
54     for (Int_t j = 0; j < 3; j++) {
55       fHistRefs[i][j] = 0x0;
56     }
57   }
58
59 }
60
61 //////////////////////////////////////////////////////////////////////////////////
62 AliTRDqaGuiClusters::AliTRDqaGuiClusters(TGWindow *parent) 
63   : TGCompositeFrame(parent, 720, 400) {
64   //
65   // main constructor
66   //
67   
68   SetLayoutManager(new TGMatrixLayout(this,2,2,1,1));
69
70   fNameList[0] = "detMap";
71   fNameList[1] = "nCls";
72   fNameList[2] = "signal";
73   fNameList[3] = "ampMPV";
74
75   for(Int_t i=0; i<4; i++) {
76     fCanvasList[i] = new TRootEmbeddedCanvas(fNameList[i], this, 480, 320);
77     AddFrame(fCanvasList[i]);
78   }
79   
80   for(Int_t i=0; i<4; i++) {
81     fHistList[i] = 0;
82   }
83
84 }
85
86 //////////////////////////////////////////////////////////////////////////////////
87 void AliTRDqaGuiClusters::SetQAFile(const char *filename) {
88   //
89   // sets the file with histograms
90   //
91
92   for(Int_t i=0; i<4; i++) {
93     if (fHistList[i]) delete fHistList[i];
94     for(Int_t j=0; j<3; j++) 
95       if (fHistRefs[i][j]) delete fHistRefs[i][j];
96   }
97   
98   const char *opt[4] = {"colz", "", "", ""};
99
100   TFile *file = new TFile(filename);
101   file->cd("TRD/RecPoints");
102   
103   for(int i=0; i<4; i++) {
104
105     fCanvasList[i]->GetCanvas()->cd();
106     gPad->SetLogy(fgkLogList[i]);
107
108     fHistList[i] = (TH1D*)gDirectory->Get(Form("qaTRD_recPoints_%s", fNameList[i]));
109     if (fHistList[i]) fHistList[i]->Draw(opt[i]);
110     // if (fgkLogList[i]) fHistList[i]->SetMinimum(0.1);
111     
112     TH1D *refHist = (TH1D*)gDirectory->Get(Form("qaTRD_recPoints_%s_%s", fNameList[i], "ref"));
113     
114     if (refHist) {
115       BuildColor(i, refHist);
116       for(Int_t j=0; j<3; j++) fHistRefs[i][j]->Draw("SAME");
117       delete refHist;
118     }
119     
120     fCanvasList[i]->GetCanvas()->Update();
121   }
122 }
123
124 //////////////////////////////////////////////////////////////////////////////////
125 /*
126 TH2D *AliTRDqaGuiClusters::BuildHisto(TH1D *ref, TH1D *data) {
127   
128   Int_t nbinsx = ref->GetNbinsX();
129   Double_t minx = ref->GetXaxis()->GetXmin();
130   Double_t maxx = ref->GetXaxis()->GetXmax();
131   Double_t min  = data->GetMinimum(); 
132   Double_t max  = data->GetMaximum();
133   
134   TH2D *pad = new TH2D("pad", "", nbinsx, minx, maxx, 1, 0.7, 0.9);
135
136   // rewriting
137   for(Int_t i=0; i<nbinsx; i++) {
138     Double_t x = ref->GetBinCenter(i+1);
139     pad->Fill(x, 0.8, ref->GetBinContent(i+1));
140   }
141   
142   pad->SetMinimum(-1);
143   return pad;
144 }
145 */
146
147 //////////////////////////////////////////////////////////////////////////////////
148
149 void AliTRDqaGuiClusters::BuildColor(Int_t i, TH1D *ref) {
150   
151   const Int_t nHist = 3;
152   const Int_t clr[nHist] = {3, 5, 2};
153
154   for(Int_t j=0; j<nHist; j++) {
155     fHistRefs[i][j] = (TH1D*)fHistList[i]->Clone(Form("%s_%d", fHistList[i]->GetName(), j));
156     fHistRefs[i][j]->SetFillColor(clr[j]);
157     fHistRefs[i][j]->SetLineColor(clr[j]);   
158   }
159   
160   for(Int_t k=0; k<ref->GetNbinsX(); k++) {
161     Double_t v = ref->GetBinContent(k+1);
162     if (v < 0.3) fHistRefs[i][1]->SetBinContent(k+1, 0);
163     if (v < 0.7) fHistRefs[i][2]->SetBinContent(k+1, 0);
164   }
165   
166 }
167
168 //////////////////////////////////////////////////////////////////////////////////