Test code for the inspiral modules. More...
Test code for the inspiral modules.
A template bank can be created either using a full range for component masses of the binary \(m_1\) and \(m_2\), that is (mMin,mMax),
OR minimum value of the component masses mMin
and maximum value of the total mass MMax.
In the first case chirptimes satisfying the constraint mMin\
\(\le m_1, m_2 \le\)\ mMax
are accepted as valid systems. In the second case chirptimes satisfying the constraint mMin\
\(\le m_1, m_2,\) and MMax
\(\le m=m_1+m_2,\) are treated as valid. Users are expected to provide both mMax
and MMax
.
For LALLIGOIPsd
the choice mMin
\(=1M_\odot\) mMax
\(=20 M_\odot\) gives 2292 templates, while the same mMin
but choosing MMax
\(=40 M_\odot\) instead gives 2512 templates – about 10% incrase. However, the extra templates are ALL short-lived templates and therefore potentially trigger a large number of false alarms, as noted by Brown in E7 analysis.
This test code creates a template bank and stores it into CoarseTest.out . Then, it creates a finer template bank around a sub-set of the original template bank and stores it in the same output file.
Input the following values of the InspiralCoarseBankIn structure to create a template bank:
Minimum component mass in solar masses. mMin = 1.0;
Maximum component mass in solar masses. mMax = 20.0;
Maximum total mass. This should be specified independently of whether or not mMax is specified.
It is used in setting up a rectangular area in the space of chirptimes where the templates will be laid. MMax = 40.0;
The type of search space. massRange = MinComponentMassMaxTotalMass;
OR massRange = MinMaxComponentMasses;
Coarse bank minimal match coarseIn->mmCoarse = 0.80;
Fine bank minimal match mmFine = 0.97;
Lower frequency cutoff fLower = 40.;
Upper frequency cutoff fUpper = 1024L;
Whether or not lso should be used as an upper frequency cutoff (Currently not used; so please specify fUpper
. coarseIn->iflso = 0;
Sampling rate tSampling = 4000L;
Space in which templates should be created: whether \((\tau_0,\tau_2)\) or \((\tau_0, \tau_3).\) coarseIn->space = Tau0Tau2;
OR coarseIn->space = Tau0Tau3;
OR
Order and type of the approximant to be used in template generation. These members will NOT be used in creating a template bank but in filling up the InspiralTemplate
structure created by the bank.
order = twoPN;
coarseIn->approximant = TaylorT1;
minimum value of eta etamin = mMin * ( MMax - mMin)/pow(MMax,2.);
shf.data
and filling it with noise PSD as also choosing the starting frequency and the frequency resolution. This test code gives an example of how one calls LALInspiralCreateCoarseBank
and LALInspiralCreateFineBank
modules.
Definition in file CoarseTest.c.