In This Package:

EsDarkDetectorAlg.cc

Go to the documentation of this file.

#include "EsDarkDetectorAlg.h"

#include "Event/SimHeader.h"
#include "Conventions/Detectors.h"
#include "CLHEP/Units/SystemOfUnits.h"

EsDarkDetectorAlg::EsDarkDetectorAlg(const std::string& name, ISvcLocator* pSvcLocator)
  : DybAlgorithm<DayaBay::SimHeader>(name,pSvcLocator),
    m_detectors()
{
    declareProperty("SimLocation",m_simLocation=DayaBay::SimHeaderLocation::Default,
                    "Location in the TES where the input SimHeader is to be found.");
    declareProperty("Detectors",m_detectors,
                    "List of detectors to simulate");
    declareProperty("SimWindow",m_simWindow=0,
                    "Length of the simulation time window");
    declareProperty("SimDeltaT",m_simDeltaT=1.0*CLHEP::second,
                    "Time between each simulation cycle");
    declareProperty("SimStartTime",m_simStartTime=0.0,
                    "Starting timestamp of the simulation");
}

EsDarkDetectorAlg::~EsDarkDetectorAlg()
{
}

StatusCode EsDarkDetectorAlg::initialize()
{
  m_currentTime = m_simStartTime;
  return StatusCode::SUCCESS;
}

StatusCode EsDarkDetectorAlg::execute()
{
    // Create output header 
    DayaBay::SimHeader* simHeader = MakeHeaderObject();
    simHeader->setTimeStamp( m_currentTime );

    // Set the simulation time range
    // Add padding to ElecHeader time window
    TimeStamp earliest = m_currentTime;
    TimeStamp latest = m_currentTime;
    latest.Add(m_simWindow);
    simHeader->setEarliest(earliest);
    simHeader->setLatest(latest);
    Context context;
    context.SetSimFlag(SimFlag::kMC);
    context.SetTimeStamp(earliest);

    // Add the sim hit header
    DayaBay::SimHitHeader* hitHeader = 
      new DayaBay::SimHitHeader(simHeader);
    simHeader->setHits(hitHeader);
    DayaBay::SimHitCollection::hit_container hits;

    bool firstDetector = true;
    // Loop over detectors and add an empty sim hit collection for each detector
    std::vector<std::string>::const_iterator detIter, end = 
      m_detectors.end();
    for(detIter = m_detectors.begin(); detIter != end; detIter++){
      DayaBay::Detector detector 
        = DayaBay::Detector::siteDetPackedFromString(*detIter);
      hitHeader->addHitCollection(
                      new DayaBay::SimHitCollection(hitHeader,detector,hits));

      // Reset context if needed
      if(firstDetector){
        // Keep track of sites/detectors in this simulation for context
        context.SetSite(detector.site());
        context.SetDetId(detector.detectorId());
        firstDetector = false;
      }
      if(context.GetSite() != detector.site()){
        // Simulation contains multiple sites, unset
        context.SetSite(Site::kUnknown);
      }
      if(context.GetDetId() != detector.detectorId()){
        // Simulation contains multiple detectors, unset
        context.SetDetId(DetectorId::kUnknown);
      }
    }
    simHeader->setContext(context);

    m_currentTime.Add(m_simDeltaT / CLHEP::second);
    return StatusCode::SUCCESS;
}

StatusCode EsDarkDetectorAlg::finalize()
{
    return StatusCode::SUCCESS;
}

---