isis> set_fit_statistic("chisqr;sigma=data"); isis> set_fit_method("lmdif"); isis> eval_counts; Parameters[Variable] = 23[5] Data bins = 1011 Chi-square = 1072.173 Reduced chi-square = 1.065778 0 isis> fit_counts; Parameters[Variable] = 23[5] Data bins = 1011 Chi-square = 935.0749 Reduced chi-square = 0.9294979 0 isis> list_par; wabs(1) * xaped(1, He_triplets(1, He_Data)) idx param tie-to freeze value min max 1 wabs(1).nH 0 0 0.9725297 0.1 10 10^22 2 He_triplets(1).density 0 1 1e+12 1e+08 1e+16 3 xaped(1).norm1 0 0 0.001693721 1e-06 0.1 4 xaped(1).norm2 0 0 0.01348842 1e-06 0.1 5 xaped(1).temperature1 0 0 1.480109e+07 2000000 2e+07 6 xaped(1).temperature2 0 0 6.792476e+07 2e+07 2e+08 7 xaped(1).density 0 1 1 0 0 8 xaped(1).vturb 0 1 30 0 100000 9 xaped(1).redshift 0 1 0 -0.1 0.1 10 xaped(1).metal_abund 0 1 1 0 10 11 xaped(1).He 0 1 1 0 0 12 xaped(1).C 0 1 1 0 0 13 xaped(1).N 0 1 1 0 0 14 xaped(1).O 0 1 1 0 0 15 xaped(1).Ne 0 1 1 0 0 16 xaped(1).Mg 0 1 1 0 0 17 xaped(1).Al 0 1 1 0 0 18 xaped(1).Si 0 1 1 0 0 19 xaped(1).S 0 1 1 0 0 20 xaped(1).Ar 0 1 1 0 0 21 xaped(1).Ca 0 1 1 0 0 22 xaped(1).Fe 0 1 0.5 0 0 23 xaped(1).Ni 0 1 1 0 0 isis> variable lo,hi; isis> (lo,hi) = conf("xaped(1).temperature1",0,0.1); Found improved fit, stat= 934.983 for param[5] = 1.35248e+07 **** Found improved fit isis> (lo,hi) = conf("xaped(1).temperature1",0,0.1); isis> print([lo,hi]); 1.32374e+07 1.54275e+07 isis> list_data; Current Spectrum List: id instrument m prt src use/nbins A R totcts exp(ksec) target 1 HETG-ACIS -1 1 1 1011/ 1016 1 1 4.5170e+03 90.990 DOAR 21 2 HETG-ACIS 1 1 1 1011/ 1016 2 2 4.2640e+03 90.990 DOAR 21 3 HETG-ACIS -1 2 1 1011/ 1016 3 3 7.3290e+03 90.990 DOAR 21 4 HETG-ACIS 1 2 1 1011/ 1016 4 4 8.3170e+03 90.990 DOAR 21 5 -1 1 1 0/ 764 5 5 8.7810e+03 90.990 DOAR 21 6 -1 2 1 0/ 1074 6 6 1.5698e+04 90.990 DOAR 21 isis> eval_counts; Parameters[Variable] = 23[5] Data bins = 1011 Chi-square = 934.9829 Reduced chi-square = 0.9294065 isis> (lo,hi) = conf("xaped(1).temperature2",0); Found improved fit, stat= 934.342 for param[6] = 6.4309e+07 **** Found improved fit isis> (lo,hi) = conf("xaped(1).temperature2",0); isis> print([lo,hi]); 6.26475e+07 6.9524e+07 isis> list_par; wabs(1) * xaped(1, He_triplets(1, He_Data)) idx param tie-to freeze value min max 1 wabs(1).nH 0 0 0.9831352 0.1 10 10^22 2 He_triplets(1).density 0 1 1e+12 1e+08 1e+16 3 xaped(1).norm1 0 0 0.001489894 1e-06 0.1 4 xaped(1).norm2 0 0 0.01380541 1e-06 0.1 5 xaped(1).temperature1 0 0 1.38829e+07 2000000 2e+07 6 xaped(1).temperature2 0 0 6.4309e+07 2e+07 2e+08 7 xaped(1).density 0 1 1 0 0 8 xaped(1).vturb 0 1 30 0 100000 9 xaped(1).redshift 0 1 0 -0.1 0.1 10 xaped(1).metal_abund 0 1 1 0 10 11 xaped(1).He 0 1 1 0 0 12 xaped(1).C 0 1 1 0 0 13 xaped(1).N 0 1 1 0 0 14 xaped(1).O 0 1 1 0 0 15 xaped(1).Ne 0 1 1 0 0 16 xaped(1).Mg 0 1 1 0 0 17 xaped(1).Al 0 1 1 0 0 18 xaped(1).Si 0 1 1 0 0 19 xaped(1).S 0 1 1 0 0 20 xaped(1).Ar 0 1 1 0 0 21 xaped(1).Ca 0 1 1 0 0 22 xaped(1).Fe 0 1 0.5 0 0 23 xaped(1).Ni 0 1 1 0 0 isis> (lo,hi) = conf("xaped(1).norm1",0); isis> print([lo,hi]); 0.00135286 0.00172229 isis> list_par; wabs(1) * xaped(1, He_triplets(1, He_Data)) idx param tie-to freeze value min max 1 wabs(1).nH 0 0 0.9831352 0.1 10 10^22 2 He_triplets(1).density 0 1 1e+12 1e+08 1e+16 3 xaped(1).norm1 0 0 0.001489894 1e-06 0.1 4 xaped(1).norm2 0 0 0.01380541 1e-06 0.1 5 xaped(1).temperature1 0 0 1.38829e+07 2000000 2e+07 6 xaped(1).temperature2 0 0 6.4309e+07 2e+07 2e+08 7 xaped(1).density 0 1 1 0 0 8 xaped(1).vturb 0 1 30 0 100000 9 xaped(1).redshift 0 1 0 -0.1 0.1 10 xaped(1).metal_abund 0 1 1 0 10 11 xaped(1).He 0 1 1 0 0 12 xaped(1).C 0 1 1 0 0 13 xaped(1).N 0 1 1 0 0 14 xaped(1).O 0 1 1 0 0 15 xaped(1).Ne 0 1 1 0 0 16 xaped(1).Mg 0 1 1 0 0 17 xaped(1).Al 0 1 1 0 0 18 xaped(1).Si 0 1 1 0 0 19 xaped(1).S 0 1 1 0 0 20 xaped(1).Ar 0 1 1 0 0 21 xaped(1).Ca 0 1 1 0 0 22 xaped(1).Fe 0 1 0.5 0 0 23 xaped(1).Ni 0 1 1 0 0 isis> (lo,hi) = conf("xaped(1).norm2",0); **** Upper confidence limit didn't converge[4]: allow wider parameter ranges? isis> print([lo,hi]); 0.0133497 0.0139204 isis> list_par; wabs(1) * xaped(1, He_triplets(1, He_Data)) idx param tie-to freeze value min max 1 wabs(1).nH 0 0 0.9831352 0.1 10 10^22 2 He_triplets(1).density 0 1 1e+12 1e+08 1e+16 3 xaped(1).norm1 0 0 0.001489894 1e-06 0.1 4 xaped(1).norm2 0 0 0.01380541 1e-06 0.1 5 xaped(1).temperature1 0 0 1.38829e+07 2000000 2e+07 6 xaped(1).temperature2 0 0 6.4309e+07 2e+07 2e+08 7 xaped(1).density 0 1 1 0 0 8 xaped(1).vturb 0 1 30 0 100000 9 xaped(1).redshift 0 1 0 -0.1 0.1 10 xaped(1).metal_abund 0 1 1 0 10 11 xaped(1).He 0 1 1 0 0 12 xaped(1).C 0 1 1 0 0 13 xaped(1).N 0 1 1 0 0 14 xaped(1).O 0 1 1 0 0 15 xaped(1).Ne 0 1 1 0 0 16 xaped(1).Mg 0 1 1 0 0 17 xaped(1).Al 0 1 1 0 0 18 xaped(1).Si 0 1 1 0 0 19 xaped(1).S 0 1 1 0 0 20 xaped(1).Ar 0 1 1 0 0 21 xaped(1).Ca 0 1 1 0 0 22 xaped(1).Fe 0 1 0.5 0 0 23 xaped(1).Ni 0 1 1 0 0 sis> list_par; wabs(1) * xaped(1, He_triplets(1, He_Data)) idx param tie-to freeze value min max 1 wabs(1).nH 0 0 0.9831352 0.1 10 10^22 2 He_triplets(1).density 0 1 1e+12 1e+08 1e+16 3 xaped(1).norm1 0 0 0.001489894 1e-06 0.1 4 xaped(1).norm2 0 0 0.01380541 1e-06 0.1 5 xaped(1).temperature1 0 0 1.38829e+07 2000000 2e+07 6 xaped(1).temperature2 0 0 6.4309e+07 2e+07 2e+08 7 xaped(1).density 0 1 1 0 0 8 xaped(1).vturb 0 1 30 0 100000 9 xaped(1).redshift 0 1 0 -0.1 0.1 10 xaped(1).metal_abund 0 1 1 0 10 11 xaped(1).He 0 1 1 0 0 12 xaped(1).C 0 1 1 0 0 13 xaped(1).N 0 1 1 0 0 14 xaped(1).O 0 1 1 0 0 15 xaped(1).Ne 0 1 1 0 0 16 xaped(1).Mg 0 1 1 0 0 17 xaped(1).Al 0 1 1 0 0 18 xaped(1).Si 0 1 1 0 0 19 xaped(1).S 0 1 1 0 0 20 xaped(1).Ar 0 1 1 0 0 21 xaped(1).Ca 0 1 1 0 0 22 xaped(1).Fe 0 1 0.5 0 0 23 xaped(1).Ni 0 1 1 0 0 isis> (lo,hi) = conf("wabs(1).nH",0); isis> print([lo,hi]); 0.963041 1.0006 isis> freeze([1,3,4,5,6]); isis> .load doar21_he_modifier_fit.sl "Density 2T eval" Parameters[Variable] = 23[1] Data bins = 104 Chi-square = 84.24766 Reduced chi-square = 0.8179385 "Density 2T fit" Parameters[Variable] = 23[1] Data bins = 104 Chi-square = 84.24766 Reduced chi-square = 0.8179384 "Density 68% conf limits" **** Parameter range endpoint 1e+06 is inside the confidence limit **** Lower confidence limit didn't converge[2]: allow wider parameter ranges? 1e+06 1.88324e+13 isis> list_par; wabs(1) * xaped(1, He_triplets(1, He_Data)) idx param tie-to freeze value min max 1 wabs(1).nH 0 1 0.9831352 0.1 10 10^22 2 He_triplets(1).density 0 0 1.218506e+08 1000000 1e+16 3 xaped(1).norm1 0 1 0.001489894 1e-06 0.1 4 xaped(1).norm2 0 1 0.01380541 1e-06 0.1 5 xaped(1).temperature1 0 1 1.38829e+07 2000000 2e+07 6 xaped(1).temperature2 0 1 5.135614e+07 2e+07 2e+08 7 xaped(1).density 0 1 1 0 0 8 xaped(1).vturb 0 1 30 0 100000 9 xaped(1).redshift 0 1 0 -0.1 0.1 10 xaped(1).metal_abund 0 1 1 0 10 11 xaped(1).He 0 1 1 0 0 12 xaped(1).C 0 1 1 0 0 13 xaped(1).N 0 1 1 0 0 14 xaped(1).O 0 1 1 0 0 15 xaped(1).Ne 0 1 1 0 0 16 xaped(1).Mg 0 1 1 0 0 17 xaped(1).Al 0 1 1 0 0 18 xaped(1).Si 0 1 1 0 0 19 xaped(1).S 0 1 1 0 0 20 xaped(1).Ar 0 1 1 0 0 21 xaped(1).Ca 0 1 1 0 0 22 xaped(1).Fe 0 1 0.5 0 0 23 xaped(1).Ni 0 1 1 0 0 isis> .load doar21_he_modifier_fit.sl "Density 2T eval" Parameters[Variable] = 23[1] Data bins = 325 Cash Statistic = 236.1056 Reduced Cash Statistic = 0.7287208 "Density 2T fit" Parameters[Variable] = 23[1] Data bins = 325 Cash Statistic = 236.1055 Reduced Cash Statistic = 0.7287208 "Density 68% conf limits" **** Parameter range endpoint 1e+06 is inside the confidence limit **** Lower confidence limit didn't converge[2]: allow wider parameter ranges? 1e+06 1.15973e+13 isis> list_par; wabs(1) * xaped(1, He_triplets(1, He_Data)) idx param tie-to freeze value min max 1 wabs(1).nH 0 1 0.9831352 0.1 10 10^22 2 He_triplets(1).density 0 0 1.315674e+08 1000000 1e+16 3 xaped(1).norm1 0 1 0.001489894 1e-06 0.1 4 xaped(1).norm2 0 1 0.01380541 1e-06 0.1 5 xaped(1).temperature1 0 1 1.38829e+07 2000000 2e+07 6 xaped(1).temperature2 0 1 5.135614e+07 2e+07 2e+08 7 xaped(1).density 0 1 1 0 0 8 xaped(1).vturb 0 1 30 0 100000 9 xaped(1).redshift 0 1 0 -0.1 0.1 10 xaped(1).metal_abund 0 1 1 0 10 11 xaped(1).He 0 1 1 0 0 12 xaped(1).C 0 1 1 0 0 13 xaped(1).N 0 1 1 0 0 14 xaped(1).O 0 1 1 0 0 15 xaped(1).Ne 0 1 1 0 0 16 xaped(1).Mg 0 1 1 0 0 17 xaped(1).Al 0 1 1 0 0 18 xaped(1).Si 0 1 1 0 0 19 xaped(1).S 0 1 1 0 0 20 xaped(1).Ar 0 1 1 0 0 21 xaped(1).Ca 0 1 1 0 0 22 xaped(1).Fe 0 1 0.5 0 0 23 xaped(1).Ni 0 1 1 0 0 isis> list_par; wabs(1) * xaped(1, He_triplets(1, He_Data)) idx param tie-to freeze value min max 1 wabs(1).nH 0 0 0.9831352 0.1 10 10^22 2 He_triplets(1).density 0 1 1.315674e+08 1000000 1e+16 3 xaped(1).norm1 0 0 0.001489894 1e-06 0.1 4 xaped(1).norm2 0 0 0.01380541 1e-06 0.1 5 xaped(1).temperature1 0 0 1.38829e+07 2000000 2e+07 6 xaped(1).temperature2 0 0 5.135614e+07 2e+07 2e+08 7 xaped(1).density 0 1 1 0 0 8 xaped(1).vturb 0 1 30 0 100000 9 xaped(1).redshift 0 1 0 -0.1 0.1 10 xaped(1).metal_abund 0 1 1 0 10 11 xaped(1).He 0 1 1 0 0 12 xaped(1).C 0 1 1 0 0 13 xaped(1).N 0 1 1 0 0 14 xaped(1).O 0 1 1 0 0 15 xaped(1).Ne 0 1 1 0 0 16 xaped(1).Mg 0 1 1 0 0 17 xaped(1).Al 0 1 1 0 0 18 xaped(1).Si 0 1 1 0 0 19 xaped(1).S 0 1 1 0 0 20 xaped(1).Ar 0 1 1 0 0 21 xaped(1).Ca 0 1 1 0 0 22 xaped(1).Fe 0 1 0.5 0 0 23 xaped(1).Ni 0 1 1 0 0 isis> rebin_combined( hmc, min_counts ) ; isis> list_data; Current Spectrum List: id instrument m prt src use/nbins A R totcts exp(ksec) target 1 HETG-ACIS -1 1 1 1016/ 1016 1 1 4.5170e+03 90.990 DOAR 21 2 HETG-ACIS 1 1 1 1016/ 1016 2 2 4.2640e+03 90.990 DOAR 21 3 HETG-ACIS -1 2 1 1016/ 1016 3 3 7.3290e+03 90.990 DOAR 21 4 HETG-ACIS 1 2 1 1016/ 1016 4 4 8.3170e+03 90.990 DOAR 21 5 -1 1 1 0/ 764 5 5 8.7810e+03 90.990 DOAR 21 6 -1 2 1 0/ 1074 6 6 1.5698e+04 90.990 DOAR 21 isis> print(min_counts); 20 isis> set_fit_statistic("chisqr;sigma=data"); isis> set_fit_method("lmdif"); isis> thaw("xaped(1).Fe"); isis> fit_counts; Parameters[Variable] = 23[6] Data bins = 1016 Chi-square = 935.7205 Reduced chi-square = 0.926456 isis> list_par; wabs(1) * xaped(1, He_triplets(1, He_Data)) idx param tie-to freeze value min max 1 wabs(1).nH 0 0 0.9711006 0.1 10 10^22 2 He_triplets(1).density 0 1 1.315674e+08 1000000 1e+16 3 xaped(1).norm1 0 0 0.001472058 1e-06 0.1 4 xaped(1).norm2 0 0 0.01395153 1e-06 0.1 5 xaped(1).temperature1 0 0 1.372742e+07 2000000 2e+07 6 xaped(1).temperature2 0 0 6.283632e+07 2e+07 2e+08 7 xaped(1).density 0 1 1 0 0 8 xaped(1).vturb 0 1 30 0 100000 9 xaped(1).redshift 0 1 0 -0.1 0.1 10 xaped(1).metal_abund 0 1 1 0 10 11 xaped(1).He 0 1 1 0 0 12 xaped(1).C 0 1 1 0 0 13 xaped(1).N 0 1 1 0 0 14 xaped(1).O 0 1 1 0 0 15 xaped(1).Ne 0 1 1 0 0 16 xaped(1).Mg 0 1 1 0 0 17 xaped(1).Al 0 1 1 0 0 18 xaped(1).Si 0 1 1 0 0 19 xaped(1).S 0 1 1 0 0 20 xaped(1).Ar 0 1 1 0 0 21 xaped(1).Ca 0 1 1 0 0 22 xaped(1).Fe 0 0 0.3493926 0.1 2 23 xaped(1).Ni 0 1 1 0 0 isis> list_data; Current Spectrum List: id instrument m prt src use/nbins A R totcts exp(ksec) target 1 HETG-ACIS -1 1 1 1016/ 1016 1 1 4.5170e+03 90.990 DOAR 21 2 HETG-ACIS 1 1 1 1016/ 1016 2 2 4.2640e+03 90.990 DOAR 21 3 HETG-ACIS -1 2 1 1016/ 1016 3 3 7.3290e+03 90.990 DOAR 21 4 HETG-ACIS 1 2 1 1016/ 1016 4 4 8.3170e+03 90.990 DOAR 21 5 -1 1 1 0/ 764 5 5 8.7810e+03 90.990 DOAR 21 6 -1 2 1 0/ 1074 6 6 1.5698e+04 90.990 DOAR 21 isis> xnotice(hm,1.8,18.); isis> (lo,hi) = conf("xaped(1).Fe",0); isis> print([lo,hi]); 0.310759 0.396759 isis> set_par("xaped(1).Si",0.70,0,0.1,2.0); isis> fit_counts; Parameters[Variable] = 23[6] Data bins = 1011 Chi-square = 897.7372 Reduced chi-square = 0.8932708 0 isis> list_par; wabs(1) * xaped(1, He_triplets(1, He_Data)) idx param tie-to freeze value min max 1 wabs(1).nH 0 0 1.006647 0.1 10 10^22 2 He_triplets(1).density 0 1 1.315674e+08 1000000 1e+16 3 xaped(1).norm1 0 0 0.002423091 1e-06 0.1 4 xaped(1).norm2 0 0 0.01364595 1e-06 0.1 5 xaped(1).temperature1 0 0 1.503305e+07 2000000 2e+07 6 xaped(1).temperature2 0 0 6.500116e+07 2e+07 2e+08 7 xaped(1).density 0 1 1 0 0 8 xaped(1).vturb 0 1 30 0 100000 9 xaped(1).redshift 0 1 0 -0.1 0.1 10 xaped(1).metal_abund 0 1 1 0 10 11 xaped(1).He 0 1 1 0 0 12 xaped(1).C 0 1 1 0 0 13 xaped(1).N 0 1 1 0 0 14 xaped(1).O 0 1 1 0 0 15 xaped(1).Ne 0 1 1 0 0 16 xaped(1).Mg 0 1 1 0 0 17 xaped(1).Al 0 1 1 0 0 18 xaped(1).Si 0 0 0.7022756 0.1 2 19 xaped(1).S 0 1 1 0 0 20 xaped(1).Ar 0 1 1 0 0 21 xaped(1).Ca 0 1 1 0 0 22 xaped(1).Fe 0 1 0.3493926 0.1 2 23 xaped(1).Ni 0 1 1 0 0 isis> (lo,hi) = conf("xaped(1).Si",0); isis> print([lo,hi]); 0.64686 0.760835 isis> fit_counts; Parameters[Variable] = 23[6] Data bins = 1011 Chi-square = 896.9453 Reduced chi-square = 0.8924829 0 isis> list_par; wabs(1) * xaped(1, He_triplets(1, He_Data)) idx param tie-to freeze value min max 1 wabs(1).nH 0 0 1.010074 0.1 10 10^22 2 He_triplets(1).density 0 1 1.315674e+08 1000000 1e+16 3 xaped(1).norm1 0 0 0.002472685 1e-06 0.1 4 xaped(1).norm2 0 0 0.01364138 1e-06 0.1 5 xaped(1).temperature1 0 0 1.501804e+07 2000000 2e+07 6 xaped(1).temperature2 0 0 6.484358e+07 2e+07 2e+08 7 xaped(1).density 0 1 1 0 0 8 xaped(1).vturb 0 0 130.9176 0 3000 9 xaped(1).redshift 0 1 0 -0.1 0.1 10 xaped(1).metal_abund 0 1 1 0 10 11 xaped(1).He 0 1 1 0 0 12 xaped(1).C 0 1 1 0 0 13 xaped(1).N 0 1 1 0 0 14 xaped(1).O 0 1 1 0 0 15 xaped(1).Ne 0 1 1 0 0 16 xaped(1).Mg 0 1 1 0 0 17 xaped(1).Al 0 1 1 0 0 18 xaped(1).Si 0 1 0.7022756 0.1 2 19 xaped(1).S 0 1 1 0 0 20 xaped(1).Ar 0 1 1 0 0 21 xaped(1).Ca 0 1 1 0 0 22 xaped(1).Fe 0 1 0.3493926 0.1 2 23 xaped(1).Ni 0 1 1 0 0 isis> (lo,hi) = conf("xaped(1).vturb",0); **** Parameter range endpoint 0 is inside the confidence limit **** Lower confidence limit didn't converge[8]: allow wider parameter ranges? isis> print([lo,hi]); 0 200.654