usage of apedx model

HI:

I'm trying to use apedx to model an spectral line. I'm following the example on LSF_script.sl from David Huenemoerder:


%%%%%%%%%%%% Fit with APED plasma %%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%-
%  These models take longer to evaluate, so depending on the speed of
%  the computer, it may be a few seconds per fit, or much longer.
%
% define the model functions.
%
% aped_fit_models.sl is from dph and defines several plasma functions;
% the one used here is a single temperature plasma, and has parameters
% for temperature, radial velocity, turbulent broadening, and
% normalization.   
%
% Available from:
%
% http://space.mit.edu/CXC/ISIS/contrib/Aped_fit/aped_fit_models.sl
% (and in the isis contrib directory, with a bug)

()=evalfile("../../../isis_fun/aped_fit_models.sl");   % the aped model definitions

% define a single-temperature plasma model
%  (list_par will show the params)

fit_fun("apedx(1)");    % (renamed apedx because of
                        % collision w/ aped variable for database.)

%%%
% the following are the result of an ad hoc fit to Ne X 12A:
%
set_par("apedx(1).norm", 0.048, 0, 0.0001, 0.9);
set_par("apedx(1).temp", 10.^6.8, 1, 1.e6, 3.e8);   % a good enough temperature

% The lines do seem to need more-than-instrumental broadening.  Is
% thermal broadening  significant for Ne X at 12A in the HEG?

use_thermal_profile;                          % turn on broadening, 
set_par("apedx(1).vturb", 0, 1, 0.0, 200);    %  but freeze vturb to 0.


% Since there seems to be a problem with the wavelength scale, we will 
% work around it by letting the radial velocity be free.  (This is OK
% for single orders, but not for joint fits w/ different offsets.)

set_par("apedx(1).vrad", 100, 0, 0, 200);   % vrad

% ignore all but HEG -1 for this test:
%
xlo=18.7;   xhi=19.7;
xnotice(3, xlo, xhi);
ignore([1,2,4,5]);

yrange(0,450);

() = renorm_counts;   % adjust the normalization to the data.

list_data;   % inspect lists
list_par;

xrange(xlo-.1, xhi+.1);         % select Ne X + Fe  region
label("Wavelength [\\A]", "Flux [counts/bin]", "ObsID 1103/Capella/HEG -1");
plot_data_counts(3,red);
oplot_model_counts(3,green);

rplot_counts(3);             % residual plot

() = fit_counts;

rplot_counts(3);

% the profile isn't really as bad as it looks. (but the model isn't
%   perfect; residuals aren't quite random);
list_par;


% compute contour map for data fit
%  dcv = conf_grid("apedx(1).vrad", 0.0, 150.0, 20);
%  dcf = conf_grid("apedx(1).norm", 0.012, 0.023, 20);
%  dfvmap = conf_map_counts(dcf, dcv);


But I'm getting only continuum fitting and no lines.

If I ask isis for the value of Aped_Filter I got:

isis> print(Aped_Filter);
{type=NULL, lines=NULL}


which I guess that are the values in model_spectrum to have spectral lines (right?), if not, how do I change it?


Thanks
Juan
########################################################
GERARDO JUAN MANUEL LUNA
Center for Astrophysics/Smithsonian Astrophysical Observatory
60 Garden st. MS 15
Cambridge-MA-USA
02138
########################################################



      ____________________________________________________________________________________
¡Buscá desde tu celular!

Yahoo! oneSEARCH ahora está en Claro

http://ar.mobile.yahoo.com/onesearch

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