ROsSequencerAlg.cc

Go to the documentation of this file.

#include "ROsSequencerAlg.h"

#include "Event/SimReadoutHeader.h"
#include "Event/ReadoutPmtCrate.h"
#include "Event/ReadoutRpcCrate.h"

ROsSequencerAlg::ROsSequencerAlg(const std::string& name, 
                                 ISvcLocator* pSvcLocator)
  : DybAlgorithm<DayaBay::ReadoutHeader>(name,pSvcLocator)
{
    declareProperty("SimReadoutLocation",m_simRoLocation=DayaBay::SimReadoutHeaderLocation::Default,
                    "Location in the TES where the input SimReadoutHeader is to be found.");
    declareProperty("ReadoutLocation",m_roLocation=DayaBay::ReadoutHeaderLocation::Default,
                    "Location in the TES where the output ReadoutHeader is to be put.");

}

ROsSequencerAlg::~ROsSequencerAlg()
{
}

StatusCode ROsSequencerAlg::initialize()
{
  return this->GaudiAlgorithm::initialize();
}

StatusCode ROsSequencerAlg::execute()
{
    // Look for pre-existing header object or make new one
    DayaBay::ReadoutHeader* roHeader = MakeHeaderObject();
    
    DayaBay::SimReadoutHeader* simReadoutHeader 
      = getTES<DayaBay::SimReadoutHeader>(m_simRoLocation);

    TimeStamp earliest = simReadoutHeader->earliest();
    TimeStamp latest = simReadoutHeader->latest();
    
    roHeader->setEarliest(earliest);
    roHeader->setLatest(latest);

    Context context = simReadoutHeader->context();
    
    
    
    const DayaBay::SimReadoutHeader::SimReadoutContainer& simReadouts
      = simReadoutHeader->readouts();
    DayaBay::SimReadoutHeader::SimReadoutContainer::const_iterator simReadoutIter 
      = simReadouts.begin();

    if(simReadouts.size() >= 1){
      const DayaBay::SimReadout* simReadout = *simReadoutIter;
      DayaBay::Detector det(simReadout->readout()->detector());
      context.SetSite(det.site());
      context.SetDetId(det.detectorId());

      if(det.detectorId() == DetectorId::kAD1
         || det.detectorId() == DetectorId::kAD2
         || det.detectorId() == DetectorId::kAD3
         || det.detectorId() == DetectorId::kAD4
         || det.detectorId() == DetectorId::kIWS
         || det.detectorId() == DetectorId::kOWS){
        const DayaBay::ReadoutPmtCrate* pmtReadout 
          = dynamic_cast<const DayaBay::ReadoutPmtCrate*>(simReadout->readout());
        if(!pmtReadout){
          error() << "Readout from detector " << det 
                  << " is not a PMT Crate readout." << endreq;
          return StatusCode::FAILURE;
        }
        DayaBay::ReadoutPmtCrate* readout 
          = new DayaBay::ReadoutPmtCrate(*pmtReadout);
        readout->setHeader(roHeader);
        roHeader->setReadout(readout);
      }else if(det.detectorId() == DetectorId::kRPC){
        const DayaBay::ReadoutRpcCrate* rpcReadout 
          = dynamic_cast<const DayaBay::ReadoutRpcCrate*>(simReadout->readout());
        if(!rpcReadout){
          error() << "Readout from detector " << det 
                  << " is not an RPC Crate readout." << endreq;
          return StatusCode::FAILURE;
        }
        DayaBay::ReadoutRpcCrate* readout 
          = new DayaBay::ReadoutRpcCrate(*rpcReadout);
        readout->setHeader(roHeader);
        roHeader->setReadout(readout);
      }else{
        error() << "Unknown detector " << det 
                << " for readout " << simReadout << endreq;
        return StatusCode::FAILURE;
      }
      if(simReadouts.size() > 1){
        warning() << "Simulation produced " << simReadouts.size()
                  <<" readouts, ignoring all but first."  << endreq;
      }
    }else{
      roHeader->setReadout(0);
      info() << "No Readouts produced this cycle." << endreq;
      context.SetSite(Site::kUnknown);
      context.SetDetId(DetectorId::kUnknown);
    }
    roHeader->setContext(context);
    return StatusCode::SUCCESS;
}

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