Skip to content

IMRPhenomX settings and nestcheck test

Maite Mateu-Lucena requested to merge maite.mateu-lucena/bilby:imrphenomx-nestcheck into master

This MR adds the nestcheck sampler setting to dynesty. If it is nestcheck = True, a pickle file is created for each run in the output directory storing a nestcheck.data_processing.process_dynesty_run object, (nestcheck needs to be installed). That object can be used then in the post processing in order to compute the implementation and bootstrap errors explained in “Sampling Errors In Nested Sampling Parameter Estimation” (Higson et al. 2018).

Other waveform_kwargs of the IMRPhenomX family are also added in that MR:

phenomXHMMband: Threshold for the multibanding algorithm in IMRPhenomXHM.
                If option not specified use 10^-3. If set to 0 then multibanding is disabled.
phenomXPHMMband: Threshold for the multibanding algorithm for the Euler angles in IMRPhenomXPHM.
                If option not specified use 0. If set to 0 then multibanding is disabled.
phenomXPrecVersion: Version for the Euler angles of the twisting-up in IMRPhenomXP/IMRPhenomXPHM.
                101 : NNLO PN Euler angles and a 2PN non-spinning approximation to L
                102 : NNLO PN Euler angles and a 3PN spinning approximation to L
                103 : NNLO PN Euler angles and a 4PN spinning approximation to L
                104 : same than 103 augmented with leading PN order at all order in spin terms.
                220 : MSA Euler angles and a 3PN spinning approximation to L. Falls back to NNLO angles in case of fail.
                221 : MSA Euler angles and a 3PN spinning approximation to L. Throw an error message in case of fail.
phenomXPConvention: Version for computing the offsets of the Euler Angles.
                2: recover aligned spin limit.
                7: PhenomPv3HM style needs to optimize to recover aligned spin limit.
phenomXPFinalSpinMod: Version for the final spin in the precessing waveform.
                0: use chip
                1: use chi1x
                2: use chi_perp
phenomXPHMTwistPhenomHM: Twist PhenomHM instead of PhenomXHM
                0: twist PhenomXHM
                1: twist PhenomHM

Finally numerical relativity simulations can be injected using numrelfile: path to LVCNR h5.

Merge request reports