#include <stdio.h>#include <stdlib.h>#include <config.h>#include <math.h>#include <ctype.h>#include <assert.h>#include <string.h>#include <getopt.h>#include <time.h>#include <lalapps.h>#include <processtable.h>#include <lal/LALStdio.h>#include <lal/LALStdlib.h>#include <lal/LALConstants.h>#include <lal/LIGOMetadataTables.h>#include <lal/LIGOMetadataUtils.h>#include <lal/LIGOLwXML.h>#include <lal/Units.h>#include <lal/Date.h>#include <lal/Inject.h>#include <lal/InspiralInjectionParams.h>#include <lal/GenerateInspiral.h>#include <lal/GenerateInspRing.h>#include <lal/FindChirp.h>#include <lal/GenerateRing.h>#include <lal/Ring.h>#include <lal/LALNoiseModels.h>#include <lal/RealFFT.h>#include <lal/FrequencySeries.h>#include <lal/TimeSeries.h>#include <lal/TimeFreqFFT.h>#include <lal/VectorOps.h>Include dependency graph for blindinj.c:
Go to the source code of this file.
Data Structures | |
| struct | actuationparameters |
Defines | |
| #define | CVS_ID_STRING "$Id: blindinj.c,v 1.17 2008/07/03 23:50:25 dfazi Exp $" |
| #define | CVS_NAME_STRING "$Name: $" |
| #define | CVS_REVISION "$Revision: 1.17 $" |
| #define | CVS_SOURCE "$Source: /usr/local/cvs/lscsoft/lalapps/src/inspiral/blindinj.c,v $" |
| #define | CVS_DATE "$Date: 2008/07/03 23:50:25 $" |
| #define | PROGRAM_NAME "blindinj" |
| #define | USAGE |
Typedefs | |
| typedef actuationparameters | ActuationParameters |
Enumerations | |
| enum | ResponseFunction { noResponse, unityResponse, LIGOdesign, actuationX, actuationY } |
Functions | |
| RCSID ("$Id: blindinj.c,v 1.17 2008/07/03 23:50:25 dfazi Exp $") | |
| 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 | |
| int | vrbflg |
| ActuationParameters | actuationParams [LAL_NUM_IFO] |
| REAL8 | dynRange = 1.0/3.0e-23 |
| INT4 | duration = 100.0 |
| INT4 | sampleRate = 16384 |
| REAL4 | fLower = 30 |
| REAL4 | mergerLength = 10 |
| REAL8 | longestSignal = 95.0 |
| REAL8 | timeWindow = 4 |
| UINT4 | numInjections = 0 |
| INT4 | randSeed = 1 |
| REAL4 | minNSMass = 1.0 |
| REAL4 | maxNSMass = 2.0 |
| REAL4 | minBHMass = 2.0 |
| REAL4 | maxBHMass = 30.0 |
| REAL4 | minTotalMass = 0.0 |
| REAL4 | maxTotalMass = 35.0 |
| REAL4 | minNSSpin = 0.0 |
| REAL4 | maxNSSpin = 0.2 |
| REAL4 | minBHSpin = 0.0 |
| REAL4 | maxBHSpin = 1.0 |
| REAL4 | BNSfrac = 0.35 |
| REAL4 | BBHfrac = 0.35 |
| REAL4 | bnsSnrMean = 9.0 |
| REAL4 | bnsSnrStd = 0.5 |
| REAL4 | snrMean = 12.0 |
| REAL4 | snrStd = 1.0 |
| REAL4Vector * | normalDev |
| #define CVS_ID_STRING "$Id: blindinj.c,v 1.17 2008/07/03 23:50:25 dfazi Exp $" |
Definition at line 66 of file blindinj.c.
| #define CVS_NAME_STRING "$Name: $" |
Definition at line 67 of file blindinj.c.
| #define CVS_REVISION "$Revision: 1.17 $" |
Definition at line 68 of file blindinj.c.
| #define CVS_SOURCE "$Source: /usr/local/cvs/lscsoft/lalapps/src/inspiral/blindinj.c,v $" |
Definition at line 69 of file blindinj.c.
| #define CVS_DATE "$Date: 2008/07/03 23:50:25 $" |
Definition at line 70 of file blindinj.c.
| #define PROGRAM_NAME "blindinj" |
Definition at line 71 of file blindinj.c.
| #define USAGE |
Value:
"lalapps_blindinj [options]\n""\nDefaults are shown in brackets\n\n" " --help display this message\n"" --version print version information and exit\n"" --verbose be verbose\n"" --gps-start-time TIME start time of injection\n"" --injection-type TYPE type of injection, must be one of \n"" (strain, etmx, etmy)\n"" --seed SEED seed random number generator with SEED (1)\n"" --debug-level LEVEL set the LAL debug level to LEVEL\n""\n"
Definition at line 73 of file blindinj.c.
| typedef struct actuationparameters ActuationParameters |
| enum ResponseFunction |
Definition at line 88 of file blindinj.c.
| RCSID | ( | "$Id: blindinj. | c, | |
| v 1.17 2008/07/03 23:50:25 dfazi Exp $" | ||||
| ) |
| 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, ¶ms ); 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, ¶m ); 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 426 of file blindinj.c.
| ActuationParameters actuationParams[LAL_NUM_IFO] |
Definition at line 126 of file blindinj.c.
Definition at line 127 of file blindinj.c.
Definition at line 128 of file blindinj.c.
| INT4 sampleRate = 16384 |
Definition at line 129 of file blindinj.c.
Definition at line 132 of file blindinj.c.
| REAL4 mergerLength = 10 |
Definition at line 133 of file blindinj.c.
| REAL8 longestSignal = 95.0 |
Definition at line 134 of file blindinj.c.
| REAL8 timeWindow = 4 |
Definition at line 135 of file blindinj.c.
| UINT4 numInjections = 0 |
Definition at line 137 of file blindinj.c.
Definition at line 138 of file blindinj.c.
Definition at line 139 of file blindinj.c.
Definition at line 140 of file blindinj.c.
Definition at line 141 of file blindinj.c.
Definition at line 142 of file blindinj.c.
| REAL4 minTotalMass = 0.0 |
Definition at line 143 of file blindinj.c.
| REAL4 maxTotalMass = 35.0 |
Definition at line 144 of file blindinj.c.
Definition at line 146 of file blindinj.c.
Definition at line 147 of file blindinj.c.
Definition at line 148 of file blindinj.c.
Definition at line 149 of file blindinj.c.
Definition at line 151 of file blindinj.c.
Definition at line 152 of file blindinj.c.
| REAL4 bnsSnrMean = 9.0 |
Definition at line 155 of file blindinj.c.
Definition at line 156 of file blindinj.c.
Definition at line 157 of file blindinj.c.
Definition at line 158 of file blindinj.c.
Definition at line 160 of file blindinj.c.
1.5.2