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 **************************************************************************/
17 // *************************************************************
18 // Checks the quality assurance
19 // by comparing with reference data
21 // -------------------------------------------------------------
22 // W. Ferrarese + P. Cerello INFN Torino Feb 2008
23 // M. Nicassio D. Elia INFN Bari April 2008
24 // maria.nicassio@ba.infn.it
26 // --- ROOT system ---
30 // --- Standard library ---
32 // --- AliRoot header files ---
34 #include "AliITSQADataMakerSim.h"
35 #include "AliITSQASPDDataMakerSim.h"
37 #include "AliQAChecker.h"
38 #include "AliITSdigit.h"
39 #include "AliITSdigitSPD.h"
41 #include "AliITSmodule.h"
42 #include "AliITShit.h"
43 #include "AliITSLoader.h"
44 #include "AliRunLoader.h"
46 ClassImp(AliITSQASPDDataMakerSim)
48 //____________________________________________________________________________
49 AliITSQASPDDataMakerSim::AliITSQASPDDataMakerSim(AliITSQADataMakerSim *aliITSQADataMakerSim) :
51 fAliITSQADataMakerSim(aliITSQADataMakerSim),
59 //ctor used to discriminate OnLine-Offline analysis
62 //____________________________________________________________________________
63 AliITSQASPDDataMakerSim::AliITSQASPDDataMakerSim(const AliITSQASPDDataMakerSim& qadm) :
65 fAliITSQADataMakerSim(qadm.fAliITSQADataMakerSim),
66 fSPDhDigits(qadm.fSPDhDigits),
67 fSPDhSDigits(qadm.fSPDhSDigits),
68 fSPDhHits(qadm.fSPDhHits),
69 fDigitsOffset(qadm.fDigitsOffset),
70 fSDigitsOffset(qadm.fSDigitsOffset),
71 fHitsOffset(qadm.fHitsOffset)
74 fAliITSQADataMakerSim->SetName((const char*)qadm.fAliITSQADataMakerSim->GetName()) ;
75 fAliITSQADataMakerSim->SetTitle((const char*)qadm.fAliITSQADataMakerSim->GetTitle());
78 //__________________________________________________________________
79 AliITSQASPDDataMakerSim& AliITSQASPDDataMakerSim::operator = (const AliITSQASPDDataMakerSim& qac )
82 this->~AliITSQASPDDataMakerSim();
83 new(this) AliITSQASPDDataMakerSim(qac);
87 //____________________________________________________________________________
88 void AliITSQASPDDataMakerSim::StartOfDetectorCycle()
90 //Detector specific actions at start of cycle
91 AliDebug(1,"AliITSQADM::Start of SPD Cycle\n");
94 //____________________________________________________________________________
95 void AliITSQASPDDataMakerSim::EndOfDetectorCycle(AliQA::TASKINDEX_t /*task*/, TObjArray* /*list*/)
97 // launch the QA checking
98 AliDebug(1,"AliITSDM instantiates checker with Run(AliQA::kITS, task, list)\n");
100 //AliQAChecker::Instance()->Run( AliQA::kITS , task, list);
103 //____________________________________________________________________________
104 void AliITSQASPDDataMakerSim::InitDigits()
106 // Initialization for DIGIT data - SPD -
108 fDigitsOffset = (fAliITSQADataMakerSim->fDigitsQAList)->GetEntries();
109 //fSPDhDigits must be incremented by one unit every time a histogram is ADDED to the QA List
114 TH1F *hlayer = new TH1F("LayPattern_SPD","Layer map - SPD",6,0.,6.);
115 hlayer->GetXaxis()->SetTitle("Layer number");
116 hlayer->GetYaxis()->SetTitle("Entries");
117 fAliITSQADataMakerSim->Add2DigitsList(hlayer,1+fDigitsOffset);
120 TH1F **hmod = new TH1F*[2];
121 for (Int_t iLay=0; iLay<2; iLay++) {
122 sprintf(name,"ModPattern_SPD%d",iLay+1);
123 sprintf(title,"Module map - SPD Layer %d",iLay+1);
124 hmod[iLay]=new TH1F(name,title,240,0,240);
125 hmod[iLay]->GetXaxis()->SetTitle("Module number");
126 hmod[iLay]->GetYaxis()->SetTitle("Entries");
127 fAliITSQADataMakerSim->Add2DigitsList(hmod[iLay],2+iLay+fDigitsOffset);
131 TH1F *hcolumns = new TH1F("Columns_SPD","Columns - SPD",160,0.,160.);
132 hcolumns->GetXaxis()->SetTitle("Column number");
133 hcolumns->GetYaxis()->SetTitle("Entries");
134 fAliITSQADataMakerSim->Add2DigitsList(hcolumns,4+fDigitsOffset);
137 TH1F *hrows = new TH1F("Rows_SPD","Rows - SPD",256,0.,256.);
138 hrows->GetXaxis()->SetTitle("Row number");
139 hrows->GetYaxis()->SetTitle("Entries");
140 fAliITSQADataMakerSim->Add2DigitsList(hrows,5+fDigitsOffset);
143 TH1F** hMultSPDdigits = new TH1F*[2];
144 for (Int_t iLay=0; iLay<2; ++iLay) {
145 sprintf(name,"DigitMultiplicity_SPD%d",iLay+1);
146 sprintf(title,"Digit multiplicity - SPD Layer %d",iLay+1);
147 hMultSPDdigits[iLay]=new TH1F(name,title,200,0.,200.);
148 hMultSPDdigits[iLay]->GetXaxis()->SetTitle("Digit multiplicity");
149 hMultSPDdigits[iLay]->GetYaxis()->SetTitle("Entries");
150 fAliITSQADataMakerSim->Add2DigitsList(hMultSPDdigits[iLay], 6+iLay+fDigitsOffset);
154 TH2F *hMultSPDdig2MultSPDdig1 = new TH2F("DigitMultCorrelation_SPD","Digit multiplicity correlation - SPD",200,0.,200.,200,0.,200.);
155 hMultSPDdig2MultSPDdig1->GetXaxis()->SetTitle("Digit multiplicity (Layer 1)");
156 hMultSPDdig2MultSPDdig1->GetYaxis()->SetTitle("Digit multiplicity (Layer 2)");
157 fAliITSQADataMakerSim->Add2DigitsList(hMultSPDdig2MultSPDdig1,8+fDigitsOffset);
160 AliDebug(1,Form("%d SPD Digits histograms booked\n",fSPDhDigits));
164 //____________________________________________________________________________
165 void AliITSQASPDDataMakerSim::MakeDigits(TTree *digits)
167 // Fill QA for DIGIT - SPD -
168 AliITS *fITS = (AliITS*)gAlice->GetModule("ITS");
169 fITS->SetTreeAddress();
170 TClonesArray *iITSdigits = fITS->DigitsAddress(0); // 0->SPD
175 for (Int_t imod=0; imod<240; ++imod){
176 digits->GetEvent(imod);
177 Int_t ndigits = iITSdigits->GetEntries();
179 fAliITSQADataMakerSim->GetDigitsData(1+fDigitsOffset)->Fill(0.5,ndigits);
180 fAliITSQADataMakerSim->GetDigitsData(2+fDigitsOffset)->Fill(imod,ndigits);
184 fAliITSQADataMakerSim->GetDigitsData(1+fDigitsOffset)->Fill(1,ndigits);
185 fAliITSQADataMakerSim->GetDigitsData(3+fDigitsOffset)->Fill(imod,ndigits);
188 for (Int_t idig=0; idig<ndigits; ++idig) {
189 AliITSdigit *dig=(AliITSdigit*)iITSdigits->UncheckedAt(idig);
190 Int_t col=dig->GetCoord1(); // cell number z
191 Int_t row=dig->GetCoord2(); // cell number x
192 fAliITSQADataMakerSim->GetDigitsData(4+fDigitsOffset)->Fill(col);
193 fAliITSQADataMakerSim->GetDigitsData(5+fDigitsOffset)->Fill(row);
196 fAliITSQADataMakerSim->GetDigitsData(6+fDigitsOffset)->Fill(nDigitsL1);
197 fAliITSQADataMakerSim->GetDigitsData(7+fDigitsOffset)->Fill(nDigitsL2);
198 fAliITSQADataMakerSim->GetDigitsData(8+fDigitsOffset)->Fill(nDigitsL1,nDigitsL2);
201 //____________________________________________________________________________
202 void AliITSQASPDDataMakerSim::InitSDigits()
204 // Initialization for SDIGIT data - SPD -
205 fSDigitsOffset = (fAliITSQADataMakerSim->fSDigitsQAList)->GetEntries();
207 //fSPDhSDigits must be incremented by one unit every time a histogram is ADDED to the QA List
211 TH1F *hlayer = new TH1F("LayPattern_SPD","Layer map - SPD",6,0.,6.);
212 hlayer->GetXaxis()->SetTitle("Layer number");
213 hlayer->GetYaxis()->SetTitle("Entries");
214 fAliITSQADataMakerSim->Add2SDigitsList(hlayer,1+fSDigitsOffset);
217 TH1F **hmod = new TH1F*[2];
218 for (Int_t iLay=0; iLay<2; ++iLay) {
219 sprintf(name,"ModPattern_SPD%d",iLay+1);
220 sprintf(title,"Module map - SPD Layer %d",iLay+1);
221 hmod[iLay]=new TH1F(name,title,240,0,240);
222 hmod[iLay]->GetXaxis()->SetTitle("Module number");
223 hmod[iLay]->GetYaxis()->SetTitle("Entries");
224 fAliITSQADataMakerSim->Add2SDigitsList(hmod[iLay],2+iLay+fSDigitsOffset);
229 AliDebug(1,Form("%d SPD SDigits histograms booked\n",fSPDhSDigits));
233 //____________________________________________________________________________
234 void AliITSQASPDDataMakerSim::MakeSDigits(TTree *sdigits)
236 // Fill QA for SDIGIT - SPD -
237 TBranch *brchSDigits = sdigits->GetBranch("ITS");
238 for (Int_t imod=0; imod<240; ++imod){
239 TClonesArray * sdig = new TClonesArray( "AliITSpListItem",1000 );
240 brchSDigits->SetAddress( &sdig );
241 brchSDigits->GetEvent(imod);
242 Int_t nsdig=sdig->GetEntries();
244 fAliITSQADataMakerSim->GetSDigitsData(1+fSDigitsOffset)->Fill(0.5,nsdig);
245 fAliITSQADataMakerSim->GetSDigitsData(2+fSDigitsOffset)->Fill(imod,nsdig);
248 fAliITSQADataMakerSim->GetSDigitsData(1+fSDigitsOffset)->Fill(1,nsdig);
249 fAliITSQADataMakerSim->GetSDigitsData(3+fSDigitsOffset)->Fill(imod,nsdig);
256 //____________________________________________________________________________
257 void AliITSQASPDDataMakerSim::InitHits()
259 // Initialization for HITS data - SPD -
260 fHitsOffset = (fAliITSQADataMakerSim->fHitsQAList)->GetEntries();
262 //fSPDhHits must be incremented by one unit every time a histogram is ADDED to the QA List
266 TH1F *hlayer = new TH1F("LayPattern_SPD","Layer map - SPD",6,0.,6.);
267 hlayer->GetXaxis()->SetTitle("Layer number");
268 hlayer->GetYaxis()->SetTitle("Entries");
269 fAliITSQADataMakerSim->Add2HitsList(hlayer,1+fHitsOffset);
272 TH1F **hmod = new TH1F*[2];
273 for (Int_t iLay=0; iLay<2; ++iLay) {
274 sprintf(name,"ModPattern_SPD%d",iLay+1);
275 sprintf(title,"Module map - SPD Layer %d",iLay+1);
276 hmod[iLay]=new TH1F(name,title,240,0,240);
277 hmod[iLay]->GetXaxis()->SetTitle("Module number");
278 hmod[iLay]->GetYaxis()->SetTitle("Entries");
279 fAliITSQADataMakerSim->Add2HitsList(hmod[iLay],2+iLay+fHitsOffset);
283 TH1F *hhitlenght = new TH1F("Lenght_SPD","Hit lenght along y_{loc} coord",210,0.,210.);
284 hhitlenght->GetXaxis()->SetTitle("Hit lenght [#mum]");
285 hhitlenght->GetYaxis()->SetTitle("# hits");
286 fAliITSQADataMakerSim->Add2HitsList(hhitlenght,4+fHitsOffset);
289 TH1F *hEdepos = new TH1F("EnergyDeposit_SPD","Deposited energy distribution (y_{loc}>180 #mum)",150,0.,300.);
290 hEdepos->GetXaxis()->SetTitle("Deposited energy [keV]");
291 hEdepos->GetYaxis()->SetTitle("# hits");
292 fAliITSQADataMakerSim->Add2HitsList(hEdepos,5+fHitsOffset);
295 AliDebug(1,Form("%d SPD Hits histograms booked\n",fSPDhHits));
299 //____________________________________________________________________________
300 void AliITSQASPDDataMakerSim::MakeHits(TTree *hits)
302 // Fill QA for HITS - SPD -
303 AliITS *fITS = (AliITS*)gAlice->GetModule("ITS");
304 fITS->SetTreeAddress();
306 fITS->InitModules(-1,nmodules); //-1->number of modules taken from AliITSgeom class kept in fITSgeom
309 fITS->FillModules(hits,0);
311 for (Int_t imod=0; imod<240; ++imod){
312 AliITSmodule *module = fITS->GetModule(imod);
313 TObjArray *arrHits = module->GetHits();
314 Int_t nhits = arrHits->GetEntriesFast();
316 fAliITSQADataMakerSim->GetHitsData(1+fHitsOffset)->Fill(0.5,nhits);
317 fAliITSQADataMakerSim->GetHitsData(2+fHitsOffset)->Fill(imod,nhits);
319 fAliITSQADataMakerSim->GetHitsData(1+fHitsOffset)->Fill(1,nhits);
320 fAliITSQADataMakerSim->GetHitsData(3+fHitsOffset)->Fill(imod,nhits);
322 //printf("--w--AliITSQASPDDataMakerSim::MakeHits nhits = %d\n",nhits);
323 for (Int_t iHit=0; iHit<nhits; ++iHit) {
324 AliITShit *hit = (AliITShit*) arrHits->At(iHit);
325 Double_t xl,yl,zl,xl0,yl0,zl0;
327 hit->GetPositionL(xl,yl,zl,tof);
328 hit->GetPositionL0(xl0,yl0,zl0,tof0);
329 Float_t dyloc=TMath::Abs(yl-yl0)*10000.;
330 fAliITSQADataMakerSim->GetHitsData(4+fHitsOffset)->Fill(dyloc);
331 Float_t edep=hit->GetIonization()*1000000;
333 fAliITSQADataMakerSim->GetHitsData(5+fHitsOffset)->Fill(edep);