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
26 // --- ROOT system ---
28 // --- Standard library ---
30 // --- AliRoot header files ---
31 #include "AliITSQASDDDataMakerSim.h"
34 #include "AliQAChecker.h"
35 #include "AliQADataMakerSim.h"
36 #include "AliITSQADataMakerSim.h"
37 #include "AliRawReader.h"
38 #include "AliITSdigit.h"
40 #include "AliITSmodule.h"
41 #include "AliITShit.h"
42 #include "AliITSLoader.h"
43 #include "AliRunLoader.h"
45 #include "AliITSsegmentationSDD.h"
46 #include "AliITSpList.h"
48 ClassImp(AliITSQASDDDataMakerSim)
50 //____________________________________________________________________________
51 AliITSQASDDDataMakerSim::AliITSQASDDDataMakerSim(AliITSQADataMakerSim *aliITSQADataMakerSim) :
53 fAliITSQADataMakerSim(aliITSQADataMakerSim),
61 //ctor used to discriminate OnLine-Offline analysis
62 fGenOffsetH= new Int_t[AliRecoParam::kNSpecies];
63 fGenOffsetS= new Int_t[AliRecoParam::kNSpecies];
64 fGenOffsetD= new Int_t[AliRecoParam::kNSpecies];
65 for(Int_t i=0; i<AliRecoParam::kNSpecies; i++)
73 //____________________________________________________________________________
74 AliITSQASDDDataMakerSim::AliITSQASDDDataMakerSim(const AliITSQASDDDataMakerSim& qadm) :
76 fAliITSQADataMakerSim(qadm.fAliITSQADataMakerSim),
77 fSDDhHTask(qadm.fSDDhHTask),
78 fSDDhSTask(qadm.fSDDhSTask),
79 fSDDhDTask(qadm.fSDDhDTask),
80 fGenOffsetH(qadm.fGenOffsetH),
81 fGenOffsetS(qadm.fGenOffsetS),
82 fGenOffsetD(qadm.fGenOffsetD)
85 fAliITSQADataMakerSim->SetName((const char*)qadm.fAliITSQADataMakerSim->GetName()) ;
86 fAliITSQADataMakerSim->SetTitle((const char*)qadm.fAliITSQADataMakerSim->GetTitle());
90 //__________________________________________________________________
91 AliITSQASDDDataMakerSim& AliITSQASDDDataMakerSim::operator = (const AliITSQASDDDataMakerSim& qac )
94 this->~AliITSQASDDDataMakerSim();
95 new(this) AliITSQASDDDataMakerSim(qac);
99 //____________________________________________________________________________
100 void AliITSQASDDDataMakerSim::StartOfDetectorCycle()
102 //Detector specific actions at start of cycle
103 AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM::Start of SDD Cycle\n");
106 //____________________________________________________________________________
107 void AliITSQASDDDataMakerSim::EndOfDetectorCycle(AliQAv1::TASKINDEX_t /*task*/, TObjArray** /*list*/)
109 // launch the QA checking
111 AliDebug(AliQAv1::GetQADebugLevel(),"AliITSDM instantiates checker with Run(AliQAv1::kITS, task, list)\n");
112 //AliQAChecker::Instance()->Run( AliQAv1::kITS , task, list);
115 //____________________________________________________________________________
116 Int_t AliITSQASDDDataMakerSim::InitDigits()
118 // Initialization for DIGIT data - SDD -
119 const Bool_t expert = kTRUE ;
120 const Bool_t image = kTRUE ;
122 //fGenOffsetD = (fAliITSQADataMakerSim->fDigitsQAList[AliRecoParam::kDefault])->GetEntries();
123 //fSDDhTask must be incremented by one unit every time a histogram is ADDED to the QA List
124 TH1F* h0=new TH1F("SDD DIGITS Module Pattern","SDD DIGITS Module Pattern",260,239.5,499.5); //hmod
125 h0->GetXaxis()->SetTitle("SDD Module Number");
126 h0->GetYaxis()->SetTitle("# DIGITS");
127 rv = fAliITSQADataMakerSim->Add2DigitsList(h0,fGenOffsetD[fAliITSQADataMakerSim->GetEventSpecie()], !expert, image);
129 TH1F* h1=new TH1F("SDD Anode Distribution","SDD DIGITS Anode Distribution",512,-0.5,511.5); //hanocc
130 h1->GetXaxis()->SetTitle("Anode Number");
131 h1->GetYaxis()->SetTitle("# DIGITS");
132 rv = fAliITSQADataMakerSim->Add2DigitsList(h1,1+fGenOffsetD[fAliITSQADataMakerSim->GetEventSpecie()], !expert, image);
134 TH1F* h2=new TH1F("SDD Tbin Distribution","SDD DIGITS Tbin Distribution",256,-0.5,255.5); //htbocc
135 h2->GetXaxis()->SetTitle("Tbin Number");
136 h2->GetYaxis()->SetTitle("# DIGITS");
137 rv = fAliITSQADataMakerSim->Add2DigitsList(h2,2+fGenOffsetD[fAliITSQADataMakerSim->GetEventSpecie()], !expert, image);
139 TH1F* h3=new TH1F("SDD ADC Counts Distribution","SDD DIGITS ADC Counts Distribution",200,0.,1024.); //hsig
140 h3->GetXaxis()->SetTitle("ADC Value");
141 h3->GetYaxis()->SetTitle("# DIGITS");
142 rv = fAliITSQADataMakerSim->Add2DigitsList(h3,3+fGenOffsetD[fAliITSQADataMakerSim->GetEventSpecie()], !expert, image);
144 AliDebug(AliQAv1::GetQADebugLevel(),Form("%d SDD Digits histograms booked\n",fSDDhDTask));
148 //____________________________________________________________________________
149 Int_t AliITSQASDDDataMakerSim::MakeDigits(TTree * digits)
152 // Fill QA for DIGIT - SDD -
155 AliITS *fITS = (AliITS*)gAlice->GetModule("ITS");
156 fITS->SetTreeAddress();
157 TClonesArray *iITSdigits = fITS->DigitsAddress(1);
158 for(Int_t i=0; i<260; i++){
160 digits->GetEvent(nmod);
161 Int_t ndigits = iITSdigits->GetEntries();
162 fAliITSQADataMakerSim->FillDigitsData(fGenOffsetD[fAliITSQADataMakerSim->GetEventSpecie()],nmod,ndigits);
163 for (Int_t idig=0; idig<ndigits; idig++) {
164 AliITSdigit *dig=(AliITSdigit*)iITSdigits->UncheckedAt(idig);
165 Int_t iz=dig->GetCoord1(); // cell number z
166 Int_t ix=dig->GetCoord2(); // cell number x
167 Int_t sig=dig->GetSignal();
168 fAliITSQADataMakerSim->FillDigitsData(1+fGenOffsetD[fAliITSQADataMakerSim->GetEventSpecie()],iz);
169 fAliITSQADataMakerSim->FillDigitsData(2+fGenOffsetD[fAliITSQADataMakerSim->GetEventSpecie()],ix);
170 fAliITSQADataMakerSim->FillDigitsData(3+fGenOffsetD[fAliITSQADataMakerSim->GetEventSpecie()],sig);
177 //____________________________________________________________________________
178 Int_t AliITSQASDDDataMakerSim::InitSDigits()
180 // Initialization for SDIGIT data - SDD -
181 const Bool_t expert = kTRUE ;
182 const Bool_t image = kTRUE ;
184 //fGenOffsetS = (fAliITSQADataMakerSim->fSDigitsQAList[AliRecoParam::kDefault])->GetEntries();
185 //fSDDhTask must be incremented by one unit every time a histogram is ADDED to the QA List
186 TH1F* h0=new TH1F("SDD SDIGITS Module Pattern","SDIGITS SDD Module Pattern",260,239.5,499.5); //hmod
187 h0->GetXaxis()->SetTitle("SDD Module Number");
188 h0->GetYaxis()->SetTitle("# SDIGITS");
189 rv = fAliITSQADataMakerSim->Add2SDigitsList(h0,fGenOffsetS[fAliITSQADataMakerSim->GetEventSpecie()], !expert, image);
191 TH1F* h1=new TH1F("SDD Anode Distribution","SDIGITS SDD Anode Distribution",512,-0.5,511.5); //hanocc
192 h1->GetXaxis()->SetTitle("Anode Number");
193 h1->GetYaxis()->SetTitle("# SDIGITS");
194 rv = fAliITSQADataMakerSim->Add2SDigitsList(h1,1+fGenOffsetS[fAliITSQADataMakerSim->GetEventSpecie()], !expert, image);
196 TH1F* h2=new TH1F("SDD Tbin Distribution","SDIGITS SDD Tbin Distribution",256,-0.5,255.5); //htbocc
197 h2->GetXaxis()->SetTitle("Tbin Number");
198 h2->GetYaxis()->SetTitle("# SDIGITS");
199 rv = fAliITSQADataMakerSim->Add2SDigitsList(h2,2+fGenOffsetS[fAliITSQADataMakerSim->GetEventSpecie()]);
201 TH1F* h3=new TH1F("SDD ADC Counts Distribution","SDIGITS SDD ADC Counts Distribution",200,0.,1024.); //hsig
202 h3->GetXaxis()->SetTitle("ADC Value");
203 h3->GetYaxis()->SetTitle("# SDIGITS");
204 rv = fAliITSQADataMakerSim->Add2SDigitsList(h3,3+fGenOffsetS[fAliITSQADataMakerSim->GetEventSpecie()], !expert, image);
207 AliDebug(AliQAv1::GetQADebugLevel(),Form("%d SDD SDigits histograms booked\n",fSDDhSTask));
211 //____________________________________________________________________________
212 Int_t AliITSQASDDDataMakerSim::MakeSDigits(TTree * sdigits)
214 // Fill QA for SDIGIT - SDD -
217 AliITSsegmentationSDD* seg = new AliITSsegmentationSDD();
218 Int_t nan=seg->Npz();
219 Int_t ntb=seg->Npx();
221 AliITSpList* list=new AliITSpList(nan,ntb*scaleSize);
223 //AliITS *fITS = (AliITS*)gAlice->GetModule("ITS");
224 //fITS->SetTreeAddress();
225 //TClonesArray *ITSdigits = fITS->DigitsAddress(1);
226 //TFile *sper = new TFile("sper.root","CREATE"); //agginto a mano x prova
231 TBranch *brchSDigits = sdigits->GetBranch("ITS");
232 static TClonesArray * sdig ;
234 sdig = new TClonesArray( "AliITSpListItem",1000 );
235 for(Int_t id=0; id<260; id++){
237 brchSDigits->SetAddress( &sdig );
238 brchSDigits->GetEvent(nmod);
239 Int_t nsdig=sdig->GetEntries();
240 fAliITSQADataMakerSim->FillSDigitsData(fGenOffsetS[fAliITSQADataMakerSim->GetEventSpecie()],nmod,nsdig);
241 for(Int_t i=0;i<nsdig;i++){
242 AliITSpListItem *cell=(AliITSpListItem*)sdig->At(i);
243 Float_t sig=cell->GetSignal();
244 Int_t idx=cell->GetIndex();
246 list->GetCell(idx,ia,it);
247 fAliITSQADataMakerSim->FillSDigitsData(1+fGenOffsetS[fAliITSQADataMakerSim->GetEventSpecie()],ia);
248 fAliITSQADataMakerSim->FillSDigitsData(2+fGenOffsetS[fAliITSQADataMakerSim->GetEventSpecie()],it);
249 fAliITSQADataMakerSim->FillSDigitsData(3+fGenOffsetS[fAliITSQADataMakerSim->GetEventSpecie()],sig);
257 //____________________________________________________________________________
258 Int_t AliITSQASDDDataMakerSim::InitHits()
261 // Initialization for HITS data - SDD -
262 const Bool_t expert = kTRUE ;
263 const Bool_t image = kTRUE ;
266 //fGenOffsetH = (fAliITSQADataMakerSim->fHitsQAList[AliRecoParam::kDefault])->GetEntries();
267 //fSDDhTask must be incremented by one unit every time a histogram is ADDED to the QA List
268 //printf("AliITSQASDDDataMakerSim::InitHits called \n");
269 TH1F *h0=new TH1F("SDD HITS Module Pattern","SDD HITS Module Pattern",260,239.5,499.5);
270 h0->GetXaxis()->SetTitle("SDD Module Number");
271 h0->GetYaxis()->SetTitle("# HITS");
272 rv = fAliITSQADataMakerSim->Add2HitsList(h0,fGenOffsetH[fAliITSQADataMakerSim->GetEventSpecie()], !expert, image);
274 TH1F *h1=new TH1F("SDD HIT lenght along local Y Coord","SDD HIT lenght along local Y Coord",200,0.,350.);
275 h1->GetXaxis()->SetTitle("HIT lenght (um)");
276 h1->GetYaxis()->SetTitle("# HITS");
277 rv = fAliITSQADataMakerSim->Add2HitsList(h1,1+fGenOffsetH[fAliITSQADataMakerSim->GetEventSpecie()], !expert, image);
279 TH1F *h2=new TH1F("SDD HIT lenght along local Y Coord - Zoom","SDD HIT lenght along local Y Coord - Zoom",200,250.,350.);
280 h2->GetXaxis()->SetTitle("HIT lenght (um)");
281 h2->GetYaxis()->SetTitle("# HITS");
282 rv = fAliITSQADataMakerSim->Add2HitsList(h2,2+fGenOffsetH[fAliITSQADataMakerSim->GetEventSpecie()], !expert, image);
284 TH1F *h3=new TH1F("SDD Deposited Energy Distribution (loc Y > 200um)","SDD HITS Deposited Energy Distribution (loc Y > 200um)",200,0.,350.);
285 h3->GetXaxis()->SetTitle("ADC counts ");
286 h3->GetYaxis()->SetTitle("# HITS");
287 rv = fAliITSQADataMakerSim->Add2HitsList(h3,3+fGenOffsetH[fAliITSQADataMakerSim->GetEventSpecie()], !expert, image);
289 AliDebug(AliQAv1::GetQADebugLevel(),Form("%d SDD Hits histograms booked\n",fSDDhHTask));
293 //____________________________________________________________________________
294 Int_t AliITSQASDDDataMakerSim::MakeHits(TTree * hits)
296 // Fill QA for HITS - SDD -
299 AliITS *fITS = (AliITS*)gAlice->GetModule("ITS");
300 fITS->SetTreeAddress();
302 if(!(fITS->InitModules(-1,nmodules))){
303 AliError("ITS geometry not available - nothing done");
307 fITS->FillModules(hits,0);
309 for(Int_t i=0; i<260; i++){
311 AliITSmodule *modu = fITS->GetModule(nmod);
312 TObjArray *arrHits = modu->GetHits();
313 Int_t nhits = arrHits->GetEntriesFast();
314 ////printf("--w--AliITSQASDDDataMakerSim::MakeHits nhits = %d\n",nhits);
315 for (Int_t iHit=0;iHit<nhits;iHit++) {
316 AliITShit *hit = (AliITShit*) arrHits->At(iHit);
317 fAliITSQADataMakerSim->FillHitsData(fGenOffsetH[fAliITSQADataMakerSim->GetEventSpecie()],nmod);
318 Double_t xl,yl,zl,xl0,yl0,zl0;
320 hit->GetPositionL(xl,yl,zl,tof);
321 hit->GetPositionL0(xl0,yl0,zl0,tof0);
322 Float_t dyloc=TMath::Abs(yl-yl0)*10000.;
323 fAliITSQADataMakerSim->FillHitsData(1+fGenOffsetH[fAliITSQADataMakerSim->GetEventSpecie()],dyloc);
324 Float_t edep=hit->GetIonization()*1000000;
326 fAliITSQADataMakerSim->FillHitsData(2+fGenOffsetH[fAliITSQADataMakerSim->GetEventSpecie()],edep);
327 fAliITSQADataMakerSim->FillHitsData(3+fGenOffsetH[fAliITSQADataMakerSim->GetEventSpecie()],dyloc);
335 //_______________________________________________________________
337 Int_t AliITSQASDDDataMakerSim::GetOffset(AliQAv1::TASKINDEX_t task, Int_t specie){
338 // Returns histogram offset according to the specified task
340 if( task == AliQAv1::kHITS){
341 offset=fGenOffsetH[specie];
343 else if( task == AliQAv1::kSDIGITS) {
344 offset=fGenOffsetS[specie];
346 else if( task == AliQAv1::kDIGITS) {
347 offset=fGenOffsetD[specie];
350 AliInfo("No task has been selected. TaskHisto set to zero.\n");
356 //____________________________________________________________________________
357 void AliITSQASDDDataMakerSim::SetOffset(AliQAv1::TASKINDEX_t task, Int_t offset,Int_t specie ){
358 // Returns histogram offset according to the specified task
359 if( task == AliQAv1::kHITS){
360 fGenOffsetH[specie] = offset;
362 else if( task == AliQAv1::kSDIGITS) {
363 fGenOffsetS[specie] = offset;
365 else if( task == AliQAv1::kDIGITS) {
366 fGenOffsetD[specie] = offset;
369 AliInfo("No task has been selected. TaskHisto set to zero.\n");
373 //_______________________________________________________________
375 Int_t AliITSQASDDDataMakerSim::GetTaskHisto(AliQAv1::TASKINDEX_t task) {
376 // Returns the number of booked histograms for the selected task
378 if( task == AliQAv1::kHITS) {
379 histotot=fSDDhHTask ;
381 else if( task == AliQAv1::kSDIGITS) {
384 else if( task == AliQAv1::kDIGITS) {
385 histotot=fSDDhDTask ;
388 AliInfo("No task has been selected. TaskHisto set to zero.\n");