we are poised on the threshold of unprecedented technical growth in wide - field time domain astronomy , where ground - based observations yield very precise measurements of stellar brightness from high - volume data streams .
so far , wide - field time - series surveys has been spearheaded by relatively small telescopes since they are supported by large field of view ( fov ) instruments operating with high duty cycle ( see @xcite for a summary of optical variability surveys ) . within the last decade ,
the advent of large mosaic ccds has facilitated the coverage of large sky area even for large - aperture telescopes ( e.g. , mmt megacam : @xcite ; eso very large telescope omegacam : @xcite ; subaru suprime - cam : @xcite ; chft megacam : @xcite ; iptf : @xcite ) .
although these facilities are generally devoted to imaging surveys , researchers are attempting to utilize them for short- and long - term variability surveys with short - cadence exposures ( e.g. , @xcite ) .
such wide - field imaging systems have enabled us to observe hundred of thousands of target stars simultaneously and also to detect various variability phenomena .
a remarkable thing about these surveys is that the fraction of variable sources increases as the photometric precision of the survey improves .
for this reason , it is important to improve the accuracy in photometry .
another key issue in wide - field time - series photometry is the removal of temporal systematics from a single image frame or several consecutive image frames .
it has recently become known that systematic trends in time - series data can be different and localized within the image frame when the fov is large . such spatially localized patterns may be related to subtle point spread function ( psf ) differences and sky condition within the detector fov ( e.g. , @xcite ) .
as these patterns change in time , we can see how the temporal variations of systematic trends affect the brightness and shape of light curves directly .
the time - scale of systematic variation is sometimes comparable to short - term variability , such as transits or eclipses , and in some cases even long - term variability .
thus , it is often difficult to identify and characterize true variabilities .    in this paper
, we introduce a new photometry procedure , called multi - aperture indexing , which is suited to analyzing well - sampled wide - field images of non - crowded fields with a highly varying psf , such as those produced by wide - field mosaic imagers on large telescopes .
we apply this procedure to archival imaging data from the mmt / megacam transit survey of the open cluster m37 ( hartman et al .
2008a ) , demonstrating a substantial improvement over the existing photometry .
section 2 describes the mmt imaging database and identifies problems in the existing photometry which motivated the development of our new methods .
section 3 describes the multi - aperture photometry that utilize newly defined contamination index and carefully tuned calibration procedures , including the results of the basic tests to validate our approach .
section 4 gives an in - depth discussion about systematic trends in time - series data and suggests an efficient way for identifying , measuring , and removing spatio - temporal trends .
section 5 describes the effects of new calibration on period search , and we summarize our main results in the last section .
@xcite have conducted a study to find neptune - sized planets transiting solar - like stars in the rich open cluster m37 .
the observing strategy was carefully designed for a transiting planet search by several considerations ( e.g. , the reliability of exposure time per frame , the effects of pixel - to - pixel sensitivity variations , and sensitivity of filter ) .
their work did not reveal any transiting planets , but it did provide a rare opportunity to explore photometric variability at relatively high temporal resolution with 3090 s. @xcite discovered 1430 new variable stars , including very short - period eclipsing binaries ( e.g. , v37 , v706 , v1160 ) and @xmath1 sct - type pulsating stars ( e.g. , v397 , v744 , v1412 ) .
we used the same data set on the open cluster m37 .
a detailed discussion of the observations , original data reduction , and light curve production is described in @xcite .
the data archive consists of approximately 5000 @xmath2-filter images taken over 24 nights with the wide - field mosaic imager ( megacam ) mounted at the @xmath3 cassegrain focus of the 6.5 m mmt telescope .
note that megacam is made up of 36 2048 @xmath4 4608 pixel ccd chips in a @xmath5 pattern , covering a 24@xmath6@xmath424@xmath6 fov @xcite .
this instrument has an unbinned pixel scale of @xmath7 , but it was used in @xmath8 binning mode for readout .
the observation logs are summarized in table 1 of @xcite . in brief , the @xmath9-band time - series observations were undertaken between 2005 december 21 and 2006 january 21 , with a median fwhm of @xmath10 arcsec .
exposure times are chosen to keep an @xmath11 mag star as close to the saturation limit , which is expressed as a function of seeing conditions . with an average seeing @xmath00@xmath12.89 on images ,
the quality of the images is good to achieve high - precision light curves ( less than 1% rms value ) down to 20 .
in addition to the imaging data set , this database includes light curve data sets for a total of 23,790 sources detected in a co - added reference image .
theses light curves are obtained by the image subtraction technique using a modified version of isis software ( @xcite ; @xcite ) .    as shown in figure [ fig : fig1 ] , however , the raw light curves from the original image subtraction procedures exhibit many unusual outliers , and more than @xmath015% of data get rejected by a simple filtering algorithm after cleaning procedures . in practice , brutal filtering that is often applied to remove outlying data points can result in the loss of vital data , with seriously negative impact to short - term variations such as flares and deep eclipses .
we also find that the image subtraction technique often resulted in measurement failures from several frames due to poor seeing or tracking problem . after removing these bad frames
, it leads to loss of additional @xmath05% data points from most light curves . in order to overcome this problem ,
we have re - processed the entire image database with new photometric reduction procedures .
3-@xmath13 control limits ( dotted lines ) .
these light curves contain outlier points that significantly increase the rms scatter of the raw light curves . ]
we followed the standard ccd reduction procedures of the bias correction , overscan trimming , dark correction , and flat - fielding as described in @xcite .
the individual ccd frames were calibrated in iraf , using the mosaic data reduction package megared .
the first step is to correct the pixel - to - pixel zero - point differences that are usually described by the sum of a mean bias level and a bias structure .
as the bias frames were not separately taken during the time of the observations , the mean bias level was subtracted from each image extension using an overscan correction and so we can not remove any remaining bias structure from all overscan - subtracted data frames . according to description in matt ashby s megacam reduction guide , the bias structure can be very significant in some small regions such as the portions of the arrays close to the readout leads .
the dark currents are normally insignificant for megacam so that corrections are not needed even for long exposures .
the next step is to correct pixel - to - pixel variations in the sensitivity of the ccd ; we used the program ` domegacamflat2 ` .
this program determines the scaling factors to correct the gain difference between the two amplifiers of each chip by finding the mode in the quotient of pixels to the left and right of the amplifier boundary , and then flattens each of the frames with a master flat field frame .
it is worth mentioning that the sky conditions were rarely photometric during the observing run , with persistent light cirrus for most of the nights .
therefore it was only possible to obtain twilight sky flats on a handful of nights ( dome flats were not possible ) .
we removed bad pixels using the megacam bad pixel masks distributed with the megared package .
the values of bad pixels are replaced with interpolated value of the surrounding pixels using the iraf task ` fixpix ` .
the numerous single pixel events ( cosmic - rays ) were identified and removed using the lacosmic package @xcite .
the megacam data already have a rough world coordinates system ( wcs ) solution that is based on a single value of the telescope pointing .
to update these with a more precise solution , we applied astrometric correction to each ccd in the mosaic using the wcstools imwcs program @xcite .
the new solution is derived by minimizing the differences between the r.a . and decl .
positions of sources in a single ccd chip and their positions listed in the 2mass point source catalog @xcite .
the resulting astrometric accuracy is typically better than 0@xmath14.1 rms in both r.a . and decl .
typically , a point or extended source detection algorithm is applied to each frame independently and it always requires criteria for what should be regarded as a true detection . in obtaining the pixel coordinates for all objects in the m37 fields , this procedure often misses some objects when the detection threshold approaches the noise level .
also it needs a substantial effort to match the objects that are detected in only some of the frames .
our approach is as follows : a complete list of all objects is obtained from a co - added reference frame , and then the photometry is performed for each frame using the fixed positions of the sources detected on the reference . since the relative centroid positions of all objects are the same for all frames in the time series , we can easily place an aperture on each target and measure the flux even for the stars at the faint magnitude end .
we constructed the reference frame for each chip from the best seeing frames using the swarp software .
benefiting from a highly accurate astrometric calibration of input frames , we were able to improve the quality of co - added images . in the swarp implementation
, the pixels of each frame were resampled using the lanczos3 convolution kernel , then combined into the reference frame by taking a median or average .
after this was done , sources were detected and extracted on the reference frame using the sextractor software @xcite .
when configured with a lower detection threshold , sextractor extracts the number of spurious detections ( e.g. , diffraction spikes around bright stars , or outer features of bright galaxies ) .
these false detections were removed by careful visual inspection for each chip .
the final catalog contains a total of 30,294 objects including both point and extended sources .
prior to the photometry , the initial centroid coordinates of the target objects for each frame were computed by using the wcstools sky2xy routine @xcite .
the stored world coordinate system for each frame is used to convert the ( r.a .
, decl . ) coordinates from the master source catalog to the @xmath15 pixel locations .
however , actual positions of objects for each frame can be slightly moved from its original locations depending on the focus condition of instrument , seeing condition , and the signal - to - noise ratio ( s / n ) in the individual observations .
these types of positioning errors ( i.e. , centroid noise ) will lead to the internal error for a photometric measurement that results from placement of the measuring aperture on the object being measured .
the situation gets worse for faint stars because the centroid position of them is itself subject to some uncertainty . the fractional error in the measured flux as a result of mis - centering
is given by:@xmath16 where @xmath17 is the positioning error , @xmath18 is the radius of aperture , @xmath13 is the profile width for a source with a gaussian psf , and total flux @xmath19 ( see appendix a in irwin et al . , 2007 for details ) .
this expression shows that if we set the aperture radius equal to the fwhm ( @xmath18=2.35-@xmath13 ) , even small differences in placement of the aperture ( e.g. , @xmath17=0.1-@xmath13 ) may increase the uncertainty in the flux measurements ( @xmath20 1 mmag ) .
thus , accurate centroid determination is important to achieve the high - precision photometry .    following the windowed centroid procedure in the sextractor
, a refined centroid of each object is calculated iteratively .
on average , the rms uncertainty in the coordinate transformation using the wcs information was @xmath21 pixels for the bright reference stars .
after the coordinate transformation from sky to @xmath22 , however , the centroid coordinates ( @xmath23 , @xmath24 ) are slightly misaligned from their actual ones ( @xmath25 , @xmath26 ) .
the refined centroid values ( @xmath27 , @xmath28 ) are used only if the maximum displacement is at least less than 1.5 pixels .
this condition prevents arbitrary shifting of a source centroid , especially for faint stars .
we estimated a local sky background by measuring the mode of the histogram of pixel values within a local annulus around each object , which is suitable choice for our uncrowded field ( less than @xmath01000 stars per chip ) .
this process is a combination of @xmath29-@xmath13 clipping and mode estimation .
the background histogram is clipped iteratively at @xmath303-@xmath13 around its median , and then the mode value is taken as : @xmath31 it represents the most probable sky value of a randomly chosen pixel in the sample of sky pixels @xcite .
for relatively crowded regions , we utilized a background map created by sextractor package using a mesh of @xmath32 pixels and a median filter box of @xmath33 pixels .
this map is used to confirm the properness of individual sky values from annulus estimates .
modern data reduction techniques aim to reach photon noise limit and minimize systematic effects .
for example , differential photometry technique can be achieved better than 1% precision for brighter stars ( e.g. , @xcite ) , and the deconvolution - based photometry algorithm leads to the minimization of systematic effects in very crowded fields ( e.g. , @xcite ) . however
, conventional data reduction methods often fail to handle various artifacts in wide - field survey data .
we present below a new photometric reduction method for precise time - series photometry of non - crowded fields , without the need to involve complicated and cpu intensive process ( e.g. , psf fitting or difference image analysis ) .
our photometry is similar to standard aperture photometry , except in that we compute the flux in a sequence of several apertures and then determine the optimum aperture individually to each object at each epoch .
this multi - aperture photometry is an efficient way to determine the optimum aperture size that gives the maximum s / n for a flux measurement . the maximum s / n is not necessarily at the same aperture for all objects , and it can be obtained from a relatively small aperture @xcite .
this _ photometric _
aperture is to achieve the optimal balance between flux loss and noises based on a relationship derived from the ccd equation ( see @xcite ) .
figure [ fig : fig2 ] shows how the optimum apertures vary with the stellar magnitude .
there is an obvious trend of decreasing aperture sizes with increasing magnitudes down to the faint magnitude limit in the example frame .
once we measure the flux of each object with the optimum aperture , we need to apply the aperture correction for small apertures .
the aperture correction terms are estimated from the growth curve analysis of selected isolated bright stars ( i.e. , reference stars ) .
the average curve - of - growth for each frame is calculated by measuring the difference in magnitude between different pairs of apertures ( up to 10 pixels aperture radius ) and then an automatic correction is applied to all objects for each photometric aperture .
the use of a common aperture correction for each ccd assumes that there is no variation in the correction across the ccd .
this flux correction method gives nearly the same brightness within the measurement uncertainties for all apertures .
any psf variation across the ccd causes systematic errors , however , and we deal with this in section 3.4 and section 4 .    100 thumbnail images of the target star . _
middle panels _ : normalized s / n as a function of aperture size . _ bottom panels _ : aperture corrected magnitude as a function of aperture size .
the arrows represent the peak locations in the aperture - s / n diagram ( see text for details ) . ]
we performed the multi - aperture photometry based on the concentric aperture photometry algorithm in daophot package @xcite , using several circular apertures ( up to 10 pixels aperture radius ) with a fixed sky annulus from 35 to 45 pixels .
the initial results of multi - aperture photometry are stored in ascii - format photometry tables , including the date of the observations ( mjd ) , the pixel @xmath15 coordinates , the aperture - corrected magnitudes with errors for each aperture , the sky values and its errors .
figure [ fig : fig3 ] shows the details of the multi - aperture photometry for one star at different epochs . in the former two epochs , the photometric apertures can be properly selected by the s / n cuts , while in the latter two epochs , s / n increases for lager apertures .
this unusual behavior is due to contamination by a moving object .
we automatically identifies similar unusual cases by the method of aperture indexing .
-filter light curve of same star ( id=10213 ) as shown in figure [ fig : fig3 ] .
bottom panels show the multi - aperture indexing scheme .
the @xmath34-axis is the aperture size and @xmath35 @xmath36-axis is the differential magnitude between pairs of apertures @xmath37(=@xmath38 ) .
we can see whether and at what aperture the differential magnitude ( solid lines ) begins to deviate from the model curve ( dashed lines ) for each epoch . ]
our multi - aperture indexing method is similar to the basic concept of the discrete curve - of - growth method @xcite .
each object is indexed based on the difference in aperture - corrected magnitude between pairs of apertures @xmath37(=@xmath38 ) with mean trend for stars of similar brightness ( see solid and dotted lines in bottom panels of figure [ fig : fig4 ] , respectively ) .
the aperture with a 10 pixel radius is used as the fixed reference aperture .
the mean trend is determined by computing the rms curve of the aperture correction values for all measured apertures , and used to evaluate whether magnitude at a given aperture significantly differs from the mean trend .
since the rms value depends on the chosen magnitude interval , all stars are divided into groups according to their brightness in the individual frames .
we determine the rms curve for each magnitude group using an iterative @xmath13-clipping until convergence is reached .
objects lying within @xmath303-@xmath13 of the model curve are indexed as contamination - free , and those above @xmath303-@xmath13 as a contaminated source .
figure [ fig : fig4 ] shows that multi - aperture indexing guides us to throw out some photometric measurements if they are discrepant from the mean trend .
this approach also gives us a chance to recover a measurement that would be otherwise thrown out .
the problematic aperture can be simply replaced by one of the smaller apertures if it is indexed as contamination - free .
this help us make a full use of the information offered by the data .
we present a new photometric calibration to convert the instrumental magnitudes onto the standard system , including a relative flux correction of the left and right half - region of each ccd chip .
as mentioned in the section 3.1 , mmt / megacam shows the temporal variations in the gain between two amplifiers on each ccd , as well as between ccds that are part of the same mosaic .
it may have been caused by unstable bias voltage of the ccd output drain which has a profound impact on the gain of the output amplifier .
the level of readout noise is also unstable between two amplifiers . to correct for this effect
, the photometric calibration needs to be performed individually for each amplifier region .
we use a sufficient number of ( pre - selected ) bright isolated stars as standard stars and compute the relative flux correction terms .
these terms were derived for each frame using the mean magnitude offset ( @xmath39 ) of standard stars ( @xmath40 ) with respect to corresponding magnitudes ( @xmath41 ) in the master frame chosen as an internal photometric reference    @xmath42    where @xmath43 is the aperture - corrected instrumental magnitudes in other frames and @xmath44 is the photometric zero - points for the left- and right - side of each chip , respectively . to calculate the zero - points , we solve
a linear calibration relation of the form : @xmath45 where @xmath46 and @xmath47 are standard magnitudes from the photometric catalog of m37 @xcite and @xmath48 is an airmass term .
the fit is performed iteratively using a sigma - clipping method .
figure [ fig : fig5 ] shows the difference in the magnitude offsets between the left- and right - side of each chip ( @xmath49 ) for the whole data set , which is within @xmath50 magnitude level for all 36 ccd chips .
the histograms are normalized by the total number of data frames @xmath51 and are described by a gaussian function with slightly different mean values and shapes ( dashed line ) .
we clearly see a significant variation in difference between a pair of magnitude offsets for all ccds .      the photometric calibration for wide - field imaging systems is also affected by position - dependent systematic errors due to a psf variation across the fov ( e.g. , @xcite ) .
we derive psf variations across the fov with the sextractor package .
the change of the psf shapes in the image plane is represented by spatial distribution of psf fwhm values for several bright stars , with parameters of ` class\_star ` @xmath52 0.9 , ` magerr\_auto ` @xmath53 0.01 mag , and ` flags ` = 0 ( i.e. , isolated point sources with no contamination ) .
note that the psf fwhm values are defined as the diameter of the disk that contains half of the object flux based on a circular gaussian kernel .
figure [ fig : fig6 ] presents the variation of the psf fwhm as a function of distance from the image center for various seeing conditions . for each quadrant ,
the dashed lines represent the weighted spline approximation of the median value of each distance bin ( 1 arcmin ) .
the result shows that the psf fwhm varies significantly as a function of position on the single - epoch image frames and variations are at the level of @xmath010% to 20% ( @xmath54 ) across the fov . as the field distortion is not negligible from the center of field to its edges ,
such variations limit the accuracy of stellar photometry .        to address this issue , we perform a 2d polynomial fitting technique . for each frame , the correction terms are described by a linear or quadratic polynomial depending on the position @xmath55 only .
@xmath56    where @xmath57 are the pixel coordinates of @xmath58 bright isolated stars , @xmath59 is the statistical weight in the fitting procedure , @xmath60 are the sets of polynomial coefficients for each aperture size , and @xmath61 are the difference in magnitude between the reference aperture and @xmath62 aperture , @xmath63 , at the position @xmath55 for each chip @xmath64 .
we derived the optimal parameter values from a nonlinear least - squares fit using the levenberg
marquardt algorithm and automatic differentiation , and choose between two models that best fit the data .
figure [ fig : fig7 ] shows the field - dependent magnitude offsets and the distortion correction by 2d polynomial fitting method for one example mosaic ccd . for the outer and the central region of the mosaic , we compare the magnitude offsets between the reference aperture and the relatively smaller apertures as a function of @xmath55 coordinates . here
@xmath34-axis is in the declination direction and @xmath36-axis is opposite to the right ascension direction .
we find that the magnitude difference depends on position @xmath15 and is most discrepant in the outer part of the fov .
this effect is usually more significant in the @xmath36-direction than in the @xmath34-direction , especially for the case of aperture photometry performed with small apertures .
the correction for field - dependent psf variation reduces the initial @xmath010% variation ( gray lines ) to less than @xmath01% ( black lines ) .    .
]        we compare the photometric performance of the re - calibrated light curves with the non - de - trended archival light curves by means of the two representative measures : ( i ) the rms photometric precision @xmath65 , and ( ii ) the data recovery rate @xmath66 .
the former is defined as the standard deviation of light curves around the mean value as a function of @xmath46 magnitude : @xmath67 where @xmath68 is the number of data points in each light curve , @xmath69 is the observed magnitude , and @xmath70 is the mean magnitude of the object @xmath71 , and the latter refers to the number of _ analyzed _ data frames normalized by the total number of observed data frames @xmath58 for each object . in typical cases , the data recovery rate should be near unity in the bright magnitude regime and decreases with magnitude for fainter objects . for comparison , we decided to select light curve samples which show either no significant variability or seeing - correlated variations induced by image blending .
we remove all known variable stars from the sample list based on a new catalog of variable stars in m37 field ( s .- w .
chang et al .
2015 , hereafter paper ii ) . to remove the light curves of blended objects
, we use an empirical statistical technique to quantify the level of blending by looking for seeing - correlated shifts of the object from its median magnitude @xcite .
@xmath72 where @xmath73 for light curve points @xmath69 with uncertainties @xmath74 , and @xmath75 is the same statistic measured with respect to a fourth order polynomial in fwhm fitted to the data .
we adopt the value @xmath76 for the selecting light curves with no blending .
the last selection criterion is that the light curves must exit both in the archive and our database .    in the bottom panels of figure [ fig : fig8 ] and figure [ fig : fig9 ] , we plot the rms photometric precision of light curves for the two megacam ccd chips in the outer ( ccd 1 ) and central ( ccd 21 ) part of the fov , respectively .
the black points show the rms values of the re - calibrated light curves , while the gray points are for the raw and filtered light curves in archive .
the first impression from this comparison is that the typical rms scatter is overestimated from the raw light curves because of many outliers in the photometric data ( bottom left panels ) . for the better results ,
these light curves were filtered out in two steps : ( i ) clipping 5-@xmath13 outliers from each light curve and ( ii ) removing every data points that are outliers in a large number of light curves @xcite . in the second step ,
the outlier candidates are estimated by choosing a cutoff value for each ccd chip .
the cutoff value is defined as the fraction of light curves for which a given image is a 3-@xmath13 outlier .
this filtering was applied to remove bad measurements due to image artifacts or poor conditions , which were previously thought to be unrecoverable , but resulting data loss is up to 20% of the total number of data points ( bottom right panels ) . as shown in the top panels of the two figures ,
the data recovery rate for the re - calibrated light curves is close to 100@xmath77 over a wide range of magnitude and it appears to be more complete compared with the raw ( top left panels ) and filtered ( top right panels ) light curves . at bright magnitudes ( @xmath78 )
the data recovery rate does not reach 100@xmath77 because the exposure time was chosen to be saturated at a magnitude of @xmath11 .
this comparison proves that our approach is a powerful strategy for improving overall photometric accuracy without the need to throw out many outlier data points .
finally , we compare the light curves themselves for selected variable stars between the archive and our own .
this comparison serves to illustrate how the photometric precision and data recovery rate of the time - series data affect the ability to address a variety of variability characteristics .
figure [ fig : fig10 ] shows a direct comparison with the filtered light curves ( top panels ) and our re - calibrated light curves ( bottom panels ) for four variable stars .
it is shown that our method recovers more data points ( black ) from the same data set of images .
and @xmath36 coordinates , respectively , for selected ccd chips .
the filled squares are relatively bright stars with @xmath79 mag , while the open squares are stars with @xmath80 mag . ]      in order to further investigate possible systematics in our approach , we conducted psf - fitting photometry with daophot ii and allstar @xcite .
for each mosaic frame , we select bright , isolated , and unsaturated stars to make the psf model varying quadratically with @xmath55 coordinates .
after psf modeling , we run allstar to perform iterative psf photometry of all detected sources in the frame with initial centroids set to the same values used for our own photometry .
we then calculated aperture corrections using the package daogrow @xcite after subtraction of all but psf stars , which creates aperture growth curves for each frame and then integrates them out to infinity to obtain a total magnitude for each psf star .
the final step is to convert the instrumental magnitudes into the standard photometric system . for each frame , the initial zero - point correction is applied by correcting the magnitude offset with respect to the master frame .
this places photometry for all frames on a common instrumental system .
following the same procedure in section 3.3 , the photometric calibration is performed individually for each amplifier region .
figure [ fig : fig11 ] shows the residual magnitudes and sky values between our multi - aperture photometry and the psf - fitting photometry as a function of position in the selected ccd chips .
there are no position - dependent trends in the magnitude residuals . for the brighter stars with @xmath79 mag ,
the rms magnitude difference between the two methods is very small ( @xmath81 = @xmath82 and @xmath83 = @xmath84 , respectively ) , while for the relatively faint stars the rms difference is somewhat larger ( @xmath81 = @xmath85 and @xmath83 = @xmath86 , respectively ) .
the results of this example indicate that we can reliably correct for the psf variations by our calibration procedures
. meanwhile , our sky values are slightly higher than the allstar sky values , but not to the degree that can seriously affect photometric measurements .
magnitude in the central region of the open cluster m37 ( left panel ) .
the arrows indicate the three different magnitude levels from bright to faint in our sample shown in the right panels .
note that the variable object ids are taken from the new variable catalog of the m37 ( see paper ii ) . ]
figure [ fig : fig12 ] shows a comparison of the rms dispersion of the light curves obtained with our photometry with respect to the that of the psf - fitting photometry in the central region of the open cluster m37 .
we only compare the light curves of non - blended objects as described in section 3.5 .
our multi - aperture photometry does not reach the same level of precision as psf - fitting photometry for the faintest stars , while the psf - fitting approach results in poorer photometry for bright stars .
as shown in the right panels of figure [ fig : fig12 ] , it is clear that our photometry tend to have smaller measurement errors with respect to the psf photometry for the bright stars .
-axis indicate the corresponding timestamps in each frame .
the deviations from the mean value , @xmath87 , are less than @xmath300.01 mag level ( with a rms values of 0.0021 , 0.0025 , and 0.0033 mag from the top panel down ) .
these stars show a similar pattern of light variations over the observation span . ]
from a visual inspection of the re - calibrated light curves in the same ccd chip , we found that some light curves tend to have the same pattern of variations over the observation span ( figure [ fig : fig13 ] ) .
this kind of systematic variation ( i.e. , _ trend _ ) is often noticed in other studies .
for example , the importance of minimizing known ( or unknown ) systematics have been recognized by several exo - planet surveys because planet detection performance can be easily damaged by them ( e.g.,@xcite ) .
also space - based time - series data ( e.g. , corot and kepler ) are no exception to this behavior although it is completely free from systematics caused by the turbulent atmosphere .
most of the raw light curves are affected by a secular ( or a sudden ) variation of flux without any obvious physical reason @xcite .    in order to check the properties of temporal systematics
, we examined the correlation coefficients as measure of similarity between two light curves @xmath47 and @xmath71 obtained from a single ccd chip ( ccd 1 ) .
@xmath88 where @xmath89 is the flux of each star at time @xmath90 , @xmath58 is the total number of measurements , @xmath91 is the mean flux of each star , and @xmath13 is the standard deviation of @xmath89 .
this comparison is a point - by - point comparison and is done for every pair of light curves in the data set .
the resultant similarity matrix can be used to identify correlated pairs of light curves and to determine which light curve is least like all other light curves ( e.g. , @xcite ) .
after that , we selected stars showing most systematics based on a hierarchical clustering method with the correlation coefficients ( see @xcite , for more details ) .
figure [ fig : fig14 ] represents spatial distribution of the most prominent trend groups on the ccd plane ( top panel ) and its strongly correlated features determined by the weighted sum of normalized light curves ( bottom panel ) .
there are two interesting features in this figure : the first one is that each trend covers only a certain part of the sky area and the second one is that some portions of neighboring trends show different variation patterns even at the same moment in time ( shaded gray region in the figure ) . in particular , we found an anti - correlated variation for the trends between the group 1 ( @xmath92 ) and the group 2 ( @xmath93 ) , so we might expect to find possible noise sources that are responsible for these discrepancies .
why the trends are different and localized within a single ccd frame is a subject of further study , but it is probably related to subtle changes in point spread function and sky condition within the detector fov .     and @xmath94 ) of stars .
the numbers on the @xmath34-axis indicate the corresponding timestamps between the frame 1 and the frame 400 . from the top to bottom , the panels show the average difference of trend , flux concentrations in which we applied a @xmath95 mag shift to the @xmath96 value , sky level , and psf fwhm value between the two groups . ]    , for the original multi - aperture magnitude measurements as a function of aperture size , which is marked with a gray dots in all panels .
a positive @xmath97 indicates that the photometric measurements with the corresponding aperture size , @xmath98 , are brighter than those of reference aperture @xmath99 , while a negative @xmath97 indicates vise versa .
the dashed lines are the rms model profiles introduced in section 3.2.2 .
there is a noticeable distinction between the group 1 ( _ top panels _ ) and group 2 ( _ bottom panels _ ) when looked at different concentration levels .
but the trends are not fully explained by the difference patterns in @xmath97 because those correlated variations become small after the correction for the psf variation ( black dots ) . ]    for these two groups , we consider a possible causal relationship between the systematic trends and average object / image properties .
figure [ fig : fig15 ] shows the differences in trend , differential magnitudes , sky level , and psf fwhm between the two groups , respectively . in the top panel , we plot the magnitude difference in trends , which shows variations in the range of @xmath50 mag .
we suspect that this may be due to the different concentration of star light between these two groups
. it can be checked by using the magnitude difference @xmath100 , where @xmath99 and @xmath98 are the reference aperture and the relatively small aperture , respectively .
in fact , we already know that there is a magnitude variation in @xmath37 depend on the aperture size due to the field - dependent psf variation ( see section 3.4 ) . for example , figure [ fig : fig16 ] shows the response of multi - aperture photometry for the two groups of stars at one epoch ( mjd = 53726.14817 ) before and after applying the distortion corrections .
although the magnitude variation between the group 1 and the group 2 seems to have different behavior as a function of aperture size ( gray points ) , it is negligible after the removal of the psf variation ( black points ) .
we also check that the possible contribution of sky level ( @xmath101 ) and psf fwhm differences ( @xmath102 ) to the systematic trends on the re - calibrated light curves . as mentioned by @xcite , sky over - subtraction
may lead to the systematic trends as a function of the psf fwhm , the amplitude of which increase for fainter stars .
the third and forth panels of figure [ fig : fig15 ] show the variation of the mean @xmath101 and @xmath102 , respectively . in our case , however , the form and amplitude of trends seem independent of sky level and psf fwhm .
some other possible sources that may contribute to the observed systematic trends include : higher order variations in the psf shape beyond just the fwhm ; cross - talk from other amplifiers , or ghosts from bright stars undergoing multiple reflections within the optics ; non - uniform variations in the gain ; and unmodeled temporal atmospheric variations that are dominated by rayleigh scattering , molecular absorption by ozone and water vapor , and aerosol scattering @xcite .
while we find a clear presence of trends that should be removed , we are not able to identify their exact cause .      in order to reduce systematic effects in photometric time - series data , several methods were introduced ( e.g. , tfa : @xcite ; sys - rem : @xcite ; pdt : @xcite ; cda : @xcite ; pdc : @xcite ) .
all of these algorithms share a common advantage that they work without any prior knowledge of the systematic effects .
we use the pdt algorithm , which has been designed to detect and remove spatially localized patterns . by default
, this algorithm works with a set of light curves that contain the same number of data points distributed in the same series of epoch . in many cases , however ,
missing data occur when no photometric measurements are available for some stars in a given observed frame .
these missing data can be simply replaced by means , medians , or the values from the interpolation of adjacent data points in each light curve ( e.g. , @xcite ) .
although using the replaced value is the easiest way to reconstruct the light curve to be analyzed , it is not appropriate if the time separation between two subsequent observations is too large .
instead we use more straightforward approach by applying the pdt algorithm in two separate steps : ( i ) we construct the master trends from the subset of bright stars , and ( ii ) de - trend light curves of all stars with most similar master trend and matching time line .    -axis indicate the corresponding timestamps between the frame 1 and the frame 2000 .
the @xmath36-axis is @xmath46-filter magnitude ( normalized by its mean value ) .
while the morphology of the two light curves in the left panels appear to be variable stars of some kind , these turn out to be non - variable after applying the photometric de - trending method . in the case of the right panels ,
all true variabilities are preserved from the raw light curves . from upper left
to lower right : id=10032 , id=10039 , id=170088 , and id=170108 . ]
we briefly describe the main procedure of our de - trending process following the algorithm derived by @xcite .
we first select the template light curves from bright stars that show the highest correlation in the light - curve features .
the total length of template light curves should be long enough to cover the whole time span of observations . in this step , we take a sequence of data points , @xmath103 as the reference time line . using the correlation matrix calculated from equation ( 7 ) , we extract all subset of light curves that show spatio - temporally correlated features ( i.e. , clusters ) .
each cluster is determined by hierarchical tree clustering algorithm based on the degree of similarity .
next , we obtain master trends @xmath104 for each cluster by weighted average of the normalized differential light curves , @xmath105 : @xmath106 @xmath107 @xmath108 where @xmath109 is the total number of light curves in each cluster @xmath64 , @xmath110 is mean value of _ i light curve , and @xmath111 is the standard deviation of @xmath105 . after determining the master trends , we de - trend the light curves of all stars with matching master trend and time line .
we adjust the temporal sequence of measurements for the master trends @xmath112 by that of individual light curves to be de - trended @xmath113 . because each light curve is assumed as a linear combination of master trends and noise
, we can determine the optimal solution by minimizing the residual between the master trends and the light curve : @xmath114 @xmath115 where @xmath116 is the total number of master trends and @xmath117 are free parameters to be determined by means of minimization of noise term @xmath118 .
_    figure [ fig : fig17 ] shows examples of our light curves before and after removing the systematic trends .
the algorithm we used for de - trending removes only the systematic variations that are shared by light curves of stars in the adjacent sky regions ( left panels ) , while all kinds of true variabilities are preserved ( right panels ) .
( indicated by arrow ) .
the phased diagram of this candidate frequency is shown in the low - right panel .
we see that the model fits well the overall pulsating variability ( red line ) . ]    , but the subfigure of low - right panel is 200@xmath4200 thumbnail image of the target star .
there is no potential sources of contamination to hamper the interpretation of the power spectrum . ]
the usefulness of our photometry is tested for a set of variable stars .
we immediately find abundant cases of improvements in the following three aspects : ( i ) refinement of the derived period , ( ii ) detection of a new significant peak in the periodogram , and ( iii ) separation of non - variable candidates where systematics in the light curves were mistaken for true variability .
for each case , we compare light curves and power spectra for archival data and our data .
for the first example , we show that a new photometric measurement and calibration allowed us to derive a much improved refinement of the light curves and of the derived periods ( figure [ fig : fig18 ] ) . we performed a lomb - scargle ( l - s : @xcite ) search of both archival and new light curves for periodic variable star ( v427 ) .
the light curves are folded by the best - fit period of 5.4615 and 4.4158 days , respectively .
we also calculated the false - alarm probability ( @xmath119 fap ) for each peak and its signal to noise ratio ( s / n ) : @xmath119 fap@xmath120 = @xmath121 , s / n@xmath120 = 43.2 for the archival data and @xmath119 fap@xmath122 = @xmath123 , s / n@xmath122 = 84.5 for the new data .
since we can get better estimation with much lowered minimum fap value , our new period is the most likely result . in the bottom panels ,
the resulting amplitude spectrum was calculated with ` period04 ` package .
since the archival data is more noisy than the new one , it is rather complicated to interpret the peaks of its power spectrum .    for the second example , we show the newly discovered low - amplitude pulsating variable star ( figure [ fig : fig19 ] ) .
we used the ` period04 ` package to find multiple pulsation periods .
the whole process of identifying , fitting , and pre - whitening successive frequencies was repeated until no significant frequencies were found .
we adopt a conservative approach in selecting the statistical significant peaks from the amplitude spectrum
. a s / n amplitude ratio of 4.0 is a good criterion for independent frequencies , equivalent to 99.9% certainty of variability @xcite .
while no clear periodicity was found in the archive data , our amplitude spectrum shows a clear excess of power centered at 22.3979 days@xmath124 with peak amplitudes of about 1 mmag ( s / n@xmath122 = 9.02 ) .
the last example is the opposite case of the second .
figure [ fig : fig20 ] shows that this object is unlikely to be a variable source because there is no evidence for any significant peaks , which indicates that the variations are mostly noise .
extensive study on variabilities will be presented in paper ii .
in this paper , we introduce a new time - series photometry with multi - aperture indexing and spatio - temporal de - trending techniques , together with complex corrections to minimize instrumental biases .
we used the archival , high - temporal time - series data from one - month long mmt / megacam transit survey program .
the re - calibration of the archival data has made several improvements as follows : ( i ) the photometric information derived from the multi - aperture indexing measurements is useful to obtain the best s / n measurement , but also to diagnose whether or not the targets are contaminated ; ( ii ) the resulting light curves utilize nearly 100% of available data and reach precisions down to sub mmag level at the bright magnitude end without the need to throw out many outlier data points , which makes it possible to preserve data points that show intrinsic sudden variations such as flare events ; ( iii ) corrections for position - dependent psf variations and de - trending of spatio - temporal systematic trends improve the quality of light curves ; and ( iv ) new photometry enables us to determine the variability nature and period estimate more accurately .    while this study deals with a particular set of data from mmt
, we find our approach has a potential for other wide - field time - series observations .
multi - aperture indexing measurement is a powerful tool in isolating and even correcting various contaminations .
spatio - temporal de - trending is also very useful in removing systematics caused by psf variation and even non - uniform extinction of thin clouds across the fov .
@xcite proved this for different sets of archival survey data .
this research was supported by basic science research program of the national research foundation of korea ( 2011 - 0030875 ) .
i.b . is grateful for support from kasi - yonsei drc program of korea research council of fundamental science and technology ( drc-12 - 2-kasi ) .
we thank the mmt / m37 survey team for the kind provision of raw image data .
kim , d .- w . helped us to test the modified photometric de - trending algorithm .
additionally , we would like to thank the anonymous referee for many helpful comments , including the suggestion to use psf photometry to further justify the successful removal of systematics ( section 3.5.2 ) .
alard , c. 2000 , , 144 , 363 alard , c. & lupton , r. h. 1998 , , 503 , 325 becker , a. c. , wittman , d. m. , boeshaar , p. c. , et al .
2004 , , 611 , 418 bertin , e. & arnouts , s. 1996 , , 117 , 393 bianco , f. b. , protopapas , o. , mcleod , b. a. , et al .
2009 , , 138 , 568 boulade , o. , charlot , x. , abbon , p. , et al .
2003 , , 4841 , 72 bramich , d. m. & freudling , w. 2012 , , 424 , 1584 breger , m. , stich , j. , garrido , r. , et al .
1993 , , 271 , 482 chang , s .- w . & byun , y .- i .
2013 , in astronomical society of the pacific conference series , vol .
475 , astronomical data analysis software and systems xxii , ed .
d. n. freidel , 41 everett , m. e. & howell , s. b. 2001 , , 113 , 1428 gillon , m. , pont , f. , moutou , c. , et al .
2007 , , 466 , 743 hartman , j. d. , stanek , k. z. , gaudi , b. s. , holman , m. j. , & mcleod , b. a. 2005 , , 130 , 2241 hartman , j. d. , gaudi , b. s. , holman , m. j. , et al . 2008a , , 675 , 1233 hartman , j. d. , gaudi , b. s. , holman , m. j. , et al . 2008b , , 675 , 1254 hartman , j. d. , gaudi , b. s. , holman , m. j. , et al .
2009 , , 695 , 336 hodgkin , s. t. , irwin , m. j. , hewett , p. c. , & warren , s. j. 2009 , , 394 , 675 howell , s. b. 1989 , , 101 , 616 irwin , j. , irwin , m. , aigrain , s. , et al .
2007 , , 375 , 1449 ivezi , z. , smith , j. a. , miknaitis , g. , et al .
2007 , , 134 , 973 jenkins , j. m. , caldwell , douglas a. , chandrasekaran , h. , et al .
2010 , , 713 , l87 kim , d .- w . ,
protopapas , p. , alcock , c. , byun , y .-
i . , & bianco , f. 2009 , , 397 , 558 kim , d .- w . , protopapas , p. , alcock , c. , byun , y .- i .
, et al . 2010 , , 139 , 757 kovcs , g. , bakos , g. , & noyes , r. w. 2005 , , 356 , 557 kuijken , k. , bender , r. , cappellaro , e. , et al .
2002 , the eso messenger , 110 , 15 kulkarni , s. r. 2013 , the astron .
telegram , 4807 , 1 magain , p. , courbin , f. , gillon , m. , et al .
2007 , , 461 , 373 mazeh , t. , guterman , p. , aigrain , s. , et al .
2009 , , 506 , 431 mcleod , b. a. , conroy , m. , gauron , t. m. , geary , j. c. , & ordway , m. p. 2000 , in further developments in scientific optical imaging , ed .
m. b. denton ( cambridge : royal society of chemistry ) , 11 merline , w. & howell , s. b. 1995 , exp .
6 , 163 . mink , d.j .
2002 in astronomical society of the pacific conference series , vol .
281 , astronomical data analysis software and systems xi , ed .
d. a. bohlender , d. durand , & t. h. handley , 169 mislis , d. , schmitt , j. h. m. m. , carone , l. , guenther , e. w. , & ptzold , m. 2010 , , 522 , 86 miyazaki , s. , komiyama , y. , sekiguchi , m. , et al .
2002 , pasj , 54 , 833 padmanabhan , n. , schlegel , d. j. , finkbeiner , d. p. , et al .
2008 , , 674 , 1217 pepper , j. , stanek , k. z. , pogge , r. w. , et al .
2008 , , 135 , 907 pietrukowicz , p. , minniti , d. , fernndez , j. m. , et al .
2009 , , 503 , 651 pont , f. , zucker , s. , & queloz , d. 2006 , , 373 , 231 protopapas , p. , giammarco , j. m. , faccioli , l. , et al .
2006 , , 369 , 677 randall , s. k. , calamida , a. , fontaine , g. , et al .
2011 , , 737 , 27 scargle , j. d. 1982 , , 263 , 835 skrutskie , m. f. , cutri , r. m. , stiening , r. , et al .
2006 , , 131 , 1163 stetson , p. b. 1987 , , 99 , 191 stetson , p. b. 1990 , , 102 , 932 tamuz , o. , mazeh , t. , & zucker , s. 2005 , , 356 , 1466 twicken , j. d. , chandrasekaran , h. , jenkins , j. m. , et al .
2010 , , 7740 , 62 van dokkum , p. g. 2001 , , 113 , 1420