Function: splitup_freq()

0 float splitup_freq(float c0, float c90, float *chirp0, float *chirp90, float norm, float* twice_inv_noise, int n, int offset, int p, int* indices, float* stats, float* working, float* htilde)
This routine returns the value of the statistic $r^2=\sum_{i=1}^p (\Delta S_i)^2$. This is a less-efficient version, which internally constructs filters for each of the different frequency subintervals, and then filters the metric perturbation through those filters. It is useful to understand how the different frequency components behave in the time domain, after filtering.

The arguments are:

c0: Input. The coefficient of the 0-phase template.

c90: Input. The coefficient of the $90^\circ$-phase template. Note that $c_0^2 + c_{90}^2$ should be 1.

chirp0: Input. An array chirp0[0..n-1] containing the FFT of the 0-phase chirp.

chirp90: Input. An array chirp90[0..n-1] containing the FFT of the $90^\circ$-phase chirp.

norm: Input. The normalization of the 0-phase chirp.

twice_inv_noise: Input. The array twice_inv_noise[0..n/2] contains $2/S_h(f)$, as described previously. The array element twice_inv_noise[0] contains the DC value, and the array element twice_inv_noise[n/2] contains the value at the Nyquist frequency.

n: Input. Defines the lengths of the previous arrays.

offset: Input. The offset of the moment of maximum signal in the filter output.

p: Input. The number of frequency bands $p$ for the vetoing test.

indices: Output. An array indices[0..p-1] used for internal storage of the frequency subintervals (see splitup().

stats: Output. An array stats[0..p-1] containing the values of the $S_i$ for $i=1,\cdots,p$.

working: Output. An array working[0..n-1] used for internal storage.

htilde: Input. An array htilde[0..n-1] containing the positive frequency part of $\tilde h(f)$.

Author: Bruce Allen, ballen@dirac.phys.uwm.edu

Comments: None.