In This Package:

CalibHistoCheck.cc

Go to the documentation of this file.

#include "CalibHistoCheck.h"

#include "StatisticsSvc/IStatisticsSvc.h"

#include "Event/CalibReadoutHeader.h"
#include "Event/CalibReadout.h"
#include "Event/CalibReadoutPmtCrate.h"
#include "Event/CalibReadoutPmtChannel.h"

#include "TH1F.h"
#include "TH2D.h"

#include <vector>

using namespace std;
using namespace DayaBay;

CalibHistoCheck::CalibHistoCheck(const std::string& name,
                                 ISvcLocator* pSvcLocator) : 
  GaudiAlgorithm(name,pSvcLocator)
{
}

StatusCode CalibHistoCheck::initialize()
{
  if ( service("StatisticsSvc", p_statsSvc).isFailure() ) {
    error()<<" No StatisticsSvc available ! "<<endreq;
    return StatusCode::FAILURE;
  } 
  
  TH1F*  h1 = new TH1F("Tdc","Tdc", 500, -500., 0. );
  h1->GetXaxis()->SetTitle("ns");
  if ( p_statsSvc->put("/File1/Calib/CalibPmt/Tdc", h1).isFailure() ) {
    error()<<"Could't register Tdc "<<endreq; 
  }
  
  TH1F*  h2 = new TH1F("Adc","Adc", 400, 0., 40. );
  h2->GetXaxis()->SetTitle("p.e");
  if ( p_statsSvc->put("/File1/Calib/CalibPmt/Adc", h2).isFailure() ) {
    error()<<"Could't register Adc "<<endreq; 
  }

  info() << "Finished booking Histograms" << endreq;
  return StatusCode::SUCCESS;
}

StatusCode CalibHistoCheck::execute()
{
  
  TH1* tdcH = p_statsSvc->getTH1F("/File1/Calib/CalibPmt/Tdc");
  TH1* adcH = p_statsSvc->getTH1F("/File1/Calib/CalibPmt/Adc");
  if (tdcH && adcH) {
    
    const  DayaBay::CalibReadoutHeader* calibHeader 
      = get<DayaBay::CalibReadoutHeader>(
                                 DayaBay::CalibReadoutHeaderLocation::Default);
    
    if(!calibHeader) {
      error()<< "Could not find calibHeader" << endreq;
      return StatusCode::FAILURE; 
    }else {
      if(!calibHeader->calibReadout()) {
        info() << "No CalibReadout in calibHeader. No tdc/adc will be filled." << endreq;
        return StatusCode::SUCCESS;
      }
      
      const DayaBay::CalibReadoutPmtCrate* calibCrate 
        = dynamic_cast<const DayaBay::CalibReadoutPmtCrate*>(
                                                  calibHeader->calibReadout());
      const DayaBay::CalibReadoutPmtCrate::PmtChannelReadouts calib_chVec 
        = calibCrate->channelReadout();

      DayaBay::CalibReadoutPmtCrate::PmtChannelReadouts::const_iterator cpcrIter
        = calib_chVec.begin();

      for(; cpcrIter!=calib_chVec.end(); ++cpcrIter) {
        
        // Get charge on this PMT
        double charge = cpcrIter->maxCharge();
        adcH->Fill(charge);

        // Get time on this PMT
        vector<double> time = cpcrIter->time();
        for(size_t len=0; len<time.size(); len++) tdcH->Fill(time[len]);
      }
        
    }
  } else {
    error() <<" Could not retrieve Histograms "<<endreq;
    return StatusCode::FAILURE;
  } 
  
  return StatusCode::SUCCESS;
}

StatusCode CalibHistoCheck::finalize()
{
  debug() << "Finalizing..." << endreq;
  
  return StatusCode::SUCCESS;
}

---