fit_uncertainty_ellipse

dysmalpy.utils.fit_uncertainty_ellipse(chain_x, chain_y, bins=50)[source]

Get the uncertainty ellipse of the sample for each photocenter. (Modified from code from Jinyi Shangguan to do photocenter uncertainty ellipses)

Parameters:

  • pos_list (2D array) – Sampler chains for photocenter positions

  • bins (integer) – The number of bins to use in the 2D histogram

Returns:

  • PA (float) – angle of ellipse, in degrees

  • stddev_x (float) – stddev of the “x” axis of the 2D gaussian; double to get the full “width” of a 1sig ellipse for matplotlib.Ellipse

  • stddev_y (float) – stddev of the “y” axis of the 2D gaussian