FMD/AliFMDQADataMakerSim.cxx

   1 /**************************************************************************
   2  * Copyright(c) 2004, 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 // --- ROOT system ---
  16 #include <iostream>
  17 #include <TClonesArray.h>
  18 #include <TFile.h>
  19 #include <TH1F.h>
  20 #include <TH1I.h>
  21
  22 // --- AliRoot header files ---
  23 #include "AliESDEvent.h"
  24 #include "AliLog.h"
  25 #include "AliFMDQADataMakerSim.h"
  26 #include "AliFMDDigit.h"
  27 #include "AliFMDHit.h"
  28 #include "AliQAChecker.h"
  29 #include "AliFMDParameters.h"
  30
  31 //_____________________________________________________________________
  32 // This is the class that collects the QA data for the FMD during simulation.
  33 // The following data types are picked up:
  34 // - hits
  35 // - digits
  36 // The following data types are not supported (yet):
  37 // - raws
  38 // - sdigits
  39 // Author : Hans Hjersing Dalsgaard, Niels Bohr Institute, hans.dalsgaard@cern.ch
  40 //_____________________________________________________________________
  41
  42 ClassImp(AliFMDQADataMakerSim)
  43 #if 0
  44 ; // This line is for Emacs - do not delete!
  45 #endif
  46 //_____________________________________________________________________
  47 AliFMDQADataMakerSim::AliFMDQADataMakerSim()
  48   :  AliQADataMakerSim(AliQA::GetDetName(AliQA::kFMD),
  49                        "FMD Quality Assurance Data Maker"),
  50      fDigitsArray(0),
  51      fHitsArray(0)
  52 {
  53   // ctor
  54   fDigitsArray = new TClonesArray("AliFMDDigit", 1000) ;
  55   fHitsArray   = new TClonesArray("AliFMDHit", 10);
  56 }
  57
  58 //_____________________________________________________________________
  59 AliFMDQADataMakerSim::AliFMDQADataMakerSim(const AliFMDQADataMakerSim& qadm)
  60   : AliQADataMakerSim(),
  61     fDigitsArray(qadm.fDigitsArray),
  62     fHitsArray(qadm.fHitsArray)
  63 {
  64   //copy ctor
  65
  66   // Parameters:
  67   //    qadm    Object to copy from
  68
  69 }
  70 //_____________________________________________________________________
  71 AliFMDQADataMakerSim& AliFMDQADataMakerSim::operator = (const AliFMDQADataMakerSim& qadm )
  72 {
  73   fDigitsArray = qadm.fDigitsArray;
  74   fHitsArray = qadm.fHitsArray;
  75
  76   return *this;
  77 }
  78 //_____________________________________________________________________
  79 AliFMDQADataMakerSim::~AliFMDQADataMakerSim()
  80 {
  81   delete fDigitsArray;
  82   delete fHitsArray;
  83 }
  84
  85 //_____________________________________________________________________
  86 void AliFMDQADataMakerSim::EndOfDetectorCycle(AliQA::TASKINDEX_t task,
  87                                               TObjArray * list)
  88 {
  89   //Detector specific actions at end of cycle
  90   // do the QA checking
  91   AliLog::Message(5,"FMD: end of detector cycle",
  92                   "AliFMDQADataMakerSim","AliFMDQADataMakerSim",
  93                   "AliFMDQADataMakerSim::EndOfDetectorCycle",
  94                   "AliFMDQADataMakerSim.cxx",83);
  95   AliQAChecker::Instance()->Run(AliQA::kFMD, task, list) ;
  96
  97 }
  98
  99 //____________________________________________________________________
 100 void AliFMDQADataMakerSim::InitHits()
 101 {
 102   // create Digits histograms in Digits subdir
 103   TH1F* hEnergyOfHits = new TH1F("hEnergyOfHits","Energy distribution",100,0,3);
 104   hEnergyOfHits->SetXTitle("Edep");
 105   hEnergyOfHits->SetYTitle("Counts");
 106   Add2HitsList(hEnergyOfHits, 0);
 107 }
 108
 109 //_____________________________________________________________________
 110 void AliFMDQADataMakerSim::InitDigits()
 111 {
 112   // create Digits histograms in Digits subdir
 113   TH1I* hADCCounts = new TH1I("hADCCounts","Dist of ADC counts",1024,0,1024);
 114   hADCCounts->SetXTitle("ADC counts");
 115   Add2DigitsList(hADCCounts, 0);
 116 }
 117
 118 //_____________________________________________________________________
 119 void AliFMDQADataMakerSim::MakeHits(TClonesArray * hits)
 120 {
 121   TIter next(hits);
 122   AliFMDHit * hit;
 123   while ((hit = static_cast<AliFMDHit *>(next())))
 124     GetHitsData(0)->Fill(hit->Edep()/hit->Length()*0.032);
 125 }
 126
 127 //_____________________________________________________________________
 128 void AliFMDQADataMakerSim::MakeHits(TTree * hitTree)
 129 {
 130   // make QA data from Hit Tree
 131
 132   TBranch * branch = hitTree->GetBranch("FMD") ;
 133   if (!branch) {
 134     AliWarning("FMD branch in Hit Tree not found") ;
 135     return;
 136   }
 137   fHitsArray->Clear();
 138   branch->SetAddress(&fHitsArray) ;
 139
 140   for (Int_t ientry = 0 ; ientry < branch->GetEntries() ; ientry++) {
 141     branch->GetEntry(ientry);
 142     MakeHits(fHitsArray);   //tmp);
 143   }
 144
 145 }
 146
 147 //_____________________________________________________________________
 148 void AliFMDQADataMakerSim::MakeDigits(TClonesArray * digits)
 149 {
 150   // makes data from Digits
 151   if(!digits) return;
 152
 153   for(Int_t i = 0 ; i < fDigitsArray->GetEntriesFast() ; i++) {
 154     //Raw ADC counts
 155     AliFMDDigit* digit = static_cast<AliFMDDigit*>(digits->At(i));
 156     GetDigitsData(0)->Fill(digit->Counts());
 157   }
 158 }
 159
 160 //_____________________________________________________________________
 161 void AliFMDQADataMakerSim::MakeDigits(TTree * digitTree)
 162 {
 163
 164   fDigitsArray->Clear();
 165   TBranch * branch = digitTree->GetBranch("FMD") ;
 166   if (!branch)    {
 167       AliWarning("FMD branch in Digit Tree not found") ;
 168       return;
 169   }
 170   branch->SetAddress(&fDigitsArray) ;
 171   branch->GetEntry(0) ;
 172   MakeDigits(fDigitsArray) ;
 173 }
 174
 175 //_____________________________________________________________________
 176 void AliFMDQADataMakerSim::StartOfDetectorCycle()
 177 {
 178
 179 }
 180 //_____________________________________________________________________
 181 //
 182 // EOF
 183 //