Bfit_Utils: Utilities for baseline optimal fit

tesfdmtools.methods.Bfit_Utils.fitgauss(pltdata, iplt, hdf, channel, freq, spectrum, bins, debug=False)

Fit gaussian function to the optimal fit spectrum

Args:

• pltdata = dictionary for the fit results

• iplt = number of the plot to process

• hdf = HDF5 input file object

• channel = channel number being processed

• freq = frequency number (pixel) being processed

• spectrum = spectrum of the pulse fit parameters

• bins = centers of the spectrum bins

Kwargs:

• debug = if True, plot debug information [default: False]

Returns:

• afit = average fitted position

• fitpars = fit parameters

  • fitpars[0] = Gaussian scale

  • fitpars[1] = Gaussian position

  • fitpars[2] = Gaussian position width

  • fitpars[3-5] = Errors for fitpars[0-2]

tesfdmtools.methods.Bfit_Utils.nfilt(hdf, channel, freq, ptp, ofilter, norm, debug=False, absolute=False, rotate=False, freqcutoff=None, prlen=None, threshold=1000)

Perform optimal filter for baseline records

Args:

• hdf = HDF5 input file object

• channel = channel number being processed

• freq = frequency number (pixel) being processed

• ptp = delta difference within records

• ofilter = the optimal filter

• norm = optimal filter scaling factor

Kwargs:

• debug = if True, plot various fit parameters [default: False]

• absolute = use sqrt(I^2+Q^2) signal [default: False]

• rotate = rotation angle (radians) for I/Q pulse record

• freqcutoff = cutoff frequency (Hz) for optimal filtering [default: None]

• prlen = length of record to process (None = entire record) [default: None]

• threshold = treshold for event detection

Returns:

• ofits = list with optimal fits for noise records

• indx = index of records with the optimal fits

tesfdmtools.methods.Bfit_Utils.optfilt(hdf, channel, freq, indx, base1, base2, noisspec, debug=False, freqcutoff=None, rotate=False, bsec=0.05, absolute=False, prlen=None, avpulsepos=None)

compute optimal filter for pulses

Args:

• hdf = HDF5 input file object

• channel = channel number being processed

• freq = frequency number (pixel) being processed

• indx = index of selected events, to be processed

• base1 = baseline level at start of record

• base2 = baseline level at end of record

• noisspec = noise spectrum

Kwargs:

• debug = if True, plot various fit parameters [default: False]

• freqcutoff = cutoff frequency (Hz) for optimal filtering [default: None]

• rotate = rotation angle (radians) for I/Q pulse record

• bsec = section of record to take for background (only used for ‘rotate’ or ‘absolute’) [default: 0.05]

• absolute = use sqrt(I^2+Q^2) signal [default: False]

• prlen = length of record to process (None = entire record) [default: None]

• avpulsepos = average position of pulse maximum; used for limited prlen. [default: None]

Returns:

• optf = optimal filter

• norm = optimal filter scaling factor