IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Changeset 41402 for trunk


Ignore:
Timestamp:
Aug 15, 2020, 10:43:08 AM (6 years ago)
Author:
eugene
Message:

updates to the analysis paper

Location:
trunk/doc/release.2015/ps1.analysis
Files:
5 edited

Legend:

Unmodified
Added
Removed
  • trunk/doc/release.2015/ps1.analysis/Makefile

    r41324 r41402  
    44#
    55DO_PDFLATEX = 1
    6 DO_BIBTEX = 1
     6DO_BIBTEX = 0
    77
    88help:
     
    2424PDFPICS = \
    2525pics/peaks.pdf \
    26 pics/galplanecut.pdf \
    2726pics/FWHM.smooth.trend.v1.ps1.pdf \
     27pics/iq.exposures.pdf \
    2828pics/radial.profiles.pdf \
    2929pics/moment.class.pdf \
    30 pics/mag.resid.psf.pdf \
    31 pics/mag.resid.aper.pdf \
     30pics/completion.ppsim.pdf \
     31pics/psphot.complete.pv3.pdf \
     32pics/mag.resid.psf.v1.pdf \
     33pics/mag.resid.aper.v1.pdf \
     34pics/bright.mag.resid.pdf \
     35pics/zpt.mjd.v0.i.pdf \
     36pics/zptres.hist.v0.i.pdf \
     37pics/galplanecut.pdf \
     38pics/petrosians.mags.pdf \
    3239pics/galaxy.exp.complete.pdf \
    3340pics/galaxy.dev.complete.pdf \
    3441pics/galaxy.exp.params.pdf \
    35 pics/galaxy.dev.params.pdf
     42pics/galaxy.dev.params.pdf \
     43pics/compare.mags.pdf
    3644
    3745PNGPICS = \
  • trunk/doc/release.2015/ps1.analysis/analysis.tex

    r41347 r41402  
    2929
    3030%\def\picdir{/home/eugene/chipresid.20140404}
    31 \def\picdir{pics}
    32 %\def\picdir{.}
     31%\def\picdir{pics}
     32\def\picdir{.}
    3333
    3434% Pick a terse version of the title here;
     
    184184  the image data products and a hierachical database of measurements
    185185  using a system developed specifically for the Pan-STARRS dataset.
    186   Development of this database systems was the product of a
     186  Development of this database system was the product of a
    187187  collaboration between the Pan-STARRS Project and Alex Szalay's
    188188  database development group at The Johns Hopkins University (JHU)
     
    375375  is performed in parallel on each of the individual CCDs in the
    376376  camera.  This so-called \ippstage{chip} stage analysis includes the
    377   detrending of CCD image as well as the detection and analysis of
     377  detrending of CCD image (instrumental signature removal) as well as the detection and analysis of
    378378  sources in the image using the basic version of \ippprog{psphot}.
    379379  The next stage of the analysis, the \ippstage{camera} stage,
     
    388388  a large regular pixel grid is defined, and then subdivided along
    389389  pixel boundaries into smaller units which are well-matched to the
    390   memory footprint of our processing computesr.  These smaller images,
     390  memory footprint of our processing computers.  These smaller images,
    391391  called `skycells' are defined with 1 arcminute of overlap with their
    392392  neighbors to that any modest-sized object can be analysed entirely
     
    525525  and astrometry accuracy at the level of our goals, not only must the
    526526  measurement of the astronomical detections be precise, but it is
    527   necessary for the detrending (instrumental signature remove) and
     527  necessary for the detrending and
    528528  calibration processes to correct for a wide variety of systematic
    529529  effects and it is also necessary for the observations to be
    530530  performed in such a way that the data can be calibrated well.  These
    531   others aspects of the process are discussed in detail elsewhere
     531  other aspects of the process are discussed in detail elsewhere
    532532  (Papers I, III, V).  In the end, the goals were largely achieved for
    533533  the Pan-STARRS\,1 $3\pi$ survey. As reported in Paper V, the
     
    12271227  classification of the sources.  As discussed below, the second
    12281228  moments are used to select candidate stellar sources to be used in
    1229   modeling the PSF and the exclude `cosmic rays' and extended sources.
     1229  modeling the PSF and to identify `cosmic rays' and extended sources.
    12301230  The radial moment is used in the measurement of the Kron magnitudes \citep{1980ApJS...43..305K}.
    12311231  The higher-order moments are provided primarily for image quality
     
    17361736radial moment as the major axis size for the Gaussian ($\sigma_a$), retaining
    17371737the position angle and axial ratio from the calculation above.  We use
    1738 these guess parameters for all version of the PSF analytical models,
     1738these guess parameters for all versions of the PSF analytical models,
    17391739despite the fact that for the versions which are not approximations of
    17401740Gaussians these guess values will be systematically incorrect. 
     
    21462146  source model parameters (position in $X$ and $Y$ and flux
    21472147  normalization) are allowed to vary in the fit.  Note that we do {\em
    2148     not} allow the local sky to be fitted as a free parameters.  Since
     2148    not} allow the local sky to be fitted as a free parameter.  Since
    21492149  we have subtracted a model for the background, allowing the sky to
    2150   be again at this stage is redundant.  In fact, in our testing, we
    2151   found that allowing the sky to float resulted in higher scatter for
     2150  be fitted again at this stage is redundant.  In fact, in our testing, we
     2151  found that allowing the sky background value to float resulted in higher scatter for
    21522152  the flux normalizations.  For the non-linear fitting,
    21532153\ippprog{psphot} again uses the Levenberg-Marquardt technique.}  The
     
    26222622crowding and confusion.  Since the injection and recovery analysis of
    26232623the fake sources operates on the source-subtracted image and does not
    2624 attempt to fully discovery the sources, this analysis over-estimates
     2624attempt to fully discover the sources, this analysis over-estimates
    26252625the completeness in crowded fields.  To explore the completeness in
    26262626crowded field images, we generate a series of simulated images using a
     
    34803480least the smallest 4 apertures.  Sources for which photometry in these
    34813481fixed aperture are calculated have the flag bit
    3482 \code{PM_SOURCE_MODE_RADIAL_FLUX} set.  \textadd{Although these aperture are
     3482\code{PM_SOURCE_MODE_RADIAL_FLUX} set.  \textadd{Although these apertures are
    34833483chosen to match the SDSS apertures, the SDSS images are measured on
    34843484unconvolved images.  Since the median seeing for the SDSS images is
     
    37693769 \includegraphics[width=\hsize,clip]{\picdir/{compare.mags}.pdf}
    37703770  \caption{\label{fig:compare.mags} Comparison of {\tt psphot} average
    3771     chip photometry, average forced-warp photometry, and stack
    3772     photometry from $3\pi$ Survey data to average forced-warp
     3771    chip photometry (panel a), average forced-warp photometry (panel b), and stack
     3772    photometry (panel c) from $3\pi$ Survey data to average forced-warp
    37733773    photometry from the Pan-STARRS\,1 Medium-Deep Survey field MD06
    37743774    At bright magnitudes, average chip photometry is the most
     
    37953795Paper V).}
    37963796
    3797 {\TEXTADD As can be clearly seen in the figure, the average from the forced-warp
    3798 photometry is slightly worse than the chip photometry, while the stack
    3799 PSF photometry is significantly degraded.  We attribute the latter
    3800 effect to the highly-textured PSF observed in the stack images due to
    3801 the combination of variable PSFs in each exposure and significant
    3802 masking fraction in the PS1 camera.  At the faint end, the chip
    3803 photometry is significantly worse that both average warp and stack
    3804 photometry.  First, in order to have a measurement, a source must be
    3805 detected above the detection threshold in at least one of the
    3806 exposures, limiting the depth possible of the average chip
    3807 photometry. Second, at the faint end, only bright fluctuations will be
    3808 detected, resulting in a bright bias. This latter effect is clearly
    3809 seen in Figure~\ref{fig:compare.mags} as the average chip magnitudes
    3810 diverge from the deeper Medium Deep photometry measurements.  As has
    3811 been noted elsewhere \citep{2018ApJS..234....1B}, the warp and stack
    3812 photometry is also degraded for objects which have significant proper
    3813 motion over the course of the $3\pi$ Survey since the position is held
    3814 constant for all epochs, while the average chip photometry is
    3815 calculated on detections which are cross-matched in the database.
    3816 Thus, warp and stack photometry should be avoided for sources with
    3817 proper motion greater than roughly 100 milliarcseconds per year.}
     3797{\TEXTADD As can be clearly seen in the figure, the average from the
     3798  forced-warp photometry is slightly worse than the chip photometry,
     3799  while the stack PSF photometry is significantly degraded.  We
     3800  attribute the latter effect to the highly-textured PSF observed in
     3801  the stack images due to the combination of variable PSFs in each
     3802  exposure and significant masking fraction in the PS1 camera.  At the
     3803  faint end, the chip photometry is significantly worse that both
     3804  average warp and stack photometry.  First, in order to have a
     3805  measurement, a source must be detected above the detection threshold
     3806  in at least one of the exposures, limiting the depth possible of the
     3807  average chip photometry. Second, at the faint end, only bright
     3808  fluctuations will be detected, resulting in a bright bias, a form of
     3809  Eddington bias \citep{1913MNRAS..73..359E}. This latter effect is
     3810  clearly seen in Figure~\ref{fig:compare.mags} as the average chip
     3811  magnitudes diverge from the deeper Medium Deep photometry
     3812  measurements.  As has been noted elsewhere
     3813  \citep{2018ApJS..234....1B}, the warp and stack photometry is also
     3814  degraded for objects which have significant proper motion over the
     3815  course of the $3\pi$ Survey since the position is held constant for
     3816  all epochs, while the average chip photometry is calculated on
     3817  detections which are cross-matched in the database.  Thus, warp and
     3818  stack photometry should be avoided for sources with proper motion
     3819  greater than roughly 100 milliarcseconds per year.}
    38183820
    38193821\subsection{Forced Galaxy Models}
     
    40894091\cite{2008ApJ...677..808Y}.  The analysis of the sources detected in
    40904092these difference images uses a portion of the \ippprog{psphot} code
    4091 embedded in the program, \ippprog{ppSub}, which generates those image.
     4093embedded in the program, \ippprog{ppSub}, which generates those images.
    40924094Difference images are generated from three different possible image
    40934095combinations: 1) pairs of individual exposures are differenced using
     
    42874289
    42884290\bibliographystyle{apj}
    4289 \bibliography{lib}{}
    4290 %\input{analysis.bbl}
     4291%\bibliography{lib}{}
     4292\input{analysis.bbl}
    42914293
    42924294\end{document}
  • trunk/doc/release.2015/ps1.analysis/response.txt

    r41347 r41402  
     1
     2We thank the referee for the many detailed suggestions.  We have
     3incorportated these in almost all cases into the article.  The
     4result is greatly improved in terms of clarity and context.  Below are
     5our specific responses the the suggestions of the referee.
     6
     7(Our responses to the comments start with "**" and are indented).
    18
    29Referee Report
     
    6774mask values appear in the DR1 and DR2 data releases. For example,
    6875Tables 1-4 must refer to some named quantities in the data releases.
    69 This paper should also clearly indicate which quantities are recommended
    70 to be most reliable for point source astrometry, fluxes and colors.
    71 (My guess would be the astrometry from the stacked images in Sec 4.7, since
    72 those aren't recomputed; fluxes from the averaged forced photometry
    73 in Sec 6; and colors from either aperture photometry on the stacked
    74 images or from the averaged forced photometry, and the authors must
    75 know which is demonstrated to be more reliable). The paper should
    76 state the same for galaxy astrometry, fluxes and colors.
     76This paper should also clearly indicate which quantities are
     77recommended to be most reliable for point source astrometry, fluxes
     78and colors.  (My guess would be the astrometry from the stacked images
     79in Sec 4.7, since those aren't recomputed; fluxes from the averaged
     80forced photometry in Sec 6; and colors from either aperture photometry
     81on the stacked images or from the averaged forced photometry, and the
     82authors must know which is demonstrated to be more reliable). The
     83paper should state the same for galaxy astrometry, fluxes and colors.
    7784
    7885** for each section, we have added a summary of where the values may
     
    97104that the photometric goals are achieved
    98105
    99 ** added comparion discussion of chip, warp, stack photometry at the end of Sec 6.1
     106** added comparion discussion of chip, warp, stack photometry at the
     107   end of Sec 6.1
    100108
    101109- Sec 7, where the image differencing detections and photometry is used
     
    643651   are completely constrained.
    644652
    645 We have renamed Section 6 as Forced Warp Analysis and split out the
    646 PSF vs extended source analysis sections as recommended.  We added
    647 some explanation at the end of the section to explain what we mean by
    648 'forced'.
     653** We have renamed Section 6 as Forced Warp Analysis and split out the
     654   PSF vs extended source analysis sections as recommended.  We added
     655   some explanation at the end of the section to explain what we mean
     656   by 'forced'.
    649657
    650658- The general description of the section should end with "variant of psphot",
     
    665673- The terms "skycell" and "warp image" are first used here without
    666674definition. Are warp images the same as CAMERA and CHIP?
     675
    667676** Updated Section 2 to outline the relevant processing stages and
    668677define 'warp', 'skycell', and 'stack' more cleanly.
Note: See TracChangeset for help on using the changeset viewer.