DecayConvolution.h

Go to the documentation of this file.

#ifndef TTTRLIB_FSCONV_H
#define TTTRLIB_FSCONV_H
 
#include <cmath>  /* std::ceil */
#include <numeric> /* std::accumulate */
#include <iostream>
#include <vector>
#include <algorithm> /* std::max */
#include <string.h> /* strcmp */
 
#if defined(_MSC_VER)
  /* Microsoft C/C++-compatible compiler */
  #include <intrin.h>
  #include <immintrin.h>
#endif
 
#if defined(__AVX__)
  #if (defined(__GNUC__) || defined(__clang__))
    #include <immintrin.h>
    #include <avxintrin.h>
  #endif
  #if !defined(__FMA__)
    #define __FMA__ 1
  #endif
#endif 
 
 
void rescale(double *fit, double *decay, double *scale, int start, int stop);
 
 
void rescale_w(double *fit, double *decay, double *w_sq, double *scale, int start, int stop);
 
 
void rescale_w_bg(double *fit, double *decay, double *w_sq, double bg, double *scale, int start, int stop);
 
 
void fconv(double *fit, double *x, double *lamp, int numexp, int start, int stop, double dt=0.05);
 
 
void fconv_avx(double *fit, double *x, double *lamp, int numexp, int start, int stop, double dt=0.05);
 
 
void fconv_per(
        double *fit, double *x, double *lamp, int numexp, int start, int stop,
        int n_points, double period, double dt=0.05
);
void fconv_per_avx(
        double *fit, double *x, double *lamp, int numexp, int start, int stop,
        int n_points, double period, double dt=0.05
);
 
 
void fconv_per_cs(double *fit, double *x, double *lamp, int numexp, int stop,
                  int n_points, double period, int conv_stop, double dt);
 
 
void fconv_ref(double *fit, double *x, double *lamp, int numexp, int start, int stop, double tauref, double dt=0.05);
 
 
void sconv(double *fit, double *p, double *lamp, int start, int stop);
 
 
void shift_lamp(double *lampsh, double *lamp, double ts, int n_points, double out_value=0.0);
 
 
void add_pile_up_to_model(
        double* model, int n_model,
        double* data, int n_data,
        double repetition_rate,
        double instrument_dead_time,
        double measurement_time,
        std::string pile_up_model = "coates",
        int start = 0,
        int stop = -1
);
 
 
void discriminate_small_amplitudes(
    double* lifetime_spectrum, int n_lifetime_spectrum,
    double amplitude_threshold
);
 
 
void fconv_per_cs_time_axis(
    double *model, int n_model,
    double *time_axis, int n_time_axis,
    double *irf, int n_irf,
    double *lifetime_spectrum, int n_lifetime_spectrum,
    int convolution_start = 0,
    int convolution_stop = -1,
    double period = 100.0
);
 
 
void fconv_cs_time_axis(
    double *inplace_output, int n_output,
    double *time_axis, int n_time_axis,
    double *instrument_response_function, int n_instrument_response_function,
    double *lifetime_spectrum, int n_lifetime_spectrum,
    int convolution_start = 0,
    int convolution_stop = -1
);
 
 
 
#endif //TTTRLIB_FSCONV_H