MUON/AliMUONQAChecker.cxx

   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 // $Id$
  17
  18 #include "AliMUONQAChecker.h"
  19
  20 /// \class AliMUONQAChecker
  21 ///
  22 /// Implementation of AliQACheckerBase for MCH and MTR
  23 ///
  24 /// For the moment we only implement the checking of raw data QA for the tracker
  25 /// by looking at the occupancy at the manu level.
  26 /// We count the number of manus above a given occupancy threshold, and
  27 /// depending on that number, the resulting QA flag is warning, error or fatal.
  28 /// (there's no "info" type in this case).
  29 ///
  30 /// \author Laurent Aphecetche, Subatech
  31
  32 #include "AliLog.h"
  33 #include "AliMUONVTrackerData.h"
  34 #include "AliMpManuIterator.h"
  35 #include "AliQAv1.h"
  36 #include <TDirectory.h>
  37 #include <TH1.h>
  38
  39 /// \cond CLASSIMP
  40 ClassImp(AliMUONQAChecker)
  41 /// \endcond
  42
  43 //__________________________________________________________________
  44 AliMUONQAChecker::AliMUONQAChecker() :
  45     AliQACheckerBase("MUON","MUON Quality Assurance Data Maker")
  46 {
  47         /// ctor
  48 }
  49
  50 //__________________________________________________________________
  51 AliMUONQAChecker::~AliMUONQAChecker()
  52 {
  53         /// dtor
  54 }
  55
  56 //__________________________________________________________________
  57 AliMUONQAChecker::AliMUONQAChecker(const AliMUONQAChecker& qac) :
  58     AliQACheckerBase(qac.GetName(), qac.GetTitle())
  59 {
  60         /// copy ctor
  61 }
  62
  63 //______________________________________________________________________________
  64 Double_t *
  65 AliMUONQAChecker::Check(AliQAv1::ALITASK_t index, TObjArray ** list)
  66 {
  67   /// Check objects in list
  68
  69   if ( index == AliQAv1::kRAW )
  70   {
  71     return CheckRaws(list);
  72   }
  73
  74   if ( index == AliQAv1::kREC)
  75   {
  76     return CheckRecPoints(list);
  77   }
  78
  79   if ( index == AliQAv1::kESD )
  80   {
  81     return CheckESD(list);
  82   }
  83
  84   AliWarning(Form("Checker for task %d not implement for the moment",index));
  85   return NULL;
  86 }
  87
  88 //______________________________________________________________________________
  89 TH1*
  90 AliMUONQAChecker::GetHisto(TObjArray* list, const char* hname, Int_t specie) const
  91 {
  92   /// Get a given histo from the list
  93   TH1* h = static_cast<TH1*>(list->FindObject(Form("%s_%s",AliRecoParam::GetEventSpecieName(specie),hname)));
  94   if (!h)
  95   {
  96     AliError(Form("Did not find expected histo %s",hname));
  97   }
  98   return h;
  99 }
 100
 101 //______________________________________________________________________________
 102 Double_t *
 103 AliMUONQAChecker::CheckRecPoints(TObjArray ** list)
 104 {
 105   /// Check rec points
 106   /// Very binary check for the moment.
 107
 108   Double_t * rv = new Double_t[AliRecoParam::kNSpecies] ;
 109   for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++)
 110     rv[specie] = 1.0 ;
 111
 112   for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
 113     TH1* h = GetHisto(list[specie],"hTrackerNumberOfClustersPerDE",specie);
 114
 115     if ( !h ) rv[specie] =  0.75; // only a warning if histo not found, in order not to kill anything because QA failed...
 116
 117     else if ( h->GetMean() == 0.0 ) rv[specie] =  MarkHisto(*h,0.0);
 118   }
 119   return rv;
 120 }
 121
 122 //______________________________________________________________________________
 123 Double_t
 124 AliMUONQAChecker::MarkHisto(TH1& histo, Double_t value) const
 125 {
 126   /// Mark histo as originator of some QA error/warning
 127
 128   if ( value != 1.0 )
 129   {
 130     histo.SetBit(AliQAv1::GetQABit());
 131   }
 132
 133   return value;
 134 }
 135
 136 //______________________________________________________________________________
 137 Double_t *
 138 AliMUONQAChecker::CheckESD(TObjArray ** list)
 139 {
 140   /// Check ESD
 141
 142   Double_t * rv = new Double_t[AliRecoParam::kNSpecies] ;
 143   for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++)
 144     rv[specie] = 1.0 ;
 145
 146   for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
 147
 148     TH1* h = GetHisto(list[specie],"hESDnTracks",specie);
 149
 150     if (!h) rv[specie] = 0.75;
 151
 152     else if ( h->GetMean() == 0.0 ) rv[specie] =  MarkHisto(*h,0.0); // no track -> fatal
 153
 154     h = GetHisto(list[specie],"hESDMatchTrig",specie);
 155
 156     if (!h) rv[specie] =  0.75;
 157
 158     else if (h->GetMean() == 0.0 ) rv[specie] = MarkHisto(*h,0.25); // no trigger matching -> error
 159   }
 160   return rv;
 161 }
 162
 163 //______________________________________________________________________________
 164 Double_t *
 165 AliMUONQAChecker::CheckRaws(TObjArray ** list)
 166 {
 167   /// Check raws
 168
 169   Double_t * rv = new Double_t[AliRecoParam::kNSpecies] ;
 170   for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++)
 171     rv[specie] = 1.0 ;
 172
 173   for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
 174     TIter next(list[specie]);
 175     TObject* object;
 176     AliMUONVTrackerData* data(0x0);
 177
 178     while ( (object=next()) && !data )
 179       {
 180         if (object->InheritsFrom("AliMUONVTrackerData"))
 181           {
 182             data = static_cast<AliMUONVTrackerData*>(object);
 183           }
 184       }
 185
 186     if ( !data )
 187       {
 188         AliError("Did not find TrackerData in the list !");
 189         return NULL;
 190       }
 191
 192     AliMpManuIterator it;
 193     Int_t detElemId;
 194     Int_t manuId;
 195     Int_t n50(0); // number of manus with occupancy above 0.5
 196     Int_t n75(0); // number of manus with occupancy above 0.75
 197     Int_t n(0); // number of manus with some occupancy
 198
 199     while ( it.Next(detElemId,manuId) )
 200       {
 201         Float_t occ = data->Manu(detElemId,manuId,2);
 202         if (occ > 0 ) ++n;
 203         if (occ >= 0.5 ) ++n50;
 204         if (occ >= 0.75 ) ++n75;
 205       }
 206
 207     AliInfo(Form("n %d n50 %d n75 %d",n,n50,n75));
 208
 209     if ( n == 0 )
 210       {
 211         AliError("Oups. Got zero occupancy in all manus ?!");
 212         rv[specie] =  0.0;
 213       }
 214
 215     if ( n75 )
 216       {
 217         AliError(Form("Got %d manus with occupancy above 0.75",n75));
 218         rv[specie] =  0.1;
 219       }
 220
 221     if ( n50 )
 222       {
 223         AliWarning(Form("Got %d manus with occupancy above 0.5",n50));
 224         rv[specie] =  0.9;
 225       }
 226   }
 227   return rv;
 228 }
 229
 230 //______________________________________________________________________________
 231 void AliMUONQAChecker::Init(const AliQAv1::DETECTORINDEX_t det)
 232 {
 233   /// intialises QA and QA checker settings
 234   AliQAv1::Instance(det) ;
 235   Float_t hiValue[AliQAv1::kNBIT] ;
 236   Float_t lowValue[AliQAv1::kNBIT] ;
 237   lowValue[AliQAv1::kINFO]      = 0.999   ;
 238   hiValue[AliQAv1::kINFO]       = 1.0 ;
 239   hiValue[AliQAv1::kWARNING]    = 0.99 ;
 240   lowValue[AliQAv1::kWARNING]   = 0.5 ;
 241   lowValue[AliQAv1::kERROR]     = 0.0   ;
 242   hiValue[AliQAv1::kERROR]      = 0.5 ;
 243   lowValue[AliQAv1::kFATAL]     = -1.0   ;
 244   hiValue[AliQAv1::kFATAL]      = 0.0 ;
 245   SetHiLo(&hiValue[0], &lowValue[0]) ;
 246 }
 247
 248 //______________________________________________________________________________
 249 void
 250 AliMUONQAChecker::SetQA(AliQAv1::ALITASK_t index, Double_t * value) const
 251 {
 252   /// sets the QA according the return value of the Check
 253
 254   AliQAv1 * qa = AliQAv1::Instance(index);
 255
 256   for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
 257     qa->UnSet(AliQAv1::kFATAL, specie);
 258     qa->UnSet(AliQAv1::kWARNING, specie);
 259     qa->UnSet(AliQAv1::kERROR, specie);
 260     qa->UnSet(AliQAv1::kINFO, specie);
 261
 262     if ( ! value ) { // No checker is implemented, set all QA to Fatal
 263       qa->Set(AliQAv1::kFATAL, specie) ;
 264     } else {
 265       if ( value[specie] >= fLowTestValue[AliQAv1::kFATAL] && value[specie] < fUpTestValue[AliQAv1::kFATAL] )
 266         qa->Set(AliQAv1::kFATAL, specie) ;
 267       else if ( value[specie] > fLowTestValue[AliQAv1::kERROR] && value[specie] <= fUpTestValue[AliQAv1::kERROR]  )
 268         qa->Set(AliQAv1::kERROR, specie) ;
 269       else if ( value[specie] > fLowTestValue[AliQAv1::kWARNING] && value[specie] <= fUpTestValue[AliQAv1::kWARNING]  )
 270         qa->Set(AliQAv1::kWARNING, specie) ;
 271       else if ( value[specie] > fLowTestValue[AliQAv1::kINFO] && value[specie] <= fUpTestValue[AliQAv1::kINFO] )
 272         qa->Set(AliQAv1::kINFO, specie) ;
 273     }
 274   }
 275 }