In This Package:

CalibAlg.cc

Go to the documentation of this file.

#include "CalibAlg.h"
#include "CalibAlg/ICalibTool.h"

#include "Event/CalibReadoutPmtCrate.h"
#include "Event/CalibReadoutRpcCrate.h"

#include <vector>

CalibAlg::CalibAlg(const std::string& name, ISvcLocator* pSvcLocator)
  : DybAlgorithm<DayaBay::CalibReadoutHeader>(name,pSvcLocator)
  , m_calibRpcTool(0)
  , m_calibPmtTool(0)
{
  declareProperty("ReadoutLocation", 
                  m_roLocation=DayaBay::ReadoutHeaderLocation::Default,
                  "ReadoutHeader location in the TES.");

  declareProperty("CalibReadoutLocation", 
                 m_calibRoLocation=DayaBay::CalibReadoutHeaderLocation::Default,
                 "CalibHeader location in the TES.");

  declareProperty("RpcCalibTool", m_calibRpcToolName="SimpRpcCalibTool",
                  "Tool to calibrate RPC readouts");

  declareProperty("PmtCalibTool", m_calibPmtToolName="SimpPmtCalibTool",
                  "Tool to calibrate PMT readouts");
}

CalibAlg::~CalibAlg()
{
}

StatusCode CalibAlg::initialize()
{
  // Get RPC calibration Tool
  try {
    m_calibRpcTool = tool<ICalibTool>(m_calibRpcToolName);
  }
  catch(const GaudiException& exg) {
    fatal() << "Failed to get Calib Tool: \"" << m_calibRpcToolName << "\"" << endreq;
    return StatusCode::FAILURE;
  }
  info () << "Added calib tool " << m_calibRpcToolName << endreq;

  try {
    m_calibPmtTool = tool<ICalibTool>(m_calibPmtToolName);
  }
  catch(const GaudiException& exg) {
    fatal() << "Failed to get Calib Tool: \"" << m_calibPmtToolName << "\"" << endreq;
    return StatusCode::FAILURE;
  }
  info () << "Added calib tool " << m_calibPmtToolName << endreq;

  return StatusCode::SUCCESS;
}

StatusCode CalibAlg::execute()
{
  StatusCode sc;

  // Create output calib header
  DayaBay::CalibReadoutHeader* calibHeader = MakeHeaderObject();
  
  // Get ReadoutHeader
  DayaBay::ReadoutHeader* roHeader 
    = getTES<DayaBay::ReadoutHeader>(m_roLocation);
  if(!roHeader) { 
    error() << "Get no ReadoutHeader from TES!!" << endreq;
    return StatusCode::FAILURE;
  }

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

  // Get Readout  
  const DayaBay::DaqCrate* readout = roHeader->daqCrate();
  if(!readout){
    info() << "No Readouts retrieved this cycle." << endreq;
    calibHeader->setCalibReadout(0);
    return StatusCode::SUCCESS;
  }

  const DayaBay::Detector& detector = readout->detector();
  unsigned int eventNumber = readout->eventNumber();
  const TimeStamp& triggerTime = readout->triggerTime();
  const DayaBay::Trigger::TriggerType_t& triggerType = readout->triggerType();

  debug() << "detector " << detector << endreq;
  debug() << "eventNumber " << eventNumber << endreq;
  debug() << "triggerTime " << triggerTime << endreq;
  debug() << "triggerType " << triggerType << endreq;

  if(detector.isAD() || detector.isWaterShield()){
    // Pmt based detectors calibration
    debug() << "create CalibReadoutPmtCrate " << endreq;
    DayaBay::CalibReadoutPmtCrate* calibReadout = 
      new DayaBay::CalibReadoutPmtCrate(detector, eventNumber, 
                                        triggerTime, triggerType);
    calibHeader->setCalibReadout(calibReadout);
    calibReadout->setHeader(calibHeader);
    debug() << "start calibration readout " << endreq;
    sc = m_calibPmtTool->calibration(readout, calibReadout);
    debug() << "finish calibration readout " << endreq;

  }else if(detector.detectorId() == DetectorId::kRPC){
    // RPC Based detectors Calibration
    DayaBay::CalibReadoutRpcCrate* calibReadout = 
      new DayaBay::CalibReadoutRpcCrate(detector,eventNumber, 
                                        triggerTime, triggerType);
    calibHeader->setCalibReadout(calibReadout);
    calibReadout->setHeader(calibHeader);
    sc = m_calibRpcTool->calibration(readout, calibReadout);

  }else {
    error() << "Unknown detector " << detector << endreq;
    return StatusCode::FAILURE;
  }

  return StatusCode::SUCCESS;
}

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

---