T0/AliT0QAChecker.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 //  Checks the quality assurance.
  17 //  By comparing with reference data
  18 //  Skeleton for T0
  19 //---------------------------------------------
  20 //checkig without reference data:
  21 //for RAW QA all histograms should have approximatly the same
  22 //number of entries as RefPoint
  23 //for Rec Points checks
  24 //  - amplitude measured by 2 methos
  25 //  - online and offline T0 measurements
  26 // for ESD quality of reconstruction ( and measurements):
  27 // RMS of vertex and T0 less than 75ps
  28 //
  29 // Alla.Maevskaya@cern.ch
  30 //...
  31
  32 // --- ROOT system ---
  33 #include <Riostream.h>
  34 #include <TClass.h>
  35 #include <TH1F.h>
  36 #include <TH2.h>
  37 #include <TIterator.h>
  38 #include <TKey.h>
  39 #include <TFile.h>
  40 #include <TMath.h>
  41 #include <TString.h>
  42
  43 // --- Standard library ---
  44
  45 // --- AliRoot header files ---
  46 #include "AliLog.h"
  47 #include "AliQAv1.h"
  48 #include "AliQAChecker.h"
  49 #include "AliCDBEntry.h"
  50 #include "AliQAManager.h"
  51 #include "AliT0QAChecker.h"
  52
  53 ClassImp(AliT0QAChecker)
  54
  55 //__________________________________________________________________
  56 Double_t * AliT0QAChecker::Check(AliQAv1::ALITASK_t index, TObjArray ** list, const AliDetectorRecoParam * /*recoParam*/)
  57 {
  58
  59   // Super-basic check on the QA histograms on the input list:
  60   // look whether they are empty!
  61
  62   Double_t * test = new Double_t[AliRecoParam::kNSpecies] ;
  63
  64   char * detOCDBDir = Form("T0/%s/%s", AliQAv1::GetRefOCDBDirName(), AliQAv1::GetRefDataDirName()) ;
  65
  66   AliCDBEntry *QARefRec = AliQAManager::QAManager()->Get(detOCDBDir);
  67   //  QARefRec->Dump();
  68   if( !QARefRec){
  69     AliInfo("QA reference data NOT retrieved for Reconstruction check. No T0 reference distribution");
  70   }
  71
  72    for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++)
  73     test[specie]    = 10.0 ;
  74
  75   Double_t nent[500];
  76   TString hname[500];
  77   const char *cname;
  78   memset(nent,0,500*sizeof(Double_t));
  79   Double_t w[500];
  80   memset(w,1,500*sizeof(Double_t));
  81   TH2 *fhRecDiff[3];
  82   TH1 *fhRawEff[500];
  83   TH2 *fhRawTime[500];
  84   TH1 *fhESD[2];
  85   TH2 *fhHits[10];
  86
  87   for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
  88     //  TString dataType = AliQAv1::GetAliTaskName(index);
  89     if (list[specie]->GetEntries() == 0){
  90       test[specie] = 1. ; // nothing to check
  91     }
  92     else {
  93       TH1 * hdata ;
  94       TH2 * h ;
  95       for (Int_t ir=0; ir<list[specie]->GetEntries(); ir++) {
  96
  97         //hits
  98
  99         if(index == AliQAv1::kSIM && AliQAv1::GetTaskName(AliQAv1::kHITS)){
 100           h =  (TH2*) list[specie]->UncheckedAt(ir);
 101           cname = h->GetName();
 102           hname[ir] = cname;
 103           fhHits[ir] = h;
 104         }
 105
 106         //raw
 107       if(index == AliQAv1::kRAW ){
 108         if(ir > 204 && ir<208 ) {
 109           hdata = (TH1*) list[specie]->UncheckedAt(ir);
 110           if(hdata) {
 111             cname = hdata->GetName();
 112             hname[ir] = cname;
 113             fhRawEff[ir] = hdata;
 114           }
 115         }
 116         if((ir>207 && ir<210)  && specie == AliRecoParam::kCalib) {
 117           h =  (TH2*) list[specie]->UncheckedAt(ir);
 118           if(h) {
 119             cname = h->GetName();
 120             hname[ir] = cname;
 121             fhRawTime[ir] = h;
 122           AliDebug(AliQAv1::GetQADebugLevel(), Form("count %i %s ",ir, hname[ir].Data()) );
 123           }
 124         }
 125       }
 126       //rec
 127       if(index == AliQAv1::kREC){
 128         h = (TH2*) list[specie]->UncheckedAt(ir);
 129         if(h) {
 130           cname = h->GetName();
 131           hname[ir] = cname;
 132            fhRecDiff[ir] = h;
 133         }
 134       }
 135       //esd
 136       if(index ==  AliQAv1::kESD){
 137         hdata = (TH1*) list[specie]->UncheckedAt(ir);
 138         if(hdata){
 139           fhESD[ir] = hdata;
 140           AliDebug(AliQAv1::GetQADebugLevel(), Form("count %i %s ",ir, hname[ir].Data()) );
 141         }
 142       }
 143     }
 144
 145
 146
 147     //raw data
 148
 149       if (index == AliQAv1::kRAW && specie == AliRecoParam::kCalib){
 150         test[specie] = CheckRaw(list[specie],dynamic_cast<TObjArray*>(dynamic_cast<TList *>(QARefRec->GetObject())->First()));
 151       }
 152
 153       if(index == AliQAv1::kREC){
 154         //rec points
 155         for (Int_t icase=0; icase<2; icase++) {
 156           for (Int_t idet=0; idet<24; idet++) {
 157             Double_t mean = fhRecDiff[icase]->
 158             ProjectionY(Form("%s_py", fhRecDiff[icase]->GetName()),
 159                         idet,idet+1)->GetMean();
 160             Double_t rms= fhRecDiff[icase]->
 161             ProjectionY(Form("%s_py", fhRecDiff[icase]->GetName()),
 162                         idet,idet+1)->GetRMS();
 163             AliDebug(AliQAv1::GetQADebugLevel(), Form("name %s icase %i idet %i mean %f, rms %f\n",
 164                          fhRecDiff[icase]->GetName(), icase, idet, mean,rms));
 165
 166             if(TMath::Abs(mean) >1.5 || rms >1){
 167               AliDebug(AliQAv1::GetQADebugLevel(), Form(" calibration is nor perfect; test=%f", test)) ;
 168               test[specie]=0.25;
 169             }
 170             if(mean>3 || rms >5) {
 171               test[specie] = 0.1;
 172               AliDebug(AliQAv1::GetQADebugLevel(), Form(" wrong calibration test=%f", test[specie])) ;
 173             }
 174           }
 175         }
 176       }
 177       if (index == AliQAv1::kESD) {
 178         //ESD
 179         for (Int_t icase=0; icase<2; icase++) {
 180           Double_t rmsVertex = fhESD[icase]->GetRMS();
 181           Double_t meanVertex = fhESD[icase]->GetMean();
 182           test[specie]=1;
 183           AliDebug(AliQAv1::GetQADebugLevel(), Form("numentries %d meanVertex %f rmsVertex %f", fhESD[icase]->GetEntries(), meanVertex, rmsVertex));
 184           if (TMath::Abs(rmsVertex)>3) {
 185             test[specie]=0.25;
 186             AliDebug(AliQAv1::GetQADebugLevel(), Form("Vertex position resolution not good  , rms= %f test=%f",
 187                             rmsVertex, test[specie])) ;
 188           }
 189           if (TMath::Abs(meanVertex)>3) {
 190             test[specie]=0.25;
 191             AliDebug(AliQAv1::GetQADebugLevel(), Form("Vertex position bad calibrated  , Mean= %f test=%f",
 192                             meanVertex, test[specie])) ;
 193           }
 194         }
 195       }
 196     } //  if (list->GetEntries() != 0
 197     AliDebug(AliQAv1::GetQADebugLevel(), Form("Test Result = %f", test[specie])) ;
 198   }
 199
 200   return test ;
 201 }
 202
 203 //--------------------------------------------------------------------------
 204 Double_t AliT0QAChecker::CheckRaw(TObjArray *listrec , TObjArray *listref) const
 205 {
 206
 207     //TH1 *fhRawEff;
 208     //TH1 *fhRawRef;
 209   TH2 *fhRawRec2d;
 210   TH2 *fhTime;
 211
 212   TIter next(listref) ;
 213     //Int_t counter=0;
 214   Float_t refmean[50][25];
 215   Float_t refrms[50][25];
 216   Float_t checkr = 0;
 217
 218     //Int_t nref = listref->GetEntries();
 219     //Int_t nrec = listrec->GetEntries();
 220
 221   for (Int_t iii=4; iii<6; iii++){
 222     fhRawRec2d =(TH2*) listref->At(iii);
 223     for (Int_t idet=1; idet<25; idet++) {
 224
 225       refmean[iii-4][idet] = fhRawRec2d->
 226         ProjectionY(Form("%s_py_%i_%i",                                                              fhRawRec2d ->GetName(), idet,iii-4),
 227                     idet,idet+1)->GetMean();
 228
 229       refrms[iii-4][idet] = fhRawRec2d->
 230         ProjectionY(Form("%s_py%i_%i",
 231                          fhRawRec2d ->GetName(), idet,iii-4),
 232                     idet,idet+1)->GetRMS();
 233
 234       }
 235   }
 236
 237
 238   TString nameDev[2] = {"CDF", "LED"};
 239   for (Int_t icase=208; icase<210; icase++) {
 240     fhTime = (TH2*) listrec->At(icase);
 241     for (Int_t idet=1; idet<25; idet++) {
 242       Double_t binmean = fhTime->
 243       ProjectionY(Form("%s_py_%i_%i",                                                              fhTime ->GetName(), idet,icase),
 244                     idet,idet+1)->GetMean();
 245 //        Double_t rms= fhTime ->ProjectionY(Form("%s_py%i_%i",
 246 //                                              fhTime ->GetName(), idet,icase),
 247 //                                         idet,idet+1)->GetRMS();
 248       Double_t diffmean = binmean-refmean[icase-208][idet];
 249
 250       if (TMath::Abs(diffmean) < 2 ) {
 251         checkr = 1;
 252         //      printf(" Laser calibration signal sits on its place %f for PMT %s %i : check = %f\n",  diffmean, nameDev[icase-208].Data() ,idet, checkr);
 253         AliDebug(AliQAv1::GetQADebugLevel(),
 254                  Form(" Laser calibration signal sits on its place %f for PMT %s %i : check = %f\n",  diffmean, nameDev[icase-208].Data(),idet, checkr)) ;
 255      }
 256      if (TMath::Abs(diffmean) <= 5 && TMath::Abs(diffmean) >= 2 ){
 257        checkr = 0.5;
 258        // printf(" Laser calibration signal shifted by  %f ps for PMT %s %i : check = %f\n",  diffmean*24.4, nameDev[icase-208].Data(),idet, checkr);
 259        AliDebug(AliQAv1::GetQADebugLevel(),
 260                 Form(" Laser calibration signal shifted by  %f ps (%f channels) for PMT %s %i : check = %f\n",  diffmean*24.4 ,diffmean , nameDev[icase-208].Data(),idet, checkr)) ;
 261       }
 262      if (TMath::Abs(diffmean) > 5) {
 263        checkr = 0.25;
 264        //   printf(" Big problems :laser calibration signal shifted by  %f ps (%f channels) for PMT %s %i : check = %f\n",  diffmean*24.4, diffmean, nameDev[icase-208].Data(),idet, checkr);
 265        AliDebug(AliQAv1::GetQADebugLevel(),
 266                 Form(" Big problems :laser calibration signal shifted by  %f ps (%f channels) for PMT %s %i : check = %i\n",  diffmean*24.4, diffmean, nameDev[icase-208].Data(),idet, checkr)) ;
 267
 268      }
 269
 270    }
 271  }
 272
 273   return checkr;
 274 }
 275