search

jbwhit / LargeProgramme

This is a place for the results from the UVES Large Programme 185.A-0745 collaboration.
MIT License
0 stars 0 forks source link

data re-reduction #11

Closed jbwhit closed 10 years ago

jbwhit commented 10 years ago

hem1341m1020/546_04_l hem1341m1020/564_10_l hem1341m1020/564_02_l

one transition at 5250 to check.

jbwhit commented 10 years ago

Starting situation: 

# /lustre/projects/p028_swin/jwhitmore/data/UVES/185.A-0745/qso/J134427-103541/hem1341m1020
# sstar001 $ uves_wavcheck.csh
Checking for errors reported in reduction scripts ... none.
390_01_b  1x1  TOL =  0.050  n_av =  14.05  rms =  40.87 m/s
390_02_b  1x1  TOL =  0.045  n_av =  16.00  rms =  38.75 m/s
390_03_b  1x1  TOL =  0.045  n_av =  14.47  rms =  37.72 m/s
390_04_b  1x1  TOL =  0.045  n_av =  13.43  rms =  36.63 m/s
390_05_b  1x1  TOL =  0.045  n_av =  14.10  rms =  37.48 m/s
390_06_b  1x1  TOL =  0.045  n_av =  14.10  rms =  37.60 m/s
390_07_b  1x1  TOL =  0.045  n_av =  13.88  rms =  38.25 m/s
390_08_b  1x1  TOL =  0.050  n_av =  15.20  rms =  40.52 m/s
390_09_b  1x1  TOL =  0.050  n_av =  15.03  rms =  40.92 m/s
390_10_b  1x1  TOL =  0.045  n_av =  13.00  rms =  37.33 m/s
564_01_l  1x1  TOL =  0.055  n_av =  18.50  rms =  29.98 m/s
564_01_u  1x1  TOL =  0.065  n_av =  16.19  rms =  36.59 m/s
564_02_l  1x1  TOL =  0.055  n_av =  18.58  rms =  28.94 m/s
564_02_u  1x1  TOL =  0.060  n_av =  26.56  rms =  33.48 m/s
564_03_l  1x1  TOL =  0.060  n_av =  25.17  rms =  30.06 m/s
564_03_u  1x1  TOL =  0.060  n_av =  27.25  rms =  34.42 m/s
564_04_l  1x1  TOL =  0.055  n_av =  18.38  rms =  29.08 m/s
564_04_u  1x1  TOL =  0.060  n_av =  21.06  rms =  34.42 m/s
564_05_l  1x1  TOL =  0.050  n_av =  22.12  rms =  28.84 m/s
564_05_u  1x1  TOL =  0.060  n_av =  26.69  rms =  32.53 m/s
564_06_l  1x1  TOL =  0.050  n_av =  17.75  rms =  28.23 m/s
564_06_u  1x1  TOL =  0.060  n_av =  24.25  rms =  34.20 m/s
564_07_l  1x1  TOL =  0.055  n_av =  23.83  rms =  30.11 m/s
564_07_u  1x1  TOL =  0.060  n_av =  25.31  rms =  34.33 m/s
564_08_l  1x1  TOL =  0.055  n_av =  24.29  rms =  29.72 m/s
564_08_u  1x1  TOL =  0.065  n_av =  24.75  rms =  34.62 m/s
564_09_l  1x1  TOL =  0.055  n_av =  18.96  rms =  29.62 m/s
564_09_u  1x1  TOL =  0.060  n_av =  26.94  rms =  33.19 m/s
564_10_l  1x1  TOL =  0.050  n_av =  19.08  rms =  28.89 m/s
564_10_u  1x1  TOL =  0.060  n_av =  27.75  rms =  33.16 m/s
jbwhit commented 10 years ago

There is no hem1341m1020/546_04_l -- must refer to 564.

upped-tolerance for lower 564 chip on exposures 564_02, 564_04, 564_10; 

# /lustre/projects/p028_swin/jwhitmore/data/UVES/185.A-0745/qso/J134427-103541/hem1341m1020
# sstar001 $ uves_wavcheck.csh
Checking for errors reported in reduction scripts ... none.
390_01_b  1x1  TOL =  0.050  n_av =  14.05  rms =  40.87 m/s
390_02_b  1x1  TOL =  0.045  n_av =  16.00  rms =  38.75 m/s
390_03_b  1x1  TOL =  0.045  n_av =  14.47  rms =  37.72 m/s
390_04_b  1x1  TOL =  0.045  n_av =  13.43  rms =  36.63 m/s
390_05_b  1x1  TOL =  0.045  n_av =  14.10  rms =  37.48 m/s
390_06_b  1x1  TOL =  0.045  n_av =  14.10  rms =  37.60 m/s
390_07_b  1x1  TOL =  0.045  n_av =  13.88  rms =  38.25 m/s
390_08_b  1x1  TOL =  0.050  n_av =  15.20  rms =  40.52 m/s
390_09_b  1x1  TOL =  0.050  n_av =  15.03  rms =  40.92 m/s
390_10_b  1x1  TOL =  0.045  n_av =  13.00  rms =  37.33 m/s
564_01_l  1x1  TOL =  0.055  n_av =  18.50  rms =  29.98 m/s
564_01_u  1x1  TOL =  0.065  n_av =  16.19  rms =  36.59 m/s
564_02_l  1x1  TOL =  0.065  n_av =  24.29  rms =  31.78 m/s
564_02_u  1x1  TOL =  0.060  n_av =  26.56  rms =  33.48 m/s
564_03_l  1x1  TOL =  0.060  n_av =  25.17  rms =  30.06 m/s
564_03_u  1x1  TOL =  0.060  n_av =  27.25  rms =  34.42 m/s
564_04_l  1x1  TOL =  0.065  n_av =  19.96  rms =  34.35 m/s
564_04_u  1x1  TOL =  0.060  n_av =  21.06  rms =  34.42 m/s
564_05_l  1x1  TOL =  0.050  n_av =  22.12  rms =  28.84 m/s
564_05_u  1x1  TOL =  0.060  n_av =  26.69  rms =  32.53 m/s
564_06_l  1x1  TOL =  0.050  n_av =  17.75  rms =  28.23 m/s
564_06_u  1x1  TOL =  0.060  n_av =  24.25  rms =  34.20 m/s
564_07_l  1x1  TOL =  0.055  n_av =  23.83  rms =  30.11 m/s
564_07_u  1x1  TOL =  0.060  n_av =  25.31  rms =  34.33 m/s
564_08_l  1x1  TOL =  0.055  n_av =  24.29  rms =  29.72 m/s
564_08_u  1x1  TOL =  0.065  n_av =  24.75  rms =  34.62 m/s
564_09_l  1x1  TOL =  0.055  n_av =  18.96  rms =  29.62 m/s
564_09_u  1x1  TOL =  0.060  n_av =  26.94  rms =  33.19 m/s
564_10_l  1x1  TOL =  0.065  n_av =  20.21  rms =  34.84 m/s
564_10_u  1x1  TOL =  0.060  n_av =  27.75  rms =  33.16 m/s
jbwhit commented 10 years ago

Specifying min/max lines seems to be the best -- as done in 564_10_l: 

# /lustre/projects/p028_swin/jwhitmore/data/UVES/185.A-0745/qso/J134427-103541/hem1341m1020
# sstar001 $ uves_wavcheck.csh
Checking for errors reported in reduction scripts ... none.
390_01_b  1x1  TOL =  0.050  n_av =  14.05  rms =  40.87 m/s
390_02_b  1x1  TOL =  0.045  n_av =  16.00  rms =  38.75 m/s
390_03_b  1x1  TOL =  0.045  n_av =  14.47  rms =  37.72 m/s
390_04_b  1x1  TOL =  0.045  n_av =  13.43  rms =  36.63 m/s
390_05_b  1x1  TOL =  0.045  n_av =  14.10  rms =  37.48 m/s
390_06_b  1x1  TOL =  0.045  n_av =  14.10  rms =  37.60 m/s
390_07_b  1x1  TOL =  0.045  n_av =  13.88  rms =  38.25 m/s
390_08_b  1x1  TOL =  0.050  n_av =  15.20  rms =  40.52 m/s
390_09_b  1x1  TOL =  0.050  n_av =  15.03  rms =  40.92 m/s
390_10_b  1x1  TOL =  0.045  n_av =  13.00  rms =  37.33 m/s
564_01_l  1x1  TOL =  0.065  n_av =  19.04  rms =  33.00 m/s
564_01_u  1x1  TOL =  0.065  n_av =  16.19  rms =  36.59 m/s
564_02_l  1x1  TOL =  0.065  n_av =  24.29  rms =  31.78 m/s
564_02_u  1x1  TOL =  0.060  n_av =  26.56  rms =  33.48 m/s
564_03_l  1x1  TOL =  0.060  n_av =  25.17  rms =  30.06 m/s
564_03_u  1x1  TOL =  0.060  n_av =  27.25  rms =  34.42 m/s
564_04_l  1x1  TOL =  0.065  n_av =  19.96  rms =  34.35 m/s
564_04_u  1x1  TOL =  0.060  n_av =  21.06  rms =  34.42 m/s
564_05_l  1x1  TOL =  0.050  n_av =  22.12  rms =  28.84 m/s
564_05_u  1x1  TOL =  0.060  n_av =  26.69  rms =  32.53 m/s
564_06_l  1x1  TOL =  0.065  n_av =  19.17  rms =  33.92 m/s
564_06_u  1x1  TOL =  0.060  n_av =  24.25  rms =  34.20 m/s
564_07_l  1x1  TOL =  0.055  n_av =  23.83  rms =  30.11 m/s
564_07_u  1x1  TOL =  0.060  n_av =  25.31  rms =  34.33 m/s
564_08_l  1x1  TOL =  0.055  n_av =  24.29  rms =  29.72 m/s
564_08_u  1x1  TOL =  0.065  n_av =  24.75  rms =  34.62 m/s
564_09_l  1x1  TOL =  0.065  n_av =  23.08  rms =  34.14 m/s
564_09_u  1x1  TOL =  0.060  n_av =  26.94  rms =  33.19 m/s
564_10_l  1x1  TOL =  0.065  n_av =  31.71  rms =  37.58 m/s
564_10_u  1x1  TOL =  0.060  n_av =  27.75  rms =  33.16 m/s
jbwhit commented 10 years ago

Changing all 564_*_l files via: 

for infile in reduce_564_*.cpl 
do
    awk '{if (/redl/ && /extract/) print $1, $2, $3, $4, $5, $6, $7, $8, $9, "--tolerance=0.065 --minlines=3000 --maxlines=5500", $13; else print}' $infile > temp
    mv temp $infile
done

Rerunning all reductions with the min/max lines changed. Will post results when finished.

jbwhit commented 10 years ago

Final results: 

# /lustre/projects/p028_swin/jwhitmore/data/UVES/185.A-0745/qso/J134427-103541/hem1341m1020
# sstar001 $ uves_wavcheck.csh
Checking for errors reported in reduction scripts ... none.
390_01_b  1x1  TOL =  0.050  n_av =  14.05  rms =  40.87 m/s
390_02_b  1x1  TOL =  0.045  n_av =  16.00  rms =  38.75 m/s
390_03_b  1x1  TOL =  0.045  n_av =  14.47  rms =  37.72 m/s
390_04_b  1x1  TOL =  0.045  n_av =  13.43  rms =  36.63 m/s
390_05_b  1x1  TOL =  0.045  n_av =  14.10  rms =  37.48 m/s
390_06_b  1x1  TOL =  0.045  n_av =  14.10  rms =  37.60 m/s
390_07_b  1x1  TOL =  0.045  n_av =  13.88  rms =  38.25 m/s
390_08_b  1x1  TOL =  0.050  n_av =  15.20  rms =  40.52 m/s
390_09_b  1x1  TOL =  0.050  n_av =  15.03  rms =  40.92 m/s
390_10_b  1x1  TOL =  0.045  n_av =  13.00  rms =  37.33 m/s
564_01_l  1x1  TOL =  0.065  n_av =  20.62  rms =  36.87 m/s
564_01_u  1x1  TOL =  0.065  n_av =  16.19  rms =  36.59 m/s
564_02_l  1x1  TOL =  0.065  n_av =  26.67  rms =  36.59 m/s
564_02_u  1x1  TOL =  0.060  n_av =  26.56  rms =  33.48 m/s
564_03_l  1x1  TOL =  0.065  n_av =  27.42  rms =  37.12 m/s
564_03_u  1x1  TOL =  0.060  n_av =  27.25  rms =  34.42 m/s
564_04_l  1x1  TOL =  0.065  n_av =  26.88  rms =  36.62 m/s
564_04_u  1x1  TOL =  0.060  n_av =  21.06  rms =  34.42 m/s
564_05_l  1x1  TOL =  0.065  n_av =  31.96  rms =  36.69 m/s
564_05_u  1x1  TOL =  0.060  n_av =  26.69  rms =  32.53 m/s
564_06_l  1x1  TOL =  0.065  n_av =  28.00  rms =  37.22 m/s
564_06_u  1x1  TOL =  0.060  n_av =  24.25  rms =  34.20 m/s
564_07_l  1x1  TOL =  0.065  n_av =  32.25  rms =  36.13 m/s
564_07_u  1x1  TOL =  0.060  n_av =  25.31  rms =  34.33 m/s
564_08_l  1x1  TOL =  0.065  n_av =  23.88  rms =  38.16 m/s
564_08_u  1x1  TOL =  0.065  n_av =  24.75  rms =  34.62 m/s
564_09_l  1x1  TOL =  0.065  n_av =  30.50  rms =  36.59 m/s
564_09_u  1x1  TOL =  0.060  n_av =  26.94  rms =  33.19 m/s
564_10_l  1x1  TOL =  0.065  n_av =  31.71  rms =  37.58 m/s
564_10_u  1x1  TOL =  0.060  n_av =  27.75  rms =  33.16 m/s
jbwhit commented 10 years ago

Still have some problems somehow with: 

015_FLUX = fxb_hem1341m1020_sci_564_03_l.fits
025_FLUX = fxb_hem1341m1020_sci_564_08_l.fits

screenshot 2014-04-07 19 04 18

jbwhit commented 10 years ago

Running: 

tcsh
foreach i ( `\ls reduce_564_*.cpl | grep -v _01 | grep -v _06` )
sed 's/--tolerance=0.065 --minlines=3000 --maxlines=5500/--tolerance=0.065 --minlines=3000 --maxlines=4000/' $i > temp.cpl; \mv -f temp.cpl $i
end
exit

to update all (except 1, and 6) exposures to an updated min/max lines.

jbwhit commented 10 years ago

Have to relink names. Couldn't figure out uves_changelinks.csh

cd /lustre/projects/p028_swin/jwhitmore/data/UVES/185.A-0745/qso/J134427-103541/hem1341m1020
newdir="/nfs/cluster/qso/mmurphy/data1/uves_qso_archive/185.A-0745/raw"
for brokenfile in $(find . -type l ! -exec test -e {} \; -print); 
do 
    rootname=$(ls -l $brokenfile | awk '{print $NF}' | awk 'BEGIN{FS="/"}{print $NF}')
    \rm $brokenfile
    ln -s $newdir/$rootname $brokenfile 
done