Re: output plots

From: David P. Huenemoerder <dph>
Date: Mon, 4 Nov 2002 17:13:06 -0500
Hi Saku - 

    Saku> how do I plot data+model and data/model ratio?  the command
    Saku> rplot_counts (1); gives data+model and delta chi?  Also can
    Saku> I make the xaxis be wavelength instead of energy?

 plot_unit("Angstrom");

should do it. For grating spectra, wavelength is the default, I never
have to set this, but it could come from the bin_lo, bin_hi columns of
the data).


rplot_counts();  % yes - data, model , and delta chi in lower plot,
                   but optionally other quantities:

isis> help rplot_counts;
rplot_counts

 SYNOPSIS
    Plot counts data and model with residuals

 USAGE
    rplot_counts (hist_index)

 DESCRIPTION
    This function generates a two-paned plot, the upper pane
    showing the counts data with model overlaid and the lower pane
    showing the residuals.  The type of residuals plotted is
    determined by the value of the intrinsic variable
    Isis_Residual_Plot_Type; supported values are:

    ___Value____Definition_________
       DELCHI   (data-model)/sigma
       DIFF     (data-model)
       RATIO    (data/model)

    The default is Isis_Residual_Plot_Type=DELCHI.

 SEE ALSO
    rplot_flux, plot_data_counts, title



So you do this:

isis> Isis_Residual_Plot_Type=RATIO;
isis> rplot_counts(1);


Or you could define a function if you flip back and forth a lot:

define rpr(n)
{
  variable save_type = Isis_Residual_Plot_Type;
  Isis_Residual_Plot_Type=RATIO;

  rplot_counts(n);

  Isis_Residual_Plot_Type = save_type;

}
  

isis> rpr(1);





(FYI on the HEAD network, the latest isis is usually in
  /soft/isis/isis/bin/isis.)



-- Dave

David Huenemoerder (617-253-4283; fax: 253-8084)
MIT Center for Space Research/Chandra X-ray Center
70 Vassar St., NE80-6023,
Cambridge, MA  02139
[Admin. Asst.: Deborah Gage 617-253-0228, dgage_at_email.domain.hidden----
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
Received on Mon Nov 04 2002 - 17:13:11 EST

This archive was generated by hypermail 2.2.0 : Thu Mar 15 2007 - 08:45:50 EDT