CorrelatorPhotonStream.h

Go to the documentation of this file.

#ifndef TTTRLIB_CORRELATORPHOTONSTREAM_H
#define TTTRLIB_CORRELATORPHOTONSTREAM_H
 
#include <vector>
#include <numeric>  /* std::accumulate */
 
#include "TTTR.h"
 
 
 
class CorrelatorPhotonStream{
 
    friend class Correlator;
 
protected:
 
    double time_axis_calibration = 1.0;
 
    std::shared_ptr<TTTR> tttr = nullptr;
 
public:
 
    CorrelatorPhotonStream() = default;
 
    CorrelatorPhotonStream(const CorrelatorPhotonStream& a):
    times(a.times), weights(a.weights) {};
 
    ~CorrelatorPhotonStream() = default;
 
    std::vector<unsigned long long> times;
 
    std::vector<double> weights;
 
    bool empty() const{
        return times.empty() && weights.empty();
    }
 
    size_t size() const{
        return times.size();
    }
 
    void clear(){
        times.clear();
        weights.clear();
    }
 
    void resize(size_t n, double x = 1.0){
        times.resize(n);
        weights.resize(n, x);
    }
 
    static void make_fine_times(
            unsigned long long *t, unsigned int n_times,
            unsigned short *tac,
            unsigned int n_tac
    );
 
    void make_fine(
            unsigned short *tac, int n_tac,
            unsigned int number_of_microtime_channels
    );
 
    void set_weights(
        const std::map<short, std::vector<double>>& filter,
        std::vector<unsigned int> micro_times = std::vector<unsigned int>(),
        std::vector<signed char> routing_channels = std::vector<signed char>()
    );
 
    void set_events(
        unsigned long long  *t, int n_t,
        double* weight, int n_weight
    ){
        resize(std::min(n_t, n_weight));
        for(size_t i = 0; i < size();i++){
            times[i] = t[i];
            weights[i] = weight[i];
        }
    }
 
    void coarsen();
 
    unsigned long long dt();
 
    double sum_of_weights();
 
    double mean_count_rate();
 
    void set_time_axis_calibration(double v);
 
    double get_time_axis_calibration() const{
        return time_axis_calibration;
    }
 
    void set_tttr(std::shared_ptr<TTTR> tttr, bool make_fine = false);
 
    std::shared_ptr<TTTR> get_tttr() const {
        return tttr;
    }
 
};
 
#endif //TTTRLIB_CORRELATORPHOTONSTREAM_H