Re: using rplot_counts with get_combined

On Apr 30, 2007, at 4:39 PM, David P. Huenemoerder wrote:

> Mike Nowak has written a suite of plotting functions which handles the
> combined indices nicely and you can find them here (his .isisrc_plots
> file):
>
>    http://space.mit.edu.ezproxyberklee.flo.org/home/mnowak/isis_vs_xspec/download.html

>     juan> -Also, is the residual from rplot_counts defined
>     juan> equally to XSPEC delchi or chi ?

Thanks for the shout out, Dave.

Basically, there I am doing the combined data a little bit "by hand", 
rather than relying on ISIS intrinsics.  That's partly so I can do the 
"combined" plotting whether the data are combined or not, and also so I 
can separately combine the data, or just their residuals.

Less well-documented on that page:

plot_bsub_all([1,-2],popt);
plot_bsub_all([-1,-2],popt);
plot_bsub_all([-1,2],popt);

*or*

plot_bsub_all([-1],popt);

all produce the same combined data plot, if 1 & 2 share a common grid 
(and noticed range), and the latter if 1 refers to combined group #1 
(which in Dave's example, g==1, and represents data sets 1 & 2; i.e., 
plot_bsub_all([-g],popt), if there is only one index, gives you the 
combined group, plot_bsub_all([#,#,#-#,...],popt) combines all the # if 
there are more than one, and *any* are negative.).

I have put toggles in so that you can get equivalents of either delchi 
or chi in XSPEC terms.  Furthermore, the code is semi-intelligent about 
figuring out whether that chi square is referred to model counts, or 
the Gehrels definition of chi^2, etc.

One bit I haven't I figured out and corrected yet:  the code can give 
funny answers when plotting flux-corrected data *and* mixing/matching 
heg/meg data to which you've applied match_dataset_grid.  (Non-flux 
corrected data, or flux corrected heg viewed by itself, or 
flux-corrected meg viewed by itself, are all fine.)

-Mike

----
You received this message because you are
subscribed to the isis-users list.
To unsubscribe, send a message to
isis-users-request_at_email.domain.hiddenwith the first line of the message as:
unsubscribe