Tips for sky subtraction¶
Not using default sky subtraction in hscPipe¶
Sky subtraction is carried out in singleFrameDriver.py and coaddDriver.py. The sky subtraction performed in singleFrameDriver.py can be recovered in coaddDriver.py.
# Restore the sky background subtracted in singleFrameDriver.py, then run coaddDriver.py
coaddDriver.py $home/HSC --calib=$home/HSC/CALIB --rerun test --id filter=HSC-I tract=0 --selectId visit=124776..124790:2 ccd=0..103 --config makeCoaddTempExp.bgSubtracted=False
detectCoaddSources.detection.reEstimateBackground=False
# Options
# --config makeCoaddTempExp.bgSubtracted=False:It is set to True by default.
# detectCoaddSources.detection.reEstimateBackground=False : If it is set to True, the background estimation have a problem.
You can check the restored data as ~/HSC/rerun/[rerun]/deepCoadd/[filter]/[tract]/[patch]/warp-[filter]-[tract]-[patch]-[visit].fits. Regarding calexp data, the difference of counts between CCDs sometimes looks clear.
Warning
If you have already done coaddDriver.py and created deepCoadd directory, hscPipe loads warp images under the existing deepCoadd directory when you run the second coaddDriver.py. So you need to rename the [filter] directory name, such as HSC-I to HSC-I_old, then the task generates new warp images which restored sky background.
Changing the size of region used for sky estimation in singleFrameDriver.py¶
You can change the size of region used for sky estimation with the following option.
# Change the size to 256 × 256
singleFrameDriver.py $home/hsc --calib=$home/hsc/CALIB --rerun test --id filter=HSC-I visit=902798..902808:2 --config processCcd.isr.doFringe=False processCcd.detection.background.binSize=256 processCcd.detection.background.binSize=256
# Option
# --config processCcd.detection.background.binSize=256 processCcd.calibrate.background.binSize=256
Running coaddDriver.py without skycorrection.py¶
Run coaddDriver.py with “–config makeCoaddTempExp.doApplySkyCorr=False”
(https://community.lsst.org/t/sky-subtraction/2415)
# Running coaddDriver.py without skycorrection.py
coaddDriver.py $home/hsc --calib=$home/hsc/CALIB --rerun test --id filter=HSC-I visit=902798..902808:2 --config makeCoaddTempExp.doApplySkyCorr=False
# Option
# --config makeCoaddTempExp.doApplySkyCorr=False