BFieldRange = [2.45 2.44 2.43 2.42 2.40 2.39 2.38 2.37 2.36 2.35 2.34 2.32 2.30 2.28 2.25 2.20 2.15 2.10 2.05 2.00 1.95 1.90 1.85 1.80];
AlphaRange  = [180 175 170 168 166 164 162 160 158 156 154 152 150 145 140];

% Z must be 15x24 to match y-by-x
Z = zeros([numel(BFieldRange), numel(AlphaRange)]);... % your data

[X,Y] = meshgrid(x, y);

%% ------------------ Plot Heatmaps ------------------

% Heatmap of radial_spectral_contrast
Plotter.plotHeatmap(compiled_results, scan_reference_values, BFields, 'radial_spectral_contrast', ...
    'Colormap', @sky, ...
    'CLim', [0 0.008], ...
    'XLabel', '\alpha (degrees)', ...
    'YLabel', 'BField (G)', ...
    'Title', 'Radial Spectral Contrast', ...
    'FigNum', 10, ...
    'SaveFileName', 'Heatmap_RadialSpectralContrast.fig', ...
    'SaveDirectory', options.saveDirectory);

% Heatmap of mean_max_g2_values
Plotter.plotHeatmap(compiled_results, scan_reference_values, BFields, 'mean_max_g2_values', ...
    'Colormap', @sky, ...
    'CLim', [0 1], ...
    'XLabel', '\alpha (degrees)', ...
    'YLabel', 'BField (G)', ...
    'Title', '$\mathrm{max}[g^{(2)}_{[50,70]}(\delta\theta)]$', ...
    'FigNum', 9, ...
    'SaveFileName', 'Heatmap_MaxG2.fig', ...
    'SaveDirectory', options.saveDirectory);