Name = dataheader.get('OBJECT')
headerKWs['OBJECT'] = Name
headerKWs["DATE-OBS"] = dataheader.get('DATE-OBS')
headerKWs["INSTRUMENT"] = dataheader.get('INSTRUME')
headerKWs["OBSERVER"] = dataheader.get('OBSERVER')
headerKWs["EXPTIME"] = dataheader.get('EXPTIME')
Score = Moog960.Score()
Melody = Moog960.ObservedMelody(Score=Score, name=Name)
for wl, I, sn, v in zip(wlsol, data, snr, var):
header = pyfits.Header()
header.set('NAME', Name)
header.set("WLSTART", numpy.min(wl))
header.set("WLSTOP", numpy.max(wl))
Melody.addPhrase(Moog960.ObservedPhrase(observedData=[SpectralTools.Spectrum(wl=wl,
I=I, dI=I/sn, header=header, spectrum_type="OBSERVED")], Melody=Melody))
ID = Melody.phrases[-1].ID
parameters = {}
parameters["PHRASE"] = ID
parameters["MELODY"] = Melody.ID
parameters["SCORE"] = Melody.Score.ID
parameters["WLSTART"] = header.get('WLSTART')
parameters["WLSTOP"] = header.get('WLSTOP')
parameters["SELECTED"] = False
if not(Melody.ID in Melody.Score.observed_labels.keys()):
Melody.Score.observed_labels[Melody.ID] = {}
Melody.Score.observed_labels[Melody.ID][ID] = []
Melody.Score.observed_labels[Melody.ID][ID].append(Moog960.Label(parameters,
reference=Melody.phrases[-1].observedData[-1], Phrase=Melody.phrases[-1],
headerKWs["OBSERVER"] = dataheader.get('OBSERVER')
headerKWs["EXPTIME"] = dataheader.get('EXPTIME')
Score = Moog960.Score()
Melody = Moog960.ObservedMelody(Score=Score, name=Name)
for wl, I, sn, v in zip(wlsol, data, snr, var):
header = pyfits.Header()
header.set('NAME', Name)
header.set("WLSTART", numpy.min(wl))
header.set("WLSTOP", numpy.max(wl))
Melody.addPhrase(
Moog960.ObservedPhrase(observedData=[
SpectralTools.Spectrum(wl=wl,
I=I,
dI=I / sn,
header=header,
spectrum_type="OBSERVED")
],
Melody=Melody))
ID = Melody.phrases[-1].ID
parameters = {}
parameters["PHRASE"] = ID
parameters["MELODY"] = Melody.ID
parameters["SCORE"] = Melody.Score.ID
parameters["WLSTART"] = header.get('WLSTART')
parameters["WLSTOP"] = header.get('WLSTOP')
parameters["SELECTED"] = False
if not (Melody.ID in Melody.Score.observed_labels.keys()):
Melody.Score.observed_labels[Melody.ID] = {}
