<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd"> <html xmlns="http://www.w3.org/1999/xhtml"> <head> <meta http-equiv="Content-Type" content="text/html; charset=utf-8" /> <title>Discovery case — TSR 2009-02 2009-02 documentation</title> <link rel="stylesheet" href="_static/stsci_sphinx.css" type="text/css" /> <link rel="stylesheet" href="_static/pygments.css" type="text/css" /> <script type="text/javascript"> var DOCUMENTATION_OPTIONS = { URL_ROOT: '', VERSION: '2009-02', COLLAPSE_INDEX: false, FILE_SUFFIX: '.html', HAS_SOURCE: true }; </script> <script type="text/javascript" src="_static/jquery.js"></script> <script type="text/javascript" src="_static/underscore.js"></script> <script type="text/javascript" src="_static/doctools.js"></script> <link rel="top" title="TSR 2009-02 2009-02 documentation" href="index.html" /> <link rel="next" title="Appendix: Commands used to generate plots" href="commands.html" /> <link rel="prev" title="TSR 2009-02: Choppy source plot in ETC 18.0" href="index.html" /> </head> <body> <div class="related"> <h3>Navigation</h3> <ul> <li class="right" style="margin-right: 10px"> <a href="genindex.html" title="General Index" accesskey="I">index</a></li> <li class="right" > <a href="commands.html" title="Appendix: Commands used to generate plots" accesskey="N">next</a> |</li> <li class="right" > <a href="index.html" title="TSR 2009-02: Choppy source plot in ETC 18.0" accesskey="P">previous</a> |</li> <li><a href="index.html">TSR 2009-02 2009-02 documentation</a> »</li> </ul> </div> <div class="sphinxsidebar"> <div class="sphinxsidebarwrapper"> <p class="logo"><a href="index.html"> <img class="logo" src="_static/stsci_logo.png" alt="Logo"/> </a></p> <h3><a href="index.html">Table Of Contents</a></h3> <ul> <li><a class="reference internal" href="#">Discovery case</a></li> <li><a class="reference internal" href="#investigation">Investigation</a></li> <li><a class="reference internal" href="#conclusions">Conclusions</a></li> </ul> <h4>Previous topic</h4> <p class="topless"><a href="index.html" title="previous chapter">TSR 2009-02: Choppy source plot in ETC 18.0</a></p> <h4>Next topic</h4> <p class="topless"><a href="commands.html" title="next chapter">Appendix: Commands used to generate plots</a></p> <h3>This Page</h3> <ul class="this-page-menu"> <li><a href="_sources/main.txt" rel="nofollow">Show Source</a></li> </ul> <div id="searchbox" style="display: none"> <h3>Quick search</h3> <form class="search" action="search.html" method="get"> <input type="text" name="q" /> <input type="submit" value="Go" /> <input type="hidden" name="check_keywords" value="yes" /> <input type="hidden" name="area" value="default" /> </form> <p class="searchtip" style="font-size: 90%"> Enter search terms or a module, class or function name. </p> </div> <script type="text/javascript">$('#searchbox').show(0);</script> </div> </div> <div class="document"> <div class="documentwrapper"> <div class="bodywrapper"> <div class="body"> <div class="section" id="discovery-case"> <h1>Discovery case<a class="headerlink" href="#discovery-case" title="Permalink to this headline">¶</a></h1> <p>This issue was reported by Cristina Oliveira during user testing of ETC 18.0 (etc.stsci.edu) on 26 October 2009, for exposure ID COS.A192124 which specifies LDS749B (crcalspec$lds749b_mod_001.fits) as the target spectrum.</p> <div class="figure"> <img alt="Source spectrum plot, in counts, as rendered in ETC 18.0" src="_images/src_A20631.png" /> <p class="caption">Spectrum, in counts, as rendered in ETC 18.0</p> </div> <p>She observed that the file in CDBS looked much smoother.</p> <div class="figure"> <img alt="Corresponding spectrum, in flam, when plotted manually" src="_images/cdbs_flam.png" /> <p class="caption">Corresponding spectrum, in flam, when plotted manually</p> </div> </div> <div class="section" id="investigation"> <h1>Investigation<a class="headerlink" href="#investigation" title="Permalink to this headline">¶</a></h1> <p>Using pysynphot v0.61, the observed ETC behavior was duplicated. Further examination of the input spectrum revealed that its wavelength table is far from evenly spaced. Its bin widths range from 0.00098A to 50A. In the specific domain of interest (1100A to 1500A), the bin size ranges from 0.00195 to 0.477.</p> <div class="figure"> <img alt="_images/binwidth.png" src="_images/binwidth.png" /> <p class="caption">Scatter plot showing bin width as a function of wavelength</p> </div> <div class="figure"> <img alt="_images/binwidth_zoom.png" src="_images/binwidth_zoom.png" /> <p class="caption">Scatter plot zoomed to x(1100:1500) and y(0:0.6)</p> </div> <p>Unlike most of the other units supported by (py)synphot, “counts” is not a flux density. When a spectrum is plotted in counts, the y value per bin is obtained by integrating the flux density over the bin width, then converting that integrated flux into counts. Therefore, the appearance of a plotted spectrum in counts is critically dependent on the underlying wavelength array in which it is presented.</p> <p>This was not apparent in SYNPHOT plots (and, therefore, not in versions of the ETC which used synphot for plotting) because SYNPHOT uses a logarithmically spaced wavelength array when combining spectra (and SYNPHOT always combines spectra with Vega, even if they’re just being read in from a file and then plotted).</p> <div class="figure"> <img alt="_images/syn_binwidth.png" src="_images/syn_binwidth.png" /> <p class="caption">bin width spacing of SYNPHOT spectrum</p> </div> <div class="figure"> <img alt="_images/syn_binwidth_zoom.png" src="_images/syn_binwidth_zoom.png" /> <p class="caption">SYNPHOT bin width zoomed to x(1100:1500), y(0:0.6)</p> </div> <p>When pysynphot is used to sample the LDS749b spectrum, in counts, on the wavelength array used by SYNPHOT, the appearance is much more “reasonable”.</p> <div class="figure"> <img alt="_images/syn_sampled_counts.png" src="_images/syn_sampled_counts.png" /> <p class="caption">pysynphot LDS749B spectrum sampled on SYNPHOT wavelength array</p> </div> </div> <div class="section" id="conclusions"> <h1>Conclusions<a class="headerlink" href="#conclusions" title="Permalink to this headline">¶</a></h1> <p>The peculiar appearance of the plot is an artifact of the irregular wavelength array in the spectrum being plotted, and the fact that a spectrum presented in counts has been integrated over wavelength in each bin.</p> <p>Presenting a source spectrum in counts is not a useful plot, and the ETC should probably not be presenting this plot in counts.</p> <p>This was a useful case study in differences between SYNPHOT and pysynphot which will be added to the pysynphot documentation for SYNPHOT users.</p> </div> </div> </div> </div> <div class="clearer"></div> </div> <div class="related"> <h3>Navigation</h3> <ul> <li class="right" style="margin-right: 10px"> <a href="genindex.html" title="General Index" >index</a></li> <li class="right" > <a href="commands.html" title="Appendix: Commands used to generate plots" >next</a> |</li> <li class="right" > <a href="index.html" title="TSR 2009-02: Choppy source plot in ETC 18.0" >previous</a> |</li> <li><a href="index.html">TSR 2009-02 2009-02 documentation</a> »</li> </ul> </div> <div class="footer"> © Copyright 2009, Victoria G. Laidler. Created using <a href="http://sphinx.pocoo.org/">Sphinx</a> 1.1.3. </div> </body> </html>