Classic Cam
0.60"/pix
7.8e-/ADU
40e- r.o.n.'
12500 ADU full well (100000 e-)
Recommends keeping below 10000ADU
Mark has reduced the IR data through flat-fielding and sky subtraction.
For the program images, he has combined the data.
cp /uw50/nick/daophot/irafcl/lco/setuplco .
cp /uw50/nick/daophot/irafcl/lco/*.cl .
0. After dumping the data, clean up all the useless files. Create
mkdir crap
mkdir data
and run
cl < crap.cl
cl < data.cl
1. Edit in vital information. First:
hsel *.fits $I,date-obs,ut yes
Run the task
task lcohead = home$nickcl/lcohead.cl
epar lcohead
lcohead is*.fits
Enter the correct DATE-OBS as DD/MM/YY (if 1900-1999) or YYYY-MM-DD if 2000+ (this also works for all dates) into the param file. Note that this task is not smart enough to realize that the date changed during the night. The DATE-OBS must be corrected if the observations start at 23h UT.
Mark has named the directories as the UT date (ie 1/2 Nov 1999 is 2nov99).
Note that the DATE_OBS and the UT may or may not contain extraneous blanks which may screw upt the program. On the 40", the UT had blanks between hh mm ss. On magellan it did not. You will have to modify the lcohead program. The code is there - it just has to be commented out.
2. Correct any DATE_OBS now.
hsel is*.fits $I,ut,date-obs yes > in1
emacs in1 (just the image names for the date changes)
hedit @in1 date-obs "01/11/99" up+ ver-
hsel is*.fits $I,ut,date-obs yes
3. Run
setjd i*.imh hjd="" epoch="EQUINOX"
Can't run setairmass because there is no ST.
4. Copy all the individual SN image to the data directory to clean up the working directory.
hsel is*.fits $I,title yes | sort - col=2 > in2
edit in2 and copy the indivdual SN data to data/
imren @in2 data
5. translate fits to imh
task cpimh = home$nickcl/cpimh.cl
cpimh *.fits delin+
6. Copy in the missing information into the SN mosaicked data.
files SN*.imh > in3
task lcomosaic = home$nickcl/lcomosaic.cl
Edit in3 as:
lcom SN1999ebHn_1656_1675_02nov99 data/is1656
lcom SN1999ebHn_1656_1675_02nov99_med data/is1656
lcom SN1999ebXn_1726_1745_02nov99 data/is1726
etc.
7. Copy the SN*med data to data/
imrename *med.imh data/
8. Rename the SN*.imh data to something shorter.
MASK IMAGE
We need to identify the bad pixels.
mask1.cl:
imdel mask.imh,mask.pl
imcomb is*.imh mask.imh comb=median reject- scale=exposure zero=median
imhist mask.imh z1=-10 z2=10
Look at mask.imh and get low and high values
mask2.cl:
#
# good pix are 0, bad are 1 for IRAF mask
# the values -2 and 6 need to be checked on the histogram
# each time you make the mask.
#
real lll, uuu
lll = -1
uuu = 5
displ mask.imh 1
imcopy mask.imh temp.imh
imrep("mask.imh", lower=lll, upper=uuu, val=0)
imrep("mask.imh", lower=INDEF, upper=lll, val=-1 )
imrep("mask.imh", lower=uuu, upper=INDEF, val=-1)
imar mask.imh * mask.imh mask.imh
imcopy mask.imh mask.pl
# make DAOPHOT mask where bad pix are 0 and good are 1
imrename mask.imh maskdao
imar maskdao - 1 maskdao
imar maskdao * -1 maskdao
hedit maskdao.imh title "mask" up+ ver-
hedit mask.pl title "mask" up+ ver-
display mask.pl 1 zs- zr- z1=0 z2=1
display maskdao 1 zs- zr- z1=0 z2=1
imar is*.imh / maskdao is*.imh divzero=32766 calc=real pix=short
1. Copy daophot.opt,allstar.opt,photo.opt,jhk_short.lib,filters.dat to the directory.
cp ../*.opt .
cp ../jhk_short.lib .
cp ../filters.dat .
cp ../jhk.tfm novxx.tfm
cp ../*.cl .
daophot.opt:
| Read noise | = | 5.1 | |
| Gain | = | 7.8 | |
| FWHM | = | 2.0 | |
| Fitting radius | = | 2.5 | |
| PSF radius | = | 12 | |
| Analytic model PSF | = | -3 | |
| Variable PSF | = | 0 | |
| Extra PSF cleaning passes | = | 5 | |
| High good datum | = | 10000 | |
| Watch progess | = | -2 | |
| Threshold | = | 7 |
allstar.opt:
| Fitting radius | = | 3.3 | |
| IS (Inner sky radius) | = | 2 | |
| OS (Outer sky radius) | = | 16 | |
| Redetermine Centroids | = | 1 |
photo.opt:
| A1 = | 3.0000 | |
| A2 = | 3.1429 | |
| A3 = | 3.3571 | |
| A4 = | 3.6429 | |
| A5 = | 4.0000 | |
| A6 = | 4.4286 | |
| A7 = | 4.9286 | |
| A8 = | 5.5000 | |
| A9 = | 6.1429 | |
| AA = | 6.8571 | |
| AB = | 7.6429 | |
| AC = | 8.5000 | |
| IS = | 8.5000 | |
| OS = | 17.0000 |
2. Make the *.inf file
del in*
files *%.imh%% > in1
hsel @in1 $I,filter,ut,airmass,exptime,jd,title yes > junk.dat
!$myprog/prog3lco junk.dat
filters.dat:
'Jshort' 1
'H' 2
'Kshort' 3
'X' 4
3. Run BFIND on the is*.imh files. Seems to run without problems or the need to input FWHM. Use threshold of 25.
4. For the program fields, first get the FWHM by running:
loadit SNxxx dis- ; imexam keep+
etc.
$myprog/prog39 junk.dat
Edit fwhm.dat. Most of the data will use "0" variance.
Then run BPASS2 on the data. I used a threshold of 10
5. Check the psf's and run final2. Check data for new stars. Run:
loadit SN1999emH dis+ ; addals SN1999emH
loadit SN1999emj dis+ ; addals SN1999emj
etc.
5.5 Edit the *.ap files if they are corrupted by a bright galaxy or star.
loadit SN1999emH ; mark SN1999emH.kap 1 ; emacs SN1999emH.kap
etc.
I had trouble with SN1999em and p9181
6. Run DAOGROW on the data. I used 0.025mag error and 4 variables.
This worked okay only after I got rid of the bad pixels (with the mask). If I included SN data, I used up to 0.05.
In some cases, the brightest star was below 10000ADU but was clearly saturated. This threw off DAOGROW, and the some of the curves were crazy. If you look at each curve quickly you can find the bad ones. Alternatively, you can do:
grep Ro nov06.gro | sort +4
and the bad ones will have *huge* Ro values. You have to edit out the *.ap files for the saturated stars - typically there will be a few ap values on the star.
7. DAOMATCH/DAOMASTER
For simple shifts, you can run
hsel is*.imh $I,title yes | sort - col=2 | grep 9115 | sort > in9115
hsel is*.imh $I,title yes | sort - col=2 | grep 9181 | sort > in9181
hsel is*.imh $I,title yes | sort - col=2 | grep 9106 | sort > in9106
hsel is*.imh $I,title yes | sort - col=2 | grep 9138 | sort > in9138
hsel is*.imh $I,title yes | sort - col=2 | grep 9136 | sort > in9136
hsel is*.imh $I,title yes | sort - col=2 | grep 9137 | sort > in9137
hsel is*.imh $I,title yes | sort - col=2 | grep 9155 | sort > in9155
hsel is*.imh $I,title yes | sort - col=2 | grep 9157 | sort > in9157
hsel is*.imh $I,title yes | sort - col=2 | grep 9164 | sort > in9164
hsel is*.imh $I,title yes | sort - col=2 | grep 99em | sort > in99em
task irtest = home$nickcl/irtest.cl
Now mark the star.
irtest @inx usedither-
Check for bad centers.
sort junk.dat col=5
!$myprog/prog52 junk.dat
head temp.mch nlines=1
mv temp.mch
Copy temp.mch to another .mch file.
Run DAOMASTER on all the data. I used a simple translation only.
8. Run fetch to get the *.fet files.
Note that for SN1999ei, there were two different positions, and not all the stars are on the second position. To make the *.fet file, you have to add the stars that are missing from the first frame. You can do this two ways.
a. find the star in the second frame:
2 145.132 381.642 15.879 0.0030 6.653 4. 0.129 -0.003
get the *.mch file:
' SN1999eiH1.als ' 0.000 0.000 1.00000 0.00000 0.00
000 1.00000 0.000 0.0331
' SN1999eiH2.als ' -15.432 65.735 1.00000 0.00000 0.00
000 1.00000 0.239 0.1024
Add the offsets for the second frame to get to the master frame:
145.132 + -15.432
381.642 + 65.735
and put these in the *.fet file.
b. Run daomaster and output the *.tfr table. Find the star in the table:
200005 129.700 447.376 0 1
and edit the *.fet file.
9. Run collect to make the *.obs file.
$myprog/prog43
10. Run CCDSTD and CCDAVE