Macros used to produce the OCDB entries for MC. One for default (p-p) and one for...
[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() :
4e25ac79 44 AliQADataMakerSim(AliQAv1::GetDetName(AliQAv1::kVZERO), "VZERO Quality Assurance Data Maker")
508b9fc0 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//____________________________________________________________________________
4e25ac79 71void AliVZEROQADataMakerSim::EndOfDetectorCycle(AliQAv1::TASKINDEX_t task, TObjArray ** list)
508b9fc0 72{
73 //Detector specific actions at end of cycle
74 // do the QA checking
4e25ac79 75 AliQAChecker::Instance()->Run(AliQAv1::kVZERO, task, list) ;
508b9fc0 76}
77
78
79//____________________________________________________________________________
80void AliVZEROQADataMakerSim::InitHits()
81{
508b9fc0 82
83 // create Hits histograms in Hits subdir
7d297381 84 const Bool_t expert = kTRUE ;
85 const Bool_t image = kTRUE ;
095fab41 86
db72ff3b 87 TH1I * h0 = new TH1I("hHitMultiplicity", "Hit multiplicity distribution in VZERO;# of Hits;Entries", 300, 0, 299) ;
508b9fc0 88 h0->Sumw2() ;
7d297381 89 Add2HitsList(h0, 0, !expert, image) ;
095fab41 90
db72ff3b 91 TH1I * h1 = new TH1I("hHitCellNumber", "Hit cell distribution in VZERO;# of Hits;Entries", 80, 0, 79) ;
095fab41 92 h1->Sumw2() ;
7d297381 93 Add2HitsList(h1, 1, !expert, image) ;
508b9fc0 94
508b9fc0 95}
96
97//____________________________________________________________________________
98void AliVZEROQADataMakerSim::InitDigits()
99{
100 // create Digits histograms in Digits subdir
7d297381 101 const Bool_t expert = kTRUE ;
102 const Bool_t image = kTRUE ;
508b9fc0 103
a3f8715e 104 char TDCname[100];
105 char ADCname[100];
095fab41 106 TH1I *fhDigTDC[64];
107 TH1I *fhDigADC[64];
a3f8715e 108 char texte[100];
508b9fc0 109
110 // create Digits histograms in Digits subdir
db72ff3b 111 TH1I * h0 = new TH1I("hDigitMultiplicity", "Digits multiplicity distribution in VZERO;# of Digits;Entries", 100, 0, 99) ;
508b9fc0 112 h0->Sumw2() ;
7d297381 113 Add2DigitsList(h0, 0, !expert, image) ;
508b9fc0 114
115 for (Int_t i=0; i<64; i++)
116 {
117 sprintf(TDCname, "hDigitTDC%d", i);
db72ff3b 118 sprintf(texte,"Digit TDC in cell %d; TDC value;Entries",i);
095fab41 119 fhDigTDC[i] = new TH1I(TDCname,texte,300,0.,149.);
508b9fc0 120
121 sprintf(ADCname,"hDigitADC%d",i);
db72ff3b 122 sprintf(texte,"Digit ADC in cell %d;ADC value;Entries",i);
095fab41 123 fhDigADC[i]= new TH1I(ADCname,texte,1024,0.,1023.);
508b9fc0 124
7d297381 125 Add2DigitsList(fhDigTDC[i],i+1, !expert, image);
126 Add2DigitsList(fhDigADC[i],i+1+64, !expert, image);
508b9fc0 127 }
128}
129
130
131//____________________________________________________________________________
6252ceeb 132void AliVZEROQADataMakerSim::MakeHits()
508b9fc0 133{
134 //make QA data from Hits
135
6252ceeb 136 Int_t nhits = fHitsArray->GetEntriesFast();
137 GetHitsData(0)->Fill(nhits) ; // fills Hit multiplicity
508b9fc0 138 for (Int_t ihit=0;ihit<nhits;ihit++)
139 {
6252ceeb 140 AliVZEROhit * VZEROHit = (AliVZEROhit*) fHitsArray->UncheckedAt(ihit);
508b9fc0 141 if (!VZEROHit) {
142 AliError("The unchecked hit doesn't exist");
143 break;
144 }
095fab41 145 GetHitsData(1)->Fill(VZEROHit->Cell());
508b9fc0 146 }
147}
148
149
150//____________________________________________________________________________
151
152void AliVZEROQADataMakerSim::MakeHits(TTree *hitTree)
153{
154 //fills QA histos for Hits
6252ceeb 155 if (fHitsArray)
156 fHitsArray->Clear() ;
157 else
158 fHitsArray = new TClonesArray("AliVZEROhit", 1000);
508b9fc0 159
160 TBranch * branch = hitTree->GetBranch("VZERO") ;
161 if ( ! branch ) {
162 AliWarning("VZERO branch in Hit Tree not found") ;
163 } else {
164
165 if (branch) {
6252ceeb 166 branch->SetAddress(&fHitsArray);
508b9fc0 167 }else{
168 AliError("Branch VZERO hit not found");
169 exit(111);
170 }
eca4fa66 171 // Check id histograms already created for this Event Specie
172 if ( ! GetHitsData(0) )
173 InitHits() ;
174
508b9fc0 175 Int_t ntracks = (Int_t) hitTree->GetEntries();
176
177 if (ntracks<=0) return;
178 // Start loop on tracks in the hits containers
179 for (Int_t track=0; track<ntracks;track++) {
180 branch->GetEntry(track);
6252ceeb 181 Int_t nhits = fHitsArray->GetEntriesFast();
182 GetHitsData(0)->Fill(nhits) ; // fills Hit multiplicity
508b9fc0 183 for (Int_t ihit=0;ihit<nhits;ihit++)
184 {
6252ceeb 185 AliVZEROhit * VZEROHit = (AliVZEROhit*) fHitsArray->UncheckedAt(ihit);
508b9fc0 186 if (!VZEROHit) {
187 AliError("The unchecked hit doesn't exist");
188 break;
189 }
095fab41 190 GetHitsData(1)->Fill(VZEROHit->Cell());
508b9fc0 191 }
192 }
193 }
194}
195
196
197//____________________________________________________________________________
6252ceeb 198void AliVZEROQADataMakerSim::MakeDigits()
508b9fc0 199{
200 // makes data from Digits
201
6252ceeb 202 GetDigitsData(0)->Fill(fDigitsArray->GetEntriesFast()) ;
203 TIter next(fDigitsArray) ;
508b9fc0 204 AliVZEROdigit *VZERODigit ;
205 while ( (VZERODigit = dynamic_cast<AliVZEROdigit *>(next())) ) {
206 Int_t PMNumber = VZERODigit->PMNumber();
207 GetDigitsData(PMNumber +1)->Fill( VZERODigit->Time()) ; // in 100 of picoseconds
208 GetDigitsData(PMNumber +1+64)->Fill( VZERODigit->ADC()) ;
209 }
210}
211
212
213//____________________________________________________________________________
214void AliVZEROQADataMakerSim::MakeDigits(TTree *digitTree)
215{
216 // makes data from Digit Tree
217
6252ceeb 218 if (fDigitsArray)
219 fDigitsArray->Clear() ;
220 else
221 fDigitsArray = new TClonesArray("AliVZEROdigit", 1000) ;
508b9fc0 222
095fab41 223 TBranch * branch = digitTree->GetBranch("VZERODigit") ;
508b9fc0 224 if ( ! branch ) {
225 AliWarning("VZERO branch in Digit Tree not found") ;
226 } else {
6252ceeb 227 branch->SetAddress(&fDigitsArray) ;
508b9fc0 228 branch->GetEntry(0) ;
6252ceeb 229 MakeDigits() ;
508b9fc0 230 }
231}
232
233
234//____________________________________________________________________________
235void AliVZEROQADataMakerSim::StartOfDetectorCycle()
236{
237 //Detector specific actions at start of cycle
238
239}