// --- AliRoot header files ---
#include "AliLog.h"
-#include "AliQA.h"
+#include "AliQAv1.h"
#include "AliQAChecker.h"
#include "AliVZEROQAChecker.h"
+#include "AliVZEROQADataMakerRec.h"
//#include "AliCDBEntry.h"
//#include "AliCDBManager.h"
ClassImp(AliVZEROQAChecker)
+//____________________________________________________________________________
+Double_t * AliVZEROQAChecker::Check(AliQAv1::ALITASK_t /*index*/)
+{
+ Double_t * rv = new Double_t[AliRecoParam::kNSpecies] ;
+ for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++)
+ rv[specie] = 0.0 ;
+ return rv ;
+}
+
//__________________________________________________________________
-Double_t AliVZEROQAChecker::Check(AliQA::ALITASK_t index, TObjArray * list)
+Double_t * AliVZEROQAChecker::Check(AliQAv1::ALITASK_t index, TObjArray ** list)
{
// Main check function: Depending on the TASK, different checks will be applied
// }
// Check that histos are filled - (FATAL) set if empty
- if(CheckEntries(list) == 0.0){
- return CheckEntries(list);
- }
-
-// Check for one disk missing (FATAL) or one ring missing (ERROR) in ESDs
- if(index == AliQA::kESD) {
- return CheckEsds(list);
- }
-
- return 1.0;
+ Double_t * check = new Double_t[AliRecoParam::kNSpecies] ;
+ for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
+ check[specie] = 1.0 ;
+ if ( !AliQAv1::Instance()->IsEventSpecieSet(specie) )
+ continue ;
+ if(CheckEntries(list[specie]) == 0.0){
+ check[specie] = CheckEntries(list[specie]);
+ } else {
+ // Check for one disk missing (FATAL) or one ring missing (ERROR) in ESDs
+ if(index == AliQAv1::kESD)
+ check[specie] = CheckEsds(list[specie]);
+ }
+ }
+ return check;
}
+
//_________________________________________________________________
Double_t AliVZEROQAChecker::CheckEntries(TObjArray * list) const
{
-
+
// check on the QA histograms on the input list: list
// list->Print();
Int_t count = 0 ;
if (list->GetEntries() == 0){
- test = 1.0;
- AliInfo(Form("There are NO ENTRIES to be checked..."));
- }
- else {
- TIter next(list) ;
- TH1 * hdata ;
- count = 0 ;
- while ( (hdata = dynamic_cast<TH1 *>(next())) ) {
- if (hdata) {
- Double_t rv = 0.0;
- if(hdata->GetEntries()>0) rv=1.0;
+ test = 1.0;
+ AliInfo(Form("There are NO ENTRIES to be checked..."));
+ } else {
+ TIter next(list) ;
+ TH1 * hdata ;
+ count = 0 ;
+ while ( (hdata = dynamic_cast<TH1 *>(next())) ) {
+ if (hdata) {
+ Double_t rv = 0.0;
+ if(hdata->GetEntries()>0) rv=1.0;
// AliInfo(Form("%s -> %f", hdata->GetName(), rv));
- count++ ; // number of histos
- test += rv ; // number of histos filled
- }
- else{
- AliError(Form("Data type cannot be processed"));
+ count++ ; // number of histos
+ test += rv ; // number of histos filled
+ }else{
+ AliError(Form("Data type cannot be processed"));
}
- }
- if (count != 0) {
- if (test==0.0) {
- AliWarning(Form("Histograms are BOOKED for this specific task, but they are all EMPTY"));
- }
- else {
- test /= count;
- }
- }
+ }
+ if (count != 0) {
+ if (test==0.0) {
+ AliWarning(Form("Histograms are BOOKED for this specific task, but they are all EMPTY"));
+ } else {
+ test /= count;
+ }
+ }
}
return test ;
}
TH1 * hdata ;
while ( (hdata = dynamic_cast<TH1 *>(next())) ) {
- if (hdata) {
- histo_nb++;
- if (histo_nb == 3) {
- Mult_V0A = hdata->GetMean();
-// printf(" Histo ESD number %d; Mean Mult on V0A = %f\n",histo_nb, Mult_V0A);
- }
- if (histo_nb == 4) {
- Mult_V0C = hdata->GetMean();
-// printf(" Histo ESD number %d; Mean Mult on V0C = %f\n",histo_nb, Mult_V0C);
- }
- if (histo_nb == 6) {
- for (Int_t i=0; i<4; i++) {
- V0A_BB_Ring[i] = hdata->Integral((i*8)+1, (i*8) +8);
-// printf(" Histo ESD number %d; Ring = %d; BB A %f\n",histo_nb, i, V0A_BB_Ring[i]);
- }
- }
- if (histo_nb == 7) {
- for (Int_t i=0; i<4; i++) {
- V0A_BG_Ring[i] = hdata->Integral((i*8)+1, (i*8) +8);
-// printf(" Histo ESD number %d; Ring = %d; BG A %f\n",histo_nb, i, V0A_BG_Ring[i]);
- }
- }
- if (histo_nb == 8) {
- for (Int_t i=0; i<4; i++) {
- V0C_BB_Ring[i] = hdata->Integral((i*8)+1, (i*8) +8);
-// printf(" Histo ESD number %d; Ring = %d; BB C %f\n",histo_nb, i, V0C_BB_Ring[i]);
- }
- }
- if (histo_nb == 9) {
- for (Int_t i=0; i<4; i++) {
- V0C_BG_Ring[i] = hdata->Integral((i*8)+1, (i*8) +8);
-// printf(" Histo ESD number %d; Ring = %d; BG C %f\n",histo_nb, i, V0C_BG_Ring[i]);
- }
- }
- }
+ if (hdata) {
+ switch (histo_nb) {
+ case AliVZEROQADataMakerRec::kCellMultiV0A:
+ Mult_V0A = hdata->GetMean();
+ break;
+ case AliVZEROQADataMakerRec::kCellMultiV0C:
+ Mult_V0C = hdata->GetMean();
+ break;
+ case AliVZEROQADataMakerRec::kBBFlag:
+ for (Int_t i=0; i<8; i++) {
+ if(i<4) V0C_BB_Ring[i] = hdata->Integral((i*8)+1, (i*8) +8);
+ else V0A_BB_Ring[i-4] = hdata->Integral((i*8)+1, (i*8) +8);
+ }
+ break;
+ case AliVZEROQADataMakerRec::kBGFlag:
+ for (Int_t i=0; i<8; i++) {
+ if(i<4) V0C_BG_Ring[i] = hdata->Integral((i*8)+1, (i*8) +8);
+ else V0A_BG_Ring[i-4] = hdata->Integral((i*8)+1, (i*8) +8);
+ }
+ break;
+ }
+ }
+ histo_nb++;
}
if(Mult_V0A == 0.0 || Mult_V0C == 0.0) {
return test ;
}
-
+
//______________________________________________________________________________
-void AliVZEROQAChecker::SetQA(AliQA::ALITASK_t index, const Double_t value) const
+void AliVZEROQAChecker::Init(const AliQAv1::DETECTORINDEX_t det)
{
-// sets the QA word according to return value of the Check
+ // intialises QA and QA checker settings
+ AliQAv1::Instance(det) ;
+ Float_t * hiValue = new Float_t[AliQAv1::kNBIT] ;
+ Float_t * lowValue = new Float_t[AliQAv1::kNBIT] ;
+ lowValue[AliQAv1::kINFO] = 0.5 ;
+ hiValue[AliQAv1::kINFO] = 1.0 ;
+ lowValue[AliQAv1::kWARNING] = 0.2 ;
+ hiValue[AliQAv1::kWARNING] = 0.5 ;
+ lowValue[AliQAv1::kERROR] = 0.0 ;
+ hiValue[AliQAv1::kERROR] = 0.2 ;
+ lowValue[AliQAv1::kFATAL] = -1.0 ;
+ hiValue[AliQAv1::kFATAL] = 0.0 ;
+ SetHiLo(hiValue, lowValue) ;
+}
- AliQA * qa = AliQA::Instance(index);
-
- qa->UnSet(AliQA::kFATAL);
- qa->UnSet(AliQA::kWARNING);
- qa->UnSet(AliQA::kERROR);
- qa->UnSet(AliQA::kINFO);
-
- if ( value <= 0.0 ) qa->Set(AliQA::kFATAL) ;
- else if ( value > 0.0 && value <= 0.2 ) qa->Set(AliQA::kERROR) ;
- else if ( value > 0.2 && value <= 0.5 ) qa->Set(AliQA::kWARNING) ;
- else if ( value > 0.5 && value < 1.0 ) qa->Set(AliQA::kINFO) ;
+//______________________________________________________________________________
+void AliVZEROQAChecker::SetQA(AliQAv1::ALITASK_t index, Double_t * value) const
+{
+// sets the QA word according to return value of the Check
+ AliQAv1 * qa = AliQAv1::Instance(index);
+ for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
+ qa->UnSet(AliQAv1::kFATAL, specie);
+ qa->UnSet(AliQAv1::kWARNING, specie);
+ qa->UnSet(AliQAv1::kERROR, specie);
+ qa->UnSet(AliQAv1::kINFO, specie);
+ if ( ! value ) { // No checker is implemented, set all QA to Fatal
+ qa->Set(AliQAv1::kFATAL, specie) ;
+ } else {
+ if ( value[specie] >= fLowTestValue[AliQAv1::kFATAL] && value[specie] < fUpTestValue[AliQAv1::kFATAL] )
+ qa->Set(AliQAv1::kFATAL, specie) ;
+ else if ( value[specie] > fLowTestValue[AliQAv1::kERROR] && value[specie] <= fUpTestValue[AliQAv1::kERROR] )
+ qa->Set(AliQAv1::kERROR, specie) ;
+ else if ( value[specie] > fLowTestValue[AliQAv1::kWARNING] && value[specie] <= fUpTestValue[AliQAv1::kWARNING] )
+ qa->Set(AliQAv1::kWARNING, specie) ;
+ else if ( value[specie] > fLowTestValue[AliQAv1::kINFO] && value[specie] <= fUpTestValue[AliQAv1::kINFO] )
+ qa->Set(AliQAv1::kINFO, specie) ;
+ }
+ }
}