ACORDE/AliACORDEQADataMakerSim.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
  18 //---
  19 //  Produces the data needed to calculate the quality assurance.
  20 //  All data must be mergeable objects.
  21
  22 //  Authors:
  23 //
  24 //  Luciano Diaz Gonzalez <luciano.diaz@nucleares.unam.mx> (ICN-UNAM)
  25 //  Mario Rodriguez Cahuantzi <mrodrigu@mail.cern.ch> (FCFM-BUAP)
  26 //  Arturo Fernandez Tellez <afernan@mail.cern.ch (FCFM-BUAP)
  27 //
  28 //  Created: June 13th 2008
  29 //---
  30
  31 // --- ROOT system ---
  32 #include <TClonesArray.h>
  33 #include <TFile.h>
  34 #include <TH1F.h>
  35 #include <TH2F.h>
  36 #include <TDirectory.h>
  37 // --- Standard library ---
  38
  39 // --- AliRoot header files ---
  40 #include "AliESDEvent.h"
  41 #include "AliLog.h"
  42 #include "AliACORDEdigit.h"
  43 #include "AliACORDEhit.h"
  44 #include "AliACORDEQADataMakerSim.h"
  45 #include "AliQAChecker.h"
  46 #include "AliACORDERawReader.h"
  47 ClassImp(AliACORDEQADataMakerSim)
  48
  49 //____________________________________________________________________________
  50 AliACORDEQADataMakerSim::AliACORDEQADataMakerSim():AliQADataMakerSim(AliQA::GetDetName(AliQA::kACORDE), "ACORDE Quality Assurance Data Maker")
  51 {
  52 }
  53 //____________________________________________________________________________
  54 AliACORDEQADataMakerSim::AliACORDEQADataMakerSim(const AliACORDEQADataMakerSim& qadm) :
  55   AliQADataMakerSim()
  56 {
  57   SetName((const char*)qadm.GetName()) ;
  58   SetTitle((const char*)qadm.GetTitle());
  59 }
  60 //__________________________________________________________________
  61 AliACORDEQADataMakerSim& AliACORDEQADataMakerSim::operator = (const AliACORDEQADataMakerSim& qadm )
  62 {
  63   // Equal operator.
  64   this->~AliACORDEQADataMakerSim();
  65   new(this) AliACORDEQADataMakerSim(qadm);
  66   return *this;
  67 }
  68 //____________________________________________________________________________
  69 void AliACORDEQADataMakerSim::EndOfDetectorCycle(AliQA::TASKINDEX_t task, TObjArray * list)
  70 {
  71   //Detector specific actions at end of cycle
  72   // do the QA checking
  73    printf("ACORDE---->Detector specific actions at END of cycle\n................\n");
  74
  75   AliQAChecker::Instance()->Run(AliQA::kACORDE, task, list) ;
  76 }
  77 //____________________________________________________________________________
  78 void AliACORDEQADataMakerSim::StartOfDetectorCycle()
  79 {
  80   //Detector specific actions at start of cycle
  81   printf("ACORDE---->Detector specific actions at START of cycle\n................\n");
  82 }
  83 //____________________________________________________________________________
  84 void AliACORDEQADataMakerSim::InitHits()
  85 {
  86   // create Hits histograms in Hits subdir
  87
  88         TH1F *   fHitsACORDE;
  89         fHitsACORDE = new TH1F("hACORDEBitPattern","Distribution of fired modules",60,0,60);
  90         Add2HitsList(fHitsACORDE,0);
  91 }
  92 //____________________________________________________________________________
  93 void AliACORDEQADataMakerSim::InitDigits()
  94 {
  95   // create Digits histograms in Digits subdir
  96
  97    TH1F *    fhDigitsModule;
  98    TString   modulename;
  99    modulename = "hDigitsModule";
 100    fhDigitsModule = new TH1F(modulename.Data(),"hDigitsModuleSingle",60,0,60);
 101    Add2DigitsList( fhDigitsModule,0);
 102
 103 }
 104 //____________________________________________________________________________
 105
 106 void AliACORDEQADataMakerSim::MakeHits(TTree *hitTree)
 107 {
 108   // Here we fill the QA histos for Hits declared above
 109
 110         TClonesArray * hits = new TClonesArray("AliACORDEhit",1000);
 111         TBranch * branch = hitTree->GetBranch("ACORDE");
 112         if (!branch)
 113         {
 114                 AliWarning("ACORDE branch in Hit Tree not found");
 115         }else
 116         {
 117                 if (branch)
 118                 {
 119                         branch->SetAddress(&hits);
 120                 }else
 121                 {
 122                         AliError("Branch ACORDE hit not found");
 123                         exit(111);
 124                 }
 125                 Int_t ntracks = (Int_t)hitTree->GetEntries();
 126                 if (ntracks<=0) return;
 127                 for(Int_t track=0;track<ntracks;track++)
 128                 {
 129                         branch->GetEntry(track);
 130                         Int_t nhits = hits->GetEntriesFast();
 131                         for(Int_t ihit=1;ihit<=nhits;ihit++)
 132                         {
 133                                 AliACORDEhit *AcoHit = (AliACORDEhit*) hits->UncheckedAt(ihit);
 134                                 if(!AcoHit)
 135                                 {
 136                                         AliError("The unchecked hit doesn't exist");
 137                                         break;
 138                                 }
 139                                 GetHitsData(0)->Fill(AcoHit->GetModule()-1);
 140                         }
 141                 }
 142         }
 143
 144 }
 145 //____________________________________________________________________________
 146 void AliACORDEQADataMakerSim::MakeDigits( TTree *digitsTree)
 147 {
 148   //fills QA histos for Digits
 149
 150
 151         TClonesArray * digits = new TClonesArray("AliACORDEdigit",1000);
 152         TBranch * branch = digitsTree->GetBranch("ACORDEdigit");
 153         if (!branch)
 154         {
 155                 AliWarning("ACORDE branch in Digits Tree not found");
 156         }else
 157         {
 158                 if (branch)
 159                 {
 160                         branch->SetAddress(&digits);
 161                 }else
 162                 {
 163                         AliError("Branch ACORDE digit not found");
 164                         exit(111);
 165                 }
 166                 Int_t ntracks = (Int_t)digitsTree->GetEntries();
 167                 if (ntracks<=0) return;
 168                 printf("Entries in DigitsTree:%d\n",ntracks);
 169                 for(Int_t track=0;track<ntracks;track++)
 170                 {
 171                         branch->GetEntry(track);
 172                         Int_t ndigits = digits->GetEntriesFast();
 173                         for(Int_t idigit=0;idigit<ndigits;idigit++)
 174                         {
 175                                 AliACORDEdigit *AcoDigit = (AliACORDEdigit*) digits->UncheckedAt(idigit);
 176                                 if (!AcoDigit)
 177                                 {
 178                                         AliError("The unchecked digit doesn't exist");
 179                                         break;
 180                                 }
 181                                 GetDigitsData(0)->Fill(AcoDigit->GetModule()-1);
 182                         }
 183                 }
 184
 185
 186         }
 187
 188 }