ReadTBLArCalibDigits.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2021 CERN for the benefit of the ATLAS collaboration
*/
 
#include "TBCnv/ReadTBLArCalibDigits.h"
#include "LArIdentifier/LArOnlineID.h"
#include <vector>
 
ReadTBLArCalibDigits::ReadTBLArCalibDigits(const std::string& name, ISvcLocator* pSvcLocator)
  : AthAlgorithm(name, pSvcLocator),
    m_count(0),
    m_onlineHelper(0)
{
  declareProperty("ContainerKey",m_containerKey="");
  declareProperty("DumpFile",m_dumpFile="");
}
 
ReadTBLArCalibDigits::~ReadTBLArCalibDigits() 
{
}
 
StatusCode ReadTBLArCalibDigits::initialize()
{
  MsgStream log(msgSvc(), name());
  log << MSG::INFO << "Initialize" << endmsg;
  
  StatusCode sc = detStore()->retrieve(m_onlineHelper, "LArOnlineID");
  if (sc.isFailure()) {
    log << MSG::ERROR << "Could not get LArOnlineID helper" << endmsg;
    return StatusCode::FAILURE;
  }
  
  if (m_dumpFile.size()>0)
    m_outfile.open(m_dumpFile.c_str(),std::ios::out);
 
  m_count=0;
  
  log << MSG::INFO << "======== ReadTBLArCalibDigits initialize successfully ========" << endmsg;
  return StatusCode::SUCCESS;
}
 
 
StatusCode ReadTBLArCalibDigits::execute(const EventContext& /*ctx*/) {
  MsgStream log(msgSvc(), name());
  
  m_count++;
 
  int cellCounter = 0;
  
  StatusCode sc;
  
  log << MSG::DEBUG << "======== executing event "<< m_count << " ========" << endmsg;
  log << MSG::DEBUG << "Retrieving TBLArCalibDigitContainer. Key= " << m_containerKey << endmsg; 
  TBLArCalibDigitContainer* larCalibDigitCont;
  sc = evtStore()->retrieve(larCalibDigitCont ,m_containerKey);
 
  if (sc.isFailure()) {
    log << MSG::FATAL << " Cannot read TBLArCalibDigitContainer from StoreGate! key=" << m_containerKey << endmsg;
    return StatusCode::FAILURE;
  }
  
  log << MSG::VERBOSE << "Now loop over container " << endmsg;
 
  for (const LArCalibDigit* digit : *larCalibDigitCont) {
 
    log << MSG::VERBOSE << " Get hardware ID " << endmsg;
    
    HWIdentifier chid=digit->hardwareID();
    
    log << MSG::VERBOSE << "  chid =  " << chid << endmsg;
 
    const std::vector<short>& vSamples=digit->samples();    
    int nSamples=vSamples.size();
    
    log << MSG::VERBOSE << "  nSamples =  " << nSamples << endmsg;
 
    if (m_outfile.is_open()) {
      m_outfile << "Evt="<< m_count << " ";
      m_outfile << "FebId= 0x" << std::hex << m_onlineHelper->feb_Id(chid).get_compact()
                << std::dec << " Ch= " << m_onlineHelper->channel(chid) << " ";
      for(int i=0; i<nSamples; i++) {
        m_outfile << " " << vSamples[i];
      }
      m_outfile << " G=" << digit->gain() << std::endl;
    }
 
    cellCounter++;
 
    log << MSG::VERBOSE << "  cellCounter =  " << cellCounter << endmsg;
    
  }
  
  log << MSG::DEBUG  << "Event " << m_count << " contains " << cellCounter << " channels" << endmsg;
  
  return StatusCode::SUCCESS;
  
}
 
StatusCode ReadTBLArCalibDigits::finalize()
{
  MsgStream log(msgSvc(), name());
  if (m_outfile.is_open()) 
    m_outfile.close();
  log << MSG::INFO << "finalize ReadTBLArCalibDigits" << endmsg;
  return StatusCode::SUCCESS;
}