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 RAW data
  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
 109   if(fSubDetector == 0 || fSubDetector == 1) fSPDDataMaker->StartOfDetectorCycle();
 110   if(fSubDetector == 0 || fSubDetector == 2) fSDDDataMaker->StartOfDetectorCycle();
 111   if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->StartOfDetectorCycle();
 112 }
 113
 114 //____________________________________________________________________________
 115 void AliITSQADataMakerSim::EndOfDetectorCycle(AliQAv1::TASKINDEX_t task, TObjArray** list)
 116 {
 117   // launch the QA checking
 118   for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
 119     AliDebug(AliQAv1::GetQADebugLevel(),"AliITSDM instantiates checker with Run(AliQAv1::kITS, task, list)\n");
 120     if(fSubDetector == 0 || fSubDetector == 1) fSPDDataMaker->EndOfDetectorCycle(task, list[specie]);
 121     if(fSubDetector == 0 || fSubDetector == 2) fSDDDataMaker->EndOfDetectorCycle(task, list[specie]);
 122     if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->EndOfDetectorCycle(task, list[specie]);
 123
 124
 125     AliQAChecker *qac = AliQAChecker::Instance();
 126     AliITSQAChecker *qacb = (AliITSQAChecker *) qac->GetDetQAChecker(0);
 127     Int_t subdet=GetSubDet();
 128     qacb->SetSubDet(subdet);
 129
 130     if(subdet== 0 ){
 131       qacb->SetTaskOffset(fSPDDataMaker->GetOffset(task),fSDDDataMaker->GetOffset(task),fSSDDataMaker->GetOffset(task)); //Setting the offset for the QAChecker list
 132     }
 133     else
 134       if(subdet!=0){
 135         Int_t offset=GetDetTaskOffset(subdet, task);
 136         qacb->SetDetTaskOffset(subdet,offset);
 137       }
 138     qac->Run( AliQAv1::kITS , task, list);
 139   }
 140 }
 141
 142 //____________________________________________________________________________
 143 void AliITSQADataMakerSim::InitDigits()
 144 {
 145   // Initialization for RAW data
 146         if(fSubDetector == 0 || fSubDetector == 1) {
 147           AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SPD InitDigits\n");
 148           fSPDDataMaker->InitDigits();
 149         }
 150         if(fSubDetector == 0 || fSubDetector == 2) {
 151           AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SDD InitDigits\n");
 152           fSDDDataMaker->InitDigits();
 153         }
 154         if(fSubDetector == 0 || fSubDetector == 3) {
 155           AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SSD InitDigits\n");
 156           fSSDDataMaker->InitDigits();
 157         }
 158 }
 159
 160 //____________________________________________________________________________
 161 void AliITSQADataMakerSim::MakeDigits(TClonesArray * digits)
 162 {
 163   // Fill QA for digits
 164   if(fSubDetector == 0 || fSubDetector == 1) {
 165     Int_t rv = fSPDDataMaker->MakeDigits(digits) ;
 166     if ( rv != 0 )
 167       fSDDDataMaker->SetOffset(AliQAv1::kDIGITS, fDigitsQAList [AliRecoParam::AConvert(fEventSpecie)]->GetEntries());
 168   }
 169
 170   if(fSubDetector == 0 || fSubDetector == 2) {
 171     Int_t rv = fSDDDataMaker->MakeDigits(digits) ;
 172     if ( rv != 0 )
 173       fSSDDataMaker->SetOffset(AliQAv1::kDIGITS, fDigitsQAList [AliRecoParam::AConvert(fEventSpecie)]->GetEntries());
 174   }
 175
 176   if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeDigits(digits);
 177 }
 178
 179 //____________________________________________________________________________
 180 void AliITSQADataMakerSim::MakeDigits(TTree * digits)
 181 {
 182   // Fill QA for digits
 183   if(fSubDetector == 0 || fSubDetector == 1) {
 184     Int_t rv = fSPDDataMaker->MakeDigits(digits) ;
 185     if ( rv != 0 )
 186       fSDDDataMaker->SetOffset(AliQAv1::kDIGITS, fDigitsQAList [AliRecoParam::AConvert(fEventSpecie)]->GetEntries());
 187   }
 188
 189   if(fSubDetector == 0 || fSubDetector == 2){
 190     Int_t rv = fSDDDataMaker->MakeDigits(digits) ;
 191     if ( rv != 0 )
 192       fSSDDataMaker->SetOffset(AliQAv1::kDIGITS, fDigitsQAList [AliRecoParam::AConvert(fEventSpecie)]->GetEntries());
 193   }
 194
 195   if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeDigits(digits);
 196 }
 197
 198 //____________________________________________________________________________
 199 void AliITSQADataMakerSim::InitSDigits()
 200 {
 201   // Initialization for RECPOINTS
 202   if(fSubDetector == 0 || fSubDetector == 1) {
 203     AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SPD InitSDigits\n");
 204     fSPDDataMaker->InitSDigits();
 205   }
 206   if(fSubDetector == 0 || fSubDetector == 2) {
 207     AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SDD InitSDigits\n");
 208         fSDDDataMaker->InitSDigits();
 209   }
 210   if(fSubDetector == 0 || fSubDetector == 3) {
 211     AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SSD InitSDigits\n");
 212         fSSDDataMaker->InitSDigits();
 213   }
 214 }
 215
 216 //____________________________________________________________________________
 217 void AliITSQADataMakerSim::MakeSDigits(TClonesArray * sdigits)
 218 {
 219   // Fill QA for recpoints
 220   if(fSubDetector == 0 || fSubDetector == 1){
 221     Int_t rv = fSPDDataMaker->MakeSDigits(sdigits) ;
 222     if ( rv != 0 )
 223       fSDDDataMaker->SetOffset(AliQAv1::kSDIGITS, fSDigitsQAList [AliRecoParam::AConvert(fEventSpecie)]->GetEntries());
 224   }
 225
 226   if(fSubDetector == 0 || fSubDetector == 2) {
 227     Int_t rv = fSDDDataMaker->MakeSDigits(sdigits) ;
 228     if ( rv != 0 )
 229       fSSDDataMaker->SetOffset(AliQAv1::kSDIGITS, fSDigitsQAList [AliRecoParam::AConvert(fEventSpecie)]->GetEntries());
 230   }
 231
 232   if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeSDigits(sdigits);
 233 }
 234
 235 //____________________________________________________________________________
 236 void AliITSQADataMakerSim::MakeSDigits(TTree * sdigits)
 237 {
 238   // Fill QA for recpoints
 239   if(fSubDetector == 0 || fSubDetector == 1){
 240     Int_t rv = fSPDDataMaker->MakeSDigits(sdigits) ;
 241     if ( rv != 0 )
 242       fSDDDataMaker->SetOffset(AliQAv1::kSDIGITS, fSDigitsQAList [AliRecoParam::AConvert(fEventSpecie)]->GetEntries());
 243   }
 244
 245   if(fSubDetector == 0 || fSubDetector == 2){
 246     Int_t rv = fSDDDataMaker->MakeSDigits(sdigits) ;
 247     if ( rv != 0 )
 248       fSSDDataMaker->SetOffset(AliQAv1::kSDIGITS, fSDigitsQAList [AliRecoParam::AConvert(fEventSpecie)]->GetEntries());
 249   }
 250
 251   if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeSDigits(sdigits);
 252 }
 253
 254 //____________________________________________________________________________
 255 void AliITSQADataMakerSim::InitHits()
 256 {
 257   // Initialization for hits
 258   if(fSubDetector == 0 || fSubDetector == 1) {
 259     AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SPD InitHits\n");
 260     fSPDDataMaker->InitHits();
 261   }
 262   if(fSubDetector == 0 || fSubDetector == 2) {
 263     AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SDD InitHits\n");
 264         fSDDDataMaker->InitHits();
 265   }
 266   if(fSubDetector == 0 || fSubDetector == 3) {
 267     AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM:: SSD InitHits\n");
 268         fSSDDataMaker->InitHits();
 269   }
 270 }
 271
 272 //____________________________________________________________________________
 273 void AliITSQADataMakerSim::MakeHits(TClonesArray * hits)
 274 {
 275   // Fill QA for hits
 276   if(fSubDetector == 0 || fSubDetector == 1) {
 277     Int_t rv = fSPDDataMaker->MakeHits(hits) ;
 278     if ( rv != 0 )
 279       fSDDDataMaker->SetOffset(AliQAv1::kHITS, fHitsQAList [AliRecoParam::AConvert(fEventSpecie)]->GetEntries());
 280   }
 281
 282   if(fSubDetector == 0 || fSubDetector == 2) {
 283     Int_t rv = fSDDDataMaker->MakeHits(hits) ;
 284     if ( rv != 0 )
 285       fSSDDataMaker->SetOffset(AliQAv1::kHITS, fHitsQAList [AliRecoParam::AConvert(fEventSpecie)]->GetEntries());
 286   }
 287
 288   if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeHits(hits);
 289 }
 290
 291 //____________________________________________________________________________
 292 void AliITSQADataMakerSim::MakeHits(TTree * hits)
 293 {
 294   // Fill QA for hits
 295   if(fSubDetector == 0 || fSubDetector == 1) {
 296     Int_t rv = fSPDDataMaker->MakeHits(hits) ;
 297     if ( rv != 0 )
 298       fSDDDataMaker->SetOffset(AliQAv1::kHITS, fHitsQAList [AliRecoParam::AConvert(fEventSpecie)]->GetEntries());
 299   }
 300   if(fSubDetector == 0 || fSubDetector == 2) {
 301     Int_t rv = fSDDDataMaker->MakeHits(hits) ;
 302     if ( rv != 0 )
 303       fSSDDataMaker->SetOffset(AliQAv1::kHITS, fHitsQAList [AliRecoParam::AConvert(fEventSpecie)]->GetEntries());
 304   }
 305
 306   if(fSubDetector == 0 || fSubDetector == 3) fSSDDataMaker->MakeHits(hits);
 307 }
 308
 309 //_________________________________________________________________
 310 Int_t AliITSQADataMakerSim::GetDetTaskOffset(Int_t subdet,AliQAv1::TASKINDEX_t task)
 311 {
 312   switch(subdet)
 313     {
 314
 315       Int_t offset;
 316     case 1:
 317       offset=fSPDDataMaker->GetOffset(task);
 318       return offset;
 319       break;
 320     case 2:
 321       offset=fSDDDataMaker->GetOffset(task);
 322       return offset;
 323       break;
 324     case 3:
 325       offset=fSSDDataMaker->GetOffset(task);
 326       return offset;
 327       break;
 328     default:
 329       AliWarning("No specific subdetector (SPD, SDD, SSD) selected!! Offset set to zero \n");
 330       offset=0;
 331       return offset;
 332       break;
 333     }
 334   //return offset;
 335 }