Re: systematic errors

Hi,

>> I wasn't hard at work, but rather on vaccation.
>> (But I'll only be left like *gone* tomorrow...)
>
> Thanks! It would be really nice to have the patch
> on Mike's plotting routines, I use them a lot.

Argh! I just realized that all my efforts to provide the patch from a 
*very* bad internet connection before I had left were useless, as I was 
so stupid to somehow reply to myself only...

Cheers,

Manfred,
who'll most likely be online properly again next week


----- Forwarded Message -----
   Date: Mon, 3 Aug 2009 19:33:28 +0200
   From: Manfred Hanke <Manfred.Hanke_at_email.domain.hiddenSubject: Re: systematic errors for fancy plots
     To: Manfred Hanke <Manfred.Hanke_at_email.domain.hidden
Hi,

> I once had prepared a patch to an elder version of your fancy 
> plotting routines that incorporates ISIS' systematic errors. I can 
> try to update this and send the new version,
Huh, the scripts have changed a lot in recent times (well, during the 
last year or so...). So I wasn't able to learn too much from some quick 
diff's.
However, as far as I could figure out, the patch for the fancy plot 
functions might indeed be three lines only:

--- isis_fancy_plots.sl
+++ isis_fancy_plots.sl	2009-08-03 19:00:01 +0200
_at_email.domain.hiddenimplements("isis_fancy_plots");

% Most recent update: March 3, 2009
+% 2009-08-03 (vacation time):
+%   Manfred Hanke added syst_err_frac-support to make_data

%  Soon to undergo a major overhaul, which will break backward
%  compatability, but will make the functions somewhat simpler, and
%  make maintaining them much easier.
_at_email.domain.hidden       indx = ();

       data = get_data_counts(indx);
+      variable sys_err_frac = get_sys_err_frac(indx);
+      if(length(sys_err_frac))
+        data.err = sqrt( data.err^2 + (sys_err_frac*data.value)^2 );
+           dnote = get_data_info(indx).notice;

       if(modc > 0)


The same might be required after the
       data = get_data_counts(indx);
assignment in make_counts.


I suggest another patch to the plot color defaults, which makes 
different data sets shown in different colors, with residuals 
corresponding to the data:

--- isis_fancy_plots.sl
+++ isis_fancy_plots.sl	2009-08-03 19:00:02 +0200
_at_email.domain.hidden

style_df = 4;  % Default symbol is a circle, if current default is a dot
-mcol_df = 2;   % Default model color red
+mcol_df = 1;   % Default model color black
res_df = 0;    % Default residual is none at all!
power_df = 1;  % Default unfolded spectrum is Photons/cm^2/sec/Unit
oplt_df = 0;   % Default is to erase the old plot
_at_email.domain.hidden{
    plopt = get_plot_options;
    col_df = plopt.color;
    if(plopt.point_style==1){ style_df = plopt.point_style; }
    bkg = Integer_Type[indx_ln]; bkg[*] = bkg_df;
-   dcol = Integer_Type[indx_ln]; dcol[*] = col_df;
-   decol = Integer_Type[indx_ln]; decol[*] = col_df;
+   dcol = Integer_Type[indx_ln]; dcol[*] = [2:indx_ln+1];  % different 
colors for the different data sets
+   decol = Integer_Type[indx_ln]; decol[*] = dcol;  % same color as for data
    mcol = Integer_Type[indx_ln];    if( get_fit_fun() != NULL and 
get_fit_fun() != "bin_width(1)" )
    {
_at_email.domain.hidden    {
       mcol[*] = 0;
    }
-   rcol = Integer_Type[indx_ln]; rcol[*] = col_df;
-   recol = Integer_Type[indx_ln]; recol[*] = col_df;
+   rcol = Integer_Type[indx_ln]; rcol[*] = dcol;  % same color as for data
+   recol = Integer_Type[indx_ln]; recol[*] = dcol;  % same color as for data
    dsym = Integer_Type[indx_ln]; dsym[*] = style_df;
    rsym = Integer_Type[indx_ln]; rsym[*] = style_df;
    yrng = [NULL];



The other change I once made was to have reasonable defaults for 
x/yrng, considering all data sets to be plotted. However, what worked 
with an .isisrc_plots from 2008, it doesn't work anymore with a newer 
version of isis_fancy_plots.sl, because the interal structures seem to 
have changed. But Mike has already promised to clean them up on a 
medium time scale.

Cheers,

Manfred

P.S. As I believe that the mailing list does not forward attachments, 
I've put my changed version of isis_fancy_plots.sl on 
http://www.sternwarte.uni-erlangen.de/~hanke/ISIS/isis_fancy_plots.sl
or
http://www.sternwarte.uni-erlangen.de/~hanke/ISIS/isis_fancy_plots.html


----------------------------------------------------------------
This message was sent using IMP, the Internet Messaging Program.


================================================================================================
This message and any attachments are intended for the use of the addressee or addressees only. The
unauthorised disclosure, use, dissemination or copying (either in whole or in part) of its content
is prohibited. If you received this message in error, please delete it from your system and notify
the sender. E-mails can be altered and their integrity cannot be guaranteed. ESA shall not be liable
for any e-mail if modified.
=================================================================================================

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