EPSearch.c File Reference

Prototypes
static double	min (double a, double b)
static double	max (double a, double b)
static LALExcessPowerFilterBank *	XLALCreateExcessPowerFilterBank (double filter_deltaF, double flow, double channel_bandwidth, int n_channels, const REAL8FrequencySeries *psd, const REAL8Sequence *two_point_spectral_correlation)
	From the power spectral density function, generate the comb of channel filters for the time-frequency plane — an excess power filter bank. More...
static void	XLALDestroyExcessPowerFilterBank (LALExcessPowerFilterBank *bank)
	Destroy and excess power filter bank. More...
int	XLALOverlappedSegmentsCommensurate (int target_length, int segment_length, int segment_shift)
	Round target_length down so that an integer number of intervals of length segment_length, each shifted by segment_shift with respect to the interval preceding it, fits into the result. More...
int	XLALEPGetTimingParameters (int window_length, int max_tile_length, double fractional_tile_shift, int *psd_length, int *psd_shift, int *window_shift, int *window_pad, int *tiling_length)
	Compute and return the timing parameters for an excess power analysis. More...
static REAL8TimeFrequencyPlane *	XLALCreateTFPlane (unsigned tseries_length, double tseries_deltaT, double flow, double bandwidth, double tiling_fractional_stride, double max_tile_bandwidth, double max_tile_duration, const REAL8FFTPlan *plan)
	Create and initialize a time-frequency plane object. More...
static void	XLALDestroyTFPlane (REAL8TimeFrequencyPlane *plane)
	Free a time-frequency plane object. More...
static COMPLEX16Sequence *	apply_filter (COMPLEX16Sequence *outputseq, const COMPLEX16FrequencySeries *inputseries, const COMPLEX16FrequencySeries *filterseries)
static int	XLALFreqSeriesToTFPlane (REAL8TimeFrequencyPlane *plane, const LALExcessPowerFilterBank *filter_bank, const COMPLEX16FrequencySeries *fseries, const REAL8FFTPlan *reverseplan)
static double	compute_unwhitened_mean_square (const LALExcessPowerFilterBank *filter_bank, unsigned channel, unsigned channels)
static SnglBurst *	XLALTFTileToBurstEvent (const REAL8TimeFrequencyPlane *plane, double tile_start, double tile_length, double f_centre, double bandwidth, double h_rss, double E, double d, double confidence)
static SnglBurst *	XLALComputeExcessPower (const REAL8TimeFrequencyPlane *plane, const LALExcessPowerFilterBank *filter_bank, SnglBurst *head, double confidence_threshold)
SnglBurst *	XLALEPSearch (LIGOLwXMLStream *diagnostics, const REAL8TimeSeries *tseries, REAL8Window *window, double flow, double bandwidth, double confidence_threshold, double fractional_stride, double maxTileBandwidth, double maxTileDuration)
	Generate a linked list of burst events from a time series. More...

Go to the source code of this file.

Data Structures
struct	ExcessPowerFilter
struct	LALExcessPowerFilterBank
struct	REAL8TimeFrequencyPlaneTiles
	A time-frequency plane. More...
struct	REAL8TimeFrequencyPlane

Function Documentation

min()

static double min ( double  a, double  b  )

static

Definition at line 58 of file EPSearch.c.

max()

static double max ( double  a, double  b  )

static

Definition at line 64 of file EPSearch.c.

XLALCreateExcessPowerFilterBank()

static LALExcessPowerFilterBank* XLALCreateExcessPowerFilterBank ( double  filter_deltaF, double  flow, double  channel_bandwidth, int  n_channels, const REAL8FrequencySeries *  psd, const REAL8Sequence *  two_point_spectral_correlation  )

static

From the power spectral density function, generate the comb of channel filters for the time-frequency plane — an excess power filter bank.

Definition at line 97 of file EPSearch.c.

XLALDestroyExcessPowerFilterBank()

static void XLALDestroyExcessPowerFilterBank ( LALExcessPowerFilterBank *  bank )

static

Destroy and excess power filter bank.

Definition at line 158 of file EPSearch.c.

XLALCreateTFPlane()

static REAL8TimeFrequencyPlane* XLALCreateTFPlane ( unsigned  tseries_length, double  tseries_deltaT, double  flow, double  bandwidth, double  tiling_fractional_stride, double  max_tile_bandwidth, double  max_tile_duration, const REAL8FFTPlan *  plan  )

static

Create and initialize a time-frequency plane object.

Parameters

tseries_length	length of time series from which TF plane will be computed
tseries_deltaT	sample rate of time series
flow	minimum frequency to search for
bandwidth	bandwidth of TF plane
tiling_fractional_stride	overlap of adjacent tiles
max_tile_bandwidth	largest tile's bandwidth
max_tile_duration	largest tile's duration
plan	forward plan whose length is tseries_length

Definition at line 396 of file EPSearch.c.

XLALDestroyTFPlane()

static void XLALDestroyTFPlane ( REAL8TimeFrequencyPlane *  plane )

static

Free a time-frequency plane object.

Definition at line 566 of file EPSearch.c.

apply_filter()

static COMPLEX16Sequence* apply_filter ( COMPLEX16Sequence *  outputseq, const COMPLEX16FrequencySeries *  inputseries, const COMPLEX16FrequencySeries *  filterseries  )

static

Definition at line 598 of file EPSearch.c.

XLALFreqSeriesToTFPlane()

static int XLALFreqSeriesToTFPlane ( REAL8TimeFrequencyPlane *  plane, const LALExcessPowerFilterBank *  filter_bank, const COMPLEX16FrequencySeries *  fseries, const REAL8FFTPlan *  reverseplan  )

static

Definition at line 635 of file EPSearch.c.

compute_unwhitened_mean_square()

static double compute_unwhitened_mean_square ( const LALExcessPowerFilterBank *  filter_bank, unsigned  channel, unsigned  channels  )

static

Definition at line 735 of file EPSearch.c.

XLALTFTileToBurstEvent()

static SnglBurst* XLALTFTileToBurstEvent ( const REAL8TimeFrequencyPlane *  plane, double  tile_start, double  tile_length, double  f_centre, double  bandwidth, double  h_rss, double  E, double  d, double  confidence  )

static

Definition at line 756 of file EPSearch.c.

XLALComputeExcessPower()

static SnglBurst* XLALComputeExcessPower ( const REAL8TimeFrequencyPlane *  plane, const LALExcessPowerFilterBank *  filter_bank, SnglBurst *  head, double  confidence_threshold  )

static

Definition at line 797 of file EPSearch.c.