mosaic_findpeak

 peak=mosaic_findpeak(ad,mos,sign,o)

 Find max peak/trough of mosaic significance map

 sign=+1 for peak, -1 for trough

0001 function [peak noiseVal]=mosaic_findpeak(ad,mos,sign,o)
0002 % peak=mosaic_findpeak(ad,mos,sign,o)
0003 %
0004 % Find max peak/trough of mosaic significance map
0005 %
0006 % sign=+1 for peak, -1 for trough
0007 
0008 % check to see if we have to cut out holes in our data
0009 if(isfield(mos, 'sigmaHole'))
0010     if(isempty(mos.sigmaHole))
0011         % do nothing;
0012     else
0013         mos.sigma = mos.sigmaHole;
0014     end
0015 end
0016 
0017 % Find the peak significance value
0018 [maxsig,maxind]=maxmd(mos.sigma*sign);
0019 
0020 peak.sigma=maxsig*sign;
0021 peak.flux=mos.signal(maxind(1),maxind(2));
0022 noiseVal = mos.noise(maxind(1), maxind(2));
0023 
0024 % Return pixel values, too
0025 peak.xpix=maxind(2);
0026 peak.ypix=maxind(1);
0027 
0028 % Convert pixel numbers to deg
0029 if(isfield(ad, 't_val_x_deg'))
0030     peak.x = ad.t_val_x_deg(maxind(2));
0031 else
0032     peak.x=ad.t_val_deg(maxind(2));
0033 end
0034 peak.y=ad.t_val_deg(maxind(1));
0035 
0036 % Convert degrees to RA/DEC
0037 [peak.RA peak.DEC] = fieldoff_to_radec(peak.x, peak.y, o.rafc*15,o.decfc);
0038 [peak.RAstr peak.DECstr] = radec2str(peak.RA/15, peak.DEC);
0039 if exist('o')
0040   [peak.RA peak.DEC] = fieldoff_to_radec(peak.x, peak.y, o.rafc*15, ...
0041                                          o.decfc);
0042   [peak.RAstr peak.DECstr] = radec2str(peak.RA/15, peak.DEC);
0043 end
0044 
0045 return