wave

Wavelength Calibration

wave.py

apogee.apred.wave.allplots()

Routine to put together master web page for summary plots from all years

apogee.apred.wave.ditherplots(planfile, vers=None, inst=None)

Make HTML file for dither plots for a visit

apogee.apred.wave.findlines(frame, rows, waves, lines, out=None, verbose=False, estsig=2)

Determine positions of lines from input file in input frame for specified rows

Parameters:

• frame (dict) – dictionary with [‘a’,’b’,’c’] keys for each chip containing HDULists with flux, error, and mask
• rows (list) – list of rows to look for lines in
• waves (list) – list of wavelength arrays to be used to get initial pixel guess for input lines
• lines – table with desired lines, must have at least CHIPNUM and WAVE tags
• out= (str) – optional name of output ASCII file for lines (default=None)

Returns :

structure with identified lines, with tags chip, row, wave, peak, pixrel, dpixel, frameid

apogee.apred.wave.func_multi_poly(x, *pars, **kwargs)

Convert pixel to wavelength using wavecal parameters

w = poly(x + offset(group,chip)) pars = [npoly coefficients, ngroup*3 chip offsets]

Parameters:x – [3,npts] array of (pixel,chip,group)

Returns :

wave (float) : wavelength array for input pixel(s), parameters

apogee.apred.wave.gauss(x, a, x0, sig)

Evaluate Gaussian function

apogee.apred.wave.gaussbin(x, a, x0, sig)

Evaluate integrated Gaussian function

apogee.apred.wave.getgroup(groups)

Given input list of group ids that may not be consecutive, return consecutive list

apogee.apred.wave.getskywave(frame, waveid, group=-1, vers='test', telescope='apo25m', plugmap=None)

Given input frame and waveid/group for frame taken without dither move to waveid, return skyline wavelengths

apogee.apred.wave.myerf(t)

Evaluate function that integrates Gaussian from -inf to t

apogee.apred.wave.peakfit(spec, pix0, estsig=5, sigma=None, mask=None, plot=False, func=<function gaussbin>)

Return integrated-Gaussian centers near input pixel center

Parameters:

• spec (float) – data spectrum arraty
• pix0 (float) – initial pixel guess
• estsig (float) – initial guess for window width=5*estsig (default=5)
• sigma (float) – uncertainty array (default=None)
• mask (float) – mask array (default=None), NOT CURRENTLY IMPLEMENT
• plot (bool) – plot spectrum and fit in current plot window (default=False)
• func (function) – user-supplied function to use to fit (default=gaussbin)

apogee.apred.wave.pix2wave(pix, wave0)

convert pixel(s) to wavelength(s) Args :

pix : pixel(s) to get wavelength at wave0 : array with wavelength as a function of pixel number

Returns :

wavelength(s) in AA

apogee.apred.wave.plot_apWave(nums, apred='current', inst='apogee-n', out=None, hard=False)

Diagnostic plots of wavecal

apogee.apred.wave.refine(oldpars, npoly=4)

Refine wavelength solution by averaging over groups,and smoothing over rows

apogee.apred.wave.save_apWave(pars, out=None, group=0, rows=array([ 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108, 109, 110, 111, 112, 113, 114, 115, 116, 117, 118, 119, 120, 121, 122, 123, 124, 125, 126, 127, 128, 129, 130, 131, 132, 133, 134, 135, 136, 137, 138, 139, 140, 141, 142, 143, 144, 145, 146, 147, 148, 149, 150, 151, 152, 153, 154, 155, 156, 157, 158, 159, 160, 161, 162, 163, 164, 165, 166, 167, 168, 169, 170, 171, 172, 173, 174, 175, 176, 177, 178, 179, 180, 181, 182, 183, 184, 185, 186, 187, 188, 189, 190, 191, 192, 193, 194, 195, 196, 197, 198, 199, 200, 201, 202, 203, 204, 205, 206, 207, 208, 209, 210, 211, 212, 213, 214, 215, 216, 217, 218, 219, 220, 221, 222, 223, 224, 225, 226, 227, 228, 229, 230, 231, 232, 233, 234, 235, 236, 237, 238, 239, 240, 241, 242, 243, 244, 245, 246, 247, 248, 249, 250, 251, 252, 253, 254, 255, 256, 257, 258, 259, 260, 261, 262, 263, 264, 265, 266, 267, 268, 269, 270, 271, 272, 273, 274, 275, 276, 277, 278, 279, 280, 281, 282, 283, 284, 285, 286, 287, 288, 289, 290, 291, 292, 293, 294, 295, 296, 297, 298, 299]), npoly=4, frames=[], rms=None, sig=None, allpars=None)

Write the apWave files in standard format given the wavecal parameters

apogee.apred.wave.scalarDecorator(func)

Decorator to return scalar outputs for wave2pix and pix2wave

apogee.apred.wave.skycal(planfile, out=None, inst=None, waveid=None, group=-1, skyfile='airglow', vers=None, nosky=False)

Determine positions of skylines for all frames in input planfile

apogee.apred.wave.skywaves()

get skyline wavelengths from some particular frames, 16390029 and 22430033

apogee.apred.wave.test()

test routine for peakfity

apogee.apred.wave.wave2pix(wave, wave0)

convert wavelength to pixel given wavelength array Args :

wave(s) : wavelength(s) (AA) to get pixel of wave0 : array with wavelength as a function of pixel number

Returns :

pixel(s) in the chip

apogee.apred.wave.wavecal(nums=[2420038], name=None, vers='current', inst='apogee-n', rows=[150], npoly=4, reject=3, pixellim=[3, 2043], plot=False, hard=True, verbose=False, clobber=False, init=False, nofit=False, test=False)

APOGEE wavelength calibration

Solves for wavelength calibration given input frameid(s) allowing for a single polynomial wavelength solution with offsets for each chip and each group of wavecals

Parameters:

• nums (list of ints) – list of input frame ids
• name (int) – name for output ID (default = first frame from list of nums)
• vers (str) – apred version (defaut=’current’)
• inst (str) – instrument (default=’apogee-n’)
• rows (list of ints) – list of rows to get solutions at
• npoly (int) – order of polynomial fit (default=4)
• reject (float) – factor for line rejection
• verbose (bool) –
• verbose – plot fits (default=False)
• clobber (bool) – forces remeasuring lines
• init (bool) – if true, use simple quadratic estimate for first guess (default=False)
• nofit (bool) – only find lines, skip fit (default=False)
• test (bool) – if True, use polynomial from first set, only let centers shift for all other sets