Subject: 6-7 A region, DoAr 21, MEG +1 From: Marc Gagne Date: Fri, 11 Jul 2008 17:54:50 -0400 To: David Cohen In this little ISIS code stub, I consider only the MEG +1 data (unbinned) from 6-7 A, then exclude 0.034 A around the Si XIII r line – the r line is over predicted by the APEC model. Using this data only, we find 68% upper and lower limits on ne: 1.5358e+13 and 1.39124e+14. The 90% upper limit: 2.42819e+14. The temperature model from the S XIV/XIII region is very close to the model from the global fit using all data. isis> ignore(4,6.632,6.666); isis> eval_counts; Parameters[Variable] = 23[2] Data bins = 146 Chi-square = 94.76048 Reduced chi-square = 0.6580589 0 isis> print("Density 68% conf limits"); "Density 68% conf limits" isis> variable lo,hi; isis> (lo,hi) = conf("He_triplets(1).density",0); isis> print([lo,hi]); 1.5358e+13 1.39124e+14 isis> (lo,hi) = conf("He_triplets(1).density",1); **** Parameter range endpoint 1 is inside the confidence limit **** Lower confidence limit didn't converge[2]: allow wider parameter ranges? isis> print([lo,hi]); 1 2.42819e+14 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.98 0.0001 100 10^22 2 He_triplets(1).density 0 0 6.319202e+13 1 1e+16 3 xaped(1).norm1 0 0 0.002901925 0 0 4 xaped(1).norm2 0 1 0.0125 0 0 5 xaped(1).temperature1 0 1 1e+07 0 0 6 xaped(1).temperature2 0 1 6.4e+07 0 0 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> doar21_SiXIII_2t_he_modifier.ps Content-Type: application/octet-stream Content-Encoding: base64