VZERO/AliVZEROQADataMakerSim.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
  17 /* $Id: AliVZEROQADataMakerSim.cxx 23123 2007-12-18 09:08:18Z hristov $ */
  18
  19 //---
  20 //  Produces the data needed to calculate the quality assurance.
  21 //  All data must be mergeable objects.
  22 //  Author : BC
  23 //---
  24
  25 // --- ROOT system ---
  26 #include <TClonesArray.h>
  27 #include <TFile.h>
  28 #include <TH1F.h>
  29 #include <TDirectory.h>
  30 // --- Standard library ---
  31
  32 // --- AliRoot header files ---
  33 #include "AliESDEvent.h"
  34 #include "AliLog.h"
  35 #include "AliVZEROdigit.h"
  36 #include "AliVZEROhit.h"
  37 #include "AliVZEROQADataMakerSim.h"
  38 #include "AliQAChecker.h"
  39
  40 ClassImp(AliVZEROQADataMakerSim)
  41
  42 //____________________________________________________________________________
  43   AliVZEROQADataMakerSim::AliVZEROQADataMakerSim() :
  44   AliQADataMakerSim(AliQAv1::GetDetName(AliQAv1::kVZERO), "VZERO Quality Assurance Data Maker")
  45
  46 {
  47   // constructor
  48
  49
  50 }
  51
  52 //____________________________________________________________________________
  53 AliVZEROQADataMakerSim::AliVZEROQADataMakerSim(const AliVZEROQADataMakerSim& qadm) :
  54   AliQADataMakerSim()
  55 {
  56   //copy constructor
  57
  58   SetName((const char*)qadm.GetName()) ;
  59   SetTitle((const char*)qadm.GetTitle());
  60 }
  61
  62 //__________________________________________________________________
  63 AliVZEROQADataMakerSim& AliVZEROQADataMakerSim::operator = (const AliVZEROQADataMakerSim& qadm )
  64 {
  65   // Assign operator.
  66   this->~AliVZEROQADataMakerSim();
  67   new(this) AliVZEROQADataMakerSim(qadm);
  68   return *this;
  69 }
  70 //____________________________________________________________________________
  71 void AliVZEROQADataMakerSim::EndOfDetectorCycle(AliQAv1::TASKINDEX_t task, TObjArray ** list)
  72 {
  73   //Detector specific actions at end of cycle
  74   // do the QA checking
  75   AliQAChecker::Instance()->Run(AliQAv1::kVZERO, task, list) ;
  76 }
  77
  78
  79 //____________________________________________________________________________
  80 void AliVZEROQADataMakerSim::InitHits()
  81 {
  82
  83   // create Hits histograms in Hits subdir
  84   const Bool_t expert   = kTRUE ;
  85   const Bool_t image    = kTRUE ;
  86
  87   TH1I * h0 = new TH1I("hHitMultiplicity", "Hit multiplicity distribution in VZERO;# of Hits;Entries", 300, 0, 299) ;
  88   h0->Sumw2() ;
  89   Add2HitsList(h0, 0, !expert, image) ;
  90
  91   TH1I * h1 = new TH1I("hHitCellNumber", "Hit cell distribution in VZERO;# of Hits;Entries", 80, 0, 79) ;
  92   h1->Sumw2() ;
  93   Add2HitsList(h1, 1, !expert, image) ;
  94
  95 }
  96
  97 //____________________________________________________________________________
  98 void AliVZEROQADataMakerSim::InitDigits()
  99 {
 100   // create Digits histograms in Digits subdir
 101   const Bool_t expert   = kTRUE ;
 102   const Bool_t image    = kTRUE ;
 103
 104   char TDCname[100];
 105   char ADCname[100];
 106   TH1I *fhDigTDC[64];
 107   TH1I *fhDigADC[64];
 108   char texte[100];
 109
 110   // create Digits histograms in Digits subdir
 111   TH1I * h0 = new TH1I("hDigitMultiplicity", "Digits multiplicity distribution in VZERO;# of Digits;Entries", 100, 0, 99) ;
 112   h0->Sumw2() ;
 113   Add2DigitsList(h0, 0, !expert, image) ;
 114
 115   for (Int_t i=0; i<64; i++)
 116     {
 117        sprintf(TDCname, "hDigitTDC%d", i);
 118        sprintf(texte,"Digit TDC in cell %d; TDC value;Entries",i);
 119        fhDigTDC[i] = new TH1I(TDCname,texte,300,0.,149.);
 120
 121        sprintf(ADCname,"hDigitADC%d",i);
 122        sprintf(texte,"Digit ADC in cell %d;ADC value;Entries",i);
 123        fhDigADC[i]= new TH1I(ADCname,texte,1024,0.,1023.);
 124
 125        Add2DigitsList(fhDigTDC[i],i+1, !expert, image);
 126        Add2DigitsList(fhDigADC[i],i+1+64, !expert, image);
 127      }
 128 }
 129
 130
 131 //____________________________________________________________________________
 132 void AliVZEROQADataMakerSim::MakeHits()
 133 {
 134         //make QA data from Hits
 135
 136   Int_t nhits = fHitsArray->GetEntriesFast();
 137   GetHitsData(0)->Fill(nhits) ;    // fills Hit multiplicity
 138     for (Int_t ihit=0;ihit<nhits;ihit++)
 139         {
 140            AliVZEROhit  * VZEROHit   = (AliVZEROhit*) fHitsArray->UncheckedAt(ihit);
 141            if (!VZEROHit) {
 142               AliError("The unchecked hit doesn't exist");
 143               break;
 144            }
 145            GetHitsData(1)->Fill(VZEROHit->Cell());
 146         }
 147 }
 148
 149
 150 //____________________________________________________________________________
 151
 152 void AliVZEROQADataMakerSim::MakeHits(TTree *hitTree)
 153 {
 154   //fills QA histos for Hits
 155  if (fHitsArray)
 156    fHitsArray->Clear() ;
 157   else
 158     fHitsArray = new TClonesArray("AliVZEROhit", 1000);
 159
 160   TBranch * branch = hitTree->GetBranch("VZERO") ;
 161   if ( ! branch ) {
 162     AliWarning("VZERO branch in Hit Tree not found") ;
 163   } else {
 164
 165    if (branch) {
 166       branch->SetAddress(&fHitsArray);
 167     }else{
 168       AliError("Branch VZERO hit not found");
 169       exit(111);
 170     }
 171     // Check id histograms already created for this Event Specie
 172     if ( ! GetHitsData(0) )
 173       InitHits() ;
 174
 175     Int_t ntracks    = (Int_t) hitTree->GetEntries();
 176
 177     if (ntracks<=0) return;
 178     // Start loop on tracks in the hits containers
 179     for (Int_t track=0; track<ntracks;track++) {
 180       branch->GetEntry(track);
 181       Int_t nhits = fHitsArray->GetEntriesFast();
 182       GetHitsData(0)->Fill(nhits) ;    // fills Hit multiplicity
 183       for (Int_t ihit=0;ihit<nhits;ihit++)
 184         {
 185           AliVZEROhit  * VZEROHit   = (AliVZEROhit*) fHitsArray->UncheckedAt(ihit);
 186           if (!VZEROHit) {
 187             AliError("The unchecked hit doesn't exist");
 188             break;
 189           }
 190           GetHitsData(1)->Fill(VZEROHit->Cell());
 191         }
 192     }
 193   }
 194 }
 195
 196
 197 //____________________________________________________________________________
 198 void AliVZEROQADataMakerSim::MakeDigits()
 199 {
 200   // makes data from Digits
 201
 202   GetDigitsData(0)->Fill(fDigitsArray->GetEntriesFast()) ;
 203     TIter next(fDigitsArray) ;
 204     AliVZEROdigit *VZERODigit ;
 205     while ( (VZERODigit = dynamic_cast<AliVZEROdigit *>(next())) ) {
 206          Int_t   PMNumber  = VZERODigit->PMNumber();
 207          GetDigitsData(PMNumber +1)->Fill( VZERODigit->Time()) ;    // in 100 of picoseconds
 208          GetDigitsData(PMNumber +1+64)->Fill( VZERODigit->ADC()) ;
 209     }
 210 }
 211
 212
 213 //____________________________________________________________________________
 214 void AliVZEROQADataMakerSim::MakeDigits(TTree *digitTree)
 215 {
 216     // makes data from Digit Tree
 217
 218   if (fDigitsArray)
 219     fDigitsArray->Clear() ;
 220   else
 221     fDigitsArray = new TClonesArray("AliVZEROdigit", 1000) ;
 222
 223     TBranch * branch = digitTree->GetBranch("VZERODigit") ;
 224     if ( ! branch ) {
 225          AliWarning("VZERO branch in Digit Tree not found") ;
 226     } else {
 227          branch->SetAddress(&fDigitsArray) ;
 228          branch->GetEntry(0) ;
 229          MakeDigits() ;
 230     }
 231 }
 232
 233
 234 //____________________________________________________________________________
 235 void AliVZEROQADataMakerSim::StartOfDetectorCycle()
 236 {
 237   //Detector specific actions at start of cycle
 238
 239 }