triggeralgs::TCMakerMichelElectronAlgorithm Class Reference

#include <TCMakerMichelElectronAlgorithm.hpp>

Inheritance diagram for triggeralgs::TCMakerMichelElectronAlgorithm:

Collaboration diagram for triggeralgs::TCMakerMichelElectronAlgorithm:

Classes
class	Window

Public Member Functions
void	process (const TriggerActivity &, std::vector< TriggerCandidate > &)
	The function that gets call when there is a new activity.
void	configure (const nlohmann::json &config)
Public Member Functions inherited from triggeralgs::TriggerCandidateMaker
virtual	~TriggerCandidateMaker ()=default
void	operator() (const TriggerActivity &input_ta, std::vector< TriggerCandidate > &output_tc)
virtual bool	preprocess (const TriggerActivity &)
	TA pre-processing/filtering.
virtual void	postprocess (std::vector< TriggerCandidate > &output_tc)
	Post-processing/filtering of the TCs, e.g. prescale.
virtual void	flush (timestamp_t, std::vector< TriggerCandidate > &)

Private Member Functions
TriggerCandidate	construct_tc () const
bool	check_adjacency () const
void	add_window_to_record (Window window)
void	dump_window_record ()

Private Attributes
Window	m_current_window
uint64_t	m_activity_count = 0
bool	m_trigger_on_adc = false
bool	m_trigger_on_n_channels = false
uint32_t	m_adc_threshold = 1200000
uint16_t	m_n_channels_threshold = 600
timestamp_t	m_window_length = 80000
int	tc_number = 0
std::vector< Window >	m_window_record

Additional Inherited Members
Public Attributes inherited from triggeralgs::TriggerCandidateMaker
std::atomic< uint64_t >	m_data_vs_system_time = 0
std::atomic< uint64_t >	m_initial_offset = 0
uint64_t	m_prescale = 1
	Configurable prescale factor.
uint64_t	m_tc_count = 0
	TC made count for prescaling.
TriggerCandidate::Type	m_tc_type_out = TriggerCandidate::Type::kUnknown
	Configurable TC type output.

Detailed Description

Definition at line 20 of file TCMakerMichelElectronAlgorithm.hpp.

Member Function Documentation

add_window_to_record()

void TCMakerMichelElectronAlgorithm::add_window_to_record ( Window window )

private

Definition at line 162 of file TCMakerMichelElectronAlgorithm.cpp.

{
  m_window_record.push_back(window);
  return;
}

check_adjacency()

bool TCMakerMichelElectronAlgorithm::check_adjacency ( ) const

private

Definition at line 154 of file TCMakerMichelElectronAlgorithm.cpp.

{
  // FIX ME: An adjacency check on the channels which have hits.
  return true;
}

configure()

void TCMakerMichelElectronAlgorithm::configure ( const nlohmann::json & config )

virtual

Reimplemented from triggeralgs::TriggerCandidateMaker.

Definition at line 101 of file TCMakerMichelElectronAlgorithm.cpp.

{
  TriggerCandidateMaker::configure(config);
 
  // FIX ME: Use some schema here. Also can't work out how to pass booleans.
  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 (config.contains("channel_map")) m_channel_map = config["channel_map"];
  }
  if (m_trigger_on_adc && m_trigger_on_n_channels) {
    TLOG_DEBUG(TLVL_VERY_IMPORTANT) << "[TCM:ME] 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_DEBUG_LOW) << "[TCM:ME] Both trigger flags are false. Passing TAs through 1:1.";
  }
 
  return;
}

construct_tc()

TriggerCandidate TCMakerMichelElectronAlgorithm::construct_tc ( ) const

private

Definition at line 131 of file TCMakerMichelElectronAlgorithm.cpp.

{
  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.time_end;
  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::kMichelElectron;
 
  // 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);
  }
 
  return tc;
}

dump_window_record()

void TCMakerMichelElectronAlgorithm::dump_window_record ( )

private

Definition at line 169 of file TCMakerMichelElectronAlgorithm.cpp.

{
  // 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;
}

process()

void TCMakerMichelElectronAlgorithm::process ( const TriggerActivity & activity, std::vector< TriggerCandidate > & output_tc )

virtual

The function that gets call when there is a new activity.

Implements triggeralgs::TriggerCandidateMaker.

Definition at line 25 of file TCMakerMichelElectronAlgorithm.cpp.

{
 
  std::vector<TriggerActivity::TriggerActivityData> ta_list = { static_cast<TriggerActivity::TriggerActivityData>(
    activity) };
 
  // The first time process() is called, reset window object.
  if (m_current_window.is_empty()) {
    m_current_window.reset(activity);
    m_activity_count++;
    // Trivial TC Logic:
    // If the request has been made to not trigger on number of channels or
    // total adc, simply construct a trigger candidate from any single activity.
    if ((!m_trigger_on_adc) && (!m_trigger_on_n_channels)) {
 
      // add_window_to_record(m_current_window);
      // dump_window_record();
      TLOG_DEBUG(TLVL_DEBUG_LOW) << "[TCM:ME] Constructing TC.";
 
      TriggerCandidate tc = construct_tc();
      output_tc.push_back(tc);
 
      // Clear the current window (only has a single TA in it)
      m_current_window.clear();
    }
    return;
  }
 
  // FIX ME: Only want to call this if running in debug mode.
  // add_window_to_record(m_current_window);
 
  // 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.
  if ((activity.time_start - m_current_window.time_start) < m_window_length) {
    TLOG_DEBUG(TLVL_DEBUG_HIGH) << "[TCM:ME] Window not yet complete, adding the activity to the window.";
    m_current_window.add(activity);
  }
  // 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.
  else if (m_current_window.adc_integral > m_adc_threshold && m_trigger_on_adc) {
    TLOG_DEBUG(TLVL_DEBUG_LOW) << "[TCM:ME] ADC integral in window is greater than specified threshold.";
    TriggerCandidate tc = construct_tc();
 
    output_tc.push_back(tc);
    TLOG_DEBUG(TLVL_DEBUG_HIGH) << "[TCM:ME] Resetting window with activity.";
    m_current_window.reset(activity);
  }
  // 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);
    TLOG_DEBUG(TLVL_DEBUG_INFO) << "[TCM:ME] Should not see this!";
  }
  // If it is not, move the window along.
  else {
    TLOG_DEBUG(TLVL_DEBUG_HIGH) << "[TCM:ME] Window is at required length but specified threshold not met, shifting window along.";
    m_current_window.move(activity, m_window_length);
  }
 
  //TLOG_DEBUG(TLVL_DEBUG_ALL) << "[TCM:ME] " m_current_window;
 
  m_activity_count++;
 
  //  if(m_activity_count % 500 == 0) dump_window_record();
 
  return;
}

Member Data Documentation

m_activity_count

uint64_t triggeralgs::TCMakerMichelElectronAlgorithm::m_activity_count = 0

private

Definition at line 128 of file TCMakerMichelElectronAlgorithm.hpp.

m_adc_threshold

uint32_t triggeralgs::TCMakerMichelElectronAlgorithm::m_adc_threshold = 1200000

private

Definition at line 138 of file TCMakerMichelElectronAlgorithm.hpp.

m_current_window

Window triggeralgs::TCMakerMichelElectronAlgorithm::m_current_window

private

Definition at line 127 of file TCMakerMichelElectronAlgorithm.hpp.

m_n_channels_threshold

uint16_t triggeralgs::TCMakerMichelElectronAlgorithm::m_n_channels_threshold = 600

private

Definition at line 139 of file TCMakerMichelElectronAlgorithm.hpp.

m_trigger_on_adc

bool triggeralgs::TCMakerMichelElectronAlgorithm::m_trigger_on_adc = false

private

Definition at line 136 of file TCMakerMichelElectronAlgorithm.hpp.

m_trigger_on_n_channels

bool triggeralgs::TCMakerMichelElectronAlgorithm::m_trigger_on_n_channels = false

private

Definition at line 137 of file TCMakerMichelElectronAlgorithm.hpp.

m_window_length

timestamp_t triggeralgs::TCMakerMichelElectronAlgorithm::m_window_length = 80000

private

Definition at line 140 of file TCMakerMichelElectronAlgorithm.hpp.

m_window_record

std::vector<Window> triggeralgs::TCMakerMichelElectronAlgorithm::m_window_record

private

Definition at line 148 of file TCMakerMichelElectronAlgorithm.hpp.

tc_number

int triggeralgs::TCMakerMichelElectronAlgorithm::tc_number = 0

private

Definition at line 141 of file TCMakerMichelElectronAlgorithm.hpp.

---

The documentation for this class was generated from the following files:

• /github/workspace/dunedaq/sourcecode/triggeralgs/include/triggeralgs/MichelElectron/TCMakerMichelElectronAlgorithm.hpp
• /github/workspace/dunedaq/sourcecode/triggeralgs/src/TCMakerMichelElectronAlgorithm.cpp