In This Package:

ReconAlg.cc

Go to the documentation of this file.

#include "ReconAlg.h"
#include "ReconAlg/IReconTool.h"

#include "Conventions/Reconstruction.h"
#include "Event/CalibReadoutHeader.h"
#include "Event/CalibReadout.h"

#include "Event/RecTrigger.h"

#include <vector>

ReconAlg::ReconAlg(const std::string& name, ISvcLocator* pSvcLocator)
  : DybAlgorithm<DayaBay::RecHeader>(name,pSvcLocator)
{
  declareProperty("calibLocation", 
                  m_calibLocation=DayaBay::CalibReadoutHeaderLocation::Default,
                  "CalibReadoutHeader location in the TES.");

  declareProperty("ReconTools", m_reconToolNames,
                  "List of Reconstruction tools");

}

ReconAlg::~ReconAlg()
{
}

StatusCode ReconAlg::initialize()
{
  // Setup reconstruction styles and get tools
  for(unsigned int toolIdx = 0; toolIdx < m_reconToolNames.size(); toolIdx++){
    std::string toolName = m_reconToolNames[toolIdx];
    IReconTool* currentTool = 0;
    try {
      currentTool = tool<IReconTool>(toolName);
    }
    catch(const GaudiException& exg) {
      fatal() << "Failed to get Recon tool: \"" << toolName << "\"" << endreq;
      return StatusCode::FAILURE;
    }
    m_reconTools.push_back(currentTool);
    info() << "Added recon tool " << toolName << " to path " 
           << this->Location() << endreq;
  }

  return StatusCode::SUCCESS;
}

StatusCode ReconAlg::execute()
{
  StatusCode sc;

  // Create output rec header
  DayaBay::RecHeader* recHeader = MakeHeaderObject();

  // Get CalibHeader from TES
  DayaBay::CalibReadoutHeader* calibHeader = 
    getTES<DayaBay::CalibReadoutHeader>(m_calibLocation);
  if(!calibHeader) {
    error() << "Get no CalibHeader from TES!!" << endreq;
    return StatusCode::FAILURE;
  }

  Context context = calibHeader->context();
  recHeader->setContext(context);
  recHeader->setEarliest(calibHeader->earliest());
  recHeader->setLatest(calibHeader->latest());

  DayaBay::RecTrigger& recTrigger = recHeader->recTrigger();

  // Get Calibrated Readout
  const DayaBay::CalibReadout* calibReadout = calibHeader->calibReadout();
  if(!calibReadout) {
    info() << "No CalibReadout retrieved this cycle." << endreq;
    recTrigger.setEnergyStatus(ReconStatus::kUnknown);
    recTrigger.setPositionStatus(ReconStatus::kUnknown);
    recTrigger.setDirectionStatus(ReconStatus::kUnknown);
    return StatusCode::SUCCESS;
  }
  
  // Get detector / trigger information 
  DayaBay::Detector det = calibReadout->detector();
  unsigned int triggerNumber = calibReadout->triggerNumber();
  const TimeStamp triggerTime = calibReadout->triggerTime();
  const DayaBay::Trigger::TriggerType_t triggerType 
    = calibReadout->triggerType();
  
  recTrigger.setDetector(det);
  recTrigger.setTriggerNumber(triggerNumber);
  recTrigger.setTriggerTime(triggerTime);
  recTrigger.setTriggerType(triggerType);
  
  // Run the recon tools for this style
  for(unsigned int toolIdx = 0; toolIdx < m_reconTools.size(); toolIdx++){
    sc = m_reconTools[toolIdx]->reconstruct(*calibReadout, recTrigger);
    if(sc != StatusCode::SUCCESS) return sc;
  }

  return StatusCode::SUCCESS;
}

StatusCode ReconAlg::finalize()
{
  return this->GaudiAlgorithm::finalize();
}

---