Script to create a random bad channel map.
[u/mrichter/AliRoot.git] / VZERO / AliVZEROQADataMakerSim.cxx
CommitLineData
508b9fc0 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
17/* $Id: AliVZEROQADataMakerSim.cxx 23123 2007-12-18 09:08:18Z hristov $ */
18
19//---
20// Produces the data needed to calculate the quality assurance.
21// All data must be mergeable objects.
22// Author : BC
23//---
24
25// --- ROOT system ---
26#include <TClonesArray.h>
27#include <TFile.h>
28#include <TH1F.h>
29#include <TDirectory.h>
30// --- Standard library ---
31
32// --- AliRoot header files ---
33#include "AliESDEvent.h"
34#include "AliLog.h"
35#include "AliVZEROdigit.h"
36#include "AliVZEROhit.h"
37#include "AliVZEROQADataMakerSim.h"
38#include "AliQAChecker.h"
39
40ClassImp(AliVZEROQADataMakerSim)
41
42//____________________________________________________________________________
43 AliVZEROQADataMakerSim::AliVZEROQADataMakerSim() :
44 AliQADataMakerSim(AliQA::GetDetName(AliQA::kVZERO), "VZERO Quality Assurance Data Maker")
45
46{
47 // constructor
48
49
50}
51
52//____________________________________________________________________________
53AliVZEROQADataMakerSim::AliVZEROQADataMakerSim(const AliVZEROQADataMakerSim& qadm) :
54 AliQADataMakerSim()
55{
56 //copy constructor
57
58 SetName((const char*)qadm.GetName()) ;
59 SetTitle((const char*)qadm.GetTitle());
60}
61
62//__________________________________________________________________
63AliVZEROQADataMakerSim& AliVZEROQADataMakerSim::operator = (const AliVZEROQADataMakerSim& qadm )
64{
65 // Assign operator.
66 this->~AliVZEROQADataMakerSim();
67 new(this) AliVZEROQADataMakerSim(qadm);
68 return *this;
69}
70//____________________________________________________________________________
71void AliVZEROQADataMakerSim::EndOfDetectorCycle(AliQA::TASKINDEX_t task, TObjArray * list)
72{
73 //Detector specific actions at end of cycle
74 // do the QA checking
75 AliQAChecker::Instance()->Run(AliQA::kVZERO, task, list) ;
76}
77
78
79//____________________________________________________________________________
80void AliVZEROQADataMakerSim::InitHits()
81{
508b9fc0 82
83 // create Hits histograms in Hits subdir
095fab41 84
85 TH1I * h0 = new TH1I("hHitMultiplicity", "Hit multiplicity distribution in VZERO", 300, 0, 299) ;
508b9fc0 86 h0->Sumw2() ;
87 Add2HitsList(h0, 0) ;
095fab41 88
89 TH1I * h1 = new TH1I("hHitCellNumber", "Hit cell distribution in VZERO", 80, 0, 79) ;
90 h1->Sumw2() ;
91 Add2HitsList(h1, 1) ;
508b9fc0 92
508b9fc0 93}
94
95//____________________________________________________________________________
96void AliVZEROQADataMakerSim::InitDigits()
97{
98 // create Digits histograms in Digits subdir
99
100
101 char TDCname[10];
102 char ADCname[12];
095fab41 103 TH1I *fhDigTDC[64];
104 TH1I *fhDigADC[64];
508b9fc0 105 char texte[30];
106
107 // create Digits histograms in Digits subdir
095fab41 108 TH1I * h0 = new TH1I("hDigitMultiplicity", "Digits multiplicity distribution in VZERO", 100, 0, 99) ;
508b9fc0 109 h0->Sumw2() ;
110 Add2DigitsList(h0, 0) ;
111
112 for (Int_t i=0; i<64; i++)
113 {
114 sprintf(TDCname, "hDigitTDC%d", i);
115 sprintf(texte,"Digit TDC in cell %d",i);
095fab41 116 fhDigTDC[i] = new TH1I(TDCname,texte,300,0.,149.);
508b9fc0 117
118 sprintf(ADCname,"hDigitADC%d",i);
119 sprintf(texte,"Digit ADC in cell %d",i);
095fab41 120 fhDigADC[i]= new TH1I(ADCname,texte,1024,0.,1023.);
508b9fc0 121
122 Add2DigitsList(fhDigTDC[i],i+1);
123 Add2DigitsList(fhDigADC[i],i+1+64);
124 }
125}
126
127
128//____________________________________________________________________________
129void AliVZEROQADataMakerSim::MakeHits(TClonesArray * hits)
130{
131 //make QA data from Hits
132
133 GetHitsData(0)->Fill(hits->GetEntriesFast()) ; // fills Hit multiplicity
134 Int_t nhits = hits->GetEntriesFast();
135 for (Int_t ihit=0;ihit<nhits;ihit++)
136 {
137 AliVZEROhit * VZEROHit = (AliVZEROhit*) hits->UncheckedAt(ihit);
138 if (!VZEROHit) {
139 AliError("The unchecked hit doesn't exist");
140 break;
141 }
095fab41 142 GetHitsData(1)->Fill(VZEROHit->Cell());
508b9fc0 143 }
144}
145
146
147//____________________________________________________________________________
148
149void AliVZEROQADataMakerSim::MakeHits(TTree *hitTree)
150{
151 //fills QA histos for Hits
152 TClonesArray * hits = new TClonesArray("AliVZEROhit", 1000);
153
154 TBranch * branch = hitTree->GetBranch("VZERO") ;
155 if ( ! branch ) {
156 AliWarning("VZERO branch in Hit Tree not found") ;
157 } else {
158
159 if (branch) {
160 branch->SetAddress(&hits);
161 }else{
162 AliError("Branch VZERO hit not found");
163 exit(111);
164 }
165 Int_t ntracks = (Int_t) hitTree->GetEntries();
166
167 if (ntracks<=0) return;
168 // Start loop on tracks in the hits containers
169 for (Int_t track=0; track<ntracks;track++) {
170 branch->GetEntry(track);
171 GetHitsData(0)->Fill(hits->GetEntriesFast()) ; // fills Hit multiplicity
172 Int_t nhits = hits->GetEntriesFast();
173 for (Int_t ihit=0;ihit<nhits;ihit++)
174 {
175 AliVZEROhit * VZEROHit = (AliVZEROhit*) hits->UncheckedAt(ihit);
176 if (!VZEROHit) {
177 AliError("The unchecked hit doesn't exist");
178 break;
179 }
095fab41 180 GetHitsData(1)->Fill(VZEROHit->Cell());
508b9fc0 181 }
182 }
183 }
184}
185
186
187//____________________________________________________________________________
188void AliVZEROQADataMakerSim::MakeDigits(TClonesArray * digits)
189{
190 // makes data from Digits
191
192 GetDigitsData(0)->Fill(digits->GetEntriesFast()) ;
193 TIter next(digits) ;
194 AliVZEROdigit *VZERODigit ;
195 while ( (VZERODigit = dynamic_cast<AliVZEROdigit *>(next())) ) {
196 Int_t PMNumber = VZERODigit->PMNumber();
197 GetDigitsData(PMNumber +1)->Fill( VZERODigit->Time()) ; // in 100 of picoseconds
198 GetDigitsData(PMNumber +1+64)->Fill( VZERODigit->ADC()) ;
199 }
200}
201
202
203//____________________________________________________________________________
204void AliVZEROQADataMakerSim::MakeDigits(TTree *digitTree)
205{
206 // makes data from Digit Tree
207
095fab41 208 TClonesArray * digits = new TClonesArray("AliVZEROdigit", 1000) ;
508b9fc0 209
095fab41 210 TBranch * branch = digitTree->GetBranch("VZERODigit") ;
508b9fc0 211 if ( ! branch ) {
212 AliWarning("VZERO branch in Digit Tree not found") ;
213 } else {
214 branch->SetAddress(&digits) ;
215 branch->GetEntry(0) ;
216 MakeDigits(digits) ;
217 }
218}
219
220
221//____________________________________________________________________________
222void AliVZEROQADataMakerSim::StartOfDetectorCycle()
223{
224 //Detector specific actions at start of cycle
225
226}