Skip to content

AGN jet ray sorting

Filip Husko requested to merge AGN_jet_ray_sorting into master

This MR relates to the discussion of #824 (closed). It implements a simpler sorting of particles into two rays. All particles are now sorted in both rays, going from one side of the spin axis to the other. This ensures there are no cases with one ray being empty. The lines that do this are:

      quantity_to_minimize = cosine_theta + ray_jet_correction;
      quantity_to_minimize_pos = -cosine_theta + ray_jet_correction;

Unfortunately, this only works for the case where we launch particles along the spin axis. If we choose one of the other launching models (minimum/maximum distance or minimum density), it's not clear what's best to do in the case of empty hemispheres. I've given it some thought but haven't come up with anything for now. It appears we need to still do the hemisphere approach in those other cases.

Merge request reports