RateHistogrammer.hpp

Go to the documentation of this file.

 
#ifndef H_RATEHISTOGRAMMER
#define H_RATEHISTOGRAMMER
 
 
#include "result_collection_support.hpp"
 
#include "DetectorEvent.hpp"
#include "ContextBase.hpp"
#include "StageDetectorEvent.hpp"
 
 
namespace pottu {
 
    template <class Selector>
    class RateHistogrammer : public StageDetectorEvent {
    public:
 
    RateHistogrammer( _result_tree_type results, Selector selector,
              const std::string &name,
              uint64_t divider=60*100000000ULL, uint64_t low=0, uint64_t high=30ULL*24ULL*60ULL,
              uint64_t binCount=30UL*24UL*60ULL )
        : _ctxh( ContextBase::getActive().createHandle( fmt::format("RateHistogrammer_{}", name ) ) ),
          _selector(selector),
          _divider(divider)
        {
        _h = createHistogram1D( results, name.c_str(), low, high, binCount );
        }
 
    RateHistogrammer( _result_tree_type results, Selector selector,
              const std::string &name,
              uint64_t divider=60*100000000ULL, uint64_t low=0,
              uint64_t binCount=30UL*24UL*60ULL )
        : _ctxh( ContextBase::getActive().createHandle( fmt::format("RateHistogrammer_{}", name ) ) ),
          _selector(selector),
          _divider(divider)
        {
        _h = createHistogram1D( results, name.c_str(), low, low+binCount, binCount );
        }
 
 
    virtual void process( std::vector<DetectorEvent> &data ) noexcept {
        auto proci = _ctxh->processInstance();
 
        for( auto &event : data ) {
        if( _selector(event) ) {
            fill( _h, event.time / _divider );
        }
        }
    }
 
    const Selector &getSelector() const noexcept { return _selector; }
 
    Selector &getSelector() noexcept { return _selector; }
    
    private:
    ContextHandle                *_ctxh;
 
    Selector                      _selector;
    uint64_t                      _divider;
    _histogram_1d_type           *_h;
 
 
    };
 
 
 
}
 
 
#endif