Experimental.cc

Go to the documentation of this file.

#include "Experimental.h"
 
void Experimental::RebinAll(int ngroup, float sigma_noise, unsigned long seed) {
 
    SP *= static_cast<float>(ngroup);
    SP_inv = 1. / SP;
    binNumber /= ngroup;
    float norm = 1. / static_cast<float>(ngroup);
    cout    << "\nRebinning the data to a new sampling rate of " << 1. / SP 
            << " GS/s which corresponds to a bin size of " << SP 
            << " ns and the data now has " << binNumber << " bins" << endl;
 
    #pragma omp parallel for
    for (int j = 0; j < nwf; j++) {
        vector<float> &waveform = rundata[j];
        vector<float> rebinned;
        unsigned long thread_seed = seed + j;
        TRandom3 noise(thread_seed);
 
        for (size_t i = 0; i + ngroup <= waveform.size(); i += ngroup) {
            float sum = 0.;
            for (int k = 0; k < ngroup; k++) {
                sum += waveform[i + k];
            }
            float avg = sum * norm;
 
            if (sigma_noise != 0.) {
                avg += noise.Gaus(0., sigma_noise);
            }
            rebinned.push_back(avg);
        }
        rundata[j] = rebinned;
    }
}
 
void Experimental::DerivativeAll() {
    cout << "\nForward derivative of all waveforms:" << endl;
    binNumber--; // reduce number of bins 
 
    #pragma omp parallel for
    for (int j = 0; j < nwf; j++) {
        vector<float> &waveform = rundata[j];
        vector<float> derivative;
        for (int i = 0; i < binNumber; i++) derivative.push_back(waveform[i + 1] - waveform[i]);
        rundata[j] = derivative;
    }
}