Re: Rebinning to grid specified by psextract

From: Manfred Hanke <Manfred.Hanke_at_email.domain.hidden>
Date: Tue, 22 Jun 2010 21:33:07 +0200
Hi,

> When I extracted the pi spectrum using the psextract tool from CIAO,
> I specified a minimum number of counts per bin (following the
> example, 15). However, when I load the data into ISIS, it reverts
> back to the raw grid.
I have to admit that I don't know psextract, but on the ISIS side, it
might be that you have to set
    Isis_Use_PHA_Grouping = 1;
see
    help("load_data");
which contains the paragraph:
>> By default, isis ignores the GROUPING column of the input PHA file.
>> To automatically apply the grouping on input, set the intrinsic
>> variable Isis_Use_PHA_Grouping to a non-zero value.

[One may indeed wonder why it is set to zero by default... Is that some 
relic from former times?]

> How can I rebin the data so that it has at least 15 counts per bin
You can use
    rebin_data(id, 15);
for this task. If you want to do more sophisticated rebinning (i.e. to a 
certain S/N ratio, taking the background spectrum into account, too, or 
only rebin in a certain energy range), the mighty group function (that 
was only added with ISIS 1.5., iirc) can do the job, e.g.
    group(id; min_sn=4);

> and then rebin the ARF and RMF to match these grids,
I believe that using ISIS, you will not have to (worry about) 
rebin(ning) ARFs or RMFs, as this is done internally in memory.

Cheers,

Manfred
----
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 Jun 22 2010 - 15:33:45 EDT

This archive was generated by hypermail 2.3.0 : Fri May 02 2014 - 08:35:46 EDT