Changeset 39845 for trunk/doc/release.2015/ps1.analysis/analysis.tex
- Timestamp:
- Dec 11, 2016, 11:56:48 AM (10 years ago)
- File:
-
- 1 edited
-
trunk/doc/release.2015/ps1.analysis/analysis.tex (modified) (12 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/doc/release.2015/ps1.analysis/analysis.tex
r39839 r39845 410 410 The variance image, if not supplied is constructed by default from the 411 411 flux image using the configuration supplied values of \code{GAIN} and 412 \code{READ \_NOISE} to calculate the appropriate Poisson statistics for412 \code{READ_NOISE} to calculate the appropriate Poisson statistics for 413 413 each pixel. In this case, the image is assumed to represent the 414 414 readout from a single detector, with well-defined gain and read noise … … 444 444 image. The background image and the background standard deviation 445 445 image 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. 447 447 448 448 \subsection{Initial Object Detection} … … 458 458 the covariance, if known. At this stage, the goal is only to detect 459 459 the brighter sources, above a user defined S/N limit (configuration 460 keyword: \code{PEAKS \_NSIGMA\_LIMIT}). A maximum of461 \code{PEAKS \_NMAX} are found at this stage. The detection efficiency460 keyword: \code{PEAKS_NSIGMA_LIMIT}). A maximum of 461 \code{PEAKS_NMAX} are found at this stage. The detection efficiency 462 462 for the brighter sources is not strongly dependent on the form of this 463 463 smoothing function. … … 546 546 {\em key col} for this peak (as used in topographic descriptions of a 547 547 mountain). If the key col for a given peak is less than 548 \code{FOOTPRINT \_CULL\_NSIGMA\_DELTA} (4.0) sigmas below the peak of548 \code{FOOTPRINT_CULL_NSIGMA_DELTA} (4.0) sigmas below the peak of 549 549 interest, the peak is considered to be {\em locally insignificant} and 550 550 removed from the list of possible detections. In the vicinity of a … … 581 581 to find a value of $\sigma_W$ for which $f$ is expected to be 0.65. 582 582 \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 an584 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 are583 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 588 588 measured as defined below. 589 589 … … 615 615 616 616 If the measured centroid coordinates ($x_0, y_0$) differs from the 617 peak coordinates be a large amount (\code{MOMENT \_RADIUS}), then the617 peak coordinates be a large amount (\code{MOMENT_RADIUS}), then the 618 618 peak is identified as being of poor quality and is rejected. In 619 619 both of these cases, it is likely that the `peak' was identified in a … … 638 638 limited at the low and high ends by $R_{\rm min} < M_r < R_{\rm max}$ 639 639 where $R_{\rm min}$ is the first radial moment of the PSF stars, or 640 0.75$\times$ \code{MOMENTS \_GAUSS\_SIGMA} if that cannot be641 determined. $R_{\rm max}$ is set to \code{PSF \_MOMENTS\_RADIUS}, the640 0.75$\times$ \code{MOMENTS_GAUSS_SIGMA} if that cannot be 641 determined. $R_{\rm max}$ is set to \code{PSF_MOMENTS_RADIUS}, the 642 642 size of the moments aperture. 643 643 … … 731 731 registered as part of the model function code. Another function is 732 732 then used to return the appropriate function for a specific model 733 type. For example, the \code{psModelLookup \_GetFunction} will return733 type. For example, the \code{psModelLookup_GetFunction} will return 734 734 the \code{psModelLookup} function for a given model type. This 735 735 mechanism makes it very easy to add new model functions into the … … 756 756 their peaks, as well as an approximate signal-to-noise ratio. All 757 757 objects with a S/N ratio greater than a user-defined parameter 758 (\code{PSF \_SHAPE\_NSIGMA} ???) are selected by PSPhot, though objects758 (\code{PSF_SHAPE_NSIGMA} ???) are selected by PSPhot, though objects 759 759 which have more than a certain number of saturated pixels are excluded 760 760 at this stage. PSPhot then examines the 2-D plane of $\sigma_x, … … 1014 1014 1015 1015 PSPhot will use the user-selected galaxy model to attempt the galaxy 1016 model fits. In the configuration system, the keyword \code{GAL \_MODEL}1016 model fits. In the configuration system, the keyword \code{GAL_MODEL} 1017 1017 is set to the model of interest. All suspected extended objects are 1018 1018 fitted with the model, allowing all of the parameters to float. The … … 1146 1146 value for the ApResid scatter is then used by PSPhot as the best PSF 1147 1147 model for this image. The number of models to be tested is specified 1148 by the configuration keyword \code{PSF \_MODEL\_N}. The configuration1149 variables \code{PSF \_MODEL\_0}, \code{PSF\_MODEL\_1}, through1150 \code{PSF \_MODEL\_N - 1} specify the names of the models which should be1148 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 1151 1151 tested. 1152 1152 … … 1198 1198 1199 1199 The 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}). 1200 annuli, with the radii defined in the configuration 1201 (\code{RADIAL.ANNULAR.BINS.LOWER} \& 1202 \code{RADIAL.ANNULAR.BINS.UPPER}). For each source, the resulting 1203 surface brightness profile is saved in the output cmf-file as an 1204 N-element value in the FITS table (\code{PROF_SB}). The flux at each 1205 radial position and the fill-factor (fraction of pixels used to the 1206 total possible) as also saved as equal-length vectors in the FITS 1207 table (\code{PROF_FLUX} and \code{PROF_FILL}). The values of the 1208 radial bins are saved in the cmf header (\code{RMIN_NN}, 1209 \code{RMAX_NN}). 1209 1210 1210 1211 \note{these profiles are not saved in PSPS}
Note:
See TracChangeset
for help on using the changeset viewer.
