Re: how to get the "areascale" value

On Fri, Mar 06, 2009 at 12:26 +0100, Michael Smith wrote:
> Hello,
> 
> I'm working with XMM RGS spectra, which have an "AREASCALE"
> column in their spectral files.
> 
> I'm trying to find the number of counts per bin.
> 
> To this end, I read in the data with "load_data()" and
> run "s = get_data_counts()" to find the counts per bin from the resulting
> s.value field.
> 
> However, s.value are floats, not integer values, hence they are not counts.
> It seems that, per bin "i", s.value[i] == counts[i] / areascale[i] ?

Yes.

> 
> If that is the case, how can I find the areascale[i] values (something
> equivalent to "get_data_backscale"), so I can then multiply s.value[i]
> to get the counts in bin i?
> 

Isis doesn't store the areascale value separately, so it will
be necessary to read the values from the pha file:

   a = fits_read_col ("pha.fits", "areascal");
   
Thanks,
-John
----
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