the study of galaxy formation and evolution requires detailed information about statistically significant samples of dim objects .
this , in turn , requires deep imaging and spectroscopy over wide areas of the sky . in pursuit of these data ,
several wide - deep surveys are now underway .
those covering several square degrees or more either lack spectroscopic follow - up ( e.g. noao dwfs , @xcite , and odts , @xcite ) or are restricted to the study of objects at @xmath9 ( except for quasars ) by their imaging depth e.g. the sloan digital sky survey ( @xcite and @xcite ) , the virmos - vlt deep survey ( @xcite , @xcite ) and deep2 @xcite .
other surveys target the high - redshift universe with deep hst imaging over fractions of a square degree i.e. the hubble deep fields @xcite , the hubble ultra deep field ( hudf ) , goods @xcite , and gems @xcite .
spectroscopic coverage is feasible over these areas but is presently unable to probe deeper than @xmath10 , making the added imaging depth useful only for morphological studies and photometric redshifts .
the multiwavelength survey by yale - chile ( musyc ) probes the intermediate regime of a square degree of sky to the spectroscopic limit of @xmath10 .
section 2 describes the design of our survey .
 3 reports our imaging observations for the extended hubble deep field south .
 4 describes our imaging reduction , and  5 covers our photometric calibration and photometry .
 6 gives our results for @xmath4-band number counts and the sky density and angular clustering of @xmath6-selected lyman break galaxies .
 7 concludes .
our analyses assume a standard @xmath11cdm cosmology with @xmath12=0.27,@xmath13 , and @xmath14=70km / s / mpc .
the multiwavelength survey by yale - chile ( see http://www.astro.yale.edu/musyc ) is designed to provide a fair sample of the universe for the study of the formation and evolution of galaxies and their central black holes .
the core of the survey is a deep imaging campaign in optical and near - infrared passbands of four carefully selected @xmath15 fields .
musyc is unique for its combination of depth and total area , for additional coverage at x - ray , uv , mid- and far - infrared wavelengths and for providing the @xmath16 photometry needed for high - quality photometric redshifts over a square - degree of sky .
the primary goal is to study the properties and interrelations of galaxies at a single epoch corresponding to redshift @xmath17 , using a range of selection techniques .
we chose to use the @xmath0 filter set in the optical in order to obtain six nearly - independent flux measurements with the broadest possible wavelength coverage .
lyman break galaxies at @xmath5 are selected through their dropout in @xmath18-band images combined with blue continuua in @xmath19 ( @xmath20 in the rest - frame ) typical of recent star formation @xcite .
imaging depths of @xmath18,@xmath21,@xmath22,@xmath23 were chosen to detect the lbgs , whose luminosity function has a characteristic magnitude of @xmath24 in @xmath25 , and to find their lyman break decrement in the @xmath18 filter via colors @xmath26 .
lyman @xmath27 emitters at @xmath5 are selected through additional deep narrow - band imaging using a 50 fwhm filter centered at 5000 .
these objects can be detected in narrow - band imaging and spectroscopy by their emission lines , allowing us to probe to much dimmer continuum magnitudes than possible for lyman break galaxies .
it has recently become clear that optical selection methods do not provide a full census of the galaxy population at @xmath28 , as they miss objects which are faint in the rest - frame ultraviolet @xcite . with this in mind , musyc has a comprehensive near - infrared imaging campaign .
the nir imaging comprises two components : a wide survey covering the full square degree and a deep survey of the central @xmath29 of each field .
this division between deep and wide was chosen because of the @xmath30 field - of - view of the ispi near - infrared camera on the ctio 4 m telescope .
the @xmath31 point source sensitivities of the wide and deep components are @xmath32 and @xmath33 respectively .
nir imaging over the full survey area provides a critical complement to optical imaging for breaking degeneracies in photometric redshifts and modeling star formation histories .
deeper @xmath34 imaging over @xmath29 subfields opens up an additional window into the @xmath35 universe as the @xmath36 selection technique @xcite will be used to find evolved optically - red galaxies at @xmath37 through their rest - frame balmer/4000 break .
extensive follow - up spectroscopy is being conducted over the square degree .
a subset of the color - selected lyman break galaxy candidates will turn out to be agn based on broad- or narrow - line emission features seen in follow - up spectroscopy @xcite . damped lyman @xmath27 absorption systems @xcite at @xmath38 , which comprise the neutral gas reservoir needed to form most of the stars in the universe @xcite ,
will be searched for in the spectra of the brightest color - selected lbg / agn candidates ( typically quasars at @xmath28 ) .
in addition to the optical and near - infrared , imaging campaigns at other wavelengths and follow - up spectroscopy are integral parts of musyc .
x - ray selection will be used to study agn demographics over the full range of accessible redshifts , @xmath39 , ( see * ? ? ?
* ) with spitzer imaging used to detect optically- and x - ray - obscured agn @xcite .
this also allows a census of accreting black holes at @xmath5 in the same fields to study correlations between black hole accretion and galaxy properties at this epoch .
future epochs of optical imaging will be used to conduct a proper motion survey to find white dwarfs and brown dwarfs in order to study galactic structure and the local initial mass function ( see * ? ? ?
* ) ; the additional epochs will also enable a variability study of agn .
the four survey fields ( see table  [ tab : fields ] ) were chosen to have extremely low reddening , hi column density @xcite , and 100@xmath40 m dust emission @xcite in order to facilitate satellite coverage with spitzer , hst , chandra , and xmm , to take advantage of existing multiwavelength data and to enable flexible scheduling of observing time during the year .
additionally , each field satisfies all of the following criteria : minimal bright foreground sources in the optical and radio , high galactic latitude ( @xmath41 ) to reduce stellar density , and accessibility from observatories located in chile .
the survey fields will be a natural choice for future observations with alma .
the remainder of this paper describes our optical images and catalog of e - hdfs .
the techniques used for data reduction and photometry are the same as those used for the analysis of the other three fields . optical imaging from the full survey
will be reported in e.  gawiser et al .
( 2005 , in preparation ) .
the near - infrared data will be discussed in r.  quadri et al .
( 2005 , in preparation ) .
the e - hdfs has deep public space - based observations at uv , optical , near - infrared , and far - infrared wavelengths .
the hdfs itself covers a small @xmath42 central region with wfpc2 plus stis and nicmos regions , with deep ground - based jhk coverage of the wfpc2 region available from the fires survey @xcite .
spitzer irac and mips coverage of the central @xmath43 is being performed in gto time .
the extended area around hdfs has previously been imaged by @xcite and @xcite to a depth sufficient for the study of galaxies at @xmath44 .
these images were made public and combined with deep h - band images by the las campanas infrared survey ( lcirs , @xcite ) to study red galaxies out to @xmath45 .
our survey goes about one magnitude deeper in ubvri@xmath46 to probe the @xmath2 universe .
our extended hubble deep field south ( e - hdfs ) field center ( see table  [ tab : fields ] ) was chosen to keep a bright star ( m=6.8 ) which lies just north of the wfpc2 field off of the ctio+mosaic detectors .
lcirs covered an h - shape centered on wfpc2 and thereby provides public @xmath47-band coverage of roughly half of our e - hdfs field .
optical images of e - hdfs were taken on the nights of 2002 october 6,8,10,12 and 2003 may 26,27,28 using the 8192@xmath488192 pixel mosaic ii camera consisting of 8 2048@xmath484096 ccds , each with 2 amplifiers , on the blanco 4 m telescope at ctio .
afternoon calibrations were obtained , including zero - second exposures to trace the readout bias pattern , and domeflats in each filter to be observed except for @xmath18 for which the counts from the dome lamps were insufficient .
twilight flats were therefore obtained for the @xmath18 filter .
dark exposures of comparable length to the imaging observations were obtained , but due to the negligible dark current they were not used for data reduction .
a standard dither pattern was used to fill in gaps between the 8 ccds .
the pixel scale is 0.267@xmath49/pixel , leading to coverage spanning @xmath50 of sky with each pointing .
figure [ fig : filters ] shows the filter response curves and their multiplication with the ccd quantum efficiency and atmospheric transmission at one airmass .
table [ tab : obs ] gives details of the exposure times in each filter for each run along with the approximate average seeing measured in the raw images during observations .
the seven nights were mostly cloudless , but moderate clouds affected some of the @xmath22 imaging of e - hdfs on 2002 october 10,12 and some of the @xmath4 imaging on 2003 may 28 ; our reduction methods described below allow these images to be used without biasing the photometry .
photometric standard fields from @xcite large enough to cover the full mosaic ii field of view were observed each night , and the nights of 2002 october 6 and 2003 may 26 proved photometric .
@xmath0 images from ctio4m+mosaic ii were reduced using the mscred and mscdb packages in iraf v.2.12 following the noao deep wide field survey cookbook v.7.02 .
we used custom software to work around a few difficulties in these packages , as described below .
a composite zero image is subtracted from each raw image to remove the amplifier bias level and pattern .
the resulting image is then `` flat - fielded '' by dividing by a composite domeflat , or a composite twilight flat in the case of @xmath18-band .
superskyflat _ for each filter is made by combining all of the flat - fielded , unregistered images taken in each filter each night with rejection used to remove sources .
our 2@xmath51 amplitude dither pattern was designed to eliminate the wings of bright sources from the _
superskyflat_. each flat - fielded image is then divided by the appropriate _ superskyflat _ , which offers an estimate of the pixel - by - pixel response to the spectrum of the night sky with sufficient counts to achieve 1% precision per pixel . because the _ superskyflat _ was produced using flat - fielded images , dividing by it serves to remove the original domeflat ( twilight flat ) from the reduction process .
the real influence of the original flat - fielding is to remove the illumination pattern and gross pixel - to - pixel variations before looking for cosmic rays and bright objects to reject in making the _ superskyflat_. using the _ superskyflat _ to correct for the pixel response is considered preferable to using the domeflat because the ccd response to the spectrum of the dome lamp ( or twilight sky ) can have significant systematic differences from its response to the spectrum of the night sky . for background - limited photometry , this is an important effect .
we then find an astrometric solution for each image , starting with fiducial wcs headers provided for each fits extension of the raw images and comparing the claimed positions with the known positions of stars in the usno - b catalog using the mscred routine _
msccmatch_. we found _ msccmatch _ to be finicky ; for some runs the fiducial headers were too inaccurate to be corrected with the maximum second - order terms used by the mscred package .
an iterative , non - interactive procedure of calling _ msccmatch _ multiple times proved sufficient .
the final rms astrometric errors are between 0.2@xmath49 and 0.3@xmath49 in each image , which is consistent with the uncertainties on individual usno - b stars despite having fit many more stars than free parameters , and the actual solution should be better than 0.2@xmath49 rms .    to perform aperture photometry later on ,
our final images are transformed to have a common pixel scale and tangent plane projection point .
this is accomplished by projecting each processed image ( after bias - subtraction , flat - fielding , _
superskyflat_-fielding , and correction of wcs header information ) onto the tangent plane of a common reference image using _
mscimage_. the reprojection performed by _ mscimage _ should be performed with flux conservation set to _ no _ the first time but then set to _ yes _ for any further reprojections .
this is because the pixel scales in raw mosaic ii images are a function of radius from the field center with several percent variation from center to corner of the field .
the process of flat - fielding removes the true illumination pattern caused by more photons falling in larger pixels leading to an illusion of flat background counts despite the variation in pixel scales .
mscimage _ is used for astrometric projection into a uniformly sized grid of pixels , flux conservation would re - introduce the illumination pattern , thereby causing a variable photometric zeropoint across the image . turning off
flux conservation instead produces an average of the values of initial pixels neighboring each final pixel location that offsets the error introduced by flat - fielding .
we make three different stacked `` final '' images for each filter .
unweighted versions are made first in order to determine the weights for the other two versions .
the versions denoted by _
xs _ were combined using weights optimized for surface brightness , which is the method used by the ndwfs cookbook .
the versions denoted by _
ps _ were made using weights optimized for point sources ; details can be found in the appendix . we synthesize a @xmath52__ps _
image to use as our photometric detection image by adding the @xmath21 , @xmath22 , and @xmath4 _ ps _ stacks using point - source - optimized weights derived from the signal and noise statistics of these stacks .
we decided not to add @xmath18 to @xmath52 although it is also very deep because some low - redshift objects have different morphologies in @xmath18 , the seeing is typically worse in @xmath18 , and objects at @xmath5 will be nearly invisible in @xmath18 so adding it to @xmath52 would just add noise for these objects .
note that for similar reasons objects that drop out in @xmath21 and @xmath22 are somewhat less likely to be detected in our @xmath52 image than in a @xmath53 or @xmath4 image alone .
we prefer the @xmath52 approach to the creation of a @xmath54 image advocated by @xcite because it allows us to measure object morphological parameters directly from the same image used for object detection , which is not recommended in the latter approach .
moreover , the optimally - weighted @xmath52 combination reduces the emphasis on single - filter outliers employed by the @xmath54 approach and avoids its flaw of treating both negative and positive sky fluctuations as evidence of an object .
we trim the final @xmath52__ps _ image to the maximum size region that has nearly uniform signal - to - noise ratio in all three input images @xmath21 , @xmath22 , and @xmath4 . for e - hdfs ,
the trimmed image is 7395@xmath487749 pixels or @xmath55 .
then the iraf task _ imalign _ is used to shift and trim the other images to match , giving our images the common origin and size required for aperture photometry with sextractor .
the trimmed images are normalized to an effective exposure time of one second . after the pipeline reduction
, the background is flat to better than 1% in all filters .
the remaining low level , large scale fluctuations were subtracted using sextractor @xcite with mesh size back_size=64 and median filter back_filtersize set to 6 mesh units .
the resulting final @xmath52 image is shown in figure [ fig : bvr ] with a magnified image of a small section shown in fig .
[ fig : zoom ] .
new fits headers have been added to these final images indicating input information for sextractor photometry runs and other routines , specifically : + satur_level , the empirically determined saturation level in each image , which is usually a factor of a few less than the apparent saturation level of the brightest stars , + seeing_fwhm , the mode of the seeing for a set of bright unsaturated stars , + mag_zeropoint , the ab magnitude of an object that gives 1 count per second derived from photometric calibration as described below , + flux_zeropoint , the flux in @xmath40jy corresponding to 1 count per second , + tot_exptime , the total exposure time that each final image represents , and + gain_tot , the number of photoelectrons represented by each count in the final image .
our calibration scheme is to take @xmath565 minute exposures of our fields ( referred to as `` calibration images '' ) in all optical filters on photometric nights .
this was achieved for e - hdfs on oct .
photometric calibration is performed using a version of the landolt catalog with all magnitudes and colors converted into the ab95 system of @xcite ( referred to throughout this paper as ab magnitudes ) .
we included @xmath46 magnitudes for landolt standard stars in our catalog of standards using the formula @xmath57 which was generated from the formulae of @xcite and has been corrected for a 0.02 magnitude bias found when comparing this prediction with measurements of @xmath46 magnitudes of bright landolt stars by @xcite .
the resulting 0.04 magnitude rms error was used to predict errors in these estimates .
this procedure offers a significant reduction in observing time and a tremendous increase in the number of calibration stars versus the traditional method of using spectrophotometric standard stars for filters outside the johnson - cousins system .
only a few stars per landolt field have been turned into sdss calibrators @xcite and they are typically too bright for the 4 m telescope to use without defocusing .
table [ tab : photcal ] lists the airmass and seeing for each of our calibration images along with our photometric solution for the magntiude in terms of measured counts per second , @xmath58 where @xmath59 is the magnitude zeropoint , @xmath60 is the airmass coefficient , @xmath61 is airmass , and the color term is listed in table [ tab : photcal ] .
our images of landolt standard fields occupied a range of airmasses bracketing the airmasses of these calibration images but insufficient to break the degeneracy between zeropoint and airmass coefficient .
we therefore used fixed airmass coefficients , with the resulting uncertainty reflected in the uncertainty in the zeropoint fit , which is @xmath621.5% in each filter . since ab95 was carefully calibrated for the standard johnson - cousins @xmath63 filter set , the conversion of the landolt catalog is quite accurate .
however , the filter system used at ctio is not a precise match to johnson - cousins , making the use of a color term helpful in photometric calibration .
the relatively small color coefficients shown in table [ tab : photcal ] minimize the scatter in our photometric solution by providing a better estimate of the true ab magnitude of each landolt star in our observed filter system .
we set the ab color coefficient to zero when determining photometry for our final object catalog ; this places the fluxes on our observed filter system rather than standard johnson - cousins filters .
a color correction can be calculated later for applications that require fluxes extrapolated to johnson - cousins colors . because photometric redshift codes multiply spectral templates by the atmospheric transmission , filter transmission , and ccd quantum efficiency used for the observations , they already account for the origin of the color term and
should be given fluxes in the observed filter system .
the presence of non - power - law features such as lyman and balmer breaks in these template spectra ( and real objects ) makes it desirable to avoid color transformations calibrated with stars .
we then used the calibration images to create standard stars in each filter for e - hdfs and used these new standard stars to calculate the zeropoint in our final images .
this avoids attempting to explicitly correct the photometry for the airmass and cirrus extinction in each image which contributed to our final images , allowing us to use images taken in non - photometric conditions as part of the final stack . for e - hdfs
we used a set of 10 stars , each from a different original amplifier , chosen to be bright enough to provide good signal in @xmath18 and yet to avoid saturation in @xmath4 , @xmath64 , and @xmath46 .
care was taken to choose some stars near the field center and some near the sides and corners of the field to check for systematic effects .
the photometry of the new standards was determined using the iraf routine _
phot_. we used the same @xmath65 diameter aperture for photometry of these standard stars in order to minimize aperture losses .
we repeated this process by using _
phot _ on our final science images for the same stars .
we checked for systematic patterns in the offsets but found none , and the rms scatter was 0.03 magnitudes , implying a smaller error in the mean offset .
the total exposure time , seeing , photometric zeropoints , and 5 @xmath66 point source detection depths ( ab ) of our final images are given in table [ tab : final ] .
flux zeropoints are also determined using the definition @xcite @xmath67 given the measured uncertainties in the initial zeropoints and our measured scatter in the solution for the zeropoints of the final images , we estimate a total systematic uncertainty of 3% in the zeropoint of each filter .
derivations of optimal photometric techniques and their uncertainties typically assume poisson noise which is uncorrelated between image pixels . however , the astrometric reprojection used to place all images on a common grid correlates neighboring pixels , leaving the large - scale noise properties unchanged but making a pixel - by - pixel rms a poor measure of the typical noise due to sky fluctuations in a larger aperture . as this is the method used by sextractor to estimate photometric errors
, we undertook a detailed investigation to test its accuracy , which turns out to be poor , and to find an improved method .
our empirical method for estimating photometric uncertainties as a function of aperture size is also sensitive to large - scale noise correlations from unsubtracted nearby objects and from any ccd defects that survive flat - fielding .
wings from bright objects appear to affect the noise statistics ; when we set all pixels belonging to detected objects to zero and displayed the resultant image , wings are clearly visible around the brightest objects .
we estimate the error due to background fluctuations and noise correlations in each filter via a custom idl code which places @xmath68 random apertures of a given size on the sky - subtracted image which do not overlap with any of the pixels in the segmentation map of isophotal object regions produced by sextractor .
we use circular apertures of area @xmath69 centered at integer pixels and describe them by an effective size @xmath70 .
a gaussian is then fit to the histogram of aperture fluxes to yield the rms background fluctuation as shown in fig .
[ fig : hist ] .
the histograms appear well described by the best - fit gaussians .
figure [ fig : sigmav ] shows the rms of aperture fluxes vs. n for raw , sky - subtracted raw , and final @xmath22 images .
the solid curve gives our fit to the noise properties of the final @xmath22 image using the function @xmath71 where @xmath72 , @xmath73 for the measured rms pixel noise of @xmath74=0.014 .
an alternative formula suggested by @xcite , @xmath75 yields an equally good fit , with @xmath76 , @xmath77 .
this formula shows an explicit sum between contributions from poisson noise , which is independent from one pixel to another and hence has rms proportional to @xmath78 , and correlated noise from fluctuations in background level on scales larger than the aperture which yields an rms proportional to @xmath79 .
these simple , extreme cases are shown as dotted lines starting from the value of @xmath80 in figure [ fig : sigmav ] and are seen to bracket the true behavior except at very small aperture sizes .
apertures of just a few pixels in area are affected by the small - scale noise correlations introduced by re - pixelization performed during re - projection ; this effect should actually depress @xmath80 since each final pixel is the average over several nominally independent input pixels .
we prefer the simplicity of eq .
[ eq : power ] which reflects the reality that noise exists on a range of scales leading to an effective power - law behavior intermediate between these two extremes .
table [ tab : noise ] lists our measured rms - per - pixel @xmath74 and the fit coefficients @xmath27 and @xmath81 for all of our filters ; the uncertainties in @xmath27 and @xmath81 are highly correlated but both have been determined to 5% precision , whereas the uncertainties in @xmath74 measurements are at the 1% level .
a similar analysis of our raw images shown in figure [ fig : sigmav ] reveals that the noise in the raw images ( squares ) is fully correlated due to the domination of background fluctuations . once a sky - subtraction has been performed ( diamonds ) ,
the noise properties resemble those of the final image except for being noisier by a factor roughly equal to the square root of the number of exposures .
while some of this improvement comes from flat - fielding , it appears that the correlated noise is due to background fluctuations from objects whose subtraction improves with the square root of observing time .
this is inconsistent with confusion noise from undetected sources , which is not expected to be a significant factor even in optical images this deep .
one possibility for the origin of this correlated noise is that the global sky subtraction performed by sextractor is insufficient .
however , when we mimic the backphoto_type local mode of sextractor by estimating the local sky in a square annulus around each aperture and subtracting this off we find a slight increase in the noise correlations .
hence it appears that the background contributions from sky gradients and nearby objects have been subtracted about as well as one could expect , and that the remaining contributions increase the noise significantly above the extrapolation from @xmath74 one would assume for the case of uncorrelated noise .
we wish to determine accurate photometry for a wide range of object sizes and brightnesses , with particular attention to high - redshift galaxies with half - light radii 0.2 - 0.3@xmath49 @xcite ) which are dim and typically unresolved in our @xmath82 seeing and can therefore be approximated as point sources .
indeed , @xcite showed that most galaxies at @xmath83 have half - light radii @xmath84 regardless of redshift . in traditional aperture photometry ,
the flux of an object is measured by giving constant weight to one region of the image and zero weight to the rest , with fractional weights used for pixels which fall only partially in the chosen aperture .
the simplest case is a circular aperture centered at the object centroid . for the case of uncorrelated noise dominated by the sky background and a gaussian psf ,
the circular aperture with optimal signal - to - noise is easy to derive and has a diameter equal to 1.35 times the seeing fwhm .
however , the actual psf in each filter has a near - gaussian core with broad wings i.e. more of the flux is found outside the fwhm than for a gaussian . to investigate the effect of the non - gaussian psf on photometry we measured the enclosed flux as a function of radius for a set of objects selected to be bright but unsaturated ( @xmath85 ) and unresolved according to sextractor ( @xmath86 , as discussed below ) .
[ fig : psf ] shows the results for the @xmath18 , @xmath22 , and @xmath64 bands plus a gaussian of 0.96@xmath49 fwhm that represents an excellent fit to the core of the @xmath22-band psf but only encloses 80% of the total flux .
a 14@xmath49 diameter aperture was used to define the total object flux , with rapid convergence seen in the bluer filters but only gradually in @xmath4,@xmath64 , and @xmath46 .
while this could be due in part to contributions from neighboring objects , we used the median statistic to define these profiles to avoid sensitivity to the small fraction of objects with bright neighbors .
the contribution of nearby objects should lead to a divergence with increasing aperture diameter which is not seen , so it appears that the psf wings do contain a higher fraction of the total flux in the redder bands , with perhaps a few percent missed by our 14@xmath49 diameter aperture .
figure [ fig : sn ] plots the enclosed flux , the rms background fluctuation , and the resulting signal - to - noise for point sources as a function of aperture diameter .
there are competing effects which serve to make the optimal aperture slightly smaller than the idealized result of 1.35 @xmath48 fwhm .
the non - gaussian psf places more flux outside the `` optimal '' diameter , so a larger aperture would find more signal .
however , the large - scale noise correlations mean that noise increases with aperture size faster than in the ideal poisson case .
table [ tab : apertures ] lists various seeing measurements along with the optimal circular aperture diameters and fractional flux enclosed within these apertures for each filter which is used to correct the observed fluxes to estimate the total object flux in each filter .
a systematic increase in the seeing occurs when the estimate is made using sextractor s _ fwhm _ parameter , which assumes a gaussian psf , instead of iraf s _ imexam _ routine , which allows a moffat fit .
the non - gaussian psf shape is reflected in the further increase in the diameter containing half of the object light implied by the half - light radii @xmath87 output by sextractor ; a gaussian contains half of the flux within one fwhm .
the optimal circular apertures for point sources turn out to be quite close to these values for the half - light radii .
the image quality is measured to be nearly constant across the @xmath88 field , with systematic variations in fwhm and half - light radius seen at the 10% level in the @xmath64 band but at no more than the 5% level in the other filters .    for extended sources ,
our apertures are no longer optimal , and the corrections for the fraction of enclosed flux are incorrect .
these sources are barely affected by seeing smaller than their intrinsic size , leading to a flux better described by the aperture area than by the fractional enclosed point - source flux .
the traditional solution is to use larger apertures ( @xmath89 diameter ) to measure colors of extended objects as well as point sources .
our analysis shows that the signal - to - noise of point source photometry is reduced by up to 30% by using these larger apertures , with a bias in the @xmath90 color of point sources of 20% ( 0.2 magnitudes ) introduced by the different seeing in these images .
this could be solved by convolving all of the images to the relatively poor seeing ( fwhm=1.3@xmath49 ) of the @xmath18 band ( see @xcite ) .
however , the non - gaussian psf shape requires a non - gaussian convolution kernel ; simply matching the fwhm by gaussian convolution is not sufficient to fully match the psf .
moreover , smoothing the other images to match the psf shape of the @xmath18 band would be costly in terms of decreasing the signal - to - noise for point source photometry , especially given the observed large - scale correlations in background noise .
in choosing the `` optimal '' aperture sizes for each filter in table [ tab : apertures ] , we have chosen amongst the range of apertures with signal - to - noise@xmath9195% of the optimal value in a manner that reduces the range of aperture sizes in order to reduce the errors introduced for extended sources .
this has the effect of choosing slightly larger apertures in the redder filters than would formally maximize the signal - to - noise .
this provides additional benefits , as the approximations that high - redshift galaxies are point sources , image co - registration is perfect , and object positions have no error all become less valid as the apertures decrease in size .
our simulations indicate that these combined effects will not affect object colors beyond the 10% level for the chosen aperture sizes .
our primary solution to the challenge of obtaining decent colors for both unresolved and extended objects dimmer than the night sky is to create corrected aperture fluxes ( hereafter referred to as apcorr ) .
we tried using the a_image and b_image parameters to model the flux distribution as an ellipsoid but encountered too many systematic problems with these parameters to prefer this approach .
we therefore use the half - light radius ( flux_radius ) of each object determined by sextractor in the @xmath52 detection image and assume a two - dimensional gaussian light profile with twice this value as fwhm .
the half - light radius measured in the @xmath52 detection image represents the object s intrinsic size convolved with the @xmath52 seeing , and this allows us to predict the observed size in each filter given the known seeing measured from the median half - light radii of bright point sources listed in table [ tab : apertures ] . for object @xmath92 in filter
@xmath93 we expect @xmath94 \ ; \ ; \ ; ,    \label{eq : sigma_squared_ij}\ ] ] where the last term in parentheses modifies the quadrature sum of intrinsic object size and bvr seeing to the needed sum of intrinsic object size and seeing in filter @xmath93 and the prefactor converts from half - light radius ( which equals half - width half maximum for a gaussian ) to the rms of the light profile .
the fractional flux of object @xmath92 falling within the optimal circular aperture radius @xmath95 for filter @xmath93 listed in table [ tab : apertures ] is then assumed to be @xmath96 where the first value in parentheses refers to the fractional flux for point sources listed in table [ tab : apertures ] and the second gives the fractional flux of the two - dimensional gaussian which approximates the object .
taking the minimum of these two values prevents noisy data from causing an object to be inferred to be smaller than a point source .
each object s flux measured within the optimal aperture for a given filter is then corrected by dividing by @xmath97 to infer the total flux in that filter .
another solution for attempting to handle both extended and unresolved sources is sextractor s auto aperture , which uses a kron ellipse adapted to each object s light profile to find @xmath98 of the flux for extended sources and @xmath99 for point sources using the parameters we have adopted , according to @xcite . in the @xmath18 , @xmath21 , and @xmath22 filters , we find that auto indeed measures 96% of the point source flux found in a 14@xmath49 diameter aperture , with the value decreasing to 93% in @xmath4 , 92% in @xmath64 , and 89% in @xmath46 in keeping with the broader psf discussed above .      because sextractor assumes poisson sky noise , it produces underestimated photometric uncertainties for objects dimmer than the background sky .
detection image as a source of uncertainty in object detection and the background in the measurement image as a source of uncertainty in the photometry . ]
the standard formula used by sextractor to predict photometric errors in flux as a function of aperture size is @xmath100 where the first term represents sky fluctuations assuming uncorrelated noise between pixels , @xmath101 is the object flux in counts within the aperture , and @xmath102 is the total effective gain used to turn counts in the normalized image into original photons detected by the ccd , with the fraction @xmath103 giving the poisson variance .
sextractor measures @xmath74 using the median fluctuation in its internally computed background map , and this measurement agrees with ours to within a few percent .
we modified this to make use of our fits for ( correlated ) background noise which yields @xmath104 where the first term now represents empirical sky fluctuations which include the subdominant contribution from electronic readout noise .
both noise terms will be larger in regions of the image where reduced exposure time leads to a lower effective gain and hence higher background rms .
sextractor accounts for this by using a variance map @xmath105 which is either input as a weight_image or generated internally from a map of background fluctuations in the image .
the errors produced by sextractor as a function of position are then @xmath106 the ratio of the formulas in equations [ eq : uncertainty_se ] and [ eq : uncertainty ] allows us to make a simple correction to the flux uncertainties output by sextractor to yield @xmath107 for negative fluxes , @xmath101 is set to zero in all of these error terms , since the uncertainty is entirely dominated by the background fluctuations . for auto apertures ,
the area is given by that of the kron ellipse used i.e. @xmath108 ( the so - called kron_radius is really a scaling of the a_image and b_image parameters based on the second moment of the object light distribution ) .
the resulting rms has units of counts per second and is converted to flux using the zeropoints in table [ tab : final ] . the uncertainty determined for optimized aperture fluxes
is then multiplied by the same correction factor used to estimate the total object flux i.e. the reciprocal of @xmath97 defined in eq .
[ eq : fractional_flux_ij ] .
finally , the uncertainty in the corrected aperture flux is increased by a factor @xmath109 which serves to amplify the uncertainties for extended sources to account for the uncertainty in their correction factors .
this reproduces the errors found in corrected aperture fluxes for extended objects in our simulations described in  [ sec : simulations ] .
the photometric errors derived in this manner do not include the 3% calibration uncertainty which is common to all sources in a given filter .
we used the iraf package _ artdata _ to simulate stars and galaxies with known magnitudes and positions and to add them to our observed images to get realistic crowding effects and background noise .
we ran sextractor on these `` simulated '' images in dual - image mode to simulate our full photometric pipeline .
we found object centroiding errors to be @xmath110 pixels i.e. @xmath111 at @xmath112 . in the centroiding of individual objects we found auto photometry to be nearly unbiased for both single filter fluxes and colors .
however , figure [ fig : stars ] shows that auto has larger errors for point sources than the corrected aperture fluxes due to the larger auto apertures including significantly more sky noise , so we recommend auto fluxes only for extended sources and give fluxes for both types of apertures in our catalog .
the errors seen in the @xmath18 and @xmath22 filters appear to be entirely uncorrelated .
figure [ fig : gals ] shows that auto performs similarly to corrected aperture fluxes for extended objects , although both show significant covariance between @xmath18 and @xmath22 due to misestimation of the true object light profile .
figure [ fig : sigma_stars ] shows the flux errors divided by their uncertainties for simulated point sources .
apcorr again appears unbiased but this plot reveals that the error estimates are roughly equally accurate for both cases .
the median squared value is close to one in both cases , so the reported flux uncertainties in the catalog appear trustworthy for point sources .
figure [ fig : sigma_gal ] shows the flux errors divided by their nominal uncertainties for simulated galaxies .
this plot illustrates that the estimated uncertainties are typically too small for galaxies in both cases .
figure [ fig : colors ] shows that @xmath90 colors of simulated objects appear slightly biased for auto but unbiased for apcorr , and the errors are a bit smaller for apcorr for point sources and of similar size for both types of fluxes for galaxies .
the color errors are highly correlated , especially for galaxies , implying that these errors are primarily caused by errors in the isophotal object detection by sextractor on which both fluxes depend .
these include problems caused by blending with nearby neighbors , although the median object is unaffected by neighbors in this uncrowded field .
the auto and apcorr colors of real objects are discussed in  [ sec : lbg_select ] .
we performed aperture photometry using sextractor in dual - image mode with @xmath52 as the detection image and each final image as the measurement image , leading to a common catalog for all filters .
detection was performed after filtering with a @xmath113 pixel gaussian convolution kernel with 4 pixel ( 1.07@xmath49 ) fwhm which is well - matched to the seeing in our @xmath52 image ( see appendix ) .
we used a threshold of 1.2 for both detection and analysis ; a single pixel at the 5 @xmath66 level in the original @xmath52 image would barely make the 1.2@xmath66 threshold after filtering . using the identical parameters
, we detected 318 objects in a negative of the @xmath52 image , so given the symmetrical nature of sky fluctuations and readnoise we estimate an equivalent number of spurious objects in our final catalog i.e. 0.5% .
for each filter , we used the `` optimal '' apertures and corrections described above as well as the auto aperture determined for each object by sextractor centered on the @xmath52 barycenter .
the photometry was corrected for galactic dust extinction of e(b - v)=0.03 @xcite .
the 5@xmath66 point source detection depths in table [ tab : apertures ] were determined by multiplying the @xmath66 of the measured sky noise for the optimal point source aperture sizes by 5 and then correcting for the flux missed by these apertures to turn sky fluctuations into total point source magnitudes .
many estimates of point source detection limits in the literature are based on extrapolating the rms pixel noise , @xmath74 , to the chosen aperture size assuming uncorrelated sky noise . for correlated noise like that in our images , which we expect
is typical , this will overestimate the true depth in a 2@xmath49 diameter aperture by 0.4 magnitudes and in a 3@xmath49 diameter aperture by 0.6 magnitudes .
our optical photometric catalog of 62968 objects in e - hdfs is available in full in the electronic version of the journal and from our website , with the first several lines shown in table [ tab : catalog ] .
all positions and shape parameters are measured from the final @xmath52 image .
the columns in the catalog offer the following information : + _ column 1 _ : object number starting with first object as 0 + _ column 2 _ : object name ( musyc- ) + _ column 3 _ : right ascension in decimal hours , double precision + _ column 4 _ : declination in decimal degrees , double precision + _ column 5 _ : x barycenter + _ column 6 _ : y barycenter + _ column 7 _ : stellarity classification + _ column 8 _ : half - light radius + _ column 9 _ : rms of flux distribution along major axis + _ column 10 _ : rms of flux distribution along minor axis + _ column 11 _ : object position angle counterclockwise from x - axis + _ column 12 _ : flags output by sextractor + _ column 13,15,17,19,21,23 _ : auto flux in @xmath18,@xmath21,@xmath22,@xmath4,@xmath64,@xmath46 + _ column 14,16,18,20,22,24 _ : auto flux uncertainty in @xmath18,@xmath21,@xmath22,@xmath4,@xmath64,@xmath46 from ( [ eq : correction ] ) + _ column 25,27,29,31,33,35 _ : apcorr flux in @xmath18,@xmath21,@xmath22,@xmath4,@xmath64,@xmath46 after division by correction factor based on half - light radius described in eq .
[ eq : fractional_flux_ij ] + _ column 26,28,30,32,34,36 _ : apcorr flux uncertainty in @xmath18,@xmath21,@xmath22,@xmath4,@xmath64,@xmath46 from ( [ eq : correction ] ) after division by correction factors based on half - light radius described in eq .
[ eq : fractional_flux_ij ] and multiplication by additional factor for extended sources described above + photometric measures in the catalog are given in units of flux density in @xmath40jy ( 1 @xmath40jy = @xmath114 w m@xmath7 hz@xmath115 ) .
the conversion to ab magnitudes is simple using the formula @xmath116 @xcite .
this avoids the loss of information that comes from sextractor representing objects with negative aperture fluxes as having @xmath117 and the confusion that results from flux errors consistent with a non - detection being turned into extremely large magnitudes .
photometric errors are nearly symmetrical in flux but not in magnitudes , even for objects with low signal - to - noise , making it acceptable to represent the uncertainties with a single number .
this makes flux density the natural unit to be used for calculating photometric redshifts if one wishes to use a @xmath54 function to compute the likelihood as is done by the publicly available codes .
for color selection , one can either turn color criteria into desired flux ratios between filters or turn the measured fluxes into magnitudes and subtract to obtain colors .
sextractor offers a neural network classification of objects into `` stars '' and galaxies ; @xmath86 implies an unresolved profile consistent with the psf , @xmath118 implies an extended profile , and @xmath119 represents uncertain objects which are typically too dim for this classification to work successfully .
figure [ fig : class ] shows the stellarity results for our catalog .
the classification starts to break down at @xmath120 , with brighter objects confidently split between stellar and extended profiles but many dimmer objects receiving an inconclusive value .
stellarity becomes nearly useless at @xmath121 where the majority of values fall in the uncertain regime , which is appropriate given that @xmath122 for most galaxies at @xmath123 @xcite .
however , the vast majority of objects at @xmath120 are galaxies due to the rapidly rising galaxy luminosity function and the plateau in the stellar equivalent .
this means that for many scientific goals all objects with @xmath124 may be considered galaxies .
figure [ fig : radec ] shows the locations of all objects in our catalog .
there is an overabundance of objects along the image border which is not apparent to the eye but can be eliminated by removing the 574 objects with @xmath125 , which identifies objects truncated by the image border .
the one visual blemish is a line of spurious objects detected along a heavily saturated column extending from the brightest star in the field .
this line is removed by eliminating all 1226 objects with @xmath126 , which implies saturation in at least one filter .
the sky density of objects in our catalog is greatly reduced near the brightest stars , leading us to create a mask to properly exclude these regions from analyses of sky density and clustering .
our methodology is described in detail in l.  infante et al .
( 2005 , in preparation ) . a careful analysis of the implied stellar and galaxy luminosity functions as a function of stellarity cut shows that stellar contamination is minimized with negligible loss of galaxy counts by requiring @xmath124 .
the total @xmath4-band magnitude counts ( from auto ) are shown in figure [ fig : rcounts ] and appear 90% complete to @xmath127 and 50% complete to @xmath128 .
part of the incompleteness comes from auto beginning to underestimate the flux of very dim objects , as noted by @xcite .
the fit for the number of galaxies per magnitude per square degree is @xmath129 .
our slope of 0.34@xmath1300.01 agrees well with previous measurements of 0.321@xmath1300.001 by @xcite , 0.361@xmath1300.004 by @xcite , 0.31@xmath1300.01 and 0.34@xmath1300.01 in two different fields by @xcite , and 0.39 by @xcite .
our differential counts at @xmath128 of 1.3@xmath131 agree well with the values of 1.2@xmath132 reported by @xcite .
-band galaxy counts per unit magnitude per square degree with poisson error bars and solid fit of @xmath133 based on counts at total ab magnitudes @xmath134 .
incompleteness is minimal below total magnitude @xmath127 and reaches 50% around @xmath128 .
[ fig : rcounts ] ]      various filter sets have yielded success at finding lbgs , including @xmath135 @xcite , sdss @xmath136 @xcite , @xmath137 @xcite , @xmath138 @xcite , and @xmath139 ( @xcite , @xcite ) .
figure [ fig : uvr_stars ] shows star colors expected and measured in our survey versus our adopted @xmath6 color selection criteria of @xmath140 where all magnitudes are ab . the color criteria shown as dashed lines in the figure result from the transformation of the @xmath141 @xmath142 @xmath143 criteria described by @xcite into @xmath6 . and @xmath144 . ]
although the upper branch of the stellar distribution at red @xmath145 colors represents giants which are rare at our survey depth and high galactic latitude , we shifted our criteria to avoid this region as it also contains dim dwarf stars with correspondingly large errors in color , which are the primary expected source of interlopers .
the region of the steidel et al .
criteria avoided by this shift had the highest fraction of interlopers and also shows a high interloper fraction in the more liberal @xmath6-selected subset of the @xcite lbg sample .
the lbg template curve falls within our selection region for @xmath146 . at higher redshifts ,
the lbgs rapidly become too red in @xmath145 to distinguish from lower - redshift objects and are better selected as dropouts in @xmath147 .
we allow the color selection region to extend far to the blue in @xmath145 to account for the effect of lyman @xmath27 emission lines , which add flux in @xmath22 .    for spectroscopic selection and clustering analysis , objects with @xmath18 fluxes , @xmath148 , less than their @xmath149 flux uncertainty , @xmath150 ,
were assigned an upper limit of @xmath151 to make it less likely that negative sky fluctuations were responsible for turning a dim object in @xmath18 into a formal dropout in @xmath90 .
this correction avoids lower - redshift interlopers at the cost of some incompleteness .
the vast majority of our lbg candidates are formal dropouts in the @xmath18 filter , with @xmath152 before being set to these limit values .
the limit values were used to generate the @xmath90 colors plotted in fig . [
fig : uvr ] , which shows our full catalog including @xmath35 lyman break galaxy candidates selected by their @xmath6 colors .
our selection criteria yield 1607 candidates in 1137 square arcminutes , or 1.4 arcmin@xmath7 .
@xcite found 1.2 arcmin@xmath7 at @xmath153 using stricter criteria than ours , @xcite found 1.7 arcmin@xmath7 using the extended @xmath141 @xmath142 @xmath143 criteria shown translated into @xmath6 in figure [ fig : uvr_stars ] , @xcite found 1.5 arcmin@xmath7 using ubr , and @xcite found 1.5 arcmin@xmath7 using @xmath154 criteria .
@xcite found a redshift distribution @xmath155 .
our @xmath6 criteria are expected to start selecting lbgs at redshifts higher by 0.16 because the red cutoff of the @xmath18-band filter and its effective wavelength after accounting for the mosaic ii ccd qe are both @xmath156  redder than for @xmath141 .
our criteria are expected to stop selecting lbgs at redshifts higher by up to 0.4 because the lyman alpha forest causes red @xmath145 colors at a higher redshift than in @xmath157 due to the longer effective wavelength of @xmath22 .
hence we expect @xmath158 , and this will be calibrated via spectroscopy .
figure [ fig : lbg_radec ] shows the spatial distribution of these lbg candidates .
we estimate the angular correlation function of these objects using the @xcite estimator @xmath159 where @xmath160 represents the number of pairs at separation @xmath161 in the object catalog , @xmath162 represents the appropriately - normalized number of pairs at separation @xmath161 in a poisson - distributed mock catalog , and @xmath163 represents the appropriately - normalized number of pairs at separation @xmath161 with one in the object catalog and the other in the mock catalog .
this estimator has variance @xmath164 \simeq { \frac{(1 + w(\theta))^2}{rr } } \ ; \ ; \ ; ,   \label{eq : ls_var}\ ] ] where the denominator gives the poisson error contribution i.e. the number of pairs expected to be found in the bin centered on @xmath161 given the survey geometry and the total number of lbgs in the survey . from @xcite ] the observed values of @xmath165 must then be corrected for the integral constraint caused by estimating the true sky density of lbgs from our survey to find a power - law fit @xcite , @xmath166 in practice , this correction is performed by using a modified estimator @xmath167 where the correction factor @xmath168 is given by @xmath169 a is estimated iteratively and usually converges in a few iterations .
the angular correlation function found for lyman break galaxy candidates in figure [ fig : w_lbg ] is @xmath170 , where @xmath161 is measured in arcsec , intermediate between previous results for fixed @xmath171 from @xcite of @xmath172 and @xcite of @xmath173 .
this implies a correlation length for our sample of @xmath174@xmath175mpc and an lbg bias factor of @xmath176 for the @xmath11cdm cosmology .
these results will be updated after spectroscopic calibration of our @xmath6 lbg selection and by including our full square degree survey .     ]
musyc is unique for its combination of depth and total area , for the coverage in x - ray , uv , mid- and far - infrared wavelengths and for providing the @xmath0 plus near - infrared photometry needed to produce high - quality photometric redshifts over a square - degree of sky .
this multiwavelength coverage will enable comparison of selection effects which have previously complicated the study of galaxies in the high - redshift universe .
the large area ensures good statistics for studies of clustering , luminosity functions , and surface densities .
we expect that the careful attention paid to optimized photometry methods and the estimation of accurate photometric uncertainties will allow our public data to make a significant impact for a wide range of scientific studies .
in particular , this high - quality @xmath0 data , combined with forthcoming near - infrared photometry , should yield photometric redshifts with good accuracy .
it is interesting to note that some of the standard techniques perform well despite making a number of idealized assumptions .
for instance , we found that the optimal point source apertures are remarkably close in diameter to the idealized case of 1.35 fwhm due to a partial cancellation of the non - gaussian psf preferring a larger aperture and the correlated noise preferring a smaller aperture .
the more traditional approach of using @xmath177 apertures produces biases in point - source colors . alleviating those biases by `` psf - matching '' i.e. smoothing the images to match the fwhm of the image with the worst seeing
must be done carefully given the non - gaussian psf shapes and will result in degraded signal - to - noise for point sources .
estimates of point source detection depth that assume uncorrelated noise and extrapolate from the rms pixel noise will overestimate the depth by about 0.5 magnitudes for 2@xmath49 or 3@xmath49 apertures .
the apcorr fluxes introduced in this paper perform extremely well , as seen in fig .
[ fig : uvr ] .
the nearly unbiased performance of apcorr fluxes for both point sources and extended objects is impressive , and the correlated errors between filters for extended objects are acceptable if the flux is misestimated by a constant fraction in all filters .
however , one should use caution when calculating photometric redshifts on extended objects in photometric catalogs generated using apertures which might miss a different fraction of the object flux in each filter .
the best results should be obtained by using the apcorr fluxes in our catalog for unresolved and barely resolved sources , with further investigation needed to determine whether apcorr outperforms auto for extended sources .
our @xmath0 images and catalogs of the @xmath88 extended hubble deep field south are available to the public at http://www.astro.yale.edu/musyc .
we thank buell januzzi and frank valdes for their detailed answers and trouble - shooting with questions on the mscred package and benne w. holwerda whose guide , `` source extractor for dummies , '' provided a valuable supplement to the official sextractor user s manual written by e. bertin .
we acknowledge helpful conversations with anton koekemoer , david spergel , and nick suntzeff .
we thank the additional members of the musyc collaboration , including danilo marchesini , ken rines , ned taylor , and bill van altena for their valuable input .
we are grateful for support from fundacin andes , the fondap centro de astrofsica , and the yale perspectives in science program .
this material is based upon work supported by the national science foundation under grant .
ast-0201667 , an nsf astronomy and astrophysics postdoctoral fellowship ( aapf ) awarded to e.g. s.t . acknowledges the support of the danish natural research council .
we thank the staff of cerro tololo inter - american observatory for their invaluable assistance with our observations .
this research has made use of nasa s astrophysics data system .
we wish to combine a series of processed individual images which have signal @xmath178 and noise @xmath179 .
the exact signal obviously depends on which star is measured , but if a common set of stars are considered the variations in signal will be due to differences in exposure time , cirrus extinction , and atmospheric extinction .
the noise should be sky - dominated away from bright objects but measuring the rms in `` blank '' regions of sky properly accounts for the readnoise as well . for surface - brightness - optimization
we want to average the total signal from a set of stars , as measured by the iraf routine _ phot _ or the inverse of the _ mscscale _ header output by _
mscimatch_. for point - source - optimization we want to determine the average signal that falls within our eventual optimal aperture , which will have a diameter of roughly 1.35 times the fwhm of the final image .
the seeing in the final image can be estimated or measured from an initially unweighted stack , or an iterative procedure can be used .
alternatively , a correction to the total star counts @xmath180 can be calculated by integrating an assumed gaussian psf out to the diameter of the expected optimal aperture of radius r ; if the seeing in each image is @xmath181 fwhm@xmath182 , we obtain @xmath183 \ ; \ ; \ ; . \label{eq : s_i}\ ] ] in either type of optimization , the noise value used should be the background rms for the expected aperture size , which depends on the level of noise correlation between pixels . assuming that the noise correlations have the same behavior in each individual image , it is sufficient to use the pixel - by - pixel rms i.e. @xmath74 measured in each image since this will produce the right relative weights .
one can derive the weights that optimize the @xmath184 of the combined image by noting that the combined image will have @xmath185 and @xmath186 ( the weighted sum is just an unweighted sum of new images having signal @xmath187 and noise @xmath188 ) . maximizing the function @xmath184 produces the result @xmath189 , where multiplying all the weights by a constant @xmath190 preserves the final @xmath184 . in astronomy , it is typical to first scale the images to have equal signal levels ; this accounts for differences in exposure time and extinction to provide constant photometry and makes it easier to delete outlying pixels due to cosmic rays , satellite trails , etc .
the iraf-_mscred _ routine _ mscstack _ performs this scaling by multiplying by the _ mscscale _ header value and then uses a weightfile to weight the * scaled * images .
so the weights used must be for the scaled images , which leads to the formula for surface - brightness - optimized weighting of @xmath191 which is equivalent to @xmath192 in the original image because the image being weighted is different than the one for which @xmath178 and @xmath179 were measured .
the output image is the same since instead of multiplying @xmath178 by @xmath193 and adding we are equating the signal levels and then multiplying 1 by @xmath192 , but performing the scaling first makes it easier to remove outlying pixels from this sum as mentioned above .
the weights being used correspond to @xmath194 for the scaled images ; they reduce to the familiar case of inverse - variance weighting since the signal levels have already been equalized . for point - source optimized
weighting we use @xmath195 where @xmath196 is given by @xmath197 where we have re - expressed the formula for @xmath178 from eq .
[ eq : s_i ] in terms of the measured seeing in an unweighted ( or surface - brightness optimized ) stack , @xmath198 , and that in each individual image @xmath199 .
again , these weights are only valid when used on the post - scaled images , and they differ by a factor of @xmath178 from the weights appropriate for the original images .    due to the large number of individual exposures taken in e.g. @xmath18-band and multiple observing runs , it is sometimes necessary to perform stacking as an iterative procedure . to do this , the signal , noise , and seeing in the intermediate stacked images should be measured empirically as before . then _ mscstack _ or the equivalent can be used to stack these intermediate images just as if they were individual images , scaling by the new _ mscscale _ values and calculating a new weights file using the above equations .
as long as the final seeing is estimated well , no loss of signal - to - noise should occur from the iteration .
a corollary question is when to discard an image with poor @xmath184 rather than to include it in the stack .
the formal answer appears to be never as @xmath200 for all @xmath92 .
the magic of the optimal weighting formula can be seen : @xmath201 hence the optimal weights cause @xmath184 to add in quadrature so it will never formally decrease no matter how poor an input image is .
note that having a magnitude of extinction or twice the seeing still leaves an image with significant @xmath184 ; a @xmath184 threshold of 40% would require discarding images with about 2 times worse seeing or 1 magnitude of extinction ( if the sky background is unaffected ) .
the @xmath184 of the combination of the image with a given @xmath184 and an image with 40% of that @xmath184 is 8% higher than that of the better image alone so this is somewhat useful .
a conservative approach would be to cut entirely any images with @xmath202% that of the median image as a way of reducing systematic effects not accounted for by these idealized formulae .
@xcite showed using simulations that the best performance for point source detection was achieved by filtering the image with the psf itself .
this result is derived in the sextractor manual using fourier transforms . in practice due to the constraints of computation time one
chooses to cut off the psf at a finite radius , as is standard in the convolution kernels offered as part of the sextractor package .
the general case for an optimized weighted sum of pixel values @xmath203 , where each pixel has signal @xmath204 and noise @xmath205 , can be derived in analogy to the result for adding images found above , yielding @xmath206 .
the constant background noise for sky - dominated objects thus leads to a filter of the precise shape of the psf . for objects significantly brighter than the sky , @xmath207 giving constant weight i.e. a tophat which gradually morphs into the roughly gaussian shape of the psf as you move far enough away from the object center for the sky to dominate the noise .
sextractor only allows for a fixed convolution kernel so we optimized the detection for sky - dominated objects by using the psf ( truncated to 7@xmath487 pixels ) as our filter . by the same argument
, the optimum signal - to - noise measurement of the photometry of a point source will be obtained by a weighted sum of pixel fluxes weighted using the psf shape itself centered at the barycenter of the object .
this is typically referred to as `` psf photometry '' but it is not a supported feature of sextractor photometry , just of object detection .
note that this is no longer strictly optimal when the noise is correlated between pixels as we have found in our images , since the derivation of the optimal weights for each pixel assumed uncorrelated noise .
we decided against using psf photometry for this reason , and also because many of our science objects are slightly resolved and this would overweight the fluxes in their core , increase the risk of color biases , and complicate the task of correcting the aperture fluxes of extended objects .
llrllccc e - hdfs & 22:32:35.6 & -60:47:12 & ( 328,-49 ) & ( 311,-47)&0.03 & 1.37 & 1.6e+20 + e - cdfs & 03:32:29.0 & -27:48:47 & ( 224,-54)&(41,-45)&0.01 & 0.40 & 9.0e+19 + sdss1030 + 05 & 10:30:27.1&05:24:55 & ( 239,50)&(157,-4)&0.02&1.01&2.3e+20 + cw1255 + 01 & 12:55:40&01:07:00 & ( 306,64 ) & ( 192,7 ) & 0.02 & 0.81 & 1.6e+20 +    lcrrc 2002 oct 6,8 & @xmath18 & 47 & 28200 & 1.40 + 2002 oct 6,8,10 & @xmath21 & 13 & 7500 & 1.35 + 2002 oct 6,10,12 & @xmath22 & 30 & 10440 & 0.90 + 2002 oct 6,10,12 & @xmath4 & 21 & 6300 & 0.85 + 2002 oct 6,12 & @xmath64 & 26 & 6300 & 0.85 + 2002 oct 6,12 & @xmath46 & 11 & 2700 & 0.90 + 2003 may 26,27 & @xmath18 & 25 & 15000 & 1.20 + 2003 may 27 & @xmath21 & 10 & 6000 & 1.10 + 2003 may 27,28 & @xmath4 & 20 & 6000 & 1.10 + 2003 may 26 & @xmath46 & 15 & 3600 & 1.30 +    ccclcr u & 1.34 & 2.1 & 24.15@xmath1300.012 & 0.55 & ( -0.08@xmath1300.01)@xmath208 + b & 1.38 & 1.8 & 25.24@xmath1300.009 & 0.24 & ( -0.03@xmath1300.004)@xmath209 + v & 1.40 & 1.6 & 25.64@xmath1300.004 & 0.15 & ( 0.03@xmath1300.003)@xmath210 + r & 1.42 & 1.5 & 25.95@xmath1300.005 & 0.09 & ( -0.009@xmath1300.005)@xmath211 + i & 1.44 & 1.5 & 25.40@xmath1300.006 & 0.05 & ( 0.006@xmath1300.017)@xmath212 + @xmath46 & 1.46 & 1.5 & 24.61@xmath1300.014 & 0.03 & ( -0.06@xmath1300.054)@xmath212 +        ccccccc bvr & 0.95 & 0.99 & 0.59 & 1.2 & 0.48 & 26.3 + u & 1.30 & 1.48 & 0.78 & 1.6 & 0.50 & 26.0 + b & 1.25 & 1.29 & 0.71 & 1.4 & 0.47 & 26.1 + v & 0.90 & 0.96 & 0.56 & 1.2 & 0.52 & 26.0 + r & 0.90 & 0.97 & 0.58 & 1.2 & 0.49 & 25.8 + i & 0.80 & 0.99 & 0.56 & 1.2 & 0.49 & 24.7 + @xmath46 & 0.90 & 1.06 & 0.62 & 1.2 & 0.43 & 23.6 +    lllllllll 0&j223023.68 - 604817.4&22.506580&-60.804850&7397.04&3613.83&0.490000&0.377004&0.0771630 + 0.0771630&45.0000&16 & 0.00599998&0.0123559&0.0268263&0.0133826&0.0721889&0.0189753 + 0.0535715&0.0192685&0.0584936&0.0536170&0.105203&0.142492 & 0.0153194&0.0299205&0.0582284 + 0.0300672&0.129481&0.0322477&0.0994542&0.0346933&0.0963358&0.0964758&0.129852&0.292286 + 1&j223024.28 - 604008.9&22.506747&-60.669150&7395.69&5443.56&0.380000&0.381543&0.0771630 + 0.0771630&45.0000&16 & 0.0545032&0.0134836&0.0528900&0.0133148&0.0600472&0.0181627 + 0.0531777&0.0189421&0.0346594&0.0537312&-0.0605122&0.0999073 & 0.123143&0.0319542&0.111140 + 0.0292123&0.0983857&0.0301251&0.0985063&0.0331430&0.0733320&0.0938256&-0.113113&0.287106 + 2&j223023.85 - 604558.3&22.506627&-60.766210&7397.14&4134.88&0.490000&0.555360&0.216003 + 0.133500&-0.200000&24 & 0.0909909&0.0243018&0.0549541&0.0255182&0.0440691&0.0353969 + 0.130563&0.0386372&0.159939&0.108677&0.0605380&0.286208 & 0.100151&0.0285476&0.0699377 + 0.0270144&0.0861040&0.0280836&0.125653&0.0311628&0.113170&0.0863858&0.0282574&0.261025 + 3&j223024.98 - 603027.9&22.506940&-60.507770&7394.68&7619.57&0.940000&0.328410&0.199716 + 0.120951&18.8000&24 & -0.0217973&0.0159112&0.000543353&0.0123818&0.115021&0.0169463 + 300.840&0.0511371&577.641&0.128364&-0.00643368&0.124763 & -0.0458961&0.0422158&-0.000965204 + 0.0289711&0.222445&0.0300914&601.662&0.125434&1191.94&0.301193&0.00430643&0.280111 + 4&j223023.35 - 605242.2&22.506487&-60.878400&7397.96&2622.20&0.610000&0.613032&0.369528 + 0.202653&-5.40000&24 & -0.00917560&0.0288420&0.0469781&0.0383380&0.134234&0.0542395 + 0.155328&0.0622344&-0.104166&0.130648&0.110060&0.458585 & -0.0224367&0.0274903&0.100572 + 0.0275703&0.155969&0.0290837&0.143652&0.0332477&0.0475046&0.0869854&0.677608&0.276257 +