Table model for spectropolarimetry example

Writing a table for a spectropolarimetric table requires providing three spectra for each point in the table. The initial set-up of the table descriptors now becomes:

test.setModelName("SPTest")
test.setModelUnits(" ")
test.setisRedshift(True)
test.setisAdditive(True)
test.setisError(False)
filters = ["Stokes:0","Stokes:1","Stokes:2"]
test.setFiltExps(filters)

and the loops setting the example spectra:

for i1 in range(11):
  for i2 in range(5):
    flux = empty((99))
    for k in range(3):
      for j in range(99): flux[j] = 0.2*i1+10*(4.6+0.2*i2)+j*0.1+k/3.0
      testspec = tableSpectrum()
      testspec.setParameterValues(np.array([0.2*i1,4.6+0.2*i2]))
      testspec.setFlux(flux)
      test.pushSpectrum(testspec)




HEASARC Home | Observatories | Archive | Calibration | Software | Tools | Students/Teachers/Public

Last modified: Wednesday, 28-Aug-2024 16:49:26 EDT

HEASARC Staff Scientist Position - Applications are now being accepted for a Staff Scientist with significant experience and interest in the technical aspects of astrophysics research, to work in the High Energy Astrophysics Science Archive Research Center (HEASARC) at NASA Goddard Space Flight Center (GSFC) in Greenbelt, MD. Refer to the AAS Job register for full details.