Loading [MathJax]/extensions/TeX/AMSsymbols.js
LALSimulation 6.2.0.1-ea7c608
•All Data Structures Namespaces Files Functions Variables Typedefs Enumerations Enumerator Macros Modules Pages
LALSimIMRPhenomTHM_fits.h
Go to the documentation of this file.
1#ifndef _LALSIM_IMR_PHENOMT_FITS_H
2#define _LALSIM_IMR_PHENOMT_FITS_H
3
4/*
5 * Copyright (C) 2020 Hector Estelles
6 *
7 * This program is free software; you can redistribute it and/or modify
8 * it under the terms of the GNU General Public License as published by
9 * the Free Software Foundation; either version 2 of the License, or
10 * (at your option) any later version.
11 *
12 * This program is distributed in the hope that it will be useful,
13 * but WITHOUT ANY WARRANTY; without even the implied warranty of
14 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
15 * GNU General Public License for more details.
16 *
17 * You should have received a copy of the GNU General Public License
18 * along with with program; see the file COPYING. If not, write to the
19 * Free Software Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston,
20 * MA 02110-1301 USA
21 */
22
23
24/**
25 * \author Hector Estelles
26 *
27 */
28
29#ifdef __cplusplus
30extern "C" {
31#endif
32
33#ifdef __GNUC__
34#define UNUSED __attribute__((unused))
35#else
36#define UNUSED
37#endif
38
39#include "LALSimIMRPhenomXUtilities.h"
40
41
42/********************************* IMRPhenomT 22 Frequency Fits *********************************/
43
44/*static double IMRPhenomT_MECOTime(double eta, double S, double dchi, double delta); FIXME: Promote some interesting quantities to XLAL functions*/
45
46static double IMRPhenomT_Inspiral_TaylorT3_t0(double eta, double S, double dchi, double delta); // theta = 0.45
47
48static double IMRPhenomT_Inspiral_Freq_CP1_22(double eta, double S, double dchi, double delta); // theta = 0.45
49static double IMRPhenomT_Inspiral_Freq_CP2_22(double eta, double S, double dchi, double delta);
50static double IMRPhenomT_Inspiral_Freq_CP3_22(double eta, double S, double dchi, double delta);
51static double IMRPhenomT_Inspiral_Freq_CP4_22(double eta, double S, double dchi, double delta);
52static double IMRPhenomT_Inspiral_Freq_CP5_22(double eta, double S, double dchi, double delta);
53static double IMRPhenomT_Merger_Freq_CP1_22(double eta, double S, double dchi, double delta);
54static double IMRPhenomT_PeakFrequency_22(double eta, double S, double dchi, double delta);
55static double IMRPhenomT_RD_Freq_D2_22(double eta, double S, double dchi, double delta);
56static double IMRPhenomT_RD_Freq_D3_22(double eta, double S, double dchi, double delta);
57
58/********************************* IMRPhenomT 22 Amplitude Fits *********************************/
59
60static double IMRPhenomT_Inspiral_Amp_CP1_22(double eta, double S, double dchi, double delta);
61static double IMRPhenomT_Inspiral_Amp_CP2_22(double eta, double S, double dchi, double delta);
62static double IMRPhenomT_Inspiral_Amp_CP3_22(double eta, double S, double dchi, double delta);
63static double IMRPhenomT_Merger_Amp_CP1_22(double eta, double S, double dchi, double delta);
64static double IMRPhenomT_PeakAmp_22(double eta, double S, double dchi);
65static double IMRPhenomT_RD_Amp_C3_22(double eta, double S);
66
67/********************************* IMRPhenomT 22 QNM Fits *********************************/
68
69static double evaluate_QNMfit_fring22(double finalDimlessSpin);
70
71static double evaluate_QNMfit_fdamp22(double finalDimlessSpin);
72
73static double evaluate_QNMfit_fdamp22n2(double finalDimlessSpin);
74
75/********************************* IMRPhenomTHM Fits *********************************/
76
77static double IMRPhenomT_Merger_Freq_CP1_21(double eta, double S, double dchi, double delta);
78static double IMRPhenomT_PeakFrequency_21(double eta, double S, double dchi, double delta);
79static double IMRPhenomT_Inspiral_Amp_CP1_21(double eta, double S, double dchi, double delta);
80static double IMRPhenomT_Inspiral_Amp_CP2_21(double eta, double S, double dchi, double delta);
81static double IMRPhenomT_Inspiral_Amp_CP3_21(double eta, double S, double dchi, double delta);
82static double IMRPhenomT_Merger_Amp_CP1_21(double eta, double S, double dchi, double delta);
83static double IMRPhenomT_PeakAmp_21(double eta, double S, double dchi, double delta);
84static double IMRPhenomT_RD_Amp_C3_21(double eta, double S, double dchi);
85static double IMRPhenomT_RD_Freq_D2_21(double eta, double S, double dchi, double delta);
86static double IMRPhenomT_RD_Freq_D3_21(double eta, double S, double dchi, double delta);
87
88static double IMRPhenomT_Merger_Freq_CP1_33(double eta, double S, double dchi, double delta);
89static double IMRPhenomT_PeakFrequency_33(double eta, double S, double dchi);
90static double IMRPhenomT_Inspiral_Amp_CP1_33(double eta, double S, double dchi, double delta);
91static double IMRPhenomT_Inspiral_Amp_CP2_33(double eta, double S, double dchi, double delta);
92static double IMRPhenomT_Inspiral_Amp_CP3_33(double eta, double S, double dchi, double delta);
93static double IMRPhenomT_Merger_Amp_CP1_33(double eta, double S, double dchi, double delta);
94static double IMRPhenomT_PeakAmp_33(double eta, double S, double dchi, double delta);
95static double IMRPhenomT_RD_Amp_C3_33(double eta, double S);
96static double IMRPhenomT_RD_Freq_D2_33(double eta, double S, double dchi, double delta);
97static double IMRPhenomT_RD_Freq_D3_33(double eta, double S, double dchi, double delta);
98
99static double IMRPhenomT_Merger_Freq_CP1_44(double eta, double S, double dchi, double delta);
100static double IMRPhenomT_PeakFrequency_44(double eta, double S, double dchi, double delta);
101static double IMRPhenomT_Inspiral_Amp_CP1_44(double eta, double S, double dchi, double delta);
102static double IMRPhenomT_Inspiral_Amp_CP2_44(double eta, double S, double dchi, double delta);
103static double IMRPhenomT_Inspiral_Amp_CP3_44(double eta, double S, double dchi, double delta);
104static double IMRPhenomT_Merger_Amp_CP1_44(double eta, double S, double dchi, double delta);
105static double IMRPhenomT_PeakAmp_44(double eta, double S, double dchi), double delta;
106static double IMRPhenomT_RD_Amp_C3_44(double eta, double S);
107static double IMRPhenomT_RD_Freq_D2_44(double eta, double S, double dchi, double delta);
108static double IMRPhenomT_RD_Freq_D3_44(double eta, double S, double dchi, double delta);
109
110static double IMRPhenomT_Merger_Freq_CP1_55(double eta, double S, double dchi, double delta);
111static double IMRPhenomT_PeakFrequency_55(double eta, double S, double dchi, double delta);
112static double IMRPhenomT_Inspiral_Amp_CP1_55(double eta, double S, double dchi, double delta);
113static double IMRPhenomT_Inspiral_Amp_CP2_55(double eta, double S, double dchi, double delta);
114static double IMRPhenomT_Inspiral_Amp_CP3_55(double eta, double S, double dchi, double delta);
115static double IMRPhenomT_Merger_Amp_CP1_55(double eta, double S, double dchi, double delta);
116static double IMRPhenomT_PeakAmp_55(double eta, double S, double dchi, double delta);
117static double IMRPhenomT_RD_Amp_C3_55(double eta, double S, double dchi);
118static double IMRPhenomT_RD_Freq_D2_55(double eta, double S, double dchi, double delta);
119static double IMRPhenomT_RD_Freq_D3_55(double eta, double S, double dchi, double delta);
120
121/********************************* IMRPhenomTHM QNM Fits *********************************/
122
123static double evaluate_QNMfit_fring21(double finalDimlessSpin);
124static double evaluate_QNMfit_fring33(double finalDimlessSpin);
125static double evaluate_QNMfit_fring44(double finalDimlessSpin);
126static double evaluate_QNMfit_fring55(double finalDimlessSpin);
127
128static double evaluate_QNMfit_fdamp21(double finalDimlessSpin);
129static double evaluate_QNMfit_fdamp33(double finalDimlessSpin);
130static double evaluate_QNMfit_fdamp44(double finalDimlessSpin);
131static double evaluate_QNMfit_fdamp55(double finalDimlessSpin);
132
133static double evaluate_QNMfit_fdamp21n2(double finalDimlessSpin);
134static double evaluate_QNMfit_fdamp33n2(double finalDimlessSpin);
135static double evaluate_QNMfit_fdamp44n2(double finalDimlessSpin);
136static double evaluate_QNMfit_fdamp55n2(double finalDimlessSpin);
137
138/***************************** IMRPhenomTHM Time Shifts ***********************/
139
140static double IMRPhenomT_tshift_21(double eta, double S, double dchi);
141static double IMRPhenomT_tshift_33(double eta, double S);
142static double IMRPhenomT_tshift_44(double eta, double S);
143static double IMRPhenomT_tshift_55(double eta, double S);
144
145#ifdef __cplusplus
146}
147#endif
148
149#endif // of #ifndef _LALSIM_IMR_PHENOMT_INSPIRAL_H
IMRPhenomT_Inspiral_Amp_CP1_33
static double IMRPhenomT_Inspiral_Amp_CP1_33(double eta, double S, double dchi, double delta)
IMRPhenomT_Inspiral_Freq_CP4_22
static double IMRPhenomT_Inspiral_Freq_CP4_22(double eta, double S, double dchi, double delta)
IMRPhenomT_PeakAmp_44
static double IMRPhenomT_PeakAmp_44(double eta, double S, double dchi)
IMRPhenomT_Inspiral_Amp_CP2_55
static double IMRPhenomT_Inspiral_Amp_CP2_55(double eta, double S, double dchi, double delta)
IMRPhenomT_Inspiral_Amp_CP3_33
static double IMRPhenomT_Inspiral_Amp_CP3_33(double eta, double S, double dchi, double delta)
evaluate_QNMfit_fdamp55
static double evaluate_QNMfit_fdamp55(double finalDimlessSpin)
IMRPhenomT_Inspiral_Amp_CP2_22
static double IMRPhenomT_Inspiral_Amp_CP2_22(double eta, double S, double dchi, double delta)
IMRPhenomT_RD_Amp_C3_33
static double IMRPhenomT_RD_Amp_C3_33(double eta, double S)
IMRPhenomT_PeakFrequency_22
static double IMRPhenomT_PeakFrequency_22(double eta, double S, double dchi, double delta)
IMRPhenomT_tshift_21
static double IMRPhenomT_tshift_21(double eta, double S, double dchi)
IMRPhenomT_PeakFrequency_44
static double IMRPhenomT_PeakFrequency_44(double eta, double S, double dchi, double delta)
IMRPhenomT_Inspiral_Amp_CP3_22
static double IMRPhenomT_Inspiral_Amp_CP3_22(double eta, double S, double dchi, double delta)
IMRPhenomT_RD_Amp_C3_21
static double IMRPhenomT_RD_Amp_C3_21(double eta, double S, double dchi)
IMRPhenomT_RD_Freq_D3_55
static double IMRPhenomT_RD_Freq_D3_55(double eta, double S, double dchi, double delta)
IMRPhenomT_PeakFrequency_21
static double IMRPhenomT_PeakFrequency_21(double eta, double S, double dchi, double delta)
evaluate_QNMfit_fdamp22n2
static double evaluate_QNMfit_fdamp22n2(double finalDimlessSpin)
evaluate_QNMfit_fdamp22
static double evaluate_QNMfit_fdamp22(double finalDimlessSpin)
IMRPhenomT_Inspiral_Amp_CP3_44
static double IMRPhenomT_Inspiral_Amp_CP3_44(double eta, double S, double dchi, double delta)
evaluate_QNMfit_fdamp21n2
static double evaluate_QNMfit_fdamp21n2(double finalDimlessSpin)
evaluate_QNMfit_fdamp44n2
static double evaluate_QNMfit_fdamp44n2(double finalDimlessSpin)
IMRPhenomT_Merger_Amp_CP1_44
static double IMRPhenomT_Merger_Amp_CP1_44(double eta, double S, double dchi, double delta)
IMRPhenomT_Merger_Amp_CP1_22
static double IMRPhenomT_Merger_Amp_CP1_22(double eta, double S, double dchi, double delta)
evaluate_QNMfit_fring33
static double evaluate_QNMfit_fring33(double finalDimlessSpin)
Definition: LALSimIMRPhenomXHM_qnm.c:64
IMRPhenomT_tshift_55
static double IMRPhenomT_tshift_55(double eta, double S)
IMRPhenomT_Inspiral_Amp_CP1_44
static double IMRPhenomT_Inspiral_Amp_CP1_44(double eta, double S, double dchi, double delta)
IMRPhenomT_RD_Freq_D2_21
static double IMRPhenomT_RD_Freq_D2_21(double eta, double S, double dchi, double delta)
IMRPhenomT_PeakAmp_22
static double IMRPhenomT_PeakAmp_22(double eta, double S, double dchi)
IMRPhenomT_Inspiral_Amp_CP2_21
static double IMRPhenomT_Inspiral_Amp_CP2_21(double eta, double S, double dchi, double delta)
IMRPhenomT_RD_Freq_D2_33
static double IMRPhenomT_RD_Freq_D2_33(double eta, double S, double dchi, double delta)
evaluate_QNMfit_fdamp33
static double evaluate_QNMfit_fdamp33(double finalDimlessSpin)
Definition: LALSimIMRPhenomXHM_qnm.c:82
IMRPhenomT_Merger_Freq_CP1_55
static double IMRPhenomT_Merger_Freq_CP1_55(double eta, double S, double dchi, double delta)
evaluate_QNMfit_fdamp33n2
static double evaluate_QNMfit_fdamp33n2(double finalDimlessSpin)
IMRPhenomT_RD_Amp_C3_44
static double IMRPhenomT_RD_Amp_C3_44(double eta, double S)
evaluate_QNMfit_fdamp44
static double evaluate_QNMfit_fdamp44(double finalDimlessSpin)
Definition: LALSimIMRPhenomXHM_qnm.c:151
IMRPhenomT_Merger_Amp_CP1_33
static double IMRPhenomT_Merger_Amp_CP1_33(double eta, double S, double dchi, double delta)
IMRPhenomT_RD_Freq_D3_22
static double IMRPhenomT_RD_Freq_D3_22(double eta, double S, double dchi, double delta)
IMRPhenomT_Merger_Freq_CP1_22
static double IMRPhenomT_Merger_Freq_CP1_22(double eta, double S, double dchi, double delta)
IMRPhenomT_PeakFrequency_33
static double IMRPhenomT_PeakFrequency_33(double eta, double S, double dchi)
IMRPhenomT_Inspiral_Freq_CP3_22
static double IMRPhenomT_Inspiral_Freq_CP3_22(double eta, double S, double dchi, double delta)
evaluate_QNMfit_fring22
static double evaluate_QNMfit_fring22(double finalDimlessSpin)
IMRPhenomT_Inspiral_Amp_CP2_44
static double IMRPhenomT_Inspiral_Amp_CP2_44(double eta, double S, double dchi, double delta)
IMRPhenomT_Inspiral_Amp_CP3_55
static double IMRPhenomT_Inspiral_Amp_CP3_55(double eta, double S, double dchi, double delta)
IMRPhenomT_RD_Freq_D2_22
static double IMRPhenomT_RD_Freq_D2_22(double eta, double S, double dchi, double delta)
IMRPhenomT_Merger_Freq_CP1_21
static double IMRPhenomT_Merger_Freq_CP1_21(double eta, double S, double dchi, double delta)
IMRPhenomT_Inspiral_Freq_CP2_22
static double IMRPhenomT_Inspiral_Freq_CP2_22(double eta, double S, double dchi, double delta)
IMRPhenomT_RD_Freq_D3_21
static double IMRPhenomT_RD_Freq_D3_21(double eta, double S, double dchi, double delta)
IMRPhenomT_RD_Amp_C3_55
static double IMRPhenomT_RD_Amp_C3_55(double eta, double S, double dchi)
IMRPhenomT_Merger_Amp_CP1_21
static double IMRPhenomT_Merger_Amp_CP1_21(double eta, double S, double dchi, double delta)
IMRPhenomT_PeakAmp_55
static double IMRPhenomT_PeakAmp_55(double eta, double S, double dchi, double delta)
IMRPhenomT_RD_Freq_D3_33
static double IMRPhenomT_RD_Freq_D3_33(double eta, double S, double dchi, double delta)
IMRPhenomT_tshift_44
static double IMRPhenomT_tshift_44(double eta, double S)
evaluate_QNMfit_fdamp21
static double evaluate_QNMfit_fdamp21(double finalDimlessSpin)
Definition: LALSimIMRPhenomXHM_qnm.c:47
evaluate_QNMfit_fring55
static double evaluate_QNMfit_fring55(double finalDimlessSpin)
IMRPhenomT_PeakAmp_21
static double IMRPhenomT_PeakAmp_21(double eta, double S, double dchi, double delta)
IMRPhenomT_Inspiral_Freq_CP5_22
static double IMRPhenomT_Inspiral_Freq_CP5_22(double eta, double S, double dchi, double delta)
evaluate_QNMfit_fring21
static double evaluate_QNMfit_fring21(double finalDimlessSpin)
Definition: LALSimIMRPhenomXHM_qnm.c:30
IMRPhenomT_Merger_Freq_CP1_44
static double IMRPhenomT_Merger_Freq_CP1_44(double eta, double S, double dchi, double delta)
IMRPhenomT_PeakAmp_33
static double IMRPhenomT_PeakAmp_33(double eta, double S, double dchi, double delta)
IMRPhenomT_Inspiral_Amp_CP1_21
static double IMRPhenomT_Inspiral_Amp_CP1_21(double eta, double S, double dchi, double delta)
IMRPhenomT_Merger_Freq_CP1_33
static double IMRPhenomT_Merger_Freq_CP1_33(double eta, double S, double dchi, double delta)
IMRPhenomT_RD_Freq_D2_55
static double IMRPhenomT_RD_Freq_D2_55(double eta, double S, double dchi, double delta)
IMRPhenomT_RD_Freq_D2_44
static double IMRPhenomT_RD_Freq_D2_44(double eta, double S, double dchi, double delta)
delta
static double double delta
Definition: LALSimIMRPhenomTHM_fits.h:105
IMRPhenomT_PeakFrequency_55
static double IMRPhenomT_PeakFrequency_55(double eta, double S, double dchi, double delta)
IMRPhenomT_tshift_33
static double IMRPhenomT_tshift_33(double eta, double S)
IMRPhenomT_RD_Freq_D3_44
static double IMRPhenomT_RD_Freq_D3_44(double eta, double S, double dchi, double delta)
IMRPhenomT_Inspiral_Freq_CP1_22
static double IMRPhenomT_Inspiral_Freq_CP1_22(double eta, double S, double dchi, double delta)
IMRPhenomT_Inspiral_Amp_CP1_22
static double IMRPhenomT_Inspiral_Amp_CP1_22(double eta, double S, double dchi, double delta)
IMRPhenomT_RD_Amp_C3_22
static double IMRPhenomT_RD_Amp_C3_22(double eta, double S)
evaluate_QNMfit_fring44
static double evaluate_QNMfit_fring44(double finalDimlessSpin)
Definition: LALSimIMRPhenomXHM_qnm.c:133
evaluate_QNMfit_fdamp55n2
static double evaluate_QNMfit_fdamp55n2(double finalDimlessSpin)
IMRPhenomT_Inspiral_Amp_CP2_33
static double IMRPhenomT_Inspiral_Amp_CP2_33(double eta, double S, double dchi, double delta)
IMRPhenomT_Merger_Amp_CP1_55
static double IMRPhenomT_Merger_Amp_CP1_55(double eta, double S, double dchi, double delta)
IMRPhenomT_Inspiral_TaylorT3_t0
static double IMRPhenomT_Inspiral_TaylorT3_t0(double eta, double S, double dchi, double delta)
IMRPhenomT_Inspiral_Amp_CP3_21
static double IMRPhenomT_Inspiral_Amp_CP3_21(double eta, double S, double dchi, double delta)
IMRPhenomT_Inspiral_Amp_CP1_55
static double IMRPhenomT_Inspiral_Amp_CP1_55(double eta, double S, double dchi, double delta)