22#include<lal/LALStdlib.h>
23#include<lal/LALStatusMacros.h>
24#include<lal/LALInspiral.h>
25#include<lal/LALInspiralBank.h>
26#include<lal/LIGOMetadataTables.h>
47 REAL8 *chi, *kappa, dChi, dKappa;
53 LALINSPIRALBANKH_MSGENULL );
55 LALINSPIRALBANKH_MSGENULL );
57 LALINSPIRALBANKH_MSGENUMFCUT );
91 for( cnt = 0; cnt < *ntiles; cnt++ )
96 maxfFinal = 1.0 / (6.0 * sqrt(6.0)*
LAL_PI
111 maxfFinal = 1.0 / (3.0 * sqrt(3.0)*
LAL_PI
120 maxfFinal = ( 1. - 0.63*pow(1. - 3.4641016*coarseList[cnt].
params.eta
122 , 0.3) )/( 2.*
LAL_PI*(1. - 0.057191
123 *coarseList[cnt].
params.eta - 0.498
129 maxfFinal = 1.2* ( 1. - 0.63*pow(1. - 3.4641016*coarseList[cnt].
params.eta
131 , 0.3) )/( 2.*
LAL_PI*(1. - 0.057191
132 *coarseList[cnt].
params.eta - 0.498
141 minfFinal = 1.0 / (6.0 * sqrt(6.0)*
LAL_PI
156 minfFinal = 1.0 / (3.0 * sqrt(3.0)*
LAL_PI
165 minfFinal = ( 1. - 0.63*pow(1. - 3.4641016*coarseList[cnt].
params.eta
167 , 0.3) )/( 2.*
LAL_PI*(1. - 0.057191
168 *coarseList[cnt].
params.eta - 0.498
174 minfFinal = 1.2* ( 1. - 0.63*pow(1. - 3.4641016*coarseList[cnt].
params.eta
176 , 0.3) )/( 2.*
LAL_PI*(1. - 0.057191
177 *coarseList[cnt].
params.eta - 0.498
211 if (fFinal > input->
fUpper)
223 &(coarseList[cnt].
params) );
226 &(coarseList[cnt].
params), &moments);
266 fFinal = minfFinal +
i *
267 (maxfFinal - minfFinal)/(input->
numFreqCut - 1);
268 if (fFinal > input->
fUpper)
280 &(coarseList[cnt].
params) );
283 &(coarseList[cnt].
params), &moments);
297 bank = (*first)->
next;
313 if (!
i) chi[
i] = input->
chiMin + dChi / 2.0;
314 else chi[
i] = chi[0] +
i * dChi ;
318 if (!
i) kappa[
i] = input->
kappaMin + dKappa / 2.0;
319 else kappa[
i] = kappa[0] +
i * dKappa ;
333 for ( chicnt = 0; chicnt < input->
nPointsChi; chicnt++ )
335 for( kappacnt = 0; kappacnt < input->
nPointsKappa; kappacnt++ )
337 for( cnt = 0; cnt < *ntiles; cnt++ )
340 if ( coarseList[cnt].
params.mass1 < 6.0 ||
348 numTmplts = numTmplts + 1 ;
379 bank = (*first)->
next;
#define ABORT(statusptr, code, mesg)
#define TRY(func, statusptr)
#define ATTATCHSTATUSPTR(statusptr)
#define ASSERT(assertion, statusptr, code, mesg)
#define DETATCHSTATUSPTR(statusptr)
#define INITSTATUS(statusptr)
#define RETURN(statusptr)
#define LALINSPIRALBANKH_EFCUT
Inappropriate cutoff frequency [SchwarzISCO, BKLISCO, LightRing, ERD, FRD or LRD].
#define LALINSPIRALBANKH_ECHOICE
Invalid choice for an input parameter.
#define LALINSPIRALBANKH_ENULL
Null pointer.
void LALInspiralSpinBank(LALStatus *status, SnglInspiralTable **tiles, INT4 *ntiles, InspiralCoarseBankIn *coarseIn)
This function creates a bank of BCVSpin templates to search for precessing binaries.
void LALInspiralBankGeneration(LALStatus *status, InspiralCoarseBankIn *input, SnglInspiralTable **first, INT4 *ntiles)
void LALInspiralCreateCoarseBank(LALStatus *status, InspiralTemplateList **list, INT4 *nlist, InspiralCoarseBankIn bankIn)
void LALGetInspiralMoments(LALStatus *status, InspiralMomentsEtc *moments, REAL8FrequencySeries *psd, InspiralTemplate *params)
#define LALINSPIRALBANKH_ENUMFCUT
Number of fcut must be greater or equal to 1.
#define LALINSPIRALBANKH_EMEM
Memory allocation failure.
void LALInspiralComputeMetric(LALStatus *status, InspiralMetric *metric, InspiralTemplate *params, InspiralMomentsEtc *moments)
@ FreqCut_SchwarzISCO
the innermost stable circular orbit (ISCO) for a test particle orbiting a Schwarzschild black hole
@ FreqCut_FRD
the "Fundamental RingDown" frequency which is calculated from the Berti, Cardoso and Will (arXiv:gr-q...
@ FreqCut_LightRing
the unstable circular orbit for photons orbiting a Schwarzschild black hole
@ FreqCut_ERD
an effective ringdown frequency studied in Pan et al (arXiv:0704.1964) that was found to give good fi...
@ FreqCut_BKLISCO
a mass ratio dependent ISCO derived from estimates of the final spin of a merged black found in a pap...
@ FreqCut_LRD
the "Lorentzian RingDown" frequency = 1.2*FRD which captures part of the Lorentzian tail from the dec...
Input for choosing a template bank.
INT4 nPointsKappa
UNDOCUMENTED.
INT4 nPointsChi
UNDOCUMENTED.
REAL8 kappaMax
UNDOCUMENTED.
ComputeMoments computeMoments
ComputeMoments tells whether to re-compute the moments using an upper limit defined by flso; This is ...
INT4 numFreqCut
Number of different upper frequency cutoffs (spaced evenly between minFreqCut and maxFreqCut) to use ...
REAL8FrequencySeries shf
Frequency series containing the PSD.
REAL8 chiMin
UNDOCUMENTED.
INT4 spinBank
spinBank=0:use Owen+Hanna bank, spinBank=1:use extended bank by AEI/Cardiff/Osaka,...
FreqCut minFreqCut
smallest upper frequency cutoff to use
REAL8 fUpper
Upper frequency cutoff.
Approximant approximant
Approximant of the waveform.
REAL8 kappaMin
UNDOCUMENTED.
REAL8 chiMax
UNDOCUMENTED.
FreqCut maxFreqCut
largest upper frequency cutoff to use
REAL4 Gamma[10]
3d metric co-efficients in coordinates; Gamma[6] is a vector that stores the upper triangular part o...
Parameter structure that holds the moments of the PSD and other useful constants required in the comp...
REAL8 t5
2.5 post-Newtonian chirp time in seconds (output)
REAL8 psi3
BCV parameter .
REAL8 t2
first post-Newtonian chirp time in seconds (input/output)
REAL8 eta
symmetric mass ratio (input/output)
REAL8 totalMass
total mass of the binary in solar mass (input/output)
REAL8 mass1
Mass of the primary in solar mass (input/output)
REAL8 fCutoff
upper frequency cutoff in Hz to be used in generating the waveform; If the last stable orbit frequenc...
REAL8 t4
second post-Newtonian chirp time in seconds (output)
REAL8 t3
1.5 post-Newtonian chirp time in seconds (input/output)
REAL8 psi0
BCV parameter .
REAL8 fFinal
final frequency reached, in units of Hz (output)
REAL8 t0
Newtonain chirp time in seconds (input/output)
REAL8 chirpMass
chirp mass of the binary in solar mass (output)
REAL8 tC
total chirp time seconds (output)
REAL8 mass2
Mass of the secondary in solar mass (mass1 need not be larger than mass2 (input/output)
A grid of inspiral templates (ie a template list).
InspiralTemplate params
Value of the parameters at the lattice point.
InspiralMetric metric
metric at the lattice point
struct tagSnglInspiralTable * next