Changes needed by the following commit: coding convention for type (_t) and access...
[u/mrichter/AliRoot.git] / ITS / AliITSQASDDDataMakerRec.cxx
1 /**************************************************************************\r
2  * Copyright(c) 2007-2009, ALICE Experiment at CERN, All rights reserved. *\r
3  *                                                                        *\r
4  * Author: The ALICE Off-line Project.                                    *\r
5  * Contributors are mentioned in the code where appropriate.              *\r
6  *                                                                        *\r
7  * Permission to use, copy, modify and distribute this software and its   *\r
8  * documentation strictly for non-commercial purposes is hereby granted   *\r
9  * without fee, provided that the above copyright notice appears in all   *\r
10  * copies and that both the copyright notice and this permission notice   *\r
11  * appear in the supporting documentation. The authors make no claims     *\r
12  * about the suitability of this software for any purpose. It is          *\r
13  * provided "as is" without express or implied warranty.                  *\r
14  **************************************************************************/\r
15 \r
16 /* $Id$ */\r
17 \r
18 //  *************************************************************\r
19 //  Checks the quality assurance \r
20 //  by comparing with reference data\r
21 //  contained in a DB\r
22 //  -------------------------------------------------------------\r
23 //  W. Ferrarese + P. Cerello Feb 2008\r
24 //  INFN Torino\r
25 \r
26 // --- ROOT system ---\r
27 #include <TH2D.h>\r
28 #include <TBranch.h>\r
29 #include <TTree.h>\r
30 #include <TGaxis.h>\r
31 #include <TMath.h>\r
32 // --- Standard library ---\r
33 \r
34 // --- AliRoot header files ---\r
35 #include "AliITSQASDDDataMakerRec.h"\r
36 #include "AliLog.h"\r
37 #include "AliQA.h"\r
38 #include "AliQAChecker.h"\r
39 #include "AliRawReader.h"\r
40 #include "AliITSRawStreamSDD.h"\r
41 #include "AliITSRecPoint.h"\r
42 #include "AliITSgeomTGeo.h"\r
43 \r
44 #include "AliCDBManager.h"\r
45 #include "AliCDBStorage.h"\r
46 #include "AliCDBEntry.h"\r
47 \r
48 \r
49 ClassImp(AliITSQASDDDataMakerRec)\r
50 \r
51 //____________________________________________________________________________ \r
52 AliITSQASDDDataMakerRec::AliITSQASDDDataMakerRec(AliITSQADataMakerRec *aliITSQADataMakerRec, Bool_t kMode, Short_t ldc) :\r
53 TObject(),\r
54 fAliITSQADataMakerRec(aliITSQADataMakerRec),\r
55 fkOnline(kMode),\r
56 fLDC(ldc),\r
57 fSDDhRaws(0),\r
58 fSDDhRecs(0),\r
59 fRawsOffset(0),\r
60 fRecsOffset(0),\r
61 fSDDDDLModuleMap(0)\r
62 {\r
63   //ctor used to discriminate OnLine-Offline analysis\r
64   if(fLDC < 0 || fLDC > 4) {\r
65         AliError("Error: LDC number out of range; return\n");\r
66   }\r
67   for(Int_t i=0;i<2*fgknSDDmodules;i++){\r
68     fModuleChargeMap[i] = NULL;\r
69   }\r
70 }\r
71 \r
72 //____________________________________________________________________________ \r
73 AliITSQASDDDataMakerRec::AliITSQASDDDataMakerRec(const AliITSQASDDDataMakerRec& qadm) :\r
74 TObject(),\r
75 fAliITSQADataMakerRec(qadm.fAliITSQADataMakerRec),\r
76 fkOnline(qadm.fkOnline),\r
77 fLDC(qadm.fLDC),\r
78 fSDDhRaws(qadm.fSDDhRaws),\r
79 fSDDhRecs(qadm.fSDDhRecs),\r
80 fRawsOffset(qadm.fRawsOffset),\r
81 fRecsOffset(qadm.fRecsOffset),\r
82 fSDDDDLModuleMap(0)\r
83 {\r
84   //copy ctor \r
85   fAliITSQADataMakerRec->SetName((const char*)qadm.fAliITSQADataMakerRec->GetName()) ; \r
86   fAliITSQADataMakerRec->SetTitle((const char*)qadm.fAliITSQADataMakerRec->GetTitle());\r
87 }\r
88 \r
89 //____________________________________________________________________________ \r
90 AliITSQASDDDataMakerRec::~AliITSQASDDDataMakerRec(){\r
91   // destructor\r
92 \r
93   for(Int_t i=0;i<2*fgknSDDmodules;i++){\r
94     if(fModuleChargeMap[i]) delete fModuleChargeMap[i];\r
95   }\r
96   \r
97 }\r
98 //__________________________________________________________________\r
99 AliITSQASDDDataMakerRec& AliITSQASDDDataMakerRec::operator = (const AliITSQASDDDataMakerRec& qac )\r
100 {\r
101   // Equal operator.\r
102   this->~AliITSQASDDDataMakerRec();\r
103   new(this) AliITSQASDDDataMakerRec(qac);\r
104   return *this;\r
105 }\r
106 \r
107 //____________________________________________________________________________ \r
108 void AliITSQASDDDataMakerRec::StartOfDetectorCycle()\r
109 {\r
110   //Detector specific actions at start of cycle\r
111   AliDebug(1,"AliITSQADM::Start of SDD Cycle\n");\r
112 }\r
113 \r
114 //____________________________________________________________________________ \r
115 void AliITSQASDDDataMakerRec::EndOfDetectorCycle(AliQA::TASKINDEX_t /*task*/, TObjArray* /*list*/)\r
116 {\r
117   // launch the QA checking\r
118   AliDebug(1,"AliITSDM instantiates checker with Run(AliQA::kITS, task, list)\n"); \r
119   \r
120   //AliQAChecker::Instance()->Run( AliQA::kITS , task, list);\r
121 }\r
122 \r
123 //____________________________________________________________________________ \r
124 void AliITSQASDDDataMakerRec::InitRaws()\r
125\r
126   // Initialization for RAW data - SDD -\r
127   fRawsOffset = (fAliITSQADataMakerRec->fRawsQAList)->GetEntries();\r
128 \r
129   AliCDBEntry *ddlMapSDD = AliCDBManager::Instance()->Get("ITS/Calib/DDLMapSDD");\r
130   if( !ddlMapSDD){\r
131     AliFatal("Calibration object retrieval failed! ");\r
132   }  \r
133   AliITSDDLModuleMapSDD *ddlsdd=(AliITSDDLModuleMapSDD*)ddlMapSDD->GetObject();\r
134   ddlMapSDD->SetOwner(kTRUE);\r
135   fSDDDDLModuleMap = ddlsdd;\r
136  \r
137   Int_t lay, lad, det;\r
138   Int_t LAY = -1;  //, LAD = -1;\r
139   char hname0[50];\r
140   Int_t indexlast = 0;\r
141   Int_t index1 = 0;\r
142 \r
143   if(fLDC == 1 || fLDC == 2) LAY = 2;\r
144   if(fLDC == 3 || fLDC == 4) LAY = 3;\r
145   \r
146   if(fkOnline) {\r
147     AliInfo("Book Online Histograms for SDD\n");\r
148   }\r
149   else {\r
150     AliInfo("Book Offline Histograms for SDD\n ");\r
151   }\r
152   TH1F *h0 = new TH1F("ModPattern","Modules pattern",fgknSDDmodules,-0.5,259.5);\r
153   h0->GetXaxis()->SetTitle("Module Number");\r
154   h0->GetYaxis()->SetTitle("Counts");\r
155   fAliITSQADataMakerRec->Add2RawsList(h0,0+fRawsOffset);\r
156   fSDDhRaws++;\r
157   if(fLDC==0 || fLDC==1 || fLDC==2){\r
158     TH1F *h1 = new TH1F("LadPatternL3","Ladder pattern L3",14,0.5,14.5);  \r
159     h1->GetXaxis()->SetTitle("Ladder Number on Lay3");\r
160     h1->GetYaxis()->SetTitle("Counts");\r
161     fAliITSQADataMakerRec->Add2RawsList(h1,1+fRawsOffset);\r
162     fSDDhRaws++;\r
163   }     \r
164   if(fLDC==0 || fLDC==3 || fLDC==4){\r
165     TH1F *h2 = new TH1F("LadPatternL4","Ladder pattern L4",22,0.5,22.5);  \r
166     h2->GetXaxis()->SetTitle("Ladder Number on Lay4");\r
167     h2->GetYaxis()->SetTitle("Counts");\r
168     fAliITSQADataMakerRec->Add2RawsList(h2,2+fRawsOffset);\r
169     fSDDhRaws++;\r
170   }\r
171   if(fLDC==0 || fLDC==1 || fLDC==2){\r
172     TH1D *h3[fgkLADDonLAY3] ; \r
173     for(Int_t i=1; i<=fgkLADDonLAY3; i++) {\r
174       sprintf(hname0,"ModPattern_L3_%d",i);\r
175       h3[i-1] = new TH1D(hname0,hname0,6,0.5,6.5);\r
176       h3[i-1]->GetXaxis()->SetTitle("Channel Number");\r
177       h3[i-1]->GetYaxis()->SetTitle("Counts");\r
178       fAliITSQADataMakerRec->Add2RawsList(h3[i-1],i-1+3+fRawsOffset);\r
179       fSDDhRaws++;\r
180     }\r
181   }\r
182   if(fLDC==0 || fLDC==3 || fLDC==4){\r
183     TH1D *h4[fgkLADDonLAY4] ; \r
184     for(Int_t i=1; i<=fgkLADDonLAY4; i++) {\r
185       sprintf(hname0,"ModPattern_L4_%d",i);\r
186       h4[i-1] = new TH1D(hname0,hname0,8,0.5,8.5);\r
187       h4[i-1]->GetXaxis()->SetTitle("Channel Number");\r
188       h4[i-1]->GetYaxis()->SetTitle("Counts");\r
189       fAliITSQADataMakerRec->Add2RawsList(h4[i-1],i-1+17+fRawsOffset);\r
190       fSDDhRaws++;\r
191     }\r
192   }\r
193 \r
194   Int_t indexlast1 = 0;\r
195   Int_t indexlast2 = 0;\r
196 \r
197   if(fkOnline) {\r
198     indexlast = 0;\r
199     TH1D *h5[fgknSDDmodules*2] ; \r
200     index1 = 0;\r
201 \r
202 /*\r
203     if(fLDC == 0 || fLDC == 1 || fLDC == 4){\r
204       \r
205         for(Int_t lad=0; lad<fgkLADDonLAY3; lad++){\r
206         for(Int_t ichan=0; det<3; det++){\r
207           for(Int_t iside=0; iside<fgknSide; iside++){\r
208             sprintf(hname0,"ProjYMap_L%d_%d_%d_%d",lay,lad,det+sideflag*3,iside);  //need to define "lay"\r
209             h5[index1] = new TH1D(hname0,hname0,256,-0.5,255.5);\r
210             h5[index1]->GetXaxis()->SetTitle("Anode Number");\r
211             h5[index1]->GetYaxis()->SetTitle("Counts");\r
212             fAliITSQADataMakerRec->Add2RawsList(h5[index1],index1+39+fRawsOffset);\r
213                 fSDDhRaws++;\r
214             index1++;\r
215             indexlast = index1+39;\r
216           }\r
217         }       \r
218       }\r
219     }\r
220 */  \r
221     for(Int_t moduleSDD=0; moduleSDD<fgknSDDmodules; moduleSDD++){\r
222       if((moduleSDD >= 0 && moduleSDD < 36) || (moduleSDD >= 84 && moduleSDD < 180)) {\r
223         for(Int_t iside=0; iside<fgknSide; iside++){\r
224           AliITSgeomTGeo::GetModuleId(moduleSDD+fgkmodoffset, lay, lad, det);\r
225           //if(fLDC == 0 || ((fLDC == 1 || fLDC == 2) && lay == 2) || ((fLDC == 3 || fLDC == 4) && lay == 3)\r
226           sprintf(hname0,"ProjYMap_L%d_%d_%d_%d",lay,lad,det,iside);\r
227           h5[index1] = new TH1D(hname0,hname0,256,-0.5,255.5);\r
228           h5[index1]->GetXaxis()->SetTitle("Anode Number");\r
229           h5[index1]->GetYaxis()->SetTitle("Counts");\r
230           fAliITSQADataMakerRec->Add2RawsList(h5[index1],index1+39+fRawsOffset);\r
231           fSDDhRaws++;\r
232           index1++;\r
233           indexlast = index1+39;\r
234         }\r
235       }\r
236     }\r
237     \r
238     TH1D *h6[fgknSDDmodules*8] ; \r
239     Int_t indextot = 0;\r
240     indexlast1 = 0;\r
241     for(Int_t htype=0; htype<4; htype++){\r
242       for(Int_t moduleSDD=0; moduleSDD<fgknSDDmodules; moduleSDD++){\r
243         if((moduleSDD >= 0 && moduleSDD < 36) || (moduleSDD >= 84 && moduleSDD < 180)) {         \r
244           for(Int_t iside=0; iside<fgknSide; iside++){\r
245             AliITSgeomTGeo::GetModuleId(moduleSDD+fgkmodoffset, lay, lad, det);\r
246             if(htype == 0) sprintf(hname0,"CountsVSAnode_L%d_%d_%d_%d",lay,lad,det,iside);\r
247             if(htype == 1) sprintf(hname0,"ChargeVSAnode_L%d_%d_%d_%d",lay,lad,det,iside);\r
248             if(htype == 2) sprintf(hname0,"CountsVSTbin_L%d_%d_%d_%d",lay,lad,det,iside);\r
249             if(htype == 3) sprintf(hname0,"ChargeVSTbin_L%d_%d_%d_%d",lay,lad,det,iside);\r
250             h6[indextot] = new TH1D(hname0,hname0,256,-0.5,255.5);\r
251             if(htype == 0){\r
252               h6[indextot]->GetXaxis()->SetTitle("Anode Number");\r
253               h6[indextot]->GetYaxis()->SetTitle("Counts");\r
254             }\r
255             if(htype == 1){\r
256               h6[indextot]->GetXaxis()->SetTitle("Anode Number");\r
257               h6[indextot]->GetYaxis()->SetTitle("Charge (ADC)");\r
258             }\r
259             if(htype == 2){\r
260               h6[indextot]->GetXaxis()->SetTitle("Time bin");\r
261               h6[indextot]->GetYaxis()->SetTitle("Counts");\r
262             }\r
263             if(htype == 3){\r
264               h6[indextot]->GetXaxis()->SetTitle("Time bin");\r
265               h6[indextot]->GetYaxis()->SetTitle("Charge (ADC)");\r
266             }\r
267             fAliITSQADataMakerRec->Add2RawsList(h6[indextot],indexlast + indextot + fRawsOffset);\r
268             fSDDhRaws++;\r
269             indextot++;\r
270             indexlast1 = indexlast + indextot;\r
271           }\r
272         }\r
273       }\r
274     }\r
275 \r
276 \r
277     indexlast2 = 0;\r
278     index1 = 0;\r
279     char *hname[3];\r
280     for(Int_t i=0; i<3; i++) hname[i]= new char[50];\r
281     for(Int_t moduleSDD =0; moduleSDD<fgknSDDmodules; moduleSDD++){\r
282       if((moduleSDD >= 0 && moduleSDD < 36) || (moduleSDD >= 84 && moduleSDD < 180)) {\r
283         for(Int_t iside=0;iside<fgknSide;iside++){\r
284         AliITSgeomTGeo::GetModuleId(moduleSDD+fgkmodoffset, lay, lad, det);\r
285         sprintf(hname[0],"chargeMap_L%d_%d_%d_%d",lay,lad,det,iside);\r
286         sprintf(hname[1],"TotalCharge_L%d_%d_%d_%d",lay,lad,det,iside);\r
287         sprintf(hname[2],"hmonoDMap_L%d_%d_%d_%d",lay,lad,det,iside);\r
288         fModuleChargeMap[index1] = new TH2D(hname[0],hname[1],256,-0.5,255.5,256,-0.5,255.5);\r
289         fModuleChargeMap[index1]->GetXaxis()->SetTitle("Time Bin");\r
290         fModuleChargeMap[index1]->GetYaxis()->SetTitle("Anode");\r
291         fAliITSQADataMakerRec->Add2RawsList((new TH2D(*fModuleChargeMap[index1])),indexlast1 + index1 + fRawsOffset);\r
292         fSDDhRaws++;\r
293         index1++;        \r
294         indexlast2 = indexlast1 + index1;\r
295       }\r
296     }\r
297   }\r
298   \r
299   \r
300   TH2D *h7[fgknSDDmodules*2] ;\r
301   index1 = 0;\r
302   for(Int_t moduleSDD=0; moduleSDD<fgknSDDmodules; moduleSDD++){\r
303     if((moduleSDD >= 0 && moduleSDD < 36) || (moduleSDD >= 84 && moduleSDD < 180)) {\r
304       for(Int_t iside=0; iside<fgknSide; iside++){\r
305         AliITSgeomTGeo::GetModuleId(moduleSDD+fgkmodoffset, lay, lad, det);\r
306         sprintf(hname0,"Anode_vs_Charge_L%d_%d_%d_%d",lay,lad,det,iside);\r
307         h7[index1] = new TH2D(hname0,hname0,fgknSDDmodules*2,-0.5,-0.5+fgknSDDmodules*2,256,0.5,256.5);\r
308         h7[index1]->GetXaxis()->SetTitle("Charge Value (ADC)");\r
309         h7[index1]->GetYaxis()->SetTitle("Anode");\r
310         fAliITSQADataMakerRec->Add2RawsList(h7[index1],indexlast2 + index1 + fRawsOffset);\r
311         fSDDhRaws++;\r
312         index1++;\r
313       }\r
314     }\r
315   }\r
316   /*\r
317   TH2F *h8[3];\r
318   for(Int_t htype =0; htype<3; htype++){\r
319     if(htype == 0) sprintf(hname0,"BLcompare");\r
320     if(htype == 1) sprintf(hname0,"NoiseCompare");\r
321     if(htype == 2) sprintf(hname0,"VdriftCompare");\r
322     h8[htype]= new TH2F(hname0,hname0,520,0,260,100,-0.5,99.5);\r
323     if(htype == 0){\r
324       h8[htype]->GetXaxis()->SetTitle("Module Number");\r
325       h8[htype]->GetYaxis()->SetTitle("Time #");\r
326     }\r
327     if(htype == 1){\r
328       h8[htype]->GetXaxis()->SetTitle("Module Number");\r
329       h8[htype]->GetYaxis()->SetTitle("Time #)");\r
330     }\r
331     if(htype == 2){\r
332       h8[htype]->GetXaxis()->SetTitle("Module Number");\r
333       h8[htype]->GetYaxis()->SetTitle("Time #");\r
334     }\r
335 \r
336     fAliITSQADataMakerRec->Add2RawsList(h8[htype],indexlast2 + index1 + htype+1 + fRawsOffset);\r
337     fSDDhRaws++;\r
338     printf("added histos %c, at position %d \n", hname0, indexlast2 + index1 + htype+1 + fRawsOffset);\r
339   }\r
340   */\r
341 }  // kONLINE\r
342 \r
343 \r
344   AliDebug(1,Form("%d SDD Raws histograms booked\n",fSDDhRaws));\r
345 }\r
346 \r
347 \r
348 //____________________________________________________________________________\r
349 void AliITSQASDDDataMakerRec::MakeRaws(AliRawReader* rawReader)\r
350\r
351   // Fill QA for RAW - SDD -\r
352   if(rawReader->GetType() != 7) return;  // skips non physical triggers\r
353   AliDebug(1,"entering MakeRaws\n");                 \r
354   rawReader->SelectEquipment(17,fgkeqOffset,fgkeqOffset + fgkDDLidRange); \r
355 \r
356   /*\r
357   if(rawReader->GetEventId()!=fEvtId){\r
358     TFile *DAoutput = new TFile::Open(filename);\r
359     TH1F *BLhisto;\r
360     for(Int_t imod=0; imod<nSDDmodules; imod++){ \r
361       BLhisto = (TH1F*)DAoutput->Get("BLhistoname[imod]");\r
362       mean = BLhisto->take mean;\r
363       fAliITSQADataMakerRec->GetRawsData(i+1887)->Fill(mean, imod);\r
364       Noisehisto....;\r
365       Vdrifthisto...;\r
366       <Q>histo....;\r
367     }\r
368   }\r
369   fEvtId==rawReader->GetEventId();\r
370   */\r
371 \r
372   rawReader->Reset();                         \r
373   AliITSRawStreamSDD s(rawReader); \r
374   s.SetDDLModuleMap(fSDDDDLModuleMap);\r
375   Int_t lay, lad, det; \r
376 \r
377   Int_t index=0;\r
378   if(fkOnline) {\r
379     for(Int_t moduleSDD =0; moduleSDD<fgknSDDmodules; moduleSDD++){\r
380       if((moduleSDD >= 0 && moduleSDD < 36) || (moduleSDD >= 84 && moduleSDD < 180)) {\r
381         for(Int_t iside=0;iside<fgknSide;iside++) {\r
382           if(fSDDhRaws > 39+12 * 132 + index) fAliITSQADataMakerRec->GetRawsData(39+12 * 132 + index +fRawsOffset)->Reset();\r
383           index++;\r
384         }\r
385       }\r
386     }\r
387   }\r
388   Int_t cnt = 0;\r
389   Int_t ildcID = -1;\r
390   Int_t iddl = -1;\r
391   Int_t isddmod = -1;\r
392   Int_t coord1, coord2, signal, moduleSDD, ioffset, iorder, activeModule, index1 ,nBins , iside, bin ;\r
393   while(s.Next()) {\r
394     ildcID = rawReader->GetLDCId();\r
395     iddl = rawReader->GetDDLID() - fgkDDLIDshift;\r
396     isddmod = s.GetModuleNumber(iddl,s.GetCarlosId());\r
397     if(isddmod==-1){\r
398       AliDebug(1,Form("Found module with iddl: %d, s.GetCarlosId: %d \n",iddl,s.GetCarlosId() ));\r
399       continue;\r
400     }\r
401     if(s.IsCompletedModule()) {\r
402       AliDebug(1,Form("IsCompletedModule == KTRUE\n"));\r
403       continue;\r
404     } \r
405     \r
406     coord1 = s.GetCoord1();\r
407     coord2 = s.GetCoord2();\r
408     signal = s.GetSignal();\r
409     \r
410     moduleSDD = isddmod - fgkmodoffset;\r
411     if(moduleSDD < 0 || moduleSDD>fgknSDDmodules) {\r
412       AliDebug(1,Form( "Module SDD = %d, resetting it to 1 \n",moduleSDD));\r
413       moduleSDD = 1;\r
414     }\r
415     fAliITSQADataMakerRec->GetRawsData(0 +fRawsOffset)->Fill(moduleSDD); \r
416     \r
417     AliITSgeomTGeo::GetModuleId(isddmod, lay, lad, det);\r
418     ioffset = 3;\r
419     iorder = 1;\r
420     if(lay==4) { \r
421       ioffset += 14;\r
422       iorder = 2;   \r
423     } \r
424     fAliITSQADataMakerRec->GetRawsData(iorder +fRawsOffset)->Fill(lad);\r
425     fAliITSQADataMakerRec->GetRawsData(ioffset+lad-1 +fRawsOffset)->Fill(det); //-1 because ladder# starts from 1    \r
426     \r
427     Short_t iside = s.GetChannel();\r
428     activeModule = moduleSDD;\r
429     if(moduleSDD > 35) activeModule -= 48;\r
430     index1 = activeModule * 2 + iside;\r
431     \r
432     if(index1<0){\r
433       AliDebug(1,Form("Wrong index number %d - patched to 0\n",index1));\r
434       index1 = 0;\r
435     }\r
436     \r
437     if(fkOnline) {\r
438       if(fSDDhRaws > 39+12 * 132 + index1) {\r
439         \r
440         fAliITSQADataMakerRec->GetRawsData(39+ 2 * 132 + index1 +fRawsOffset)->Fill(coord1); \r
441         fAliITSQADataMakerRec->GetRawsData(39+ 4 * 132 + index1 +fRawsOffset)->Fill(coord1,signal); \r
442         fAliITSQADataMakerRec->GetRawsData(39+ 6 * 132 + index1 +fRawsOffset)->Fill(coord2);    \r
443         fAliITSQADataMakerRec->GetRawsData(39+ 8 * 132 + index1 +fRawsOffset)->Fill(coord2,signal); \r
444         fAliITSQADataMakerRec->GetRawsData(39+12 * 132 + index1 +fRawsOffset)->Fill(signal,coord1);\r
445         ((TH2D *)(fAliITSQADataMakerRec->GetRawsData(39+10 * 132 + index1 +fRawsOffset)))->Fill(coord2, coord1, signal);\r
446 \r
447       }\r
448     }\r
449     cnt++;\r
450     if(!(cnt%10000)) AliDebug(1,Form(" %d raw digits read",cnt));\r
451   }\r
452   AliDebug(1,Form("Event completed, %d raw digits read",cnt));  \r
453 \r
454   if(fkOnline) {\r
455     TH1D *ptr = NULL;\r
456     nBins = 256;\r
457     for(moduleSDD=0; moduleSDD<fgknSDDmodules; moduleSDD++){\r
458       if((moduleSDD >= 0 && moduleSDD < 36) || (moduleSDD >= 84 && moduleSDD < 180)) {\r
459         for(iside=0; iside<fgknSide; iside++){\r
460           activeModule = moduleSDD;\r
461           if(moduleSDD > 35) activeModule -= 48;\r
462           index1 = activeModule * 2 + iside;\r
463           if(fSDDhRaws > 39 + 2 * 132 + index1) {\r
464             ptr = ((TH2D *) (fAliITSQADataMakerRec->GetRawsData(39+10 * 132 + index1 +fRawsOffset)))->ProjectionY();\r
465             for(bin=0; bin<nBins; bin++) fAliITSQADataMakerRec->GetRawsData(index1+39 +fRawsOffset)->Fill(bin,ptr->GetBinContent(bin+1) );\r
466           }  \r
467         }\r
468       }\r
469     }\r
470     for(Int_t i=0; i<fSDDhRaws; i++){\r
471       Int_t entries = static_cast<Int_t>(fAliITSQADataMakerRec->GetRawsData(i +fRawsOffset)->GetEntries());\r
472       if(entries != 0)\r
473         AliDebug(1,Form("histo %d, name %s , entries %d ",i,fAliITSQADataMakerRec->GetRawsData(i +fRawsOffset)->GetName(),entries));\r
474     }\r
475   }\r
476 }\r
477 \r
478 //____________________________________________________________________________ \r
479 void AliITSQASDDDataMakerRec::InitRecPoints()\r
480 {\r
481   // Initialization for RECPOINTS - SDD -\r
482   fRecsOffset = (fAliITSQADataMakerRec->fRecPointsQAList)->GetEntries();\r
483   \r
484   TH1F * h0 = new TH1F("Lay3TotCh","Layer 3 total charge",1000,-0.5, 499.5);\r
485   h0->GetXaxis()->SetTitle("ADC value");\r
486   h0->GetYaxis()->SetTitle("Entries");\r
487   fAliITSQADataMakerRec->Add2RecPointsList(h0, 0 +fRecsOffset);\r
488   fSDDhRecs++;\r
489  \r
490   TH1F * h1 = new TH1F("Lay4TotCh","Layer 4 total charge",1000,-0.5, 499.5);\r
491   h1->GetXaxis()->SetTitle("ADC value");\r
492   h1->GetYaxis()->SetTitle("Entries");\r
493   fAliITSQADataMakerRec->Add2RecPointsList(h1, 1 +fRecsOffset);\r
494   fSDDhRecs++;\r
495 \r
496     \r
497   char hisnam[50];\r
498   TH1F *h2[3]; \r
499   for(Int_t i=1; i<=3; i++){\r
500     sprintf(hisnam,"Charge_L3_Strip%d",i);\r
501     h2[i] = new TH1F(hisnam,hisnam,1000,-0.5, 499.5);\r
502     fAliITSQADataMakerRec->Add2RecPointsList(h2[i],i+1 +fRecsOffset);\r
503     fSDDhRecs++;\r
504   }\r
505   \r
506   TH1F *h3[4]; \r
507   for(Int_t i=1; i<=4; i++){\r
508     sprintf(hisnam,"Charge_L4_Strip%d",i);\r
509     h3[i] = new TH1F(hisnam,hisnam,1000,-0.5, 499.5);\r
510     fAliITSQADataMakerRec->Add2RecPointsList(h3[i],i+4 +fRecsOffset);\r
511     fSDDhRecs++;\r
512   }\r
513   \r
514   TH1F *h4 = new TH1F("ModPatternRP","Modules pattern RP",fgknSDDmodules,239.5,499.5); \r
515   h4->GetXaxis()->SetTitle("Module number");\r
516   h4->GetYaxis()->SetTitle("Entries");\r
517   fAliITSQADataMakerRec->Add2RecPointsList(h4,9 +fRecsOffset);\r
518   fSDDhRecs++;\r
519   TH1F *h5 = new TH1F("ModPatternL3 RP","Ladder pattern L3 RP",14,0.5,14.5);  \r
520   h5->GetXaxis()->SetTitle("Ladder #, Layer 3");\r
521   h5->GetYaxis()->SetTitle("Entries");\r
522   fAliITSQADataMakerRec->Add2RecPointsList(h5,10 +fRecsOffset);\r
523   fSDDhRecs++;\r
524   TH1F *h6 = new TH1F("ModPatternL4 RP","Ladder pattern L4 RP",22,0.5,22.5); \r
525   h6->GetXaxis()->SetTitle("Ladder #, Layer 4");\r
526   h6->GetYaxis()->SetTitle("Entries");\r
527   fAliITSQADataMakerRec->Add2RecPointsList(h6,11 +fRecsOffset);\r
528   fSDDhRecs++;\r
529   TH2F *h7 = new TH2F("Local Coord Distrib","Local Coord Distrib",1000,-4,4,1000,-4,4);\r
530   h7->GetXaxis()->SetTitle("X local coord, drift, cm");\r
531   h7->GetYaxis()->SetTitle("Z local coord, anode, cm");\r
532   fAliITSQADataMakerRec->Add2RecPointsList(h7,12 +fRecsOffset);\r
533   fSDDhRecs++;\r
534   TH2F *h8 = new TH2F("Global Coord Distrib","Global Coord Distrib",6000,-30,30,6000,-30,30);\r
535   h8->GetYaxis()->SetTitle("Y glob coord, cm");\r
536   h8->GetXaxis()->SetTitle("X glob coord, cm");\r
537   fAliITSQADataMakerRec->Add2RecPointsList(h8,13 +fRecsOffset);\r
538   fSDDhRecs++;\r
539   \r
540   TH1F *h9[2]; \r
541   for(Int_t iLay=0; iLay<=1; iLay++){\r
542     sprintf(hisnam,"hr_Layer%d",iLay+3);\r
543     h9[iLay] = new TH1F(hisnam,hisnam,100,10.,30.);\r
544     h9[iLay]->GetXaxis()->SetTitle("r (cm)");\r
545     h9[iLay]->GetXaxis()->CenterTitle();\r
546     h9[iLay]->GetYaxis()->SetTitle("Entries");\r
547     fAliITSQADataMakerRec->Add2RecPointsList(h9[iLay],iLay+14 +fRecsOffset);\r
548     fSDDhRecs++;\r
549   }\r
550   TH1F *h10[2]; \r
551   for(Int_t iLay=0; iLay<=1; iLay++){\r
552     sprintf(hisnam,"hphi_Layer%d",iLay+3);\r
553     h10[iLay] = new TH1F(hisnam,hisnam,100,-TMath::Pi(),TMath::Pi());\r
554     h10[iLay]->GetXaxis()->SetTitle("#varphi (rad)");\r
555     h10[iLay]->GetXaxis()->CenterTitle();\r
556     h10[iLay]->GetYaxis()->SetTitle("Entries");\r
557     fAliITSQADataMakerRec->Add2RecPointsList(h10[iLay],iLay+16 +fRecsOffset);\r
558     fSDDhRecs++;\r
559   }\r
560 \r
561   AliDebug(1,Form("%d SDD Recs histograms booked\n",fSDDhRecs));\r
562 }\r
563 \r
564 //____________________________________________________________________________ \r
565 void AliITSQASDDDataMakerRec::MakeRecPoints(TTree * clustersTree)\r
566 {\r
567   // Fill QA for RecPoints - SDD -\r
568   Int_t lay, lad, det; \r
569   TBranch *branchRecP = clustersTree->GetBranch("ITSRecPoints");\r
570   if (!branchRecP) { \r
571     AliError("can't get the branch with the ITS clusters !");\r
572     return;\r
573   }\r
574   TClonesArray * recpoints = new TClonesArray("AliITSRecPoint") ;\r
575   branchRecP->SetAddress(&recpoints);\r
576   Int_t npoints = 0;      \r
577   Float_t cluglo[3]={0.,0.,0.}; \r
578   for(Int_t module=0; module<clustersTree->GetEntries();module++){\r
579     branchRecP->GetEvent(module);\r
580     npoints += recpoints->GetEntries();\r
581     AliITSgeomTGeo::GetModuleId(module, lay, lad, det);\r
582     //printf("modnumb %d, lay %d, lad %d, det %d \n",module, lay, lad, det);\r
583     \r
584     for(Int_t j=0;j<recpoints->GetEntries();j++){\r
585       AliITSRecPoint *recp = (AliITSRecPoint*)recpoints->At(j);    \r
586       fAliITSQADataMakerRec->GetRecPointsData(9 +fRecsOffset)->Fill(module);\r
587       recp->GetGlobalXYZ(cluglo);\r
588       Float_t rad=TMath::Sqrt(cluglo[0]*cluglo[0]+cluglo[1]*cluglo[1]); \r
589       Float_t phi=TMath::ATan2(cluglo[1],cluglo[0]);\r
590       if(recp->GetLayer() ==2) {\r
591         fAliITSQADataMakerRec->GetRecPointsData(0 +fRecsOffset)->Fill(recp->GetQ()) ;\r
592         fAliITSQADataMakerRec->GetRecPointsData(10 +fRecsOffset)->Fill(lad);\r
593         fAliITSQADataMakerRec->GetRecPointsData(14 +fRecsOffset)->Fill(rad);\r
594         fAliITSQADataMakerRec->GetRecPointsData(16 +fRecsOffset)->Fill(phi);\r
595         fAliITSQADataMakerRec->GetRecPointsData(9 +fRecsOffset)->Fill(module);\r
596         fAliITSQADataMakerRec->GetRecPointsData(12 +fRecsOffset)->Fill(recp->GetDetLocalX(),recp->GetDetLocalZ());\r
597         fAliITSQADataMakerRec->GetRecPointsData(13 +fRecsOffset)->Fill(cluglo[0],cluglo[1]);\r
598       }\r
599       else if(recp->GetLayer() ==3) {\r
600         fAliITSQADataMakerRec->GetRecPointsData(1 +fRecsOffset)->Fill(recp->GetQ()) ;\r
601         fAliITSQADataMakerRec->GetRecPointsData(11 +fRecsOffset)->Fill(lad);\r
602         fAliITSQADataMakerRec->GetRecPointsData(15 +fRecsOffset)->Fill(rad);\r
603         fAliITSQADataMakerRec->GetRecPointsData(17 +fRecsOffset)->Fill(phi);\r
604         fAliITSQADataMakerRec->GetRecPointsData(9 +fRecsOffset)->Fill(module);\r
605         fAliITSQADataMakerRec->GetRecPointsData(12 +fRecsOffset)->Fill(recp->GetDetLocalX(),recp->GetDetLocalZ());\r
606         fAliITSQADataMakerRec->GetRecPointsData(13 +fRecsOffset)->Fill(cluglo[0],cluglo[1]);\r
607       }\r
608     }\r
609   }\r
610   recpoints->Delete();\r
611   delete recpoints;\r
612 \r
613 }\r
614 \r