It seems a bit silly to convert the flux weighted by the Gaussian to an amplitude, only to convert it back to weighted flux, when it can be done by passing a pointer directly to the `calcmom` function. It avoids having to divide and multiply by factors involving PI, although the loss in precision is negligible.