use_functions_wo_measurement.cc

#include <iomanip>
#include <TRandom.h>
 
// example how to use compiled functions in a macro
// an example how to use the implemented fit functions and other functions in python is given in cosmics-fit.ipynb
void use_functions_wo_measurement() {
 
    // generate 15 example data points with a bin size of 1.8 units
    int n_entries = 15;
    double bin_size = 1.8;
 
    // for this example an array is filled with 10/(1+(x-x0)^2) + noise:
    double* test_arr = new double[n_entries];
    double* x = new double[n_entries];
    int x0 = static_cast<int>(n_entries / 2);
 
    cout << "\noriginal: ";
    for (int i = 0; i < n_entries; i++) {
        x[i] = static_cast<double>(i) * bin_size;
        test_arr[i] = gRandom->Gaus(0, .5) + 10. / (1. + static_cast<double>((i - x0) * (i - x0)) * (bin_size * bin_size));
        cout << setprecision(2) << test_arr[i] << "\t";
    }
 
    // create graph of the original "data"
    auto orig = new TGraph(n_entries, x, test_arr);
    orig->SetTitle("original data; x title [arb.]; y title [arb.]");
 
    // now the array is smoothed with a gauss kernel (method = 2) with a sigma of 1.2 using ReadRun::SmoothArray()
    Filters::SmoothArray(test_arr, n_entries, 1.2, "Gaus", bin_size); // sigma = 1.2, method = "Gaus"
 
    // print content of smoothed array
    cout << "\n\nsmoothed: ";
    for (int i = 0; i < n_entries; i++) cout << setprecision(2) << test_arr[i] << "\t";
 
    auto smoothed = new TGraph(n_entries, x, test_arr); 
    smoothed->SetTitle("smoothed data");
    smoothed->SetLineColor(Helpers::rcolor(1));
    cout << endl;
 
    // you can also initialze a dummy of the ReadRun class to use its functions. This also works for non-static functions
    Filters::SmoothArray(test_arr, n_entries, 1, "Box"); // sigma = 1, method = 0 -> average over neigboring +-1 bins
 
    cout << "\naveraged: ";
    for (int i = 0; i < n_entries; i++) cout << setprecision(2) << test_arr[i] << "\t";
    cout << endl;
 
    auto more_smoothed = new TGraph(n_entries, x, test_arr);
    more_smoothed->SetTitle("even more smoothed data");
    more_smoothed->SetLineColor(4);
 
    // plot results
    gStyle->SetOptTitle(0); // no title
    orig->Draw("AC+");
    smoothed->Draw("C same");
    more_smoothed->Draw("C same");
    gPad->BuildLegend(.5,.7,.9,.9); // draw legend to coordinates
    gPad->SaveAs("examples/use_functions_wo_measurement.png");
 
    // clean up
    delete[] test_arr;
    delete[] x;
}