Image modulation profiles

We compare the modulation profiles predicted by the image to the observed modulation profiles to determine how good the image is (see also the viscomp... parameters for comparison of observed and expected visibilities). There are two tools - one computes the goodness of fit between the predicted and observed profiles, as well as the fit coefficients. The other performs a Monte Carlo analysis by varying the predicted profile, xp, by poisson statistics, calculating the goodness of fit of the varied xp to xp, doing this a set number of times, fitting a Gaussian to the distribution, and returning the center of the Gaussian and the FWHM.


The following control parameters are available for these tests:

profile_show_plot - If set, show plots of predicted and observed profiles while making image. Default is 0. (plots on screen or in PS or JPEG file depending on profile_ps_plot and profile_jpeg_plot)

profile_ps_plot - If set, draw profile plot in PS file (auto-generated name)

profile_jpeg_plot - If set, draw profiles in JPEG file (auto-generated name)

profile_plot_dir - Directory to write profile plot in (defaults to current directory)

profile_plot_rate - Units for profile plots. 0 means plot counts/bin; 1 means plot counts/sec. Default is 1

profile_plot_resid - If set, plot profile residuals. Default=0

profile_window - Internal variable to keep track of window to plot profiles in. If plotting on screen, and you want to plot in a new window (to preserve the previous plot instead of overwriting in that window), set to -1.

mc_ntrials - Number of trials to run in the Monte Carlo analysis

mc_show_plot - If set, show plots of Monte Carlo distribution of C-statistic while making image. Note - this only has an effect when profile_show_plot is 1, or when profile plots are requested using the hsi_image_profile_plot routine after an image was made.

To display the profile and Monte Carlo plots after an image is made, use the hsi_image_profile_plot routine with the appropriate keywords.


Info parameters resulting from profile tests:
(Note: Replace x by the algorithm prefix - clean, em, pixon, ff, ge, vf, nj, uv, vis_cs, or vis_wv)

x_profile_tot_cstat - Overall C-statistic from the comparison of expected to observed modulation profiles for all selected detectors.

x_profile_cstat - C-statistic from the comparison of each separate detector's expected and observed modulation profile.

x_profile_coeff - Constant, slope, and correlation coefficient between the expected profile and observed profile for each detector. Apply the constant and slope to the expected profile to match the observed profile. Note that the constant is in units of counts / bin.The correlation coefficient gives an overall estimate of the correlation between the expected and observed profiles.

x_profile_mc_tot_cstat - Most probable C-statistic for all selected detectors, from Monte Carlo analysis

x_profile_mc_cstat - Most probable C-statistic for each separate selected detector, from Monte Carlo analysis

x_profile_mc_tot_fwhm - FWHM of Monte Carlo distribution for all selected detectors

x_profile_mc_fwhm - FWHM of Monte Carlo distribution fo each separate selected detector


Limitations based on availability of calibrated eventlist (and therefore the profiles):

1. Image cubes for visibility-based algorithms can not make this comparison.

2. For image cubes for non-vis algorithms, profiles can only be plotted while generating the images,not afterwards.