IMRPhenomTP angular momentum bug for PN Euler angles
Hi @hector.estelles,
I think I have found a small bug in PhenomTP(HM) when using the PN Euler angles (not the default version).
When PhenomTP calls IMRPhenomXGetAndSetPrecessionVariables
it first needs to set up the powers_of_lalpi so that it computes the correct PN terms for the angular momentum, see e.g. here. Currently the terms containing powers_of_pi
are zero.
Just by adding
status = IMRPhenomX_Initialize_Powers(&powers_of_lalpi, LAL_PI);
XLAL_CHECK(XLAL_SUCCESS == status, status, "Failed to initialize useful powers of LAL_PI.");
for example here it would be fixed.