HOW TO USE ZZ.CL FOR REDUCING HARTMANN SCREEN DATA
Revision of J.Baldwin's notes from 6July98
All the updated programs are stored in /ua76/boccas/4m/to_run_shap/
SETUP:
LOAD IRAF TASKS:
- images
- proto
- lists
- digiphot
- apphot
GET THE FOLLOWING EXECUTABLES AND FILES into the directory you are
using:
- shap
- shapav9
- newshap
- map4
- setup.cl
- zz.cl
- modcheck.cl
- ann_zer.ann
- ela_ann.dat
plus a .coo file appropriate for the telescope/f-ratio you are using:
- 36 INCH: 36inch.coo
- 4m pf and f/8: model.coo
- 4m f/14: f14.coo
- 4m f/30: f30.coo
THEN:
- task setup=setup.cl
- task modcheck=modcheck.cl
- task zz=zz.cl
- task $shap=$shap9
PLUS A COUPLE OF HANDY ALIASES
- alias te 'textedit \!* &'
- alias lpflip 'enscript -r -l -fCourier9 -Pnpt0x2 \!*'
SET TELESCOPE AND CCD PARAMETERS in setup:
epar setup
TASK = setup
image = obj134 Image:
coofile = f14.coo coordinate file:
(pixsize= 0.048) CCD pixel size (mm)
(apfill = 1.0733) Hart screen fill factor, 4m=1.0733, 36in=1.0881
(mirdia = 3947.7) Mir dia (mm), 4m=3947.7, 36in=914
(fratio = 14.45) 4m=14.45,7.8,2.904; 36in=13.5
(cursor = )
(coof2 = f14.coo)
(mode = al)
TAKE A FIRST TEST FRAME as early as possible in the night.
RUN setup. Mark a pair of spots in the top and bottom rows of the
Hartmann pattern... their separation is used by Shap so you cannot just
pick any spots anywhere in the pattern. Note screwy feature where you
must mark 1st spot with cursor + any key, then get cursor back to IRAF
window and hit RETURN in order to get prompt for 1st spot number. Setup
will produce files shap.par, modshap.coo, and zz.input
setup obj134 f14.coo
Use textedit or any other editor to remove bad spots from modshap.coo
Run modcheck to be sure you have everything set up correctly.
modcheck obj134
NOW EPAR A FEW VITAL PARAMETERS FOR phot:
centerp:
cbox should be set to pacing between points
maxshift same
fitskypars:
annulus tune to dot size
dannulus tune to spacing between dots
photpars:
aperture set to same value as annulus
CHECK SHAP CALIBRATIONS, ANGLE ZERO POINTS, ETC:
- determine directions on sky using offset command.
- take shap data while forcing in in 4 mu astig at PA 0, PA 90.
MAKE SKY MAPS WITH CORRECTIONS OFF.
NOW RUN zz ON EACH IMAGE. 3 exposures 30 sec each.
- For each new object or telescope configuration:
'newshap'
(enter title. If mapping sky, enter ha, dec in hh.h dd.d format;
eg. degrees and fraction of degree, rather than deg.min)
- For each image:
zz objxxx
- as needed:
shapav9,
tail shap.work,
tail shap.sumry
where shapav9 averages together all data between successive applications of newshap and writes detailed log file in shap.sumry, averages in shap.av
TO PLOT RESULTS OF SKY MAP: use texteditor to copy relevant sections of shap.av into temporary file.
map4
Maxime Boccas, 4Feb01