TRD module
[u/mrichter/AliRoot.git] / TRD / TRDbase / AliTRDCalPadStatus.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$ */
17
18 ///////////////////////////////////////////////////////////////////////////////
19 //                                                                           //
20 //  TRD calibration class for the single pad status                          //
21 //                                                                           //
22 ///////////////////////////////////////////////////////////////////////////////
23
24 #include <TStyle.h>
25 #include <TCanvas.h>
26 #include <TH1F.h>
27 #include <TH2F.h>
28
29 #include "AliTRDCalPadStatus.h"
30 #include "AliTRDgeometry.h"
31 #include "AliTRDpadPlane.h"
32 #include "AliTRDCalSingleChamberStatus.h"  // test
33
34 ClassImp(AliTRDCalPadStatus)
35
36 //_____________________________________________________________________________
37 AliTRDCalPadStatus::AliTRDCalPadStatus()
38   :TNamed()
39 {
40   //
41   // AliTRDCalPadStatus default constructor
42   //
43
44   for (Int_t idet = 0; idet < kNdet; idet++) {
45     fROC[idet] = 0;
46   }
47
48 }
49
50 //_____________________________________________________________________________
51 AliTRDCalPadStatus::AliTRDCalPadStatus(const Text_t *name, const Text_t *title)
52   :TNamed(name,title)
53 {
54   //
55   // AliTRDCalPadStatus constructor
56   //
57
58   for (Int_t isec = 0; isec < kNsect; isec++) {
59     for (Int_t ipla = 0; ipla < kNplan; ipla++) {
60       for (Int_t icha = 0; icha < kNcham; icha++) {
61         Int_t idet = AliTRDgeometry::GetDetector(ipla,icha,isec);
62         fROC[idet] = new AliTRDCalSingleChamberStatus(ipla,icha,144);
63       }
64     }
65   }
66
67 }
68
69 //_____________________________________________________________________________
70 AliTRDCalPadStatus::AliTRDCalPadStatus(const AliTRDCalPadStatus &c)
71   :TNamed(c)
72 {
73   //
74   // AliTRDCalPadStatus copy constructor
75   //
76
77   ((AliTRDCalPadStatus &) c).Copy(*this);
78
79 }
80
81 //_____________________________________________________________________________
82 AliTRDCalPadStatus::~AliTRDCalPadStatus()
83 {
84   //
85   // AliTRDCalPadStatus destructor
86   //
87
88   for (Int_t idet = 0; idet < kNdet; idet++) {
89     if (fROC[idet]) {
90       delete fROC[idet];
91       fROC[idet] = 0;
92     }
93   }
94
95 }
96
97 //_____________________________________________________________________________
98 AliTRDCalPadStatus &AliTRDCalPadStatus::operator=(const AliTRDCalPadStatus &c)
99 {
100   //
101   // Assignment operator
102   //
103
104   if (this != &c) ((AliTRDCalPadStatus &) c).Copy(*this);
105   return *this;
106
107 }
108
109 //_____________________________________________________________________________
110 void AliTRDCalPadStatus::Copy(TObject &c) const
111 {
112   //
113   // Copy function
114   //
115
116   for (Int_t idet = 0; idet < kNdet; idet++) {
117     if (fROC[idet]) {
118       fROC[idet]->Copy(*((AliTRDCalPadStatus &) c).fROC[idet]);
119     }
120   }
121
122   TObject::Copy(c);
123
124 }
125
126 //_____________________________________________________________________________
127 Bool_t AliTRDCalPadStatus::CheckStatus(Int_t d, Int_t col, Int_t row, Int_t bitMask) const
128 {
129   //
130   // Checks the pad status
131   //
132
133   AliTRDCalSingleChamberStatus *roc = GetCalROC(d);
134   if (!roc) {
135     return kFALSE;
136   }
137   else {
138     return (roc->GetStatus(col, row) & bitMask) ? kTRUE : kFALSE;
139   }
140
141 }
142
143 //_____________________________________________________________________________
144 AliTRDCalSingleChamberStatus* AliTRDCalPadStatus::GetCalROC(Int_t p, Int_t c, Int_t s) const
145
146   //
147   // Returns the readout chamber of this pad
148   //
149
150   return fROC[AliTRDgeometry::GetDetector(p,c,s)];   
151
152 }
153
154 //_____________________________________________________________________________
155 TH1F *AliTRDCalPadStatus::MakeHisto1D()
156 {
157   //
158   // Make 1D histo
159   //
160
161   char  name[1000];
162   snprintf(name,1000,"%s Pad 1D",GetTitle());
163   TH1F * his = new TH1F(name,name,6, -0.5,5.5);
164   his->GetXaxis()->SetBinLabel(1,"Good");
165   his->GetXaxis()->SetBinLabel(2,"Masked");
166   his->GetXaxis()->SetBinLabel(3,"PadBridgedLeft");
167   his->GetXaxis()->SetBinLabel(4,"PadBridgedRight");
168   his->GetXaxis()->SetBinLabel(5,"ReadSecond");
169   his->GetXaxis()->SetBinLabel(6,"NotConnected");
170
171   for (Int_t idet = 0; idet < kNdet; idet++) 
172     {
173       if (fROC[idet])
174         {
175           for (Int_t ichannel=0; ichannel<fROC[idet]->GetNchannels(); ichannel++)
176             {
177               Int_t status = (Int_t) fROC[idet]->GetStatus(ichannel);
178               if(status==2)  status= 1;
179               if(status==4)  status= 2;
180               if(status==8)  status= 3;
181               if(status==16) status= 4;
182               if(status==32) status= 5;
183               his->Fill(status);
184             }
185         }
186     }
187
188   return his;
189
190 }
191
192 //_____________________________________________________________________________
193 TH2F *AliTRDCalPadStatus::MakeHisto2DSmPl(Int_t sm, Int_t pl)
194 {
195   //
196   // Make 2D graph
197   //
198
199   gStyle->SetPalette(1);
200   AliTRDgeometry *trdGeo = new AliTRDgeometry();
201   AliTRDpadPlane *padPlane0 = trdGeo->GetPadPlane(pl,0);
202   Double_t row0    = padPlane0->GetRow0();
203   Double_t col0    = padPlane0->GetCol0();
204
205   char  name[1000];
206   snprintf(name,1000,"%s Pad 2D sm %d pl %d",GetTitle(),sm,pl);
207   TH2F * his = new TH2F( name, name, 88,-TMath::Abs(row0),TMath::Abs(row0)
208                                    ,148,-TMath::Abs(col0),TMath::Abs(col0));
209
210   // Where we begin
211   Int_t offsetsmpl = 30*sm+pl;
212
213   for (Int_t k = 0; k < kNcham; k++){
214     Int_t det = offsetsmpl+k*6;
215     if (fROC[det]){
216       AliTRDCalSingleChamberStatus * calRoc = fROC[det];
217       for (Int_t icol=0; icol<calRoc->GetNcols(); icol++){
218         for (Int_t irow=0; irow<calRoc->GetNrows(); irow++){
219           Int_t binz     = 0;
220           Int_t kb       = kNcham-1-k;
221           Int_t krow     = calRoc->GetNrows()-1-irow;
222           Int_t kcol     = calRoc->GetNcols()-1-icol;
223           if(kb > 2) binz = 16*(kb-1)+12+krow+1+2*(kb+1);
224           else binz = 16*kb+krow+1+2*(kb+1); 
225           Int_t biny = kcol+1+2;
226           Float_t value = calRoc->GetStatus(icol,irow);
227           his->SetBinContent(binz,biny,value);
228         }
229       }
230       for(Int_t icol = 1; icol < 147; icol++){
231         for(Int_t l = 0; l < 2; l++){
232           Int_t binz     = 0;
233           Int_t kb       = kNcham-1-k;
234           if(kb > 2) binz = 16*(kb-1)+12+1+2*(kb+1)-(l+1);
235           else binz = 16*kb+1+2*(kb+1)-(l+1); 
236           his->SetBinContent(binz,icol,50.0);
237         }
238       }
239     }
240   }
241   for(Int_t icol = 1; icol < 147; icol++){
242     his->SetBinContent(88,icol,50.0);
243     his->SetBinContent(87,icol,50.0);
244   }
245   for(Int_t irow = 1; irow < 89; irow++){
246     his->SetBinContent(irow,1,50.0);
247     his->SetBinContent(irow,2,50.0);
248     his->SetBinContent(irow,147,50.0);
249     his->SetBinContent(irow,148,50.0);
250   }
251
252   his->SetXTitle("z (cm)");
253   his->SetYTitle("y (cm)");
254   his->SetMaximum(50);
255   his->SetMinimum(0.0);
256   his->SetStats(0);
257
258   return his;
259
260 }
261
262 //_____________________________________________________________________________
263 void AliTRDCalPadStatus::PlotHistos2DSm(Int_t sm, const Char_t *name)
264 {
265   //
266   // Make 2D graph
267   //
268
269   gStyle->SetPalette(1);
270   TCanvas *c1 = new TCanvas(name,name,50,50,600,800);
271   c1->Divide(3,2);
272   c1->cd(1);
273   MakeHisto2DSmPl(sm,0)->Draw("colz");
274   c1->cd(2);
275   MakeHisto2DSmPl(sm,1)->Draw("colz");
276   c1->cd(3);
277   MakeHisto2DSmPl(sm,2)->Draw("colz");
278   c1->cd(4);
279   MakeHisto2DSmPl(sm,3)->Draw("colz");
280   c1->cd(5);
281   MakeHisto2DSmPl(sm,4)->Draw("colz");
282   c1->cd(6);
283   MakeHisto2DSmPl(sm,5)->Draw("colz");
284
285 }