QC documentation system: QC procedure QCsinfo_std.py for SINFONI

This documentation is intended both for QC scientists and SciOps astronomers (who may want to ignore the technical information displayed in grey).
 
HOME
 
CRIRES
FORS1
FORS2
KMOS
NACO
 
GIRAFFE
UVES
XSHOOTER
 
ISAAC
SPHERE
VIMOS
VISIR
 
HAWKI
MUSE
SINFONI
 
AMBER
MIDI
PIONIER
 
OMEGACAM
VIRCAM
[best viewing with style sheets and javascript enabled]
QCsinfo_dark.py QCsinfo_flat.py QCsinfo_wave.py
QCsinfo_dist.py QCsinfo_line.py QCsinfo_pupil.py
QCsinfo_hcao.py QCsinfo_std.py QCsinfo_psf.py
QCsinfo_sci.py
 
bottom  plots   |   contact
NAME QCsinfo_std.py and QCsinfo_persist.py (for the persistence check)
VERSION 1.0 -- 2006-06-22 fill template by handbook
1.1 -- 2006-07-07 revised and certified
2.0 -- 2009-04-20 completely re-written in Python with added QC functions and images
SYNTAX Python
CALL measureQuality
from $DFS_PRODUCT/STD/$DATE:
QCsinfo_std.py -a $AB -i 1
INSTRUMENT SINFONI
RAWTYPE STD
PURPOSE a) generates QC reports
PROCINPUT No input is required:
a) $DATE is now read from AB
b) primary file is set in QCsinfo_std.py and is the PRO.CATG=COADD_STD product with the _0000.fits extension.
c) raw frames are implicitly read from the AB.
d) other products such as MED_COADD_STD, STD_NODDING_STACKED, STD_STAR_SPECTRA, and AO_PERFORMANCE are used and implicitly read.
QC1TABLE trending | table(s) in QC1 database:
sinfoni_std
TRENDPLOT trending | HealthCheck plot(s) associated to this procedure:
trend_report_STD_flux_HC.html
QC1PAGE trending | associated documentation:
std_QC1.html
QC1PLOTS
top

std1.png
display the STD_NODDING_STACKED full detector intermediate product on the Real Time Display
std2.png
display MED_COADD_STD product on the Real Time Display
std3.png
QC Report 1:
UL: the STD_NODDING_STACKED full detector intermediate product image with the number of saturated pixels listed.
UR: the MED_COADD_STD median cube collapsed image. A number of paramters are listed:
the Obs ID, the target name, and for the brightest source detected:
source center (pixels), the integrated flux (within 3 sigma), the isophotal area (3 sigma), the net flux per pixel,
the image FWHM (arcsec), the DIMM seeing corrected for the same airmass and filter used by this standard star,
the status of the AO loop (OPEN or CLOSED),
the position angle of the brightest source (degrees CCW from the x-axis), and its elongation and ellipticity.
The source center, its position angle, and its FWHM are marked on the image. NOTE: The conversion from pixels to arcsec
in the FWHM takes into account the fact that the reconstructed SINFONI cube has RECTANGULAR pixels such that on-sky
Delta_X = 0.5*Delta_Y.
The correct conversion, in terms of the source position angle (PA), is: FWHM(arcsec) = FWHM(pixels) * pixel_scale * sqrt[(0.5*cos(PA))**2 + (sin(PA))**2].
LL: the (scaled) difference between two cube planes. Plane 100 - Plane 2072. The wavelengths of these planes
are listed to the right of this image.
LR: a horizontal (blue line) and vertical (red line) cut through the center of the brightest detected source of the
MED_COADD_STD median cube collapsed image. The respective averages through this center are shown as dotted lines.
std4.png
QC Report 2:
UL: three single column traces of the current product (STD_NODDING_STACKED) through:
x = 500 (dark grey)
x = 1000 (middle grey)
x = 1500 (light grey)
and an averaged sum of all columns (blue).
UR: three single row traces of the current product (STD_NODDING_STACKED) through:
x = 500 (dark grey)
x = 1000 (middle grey)
x = 1500 (light grey)
and an averaged sum of all rows (blue).
LL: the background-subtracted standard star spectrum (STD_STAR_SPECTRA) extracted by the pipeline.
LC: the strehl ratio as a function of wavelength (microns). The median strehl ratio (dotted green line)
and the area denoting the 1 sigma strehl ratio (shaded green area) are shown.
LR: a logarithmic histogram of the STD_NODDING_STACKED full detector image is shown as a dark blue line, while
the reflection of the non-positive flux histogram is shown in light blue.
persist.png
Only if an over-flux is detected is another report created: QC Report 3:
top: a logarithmic histogram of the STD_NODDING_STACKED full detector image with the
persistence critical region shown as a red dotted line.
LL: the STD_NODDING_STACKED full detector intermediate product image.
LR: an estimate of the location of possible persistence critical pixels.
QC1PARAM QC1 parameters written into QC1 table (sinfoni_std):
SOURCE(header/script)       DBNAME                 value          description
QC.FRMON.MEANFLUX qc_frmon_meanflux 0.00 mean flux ON-lamp frame QC.FRMOFF.MEANFLUX qc_frmoff_meanflux 0.00 mean flux OFF-lamp frame QC.FRMDIF.MEANFLUX qc_frmdif_meanflux 16.24 mean flux ON-OFF lamp QC.FRMON.MAXFLUX qc_frmon_maxflux 0.00 max flux in ON-lamp frame QC.FRMOFF.MAXFLUX qc_frmoff_maxflux 0.00 max flux in OFF-lamp frame QC.FRMDIF.MAXFLUX qc_frmdif_maxflux 40646.00 max flux in ON-OFF lamp QC.FRMON.NPIXSAT qc_frmon_npixsat 187 number of ON-lamp saturated pixels QC.CONVFCT qc_convfct None conversion factor QC.CHECK3 qc_check3 0 check evaluation box QC.FWHM.MAJ qc_fwhm_maj 5.20 STD star fwhm along major axis (pixels) QC.FWHM.MIN qc_fwhm_min 4.48 STD star fwhm along minor axis (pixels) QC.FWHM.X qc_fwhm_x 4.67 STD star fwhm along X axis (pixels) QC.FWHM.Y qc_fwhm_y 3.87 STD star fwhm along Y axis (pixels) QC.FWHM.LLX qc_fwhm_llx 8 STD star fwhm LLX QC.FWHM.LLY qc_fwhm_lly 8 STD star fwhm LLY QC.FWHM.HBX qc_fwhm_hbx 16 STD star fwhm HBX QC.THETA qc_theta -0.85 STD star ellipsis angle theta QC.SHIFTX qc_shiftx 9.43 X centroid - image center QC.SHIFTY qc_shifty 0.51 Y centroid - image center QC.STREHL qc_strehl 0.13 strehl ratio from image QCsinfo_std.py qc_strehl_med 0.13 median strehl ratio QCsinfo_std.py qc_strehl_std 0.06 std deviation from all strehl ratios QCsinfo_std.py qc_xpos 45.69 flux-weighted X-axis centroid QCsinfo_std.py qc_ypos 31.58 flux-weighted Y-axis centroid QCsinfo_std.py qc_flux 2.444e+04 total flux (in ADU) within 3 sigma aperture QCsinfo_std.py qc_iso_area 1045 total aperture area over which flux is summed (pixels^2) QCsinfo_std.py qc_net_flux 29.17 net flux (in ADU): qc_flux/qc_iso_area QCsinfo_std.py qc_posang -0.88 STD star position angle (CCW along X-axis) QCsinfo_std.py qc_elong 1.14 STD star elongation (A_image/B_image) QCsinfo_std.py qc_ellip 0.12 STD star ellipticity (1 - B_image/A_image) QCsinfo_std.py qc_fwhm_dimm 0.55 DIMM seeing corrected to observed airmass and IR bandpass (arcsec) QCsinfo_std.py qc_fwhm_sinfo 0.90 STD star FWHM assuming a Gaussian core (arcsec) QCsinfo_std.py qc_delta_fwhm 0.35 difference (qc_fwhm_sinfo - qc_fwhm_dimm) (arcsec) QCsinfo_std.py qc_persist_danger 0 a test if source flux levels are too high (danger of persistence = 1; no danger = 0) QCsinfo_std.py qc_Npersist 0 total number of persistence critical pixels
ALGORITHM Description of algorithms: (status 2009-07-30)
RTC QC parameters: special AO software
QC.FRM parameters : simple statistics on the final product cube
QC.FWHM parameters : Gauss fit integrated into the pipeline; QC1 parameter inventory TBD
X,Y position, flux, and FWHM parameters : the post-pipeline script (QCsinfo_std.py) runs SExtractor (Bertin and Arnouts 1996 A&A 117,393)
and writes these parameters for the brightest source detected. The strehl median and sigma values are first filtered for erroneous values
produced by the pipeline.
DIMM and sinfoni FWHM parameters : The header keywords TEL.AMBI.FWHM.START and TEL.AMBI.FWHM.END store the DIMM seeing measurements
(at 0.5 microns and at the zenith) at the beginning and end of the STD exposure. The average of these two values is computed and converted to the
same airmass and wavelength as the respective telluric standard, and also include a correction for an 8m aperture using the relation given in
(Sarazin, M. 2000, Astronomical Site Monitor, VLT.SPE.ESO.17410.1174) (QCsinfo_std.py).
Persistence test : the post-pipeline script (QCsinfo_persist.py) tests the STD_NODDING_STACKED full detector image for instances of
6 or more pixels (to select against cosmic ray events) having flux levels greater than 10,000 ADU. If 10 or more such events are discovered
an over-flux alarm is given (qc_persist_danger = 1). The number of affected pixels satisfying the above criteria are given by qc_Npersist.
CERTIF Tolerated:
- multiple sources
- stars in sky frame
- strong skew of spectrum
Rejected:
- very low flux levels
- extremely high flux (saturated)
- strong persistence in SKY frame
COMMENTS None
top [2015-10-30T08:24:04]         created by qcDocu v1.1.2, a tqs tool