ITS/ITSsim/AliITSQASSDDataMakerSim.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 //  SSD QA part: P. Christakoglou - NIKHEF/UU
  26
  27 // --- ROOT system ---
  28 #include <TTree.h>
  29 #include <TH1.h>
  30 #include <TH2.h>
  31 #include <TMath.h>
  32 // --- Standard library ---
  33
  34 // --- AliRoot header files ---
  35 #include "AliITS.h"
  36 #include "AliITSmodule.h"
  37 #include "AliITShit.h"
  38 #include "AliITSdigit.h"
  39 #include "AliITSpListItem.h"
  40 #include "AliRun.h"
  41 #include "AliITSQADataMakerSim.h"
  42 #include "AliITSQASSDDataMakerSim.h"
  43 #include "AliLog.h"
  44 #include "AliQAv1.h"
  45 #include "AliQAChecker.h"
  46 #include "AliRawReader.h"
  47
  48 ClassImp(AliITSQASSDDataMakerSim)
  49
  50 //____________________________________________________________________________
  51 AliITSQASSDDataMakerSim::AliITSQASSDDataMakerSim(AliITSQADataMakerSim *aliITSQADataMakerSim) :
  52 TObject(),
  53 fAliITSQADataMakerSim(aliITSQADataMakerSim),
  54 //fSSDhTask(0),
  55 fSSDhHTask(0),
  56 fSSDhSTask(0),
  57 fSSDhDTask(0),
  58 fGenOffsetH(0),
  59 fGenOffsetS(0),
  60 fGenOffsetD(0)
  61 {
  62   //ctor used to discriminate OnLine-Offline analysis
  63   fGenOffsetH=  new Int_t[AliRecoParam::kNSpecies];
  64   fGenOffsetS=  new Int_t[AliRecoParam::kNSpecies];
  65   fGenOffsetD=  new Int_t[AliRecoParam::kNSpecies];
  66   for(Int_t i=0; i<AliRecoParam::kNSpecies; i++)
  67     {
  68       fGenOffsetH[i]= 0;
  69       fGenOffsetS[i]= 0;
  70       fGenOffsetD[i]= 0;
  71     }
  72 }
  73
  74 //____________________________________________________________________________
  75 AliITSQASSDDataMakerSim::AliITSQASSDDataMakerSim(const AliITSQASSDDataMakerSim& qadm) :
  76 TObject(),
  77 fAliITSQADataMakerSim(qadm.fAliITSQADataMakerSim),
  78 //fSSDhTask(qadm.fSSDhTask),
  79 fSSDhHTask(qadm.fSSDhHTask),
  80 fSSDhSTask(qadm.fSSDhSTask),
  81 fSSDhDTask(qadm.fSSDhDTask),
  82 fGenOffsetH(qadm.fGenOffsetH),
  83 fGenOffsetS(qadm.fGenOffsetS),
  84 fGenOffsetD(qadm.fGenOffsetD)
  85 {
  86   //copy ctor
  87   fAliITSQADataMakerSim->SetName((const char*)qadm.fAliITSQADataMakerSim->GetName()) ;
  88   fAliITSQADataMakerSim->SetTitle((const char*)qadm.fAliITSQADataMakerSim->GetTitle());
  89   }
  90
  91 //__________________________________________________________________
  92 AliITSQASSDDataMakerSim& AliITSQASSDDataMakerSim::operator = (const AliITSQASSDDataMakerSim& qac ) {
  93   // Equal operator.
  94   this->~AliITSQASSDDataMakerSim();
  95   new(this) AliITSQASSDDataMakerSim(qac);
  96   return *this;
  97 }
  98
  99 //____________________________________________________________________________
 100 void AliITSQASSDDataMakerSim::StartOfDetectorCycle() {
 101   //Detector specific actions at start of cycle
 102   AliDebug(AliQAv1::GetQADebugLevel(),"AliITSQADM::Start of SSD Cycle\n");
 103 }
 104
 105 //____________________________________________________________________________
 106 void AliITSQASSDDataMakerSim::EndOfDetectorCycle(AliQAv1::TASKINDEX_t /*task*/, TObjArray** /*list*/) {
 107   // launch the QA checking
 108   AliDebug(AliQAv1::GetQADebugLevel(),"AliITSDM instantiates checker with Run(AliQAv1::kITS, task, list)\n");
 109
 110 //  AliQAChecker::Instance()->Run( AliQAv1::kITS , task, list);
 111 }
 112
 113 //____________________________________________________________________________
 114 Int_t AliITSQASSDDataMakerSim::InitDigits() {
 115   // Initialization for DIGIT data - SSD -
 116   const Bool_t expert   = kTRUE ;
 117   const Bool_t image    = kTRUE ;
 118   Int_t rv = 0 ;
 119  // fGenOffsetD = (fAliITSQADataMakerSim->fDigitsQAList[AliRecoParam::kDefault])->GetEntries();
 120
 121   // custom code here
 122   TH1F *fHistSSDModule = new TH1F("fHistSSDDigitsModule",
 123                                   "SSD Digits Module;SSD Module Number;N_{DIGITS}",
 124                                   1698,499.5,2197.5);
 125   rv = fAliITSQADataMakerSim->Add2DigitsList(fHistSSDModule,
 126                                         fGenOffsetD[fAliITSQADataMakerSim->GetEventSpecie()] + 0, !expert, image);
 127   fSSDhDTask += 1;
 128   TH2F *fHistSSDModuleStrip = new TH2F("fHistSSDDigitsModuleStrip",
 129                                        "SSD Digits Module Strip;N_{Strip};N_{Module}",
 130                                        1540,0,1540,1698,499.5,2197.5);
 131   rv = fAliITSQADataMakerSim->Add2DigitsList(fHistSSDModuleStrip,
 132                                         fGenOffsetD[fAliITSQADataMakerSim->GetEventSpecie()] + 1, !expert, image);
 133   fSSDhDTask += 1;
 134
 135   AliDebug(AliQAv1::GetQADebugLevel(),Form("%d SSD Digits histograms booked\n",fSSDhDTask));
 136   return rv ;
 137 }
 138
 139 //____________________________________________________________________________
 140 Int_t AliITSQASSDDataMakerSim::MakeDigits(TTree *digits)
 141 {
 142   // Fill QA for DIGIT - SSD -
 143   Int_t rv = 0 ;
 144
 145   AliITS *fITS  = (AliITS*)gAlice->GetModule("ITS");
 146   fITS->SetTreeAddress();
 147   TClonesArray *iSSDdigits  = fITS->DigitsAddress(2);
 148   for(Int_t iModule = 500; iModule < 2198; iModule++) {
 149     iSSDdigits->Clear();
 150     digits->GetEvent(iModule);
 151     Int_t ndigits = iSSDdigits->GetEntries();
 152     fAliITSQADataMakerSim->FillDigitsData(fGenOffsetD[fAliITSQADataMakerSim->GetEventSpecie()] + 0,iModule,ndigits);
 153     if(ndigits != 0)
 154       AliDebug(AliQAv1::GetQADebugLevel(),Form("Module: %d - Digits: %d",iModule,ndigits));
 155
 156     for (Int_t iDigit = 0; iDigit < ndigits; iDigit++) {
 157       AliITSdigit *dig = (AliITSdigit*)iSSDdigits->UncheckedAt(iDigit);
 158       Int_t fStripNumber = (dig->GetCoord1() == 0) ? dig->GetCoord2() : dig->GetCoord2() + fgkNumberOfPSideStrips;
 159       fAliITSQADataMakerSim->FillDigitsData(fGenOffsetD[fAliITSQADataMakerSim->GetEventSpecie()] + 1,fStripNumber,iModule,dig->GetSignal());
 160     }//digit loop
 161   }//module loop
 162   //
 163   return rv ;
 164 }
 165
 166 //____________________________________________________________________________
 167 Int_t AliITSQASSDDataMakerSim::InitSDigits() {
 168   // Initialization for SDIGIT data - SSD -
 169   const Bool_t expert   = kTRUE ;
 170   const Bool_t image    = kTRUE ;
 171   Int_t rv = 0 ;
 172   //fGenOffsetS = (fAliITSQADataMakerSim->fSDigitsQAList[AliRecoParam::kDefault])->GetEntries();
 173
 174   // custom code here
 175   TH1F *fHistSSDModule = new TH1F("fHistSSDSDigitsModule",
 176                                   "SSD SDigits Module;SSD Module Number;N_{SDIGITS}",
 177                                   1698,499.5,2197.5);
 178   rv = fAliITSQADataMakerSim->Add2SDigitsList(fHistSSDModule,
 179                                          fGenOffsetS[fAliITSQADataMakerSim->GetEventSpecie()] + 0, !expert, image);
 180   fSSDhSTask += 1;
 181
 182   AliDebug(AliQAv1::GetQADebugLevel(),Form("%d SSD SDigits histograms booked\n",fSSDhSTask));
 183   //
 184   return rv ;
 185 }
 186
 187 //____________________________________________________________________________
 188 Int_t AliITSQASSDDataMakerSim::MakeSDigits(TTree *sdigits) {
 189   // Fill QA for SDIGIT - SSD -
 190   Int_t rv = 0 ;
 191
 192   static TClonesArray iSSDEmpty("AliITSpListItem",10000);
 193   iSSDEmpty.Clear();
 194   TClonesArray *iSSDsdigits = &iSSDEmpty;
 195
 196   //  AliDebug(AliQAv1::GetQADebugLevel(), Form("Trying to access the sdigits histogram: %d\n",fGenOffsetS));
 197
 198   TBranch *brchSDigits = sdigits->GetBranch("ITS");
 199   brchSDigits->SetAddress(&iSSDsdigits);
 200   for(Int_t iModule = 500; iModule < 2198; iModule++) {
 201     iSSDsdigits->Clear();
 202     sdigits->GetEvent(iModule);
 203     Int_t ndigits = iSSDsdigits->GetEntries();
 204     fAliITSQADataMakerSim->FillSDigitsData(fGenOffsetS[fAliITSQADataMakerSim->GetEventSpecie()] + 0,iModule,ndigits);
 205     if(ndigits != 0)
 206       AliDebug(AliQAv1::GetQADebugLevel(),Form("Module: %d - Digits: %d",iModule,ndigits));
 207
 208 //     for (Int_t iDigit = 0; iDigit < ndigits; iDigit++) {
 209 //       AliITSpListItem *dig=(AliITSpListItem*)iSSDsdigits->At(iDigit);
 210 //       dig=0;
 211 //     }//digit loop
 212   }//module loop
 213   //
 214   return rv ;
 215 }
 216
 217 //____________________________________________________________________________
 218 Int_t AliITSQASSDDataMakerSim::InitHits() {
 219   // Initialization for HITS data - SSD -
 220   const Bool_t expert   = kTRUE ;
 221   const Bool_t image    = kTRUE ;
 222   Int_t rv = 0 ;
 223
 224   //fGenOffsetH = (fAliITSQADataMakerSim->fHitsQAList[fEventSpecie])->GetEntries();
 225
 226   // custom code here
 227   TH1F *fHistSSDModule = new TH1F("fHistSSDHitsModule",
 228                                   "SSD Hits Module;SSD Module Number;N_{HITS}",
 229                                   1698,499.5,2197.5);
 230   rv = fAliITSQADataMakerSim->Add2HitsList(fHistSSDModule,
 231                                       fGenOffsetH[fAliITSQADataMakerSim->GetEventSpecie()] + 0, !expert, image);
 232   fSSDhHTask += 1;
 233   TH1F *fHistSSDGlobalX = new TH1F("fHistSSDHitsGlobalX",
 234                                    "SSD Hits Global X;x [cm];Entries",
 235                                    1000,-50.,50.);
 236   rv = fAliITSQADataMakerSim->Add2HitsList(fHistSSDGlobalX,
 237                                       fGenOffsetH[fAliITSQADataMakerSim->GetEventSpecie()] + 1, !expert, image);
 238   fSSDhHTask += 1;
 239   TH1F *fHistSSDGlobalY = new TH1F("fHistSSDHitsGlobalY",
 240                                    "SSD Hits Global Y;y [cm];Entries",
 241                                    1000,-50.,50.);
 242   rv = fAliITSQADataMakerSim->Add2HitsList(fHistSSDGlobalY,
 243                                       fGenOffsetH[fAliITSQADataMakerSim->GetEventSpecie()] + 2, !expert, image);
 244   fSSDhHTask += 1;
 245   TH1F *fHistSSDGlobalZ = new TH1F("fHistSSDHitsGlobalZ",
 246                                    "SSD Hits Global Z ;z [cm];Entries",
 247                                    1000,-60.,60.);
 248   rv = fAliITSQADataMakerSim->Add2HitsList(fHistSSDGlobalZ,
 249                                       fGenOffsetH[fAliITSQADataMakerSim->GetEventSpecie()] + 3, !expert, image);
 250   fSSDhHTask += 1;
 251   TH1F *fHistSSDLocalX = new TH1F("fHistSSDHitsLocalX",
 252                                   "SSD Hits Local X;x [cm];Entries",
 253                                   1000,-4.,4.);
 254   rv = fAliITSQADataMakerSim->Add2HitsList(fHistSSDLocalX,
 255                                       fGenOffsetH[fAliITSQADataMakerSim->GetEventSpecie()] + 4, !expert, image);
 256   fSSDhHTask += 1;
 257   TH1F *fHistSSDLocalY = new TH1F("fHistSSDHitsLocalY",
 258                                   "SSD Hits Local Y;y [cm];Entries",
 259                                   1000,-0.1,0.1);
 260   rv = fAliITSQADataMakerSim->Add2HitsList(fHistSSDLocalY,
 261                                       fGenOffsetH[fAliITSQADataMakerSim->GetEventSpecie()] + 5, !expert, image);
 262   fSSDhHTask += 1;
 263   TH1F *fHistSSDLocalZ = new TH1F("fHistSSDHitsLocalZ",
 264                                   "SSD Hits Local Z;z [cm];Entries",
 265                                   1000,-4.,4.);
 266   rv = fAliITSQADataMakerSim->Add2HitsList(fHistSSDLocalZ,
 267                                       fGenOffsetH[fAliITSQADataMakerSim->GetEventSpecie()] + 6, !expert, image);
 268   fSSDhHTask += 1;
 269   TH1F *fHistSSDIonization = new TH1F("fHistSSDHitsIonization",
 270                                       "SSD Hits Ionization;log(dE/dx) [KeV];N_{Hits}",
 271                                       100,-7,-2);
 272   rv = fAliITSQADataMakerSim->Add2HitsList(fHistSSDIonization,
 273                                       fGenOffsetH[fAliITSQADataMakerSim->GetEventSpecie()] + 7, !expert, image);
 274   fSSDhHTask += 1;
 275   TH2F *fHistSSDGlobalXY = new TH2F("fHistSSDHitsGlobalXY",
 276                                     "SSD Hits Global XY;x [cm];y [cm]",
 277                                     1000,-50.,50.,
 278                                     1000,-50.,50.);
 279   rv = fAliITSQADataMakerSim->Add2HitsList(fHistSSDGlobalXY,
 280                                       fGenOffsetH[fAliITSQADataMakerSim->GetEventSpecie()] + 8, !expert, image);
 281   fSSDhHTask += 1;
 282
 283   AliDebug(AliQAv1::GetQADebugLevel(),Form("%d SSD Hits histograms booked\n",fSSDhHTask));
 284   return rv ;
 285 }
 286
 287
 288 //____________________________________________________________________________
 289 Int_t AliITSQASSDDataMakerSim::MakeHits(TTree *hits) {
 290   // Fill QA for HITS - SSD -
 291   Int_t rv = 0 ;
 292
 293   AliITS *fITS  = (AliITS*)gAlice->GetModule("ITS");
 294   fITS->SetTreeAddress();
 295   Int_t nmodules;
 296   fITS->InitModules(-1,nmodules);
 297   fITS->FillModules(hits,0);
 298   for(Int_t iModule = 500; iModule < 2198; iModule++) {
 299     AliITSmodule *module = fITS->GetModule(iModule);
 300     TObjArray *arrHits = module->GetHits();
 301     Int_t nhits = arrHits->GetEntriesFast();
 302     if(nhits != 0)
 303       AliDebug(AliQAv1::GetQADebugLevel(),Form("Module: %d - Hits: %d",iModule,nhits));
 304     for (Int_t iHit = 0; iHit < nhits; iHit++) {
 305       AliITShit *hit = (AliITShit*) arrHits->At(iHit);
 306
 307       fAliITSQADataMakerSim->FillHitsData(fGenOffsetH[fAliITSQADataMakerSim->GetEventSpecie()] + 0,iModule);
 308       fAliITSQADataMakerSim->FillHitsData(fGenOffsetH[fAliITSQADataMakerSim->GetEventSpecie()] + 1,hit->GetXG());
 309       fAliITSQADataMakerSim->FillHitsData(fGenOffsetH[fAliITSQADataMakerSim->GetEventSpecie()] + 2,hit->GetYG());
 310       fAliITSQADataMakerSim->FillHitsData(fGenOffsetH[fAliITSQADataMakerSim->GetEventSpecie()] + 3,hit->GetZG());
 311       fAliITSQADataMakerSim->FillHitsData(fGenOffsetH[fAliITSQADataMakerSim->GetEventSpecie()] + 4,hit->GetXL());
 312       fAliITSQADataMakerSim->FillHitsData(fGenOffsetH[fAliITSQADataMakerSim->GetEventSpecie()] + 5,hit->GetYL());
 313       fAliITSQADataMakerSim->FillHitsData(fGenOffsetH[fAliITSQADataMakerSim->GetEventSpecie()] + 6,hit->GetZL());
 314       if(hit->GetIonization())
 315         fAliITSQADataMakerSim->FillHitsData(fGenOffsetH[fAliITSQADataMakerSim->GetEventSpecie()] + 7,TMath::Log10(hit->GetIonization()));
 316       fAliITSQADataMakerSim->FillHitsData(fGenOffsetH[fAliITSQADataMakerSim->GetEventSpecie()] + 8,hit->GetXG(),hit->GetYG());
 317     }//hit loop
 318   }//module loop
 319   return rv ;
 320 }
 321
 322 //____________________________________________________________________________
 323 Int_t AliITSQASSDDataMakerSim::GetOffset(AliQAv1::TASKINDEX_t task,Int_t specie){
 324   // Returns histogram offset according to the specified task
 325   Int_t offset=0;
 326   if( task == AliQAv1::kHITS){
 327     offset=fGenOffsetH[specie];
 328   }
 329   else if( task == AliQAv1::kSDIGITS) {
 330     offset=fGenOffsetS[specie];
 331   }
 332   else if( task == AliQAv1::kDIGITS) {
 333     offset=fGenOffsetD[specie];
 334   }
 335   else {
 336     AliInfo("No task has been selected. TaskHisto set to zero.\n");
 337   }
 338
 339   return offset;
 340 }
 341
 342
 343 //____________________________________________________________________________
 344 void AliITSQASSDDataMakerSim::SetOffset(AliQAv1::TASKINDEX_t task, Int_t offset,Int_t specie ){
 345   // Returns histogram offset according to the specified task
 346   if( task == AliQAv1::kHITS){
 347     fGenOffsetH[specie] = offset;
 348   }
 349   else if( task == AliQAv1::kSDIGITS) {
 350     fGenOffsetS[specie] = offset;
 351   }
 352   else if( task == AliQAv1::kDIGITS) {
 353     fGenOffsetD[specie] = offset;
 354   }
 355   else {
 356     AliInfo("No task has been selected. TaskHisto set to zero.\n");
 357   }
 358 }
 359
 360 //____________________________________________________________________________
 361 Int_t AliITSQASSDDataMakerSim::GetTaskHisto(AliQAv1::TASKINDEX_t task) {
 362   // Returns the number of booked histograms for the selected task
 363   Int_t histotot=0;
 364   if( task == AliQAv1::kHITS) {
 365     histotot=fSSDhHTask ;
 366   }
 367   else if( task == AliQAv1::kSDIGITS) {
 368     histotot=fSSDhSTask;
 369   }
 370   else if( task == AliQAv1::kDIGITS) {
 371     histotot=fSSDhDTask ;
 372   }
 373   else {
 374     AliInfo("No task has been selected. TaskHisto set to zero.\n");
 375   }
 376   return histotot;
 377
 378 }