Data processing version 2

Data location: /data/Reconstructed_GAGG_v2

Reprocessed only GAGG background data so far.

Data processed with a single config file for all runs, available at
Added average pulse for alpha events
Added average power spectrum for alpha events
Fixed bug in Chi2 module, which prevented to find the correct min-Chi2 for a fraction of the events. However, this bug had a very small effect, in fact the PSD was working pretty well anyway.
Added Chi2 calculation in Fourier complex space, i.e. including the phase information, to identify pileup events. From here on, we can use the chi2 computed on the power spectrum (i.e. without phase information) to perform particle discrimination, and the chi2 computed on the complex space for pileup rejection.
Updated the PSD module to fit the beta-gamma and alpha bands, and store in the output tree a list of output variables that make the application of the PSD much easier.
Implemented the auto-calibration of the beta-gamma band on the 40K peak.
Optimized the frequency range for the calculation of the Chi2. As a result, the PSD cut will change.

Baseline slope, for cutting events sitting on the tail of a previous one: fabs(baseline.slope)<1.e7
Trigger delay, for cutting (mildly) pileup events: triggerdelay.delay>0.45e-6 && triggerdelay.delay<0.47e-6
Chi2 on Fourier-transformed complex pulse for cutting (strongly) pileup events: chi2_freqsingleevent.Chi2<0.004
Energy calibrated on the 40K peak: gammacalibration.energy
PSD variable, defined as gaggpsd.psd and equal to chi2_poweraverage.amplitude/integral. The new PSD threshold for cutting beta-gamma vs alpha events could be e.g. at gaggpsd.psd=0.53. Alternatively, we can use the two variables described below.
Beta-Gamma flag (gaggpsd.betagamma), defined as true for events within +-4 sigma from the mean of the beta-gamma band, and as false otherwise
Distance from the mean of the beta-gamma band, normalized to the corresponding sigma. This variable is useful for re-defining the PSD cut more easily. For example, if we want to select the beta-gamma events within +-3 sigma from the mean of the band, we can apply the cut fabs(gaggpsd.nSigmaBetaGamma)<3