Loading [MathJax]/extensions/TeX/AMSsymbols.js
LAL 7.7.0.1-da3b9d3
All Data Structures Namespaces Files Functions Variables Typedefs Enumerations Enumerator Macros Modules Pages
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 */
50void 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