TCMakerHorizontalMuonAlgorithm.cpp

Go to the documentation of this file.

 
#include "triggeralgs/HorizontalMuon/TCMakerHorizontalMuonAlgorithm.hpp"
 
#include "TRACE/trace.h"
#define TRACE_NAME "TCMakerHorizontalMuonAlgorithm"
 
#include <math.h>
#include <vector>
 
using namespace triggeralgs;
 
using Logging::TLVL_DEBUG_ALL;
using Logging::TLVL_DEBUG_HIGH;
using Logging::TLVL_DEBUG_MEDIUM;
 
void
TCMakerHorizontalMuonAlgorithm::process(const TriggerActivity& activity,
                                                std::vector<TriggerCandidate>& output_tc)
{
 
  std::vector<TriggerActivity::TriggerActivityData> ta_list = { static_cast<TriggerActivity::TriggerActivityData>(
    activity) };
 
  // The first time operator is called, reset window object.
  if (m_current_window.is_empty()) {
    m_current_window.reset(activity);
    m_activity_count++;
 
    // TriggerCandidate tc = construct_tc();
    // output_tc.push_back(tc);
 
    // m_current_window.clear();
    // return;
  }
 
  // If the difference between the current TA's start time and the start of the window
  // is less than the specified window size, add the TA to the window.
  else if ((activity.time_start - m_current_window.time_start) < m_window_length) {
    TLOG_DEBUG(TLVL_DEBUG_HIGH) << "[TCM:HM] Window not yet complete, adding the activity to the window.";
    m_current_window.add(activity);
  }
  // If it is not, move the window along.
  else {
    TLOG_DEBUG(TLVL_DEBUG_ALL)
      << "[TCM:HM] TAWindow is at required length but specified threshold not met, shifting window along.";
    m_current_window.move(activity, m_window_length);
  }
 
  // If the addition of the current TA to the window would make it longer
  // than the specified window length, don't add it but check whether the sum of all adc in
  // the existing window is above the specified threshold. If it is, and we are triggering on ADC,
  // make a TA and start a fresh window with the current TP.
  if (m_current_window.adc_integral > m_adc_threshold && m_trigger_on_adc) {
    // TLOG_DEBUG(TRACE_NAME) << "ADC integral in window is greater than specified threshold.";
    TLOG_DEBUG(TLVL_DEBUG_MEDIUM) << "[TCM:HM] m_current_window.adc_integral " << m_current_window.adc_integral
                                  << " - m_adc_threshold " << m_adc_threshold;
    tc_number++;
    TriggerCandidate tc = construct_tc();
    TLOG_DEBUG(TLVL_DEBUG_MEDIUM) << "[TCM:HM] tc.time_start=" << tc.time_start << " tc.time_end=" << tc.time_end
                                  << " len(tc.inputs) " << tc.inputs.size();
 
    for (const auto& ta : tc.inputs) {
      TLOG_DEBUG(TLVL_DEBUG_ALL) << "[TCM:HM] [TA] ta.time_start=" << ta.time_start << " ta.time_end=" << ta.time_end
                                 << " ta.adc_integral=" << ta.adc_integral;
    }
 
    output_tc.push_back(tc);
    // m_current_window.reset(activity);
    m_current_window.clear();
  }
 
  // If the addition of the current TA to the window would make it longer
  // than the specified window length, don't add it but check whether the number of hit channels in
  // the existing window is above the specified threshold. If it is, and we are triggering on channels,
  // make a TC and start a fresh window with the current TA.
  else if (m_current_window.n_channels_hit() > m_n_channels_threshold && m_trigger_on_n_channels) {
    tc_number++;
    output_tc.push_back(construct_tc());
 
    // m_current_window.reset(activity);
    m_current_window.clear();
  }
 
  // // If it is not, move the window along.
  // else {
  //   // TLOG_DEBUG(TRACE_NAME) << "TAWindow is at required length but specified threshold not met."
  //   //                        << "shifting window along.";
  //   m_current_window.move(activity, m_window_length);
  // }
 
  m_activity_count++;
  return;
}
 
void
TCMakerHorizontalMuonAlgorithm::configure(const nlohmann::json& config)
{
  TriggerCandidateMaker::configure(config);
 
  if (config.is_object()) {
    if (config.contains("trigger_on_adc"))
      m_trigger_on_adc = config["trigger_on_adc"];
    if (config.contains("trigger_on_n_channels"))
      m_trigger_on_n_channels = config["trigger_on_n_channels"];
    if (config.contains("adc_threshold"))
      m_adc_threshold = config["adc_threshold"];
    if (config.contains("n_channels_threshold"))
      m_n_channels_threshold = config["n_channels_threshold"];
    if (config.contains("window_length"))
      m_window_length = config["window_length"];
  }
  if (m_trigger_on_adc && m_trigger_on_n_channels) {
    TLOG_DEBUG(TLVL_VERY_IMPORTANT) << "[TCM:HM] Triggering on ADC count and number of channels is not supported.";
    throw BadConfiguration(ERS_HERE, TRACE_NAME);
  }
  if (!m_trigger_on_adc && !m_trigger_on_n_channels) {
    TLOG_DEBUG(TLVL_VERY_IMPORTANT) << "[TCM:HM] Not triggering! All trigger flags are false!";
    throw BadConfiguration(ERS_HERE, TRACE_NAME);
  }
 
  return;
}
 
TriggerCandidate
TCMakerHorizontalMuonAlgorithm::construct_tc() const
{
  TriggerActivity latest_ta_in_window = m_current_window.inputs.back();
 
  TriggerCandidate tc;
  tc.time_start = m_current_window.time_start;
  tc.time_end = latest_ta_in_window.inputs.back().time_start;
  tc.time_candidate = m_current_window.time_start;
  tc.detid = latest_ta_in_window.detid;
  tc.type = m_tc_type_out;
  tc.algorithm = TriggerCandidate::Algorithm::kHorizontalMuon;
 
  // Take the list of triggeralgs::TriggerActivity in the current
  // window and convert them (implicitly) to detdataformats'
  // TriggerActivityData, which is the base class of TriggerActivity
  for (auto& ta : m_current_window.inputs) {
    tc.inputs.push_back(ta);
    if (ta.time_end > tc.time_end) {
      tc.time_end = ta.time_end;
    }
  }
 
  return tc;
}
 
bool
TCMakerHorizontalMuonAlgorithm::check_adjacency() const
{
  // FIX ME: An adjacency check on the channels which have hits.
  return true;
}
 
// Functions below this line are for debugging purposes.
void
TCMakerHorizontalMuonAlgorithm::add_window_to_record(TAWindow window)
{
  m_window_record.push_back(window);
  return;
}
 
void
TCMakerHorizontalMuonAlgorithm::dump_window_record()
{
  // FIX ME: Need to index this outfile in the name by detid or something similar.
  std::ofstream outfile;
  outfile.open("window_record_tcm.csv", std::ios_base::app);
 
  for (auto window : m_window_record) {
    outfile << window.time_start << ",";
    outfile << window.inputs.back().time_start << ",";
    outfile << window.inputs.back().time_start - window.time_start << ",";
    outfile << window.adc_integral << ",";
    outfile << window.n_channels_hit() << ",";
    outfile << window.inputs.size() << std::endl;
  }
 
  outfile.close();
 
  m_window_record.clear();
 
  return;
}
 
/*
void
TCMakerHorizontalMuonAlgorithm::flush(timestamp_t, std::vector<TriggerCandidate>& output_tc)
{
  // Check the status of the current window, construct TC if conditions are met. Regardless
  // of whether the conditions are met, reset the window.
  if(m_current_window.adc_integral > m_adc_threshold && m_trigger_on_adc){
  //else if(m_current_window.adc_integral > m_conf.adc_threshold && m_conf.trigger_on_adc){
    //TLOG_DEBUG(TRACE_NAME) << "ADC integral in window is greater than specified threshold.";
    output_tc.push_back(construct_tc());
  }
  else if(m_current_window.n_channels_hit() > m_n_channels_threshold && m_trigger_on_n_channels){
  //else if(m_current_window.n_channels_hit() > m_conf.n_channels_threshold && m_conf.trigger_on_n_channels){
    //TLOG_DEBUG(TRACE_NAME) << "Number of channels hit in the window is greater than specified threshold.";
    output_tc.push_back(construct_tc());
  }
 
  //TLOG_DEBUG(TRACE_NAME) << "Clearing the current window, on the arrival of the next input_tp, the window will be
reset."; m_current_window.clear();
 
  return;
}*/
 
REGISTER_TRIGGER_CANDIDATE_MAKER(TRACE_NAME, TCMakerHorizontalMuonAlgorithm)