TimeHistogrammer.hpp

Go to the documentation of this file.

 
#ifndef H_TIMEHISTOGRAMMER
#define H_TIMEHISTOGRAMMER
 
 
 
#include "result_collection_support.hpp"
 
#include "DetectorEvent.hpp"
#include "DAQChannels.hpp"
#include "StageDetectorEvent.hpp"
 
#include <deque>
 
 
namespace pottu {
 
 
    template <class SelectorTrigger,class SelectorTarget>
    class TimeHistogrammer : public StageDetectorEvent {
    public:
 
 
#ifndef POTTU_ROOT_SUPPORT
    
    TimeHistogrammer( SelectorTrigger selTrigger, SelectorTarget selTarget, _histogram_1d_type *histogram )
        : _ctxh( ContextBase::getActive().createHandle("TimeHistogrammer") ),
          _selTrigger(selTrigger), _selTarget(selTarget),
          _h(histogram),
          _low( _h->_binner.low() ),
          _high( _h->_binner.high() )
        {
        if( !_h )
            throw std::runtime_error( "No histogram given" );
        if( _low >= _high )
            throw std::runtime_error( "Bad range for histogram." );
        }
#endif
 
    TimeHistogrammer( SelectorTrigger selTrigger, SelectorTarget selTarget,
              _result_tree_type tree,
              const std::string &histname, double low, double high, uint32_t binc )
        : _ctxh( ContextBase::getActive().createHandle("TimeHistogrammer") ),
          _selTrigger(selTrigger), _selTarget(selTarget),
          _h( createHistogram1D( tree, histname.c_str(), low, high, binc ) ),
          _low( low ),
          _high( high )
        {
        if( !_h )
            throw std::runtime_error( "No histogram given" );
        }
    
    
    
    //template <class Cont>
    virtual void process( std::vector<DetectorEvent> &data ) noexcept {
        /*
        static_assert( std::is_same<typename Cont::value_type,DetectorEvent>::value,
               "TimeHistogrammer requires DetectorEvent as a datapacket value_type" );
        */
        for( const auto &event : data ) {
        if( event.time < _tsPrev ) {
            // Time order error! Clearing the queues
            _queueTrigger.clear();
            _queueTarget.clear();
        }
        
        if( _selTrigger(event) ) {
 
            // Popping too old items from targets
            while( _queueTarget.size() && _queueTarget.front()-event.time < _low )
            _queueTarget.pop_front();
            
            for( const auto &tTarget : _queueTarget )
            fill( _h, tTarget-event.time );
            
            _queueTrigger.push_back( event.time );
            
        }
 
        if( _selTarget(event) ) {
 
            // Popping too old items from triggers
            while( _queueTrigger.size() && event.time - _queueTrigger.front() > _high )
            _queueTrigger.pop_front();
 
            for( const auto &tTrigger : _queueTrigger )
            fill( _h, event.time-tTrigger );
 
            _queueTarget.push_back( event.time );
            
        }
        
        _tsPrev = event.time;
            
        }
    }
 
 
    ContextHandle               *_ctxh;
 
    SelectorTrigger              _selTrigger;
    SelectorTarget               _selTarget;
    
    _histogram_1d_type          *_h;
    double                       _low;
    double                       _high;
    
    int64_t                      _tsPrev{0};
    std::deque<int64_t>          _queueTrigger;
    std::deque<int64_t>          _queueTarget;
    };
 
 
 
}
 
 
#endif