1 /**************************************************************************
2 * Copyright(c) 2007-2009, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
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 **************************************************************************/
18 // *************************************************************
19 // Checks the quality assurance
20 // by comparing with reference data
22 // -------------------------------------------------------------
23 // W. Ferrarese + P. Cerello Feb 2008
25 // SSD QA part: P. Christakoglou - NIKHEF/UU
27 // --- ROOT system ---
32 // --- Standard library ---
34 // --- AliRoot header files ---
36 #include "AliITSmodule.h"
37 #include "AliITShit.h"
38 #include "AliITSdigit.h"
39 #include "AliITSpListItem.h"
41 #include "AliITSQADataMakerSim.h"
42 #include "AliITSQASSDDataMakerSim.h"
45 #include "AliQAChecker.h"
46 #include "AliRawReader.h"
48 ClassImp(AliITSQASSDDataMakerSim)
50 //____________________________________________________________________________
51 AliITSQASSDDataMakerSim::AliITSQASSDDataMakerSim(AliITSQADataMakerSim *aliITSQADataMakerSim) :
53 fAliITSQADataMakerSim(aliITSQADataMakerSim),
62 //ctor used to discriminate OnLine-Offline analysis
65 //____________________________________________________________________________
66 AliITSQASSDDataMakerSim::AliITSQASSDDataMakerSim(const AliITSQASSDDataMakerSim& qadm) :
68 fAliITSQADataMakerSim(qadm.fAliITSQADataMakerSim),
69 //fSSDhTask(qadm.fSSDhTask),
70 fSSDhHTask(qadm.fSSDhHTask),
71 fSSDhSTask(qadm.fSSDhSTask),
72 fSSDhDTask(qadm.fSSDhDTask),
73 fGenOffsetH(qadm.fGenOffsetH),
74 fGenOffsetS(qadm.fGenOffsetS),
75 fGenOffsetD(qadm.fGenOffsetD)
78 fAliITSQADataMakerSim->SetName((const char*)qadm.fAliITSQADataMakerSim->GetName()) ;
79 fAliITSQADataMakerSim->SetTitle((const char*)qadm.fAliITSQADataMakerSim->GetTitle());
82 //__________________________________________________________________
83 AliITSQASSDDataMakerSim& AliITSQASSDDataMakerSim::operator = (const AliITSQASSDDataMakerSim& qac ) {
85 this->~AliITSQASSDDataMakerSim();
86 new(this) AliITSQASSDDataMakerSim(qac);
90 //____________________________________________________________________________
91 void AliITSQASSDDataMakerSim::StartOfDetectorCycle() {
92 //Detector specific actions at start of cycle
93 AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM::Start of SSD Cycle\n");
96 //____________________________________________________________________________
97 void AliITSQASSDDataMakerSim::EndOfDetectorCycle(AliQAv1::TASKINDEX_t /*task*/, TObjArray* /*list*/) {
98 // launch the QA checking
99 AliDebug(AliQAv1::GetQADebugLevel(),"AliITSDM instantiates checker with Run(AliQAv1::kITS, task, list)\n");
101 // AliQAChecker::Instance()->Run( AliQAv1::kITS , task, list);
104 //____________________________________________________________________________
105 Int_t AliITSQASSDDataMakerSim::InitDigits() {
106 // Initialization for DIGIT data - SSD -
107 const Bool_t expert = kTRUE ;
108 const Bool_t image = kTRUE ;
110 // fGenOffsetD = (fAliITSQADataMakerSim->fDigitsQAList[AliRecoParam::kDefault])->GetEntries();
113 TH1F *fHistSSDModule = new TH1F("fHistSSDDigitsModule",
114 "SSD Digits Module;SSD Module Number;N_{DIGITS}",
116 rv = fAliITSQADataMakerSim->Add2DigitsList(fHistSSDModule,
117 fGenOffsetD + 0, !expert, image);
119 TH2F *fHistSSDModuleStrip = new TH2F("fHistSSDDigitsModuleStrip",
120 "SSD Digits Module Strip;N_{Strip};N_{Module}",
121 1540,0,1540,1698,499.5,2197.5);
122 rv = fAliITSQADataMakerSim->Add2DigitsList(fHistSSDModuleStrip,
123 fGenOffsetD + 1, !expert, image);
126 AliDebug(AliQAv1::GetQADebugLevel(),Form("%d SSD Digits histograms booked\n",fSSDhDTask));
130 //____________________________________________________________________________
131 Int_t AliITSQASSDDataMakerSim::MakeDigits(TTree *digits) {
132 // Fill QA for DIGIT - SSD -
135 AliITS *fITS = (AliITS*)gAlice->GetModule("ITS");
136 fITS->SetTreeAddress();
137 TClonesArray *iSSDdigits = fITS->DigitsAddress(2);
138 for(Int_t iModule = 500; iModule < 2198; iModule++) {
140 digits->GetEvent(iModule);
141 Int_t ndigits = iSSDdigits->GetEntries();
142 fAliITSQADataMakerSim->GetDigitsData(fGenOffsetD + 0)->Fill(iModule,ndigits);
144 AliDebug(AliQAv1::GetQADebugLevel(),Form("Module: %d - Digits: %d",iModule,ndigits));
146 for (Int_t iDigit = 0; iDigit < ndigits; iDigit++) {
147 AliITSdigit *dig = (AliITSdigit*)iSSDdigits->UncheckedAt(iDigit);
148 Int_t fStripNumber = (dig->GetCoord1() == 0) ? dig->GetCoord2() : dig->GetCoord2() + fgkNumberOfPSideStrips;
149 ((TH2F *)fAliITSQADataMakerSim->GetDigitsData(fGenOffsetD + 1))->Fill(fStripNumber,iModule,dig->GetSignal());
155 //____________________________________________________________________________
156 Int_t AliITSQASSDDataMakerSim::InitSDigits() {
157 // Initialization for SDIGIT data - SSD -
158 const Bool_t expert = kTRUE ;
159 const Bool_t image = kTRUE ;
161 //fGenOffsetS = (fAliITSQADataMakerSim->fSDigitsQAList[AliRecoParam::kDefault])->GetEntries();
164 TH1F *fHistSSDModule = new TH1F("fHistSSDSDigitsModule",
165 "SSD SDigits Module;SSD Module Number;N_{SDIGITS}",
167 rv = fAliITSQADataMakerSim->Add2SDigitsList(fHistSSDModule,
168 fGenOffsetS + 0, !expert, image);
171 AliDebug(AliQAv1::GetQADebugLevel(),Form("%d SSD SDigits histograms booked\n",fSSDhSTask));
175 //____________________________________________________________________________
176 Int_t AliITSQASSDDataMakerSim::MakeSDigits(TTree *sdigits) {
177 // Fill QA for SDIGIT - SSD -
180 static TClonesArray iSSDEmpty("AliITSpListItem",10000);
182 TClonesArray *iSSDsdigits = &iSSDEmpty;
184 AliDebug(AliQAv1::GetQADebugLevel(), Form("Trying to access the sdigits histogram: %d\n",fGenOffsetS));
186 TBranch *brchSDigits = sdigits->GetBranch("ITS");
187 brchSDigits->SetAddress(&iSSDsdigits);
188 for(Int_t iModule = 500; iModule < 2198; iModule++) {
189 iSSDsdigits->Clear();
190 sdigits->GetEvent(iModule);
191 Int_t ndigits = iSSDsdigits->GetEntries();
192 fAliITSQADataMakerSim->GetSDigitsData(fGenOffsetS + 0)->Fill(iModule,ndigits);
194 AliDebug(AliQAv1::GetQADebugLevel(),Form("Module: %d - Digits: %d",iModule,ndigits));
196 for (Int_t iDigit = 0; iDigit < ndigits; iDigit++) {
197 AliITSpListItem *dig=(AliITSpListItem*)iSSDsdigits->At(iDigit);
204 //____________________________________________________________________________
205 Int_t AliITSQASSDDataMakerSim::InitHits() {
206 // Initialization for HITS data - SSD -
207 const Bool_t expert = kTRUE ;
208 const Bool_t image = kTRUE ;
211 //fGenOffsetH = (fAliITSQADataMakerSim->fHitsQAList[fEventSpecie])->GetEntries();
214 TH1F *fHistSSDModule = new TH1F("fHistSSDHitsModule",
215 "SSD Hits Module;SDD Module Number;N_{HITS}",
217 rv = fAliITSQADataMakerSim->Add2HitsList(fHistSSDModule,
218 fGenOffsetH + 0, !expert, image);
220 TH1F *fHistSSDGlobalX = new TH1F("fHistSSDHitsGlobalX",
221 "SSD Hits Global X;x [cm];Entries",
223 rv = fAliITSQADataMakerSim->Add2HitsList(fHistSSDGlobalX,
224 fGenOffsetH + 1, !expert, image);
226 TH1F *fHistSSDGlobalY = new TH1F("fHistSSDHitsGlobalY",
227 "SSD Hits Global Y;y [cm];Entries",
229 rv = fAliITSQADataMakerSim->Add2HitsList(fHistSSDGlobalY,
230 fGenOffsetH + 2, !expert, image);
232 TH1F *fHistSSDGlobalZ = new TH1F("fHistSSDHitsGlobalZ",
233 "SSD Hits Global Z ;z [cm];Entries",
235 rv = fAliITSQADataMakerSim->Add2HitsList(fHistSSDGlobalZ,
236 fGenOffsetH + 3, !expert, image);
238 TH1F *fHistSSDLocalX = new TH1F("fHistSSDHitsLocalX",
239 "SSD Hits Local X;x [cm];Entries",
241 rv = fAliITSQADataMakerSim->Add2HitsList(fHistSSDLocalX,
242 fGenOffsetH + 4, !expert, image);
244 TH1F *fHistSSDLocalY = new TH1F("fHistSSDHitsLocalY",
245 "SSD Hits Local Y;y [cm];Entries",
247 rv = fAliITSQADataMakerSim->Add2HitsList(fHistSSDLocalY,
248 fGenOffsetH + 5, !expert, image);
250 TH1F *fHistSSDLocalZ = new TH1F("fHistSSDHitsLocalZ",
251 "SSD Hits Local Z;z [cm];Entries",
253 rv = fAliITSQADataMakerSim->Add2HitsList(fHistSSDLocalZ,
254 fGenOffsetH + 6, !expert, image);
256 TH1F *fHistSSDIonization = new TH1F("fHistSSDHitsIonization",
257 "SSD Hits Ionization;log(dE/dx) [KeV];N_{Hits}",
259 rv = fAliITSQADataMakerSim->Add2HitsList(fHistSSDIonization,
260 fGenOffsetH + 7, !expert, image);
262 TH2F *fHistSSDGlobalXY = new TH2F("fHistSSDHitsGlobalXY",
263 "SSD Hits Global XY;x [cm];y [cm]",
266 rv = fAliITSQADataMakerSim->Add2HitsList(fHistSSDGlobalXY,
267 fGenOffsetH + 8, !expert, image);
270 AliDebug(AliQAv1::GetQADebugLevel(),Form("%d SSD Hits histograms booked\n",fSSDhHTask));
275 //____________________________________________________________________________
276 Int_t AliITSQASSDDataMakerSim::MakeHits(TTree *hits) {
277 // Fill QA for HITS - SSD -
280 AliITS *fITS = (AliITS*)gAlice->GetModule("ITS");
281 fITS->SetTreeAddress();
283 fITS->InitModules(-1,nmodules);
284 fITS->FillModules(hits,0);
285 for(Int_t iModule = 500; iModule < 2198; iModule++) {
286 AliITSmodule *module = fITS->GetModule(iModule);
287 TObjArray *arrHits = module->GetHits();
288 Int_t nhits = arrHits->GetEntriesFast();
290 AliDebug(AliQAv1::GetQADebugLevel(),Form("Module: %d - Hits: %d",iModule,nhits));
291 for (Int_t iHit = 0; iHit < nhits; iHit++) {
292 AliITShit *hit = (AliITShit*) arrHits->At(iHit);
294 fAliITSQADataMakerSim->GetHitsData(fGenOffsetH + 0)->Fill(iModule);
295 fAliITSQADataMakerSim->GetHitsData(fGenOffsetH + 1)->Fill(hit->GetXG());
296 fAliITSQADataMakerSim->GetHitsData(fGenOffsetH + 2)->Fill(hit->GetYG());
297 fAliITSQADataMakerSim->GetHitsData(fGenOffsetH + 3)->Fill(hit->GetZG());
298 fAliITSQADataMakerSim->GetHitsData(fGenOffsetH + 4)->Fill(hit->GetXL());
299 fAliITSQADataMakerSim->GetHitsData(fGenOffsetH + 5)->Fill(hit->GetYL());
300 fAliITSQADataMakerSim->GetHitsData(fGenOffsetH + 6)->Fill(hit->GetZL());
301 if(hit->GetIonization())
302 fAliITSQADataMakerSim->GetHitsData(fGenOffsetH + 7)->Fill(TMath::Log10(hit->GetIonization()));
303 fAliITSQADataMakerSim->GetHitsData(fGenOffsetH + 8)->Fill(hit->GetXG(),hit->GetYG());
309 //____________________________________________________________________________
310 Int_t AliITSQASSDDataMakerSim::GetOffset(AliQAv1::TASKINDEX_t task){
311 // Returns histogram offset according to the specified task
313 if( task == AliQAv1::kHITS){
316 else if( task == AliQAv1::kSDIGITS) {
319 else if( task == AliQAv1::kDIGITS) {
323 AliInfo("No task has been selected. TaskHisto set to zero.\n");
330 //____________________________________________________________________________
331 void AliITSQASSDDataMakerSim::SetOffset(AliQAv1::TASKINDEX_t task, Int_t offset){
332 // Returns histogram offset according to the specified task
333 if( task == AliQAv1::kHITS){
334 fGenOffsetH = offset;
336 else if( task == AliQAv1::kSDIGITS) {
337 fGenOffsetS = offset;
339 else if( task == AliQAv1::kDIGITS) {
340 fGenOffsetD = offset;
343 AliInfo("No task has been selected. TaskHisto set to zero.\n");
347 //____________________________________________________________________________
348 Int_t AliITSQASSDDataMakerSim::GetTaskHisto(AliQAv1::TASKINDEX_t task) {
349 // Returns the number of booked histograms for the selected task
351 if( task == AliQAv1::kHITS) {
352 histotot=fSSDhHTask ;
354 else if( task == AliQAv1::kSDIGITS) {
357 else if( task == AliQAv1::kDIGITS) {
358 histotot=fSSDhDTask ;
361 AliInfo("No task has been selected. TaskHisto set to zero.\n");