1 /**************************************************************************
2 * Copyright(c) 1998-1999, 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 **************************************************************************/
17 /* $Id: AliCorrQADataMakerRec.cxx 27570 2008-07-24 21:49:27Z cvetan $ */
20 Produces the data needed to calculate the quality assurance.
21 All data must be mergeable objects.
22 Y. Schutz CERN July 2007
25 // --- ROOT system ---
26 #include <TClonesArray.h>
32 #include <TParameter.h>
35 // --- Standard library ---
37 // --- AliRoot header files ---
39 #include "AliCorrQADataMakerRec.h"
40 #include "AliQAChecker.h"
42 ClassImp(AliCorrQADataMakerRec)
44 //____________________________________________________________________________
45 AliCorrQADataMakerRec::AliCorrQADataMakerRec(AliQADataMaker ** qadm ) :
46 AliQADataMakerRec(AliQAv1::GetDetName(AliQAv1::kCORR), "Corr Quality Assurance Data Maker"),
52 fCorrNt = new TNtupleD *[AliRecoParam::kNSpecies] ;
53 for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++)
54 fCorrNt[specie] = NULL ;
57 //____________________________________________________________________________
58 AliCorrQADataMakerRec::AliCorrQADataMakerRec(const AliCorrQADataMakerRec& qadm) :
60 fMaxRawVar(qadm.fMaxRawVar),
65 SetName((const char*)qadm.GetName()) ;
66 SetTitle((const char*)qadm.GetTitle());
68 fVarvalue = new Double_t[fMaxRawVar] ;
70 fCorrNt = new TNtupleD *[AliRecoParam::kNSpecies] ;
71 for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++)
72 fCorrNt[specie] = qadm.fCorrNt[specie] ;
76 //__________________________________________________________________
77 AliCorrQADataMakerRec& AliCorrQADataMakerRec::operator = (const AliCorrQADataMakerRec& qadm )
80 this->~AliCorrQADataMakerRec();
81 new(this) AliCorrQADataMakerRec(qadm);
85 //____________________________________________________________________________
86 AliCorrQADataMakerRec::~AliCorrQADataMakerRec()
88 // dtor only destroy the ntuple
90 // for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
91 // if ( fCorrNt[specie] != NULL )
92 // delete fCorrNt[specie] ;
101 //____________________________________________________________________________
102 void AliCorrQADataMakerRec::EndOfDetectorCycle(AliQAv1::TASKINDEX_t task, TObjArray ** /*list*/)
104 //Detector specific actions at end of cycle
105 // do the QA checking
106 if (task == AliQAv1::kRAWS) {
107 AliQAChecker::Instance()->Run(AliQAv1::kCORR, task, fCorrNt) ;
111 //____________________________________________________________________________
112 void AliCorrQADataMakerRec::InitESDs()
114 //Create histograms to controll ESD
116 AliInfo("TO BE IMPLEMENTED") ;
120 //____________________________________________________________________________
121 void AliCorrQADataMakerRec::InitRaws()
123 // createa ntuple taking all the parameters declared by detectors
124 if (fCorrNt && fCorrNt[AliRecoParam::AConvert(fEventSpecie)])
128 fCorrNt = new TNtupleD *[AliRecoParam::kNSpecies] ; ;
129 for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++)
130 fCorrNt[specie] = NULL ;
135 delete[] fRawsQAList ; // not needed for the time being
138 TString varlist("") ;
139 for ( Int_t detIndex = 0 ; detIndex < AliQAv1::kNDET ; detIndex++ ) {
140 AliQADataMaker * qadm = fqadm[detIndex] ;
143 TList * list = qadm->GetParameterList() ;
146 TParameter<double> * p ;
147 while ( (p = static_cast<TParameter<double>*>(next()) ) ) {
148 varlist.Append(p->GetName()) ;
149 varlist.Append(":") ;
154 varlist = varlist.Strip(TString::kTrailing, ':') ;
155 if (fMaxRawVar == 0) {
156 AliWarning("NTUPLE not created") ;
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()) ;
160 fVarvalue = new Double_t[fMaxRawVar] ;
164 //____________________________________________________________________________
165 void AliCorrQADataMakerRec::InitRecPoints()
167 // create Reconstructed Points histograms in RecPoints subdir
169 AliInfo("TO BE IMPLEMENTED") ;
172 //____________________________________________________________________________
173 void AliCorrQADataMakerRec::InitRecoParams()
175 // Get the recoparam form the OCDB for every detector involved in CORR
177 AliInfo("TO BE IMPLEMENTED") ;
180 //____________________________________________________________________________
181 void AliCorrQADataMakerRec::MakeESDs(AliESDEvent * /*esd*/)
183 // make QA data from ESDs
185 AliInfo("TO BE IMPLEMENTED") ;
189 //____________________________________________________________________________
190 void AliCorrQADataMakerRec::MakeRaws(AliRawReader *)
192 //Fill prepared histograms with Raw digit properties
194 if ( !fCorrNt || ! fCorrNt[AliRecoParam::AConvert(fEventSpecie)])
197 if ( fMaxRawVar > 0 ) {
199 for ( Int_t detIndex = 0 ; detIndex < AliQAv1::kNDET ; detIndex++ ) {
200 AliQADataMaker * qadm = fqadm[detIndex] ;
203 TList * list = qadm->GetParameterList() ;
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));
212 fVarvalue[index] = p->GetVal() ;
217 static_cast<TNtupleD*>(fCorrNt[AliRecoParam::AConvert(fEventSpecie)])->Fill(fVarvalue);
221 //____________________________________________________________________________
222 void AliCorrQADataMakerRec::MakeRecPoints(TTree * /*clustersTree*/)
224 AliInfo("TO BE IMPLEMENTED") ;
227 //____________________________________________________________________________
228 void AliCorrQADataMakerRec::StartOfDetectorCycle()
230 //Detector specific actions at start of cycle