diff --git a/waveform/wrapper_mpv_waveform/mpvwaveform.py b/waveform/wrapper_mpv_waveform/mpvwaveform.py
index fa04e33101e089e8e35913c238e66279adbd88d4..1b40ad7fa90eda7c5a5a97c2b6c88d1ee43cef3b 100644
--- a/waveform/wrapper_mpv_waveform/mpvwaveform.py
+++ b/waveform/wrapper_mpv_waveform/mpvwaveform.py
@@ -25,7 +25,7 @@ def gen(**kwds):
 
     hp, hc = get_fd_waveform(approximant=kwds['baseapprox'], **kwds)
     zz = cosmology.redshift(kwds['distance'])
-    intz = integrate.quad(integrand_parityamu_mpv, 0, zz)[0]
+    intz = integrate.quad(integrand_parityamu_mpvinverse, 0, zz)[0]
     temp =  kwds['parity_mpvinverse'] * intz / 1e9 / lal.QE_SI * lal.H_SI * lal.PI * lal.PI / lal.H0_SI
 
     hp_parity = hp + hc * temp * hp.sample_frequencies **2