IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Ignore:
Timestamp:
Dec 2, 2019, 5:19:53 PM (7 years ago)
Author:
eugene
Message:

adding comments from bertrand

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/doc/release.2015/ps1.calibration/calibration.tex

    r41188 r41189  
    115115bright-star systematic error floor for individual astrometric
    116116measurements is 16 milliarcseconds.}  \textmod{The Pan-STARRS Data Release 2
    117 (DR2) astrometry is tied to the Gaia DR1 coordinate frame with a
     117(DR2) astrometric system is tied to the Gaia DR1 coordinate frame with a
    118118systematic uncertainty of $\sim 5$ milliarcseconds.}
    119119\end{abstract}
     
    132132Type Ia supernovae to measure the history of the expansion of the
    133133universe.  The majority of the time (56\%) was spent on surveying the
    134 $\frac{3}{4}$ of the sky north of $-30$ Declination with
     134$\frac{3}{4}$ of the sky north of $-30\degree$ Declination with
    135135\grizy\ filters in the so-called $3\pi$ Survey.  Another $\sim 25\%$
    136136of the time was concentrated on repeated deep observations of 10
     
    571571derived from the projection of the reference coordinates.  One caveat
    572572is that the chip reference coordinates are also degenerate with the
    573 fitted distortion.  In order to avoid being sensitive to the exact
     573fitted distortion.  \textmod{To avoid} being sensitive to the exact
    574574positions of the chips at this stage, we measure the local gradient
    575575between the focal plane and tangent plane coordinate systems.  We then
     
    635635\begin{itemize}
    636636\item \code{ZPT_REF} : the nominal zero point for this filter
    637 \item \code{ZPT_OBS} : the measured zero point for this chip /
    638   exposure
    639 \item \code{ZPT_ERR} : the measured error on \code{ZPT_OBS}
     637\item \code{ZPT_OBS} : the measured zero point for this chip / exposure
     638\item \code{ZPT_ERR} : the standard deviation of \code{ZPT_OBS}
    640639\item \code{ZPT_NREF} : the number of stars used to measure \code{ZPT_OBS}
    641640\item \code{ZPT_MIN} : minimum reference magnitude included in analysis
     
    663662
    664663The master DVO database is used to perform the full photometric and
    665 astrometric calibration of the data.  During these analysis steps, a
    666 wide variety of conditions are noted for individual measurements, for
    667 the objects (either as a whole or for specific filters) and for the
    668 images.  A set of bit-valued flags are used in the database to record
    669 these conditions.
     664astrometric calibration of the \textmod{PS1} data.  During these
     665analysis steps, a wide variety of conditions are noted for individual
     666measurements, for the objects (either as a whole or for specific
     667filters) and for the images.  A set of bit-valued flags are used in
     668the database to record these conditions.
    670669%
    671670Table~\ref{tab:measure_mask_values} lists the flags specific to
     
    868867calibrations.
    869868
    870 Photometric nights are selected and all other exposures are ignored.
     869In this first stage, the goal is to determine an initial
     870highly-reliable collection of zero points for exposures without any
     871confounding systematic error sources.  To this end, only
     872photometric nights are selected and all other exposures are ignored.
    871873Each night is allowed to have a single fitted zero point
    872874(corresponding to the sum $zp_{\rm ref} + M_{cal}$ below) and a single
     
    882884\cite{2012ApJ...756..158S} determined flat-field corrections for
    883885$2\times 2$ sub-regions of each chip in the camera and four distinct
    884 time periods (``seasons'').  Later analysis (PV2) used an $8\times8$
     886time periods (``seasons''), ranging from as short as one month to nearly
     88715 months.  Later analysis (PV2) used an $8\times8$
    885888grid of flat-field corrections to good effect.
    886889
     
    10161019corrections are applied to each of the types of measurements stored in
    10171020the database, PSF, Aperture, Kron.  The calibration math remains the
    1018 same regardless of the kind of magnitude being measured.  Also note
    1019 that for the moment, this discussion should only be considered as
    1020 relevant to the chip measurements.  Below we discuss the implications
    1021 for the stack and warp measurements.
     1021same regardless of the kind of magnitude being measured \textmod{(see
     1022  however Section~\ref{sec:phot.stack} for the difference in the stack
     1023  calibration)}.  Also note that for the moment, this discussion
     1024should only be considered as relevant to the chip measurements.  Below
     1025we discuss the implications for the stack and warp measurements.
    10221026
    10231027When the ubercal zero points and flat-field data are loaded,
     
    11051109
    11061110Similarly for images, we exclude those with more than 2 magnitudes of
    1107 extinction or for which the deviation greater of the zero points per
    1108 star are than 0.075 mags or 2$\times$ the median value, whichever is
    1109 greater.  These cuts are somewhat conservative to limit us to only
    1110 good measurements.  The images and stars rejected above are not used
    1111 to calculate the system of zero points and mean magnitudes.  These
    1112 cuts are updated several times as the iterations proceed.  After the
     1111extinction or for which the standard deviation of the zero points are
     1112more than 0.075 mags or 2$\times$ the median value, whichever is greater.
     1113These cuts are somewhat conservative to limit us to only good
     1114measurements.  The images and stars rejected above are not used to
     1115calculate the system of zero points and mean magnitudes.  These cuts
     1116are updated several times as the iterations proceed.  After the
    11131117iterations have completed, the images which have been reject are
    11141118calibrated based on their overlaps with other images.
     
    12331237they wait for the data they need to receive from their neighbors.  The
    12341238management of this stage is performed by communication between the
    1235 region host.  At the end of the iterations, the regions hosts write out
     1239region hosts.  At the end of the iterations, the regions hosts write out
    12361240their final image calibrations.  The master machine then loads the
    12371241full set of image calibrations and then applies these calibrations
     
    12751279data in DVO after the initial relphot calibration to measure the
    12761280flat-field residual with much finer resolution: 124 x 124 flat-field
    1277 values for each GPC1 chip (40x40 pixels per point).  We then used
     1281values for each GPC1 chip (40x40 pixels per point).  \textadd{For this
     1282  analysis, we did not use the entire database, but instead extracted
     1283  relatively bright, but unsaturated measurements (instrumental
     1284  magnitudes between -10.5 and -14.5) for stars with at least 8
     1285  measurements, including 3 used to measure the average photometry in
     1286  the corresponding filter.  These measurements were extracted from a
     1287  collection of 10 sky regions in both low and high-stellar density
     1288  regions covering a total of $\sim 5800$ square degrees of sky.
     1289  Unlike the lower-resolution photometric flat-fields determined in
     1290  the ubercal analysis, the photometric flat-fields calculated in this
     1291  analysis are static in time; they supplement the flats from the
     1292  ubercal analysis.  A total of 1.95 billion measurements were
     1293  extracted for this analysis.
     1294
     1295We then used
    12781296\ippprog{setphot} to apply this new flat-field correction, as well as the
    12791297ubercal flat-field corrections, to the data in the database.  At this
     
    15631581the modified weight is less than 30\% of the median weight
    15641582($\omega^\prime < 0.3 <\omega>$) then the point is treated as clipped.
    1565 The $\chi^2$ is determined from the {\em unclipped} points using the
     1583The $\chi^2$ is determined from the \textadd{remaining} {\em unclipped} points using the
    15661584standard Poisson errors.  Data points which are so excluded are marked
    15671585with bit-flags: \code{ID_MEAS_MASKED_PSF},
     
    15701588
    15711589Bootstrap-resampling analysis is used to assess the errors on the fit
    1572 parameters: A number of measurements equal to the number of {\em
     1590parameters: A number of measurements equal to the number of remaining {\em
    15731591  unclipped} data points are randomly selected from the set of
    15741592unclipped data points, with replacement after each selection.  These
     
    18741892position across the sky.  For each pixel in these images, we selected
    18751893all objects with (14.5, 14.5, 14.5, 14.0, 13.0) $<$ ($g,r,i,z,y$) $<$
    1876 (17, 17, 17, 16.5, 15.5), with at least 3 measurements in $i$-band (to
     1894(17, 17, 17, 16.5, 15.5) magnitudes, with at least 3 measurements in $i$-band (to
    18771895reject artifacts detected in a pair of exposures from the same night),
    18781896with \code{PSF_QF} $> 0.85$ (to reject excessively-masked objects),
     
    30283046community.
    30293047
     3048\note{need to add discussion of SDSS, DES, LSST, Gaia}
     3049
    30303050\acknowledgments
    30313051
Note: See TracChangeset for help on using the changeset viewer.