ReconProcessorAlg.cc

Go to the documentation of this file.

#include "ReconProcessorAlg.h"

#include "ProcessTools/IReconProcessor.h"
#include "Event/RecHeader.h"

ReconProcessorAlg::ReconProcessorAlg(const std::string& name, 
                                         ISvcLocator* pSvcLocator)
  : GaudiAlgorithm(name,pSvcLocator),
    m_tools()
{
  declareProperty("Location", m_location="/Event/Rec/AdSimple",
                  "Location of reconstructed data");
  declareProperty("Tools", m_toolNames,
                  "List of tool names to apply to the recon data");
}

ReconProcessorAlg::~ReconProcessorAlg()
{
}

StatusCode ReconProcessorAlg::initialize()
{
  // Initialize the tools you want to use to process the readout
  std::vector<std::string>::iterator toolNameIter, 
    toolNameEnd=m_toolNames.end();
  for(toolNameIter=m_toolNames.begin(); toolNameIter!=toolNameEnd; 
      toolNameIter++){
    IReconProcessor* currentTool = 0;
    try {
      currentTool = tool<IReconProcessor>(*toolNameIter);
    }
    catch(const GaudiException& exg) {
      fatal() << "Failed to get tool: \"" << *toolNameIter << "\"" << endreq;
      return StatusCode::FAILURE;
    }
    m_tools.push_back(currentTool);
    info() << "Added tool " << *toolNameIter << endreq;
  }
  return StatusCode::SUCCESS;
}

StatusCode ReconProcessorAlg::execute()
{
  // Add the current event into histograms
  DayaBay::RecHeader* reconHeader = 
    get<DayaBay::RecHeader>(m_location);
  if(!reconHeader){
    error() << "Failed to get reconstructed data header." << endreq;
    return StatusCode::FAILURE;
  }

  // Process reconstructed data using the list of tools 
  StatusCode sc;
  std::vector<IReconProcessor*>::iterator toolIter, toolEnd=m_tools.end();
  for(toolIter=m_tools.begin(); toolIter!=toolEnd; toolIter++){
    IReconProcessor* currentTool = *toolIter;
    sc = currentTool->process( reconHeader );
    if( !sc.isSuccess() ) return sc;
  }

  return StatusCode::SUCCESS;
}

StatusCode ReconProcessorAlg::finalize()
{
  // Clean up tools
  std::vector<IReconProcessor*>::iterator toolIter, toolEnd=m_tools.end();
  for(toolIter=m_tools.begin(); toolIter!=toolEnd; toolIter++){
    IReconProcessor* currentTool = *toolIter;
    if(currentTool) currentTool->release();
  }
  return StatusCode::SUCCESS;
}