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