DecayRoutines.h

Go to the documentation of this file.

 
#ifndef IMPBFF_DECAYROUTINES_H
#define IMPBFF_DECAYROUTINES_H
 
#include <IMP/bff/bff_config.h>
 
#include <cmath>  /* std::ceil */
#include <numeric> /* std::accumulate */
#include <iostream>
#include <vector>
#include <algorithm> /* std::max */
#include <string.h> /* strcmp */
 
#if defined(__AVX__)
 #if defined(_MSC_VER)
   /* Microsoft C/C++-compatible compiler */
   #include <intrin.h>
   #include <immintrin.h>
 #endif
 #if (defined(__GNUC__) || defined(__clang__))
   #include <immintrin.h>
 #endif
 #if !defined(__FMA__)
   #define __FMA__ 1
 #endif
#endif //__AVX__
 
IMPBFF_BEGIN_NAMESPACE
 
    return (n + (a % n)) % n;
}
 
IMPBFFEXPORT void decay_rescale(double *fit, double *decay, double *scale, int start, int stop);
 
 
IMPBFFEXPORT void decay_rescale_w(double *fit, double *decay, double *w_sq, double *scale, int start, int stop);
 
IMPBFFEXPORT void decay_rescale_w_bg(double *fit, double *decay, double *e_sq, double bg, double *scale, int start, int stop);
 
 
IMPBFFEXPORT void decay_fconv(double *fit, double *x, double *lamp, int numexp, int start, int stop, double dt=0.05);
 
IMPBFFEXPORT void decay_fconv_avx(double *fit, double *x, double *lamp, int numexp, int start, int stop, double dt=0.05);
 
IMPBFFEXPORT void decay_fconv_per(
        double *fit, double *x, double *lamp, int numexp, int start, int stop,
        int n_points, double period, double dt=0.05
);
 
IMPBFFEXPORT void decay_fconv_per_avx(
        double *fit, double *x, double *lamp, int numexp, int start, int stop,
        int n_points, double period, double dt=0.05
);
 
IMPBFFEXPORT void decay_fconv_per_cs(double *fit, double *x, double *lamp, int numexp, int stop,
                  int n_points, double period, int conv_stop, double dt);
 
IMPBFFEXPORT void decay_fconv_ref(double *fit, double *x, double *lamp, int numexp, int start, int stop, double tauref, double dt=0.05);
 
IMPBFFEXPORT void decay_sconv(double *fit, double *p, double *lamp, int start, int stop);
 
IMPBFFEXPORT void decay_shift_lamp(double *lampsh, double *lamp, double ts, int n_points, double out_value=0.0);
 
IMPBFFEXPORT void decay_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
);
 
IMPBFFEXPORT void discriminate_small_amplitudes(
        double* lifetime_spectrum, int n_lifetime_spectrum,
        double amplitude_threshold
);
 
IMPBFFEXPORT void decay_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
);
 
IMPBFFEXPORT void decay_fconv_cs_time_axis(
        double *inplace_output, int n_output,
        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
);
 
 
IMPBFF_END_NAMESPACE
 
#endif //IMPBFF_DECAYROUTINES_H