rippledetector.hpp

#pragma once
#include <memory>
#include <string>
 
#include "dsp/algorithms.hpp"
#include "eventdata/eventdata.hpp"
#include "iprocessor.hpp"
#include "options/options.hpp"
#include "options/units.hpp"
#include "timeseriesdata/timeseriesdata.hpp"
 
class RippleDetector : public IProcessor {
   public:
    RippleDetector();
    void CreatePorts() override;
    void CompleteStreamInfo() override;
 
    void Preprocess(ProcessingContext& context) override;
    void Process(ProcessingContext& context) override;
    void Postprocess(ProcessingContext& context) override;
 
   protected:
    double compute_value(TimeSeriesType<double>::Data* data_in, unsigned int sample);

    PortIn<TimeSeriesType<double>>* data_in_port_;
    PortOut<EventType>* event_out_port_;
    PortOut<TimeSeriesType<double>>* stats_out_port_;
 
    ProducerState<double>* threshold_;
    ProducerState<double>* signal_mean_;
    ProducerState<double>* signal_dev_;
    BroadcasterState<bool>* ripple_;
    StaticState<double>* threshold_dev_;
    StaticState<double>* detection_lockout_time_;
    FollowerState<bool>* detection_enabled_;
    StaticState<bool>* stream_events_;
    StaticState<double>* smooth_time_;
    StaticState<bool>* stats_out_;
 
    bool stats_out_enabled_;
    std::uint64_t stats_nsamples_;
    std::uint64_t block_;
    std::uint64_t burn_in_;
    double sample_rate_;
    double acc_;
    std::unique_ptr<dsp::algorithms::RunningMeanMAD> running_statistics_;
    std::unique_ptr<dsp::algorithms::ThresholdCrosser> threshold_detector_;
 
   public:
    const std::vector<std::string> STATS_LABEL = {"statistics", "threshold", "deviation"};
    const std::string THRESHOLD_DEV = "threshold_dev";
    const std::string SMOOTH_TIME = "smooth_time";
    const std::string DETECTION_LOCKOUT_TIME = "analysis_lockout_time";
    const std::string STREAM_EVENTS = "stream_events";
    const std::string STREAM_STATISTICS = "stream_statistics";
 
   protected:
    options::Double initial_threshold_dev_{6.};
    options::Measurement<double, false> initial_smooth_time_{10., "second",
                                                             options::positive<double>(true)};

    options::Measurement<double, false> initial_detection_lockout_time_{
        30., "ms", options::positive<double>(true)};

    options::Bool default_stream_events_{true};
    options::Bool initial_stats_out_{true};
    options::Measurement<double, false> stats_buffer_size_{0.5, "second",
                                                           options::positive<double>(true)};

    options::Value<unsigned int, false> stats_downsample_factor_{
        1, options::positive<unsigned int>(true)};

    options::Bool use_power_{true};
};