Comments on: Everybody needs crampons http://hea-www.harvard.edu/AstroStat/slog/2010/sherpa-cheat-sheet/ Weaving together Astronomy+Statistics+Computer Science+Engineering+Intrumentation, far beyond the growing borders Fri, 01 Jun 2012 18:47:52 +0000 hourly 1 http://wordpress.org/?v=3.4 By: vlk http://hea-www.harvard.edu/AstroStat/slog/2010/sherpa-cheat-sheet/comment-page-1/#comment-970 vlk Sun, 13 Jun 2010 03:11:44 +0000 http://hea-www.harvard.edu/AstroStat/slog/2007/sherpa-cheat-sheet/#comment-970 Thanks for the link to the Sherpa scripts page, Doug. The best way to become familiar with a new language is to see lots of simple examples. Also thanks for the other comments. I will probably include those comments by moving this page off to gWave shortly -- it is much easier to update and make changes there than in the blog format. Thanks for the link to the Sherpa scripts page, Doug. The best way to become familiar with a new language is to see lots of simple examples.

Also thanks for the other comments. I will probably include those comments by moving this page off to gWave shortly — it is much easier to update and make changes there than in the blog format.

]]>
By: Doug Burke http://hea-www.harvard.edu/AstroStat/slog/2010/sherpa-cheat-sheet/comment-page-1/#comment-969 Doug Burke Sat, 12 Jun 2010 22:16:48 +0000 http://hea-www.harvard.edu/AstroStat/slog/2007/sherpa-cheat-sheet/#comment-969 Vinay, Have you looked at the <a href="http://cxc.harvard.edu/sherpa/scripts/index.py.html" rel="nofollow">Sherpa scripts</a> page. For instance, some of the plot options are discussed in the <a href="http://cxc.harvard.edu/sherpa/threads/setplot_manual/index.py.html" rel="nofollow">plotting thread</a>, although I couldn't see the error bars explicitly addressed. Try <code> p = get_data_plot_prefs() print(p) p["yerrorbars"] = False </code> To change the units of the y error bars, use the set_analysis() routine, for instance <code> set_analysis("wave", type="counts") </code> Also, for error estimation, we suggest you use <code> conf() </code> rather than projection nowadays, and note that the error analysis will be done in parallel (so expect all your cores to get used up). And finally, <code> ,ahelp set_analysis </code> can also be used to avoid quoting the first argument by using a comma before the command. This is an ipython-ism, and can be used for other commands, e.g. <code> ,load_pha foo.pha </code> Doug Vinay,

Have you looked at the Sherpa scripts page. For instance, some of the plot options are discussed in the plotting thread, although I couldn’t see the error bars explicitly addressed. Try


p = get_data_plot_prefs()
print(p)
p["yerrorbars"] = False

To change the units of the y error bars, use the set_analysis() routine, for instance


set_analysis("wave", type="counts")

Also, for error estimation, we suggest you use


conf()

rather than projection nowadays, and note that the error analysis will be done in parallel (so expect all your cores to get used up).

And finally,


,ahelp set_analysis

can also be used to avoid quoting the first argument by using a comma before the command. This is an ipython-ism, and can be used for other commands, e.g.


,load_pha foo.pha

Doug

]]>