This documentation is intended both for QC scientists and SciOps astronomers (who may want to ignore the technical information displayed in grey).





bottom
plots

contact


NAME 
QCsinfo_wave.py


VERSION 
1.0  20060621 fill template by handbook
1.1  20060707 certified
2.0  20090101 completely rewritten (from scratch) in Python
4.0  20100301 overhauled with added QC functions and images (for e.g. spectral resolution and arc line position dispersions).


SYNTAX 
Python


CALL 
measureQuality from $DFS_PRODUCT/WAVE/$DATE:
QCsinfo_wave.py a $AB i 1


INSTRUMENT 
SINFONI


RAWTYPE 
WAVE


PURPOSE 
a) generates QC report
b) writes QC1 parameters into local database


PROCINPUT 
No input is required:
a) $DATE is now read from AB
b) primary file is set in QCsinfo_wave.py and is the PRO.CATG=WAVE_MAP product with the _0000.fits extension.
c) The remaining required raw frames are implicitly read from the AB.
d) They include products such as WAVE_MAP, WAVE_LAMP_STACKED, WAVE_COEF_SLIT, WAVE_FIT_PARAMS, SLIT_POS, and RESAMPLED_WAVE are used and implicitly read.
e) Reference WAVE_LAMP_STACKED, RESAMPLED_WAVE, WAVE_COEF_SLIT, and WAVE_MAP are read from $DFO_QC_DIR/references/WAVE.


QC1TABLE 
trending  table(s) in QC1 database:
sinfoni_wave


TRENDPLOT 
trending  HealthCheck plot(s) associated to this procedure:
trend_report_WAVE_length_HC.html


QC1PAGE 
trending  associated documentation:
wave_QC1.html


QC1PLOTS
top 
wave1.png
display WAVE_MAP (the 2d wavelength solution) on the Real Time Display
wave2.png
display WAVE_LAMP_STACKED (the fullframe arc lamp image, uncorrected) on the Real Time Display
wave3.png
display RESAMPLED_WAVE (the fullframe arc lamp image, corrected) on the RealTime Display
wave4.png
QC Report 1:
UL : the current RESAMPLED_WAVE product image.
UC : the reference RESAMPLED_WAVE product image.
UR : the difference between the reference and the current RESAMPLED_WAVE image. The full range of the difference image is
given by the colour bar which, unfortunately, reduces the size of the displayed image (this is a matplotlib bug).
LR : the current WAVE_MAP product image.
LC : the reference WAVE_MAP product image.
UR : the difference between the reference and the current WAVE_MAP image. The full range of the difference image is
given by the colour bar which, unfortunately, reduces the size of the displayed image (this is a matplotlib bug).
wave5.png
QC Report 2:
A four figure plot showing the coefficients computed for the wavelength solution. The coefficients are defined for each row
(yaxis) such that: lambda = C0 + C1*x + C2*x^2 + C3*x^3 (where x = y  (Ny1)/2). The current coefficient values are plotted
for each slitlet and depicted in blue. The reference coefficient values are shown in red.
top : 0th coefficient C0.
top center : 1st coefficient C1.
lower center : 2nd coefficient C2.
bottom : 3rd coefficient C3.
wave6.png
QC Report 3:
Plots of the extracted, calibrated arc line spectrum. The current arc spectrum is shown in blue and the reference spectrum is
plotted in red.
top: The full range arc line spectrum extracted along column pixel 1000. Green dots show arc lines detected by the SINFONI
pipeline. All emission lines contained in the arc line table are labelled.
LL: Zoomed region of the arc line spectrum. The pixel values of the zoom are given in the top left corner of the frame and
the zoom region is shown as a horizontal line in the upper fullscale plot.
LC: Zoomed region of the arc line spectrum. As above.
LR: Zoomed region of the arc line spectrum. As above.
wave7.png
QC Report 4:
top: The spectral resolution (R) as a function of wavelength. The spectral resolution is computed for each arc line
detected by the pipeline. The median spectral resolution is shown as a blue dotted line with blue shading indicating
+/1 standard deviation. The nominal spectral resolution (as given by the SINFONI User's manual pg. 65), for the given
grism and pixel scale, is indicated by a thick green line. Both are written to the top left of the plot.
center: The arc line residuals are computed from the positional differences between the computed arc line positions
and the positions given by the arc line table. The residuals are plotted in units of pixels, and the median value is shown
as a dotted blue line with blue shading indicating +/1 standard deviation.
bottom: The nonlinear wavelength conversion terms. The nonlinear terms are given by:
lambda  lambda(NL) = (C2*x^2 + C3*x^3) (where x = y  (Ny1)/2). The blue line indicates the current frame, while the
red line shows the nonlinear terms of the reference wavelength solution.


QC1PARAM 
QC1 parameters written into QC1 table (sinfoni_wave):
SOURCE(header/script) DBNAME value description
QC.WAVE.ALL qc_wave_all 15 number of identified arc lines
QC.WAVE.NPIXSAT qc_wave_npixsat 0 number of saturated pixels
QC.WAVE.MAXFLUX qc_wave_maxflux 30000 maximum flux in frame
QC.COEF0.AVG qc_coef0_avg 1.90538 zeroth order dispersion coefficient (mean shift over all columns)
QC.COEF0.MED qc_coef0_med 1.95561 zeroth order dispersion coefficient (median shift over all columns)
QC.COEF1.AVG qc_coef1_avg 0.005004 first order dispersion coefficient (mean dispersion over all columns)
QC.COEF1.MED qc_coef1_med 0.005003 first order dispersion coefficient (median dispersion over all columns)
QC.COEF2.AVG qc_coef2_avg 5.7e09 second order dispersion coefficient (mean nonlinear over all columns)
QC.COEF2.MED qc_coef2_med 5.8e09 second order dispersion coefficient (median nonlinear over all columns)
QC.FWHM.AVG qc_fwhm_avg 2.500 average of FWHM of arc lines (in pixels)
QC.FWHM.MED qc_fwhm_med 2.480 median FWHM of arc lines (in pixels)
QCsinfo_wave.py qc_slit17_x 1027.78 xposition of slitlet #17 (QC script) (pixels)
QC.WAVE.POSERR qc_wave_poserr 1.56565e07 the quality of the wavelength solution
QC.FRMON.MEANFLUX qc_frmon_meanflux 68.566 mean flux measured in lampon frame (ADU)
QC.FRMOFF.MEANFLUX qc_frmoff_meanflux 6.79 mean flux measured in lampoff frame (ADU)
QC.FRMDIF.MEANFLUX qc_frmdif_meanflux 54.677 mean flux of lampon frame minus lampoff frame (ADU)
QC.FRMON.MAXFLUX qc_frmon_maxflux 19400.00 maximum flux measured in lampon frame (ADU)
QC.FRMOFF.MAXFLUX qc_frmoff_maxflux 17800.00 maximum flux measured in lampoff frame (ADU)
QC.FRMDIF.MAXFLUX qc_frmdif_maxflux 1234.00 maximum flux measured in lampon minus lampoff frame (ADU)
QC.FRMON.NPIXSAT qc_frmon_npixsat 0 number of saturated pixels measured in lampon frame
QCsinfo_wave.py qc_spec_res_med 6028.28 median spectral resolution over full wavelength range (QC script)
QCsinfo_wave.py qc_spec_res_std 1543.36 std. dev. of spectral resolution over full wavelength range (QC script)
QCsinfo_wave.py qc_arcline_resid_med 0.0125 median residual between detected arc lines and arc line table (QC script) (pixels)
QCsinfo_wave.py qc_arcline_resid_std 0.0025 std. dev. of residuals between detected arc lines and arc line table (QC script) (pixels)


ALGORITHM 
Description of algorithms:
qc_coef0_avg .. qc_coef2_med: For each of the 2048 columns an individual wavelength dispersion solution is found, and specified by four coefficients per column.
The median and the average over the 2048 coefficients is given.
qc_fwhm_avg: The average full width at half maximum of the sample of arc lines (in pixels)
qc_fwhm_med: The median full width at half maximum of the sample of arc lines (in pixels)
qc_slit17_x: The slitlet position value of slitlet #17 taken from the SLIT_POS product table
qc_wave_poserr: The overall quality of the dispersion solution fit
qc_frmon_meanflux .. qc_frmdif_maxflux: A set of simple QC1 parameters comparing the lampon and lampoff frames
qc_frmon_npixsat: The number of saturated pixels
qc_spec_res_med: The spectral resolution is computed for each arc line detected and the median is computed
qc_spec_res_std: The standard deviation of the spectral resolution computed for each arc line detected
qc_arcline_resid_med: The median difference between the detected arc lines and the values in the arc line table
qc_arcline_resid_std: The standard deviation of the difference between the detected arc lines and the values in the arc line table


CERTIF 
Reasons for rejection:
The calibration template requests userinteraction (initialization of the grating, moving it to the desired position) in the case that a
grating timeout occurs in order to avoid a possible grating wheel stuck problem.
If the difference frames (current product minus reference product) show obvious, strong changes in current data.
The RESAMPLED_WAVE product has the wavelength correction applied to the full arc lamp frame. If there are obvious jumps between the horizontal
arc lines (at the junctions of the slitlets) then this warrants investigation.


COMMENTS 

top 
[20131011T06:32:13] created by qcDocu v1.1.2, a tqs tool
