Mscstack takes the geometrically transformed images, and "stacks" them into a single image per filter, using the values of mscscale and msczero that have been entered into the image headers by "philmatch" as intensity scaling and offset values.
PACKAGE = mscred TASK = mscstack input = @us List of images to combine output = N6822U.fits Output image (headers= ) List of header files (optional) (bpmasks= ) List of bad pixel masks (optional) (rejmask= ) List of rejection masks (optional) (nrejmas= ) List of number rejected masks (optional) (expmask= ) List of exposure masks (optional) (sigmas = ) List of sigma images (optional) (combine= average) Type of combine operation (median|average) (reject = ccdclip) Type of rejection (masktyp= goodvalue) Mask type (maskval= 0.) Mask value (blank = 0.) Value if there are no pixels (scale = !mscscale) Image scaling (zero = !msczero) Image zero point offset (weight = none) Image weights (statsec= ) Image section for computing statistics (lthresh= 1.) Lower threshold (hthresh= INDEF) Upper threshold (nlow = 1) minmax: Number of low pixels to reject (nhigh = 1) minmax: Number of high pixels to reject (nkeep = 1) Minimum to keep (pos) or maximum to reject (neg) (mclip = yes) Use median in sigma clipping algorithms? (lsigma = 3.) Lower sigma clipping factor (hsigma = 3.) Upper sigma clipping factor (rdnoise= 10) ccdclip: CCD readout noise (electrons) (gain = 2.2) ccdclip: CCD gain (electrons/DN) (snoise = 0) ccdclip: Sensitivity noise (fraction) (sigscal= 0.1) Tolerance for sigma clipping scaling corrections (pclip = -0.5) pclip: Percentile clipping parameter (grow = 0.) Radius (pixels) for neighbor rejection (mode = ql)