Re: plotting model over data gaps

Rhaana,

  It is possible to evaluate the model spectrum on a user grid (independent
of arfs or rfms) and then plot it as you wish, see
isis> help eval_fun
  Because the model spectrum is in units of photons/s/cm^2/bin
it is easiest to overplot the model on data that have been
flux corrected.  Assuming you have data and model in index 1
and are using plot_unit("keV") here's an example:

Flux correct your data and plot it:

isis> flux_corr(1);
isis> plot_bin_density;
isis> plot_data_flux(1);  % Photons/s/cm^2/keV

Evaluate the model and plot it.  The grid is in Ang's
and since I'm assuming you are plotting in keV I have made liberal
use of the constant Const_hc :

isis> (lo, hi) = linear_grid (1,30,2000);
isis> y = eval_fun (lo, hi);
isis> oplot(Const_hc/((lo+hi)/2.0), y/(Const_hc/lo - Const_hc/hi), 3);

re-plot the data on top for clarity

isis> oplot_data_flux(1,1);

This gives:




On Mar 15, 2011, at 11:45 AM, Rhaana Starling wrote:

> Dear ISIS experts,
> 
> I would like to be able to overplot my model fit on the data and interpolate/extrapolate it during data gaps. I'm just using the standard plot routines from .isisrc_plots you provide on the website and simple models such as smc*phabs*bknpl (from xspec, and smc is one I wrote into s-lang myself).
> Is it possible to do this, and does the routine/command already exist?
> 
> Also, can the model be overplotted as a smooth continuous line, ie using the values of the centres of bins only and doing a kind of dot-to-dot, rather than the default whole-bin span which in my case looks like steps?
> 
> Thanks,
> Rhaana
> ----
> 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.hidden> with the first line of the message as:
> unsubscribe

    ----------------------------------------------------------
    Daniel Dewey                           dd_at_email.domain.hidden    MIT Kavli Institute        http://space.mit.edu/home/dd
    ----------------------------------------------------------


----
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