CTIO
Published on CTIO (http://www.ctio.noao.edu/noao)

CTIO Home > Data Reduction Notes > LCO - April 2002

LCO - April 2002

LCO Classical Cam reduction notes April 2002
LCO IR Data Reduction

 

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.


FIX THE HEADERS

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

 

SETUP THE DAOPHOT FILES

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

 

back to top

 
 

Source URL (modified on 05/04/2011 - 10:57): http://www.ctio.noao.edu/noao/content/LCO-April-2002