Read parameters. Derive appropriate input and output filenames. Foreach band: Open all the files: Input image, e.g. mosaic-fovim-350-1100.fits Input Exposure Map Input Mask if applicable Input QPB, SP, SWCX Maps if applicable Determine counts for each input image Endforeach Bin the data (withbinning=T) Adaptively smooth the data (withsmoothing=T) Plot radial profile and size in QDP Open two new output files If withsmoothing=T then Write adapted and size arrays to output files else Write rate and sigma arrays to output files endif Add coord keywords Close output files.