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(1,"");
  56 }
  57
  58 //____________________________________________________________________________
  59 AliMUONQADataMakerSim::AliMUONQADataMakerSim(const AliMUONQADataMakerSim& qadm) :
  60     AliQADataMakerSim(),
  61   fHitStore(0x0),
  62   fDigitStore(0x0)
  63 {
  64   /// Copy constructor
  65
  66   AliDebug(1,"");
  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(1,"");
  86
  87     this->~AliMUONQADataMakerSim();
  88     new(this) AliMUONQADataMakerSim(qadm);
  89     return *this;
  90 }
  91
  92 //__________________________________________________________________
  93 AliMUONQADataMakerSim::~AliMUONQADataMakerSim()
  94 {
  95   /// Destructor
  96
  97   AliDebug(1,"");
  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", 1400, 100., 1500.);
 111   Add2HitsList(h0, 0, !expert, image);
 112
 113   TH1F* h1 = new TH1F("hHitPtot", "P distribution in Hits ", 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",  1400, 100, 1500);
 126   Add2SDigitsList(h0, 0, !expert, image);
 127
 128   TH1F* h1 = new TH1F("hSDigitsCharge", "Charge distribution in SDigits", 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",  1400, 100, 1500);
 141   Add2DigitsList(h0, 0, !expert, image);
 142
 143   TH1I* h1 = new TH1I("hDigitsADC", "ADC distribution in Digits", 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   if (!fHitStore)
 153     fHitStore = AliMUONVHitStore::Create(*hitsTree);
 154   fHitStore->Connect(*hitsTree, false);
 155   hitsTree->GetEvent(0);
 156
 157   TIter next(fHitStore->CreateIterator());
 158
 159   AliMUONHit* hit = 0x0;
 160
 161   while ( ( hit = static_cast<AliMUONHit*>(next()) ) )
 162   {
 163     GetHitsData(0)->Fill(hit->DetElemId());
 164     GetHitsData(1)->Fill(hit->Momentum());
 165   }
 166
 167
 168 }
 169
 170 //__________________________________________________________________
 171 void AliMUONQADataMakerSim::MakeSDigits(TTree* sdigitsTree)
 172 {
 173   /// makes data from SDigits
 174   if (!fDigitStore)
 175     fDigitStore = AliMUONVDigitStore::Create(*sdigitsTree);
 176   fDigitStore->Connect(*sdigitsTree, false);
 177   sdigitsTree->GetEvent(0);
 178
 179   TIter next(fDigitStore->CreateIterator());
 180
 181   AliMUONVDigit* dig = 0x0;
 182
 183   while ( ( dig = static_cast<AliMUONVDigit*>(next()) ) )
 184   {
 185     GetSDigitsData(0)->Fill(dig->DetElemId());
 186     GetSDigitsData(1)->Fill(dig->Charge());
 187   }
 188 }
 189
 190 //__________________________________________________________________
 191 void AliMUONQADataMakerSim::MakeDigits(TTree* digitsTree)
 192 {
 193    /// makes data from Digits
 194   if (!fDigitStore)
 195     fDigitStore = AliMUONVDigitStore::Create(*digitsTree);
 196   fDigitStore->Connect(*digitsTree, false);
 197   digitsTree->GetEvent(0);
 198
 199   TIter next(fDigitStore->CreateIterator());
 200
 201   AliMUONVDigit* dig = 0x0;
 202
 203   while ( ( dig = static_cast<AliMUONVDigit*>(next()) ) )
 204   {
 205     GetDigitsData(0)->Fill(dig->DetElemId());
 206     GetDigitsData(1)->Fill(dig->ADC());
 207   }
 208 }
 209
 210 //____________________________________________________________________________
 211 void AliMUONQADataMakerSim::EndOfDetectorCycle(AliQAv1::TASKINDEX_t task, TObjArray** list)
 212 {
 213     ///Detector specific actions at end of cycle
 214     // do the QA checking
 215     AliQAChecker::Instance()->Run(AliQAv1::kMUON, task, list) ;
 216 }
 217
 218
 219 //____________________________________________________________________________
 220 void AliMUONQADataMakerSim::StartOfDetectorCycle()
 221 {
 222     /// Detector specific actions at start of cycle
 223
 224 }