MUON/AliMUONQADataMakerSim.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 // $Id$
  17
  18 #include "AliMUONQADataMakerSim.h"
  19 #include "AliMUONHit.h"
  20 #include "AliMUONDigit.h"
  21 #include "AliMUONVHitStore.h"
  22 #include "AliMUONVDigitStore.h"
  23
  24 // --- AliRoot header files ---
  25 #include "AliLog.h"
  26 #include "AliQAChecker.h"
  27
  28 // --- ROOT system ---
  29 #include <TClonesArray.h>
  30 #include <TFile.h>
  31 #include <TH1F.h>
  32 #include <TH1I.h>
  33 #include <TH2F.h>
  34 #include <TTree.h>
  35
  36 //-----------------------------------------------------------------------------
  37 /// \class AliMUONQADataMakerSim
  38 ///
  39 /// MUON base class for quality assurance data (histo) maker
  40 ///
  41 /// \author C. Finck
  42
  43 /// \cond CLASSIMP
  44 ClassImp(AliMUONQADataMakerSim)
  45 /// \endcond
  46
  47 //____________________________________________________________________________
  48 AliMUONQADataMakerSim::AliMUONQADataMakerSim() :
  49     AliQADataMakerSim(AliQAv1::GetDetName(AliQAv1::kMUON), "MUON Quality Assurance Data Maker"),
  50     fHitStore(0x0),
  51     fDigitStore(0x0)
  52 {
  53   /// Default constructor
  54
  55       AliDebug(AliQAv1::GetQADebugLevel(),"");
  56 }
  57
  58 //____________________________________________________________________________
  59 AliMUONQADataMakerSim::AliMUONQADataMakerSim(const AliMUONQADataMakerSim& qadm) :
  60     AliQADataMakerSim(),
  61   fHitStore(0x0),
  62   fDigitStore(0x0)
  63 {
  64   /// Copy constructor
  65
  66     AliDebug(AliQAv1::GetQADebugLevel(),"");
  67
  68     if ( qadm.fHitStore )
  69     {
  70       fHitStore = static_cast<AliMUONVHitStore*>(qadm.fHitStore->Clone());
  71     }
  72     if ( qadm.fDigitStore )
  73     {
  74       fDigitStore = static_cast<AliMUONVDigitStore*>(qadm.fDigitStore->Clone());
  75     }
  76     SetName((const char*)qadm.GetName()) ;
  77     SetTitle((const char*)qadm.GetTitle());
  78 }
  79
  80 //__________________________________________________________________
  81 AliMUONQADataMakerSim& AliMUONQADataMakerSim::operator = (const AliMUONQADataMakerSim& qadm )
  82 {
  83   /// Assignment operator
  84
  85   AliDebug(AliQAv1::GetQADebugLevel(),"");
  86
  87     this->~AliMUONQADataMakerSim();
  88     new(this) AliMUONQADataMakerSim(qadm);
  89     return *this;
  90 }
  91
  92 //__________________________________________________________________
  93 AliMUONQADataMakerSim::~AliMUONQADataMakerSim()
  94 {
  95   /// Destructor
  96
  97   AliDebug(AliQAv1::GetQADebugLevel(),"");
  98
  99   delete fHitStore;
 100   delete fDigitStore;
 101 }
 102
 103 //__________________________________________________________________
 104 void AliMUONQADataMakerSim::InitHits()
 105 {
 106   /// Initialized hit spectra
 107   const Bool_t expert   = kTRUE ;
 108   const Bool_t image    = kTRUE ;
 109
 110   TH1F* h0 = new TH1F("hHitDetElem", "DetElemId distribution in Hits;Detection element Id;Counts", 1400, 100., 1500.);
 111   Add2HitsList(h0, 0, !expert, image);
 112
 113   TH1F* h1 = new TH1F("hHitPtot", "P distribution in Hits;P [erg];Counts ", 300, 0., 300.);
 114   Add2HitsList(h1, 1, !expert, image);
 115   return;
 116 }
 117
 118 //__________________________________________________________________
 119 void AliMUONQADataMakerSim::InitSDigits()
 120 {
 121   /// Initialized SDigits spectra
 122   const Bool_t expert   = kTRUE ;
 123   const Bool_t image    = kTRUE ;
 124
 125   TH1I* h0 = new TH1I("hSDigitsDetElem", "Detection element distribution in SDigits;Detection element Id;Counts",  1400, 100, 1500);
 126   Add2SDigitsList(h0, 0, !expert, image);
 127
 128   TH1F* h1 = new TH1F("hSDigitsCharge", "Charge distribution in SDigits;Charge [??];Counts", 4096, 0, 4095);
 129   Add2SDigitsList(h1, 1, !expert, image);
 130
 131 }
 132
 133 //__________________________________________________________________
 134 void AliMUONQADataMakerSim::InitDigits()
 135 {
 136   /// Initialized Digits spectra
 137   const Bool_t expert   = kTRUE ;
 138   const Bool_t image    = kTRUE ;
 139
 140   TH1I* h0 = new TH1I("hDigitsDetElem", "Detection element distribution in Digits;Detection element Id;Counts",  1400, 100, 1500);
 141   Add2DigitsList(h0, 0, !expert, image);
 142
 143   TH1I* h1 = new TH1I("hDigitsADC", "ADC distribution in Digits;ACD value;Counts", 4096, 0, 4095);
 144   Add2DigitsList(h1, 1, !expert, image);
 145
 146 }
 147
 148 //__________________________________________________________________
 149 void AliMUONQADataMakerSim::MakeHits(TTree* hitsTree)
 150 {
 151   /// makes data from Hits
 152
 153   if (!fHitStore)
 154     fHitStore = AliMUONVHitStore::Create(*hitsTree);
 155   fHitStore->Connect(*hitsTree, false);
 156   hitsTree->GetEvent(0);
 157
 158   TIter next(fHitStore->CreateIterator());
 159
 160   AliMUONHit* hit = 0x0;
 161
 162   while ( ( hit = static_cast<AliMUONHit*>(next()) ) )
 163   {
 164     GetHitsData(0)->Fill(hit->DetElemId());
 165     GetHitsData(1)->Fill(hit->Momentum());
 166   }
 167
 168
 169 }
 170
 171 //__________________________________________________________________
 172 void AliMUONQADataMakerSim::MakeSDigits(TTree* sdigitsTree)
 173 {
 174   /// makes data from SDigits
 175
 176   if (!fDigitStore)
 177     fDigitStore = AliMUONVDigitStore::Create(*sdigitsTree);
 178   fDigitStore->Connect(*sdigitsTree, false);
 179   sdigitsTree->GetEvent(0);
 180
 181   TIter next(fDigitStore->CreateIterator());
 182
 183   AliMUONVDigit* dig = 0x0;
 184
 185   while ( ( dig = static_cast<AliMUONVDigit*>(next()) ) )
 186   {
 187     GetSDigitsData(0)->Fill(dig->DetElemId());
 188     GetSDigitsData(1)->Fill(dig->Charge());
 189   }
 190 }
 191
 192 //__________________________________________________________________
 193 void AliMUONQADataMakerSim::MakeDigits(TTree* digitsTree)
 194 {
 195    /// makes data from Digits
 196
 197   if (!fDigitStore)
 198     fDigitStore = AliMUONVDigitStore::Create(*digitsTree);
 199   fDigitStore->Connect(*digitsTree, false);
 200   digitsTree->GetEvent(0);
 201
 202   TIter next(fDigitStore->CreateIterator());
 203
 204   AliMUONVDigit* dig = 0x0;
 205
 206   while ( ( dig = static_cast<AliMUONVDigit*>(next()) ) )
 207   {
 208     GetDigitsData(0)->Fill(dig->DetElemId());
 209     GetDigitsData(1)->Fill(dig->ADC());
 210   }
 211 }
 212
 213 //____________________________________________________________________________
 214 void AliMUONQADataMakerSim::EndOfDetectorCycle(AliQAv1::TASKINDEX_t task, TObjArray** list)
 215 {
 216     ///Detector specific actions at end of cycle
 217     // do the QA checking
 218     AliQAChecker::Instance()->Run(AliQAv1::kMUON, task, list) ;
 219 }
 220
 221
 222 //____________________________________________________________________________
 223 void AliMUONQADataMakerSim::StartOfDetectorCycle()
 224 {
 225     /// Detector specific actions at start of cycle
 226
 227 }