For this figure, all data is split into tope and bottom rows
For the top row, the combined symbol <N> stands for 'small_displacement'. For the bottom row, it stands for 'large_displacement'
The symbol appears in the name of every data field

Data fields:

<N>_waist_xx: The transverse coordinates of the basis centre with respect to the cavity centre in the spherical geometry scan. Units: metres
<N>_waist_ww: The central waist values of the basis in the spherical geometry scan. Units: metres
<N>_spherical_cint: The C_{int}/L_{y, eff}^{-1} values of the cavities in the spherical mirror geometry scan. Units: dimensionless
<N>_best_x_centre_wrt_cavity: The transverse coordinate of the basis whose corresponding spherical mirror gives the highest C_{int}/L_{y, eff}^{-1} of the spherical geometry scan. Units: metres
<N>_best_w_0: The central waist value of the basis whose corresponding spherical mirror gives the highest C_{int}/L_{y, eff}^{-1} of the spherical geometry scan. Units: metres
<N>_x_emitter_wrt_cavity: The transverse coordinate of the emitter (i.e. its displacement from the centre of the cavity). Units: metres
FOR <N> == 'large_displacement' only:
	'large_displacement_resonant_waists': The central waist values of the bases with noted transverse degeneracies. Units: metres

<N>_density_zz: Axial coordinates of the intensity values for the density plot. Units: metres
<N>_density_xx: Transverse x coordinates of the intensity values for the density plot. Units: metres
<N>_spherical_mirror_mode_density: Intensities of the mode with the spherical mirror that maximises C_{int}^{min} for the density plot. Units: inverse metres squared
<N>_optimised_mirror_mode_density: Intensities of the mode with the optimised for the density plot. Units: inverse metres squared 
<N>_length: The length of the cavity in the density plot. Units: metres
<N>_mirror_diameter: The diameter of the cavity mirrors in the density plot. Units: metres

<N>_surface_r: Radial coordinate values for the surface plot of the mirror. Units: metres
<N>_spherical_protrusion : The protrusion of the mirror surface that maximises C_{int} in the spherical case. Units: metres
<N>_surface_protrusion: The protrusion of the retroreflective optimised mirror. Units: metres
<N>_residuals_protrusion: The residuals of the retroreflective optimised surface from the spherical surface that maximises C_{int} (i.e. residuals_protrusion = surface_protrusion - spherical_protrusion). Units: metres


Figure plots:
Panel a:
Heat map: <N>_spherical_cint vs. <N>_waist_xx (y axis) and <N>_waist_ww (x axis)
Orange cross: <N>_best_x_centre_wrt_cavity vs <N>_best_w_0
Orange dashed line: <N>_best_x_centre_wrt_cavity
Green horizontal line: <N>_x_emitter_wrt_cavity
FOR <N> == 'large_displacement' only:
	Vertical white dashed lines: a w_0 values in large_displacement_resonant_waists

Panel b:
top:
Heat map: <N>_spherical_mirror_mode_density (colour) vs <N>_density_zz (x axis) and <N>_density_xx (y axis)
bottom:
Heat map: <N>_optimised_mirror_mode_density (colour) vs <N>_density_zz (x axis) and <N>_density_xx (y axis)
both:
Cyan cross: axial coordinate of 0, transverse coordinate of <N>_x_emitter_wrt_cavity
White lines: vertical lines of length <N>_mirror_diameter, centred on transverse coordinate of 0 with axial coordinates of 0.5 * <N>_length and -0.5 * <N>_length

Panel c:
Black line: surface_r vs spherical protrusion (identical data plotted for negative r)
Blue line: surface_r vs surface_protrusion (identical data plotted for negative r)
Red line: surface_r vs residuals_protrusion (identical data plotted for negative r)
Orange dashed vertical line: <N>_best_x_centre_wrt_cavity
Green vertical line: <N>_x_emitter_wrt_cavity
