jfactor#

jeanspy.axisymmetric_factors.jfactor

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

Return integral rho^2 ds dOmega in GeV^2 cm^-5 (gamma must be < 1.5).

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^2 cm^-5; 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<1.5 (finite central annihilation integral). 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