MUON/AliMUONTriggerQAChecker.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 "AliMUONTriggerQAChecker.h"
  19
  20 /// \class AliMUONTriggerQAChecker
  21 ///
  22 /// Implementation of QAChecker for MTR
  23 ///
  24 /// For the moment we only implement the checking of raw data QA for the trigger
  25 /// by looking at the local structure and trigger response errors.
  26 ///
  27 /// \author Diego Stocco, Subatech
  28
  29
  30 #include "AliRecoParam.h"
  31 #include "AliMUONQAIndices.h"
  32 #include "AliQAv1.h"
  33 #include "TH1.h"
  34 #include "TPaveText.h"
  35 #include "TString.h"
  36 #include "TObjArray.h"
  37 #include "TList.h"
  38
  39 /// \cond CLASSIMP
  40 ClassImp(AliMUONTriggerQAChecker)
  41 /// \endcond
  42
  43 //__________________________________________________________________
  44 AliMUONTriggerQAChecker::AliMUONTriggerQAChecker() : AliMUONVQAChecker()
  45 {
  46         /// ctor
  47 }
  48
  49 //__________________________________________________________________
  50 AliMUONTriggerQAChecker::~AliMUONTriggerQAChecker()
  51 {
  52         /// dtor
  53 }
  54
  55 //______________________________________________________________________________
  56 AliMUONVQAChecker::ECheckCode
  57 AliMUONTriggerQAChecker::MarkHisto(TH1& histo, AliMUONVQAChecker::ECheckCode value) const
  58 {
  59   /// Mark histo as originator of some QA error/warning
  60
  61   if ( value != AliMUONVQAChecker::kInfo )
  62   {
  63     histo.SetBit(AliQAv1::GetQABit());
  64   }
  65
  66   return value;
  67 }
  68
  69 //__________________________________________________________________
  70 AliMUONVQAChecker::ECheckCode*
  71 AliMUONTriggerQAChecker::CheckRaws(TObjArray** list, const AliMUONRecoParam* )
  72 {
  73   /// Check raw data
  74
  75   AliMUONVQAChecker::ECheckCode * rv = new AliMUONVQAChecker::ECheckCode[AliRecoParam::kNSpecies];
  76
  77   Int_t histoRawsPercentIndex[] = {
  78     AliMUONQAIndices::kTriggerErrorSummaryNorm,
  79     AliMUONQAIndices::kTriggerCalibSummaryNorm,
  80     AliMUONQAIndices::kTriggerReadOutErrorsNorm
  81   };
  82   const Int_t kNrawsHistos = sizeof(histoRawsPercentIndex)/sizeof(histoRawsPercentIndex[0]);
  83
  84   // BEGIN OF LIMITS
  85   // Fixme: Move me to reference histos
  86   Float_t safeFactor = 5.;
  87   Float_t warningPercentTrigAlgo[AliMUONQAIndices::kNtrigAlgoErrorBins] = {safeFactor*1., safeFactor*1., safeFactor*1., 100., 100., 100., 100., safeFactor*1., safeFactor*1., safeFactor*1.};
  88   Float_t warningPercentCalib[AliMUONQAIndices::kNtrigCalibSummaryBins] = {safeFactor*0.4, safeFactor*1., 6.2, 0.0001, safeFactor*0.4};
  89   Float_t warningPercentReadout[AliMUONQAIndices::kNtrigStructErrorBins] = {safeFactor*1., safeFactor*1., safeFactor*1., safeFactor*1.};
  90
  91   Float_t* warningPercent[kNrawsHistos] = {warningPercentTrigAlgo, warningPercentCalib, warningPercentReadout};
  92
  93   Float_t errorFactor = 30.;
  94   Float_t errorPercentTrigAlgo[AliMUONQAIndices::kNtrigAlgoErrorBins] = {errorFactor*1., errorFactor*1., errorFactor*1., 100., 100., 100., 100., errorFactor*1., errorFactor*1., errorFactor*1.};
  95   Float_t errorPercentCalib[AliMUONQAIndices::kNtrigCalibSummaryBins] = {errorFactor*0.4, errorFactor*1., 3.*6.2, 3.*0.0001, errorFactor*0.4};
  96   Float_t errorPercentReadout[AliMUONQAIndices::kNtrigStructErrorBins] = {errorFactor*1., errorFactor*1., errorFactor*1., errorFactor*1.};
  97   // END OF LIMTS
  98
  99   Float_t* errorPercent[kNrawsHistos] = {errorPercentTrigAlgo, errorPercentCalib, errorPercentReadout};
 100
 101   TObjArray messages;
 102   messages.SetOwner(kTRUE);
 103
 104   TH1* currHisto = 0x0;
 105   for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
 106     rv[specie] = AliMUONVQAChecker::kInfo;
 107
 108     TH1* hAnalyzedEvents = AliQAv1::GetData(list,AliMUONQAIndices::kTriggerRawNAnalyzedEvents,AliRecoParam::ConvertIndex(specie));
 109     Int_t nAnalyzedEvents = 0;
 110     if ( hAnalyzedEvents )
 111       nAnalyzedEvents = TMath::Nint(hAnalyzedEvents->GetBinContent(1));
 112
 113     //if ( nAnalyzedEvents == 0 ) rv[specie] = AliMUONVQAChecker::kFatal;
 114
 115     for(Int_t ihisto = 0; ihisto<kNrawsHistos; ihisto++){
 116       AliMUONVQAChecker::ECheckCode currRv = AliMUONVQAChecker::kInfo;
 117       messages.Clear();
 118       currHisto = AliQAv1::GetData(list,histoRawsPercentIndex[ihisto],AliRecoParam::ConvertIndex(specie));
 119       if ( ! currHisto ) continue;
 120
 121       Int_t nbins = currHisto->GetXaxis()->GetNbins();
 122       for (Int_t ibin = 1; ibin<=nbins; ibin++){
 123         Double_t binContent = currHisto->GetBinContent(ibin);
 124         if ( binContent > errorPercent[ihisto][ibin-1] )
 125           currRv = AliMUONVQAChecker::kError;
 126         else if ( binContent > warningPercent[ihisto][ibin-1] )
 127           currRv = AliMUONVQAChecker::kWarning;
 128         else if ( ibin == 4 && binContent > 50. && AliRecoParam::ConvertIndex(specie) == AliRecoParam::kCalib) {
 129           messages.Add(new TObjString("Do not panic:"));
 130           messages.Add(new TObjString("copy errors do not affect data"));
 131         }
 132       } // loop on bins
 133       if ( currRv != AliMUONVQAChecker::kInfo ) {
 134         switch ( histoRawsPercentIndex[ihisto] ) {
 135         case AliMUONQAIndices::kTriggerErrorSummaryNorm:
 136         case AliMUONQAIndices::kTriggerCalibSummaryNorm:
 137           messages.Add(new TObjString("Trigger algo errors"));
 138           break;
 139         case AliMUONQAIndices::kTriggerReadOutErrorsNorm:
 140           messages.Add(new TObjString("Readout errors"));
 141           break;
 142         }
 143         if ( currRv == AliMUONVQAChecker::kWarning )
 144           messages.Add(new TObjString("are a little bit high"));
 145         else if ( currRv == AliMUONVQAChecker::kError ||
 146                   currRv == AliMUONVQAChecker::kFatal )
 147           messages.Add(new TObjString("are too high"));
 148       }
 149       else if ( nAnalyzedEvents != 0 ) {
 150         switch ( histoRawsPercentIndex[ihisto] ) {
 151         case AliMUONQAIndices::kTriggerErrorSummaryNorm:
 152         case AliMUONQAIndices::kTriggerCalibSummaryNorm:
 153         case AliMUONQAIndices::kTriggerReadOutErrorsNorm:
 154           messages.Add(new TObjString("Values within limits"));
 155           break;
 156         }
 157       }
 158       if ( MarkHisto(*currHisto, currRv) < rv[specie] )
 159         rv[specie] = currRv;
 160       currHisto->GetYaxis()->SetRangeUser(0., 110.);
 161       SetupHisto(nAnalyzedEvents, messages, *currHisto, currRv);
 162     } // loop on histos
 163   } // loop on species
 164
 165   return rv;
 166 }
 167
 168 //__________________________________________________________________
 169 AliMUONVQAChecker::ECheckCode*
 170 AliMUONTriggerQAChecker::CheckRecPoints(TObjArray** , const AliMUONRecoParam* )
 171 {
 172   /// Check rec points
 173   return 0x0;
 174 }
 175
 176 //__________________________________________________________________
 177 AliMUONVQAChecker::ECheckCode*
 178 AliMUONTriggerQAChecker::CheckESD(TObjArray** , const AliMUONRecoParam* )
 179 {
 180   /// Check esd
 181   return 0x0;
 182 }
 183
 184
 185 //___________________________________________________________________
 186 void AliMUONTriggerQAChecker::SetupHisto(Int_t nevents, const TObjArray& messages, TH1& histo, AliMUONVQAChecker::ECheckCode code)
 187 {
 188   //
 189   /// Add text to histos
 190   //
 191
 192   Double_t y1 = 0.87 - (messages.GetLast()+2)*0.075;
 193   TPaveText* text = new TPaveText(0.25,y1,0.75,0.89,"NDC");
 194
 195   text->AddText(Form("MTR - Total events %i", nevents));
 196
 197   TIter next(&messages);
 198   TObjString* str;
 199
 200   while ( ( str = static_cast<TObjString*>(next()) ) ){
 201     text->AddText(str->String());
 202   }
 203
 204   TString defaultText = "";
 205   Int_t color = 0;
 206
 207   if ( nevents == 0 ) {
 208     color = AliMUONVQAChecker::kFatalColor;
 209     defaultText = "No events analyzed!";
 210   }
 211   else if ( nevents <= 20 ) {
 212     color = AliMUONVQAChecker::kWarningColor;
 213     text->AddText("Not enough events to judge");
 214     text->AddText("Please wait for more statistics");
 215   }
 216   else {
 217     switch ( code ) {
 218       case AliMUONVQAChecker::kInfo:
 219         color = AliMUONVQAChecker::kInfoColor;
 220         defaultText = "All is fine!";
 221         break;
 222       case AliMUONVQAChecker::kWarning:
 223         color = AliMUONVQAChecker::kWarningColor;
 224         defaultText = "Please keep an eye on it!";
 225         break;
 226       case AliMUONVQAChecker::kFatal:
 227         color = AliMUONVQAChecker::kFatalColor;
 228         defaultText = "PLEASE CALL MUON TRIGGER EXPERT!!!";
 229         break;
 230       default:
 231         color = AliMUONVQAChecker::kErrorColor;
 232         defaultText = "PLEASE CALL MUON TRIGGER EXPERT!";
 233         break;
 234     }
 235   }
 236
 237   text->AddText(defaultText.Data());
 238   text->SetFillColor(color);
 239
 240   histo.SetFillStyle(1001);
 241   histo.SetFillColor(color);
 242
 243   histo.SetStats(kFALSE);
 244
 245   histo.GetListOfFunctions()->Delete();
 246   histo.GetListOfFunctions()->Add(text);
 247 }