Loading [MathJax]/extensions/TeX/AMSsymbols.js
LALInspiral 5.0.3.1-8a6b96f
All Data Structures Namespaces Files Functions Variables Typedefs Enumerations Enumerator Macros Modules Pages
CoarseTest.c File Reference

Test code for the inspiral modules. More...

Detailed Description

Test code for the inspiral modules.

Author
Churches, D. K. and Sathyaprakash, B. S., Cokelaer, T

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.

Usage

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.);

  • Finally, fill the psd structure (see test code for an example). This involves allocating memory to vector shf.data and filling it with noise PSD as also choosing the starting frequency and the frequency resolution.
CoarseTest

Description

This test code gives an example of how one calls LALInspiralCreateCoarseBank and LALInspiralCreateFineBank modules.

Definition in file CoarseTest.c.

Go to the source code of this file.