dir = './tests/'

img_file        = dir+'dv3.XMMLSS_00-fd_v0-diff.fits.gz'
cat_file        = dir+'CFHTLS_D1_wth.cat'
fits_mask_file  = dir+'threshold-mask.fits.gz'
ascii_mask_file = dir+'CFHTLS_D1_XMMLSS_masks.mask'

;center of prior cat field
;needed to put objects at equator and then use small angles approximation
radec_cent = [36.659,-4.485]

;Correlation function computation parameters
theta_min = 0.005D
theta_max = 0.5D
nbins     = 15

;auto correlation of prior cat
out_prior_wtheta_file = dir+'prior.wth'
prior_wth = wtheta_fft(cat_file, out_prior_wtheta_file, $
                       ascii_mask_file, radec_cent = radec_cent, rmin = $
                       theta_min, rmax = theta_max, nbin = nbins, /verbose, $
                       /display)


;auto correlation of diff image
out_diff_wtheta_file = dir+'diff.wth'
diff_wth = wtheta_fft(img_file, out_diff_wtheta_file, $
                      fits_mask_file, rmin = theta_min, rmax = theta_max, $
                      nbin = nbins, custom_res = 1.5, xstart = 340, $
                      ystart = 340, /verbose, /display)


;cross correlation of prior cat and diff image
;here using the same mask for both inputs, but different masks are ok
out_prior_diff_xwtheta_file = dir+'prior_diff.xwth'
prior_diff_xwth = xwtheta_fft(cat_file, img_file,  $
                              out_prior_diff_xwtheta_file, fits_mask_file, $
                              fits_mask_file, radec_cent = radec_cent,$
                              rmin = theta_min, rmax = theta_max, $
                              nbin = nbins, custom_res = 1.5, xstart = 340, $
                              ystart = 340, /verbose, /display)

END