]>
Commit | Line | Data |
---|---|---|
a64b872d | 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: AliCorrQADataMakerRec.cxx 27570 2008-07-24 21:49:27Z cvetan $ */ | |
18 | ||
19 | /* | |
20 | Produces the data needed to calculate the quality assurance. | |
21 | All data must be mergeable objects. | |
22 | Y. Schutz CERN July 2007 | |
23 | */ | |
24 | ||
25 | // --- ROOT system --- | |
26 | #include <TClonesArray.h> | |
27 | #include <TFile.h> | |
28 | #include <TH1F.h> | |
29 | #include <TH1I.h> | |
30 | #include <TH2F.h> | |
31 | #include <TNtupleD.h> | |
32 | #include <TParameter.h> | |
57acd2d2 | 33 | #include <TMath.h> |
a64b872d | 34 | |
35 | // --- Standard library --- | |
36 | ||
37 | // --- AliRoot header files --- | |
38 | #include "AliLog.h" | |
39 | #include "AliCorrQADataMakerRec.h" | |
40 | #include "AliQAChecker.h" | |
41 | ||
42 | ClassImp(AliCorrQADataMakerRec) | |
43 | ||
44 | //____________________________________________________________________________ | |
45 | AliCorrQADataMakerRec::AliCorrQADataMakerRec(AliQADataMaker ** qadm ) : | |
4e25ac79 | 46 | AliQADataMakerRec(AliQAv1::GetDetName(AliQAv1::kCORR), "Corr Quality Assurance Data Maker"), |
a64b872d | 47 | fMaxRawVar(0), |
097485c2 | 48 | fqadm(qadm), |
49 | fVarvalue(NULL) | |
a64b872d | 50 | { |
51 | // ctor | |
57acd2d2 | 52 | fCorrNt = new TNtupleD *[AliRecoParam::kNSpecies] ; |
53 | for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) | |
54 | fCorrNt[specie] = NULL ; | |
a64b872d | 55 | } |
56 | ||
57 | //____________________________________________________________________________ | |
58 | AliCorrQADataMakerRec::AliCorrQADataMakerRec(const AliCorrQADataMakerRec& qadm) : | |
59 | AliQADataMakerRec(), | |
60 | fMaxRawVar(qadm.fMaxRawVar), | |
097485c2 | 61 | fqadm(qadm.fqadm), |
62 | fVarvalue(NULL) | |
a64b872d | 63 | { |
64 | //copy ctor | |
65 | SetName((const char*)qadm.GetName()) ; | |
66 | SetTitle((const char*)qadm.GetTitle()); | |
097485c2 | 67 | if ( fMaxRawVar > 0 ) |
68 | fVarvalue = new Double_t[fMaxRawVar] ; | |
4bb9542a | 69 | |
70 | fCorrNt = new TNtupleD *[AliRecoParam::kNSpecies] ; | |
71 | for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) | |
72 | fCorrNt[specie] = qadm.fCorrNt[specie] ; | |
73 | ||
a64b872d | 74 | } |
75 | ||
76 | //__________________________________________________________________ | |
77 | AliCorrQADataMakerRec& AliCorrQADataMakerRec::operator = (const AliCorrQADataMakerRec& qadm ) | |
78 | { | |
57acd2d2 | 79 | // assign operator. |
a64b872d | 80 | this->~AliCorrQADataMakerRec(); |
81 | new(this) AliCorrQADataMakerRec(qadm); | |
82 | return *this; | |
83 | } | |
57acd2d2 | 84 | |
85 | //____________________________________________________________________________ | |
86 | AliCorrQADataMakerRec::~AliCorrQADataMakerRec() | |
87 | { | |
88 | // dtor only destroy the ntuple | |
89 | if ( fCorrNt ) { | |
c6f66d2b | 90 | // for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) { |
91 | // if ( fCorrNt[specie] != NULL ) | |
92 | // delete fCorrNt[specie] ; | |
93 | // } | |
57acd2d2 | 94 | delete[] fCorrNt ; |
96dd7d86 | 95 | fCorrNt = 0x0; |
57acd2d2 | 96 | } |
097485c2 | 97 | if ( fMaxRawVar > 0 ) |
98 | delete [] fVarvalue ; | |
57acd2d2 | 99 | } |
100 | ||
a64b872d | 101 | //____________________________________________________________________________ |
4e25ac79 | 102 | void AliCorrQADataMakerRec::EndOfDetectorCycle(AliQAv1::TASKINDEX_t task, TObjArray ** /*list*/) |
a64b872d | 103 | { |
104 | //Detector specific actions at end of cycle | |
105 | // do the QA checking | |
4e25ac79 | 106 | if (task == AliQAv1::kRAWS) { |
107 | AliQAChecker::Instance()->Run(AliQAv1::kCORR, task, fCorrNt) ; | |
930e6e3e | 108 | } |
a64b872d | 109 | } |
110 | ||
111 | //____________________________________________________________________________ | |
112 | void AliCorrQADataMakerRec::InitESDs() | |
113 | { | |
114 | //Create histograms to controll ESD | |
115 | ||
116 | AliInfo("TO BE IMPLEMENTED") ; | |
117 | } | |
118 | ||
a64b872d | 119 | |
120 | //____________________________________________________________________________ | |
121 | void AliCorrQADataMakerRec::InitRaws() | |
122 | { | |
123 | // createa ntuple taking all the parameters declared by detectors | |
4bb9542a | 124 | if (fCorrNt && fCorrNt[AliRecoParam::AConvert(fEventSpecie)]) |
a64b872d | 125 | return ; |
4bb9542a | 126 | |
127 | if (!fCorrNt) { | |
128 | fCorrNt = new TNtupleD *[AliRecoParam::kNSpecies] ; ; | |
129 | for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) | |
130 | fCorrNt[specie] = NULL ; | |
131 | } | |
132 | ||
3136db6f | 133 | if ( fRawsQAList ) |
134 | { | |
135 | delete[] fRawsQAList ; // not needed for the time being | |
136 | fRawsQAList = NULL ; | |
137 | } | |
a64b872d | 138 | TString varlist("") ; |
4e25ac79 | 139 | for ( Int_t detIndex = 0 ; detIndex < AliQAv1::kNDET ; detIndex++ ) { |
a64b872d | 140 | AliQADataMaker * qadm = fqadm[detIndex] ; |
141 | if ( ! qadm ) | |
142 | continue ; | |
143 | TList * list = qadm->GetParameterList() ; | |
144 | if (list) { | |
145 | TIter next(list) ; | |
146 | TParameter<double> * p ; | |
eca4fa66 | 147 | while ( (p = static_cast<TParameter<double>*>(next()) ) ) { |
a64b872d | 148 | varlist.Append(p->GetName()) ; |
149 | varlist.Append(":") ; | |
150 | fMaxRawVar++ ; | |
151 | } | |
152 | } | |
153 | } | |
154 | varlist = varlist.Strip(TString::kTrailing, ':') ; | |
155 | if (fMaxRawVar == 0) { | |
156 | AliWarning("NTUPLE not created") ; | |
157 | } else { | |
025edd2a | 158 | char * name = Form("%s_%s", AliQAv1::GetQACorrName(), AliRecoParam::GetEventSpecieName(fEventSpecie)) ; |
159 | fCorrNt[AliRecoParam::AConvert(fEventSpecie)] = new TNtupleD(name, "Raws data correlation among detectors", varlist.Data()) ; | |
097485c2 | 160 | fVarvalue = new Double_t[fMaxRawVar] ; |
a64b872d | 161 | } |
162 | } | |
163 | ||
f14c8c46 | 164 | //____________________________________________________________________________ |
165 | void AliCorrQADataMakerRec::InitRecPoints() | |
166 | { | |
167 | // create Reconstructed Points histograms in RecPoints subdir | |
168 | ||
169 | AliInfo("TO BE IMPLEMENTED") ; | |
170 | } | |
171 | ||
172 | //____________________________________________________________________________ | |
173 | void AliCorrQADataMakerRec::InitRecoParams() | |
174 | { | |
175 | // Get the recoparam form the OCDB for every detector involved in CORR | |
176 | ||
177 | AliInfo("TO BE IMPLEMENTED") ; | |
178 | } | |
179 | ||
a64b872d | 180 | //____________________________________________________________________________ |
181 | void AliCorrQADataMakerRec::MakeESDs(AliESDEvent * /*esd*/) | |
182 | { | |
183 | // make QA data from ESDs | |
184 | ||
185 | AliInfo("TO BE IMPLEMENTED") ; | |
186 | ||
187 | } | |
188 | ||
189 | //____________________________________________________________________________ | |
930e6e3e | 190 | void AliCorrQADataMakerRec::MakeRaws(AliRawReader *) |
a64b872d | 191 | { |
192 | //Fill prepared histograms with Raw digit properties | |
025edd2a | 193 | |
4bb9542a | 194 | if ( !fCorrNt || ! fCorrNt[AliRecoParam::AConvert(fEventSpecie)]) |
025edd2a | 195 | InitRaws() ; |
196 | ||
a64b872d | 197 | if ( fMaxRawVar > 0 ) { |
a64b872d | 198 | Int_t index = 0 ; |
4e25ac79 | 199 | for ( Int_t detIndex = 0 ; detIndex < AliQAv1::kNDET ; detIndex++ ) { |
a64b872d | 200 | AliQADataMaker * qadm = fqadm[detIndex] ; |
201 | if ( ! qadm ) | |
202 | continue ; | |
203 | TList * list = qadm->GetParameterList() ; | |
4bb9542a | 204 | if (list) { |
205 | TIter next(list) ; | |
206 | TParameter<double> * p ; | |
207 | while ( (p = static_cast<TParameter<double>*>(next()) ) ) { | |
208 | if (index >= fMaxRawVar) { | |
209 | AliError(Form("Variables list size exceeded (%d) !",index)); | |
210 | break; | |
211 | } | |
212 | fVarvalue[index] = p->GetVal() ; | |
213 | index++ ; | |
214 | } | |
a64b872d | 215 | } |
216 | } | |
097485c2 | 217 | static_cast<TNtupleD*>(fCorrNt[AliRecoParam::AConvert(fEventSpecie)])->Fill(fVarvalue); |
a64b872d | 218 | } |
219 | } | |
220 | ||
221 | //____________________________________________________________________________ | |
222 | void AliCorrQADataMakerRec::MakeRecPoints(TTree * /*clustersTree*/) | |
223 | { | |
224 | AliInfo("TO BE IMPLEMENTED") ; | |
225 | } | |
226 | ||
227 | //____________________________________________________________________________ | |
228 | void AliCorrQADataMakerRec::StartOfDetectorCycle() | |
229 | { | |
230 | //Detector specific actions at start of cycle | |
930e6e3e | 231 | |
a64b872d | 232 | } |