Re: simultaneous fitting

From: Herbert Pablo <hpablo_at_email.domain.hidden>
Date: Mon, 02 Feb 2009 17:49:11 -0500
Hey,

Thanks for the reply.

First to make things a little more clear. All the combining is done in one 
script and then a fitting script is run. To that end, this is what I'm curious 
about. From what I understand, once I combine two data sets into one, there is 
only 1 piece of information now and that is what I want to fit. If that is the 
case why would I have to specify wavelength ranges for both sets of data if 
they are already combined?


On a secondary note, I think there is an issue with my script to combine 
datasets. If I run the same script using counts everything is fine and I get 
what I would expect. However, if i run the script with flux, something strange 
is happening. If I plot the plus and minus orders before I combine them, then 
they look fine. But, if I plot the plus and minus orders after I combine, one 
of them is fine and the other is empty (horizontal line at 0). I'm sure this 
is affecting my fit, but honestly I don't know what the problem is. For 
brevity, I'll just write the combining part of the code below:

variable h1_id;
variable heg1,weights;

weights=[0.5,0.5];
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


match_dataset_grids (1,2);
h1_id = combine_datasets (1,2,weights);

heg1=get_combined (h1_id, &get_data_flux);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

xrange(1.0,10.0);
variable yname = latex2pg("Flux [Photons/s/cm^{2}/bin]");
open_plot,3;
hplot(heg1);

If you see any obvious errors here let me know.

Thanks,
Bert


----
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 Mon Feb 02 2009 - 17:49:17 EST

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