LAL  7.5.0.1-8083555
LALGEOPsd.c
Go to the documentation of this file.
1 /*
2 * Copyright (C) 2007 Bernd Machenschalk, David Churches, Duncan Brown, Jolien Creighton
3 *
4 * This program is free software; you can redistribute it and/or modify
5 * it under the terms of the GNU General Public License as published by
6 * the Free Software Foundation; either version 2 of the License, or
7 * (at your option) any later version.
8 *
9 * This program is distributed in the hope that it will be useful,
10 * but WITHOUT ANY WARRANTY; without even the implied warranty of
11 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
12 * GNU General Public License for more details.
13 *
14 * You should have received a copy of the GNU General Public License
15 * along with with program; see the file COPYING. If not, write to the
16 * Free Software Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston,
17 * MA 02110-1301 USA
18 */
19 
20 #include <lal/LALNoiseModels.h>
21 
22 #ifdef __GNUC__
23 #define UNUSED __attribute__ ((unused))
24 #else
25 #define UNUSED
26 #endif
27 
28 /**
29  * \author Sathyaprakash, B. S.
30  * \ingroup LALNoiseModels_h
31  * \brief Module to calculate the expected noise power spectral density for the GEO600 detector.
32  *
33  * ### Description ###
34  *
35  * The module takes as an input a frequency \f$f\f$ in Hz, and it
36  * calculates the noise spectral density (per Hz) \f$S_{h}(f)\f$
37  * for that frequency. The noise PSD is based on data provided by
38  * J. Hough and G. Cagnoli (see T. Damour, B.R. Iyer and B.S. Sathyaprakash,
39  * Phys. Rev. D 63, 044023 (2001)) and is approximated by
40  * the following:
41  * \f{equation}{
42  * S_h(f) = 10^{-16} \left ( \frac{f}{f_0} \right)^{-30} +
43  * 34 \frac{f_0 }{ f } +
44  * \frac{20 \left [1 - (f/f_0)^2 + 0.5 (f/f_0)^4 \right ] }{ 1 + 0.5 (f/f_0)^2}
45  * \f}
46  * The returned value is scaled up by \f$s_0 = 10^{46}.\f$ In otherwords,
47  * the expected noise PSD is a factor \f$10^{46}\f$ lower.
48  *
49  */
50 
51 void
53 {
54 
55  REAL8 x, seismic, thermal, shot;
56 
57  x = f/150.;
58  seismic = pow(10.,-16.) * pow(x,-30.);
59  thermal = 34. / x;
60  shot = 20. * (1 - pow(x,2.) + 0.5 * pow(x,4.)) / (1. + 0.5 * pow(x,2.));
61  *psd = seismic + thermal + shot;
62 }
static double f(double theta, double y, double xi)
Definition: XLALMarcumQ.c:258
double REAL8
Double precision real floating-point number (8 bytes).
void LALGEOPsd(LALStatus UNUSED *status, REAL8 *psd, REAL8 f)
Module to calculate the expected noise power spectral density for the GEO600 detector.
Definition: LALGEOPsd.c:52
LAL status structure, see The LALStatus structure for more details.
Definition: LALDatatypes.h:947