<!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>Appendix: Commands used to generate plots — 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="prev" title="Discovery case" href="main.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="main.html" title="Discovery case"
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>
<h4>Previous topic</h4>
<p class="topless"><a href="main.html"
title="previous chapter">Discovery case</a></p>
<h3>This Page</h3>
<ul class="this-page-menu">
<li><a href="_sources/commands.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="appendix-commands-used-to-generate-plots">
<h1>Appendix: Commands used to generate plots<a class="headerlink" href="#appendix-commands-used-to-generate-plots" title="Permalink to this headline">ΒΆ</a></h1>
<p>The version of pysynphot used for this investigation was:</p>
<div class="highlight-python"><pre>URL: https://svn.stsci.edu/svn/ssb/astrolib/tags/pysynphot_v0.61_ETC18.0
Repository Root: https://svn.stsci.edu/svn/ssb/astrolib
Revision: 1143
Last Changed Author: laidler
Last Changed Rev: 1080
Last Changed Date: 2009-08-14 14:57:05 -0400 (Fri, 14 Aug 2009)</pre>
</div>
<p>#Examine spectrum as presented in CDBS:</p>
<div class="highlight-python"><div class="highlight"><pre><span class="n">sp</span><span class="o">=</span><span class="n">S</span><span class="o">.</span><span class="n">FileSpectrum</span><span class="p">(</span><span class="s">'/grp/hst/cdbs/calspec/lds749b_mod_001.fits'</span><span class="p">)</span>
<span class="n">clf</span><span class="p">()</span>
<span class="n">plot</span><span class="p">(</span><span class="n">sp</span><span class="o">.</span><span class="n">wave</span><span class="p">,</span><span class="n">sp</span><span class="o">.</span><span class="n">flux</span><span class="p">)</span>
<span class="n">xlabel</span><span class="p">(</span><span class="n">sp</span><span class="o">.</span><span class="n">waveunits</span><span class="p">)</span>
<span class="n">ylabel</span><span class="p">(</span><span class="n">sp</span><span class="o">.</span><span class="n">fluxunits</span><span class="p">)</span>
<span class="n">title</span><span class="p">(</span><span class="n">sp</span><span class="p">)</span>
<span class="n">xlim</span><span class="p">(</span><span class="mi">1100</span><span class="p">,</span><span class="mi">1500</span><span class="p">)</span>
<span class="n">savefig</span><span class="p">(</span><span class="s">'cdbs_flam.png'</span><span class="p">)</span>
</pre></div>
</div>
<p>#Examine bin widths of spectrum:</p>
<div class="highlight-python"><div class="highlight"><pre><span class="n">dw</span><span class="o">=</span><span class="n">sp</span><span class="o">.</span><span class="n">wave</span><span class="p">[</span><span class="mi">1</span><span class="p">:]</span><span class="o">-</span><span class="n">sp</span><span class="o">.</span><span class="n">wave</span><span class="p">[:</span><span class="o">-</span><span class="mi">1</span><span class="p">]</span>
<span class="n">plot</span><span class="p">(</span><span class="n">sp</span><span class="o">.</span><span class="n">wave</span><span class="p">[</span><span class="mi">1</span><span class="p">:],</span><span class="n">dw</span><span class="p">,</span><span class="s">'.'</span><span class="p">)</span>
<span class="n">xlabel</span><span class="p">(</span><span class="s">'angstrom'</span><span class="p">)</span>
<span class="n">ylabel</span><span class="p">(</span><span class="s">'binwidth'</span><span class="p">)</span>
<span class="n">title</span><span class="p">(</span><span class="n">sp</span><span class="p">)</span>
<span class="n">savefig</span><span class="p">(</span><span class="s">'binwidth.png'</span><span class="p">)</span>
<span class="n">xlim</span><span class="p">(</span><span class="mi">1100</span><span class="p">,</span><span class="mi">1500</span><span class="p">)</span>
<span class="n">ylim</span><span class="p">(</span><span class="mi">0</span><span class="p">,</span><span class="mf">0.6</span><span class="p">)</span>
<span class="n">savefig</span><span class="p">(</span><span class="s">'binwidth_zoom.png'</span><span class="p">)</span>
</pre></div>
</div>
<p>#Previously obtain and save SYNPHOT spectrum using calcspec task:</p>
<div class="highlight-python"><div class="highlight"><pre><span class="c">#--> lpar calcspec</span>
<span class="c"># spectrum = "crcalspec$lds749b_mod_001.fits" Spectrum to calculate</span>
<span class="c"># output = "lds749b.fits" Output table name</span>
<span class="c"># (form = "counts") Desired form of output spectrum</span>
</pre></div>
</div>
<p>#Read in this spectrum to examine its wavelength array:</p>
<div class="highlight-python"><div class="highlight"><pre><span class="n">synspec</span><span class="o">=</span><span class="n">S</span><span class="o">.</span><span class="n">FileSpectrum</span><span class="p">(</span><span class="s">'lds749b.fits'</span><span class="p">)</span>
<span class="n">syndw</span><span class="o">=</span><span class="n">synspec</span><span class="o">.</span><span class="n">wave</span><span class="p">[</span><span class="mi">1</span><span class="p">:]</span><span class="o">-</span><span class="n">synspec</span><span class="o">.</span><span class="n">wave</span><span class="p">[:</span><span class="o">-</span><span class="mi">1</span><span class="p">]</span>
<span class="n">plot</span><span class="p">(</span><span class="n">synspec</span><span class="o">.</span><span class="n">wave</span><span class="p">[</span><span class="mi">1</span><span class="p">:],</span><span class="n">syndw</span><span class="p">,</span><span class="s">'.'</span><span class="p">)</span>
<span class="n">xlabel</span><span class="p">(</span><span class="s">'angstrom'</span><span class="p">)</span>
<span class="n">ylabel</span><span class="p">(</span><span class="s">'binwidth'</span><span class="p">)</span>
<span class="n">title</span><span class="p">(</span><span class="s">'crcalspec$lds749b_mod_001.fits rendered in SYNPHOT'</span><span class="p">)</span>
<span class="n">savefig</span><span class="p">(</span><span class="s">'syn_binwidth.png'</span><span class="p">)</span>
<span class="n">xlim</span><span class="p">(</span><span class="mi">1100</span><span class="p">,</span><span class="mi">1500</span><span class="p">)</span>
<span class="n">ylim</span><span class="p">(</span><span class="mi">0</span><span class="p">,</span><span class="o">.</span><span class="mi">6</span><span class="p">)</span>
<span class="n">savefig</span><span class="p">(</span><span class="s">'syn_binwidth_zoom.png'</span><span class="p">)</span>
</pre></div>
</div>
<p>#Examine the spectrum in pysynphot on the SYNPHOT waveset:</p>
<div class="highlight-python"><div class="highlight"><pre><span class="n">sp</span><span class="o">.</span><span class="n">convert</span><span class="p">(</span><span class="s">'counts'</span><span class="p">)</span>
<span class="n">plot</span><span class="p">(</span><span class="n">synspec</span><span class="o">.</span><span class="n">wave</span><span class="p">,</span><span class="n">sp</span><span class="o">.</span><span class="n">sample</span><span class="p">(</span><span class="n">synspec</span><span class="o">.</span><span class="n">wave</span><span class="p">),</span><span class="s">'.-'</span><span class="p">)</span>
<span class="n">xlim</span><span class="p">(</span><span class="mi">1100</span><span class="p">,</span><span class="mi">1500</span><span class="p">)</span>
<span class="n">ylim</span><span class="p">(</span><span class="mi">0</span><span class="p">,</span><span class="mi">30</span><span class="p">)</span>
<span class="n">xlabel</span><span class="p">(</span><span class="n">synspec</span><span class="o">.</span><span class="n">waveunits</span><span class="p">)</span>
<span class="n">ylabel</span><span class="p">(</span><span class="n">sp</span><span class="o">.</span><span class="n">fluxunits</span><span class="p">)</span>
<span class="n">title</span><span class="p">(</span><span class="s">"</span><span class="si">%s</span><span class="s"> sampled on SYNPHOT wavelength array"</span><span class="o">%</span><span class="n">os</span><span class="o">.</span><span class="n">path</span><span class="o">.</span><span class="n">basename</span><span class="p">(</span><span class="nb">str</span><span class="p">(</span><span class="n">sp</span><span class="p">)))</span>
<span class="n">savefig</span><span class="p">(</span><span class="s">'syn_sampled_counts.png'</span><span class="p">)</span>
</pre></div>
</div>
</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="main.html" title="Discovery case"
>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>
