code to perform a coarse/fine heterodyne on LIGO or GEO data given a set of pulsar parameters. The heterodyne is a 2 stage process with the code run for a coarse heterodyne (not taking into account the SSB and BSB time delays, but using the other frequency params) and then rerun for the fine heterodyne (taking into account SSB and BSB time delays). The code can also be used to update heterodyned data using new parameters (old and new parameter files are required) i.e. it will take the difference of the original heterodyne phase and the new phase and reheterodyne with this phase difference. More...

Prototypes
static const char *	LogTimeToString (double t)

int	main (int argc, char *argv[])

void	get_input_args (InputParams inputParams, int argc, char argv[])

void	heterodyne_data (COMPLEX16TimeSeries data, REAL8Vector times, HeterodyneParams hetParams, REAL8 freqfactor, FilterResponse *filtresp)

void	get_frame_times (CHAR framefile, REAL8 gpstime, INT4 *duration)

REAL8TimeSeries *	get_frame_data (LALCache framecache, CHAR channel, REAL8 ttime, REAL8 length, INT4 duration, REAL8 samplerate, REAL8 scalefac, REAL8 highpass)

void	set_filters (Filters *iirFilters, REAL8 filterKnee, REAL8 samplerate)

void	filter_data (COMPLEX16TimeSeries data, Filters iirFilters)

COMPLEX16TimeSeries *	resample_data (COMPLEX16TimeSeries data, REAL8Vector times, INT4Vector starts, INT4Vector stops, REAL8 sampleRate, REAL8 resampleRate, INT4 hetflag)

INT4	get_segment_list (INT4Vector starts, INT4Vector stops, CHAR *seglistfile, INT4 heterodyneflag)

LALCache *	set_frame_files (INT4 starts, INT4 stops, LALCache cache, INT4 position, INT4 maxchunklength)

void	calibrate (COMPLEX16TimeSeries series, REAL8Vector datatimes, CalibrationFiles calfiles, REAL8 frequency, CHAR *channel)

void	get_calibration_values (REAL8 magnitude, REAL8 phase, CHAR *calibfilename, REAL8 frequency)

INT4	remove_outliers (COMPLEX16TimeSeries data, REAL8Vector times, REAL8 stddevthresh)

FilterResponse *	create_filter_response (REAL8 filterKnee)

void	destroy_filter_response (FilterResponse *filtresp)

Detailed Description

code to perform a coarse/fine heterodyne on LIGO or GEO data given a set of pulsar parameters. The heterodyne is a 2 stage process with the code run for a coarse heterodyne (not taking into account the SSB and BSB time delays, but using the other frequency params) and then rerun for the fine heterodyne (taking into account SSB and BSB time delays). The code can also be used to update heterodyned data using new parameters (old and new parameter files are required) i.e. it will take the difference of the original heterodyne phase and the new phase and reheterodyne with this phase difference.

Author: M. D. Pitkin

Definition in file heterodyne_pulsar.c.

Go to the source code of this file.

Macros
#define	ROUND(a) (floor(a+0.5))

#define	MAXALLOC 1048576

#define	TRACKMEMUSE 0

Variables
INT4	verbose = 0

Macro Definition Documentation

◆ ROUND

#define ROUND ( a ) (floor(a+0.5))

Definition at line 47 of file heterodyne_pulsar.c.

◆ MAXALLOC

#define MAXALLOC 1048576

Definition at line 50 of file heterodyne_pulsar.c.

◆ TRACKMEMUSE

#define TRACKMEMUSE 0

Definition at line 53 of file heterodyne_pulsar.c.

Function Documentation

◆ LogTimeToString()

const char * LogTimeToString ( double t )

static

Definition at line 2497 of file heterodyne_pulsar.c.

◆ main()

int main	(	int	argc,
		char *	argv[]
	)

Definition at line 74 of file heterodyne_pulsar.c.

◆ get_input_args()

void get_input_args	(	InputParams *	inputParams,
		int	argc,
		char *	argv[]
	)

Definition at line 812 of file heterodyne_pulsar.c.

◆ heterodyne_data()

void heterodyne_data	(	COMPLEX16TimeSeries *	data,
		REAL8Vector *	times,
		HeterodyneParams	hetParams,
		REAL8	freqfactor,
		FilterResponse *	filtresp
	)

Definition at line 1127 of file heterodyne_pulsar.c.

◆ get_frame_times()

void get_frame_times	(	CHAR *	framefile,
		REAL8 *	gpstime,
		INT4 *	duration
	)

Definition at line 1626 of file heterodyne_pulsar.c.

◆ get_frame_data()

REAL8TimeSeries * get_frame_data	(	LALCache *	framecache,
		CHAR *	channel,
		REAL8	ttime,
		REAL8	length,
		INT4	duration,
		REAL8	samplerate,
		REAL8	scalefac,
		REAL8	highpass
	)

Definition at line 1650 of file heterodyne_pulsar.c.

◆ set_filters()

void set_filters	(	Filters *	iirFilters,
		REAL8	filterKnee,
		REAL8	samplerate
	)

Definition at line 1741 of file heterodyne_pulsar.c.

◆ filter_data()

void filter_data	(	COMPLEX16TimeSeries *	data,
		Filters *	iirFilters
	)

Definition at line 1776 of file heterodyne_pulsar.c.

◆ resample_data()

COMPLEX16TimeSeries * resample_data	(	COMPLEX16TimeSeries *	data,
		REAL8Vector *	times,
		INT4Vector *	starts,
		INT4Vector *	stops,
		REAL8	sampleRate,
		REAL8	resampleRate,
		INT4	hetflag
	)

Definition at line 1797 of file heterodyne_pulsar.c.

◆ get_segment_list()

INT4 get_segment_list	(	INT4Vector *	starts,
		INT4Vector *	stops,
		CHAR *	seglistfile,
		INT4	heterodyneflag
	)

Definition at line 1973 of file heterodyne_pulsar.c.

◆ set_frame_files()

LALCache * set_frame_files	(	INT4 *	starts,
		INT4 *	stops,
		LALCache *	cache,
		INT4 *	position,
		INT4	maxchunklength
	)

Definition at line 2058 of file heterodyne_pulsar.c.

◆ calibrate()

void calibrate	(	COMPLEX16TimeSeries *	series,
		REAL8Vector *	datatimes,
		CalibrationFiles	calfiles,
		REAL8	frequency,
		CHAR *	channel
	)

Definition at line 2101 of file heterodyne_pulsar.c.

◆ get_calibration_values()

void get_calibration_values	(	REAL8 *	magnitude,
		REAL8 *	phase,
		CHAR *	calibfilename,
		REAL8	frequency
	)

Definition at line 2269 of file heterodyne_pulsar.c.

◆ remove_outliers()

INT4 remove_outliers	(	COMPLEX16TimeSeries *	data,
		REAL8Vector *	times,
		REAL8	stddevthresh
	)

Definition at line 2316 of file heterodyne_pulsar.c.

◆ create_filter_response()

FilterResponse * create_filter_response ( REAL8 filterKnee )

Definition at line 2368 of file heterodyne_pulsar.c.

◆ destroy_filter_response()

void destroy_filter_response ( FilterResponse * filtresp )

Definition at line 2489 of file heterodyne_pulsar.c.

Variable Documentation

◆ verbose

INT4 verbose = 0

Definition at line 72 of file heterodyne_pulsar.c.

Prototypes