dfactor#

jeanspy.axisymmetric_factors.dfactor

jeanspy.axisymmetric_factors.dfactor(halo, dist_pc, roi_deg, *, inclination=1.5707963267948966, n_mu=96, n_phi=96, n_radial=128)[source]#

Return integral rho ds dOmega in GeV cm^-2.

Notes

Inputs and units. halo is AxisymmetricZhaoModel; dist_pc is scalar observer distance (pc); roi_deg is scalar cone half-angle; inclination is radians; n_mu/n_phi/n_radial are integers >=16.

Returns and shape. Nonnegative Python float in GeV cm^-2; zero aperture gives zero.

Validity. Require explicit finite r_t_pc, observer distance > r_t_pc*max(1,Q), 0 <= roi_deg < 90 and gamma<2 under the AxisymmetricZhaoModel constructor domain. n_phi uses a periodic rule; refine all orders.

Errors. Invalid geometry/domain/order raises ValueError; wrong halo type raises TypeError.

Backend. NumPy/SciPy CPU postprocessing.

Differentiation. No physical-parameter automatic differentiation on this API.

Examples. examples/docs_factors.py