Experimental.cc

Go to the documentation of this file.

#include "Experimental.h"
 
void Experimental::RebinAll(int ngroup, float noise_level, unsigned long seed) {
 
    SP *= static_cast<float>(ngroup);
    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;
 
    for (int j = 0; j < nwf; j++) {
        TH1F* his = Getwf(j);
        his->Rebin(ngroup);
        his->Scale(norm);
 
        if (noise_level != 0.) {
            auto noise = new TRandom3();
            noise->SetSeed(seed);
            for (int i = 1; i <= his->GetNbinsX(); i++) his->SetBinContent(i, his->GetBinContent(i) + noise->Gaus(0, noise_level));
        }
        Helpers::PrintProgressBar(j, nwf);
    }
}
 
void Experimental::DerivativeAll() {
    cout << "\nForward derivative of all waveforms:" << endl;
    for (int j = 0; j < nwf; j++) {
        TH1F* his = Getwf(j);
        double* yvals = Helpers::gety(his);
        for (int i = 1; i <= his->GetNbinsX() - 1; i++) his->SetBinContent(i, yvals[i + 1] - yvals[i]);
        delete[] yvals;
        Helpers::PrintProgressBar(j, nwf);
    }
}