lalapps » src » inspiral

tmpltbank.c File Reference

#include <config.h>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <getopt.h>
#include <unistd.h>
#include <sys/types.h>
#include <sys/stat.h>
#include <fcntl.h>
#include <regex.h>
#include <time.h>
#include <FrameL.h>
#include <lalapps.h>
#include <series.h>
#include <processtable.h>
#include <lalappsfrutils.h>
#include <lal/LALConfig.h>
#include <lal/LALStdio.h>
#include <lal/LALStdlib.h>
#include <lal/LALError.h>
#include <lal/LALDatatypes.h>
#include <lal/AVFactories.h>
#include <lal/LALConstants.h>
#include <lal/PrintFTSeries.h>
#include <lal/FrameStream.h>
#include <lal/FrameCalibration.h>
#include <lal/Window.h>
#include <lal/TimeFreqFFT.h>
#include <lal/IIRFilter.h>
#include <lal/ResampleTimeSeries.h>
#include <lal/BandPassTimeSeries.h>
#include <lal/LIGOMetadataTables.h>
#include <lal/LIGOMetadataUtils.h>
#include <lal/LIGOLwXML.h>
#include <lal/LIGOLwXMLRead.h>
#include <lal/Date.h>
#include <lal/Units.h>
#include <lal/LALInspiral.h>
#include <lal/LALInspiralBank.h>
#include "inspiral.h"

Include dependency graph for tmpltbank.c:

Go to the source code of this file.

Defines

#define CVS_ID_STRING   "$Id: tmpltbank.c,v 1.90 2008/07/16 08:35:40 whelan Exp $"

#define CVS_NAME_STRING   "$Name: $"

#define CVS_REVISION   "$Revision: 1.90 $"

#define CVS_SOURCE   "$Source: /usr/local/cvs/lscsoft/lalapps/src/inspiral/tmpltbank.c,v $"

#define CVS_DATE   "$Date: 2008/07/16 08:35:40 $"

#define PROGRAM_NAME   "tmpltbank"

#define ADD_PROCESS_PARAM(pptype, format, ppvalue)

#define USAGE(a)

Enumerations

enum { undefined, real_4, real_8 }

enum {
  specType_mean, specType_median, specType_gaussian, specType_LIGO,
  specType_AdvLIGO, specType_undefined
}

Functions

RCSID ("$Id: tmpltbank.c,v 1.90 2008/07/16 08:35:40 whelan Exp $")

int arg_parse_check (int argc, char *argv[], MetadataTable procparams)

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

MAIN function of SideBandMCMC code Compute the posterior pdfs of the orbital and nuisance parameters of a binary signal in Fstat form.

Variables

enum { ... } calData

enum { ... } specType

int vrbflg

LIGOTimeGPS gpsStartTime = { 0, 0 }

LIGOTimeGPS gpsEndTime = { 0, 0 }

INT4 padData = 0

CHAR * fqChanName = NULL

INT4 globFrameData = 0

CHAR * frInCacheName = NULL

CHAR * frInType = NULL

INT4 numPoints = -1

INT4 numSegments = -1

CHAR ifo [3]

CHAR * channelName = NULL

INT4 inputDataLength = 0

INT4 resampFiltType = -1

INT4 sampleRate = -1

INT4 highPass = -1

REAL4 highPassFreq = 0

INT4 highPassOrder = -1

REAL4 highPassAtten = -1

REAL4 fLow = -1

CHAR * calCacheName = NULL

INT4 globCalData = 0

INT4 pointCal = 0

REAL4 dynRangeExponent = 0

REAL4 strainHighPassFreq = -1

INT4 strainHighPassOrder = -1

REAL4 strainHighPassAtten = -1

REAL4 minMass = -1

REAL4 maxMass = -1

REAL4 minTotalMass = -1

REAL4 maxTotalMass = -1

REAL4 chirpMassCutoff = -1

REAL4 psi0Min = 0

REAL4 psi0Max = 0

REAL4 psi3Min = 0

REAL4 psi3Max = 0

REAL4 alpha = 0

REAL4 betaMin = 0

REAL4 betaMax = 0

INT4 maxFcutTmplts = -1

REAL4 minMatch = -1

REAL4 fUpper = -1

REAL4 chiMin = -1

REAL4 chiMax = -1

REAL4 kappaMin = -2

REAL4 kappaMax = -2

Order order

Approximant approximant

CoordinateSpace space

INT4 haveGridSpacing = 0

INT4 computeMoments = 1

FreqCut maxFreqCut

FreqCut minFreqCut

INT4 numFreqCut = 0

GridSpacing gridSpacing = SquareNotOriented

int polygonFit = 1

INT4 unitResponse = 0

INT4 computeCandle = 0

REAL4 candleSnr = -1

REAL4 candleMass1 = -1

REAL4 candleMass2 = -1

CHAR ** tdFileNames = NULL

INT4 numTDFiles = 0

CHAR * userTag = NULL

CHAR * ifoTag = NULL

int writeRawData = 0

int writeResponse = 0

int writeSpectrum = 0

int writeStrainSpec = 0

INT4 outCompress = 0

CHAR comment [LIGOMETA_COMMENT_MAX]

Define Documentation

#define CVS_ID_STRING "$Id: tmpltbank.c,v 1.90 2008/07/16 08:35:40 whelan Exp $"

Definition at line 75 of file tmpltbank.c.

#define CVS_NAME_STRING "$Name: $"

Definition at line 76 of file tmpltbank.c.

#define CVS_REVISION "$Revision: 1.90 $"

Definition at line 77 of file tmpltbank.c.

#define CVS_SOURCE "$Source: /usr/local/cvs/lscsoft/lalapps/src/inspiral/tmpltbank.c,v $"

Definition at line 78 of file tmpltbank.c.

#define CVS_DATE "$Date: 2008/07/16 08:35:40 $"

Definition at line 79 of file tmpltbank.c.

#define PROGRAM_NAME "tmpltbank"

Definition at line 80 of file tmpltbank.c.

#define ADD_PROCESS_PARAM	(	pptype,
		format,
		ppvalue	)

Value:

this_proc_param = this_proc_param->next = (ProcessParamsTable *)   calloc( 1, sizeof(ProcessParamsTable) );  LALSnprintf( this_proc_param->program, LIGOMETA_PROGRAM_MAX, "%s",    PROGRAM_NAME );   LALSnprintf( this_proc_param->param, LIGOMETA_PARAM_MAX, "--%s",      long_options[option_index].name );     LALSnprintf( this_proc_param->type, LIGOMETA_TYPE_MAX, "%s", pptype );     LALSnprintf( this_proc_param->value, LIGOMETA_VALUE_MAX, format, ppvalue );

Definition at line 1267 of file tmpltbank.c.

#define USAGE ( a )

Definition at line 1277 of file tmpltbank.c.

Enumeration Type Documentation

anonymous enum

Enumerator:

undefined
real_4
real_8

Definition at line 85 of file tmpltbank.c.

anonymous enum

Enumerator:

specType_mean
specType_median
specType_gaussian
specType_LIGO
specType_AdvLIGO
specType_undefined

Definition at line 93 of file tmpltbank.c.

Function Documentation

RCSID	(	"$Id: tmpltbank.	c,
		v 1.90 2008/07/16 08:35:40 whelan Exp $"
	)

int arg_parse_check	(	int	argc,
		char *	argv[],
		MetadataTable	procparams
	)

Definition at line 233 of file FittingFactors.c.

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

MAIN function of SideBandMCMC code Compute the posterior pdfs of the orbital and nuisance parameters of a binary signal in Fstat form.

variable declarations

the status must be initially blank; this can be done by making it static:

static LALStatus status;

other variables:

ExampleOutput output; ExampleInput input; ExampleParams params;

parse arguments, if desired

for example:

program = *argv; while ( --argc > 0 ) { ++argv; if ( ! strcmp( *argv, "-d" ) ) { --argc; ++argv; lalDebugLevel = atoi( *argv ); continue; } if ( ! strcmp( *argv, "-v" ) ) { verbose = 1; continue; } if ( ! strcmp( *argv, "-h" ) ) { fprintf( stderr, usage, program ); return 0; } fprintf( stderr, "no such option %s\n", *argv ); fprintf( stderr, usage, program ); return FAIL_ARGS; }

/** test response to invalid data

these tests must be wrapped so they are not done when debugging is disabled:

ifndef LAL_NDEBUG if ( ! lalNoDebug ) { LALExample( &status, NULL, &input, &params ); if ( status.statusCode != EXAMPLEH_ENULL || strcmp( status.statusDescription, EXAMPLEH_MSGENULL ) ) { fprintf( stderr, "incorrect error code %d and message %s\n", status.statusCode, status.statusDescription ); fprintf( stderr, "expecting error code %d and message %s\n", EXAMPLEH_ENULL, EXAMPLEH_MSGENULL ); return FAIL_CODE; }

fputs( "PASS: Test response to invalid data\n", stderr ); } endif

test response to valid data

for example:

LALExample( &status, &output, &input, &param ); if ( status.statusCode ) { fprintf( stderr, "received error code %d and message %s\n", status.statusCode, status.statusDescription ); return FAIL_ESUB; }

(perform checks on the contents of output too!)

during default operation, output to the screen should be minimal:

for ( i = 0; i < 1048576; ++i ) printf( "%d\n", i ); !!! BAD !!!

but it is OK to indicate that tests have passed:

fputs( "PASS: Test response to valid data\n", stderr );

check for memory leaks and return success:

LALCheckMemoryLeaks(); fputs( "PASS: All tests\n" );

Definition at line 198 of file tmpltbank.c.

Variable Documentation

enum { ... } calData

enum { ... } specType

int vrbflg

Definition at line 45 of file lalapps.c.

LIGOTimeGPS gpsStartTime = { 0, 0 }

Definition at line 114 of file tmpltbank.c.

LIGOTimeGPS gpsEndTime = { 0, 0 }

Definition at line 115 of file tmpltbank.c.

INT4 padData = 0

Definition at line 116 of file tmpltbank.c.

CHAR* fqChanName = NULL

Definition at line 117 of file tmpltbank.c.

INT4 globFrameData = 0

Definition at line 118 of file tmpltbank.c.

CHAR* frInCacheName = NULL

Definition at line 119 of file tmpltbank.c.

CHAR* frInType = NULL

Definition at line 120 of file tmpltbank.c.

INT4 numPoints = -1

Definition at line 121 of file tmpltbank.c.

INT4 numSegments = -1

Definition at line 122 of file tmpltbank.c.

CHAR ifo[3]

Definition at line 123 of file tmpltbank.c.

CHAR* channelName = NULL

Definition at line 124 of file tmpltbank.c.

INT4 inputDataLength = 0

Definition at line 125 of file tmpltbank.c.

INT4 resampFiltType = -1

Definition at line 126 of file tmpltbank.c.

INT4 sampleRate = -1

Definition at line 127 of file tmpltbank.c.

INT4 highPass = -1

Definition at line 128 of file tmpltbank.c.

REAL4 highPassFreq = 0

Definition at line 129 of file tmpltbank.c.

INT4 highPassOrder = -1

Definition at line 130 of file tmpltbank.c.

REAL4 highPassAtten = -1

Definition at line 131 of file tmpltbank.c.

REAL4 fLow = -1

Definition at line 132 of file tmpltbank.c.

CHAR* calCacheName = NULL

Definition at line 133 of file tmpltbank.c.

INT4 globCalData = 0

Definition at line 134 of file tmpltbank.c.

INT4 pointCal = 0

Definition at line 135 of file tmpltbank.c.

REAL4 dynRangeExponent = 0

Definition at line 136 of file tmpltbank.c.

REAL4 strainHighPassFreq = -1

Definition at line 137 of file tmpltbank.c.

INT4 strainHighPassOrder = -1

Definition at line 138 of file tmpltbank.c.

REAL4 strainHighPassAtten = -1

Definition at line 139 of file tmpltbank.c.

REAL4 minMass = -1

Definition at line 142 of file tmpltbank.c.

REAL4 maxMass = -1

Definition at line 143 of file tmpltbank.c.

REAL4 minTotalMass = -1

Definition at line 144 of file tmpltbank.c.

REAL4 maxTotalMass = -1

Definition at line 145 of file tmpltbank.c.

REAL4 chirpMassCutoff = -1

Definition at line 146 of file tmpltbank.c.

REAL4 psi0Min = 0

Definition at line 147 of file tmpltbank.c.

REAL4 psi0Max = 0

Definition at line 148 of file tmpltbank.c.

REAL4 psi3Min = 0

Definition at line 149 of file tmpltbank.c.

REAL4 psi3Max = 0

Definition at line 150 of file tmpltbank.c.

REAL4 alpha = 0

Definition at line 151 of file tmpltbank.c.

REAL4 betaMin = 0

Definition at line 152 of file tmpltbank.c.

REAL4 betaMax = 0

Definition at line 153 of file tmpltbank.c.

INT4 maxFcutTmplts = -1

Definition at line 154 of file tmpltbank.c.

REAL4 minMatch = -1

Definition at line 155 of file tmpltbank.c.

REAL4 fUpper = -1

Definition at line 156 of file tmpltbank.c.

REAL4 chiMin = -1

Definition at line 157 of file tmpltbank.c.

REAL4 chiMax = -1

Definition at line 158 of file tmpltbank.c.

REAL4 kappaMin = -2

Definition at line 159 of file tmpltbank.c.

REAL4 kappaMax = -2

Definition at line 160 of file tmpltbank.c.

Order order

Definition at line 161 of file tmpltbank.c.

Approximant approximant

Definition at line 162 of file tmpltbank.c.

CoordinateSpace space

Definition at line 163 of file tmpltbank.c.

INT4 haveGridSpacing = 0

Definition at line 164 of file tmpltbank.c.

INT4 computeMoments = 1

Definition at line 165 of file tmpltbank.c.

FreqCut maxFreqCut

Definition at line 166 of file tmpltbank.c.

FreqCut minFreqCut

Definition at line 167 of file tmpltbank.c.

INT4 numFreqCut = 0

Definition at line 168 of file tmpltbank.c.

GridSpacing gridSpacing = SquareNotOriented

Definition at line 170 of file tmpltbank.c.

int polygonFit = 1

Definition at line 171 of file tmpltbank.c.

INT4 unitResponse = 0

Definition at line 174 of file tmpltbank.c.

INT4 computeCandle = 0

Definition at line 177 of file tmpltbank.c.

REAL4 candleSnr = -1

Definition at line 178 of file tmpltbank.c.

REAL4 candleMass1 = -1

Definition at line 179 of file tmpltbank.c.

REAL4 candleMass2 = -1

Definition at line 180 of file tmpltbank.c.

CHAR** tdFileNames = NULL

Definition at line 183 of file tmpltbank.c.

INT4 numTDFiles = 0

Definition at line 184 of file tmpltbank.c.

CHAR* userTag = NULL

Definition at line 187 of file tmpltbank.c.

CHAR* ifoTag = NULL

Definition at line 188 of file tmpltbank.c.

int writeRawData = 0

Definition at line 189 of file tmpltbank.c.

int writeResponse = 0

Definition at line 190 of file tmpltbank.c.

int writeSpectrum = 0

Definition at line 191 of file tmpltbank.c.

int writeStrainSpec = 0

Definition at line 192 of file tmpltbank.c.

INT4 outCompress = 0

Definition at line 193 of file tmpltbank.c.

CHAR comment[LIGOMETA_COMMENT_MAX]

Definition at line 196 of file tmpltbank.c.

Generated on Fri Aug 29 02:50:25 2008 for LAL by

1.5.2


Defines
#define	CVS_ID_STRING "$Id: tmpltbank.c,v 1.90 2008/07/16 08:35:40 whelan Exp $"
#define	CVS_NAME_STRING "$Name: $"
#define	CVS_REVISION "$Revision: 1.90 $"
#define	CVS_SOURCE "$Source: /usr/local/cvs/lscsoft/lalapps/src/inspiral/tmpltbank.c,v $"
#define	CVS_DATE "$Date: 2008/07/16 08:35:40 $"
#define	PROGRAM_NAME "tmpltbank"
#define	ADD_PROCESS_PARAM(pptype, format, ppvalue)
#define	USAGE(a)
Enumerations
enum	{ undefined, real_4, real_8 }
enum	{ specType_mean, specType_median, specType_gaussian, specType_LIGO, specType_AdvLIGO, specType_undefined }
Functions
	RCSID ("$Id: tmpltbank.c,v 1.90 2008/07/16 08:35:40 whelan Exp $")
int	arg_parse_check (int argc, char *argv[], MetadataTable procparams)
int	main (int argc, char *argv[])
	MAIN function of SideBandMCMC code Compute the posterior pdfs of the orbital and nuisance parameters of a binary signal in Fstat form.
Variables
enum { ... }	calData
enum { ... }	specType
int	vrbflg
LIGOTimeGPS	gpsStartTime = { 0, 0 }
LIGOTimeGPS	gpsEndTime = { 0, 0 }
INT4	padData = 0
CHAR *	fqChanName = NULL
INT4	globFrameData = 0
CHAR *	frInCacheName = NULL
CHAR *	frInType = NULL
INT4	numPoints = -1
INT4	numSegments = -1
CHAR	ifo [3]
CHAR *	channelName = NULL
INT4	inputDataLength = 0
INT4	resampFiltType = -1
INT4	sampleRate = -1
INT4	highPass = -1
REAL4	highPassFreq = 0
INT4	highPassOrder = -1
REAL4	highPassAtten = -1
REAL4	fLow = -1
CHAR *	calCacheName = NULL
INT4	globCalData = 0
INT4	pointCal = 0
REAL4	dynRangeExponent = 0
REAL4	strainHighPassFreq = -1
INT4	strainHighPassOrder = -1
REAL4	strainHighPassAtten = -1
REAL4	minMass = -1
REAL4	maxMass = -1
REAL4	minTotalMass = -1
REAL4	maxTotalMass = -1
REAL4	chirpMassCutoff = -1
REAL4	psi0Min = 0
REAL4	psi0Max = 0
REAL4	psi3Min = 0
REAL4	psi3Max = 0
REAL4	alpha = 0
REAL4	betaMin = 0
REAL4	betaMax = 0
INT4	maxFcutTmplts = -1
REAL4	minMatch = -1
REAL4	fUpper = -1
REAL4	chiMin = -1
REAL4	chiMax = -1
REAL4	kappaMin = -2
REAL4	kappaMax = -2
Order	order
Approximant	approximant
CoordinateSpace	space
INT4	haveGridSpacing = 0
INT4	computeMoments = 1
FreqCut	maxFreqCut
FreqCut	minFreqCut
INT4	numFreqCut = 0
GridSpacing	gridSpacing = SquareNotOriented
int	polygonFit = 1
INT4	unitResponse = 0
INT4	computeCandle = 0
REAL4	candleSnr = -1
REAL4	candleMass1 = -1
REAL4	candleMass2 = -1
CHAR **	tdFileNames = NULL
INT4	numTDFiles = 0
CHAR *	userTag = NULL
CHAR *	ifoTag = NULL
int	writeRawData = 0
int	writeResponse = 0
int	writeSpectrum = 0
int	writeStrainSpec = 0
INT4	outCompress = 0
CHAR	comment [LIGOMETA_COMMENT_MAX]