Loading [MathJax]/extensions/TeX/AMSsymbols.js
LAL 7.7.0.1-3a66518
All Data Structures Namespaces Files Functions Variables Typedefs Enumerations Enumerator Macros Modules Pages
LALVIRGOPsd.c
Go to the documentation of this file.
1/*
2* Copyright (C) 2007 Bernd Machenschalk, David Churches, Duncan Brown, Jolien Creighton, 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 Sathyaprakash, B. S., Cokelaer T.
30 * \ingroup LALNoiseModels_h
31 * \brief Module to calculate the noise power spectral density for the VIRGO 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-Y. Vinet and is approximated by
39 * the following:
40 * \f{equation}{
41 * S_h(f) =
42 * s_0 \left ( \frac {7.87f}{f_0} \right )^{-4.8} + \frac{6}{17} \frac{f_0}{f}
43 * + \left [1 + \left (\frac {f}{f_0} \right)^2 \right ],
44 * \f}
45 * where \f$s_0=10.2e-46\f$
46 *
47 */
48void LALVIRGOPsd (LALStatus UNUSED *status, REAL8 *psd, REAL8 f)
49{
50 REAL8 s0, x;
51
52 x = f/500.;
53/*
54 * s1 = 34.6;
55 s2 = 6.60;
56 s3 = 3.24;
57 *psd = pow(6.23*x,-5.) + 2.04/x + 1. + x*x;
58 */
59
60 /*new psds from fitted on the Design sensitivity curve from virgo web site*/
61 s0 = 10.2e-46;
62 *psd = s0*( pow(7.87*x,-4.8) + 6./17./x + 1. + x*x);
63
64}
static double f(double theta, double y, double xi)
Definition: XLALMarcumQ.c:258
double REAL8
Double precision real floating-point number (8 bytes).
void LALVIRGOPsd(LALStatus UNUSED *status, REAL8 *psd, REAL8 f)
Module to calculate the noise power spectral density for the VIRGO detector.
Definition: LALVIRGOPsd.c:48
LAL status structure, see The LALStatus structure for more details.
Definition: LALDatatypes.h:947