LAL  7.5.0.1-8083555
LALEGOPsd.c
Go to the documentation of this file.
1 /*
2 * Copyright (C) 2007 Bernd Machenschalk, Thomas Cokelaer
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 Cokelaer T.
30  * \ingroup LALNoiseModels_h
31  * \brief Function to calculate the noise power spectral density for the EGO 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  * grqc/0607092
39  * \f{equation}{
40  * S_h(f) =
41  * s_0 \left\{ x^{p_1} + a_1x^{p_2} +a_2 \frac{1+b_1x +b_2x^2+b_3x^3+b_4x^4+b_5x^5+b_6x^6}{1+c_1x+c_2x^2+c_3x^3+c_4x^4} \right\}
42  * \f}
43  * where \f$S_0=1.61e-51\f$\\
44  * \f$p_1=-4.05, p_2=-0.69\f$\\
45  * \f$a_1=185.62, a_2=232.56\f$\\
46  * \f$b_1 = 31.18, b_2=-64.72, b_3=52.24, b_4=-42.16, b_5=10.17, b_6=11.53\f$\\
47  * and \f$c_1=13.58, c_2 = -36.46, c_3=18.56, c_4=27.43\f$
48  *
49  */
50 void LALEGOPsd (LALStatus UNUSED *status, REAL8 *psd, REAL8 f)
51 {
52  REAL8 s0, x, x2, x3, x4, x5, x6;
53  REAL8 a1, a2, p1, p2, c1, c2, c3, c4, b1, b2, b3, b4, b5, b6;
54  REAL8 num,den;
55 
56  x = f/200.;
57 
58  a1 = 185.62;
59  a2 = 232.56;
60  p1 = -4.05;
61  p2 = -0.69;
62  b1 = 31.18;
63  b2 = -64.72;
64  b3 = 52.24;
65  b4 = -42.16;
66  b5 = 10.17;
67  b6 = 11.53;
68  c1 = 13.58;
69  c2 = -36.46;
70  c3 = 18.56;
71  c4 = 27.43;
72 
73  x2 = x * x;
74  x3 = x2 * x;
75  x4 = x2 * x2;
76  x5 = x3 * x2;
77  x6 = x3 * x3;
78 
79  num = 1 + b1*x + b2*x2 + b3*x3 + b4*x4 + x5*b5 + x6*b6;
80  den = 1 + c1*x + c2*x2 + c3*x3 + c4*x4;
81 
82  /*new psds from fitted on the Design sensitivity curve from virgo web site*/
83  s0 = 1.62e-51;
84  *psd = s0*( pow(x,p1) +a1*pow(x,p2) +a2*num/den);
85 
86 }
static const UINT4 c1
Definition: LALCityHash.c:146
static const UINT4 c2
Definition: LALCityHash.c:147
static double f(double theta, double y, double xi)
Definition: XLALMarcumQ.c:258
double REAL8
Double precision real floating-point number (8 bytes).
void LALEGOPsd(LALStatus UNUSED *status, REAL8 *psd, REAL8 f)
Function to calculate the noise power spectral density for the EGO detector.
Definition: LALEGOPsd.c:50
LAL status structure, see The LALStatus structure for more details.
Definition: LALDatatypes.h:947