file: fluffium_N9997.txt 1/17,18/04 -dd Make another version of Fluffium that is up to 10% more transmissive in the 0.32 to 0.48 keV range to better agree with Herman's Mkn 421 observation... Will follow the 'N9998 method, in fluffium_N9998.txt, to create data files to pass to Dave H, but first need to decide on a fluffium to use... ; - - - In ~/idl/marx/ run and adjust parameters of the twolevel model: IDL> contam_twolevel.pro Instead of the N9998 version ( which is based on a 20% thick model with transmission at 19A of ~0.92 - which then had to be raised to 1.52 power to get the desried ~0.885 transmission) Use a 15% thick-level model and get factor of 0.885 at 19A... Then use those parameters in the ISIS s/w: (iterate in ISIS so that the F-K and O-K edges don't show in the twolevel/nominal curve.) % 15%-85% 0.885 version: (*** used for N9997 ...) % the IDL contam_twolevel.pro parameters were: % thkfill = 0.15 % thkfac = 5.7 % thnfac = 0.97 % % ISIS code: set_par("aciscontam(1).OK", 0.97); set_par("aciscontam(1).FK", 0.90); fit_fun(" (0.85*(aciscontam(1))^0.97 + 0.15*(aciscontam(1))^5.7) "); in ~/dd/isis/v1.0/contam_twolevel.i to get the desired fluffium transmission at 2002.97: isis> .load contam_twolevel.040117.i Fluffium (only) transmission at t = 2002.97 1.01 1.0000 2.21 1.0000 3.40 1.0000 4.60 0.9983 5.80 0.9918 6.99 0.9835 8.19 0.9742 9.39 0.9637 10.58 0.9523 11.78 0.9407 12.97 0.9296 14.17 0.9189 15.37 0.9085 16.56 0.8994 17.76 0.8935 18.96 0.8883 20.15 0.8836 ... etc ... 39.29 0.9007 40.49 0.9091 41.69 0.9204 42.88 0.9333 44.08 0.9351 45.27 0.9313 46.47 0.9277 47.67 0.9242 48.86 0.9210 49.99 0.9181 and put these values into: acis_contam_fluf15thk.tbl Plot this fluffium version along with the previous '20thk version, in IDL... IDL> readcol, 'acis_contam_fluf20thk.tbl', fwav, f20 IDL> f20 = f20^1.52 IDL> readcol, 'acis_contam_fluf15thk.tbl', fwav, f15 IDL> plot, fwav, f20, yrange=[0.8,1.03],/ysty,xrange=[0.,65.0],/xsty,linestyle=1 IDL> oplot, fwav, f15,linestyle=2 The plot above shows a steep kinkk and drop in the 40A to 35A range : Because Herman has adjusted the Carbon near-edge structure in his model he has essentially taken out some amount of "fluffium" close to the edge, say in the 0.28 to 0.5 keV range... ... and in particular this kink and drop are likely already in his Carbon to some degree... So, level things off a little by using a linear interpolation from the Carbon edge to the minimum near 26A - spreading the change out in a less structured way... Create from '15thk the file: *** New fluffium to use: * * acis_contam_fluf15thkMOD.tbl * *** Now read in the MOD version and plot these three versions: IDL> readcol, 'acis_contam_fluf15thkMOD.tbl', fwav, f15mod IDL> plot, fwav, f20, yrange=[0.8,1.03],/ysty,xrange=[0.,65.0],/xsty,linestyle=1 IDL> oplot, fwav, f15,linestyle=2 IDL> oplot, fwav, f15mod Put this plot in: contam_f20f15f15MOD.ps and with log(E) on axis: contam_f20f15f15MOD_logE.ps It looks pretty kinky - but these are changes at the percent level... Make a plot showing the ratio of this new f15mod with f20 and put it on the same scale as Herman's 'Fractional Residual" plot for the Mkn 421 ~2002.8 data: IDL> plot, !ddhc/fwav, f15mod/f20 - 1.0, xrange=[0.2,0.5],/xsty, yrange=[-0.2,0.2],/ysty IDL> oplot, [0.2,0.5],[0.,0.],linestyle=1 put it in: contam_15thkMOD_20thk.ps OK, so that's the modified fluffium to use... and fine tune for ECS agreement. ------------------ Now continue as in the N9998 creation process - summarized here: Read in this Fluffium file as well as the "abs" file and write out a new version which has an extra Fluffium column - these steps were done in IDL as documented in: contam_add_fluffium.040118.pro Manually pre-pend header info to the IDL-written file and now we have the new file: acis_contam_abs_v4.tbl * Note version 4 is being used to match the version 4 * that will be created for the time dependance... Which looks like: # Created by DD on 1/18/04 # # file: acis_contam_abs_v4.tbl # # This file was created by adding a Fluffium column # to the existing file # acis_contam_abs_v1.tbl . # The Fluffium transmission was taken from the file: # acis_contam_fluf15thkMOD.tbl and converted to optical depths. # # See also: fluffium_N9997.txt and contam_add_fluffium.040118.pro . # # Optical depths due to each component # # wavelength C-K O-K F-K Fluffium 1.01250 0.00000 0.00000 0.00000 1.00000e-05 1.03750 0.00000 0.00000 0.00000 1.00000e-05 1.06250 0.00000 0.00000 0.00000 1.00000e-05 1.08750 0.00000 0.00000 0.00000 1.00000e-05 ... 11.9625 0.0843226 0.0186418 0.0222370 0.0629514 11.9875 0.0848125 0.0187447 0.0222988 0.0632008 12.0125 0.0853071 0.0188484 0.0223597 0.0634501 12.0375 0.0858019 0.0189503 0.0224196 0.0636995 ... 18.9875 0.292613 0.0603872 0.00365866 0.118568 19.0125 0.293583 0.0605962 0.00367481 0.118680 ... 56.2625 0.165344 0.0485167 0.0698523 0.104932 56.2875 0.165530 0.0485783 0.0699265 0.105010 56.3125 0.165716 0.0486399 0.0700007 0.105088 56.3375 0.165903 0.0487016 0.0700750 0.105165 - - - - - Optical depth at ~0.7 keV - - - - - - The Total optical depth "at ~0.7 keV" for ECS purposes is the weighted average of optical depths at the three line-group wavelengths, which for S3 Alexey gives as: Line group E QE_S3/QE_I3 S3 Flux fraction --------------------------------------------------------------- Fe alpha, beta 0.706 2.43152 41% Mn alpha, beta 0.638 3.25435 53% Fe zeta, eta 0.618 3.69358 6% --------------------------------------------------------------- Or tabulated wavelengths of: 17.5625, 18.1625, 20.0625 17.5625 0.241522 0.0499158 0.0888599 0.111526 18.1625 0.261969 0.0539561 0.00315442 0.114567 20.0625 0.336250 0.0692934 0.00439601 0.123361 For this new with-Fluffium material we have (at t = 2002.82 where the time file gives additional factors: 0.9846 1.1099 0.9157 and assuming a fluffium factor of 1.0) > print, TOTAL( [0.9846, 1.1099, 0.9157, 1.0] * [0.241522, 0.0499158, 0.0888599, 0.111526] ) 0.486099 > print, TOTAL( [0.9846, 1.1099, 0.9157, 1.0] * [0.261969 , 0.0539561 , 0.00315442 , 0.114567] ) 0.435276 > print, TOTAL( [0.9846, 1.1099, 0.9157, 1.0] * [0.336250 , 0.0692934 , 0.00439601 , 0.123361] ) 0.535367 converting to TRANSMISSION and weighting these: 17.5625 : exp(-0.486099) * 0.41 = 0.252159 18.1625 : exp(-0.435276) * 0.53 = 0.342956 20.0625 : exp(-0.535367) * 0.06 = 0.0351273 ---------- weighted trans = 0.630242 How does this compare with the Tennant equation prediction? Using the code in contam_add_fluffium.pro Tennant at 2002.82 gives tau: > print, tau_ref 0.460100 (trans = 0.631221) and the ratio if transmissions is: > print, 0.630242/0.631221 0.998449 So, this is within 0.2 % of the target goal. check. -------------------------------------- OK, now the time dependance.... We want a new version of acis_contam_tau_t_v2_2006.tbl which will have the 2002.82 line read: 2002.820 0.9846 1.1099 0.9157 1.0000 and which will be scaled from these values by the Tennant tau trend equation... So do that in contam_add_fluffium.pro and write out the new file: acis_contam_tau_t_v4.tbl and add header info to this IDL-created file... et voila. Send the 'v4 files to Dave H. - - - - -