Imrphenomx settings

This MR add 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 numerical_relativity_file: path to LVC Numerical Relativity Simulations (h5 files, https://git.ligo.org/waveforms/lvcnr-lfs).