ITS/AliITSQADataMakerSim.cxx

   1 /**************************************************************************
   2  * Copyright(c) 2007-2009, 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 //  *************************************************************
  19 //  Checks the quality assurance
  20 //  by comparing with reference data
  21 //  contained in a DB
  22 //  -------------------------------------------------------------
  23 //  W. Ferrarese + P. Cerello Feb 2008
  24 //  INFN Torino
  25
  26 // --- ROOT system ---
  27 #include <TTree.h>
  28 #include <TMath.h>
  29 // --- Standard library ---
  30
  31 // --- AliRoot header files ---
  32 #include "AliITSQADataMakerSim.h"
  33 #include "AliITSQASPDDataMakerSim.h"
  34 #include "AliITSQASDDDataMakerSim.h"
  35 #include "AliITSQASSDDataMakerSim.h"
  36 #include "AliLog.h"
  37 #include "AliQAv1.h"
  38 #include "AliQAChecker.h"
  39 #include "AliITSQAChecker.h"
  40 #include "AliRawReader.h"
  41
  42 ClassImp(AliITSQADataMakerSim)
  43
  44 //____________________________________________________________________________
  45 AliITSQADataMakerSim::AliITSQADataMakerSim(Short_t subDet) :
  46 AliQADataMakerSim(AliQAv1::GetDetName(AliQAv1::kITS), "ITS Quality Assurance Data Maker"),
  47 fSubDetector(subDet),
  48 fSPDDataMaker(NULL),
  49 fSDDDataMaker(NULL),
  50 fSSDDataMaker(NULL)
  51 {
  52   //ctor used to discriminate OnLine-Offline analysis
  53   if(fSubDetector < 0 || fSubDetector > 3) {
  54         AliError("Error: fSubDetector number out of range; return\n");
  55   }
  56
  57   // Initialization for each subdetector
  58   if(fSubDetector == 0 || fSubDetector == 1) {
  59     AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM::Create SPD DataMakerSim\n");
  60         fSPDDataMaker = new AliITSQASPDDataMakerSim(this);
  61   }
  62   if(fSubDetector == 0 || fSubDetector == 2) {
  63     AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM::Create SDD DataMakerSim\n");
  64         //printf("AliITSQADM::Create SDD DataMakerSim\n");
  65         fSDDDataMaker = new AliITSQASDDDataMakerSim(this);
  66   }
  67   if(fSubDetector == 0 || fSubDetector == 3) {
  68     AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM::Create SSD DataMakerSim\n");
  69         fSSDDataMaker = new AliITSQASSDDataMakerSim(this);
  70   }
  71 }
  72
  73 //____________________________________________________________________________
  74 AliITSQADataMakerSim::~AliITSQADataMakerSim(){
  75   // destructor
  76   if(fSPDDataMaker)delete fSPDDataMaker;
  77   if(fSDDDataMaker)delete fSDDDataMaker;
  78   if(fSSDDataMaker)delete fSSDDataMaker;
  79 }
  80
  81 //____________________________________________________________________________
  82 AliITSQADataMakerSim::AliITSQADataMakerSim(const AliITSQADataMakerSim& qadm) :
  83 AliQADataMakerSim(),
  84 fSubDetector(qadm.fSubDetector),
  85 fSPDDataMaker(NULL),
  86 fSDDDataMaker(NULL),
  87 fSSDDataMaker(NULL)
  88 {
  89   //copy ctor
  90   SetName((const char*)qadm.GetName()) ;
  91   SetTitle((const char*)qadm.GetTitle());
  92 }
  93
  94 //__________________________________________________________________
  95 AliITSQADataMakerSim& AliITSQADataMakerSim::operator = (const AliITSQADataMakerSim& qac )
  96 {
  97   // Equal operator.
  98   this->~AliITSQADataMakerSim();
  99   new(this) AliITSQADataMakerSim(qac);
 100   return *this;
 101 }
 102
 103 //____________________________________________________________________________
 104 void AliITSQADataMakerSim::StartOfDetectorCycle()
 105 {
 106   //Detector specific actions at start of cycle
 107   AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM::Start of ITS Cycle\n");
 108   ResetEventTrigClasses(); // reset triggers list to select all histos
 109   ResetEvCountCycle();
 110   //
 111   if(fSubDetector == 0 || fSubDetector == 1) fSPDDataMaker->StartOfDetectorCycle();
 112   if(fSubDetector == 0 || fSubDetector == 2) fSDDDataMaker->StartOfDetectorCycle();
 113   if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->StartOfDetectorCycle();
 114 }
 115
 116 //____________________________________________________________________________
 117 void AliITSQADataMakerSim::EndOfDetectorCycle(AliQAv1::TASKINDEX_t task, TObjArray** list)
 118 {
 119   //
 120   AliInfo(Form("End of Dedetctor Cycle called for %s\n",AliQAv1::GetTaskName(task).Data() ));
 121   ResetEventTrigClasses(); // reset triggers list to select all histos
 122   //
 123   // launch the QA checking
 124   for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
 125     if (! IsValidEventSpecie(specie, list) ) continue;
 126     SetEventSpecie(AliRecoParam::ConvertIndex(specie));
 127     //
 128     Int_t idnumber=list[specie]->GetUniqueID();
 129     //printf("specie %s \t id number == %d\n",AliRecoParam::GetEventSpecieName(specie),idnumber);
 130     if(idnumber==0) {
 131       //AliInfo(Form("No check for %s\n",AliQAv1::GetTaskName(task).Data() ))
 132       continue;
 133     } //skip kDigitsR and not filled TobjArray specie
 134     else{
 135       AliDebug(AliQAv1::GetQADebugLevel(),"AliITSDM instantiates checker with Run(AliQAv1::kITS, task, list)\n");
 136       if(fSubDetector == 0 || fSubDetector == 1) fSPDDataMaker->EndOfDetectorCycle(task, list); //[specie]);
 137       if(fSubDetector == 0 || fSubDetector == 2) fSDDDataMaker->EndOfDetectorCycle(task, list); //[specie]);
 138       if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->EndOfDetectorCycle(task, list); //[specie]);
 139       //
 140       AliQAChecker *qac = AliQAChecker::Instance();
 141       AliITSQAChecker *qacb = (AliITSQAChecker *) qac->GetDetQAChecker(0);
 142       Int_t subdet=GetSubDet();
 143       qacb->SetSubDet(subdet);
 144
 145       if(subdet== 0 ){
 146         qacb->SetTaskOffset(fSPDDataMaker->GetOffset(task,specie),fSDDDataMaker->GetOffset(task,specie),fSSDDataMaker->GetOffset(task,specie)); //Setting the offset for the QAChecker list
 147         qacb->SetHisto(fSPDDataMaker->GetTaskHisto(task), fSDDDataMaker->GetTaskHisto(task), fSSDDataMaker->GetTaskHisto(task));
 148       }
 149       else
 150         if(subdet!=0){
 151           Int_t offset=GetDetTaskOffset(subdet, task);
 152           qacb->SetDetTaskOffset(subdet,offset);
 153           Int_t histo=GetDetTaskHisto(subdet, task);
 154           qacb->SetDetHisto(subdet,histo);
 155         }
 156       qac->Run( AliQAv1::kITS , task, list);
 157     }//end else unique id
 158   }//end for
 159 }
 160
 161 //____________________________________________________________________________
 162 void AliITSQADataMakerSim::InitDigits()
 163 {
 164
 165  // Initialization for Digits data
 166   fDigitsQAList[AliRecoParam::AConvert(fEventSpecie)]->SetUniqueID(60);
 167   if(fSubDetector == 0 || fSubDetector == 1) {
 168     AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SPD InitDigits\n");
 169
 170     fSPDDataMaker->InitDigits();
 171   }
 172   if(fSubDetector == 0 || fSubDetector == 2) {
 173     AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SDD InitDigits\n");
 174
 175     fSDDDataMaker->SetOffset(AliQAv1::kDIGITS, fDigitsQAList[AliRecoParam::AConvert(fEventSpecie)]->GetEntries(),AliRecoParam::AConvert(fEventSpecie));
 176     fSDDDataMaker->InitDigits();
 177   }
 178   if(fSubDetector == 0 || fSubDetector == 3) {
 179     AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SSD InitDigits\n");
 180
 181     fSSDDataMaker->SetOffset(AliQAv1::kDIGITS, fDigitsQAList[AliRecoParam::AConvert(fEventSpecie)]->GetEntries(),AliRecoParam::AConvert(fEventSpecie));
 182     fSSDDataMaker->InitDigits();
 183   }
 184   //
 185   ClonePerTrigClass(AliQAv1::kDIGITS); // this should be the last line
 186 }
 187
 188 //____________________________________________________________________________
 189 void AliITSQADataMakerSim::MakeDigits()
 190 {
 191   // Fill QA for digits
 192   if(fSubDetector == 0 || fSubDetector == 1) fSPDDataMaker->MakeDigits();
 193   if(fSubDetector == 0 || fSubDetector == 2) fSDDDataMaker->MakeDigits();
 194   if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeDigits();
 195   //
 196   IncEvCountCycleDigits();
 197   IncEvCountTotalDigits();
 198   //
 199 }
 200
 201 //____________________________________________________________________________
 202 void AliITSQADataMakerSim::MakeDigits(TTree * digits)
 203 {
 204   // Fill QA for digits
 205   if(fSubDetector == 0 || fSubDetector == 1) fSPDDataMaker->MakeDigits(digits);
 206   if(fSubDetector == 0 || fSubDetector == 2) fSDDDataMaker->MakeDigits(digits);
 207   if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeDigits(digits);
 208   //
 209   IncEvCountCycleDigits();
 210   IncEvCountTotalDigits();
 211   //
 212 }
 213
 214 //____________________________________________________________________________
 215 void AliITSQADataMakerSim::InitSDigits()
 216 {
 217   // Initialization for SDigits
 218   fSDigitsQAList[AliRecoParam::AConvert(fEventSpecie)]->SetUniqueID(70);
 219   if(fSubDetector == 0 || fSubDetector == 1) {
 220     AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SPD InitSDigits\n");
 221
 222     fSPDDataMaker->InitSDigits();
 223   }
 224   if(fSubDetector == 0 || fSubDetector == 2){
 225     AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SDD InitSDigits\n");
 226
 227     fSDDDataMaker->SetOffset(AliQAv1::kSDIGITS, fSDigitsQAList [AliRecoParam::AConvert(fEventSpecie)]->GetEntries(),AliRecoParam::AConvert(fEventSpecie));
 228     fSDDDataMaker->InitSDigits();
 229   }
 230   if(fSubDetector == 0 || fSubDetector == 3) {
 231     AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SSD InitSDigits\n");
 232
 233     fSSDDataMaker->SetOffset(AliQAv1::kSDIGITS, fSDigitsQAList [AliRecoParam::AConvert(fEventSpecie)]->GetEntries(),AliRecoParam::AConvert(fEventSpecie));
 234     fSSDDataMaker->InitSDigits();
 235   }
 236   //
 237   ClonePerTrigClass(AliQAv1::kSDIGITS); // this should be the last line
 238 }
 239
 240 //____________________________________________________________________________
 241 void AliITSQADataMakerSim::MakeSDigits()
 242 {
 243   // Fill QA for sdigits
 244   if(fSubDetector == 0 || fSubDetector == 1) fSPDDataMaker->MakeSDigits();
 245   if(fSubDetector == 0 || fSubDetector == 2) fSDDDataMaker->MakeSDigits();
 246   if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeSDigits();
 247   //
 248   IncEvCountCycleSDigits();
 249   IncEvCountTotalSDigits();
 250   //
 251 }
 252
 253 //____________________________________________________________________________
 254 void AliITSQADataMakerSim::MakeSDigits(TTree * sdigits)
 255 {
 256   // Fill QA for recpoints
 257   if(fSubDetector == 0 || fSubDetector == 1) fSPDDataMaker->MakeSDigits(sdigits);
 258   if(fSubDetector == 0 || fSubDetector == 2) fSDDDataMaker->MakeSDigits(sdigits);
 259   if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeSDigits(sdigits);
 260   //
 261   IncEvCountCycleSDigits();
 262   IncEvCountTotalSDigits();
 263   //
 264 }
 265
 266 //____________________________________________________________________________
 267 void AliITSQADataMakerSim::InitHits()
 268 {
 269   // Initialization for hits
 270   fHitsQAList[AliRecoParam::AConvert(fEventSpecie)]->SetUniqueID(50);
 271   if(fSubDetector == 0 || fSubDetector == 1) {
 272     AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SPD InitHits\n");
 273     fSPDDataMaker->InitHits();
 274   }
 275   if(fSubDetector == 0 || fSubDetector == 2) {
 276     AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SDD InitHits\n");
 277
 278     fSDDDataMaker->SetOffset(AliQAv1::kHITS, fHitsQAList[AliRecoParam::AConvert(fEventSpecie)]->GetEntries(),AliRecoParam::AConvert(fEventSpecie));
 279     fSDDDataMaker->InitHits();
 280   }
 281   if(fSubDetector == 0 || fSubDetector == 3) {
 282     AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SSD InitHits\n");
 283
 284     fSSDDataMaker->SetOffset(AliQAv1::kHITS, fHitsQAList[AliRecoParam::AConvert(fEventSpecie)]->GetEntries(),AliRecoParam::AConvert(fEventSpecie));
 285     fSSDDataMaker->InitHits();
 286   }
 287   //
 288   ClonePerTrigClass(AliQAv1::kHITS); // this should be the last line
 289 }
 290
 291 //____________________________________________________________________________
 292 void AliITSQADataMakerSim::MakeHits()
 293 {
 294   // Fill QA for hits
 295   if(fSubDetector == 0 || fSubDetector == 1) fSPDDataMaker->MakeHits();
 296   if(fSubDetector == 0 || fSubDetector == 2) fSDDDataMaker->MakeHits();
 297   if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeHits();
 298   //
 299   IncEvCountCycleHits();
 300   IncEvCountTotalHits();
 301   //
 302 }
 303
 304 //____________________________________________________________________________
 305 void AliITSQADataMakerSim::MakeHits(TTree * hits)
 306 {
 307   // Fill QA for hits
 308   if(fSubDetector == 0 || fSubDetector == 1) fSPDDataMaker->MakeHits(hits);
 309   if(fSubDetector == 0 || fSubDetector == 2) fSDDDataMaker->MakeHits(hits);
 310   if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeHits(hits);
 311   //
 312   IncEvCountCycleHits();
 313   IncEvCountTotalHits();
 314   //
 315 }
 316
 317 //_________________________________________________________________
 318 Int_t AliITSQADataMakerSim::GetDetTaskOffset(Int_t subdet,AliQAv1::TASKINDEX_t task)
 319 {
 320
 321   //return the offset for each subdetector
 322   switch(subdet)
 323     {
 324       Int_t offset;
 325     case 1:
 326       offset=fSPDDataMaker->GetOffset(task);
 327       return offset;
 328       break;
 329     case 2:
 330       offset=fSDDDataMaker->GetOffset(task);
 331       return offset;
 332       break;
 333     case 3:
 334       offset=fSSDDataMaker->GetOffset(task);
 335       return offset;
 336       break;
 337     default:
 338       AliWarning("No specific subdetector (SPD, SDD, SSD) selected!! Offset set to zero \n");
 339       offset=0;
 340       return offset;
 341       break;
 342     }
 343   //return offset;
 344 }
 345
 346 //_________________________________________________________________
 347 Int_t AliITSQADataMakerSim::GetDetTaskHisto(Int_t subdet,AliQAv1::TASKINDEX_t task)
 348 {
 349   //return of the number of histograms for each task and for each sub detector
 350   switch(subdet)
 351     {
 352
 353       Int_t histo;
 354     case 1:
 355       histo=fSPDDataMaker->GetOffset(task);
 356       return histo;
 357       break;
 358     case 2:
 359       histo=fSDDDataMaker->GetOffset(task);
 360       return histo;
 361       break;
 362     case 3:
 363       histo=fSSDDataMaker->GetOffset(task);
 364       return histo;
 365       break;
 366     default:
 367       AliWarning("No specific subdetector (SPD, SDD, SSD) selected!! Offset set to zero \n");
 368       histo=0;
 369       return histo;
 370       break;
 371     }
 372   //return offset;
 373 }