Re: confidence limits

From: David P. Huenemoerder <dph_at_email.domain.hidden>
Date: Tue, 11 Jan 2005 12:25:19 -0500
Hi Saku - 

Just an addendum to what John said.  You might simplify a bit more.
_A() will convert the whole structure, and save you from explicitly
reversing other arrays (or worse, forgetting to).

    isis> d = get_data_counts(1);
    isis> (elo, ehi) = _A (d.bin_lo, d.bin_hi);
    isis> de_value = reverse(d.value);
    isis> writecol ("sm101dat8.ascii", elo, ehi, de_value / t / e_width);

f = _A( d ) ; 

dE = f.bin_hi - f.bin_lo ; 

writecol(  "sm101dat8.ascii", 
            f.bin_lo, 
            f.bin_hi,
            f.value/ t / dE,     % phot/cm^2/s/keV
            f.err  / t / dE ) ;  % error in same units.


-- 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.: Elaine Tirrell, 617-253-7480, egt_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 Tue Jan 11 2005 - 12:51:56 EST

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