combined_datasets questions

Hi John - 

I'm having some confusion or problems w/ combine_datasets and 
get_combined.  I loaded the HEG and MEG -1 and +1, arfs, and rmfs, and
flux_corr'd.  I matched_dataset_grids, and combined_datasets.
But when I

  f = get_combined(1, &get_data_flux);

the result is the sum of the four fluxes, not the mean.

Also,

  a = get_combined(1, &get_arf);

returns something, but it isn't a single, nor summed arf, as far as I
can tell.  get_combined isn't documented as retrieving an arf, but
there was no error.


Am I doing something wrong?  I could use weighting, but that would
probably ruin the combined counts.

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

P.S.  If you want a good test dataset - strong continuum, use

data in ~dph/d6/Proc_data/ :

h = load_data("Vela_X1_0.25/P_Standard/pha2", [1,2,3,4]);

a = load_arf("Vela_X1_0.25/P_Standard/evt1aHEG_-1_garf.fits");
a = [a, load_arf("Vela_X1_0.25/P_Standard/evt1aHEG_1_garf.fits") ] ;
a = [a, load_arf("Vela_X1_0.25/P_Standard/evt1aMEG_-1_garf.fits") ];
a = [a, load_arf("Vela_X1_0.25/P_Standard/evt1aMEG_1_garf.fits") ];

r = load_rmf("./Ard/heg_m1_rmf.fits");
r = [r, load_rmf("./Ard/heg_p1_rmf.fits") ];
r = [r, load_rmf("./Ard/meg_m1_rmf.fits") ];
r = [r, load_rmf("./Ard/meg_p1_rmf.fits") ];


for (i=1; i<5; i++) assign_rsp(i,i,i);

match_dataset_grids([3,4,1,2]);

g = combine_datasets([1,2,3,4]);


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