Update of the DA in order to cope with the dense filling scheme of 2011. Necessitates...
[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
095fab41 104 TH1I *fhDigTDC[64];
105 TH1I *fhDigADC[64];
508b9fc0 106
107 // create Digits histograms in Digits subdir
db72ff3b 108 TH1I * h0 = new TH1I("hDigitMultiplicity", "Digits multiplicity distribution in VZERO;# of Digits;Entries", 100, 0, 99) ;
508b9fc0 109 h0->Sumw2() ;
7d297381 110 Add2DigitsList(h0, 0, !expert, image) ;
508b9fc0 111
112 for (Int_t i=0; i<64; i++)
113 {
127e0d90 114 fhDigTDC[i] = new TH1I(Form("hDigitTDC%d", i),Form("Digit TDC in cell %d; TDC value;Entries",i),300,0.,149.);
115 fhDigADC[i]= new TH1I(Form("hDigitADC%d", i),Form("Digit ADC in cell %d;ADC value;Entries",i),1024,0.,1023.);
508b9fc0 116
7d297381 117 Add2DigitsList(fhDigTDC[i],i+1, !expert, image);
118 Add2DigitsList(fhDigADC[i],i+1+64, !expert, image);
508b9fc0 119 }
120}
121
122
123//____________________________________________________________________________
6252ceeb 124void AliVZEROQADataMakerSim::MakeHits()
508b9fc0 125{
126 //make QA data from Hits
127
6252ceeb 128 Int_t nhits = fHitsArray->GetEntriesFast();
129 GetHitsData(0)->Fill(nhits) ; // fills Hit multiplicity
508b9fc0 130 for (Int_t ihit=0;ihit<nhits;ihit++)
131 {
6252ceeb 132 AliVZEROhit * VZEROHit = (AliVZEROhit*) fHitsArray->UncheckedAt(ihit);
508b9fc0 133 if (!VZEROHit) {
134 AliError("The unchecked hit doesn't exist");
135 break;
136 }
095fab41 137 GetHitsData(1)->Fill(VZEROHit->Cell());
508b9fc0 138 }
139}
140
141
142//____________________________________________________________________________
143
144void AliVZEROQADataMakerSim::MakeHits(TTree *hitTree)
145{
146 //fills QA histos for Hits
6252ceeb 147 if (fHitsArray)
148 fHitsArray->Clear() ;
149 else
150 fHitsArray = new TClonesArray("AliVZEROhit", 1000);
508b9fc0 151
152 TBranch * branch = hitTree->GetBranch("VZERO") ;
153 if ( ! branch ) {
154 AliWarning("VZERO branch in Hit Tree not found") ;
155 } else {
156
157 if (branch) {
6252ceeb 158 branch->SetAddress(&fHitsArray);
508b9fc0 159 }else{
160 AliError("Branch VZERO hit not found");
161 exit(111);
162 }
eca4fa66 163 // Check id histograms already created for this Event Specie
164 if ( ! GetHitsData(0) )
165 InitHits() ;
166
508b9fc0 167 Int_t ntracks = (Int_t) hitTree->GetEntries();
168
169 if (ntracks<=0) return;
170 // Start loop on tracks in the hits containers
171 for (Int_t track=0; track<ntracks;track++) {
172 branch->GetEntry(track);
6252ceeb 173 Int_t nhits = fHitsArray->GetEntriesFast();
174 GetHitsData(0)->Fill(nhits) ; // fills Hit multiplicity
508b9fc0 175 for (Int_t ihit=0;ihit<nhits;ihit++)
176 {
6252ceeb 177 AliVZEROhit * VZEROHit = (AliVZEROhit*) fHitsArray->UncheckedAt(ihit);
508b9fc0 178 if (!VZEROHit) {
179 AliError("The unchecked hit doesn't exist");
180 break;
181 }
095fab41 182 GetHitsData(1)->Fill(VZEROHit->Cell());
508b9fc0 183 }
184 }
185 }
186}
187
188
189//____________________________________________________________________________
6252ceeb 190void AliVZEROQADataMakerSim::MakeDigits()
508b9fc0 191{
192 // makes data from Digits
193
6252ceeb 194 GetDigitsData(0)->Fill(fDigitsArray->GetEntriesFast()) ;
195 TIter next(fDigitsArray) ;
508b9fc0 196 AliVZEROdigit *VZERODigit ;
197 while ( (VZERODigit = dynamic_cast<AliVZEROdigit *>(next())) ) {
198 Int_t PMNumber = VZERODigit->PMNumber();
199 GetDigitsData(PMNumber +1)->Fill( VZERODigit->Time()) ; // in 100 of picoseconds
200 GetDigitsData(PMNumber +1+64)->Fill( VZERODigit->ADC()) ;
201 }
202}
203
204
205//____________________________________________________________________________
206void AliVZEROQADataMakerSim::MakeDigits(TTree *digitTree)
207{
208 // makes data from Digit Tree
209
6252ceeb 210 if (fDigitsArray)
211 fDigitsArray->Clear() ;
212 else
213 fDigitsArray = new TClonesArray("AliVZEROdigit", 1000) ;
508b9fc0 214
095fab41 215 TBranch * branch = digitTree->GetBranch("VZERODigit") ;
508b9fc0 216 if ( ! branch ) {
217 AliWarning("VZERO branch in Digit Tree not found") ;
218 } else {
6252ceeb 219 branch->SetAddress(&fDigitsArray) ;
508b9fc0 220 branch->GetEntry(0) ;
6252ceeb 221 MakeDigits() ;
508b9fc0 222 }
223}
224
225
226//____________________________________________________________________________
227void AliVZEROQADataMakerSim::StartOfDetectorCycle()
228{
229 //Detector specific actions at start of cycle
230
231}