Line data    Source code

       1              : /**
       2              :  * @file NaiveFrugalPedestalSubtractProcessor.cpp
       3              :  *
       4              :  * @copyright This is part of the DUNE DAQ Software Suite, copyright 2020.
       5              :  * Licensing/copyright details are in the COPYING file that you should have
       6              :  * received with this code.
       7              :  */
       8              : 
       9              : #include "tpglibs/NaiveFrugalPedestalSubtractProcessor.hpp"
      10              : 
      11              : namespace tpglibs {
      12              : 
      13           13 : REGISTER_NAIVEPROCESSOR_CREATOR("NaiveFrugalPedestalSubtractProcessor", NaiveFrugalPedestalSubtractProcessor)
      14              : 
      15              : void
      16            0 : NaiveFrugalPedestalSubtractProcessor::configure(const nlohmann::json& config, const int16_t* plane_numbers)
      17              : {
      18              : #ifdef TPGLIBS_ENABLE_STATE_MONITORING
      19              :   m_internal_state_name_registry.register_internal_state("pedestal",
      20              :                                                          std::shared_ptr<naive_array_t>(&m_pedestal, [](auto*) {}));
      21              :   m_internal_state_name_registry.register_internal_state("accum",
      22              :                                                          std::shared_ptr<naive_array_t>(&m_accum, [](auto*) {}));
      23              :   configure_internal_state_collection(config);
      24              : #endif
      25              : 
      26            0 :   m_accum_limit = config["accum_limit"];
      27            0 : }
      28              : 
      29              : NaiveFrugalPedestalSubtractProcessor::naive_array_t
      30            0 : NaiveFrugalPedestalSubtractProcessor::process(const naive_array_t& signal)
      31              : {
      32              : #ifdef TPGLIBS_ENABLE_STATE_MONITORING
      33              :   m_samples++;
      34              :   if (m_collect_internal_state_flag && (m_samples % m_sample_period == 0)) {
      35              :     m_internal_state_buffer_manager.write_to_active_buffer();
      36              :   }
      37              : #endif
      38              : 
      39            0 :   naive_array_t subtracted_signal;
      40            0 :   for (int i = 0; i < 16; i++) {
      41              :     // Increment if above.
      42            0 :     if (signal[i] > m_pedestal[i])
      43            0 :       m_accum[i]++;
      44              : 
      45              :     // Decrement if below.
      46            0 :     if (signal[i] < m_pedestal[i])
      47            0 :       m_accum[i]--;
      48              : 
      49              :     // Increment pedestal if we've hit the top limit.
      50            0 :     if (m_accum[i] > m_accum_limit) {
      51            0 :       m_pedestal[i]++;
      52            0 :       m_accum[i] = 0;
      53              :     }
      54              : 
      55              :     // Decrement pedestal if we've hit the low limit.
      56            0 :     if (m_accum[i] < -1 * m_accum_limit) {
      57            0 :       m_pedestal[i]--;
      58            0 :       m_accum[i] = 0;
      59              :     }
      60              : 
      61            0 :     subtracted_signal[i] = signal[i] - m_pedestal[i];
      62              :   }
      63              : 
      64            0 :   return NaiveProcessor::process(subtracted_signal);
      65              : }
      66              : 
      67              : } // namespace tpglibs