CorrelatorCurve.h

Go to the documentation of this file.

#ifndef TTTRLIB_CORRELATORCURVE_H
#define TTTRLIB_CORRELATORCURVE_H
 
#include <iostream>
#include <vector>
#include <algorithm>  /* std::max */
#include <cmath> /* pow */
#include <cstdint> // include this header for uint64_t
 
 
 
struct CorrelationCurveSettings{
 
    double macro_time_duration = 1.0;
 
    unsigned int n_casc = 25;
 
    unsigned int n_bins = 17;
 
    unsigned int get_ncorr() const{
        return n_casc * n_bins + 1;
    }
 
    std::string correlation_method;
 
};
 
 
class CorrelatorCurve{
 
    friend class Correlator;
 
private:
 
    void update_axis();
 
protected:
 
    std::vector<unsigned long long> x_axis;
 
    std::vector<double> correlation;
 
    std::vector<double> corr_normalized;
 
    void resize(size_t n){
        x_axis.resize(n);
        correlation.resize(n);
        corr_normalized.resize(n);
    }
 
    void clear(){
        std::fill(correlation.begin(), correlation.end(), 0.0);
    }
 
public:
 
    CorrelationCurveSettings settings;
 
    size_t size(){
        return x_axis.size();
    }
 
    void get_x_axis(double **output, int *n_output);
 
    void set_x_axis(std::vector<long long unsigned int> input);
 
    void set_n_bins(int v) {
        settings.n_bins = std::max(1, v);
        update_axis();
    }
 
    unsigned int get_n_bins() {
        return settings.n_bins;
    }
 
    void set_n_casc(int v) {
        settings.n_casc = std::max(1, v);
        update_axis();
    }
 
    unsigned int get_n_casc() const {
        return settings.n_casc;
    }
 
    void get_corr(double** output, int* n_output);
 
    void get_corr_normalized(double** output, int* n_output);
 
};
 
 
#endif //TTTRLIB_CORRELATORCURVE_H