IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Ignore:
Timestamp:
Dec 11, 2016, 11:56:48 AM (10 years ago)
Author:
eugene
Message:

updates to analysis and calibration

File:
1 edited

Legend:

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

    r39839 r39845  
    410410The variance image, if not supplied is constructed by default from the
    411411flux image using the configuration supplied values of \code{GAIN} and
    412 \code{READ\_NOISE} to calculate the appropriate Poisson statistics for
     412\code{READ_NOISE} to calculate the appropriate Poisson statistics for
    413413each pixel.  In this case, the image is assumed to represent the
    414414readout from a single detector, with well-defined gain and read noise
     
    444444image.  The background image and the background standard deviation
    445445image are kept in memory from which the values of \code{SKY} and
    446 \code{SKY\_SIGMA} are calculated for each object in the output catalog.
     446\code{SKY_SIGMA} are calculated for each object in the output catalog.
    447447
    448448\subsection{Initial Object Detection}
     
    458458the covariance, if known. At this stage, the goal is only to detect
    459459the brighter sources, above a user defined S/N limit (configuration
    460 keyword: \code{PEAKS\_NSIGMA\_LIMIT}).  A maximum of
    461 \code{PEAKS\_NMAX} are found at this stage.  The detection efficiency
     460keyword: \code{PEAKS_NSIGMA_LIMIT}).  A maximum of
     461\code{PEAKS_NMAX} are found at this stage.  The detection efficiency
    462462for the brighter sources is not strongly dependent on the form of this
    463463smoothing function.
     
    546546{\em key col} for this peak (as used in topographic descriptions of a
    547547mountain).  If the key col for a given peak is less than
    548 \code{FOOTPRINT\_CULL\_NSIGMA\_DELTA} (4.0) sigmas below the peak of
     548\code{FOOTPRINT_CULL_NSIGMA_DELTA} (4.0) sigmas below the peak of
    549549interest, the peak is considered to be {\em locally insignificant} and
    550550removed from the list of possible detections.  In the vicinity of a
     
    581581to find a value of $\sigma_W$ for which $f$ is expected to be 0.65.
    582582\note{what is the expected ratio of $\sigma_x$ to the true value?}.
    583 We call this value the \code{MOMENTS\_GAUSS\_SIGMA}.  We use an
    584 aperture with a radius of \code{PSF\_MOMENTS\_RADIUS} = 4$\times$
    585 \code{MOMENTS\_GAUSS\_SIGMA} to select the pixels for the measurement.
    586 
    587 Once \code{PSF\_MOMENTS\_SIGMA} has been determined, moments are
     583We call this value the \code{MOMENTS_GAUSS_SIGMA}.  We use an
     584aperture with a radius of \code{PSF_MOMENTS_RADIUS} = 4$\times$
     585\code{MOMENTS_GAUSS_SIGMA} to select the pixels for the measurement.
     586
     587Once \code{PSF_MOMENTS_SIGMA} has been determined, moments are
    588588measured as defined below. 
    589589
     
    615615
    616616If the measured centroid coordinates ($x_0, y_0$) differs from the
    617 peak coordinates be a large amount (\code{MOMENT\_RADIUS}), then the
     617peak coordinates be a large amount (\code{MOMENT_RADIUS}), then the
    618618peak is identified as being of poor quality and is rejected.  In
    619619both of these cases, it is likely that the `peak' was identified in a
     
    638638limited at the low and high ends by $R_{\rm min} < M_r < R_{\rm max}$
    639639where $R_{\rm min}$ is the first radial moment of the PSF stars, or
    640 0.75$\times$ \code{MOMENTS\_GAUSS\_SIGMA} if that cannot be
    641 determined.  $R_{\rm max}$ is set to \code{PSF\_MOMENTS\_RADIUS}, the
     6400.75$\times$ \code{MOMENTS_GAUSS_SIGMA} if that cannot be
     641determined.  $R_{\rm max}$ is set to \code{PSF_MOMENTS_RADIUS}, the
    642642size of the moments aperture.
    643643
     
    731731registered as part of the model function code.  Another function is
    732732then used to return the appropriate function for a specific model
    733 type.  For example, the \code{psModelLookup\_GetFunction} will return
     733type.  For example, the \code{psModelLookup_GetFunction} will return
    734734the \code{psModelLookup} function for a given model type.  This
    735735mechanism makes it very easy to add new model functions into the
     
    756756their peaks, as well as an approximate signal-to-noise ratio.  All
    757757objects with a S/N ratio greater than a user-defined parameter
    758 (\code{PSF\_SHAPE\_NSIGMA} ???) are selected by PSPhot, though objects
     758(\code{PSF_SHAPE_NSIGMA} ???) are selected by PSPhot, though objects
    759759which have more than a certain number of saturated pixels are excluded
    760760at this stage.  PSPhot then examines the 2-D plane of $\sigma_x,
     
    10141014
    10151015PSPhot will use the user-selected galaxy model to attempt the galaxy
    1016 model fits.  In the configuration system, the keyword \code{GAL\_MODEL}
     1016model fits.  In the configuration system, the keyword \code{GAL_MODEL}
    10171017is set to the model of interest.  All suspected extended objects are
    10181018fitted with the model, allowing all of the parameters to float.  The
     
    11461146value for the ApResid scatter is then used by PSPhot as the best PSF
    11471147model for this image.  The number of models to be tested is specified
    1148 by the configuration keyword \code{PSF\_MODEL\_N}.  The configuration
    1149 variables \code{PSF\_MODEL\_0}, \code{PSF\_MODEL\_1}, through
    1150 \code{PSF\_MODEL\_N - 1} specify the names of the models which should be
     1148by the configuration keyword \code{PSF_MODEL_N}.  The configuration
     1149variables \code{PSF_MODEL_0}, \code{PSF_MODEL_1}, through
     1150\code{PSF_MODEL_N - 1} specify the names of the models which should be
    11511151tested.
    11521152
     
    11981198
    11991199The surface brightness values are sampled at a number of radial
    1200 annuli, with the radii defined in the configuration ({\tt
    1201   RADIAL.ANNULAR.BINS.LOWER \& RADIAL.ANNULAR.BINS.UPPER}).  For each
    1202 source, the resulting surface brightness profile is saved in the
    1203 output cmf-file as an N-element value in the FITS table ({\tt
    1204   PROF\_SB}).  The flux at each radial position and the fill-factor
    1205 (fraction of pixels used to the total possible) as also saved as
    1206 equal-length vectors in the FITS table ({\tt PROF\_FLUX and
    1207   PROF\_FILL}).  The values of the radial bins are saved in the cmf
    1208 header ({\tt RMIN\_NN, RMAX\_NN}).
     1200annuli, with the radii defined in the configuration
     1201(\code{RADIAL.ANNULAR.BINS.LOWER} \&
     1202\code{RADIAL.ANNULAR.BINS.UPPER}).  For each source, the resulting
     1203surface brightness profile is saved in the output cmf-file as an
     1204N-element value in the FITS table (\code{PROF_SB}).  The flux at each
     1205radial position and the fill-factor (fraction of pixels used to the
     1206total possible) as also saved as equal-length vectors in the FITS
     1207table (\code{PROF_FLUX} and \code{PROF_FILL}).  The values of the
     1208radial bins are saved in the cmf header (\code{RMIN_NN},
     1209\code{RMAX_NN}).
    12091210
    12101211\note{these profiles are not saved in PSPS}
Note: See TracChangeset for help on using the changeset viewer.