% zo_1d_analysis.i % % Read in the isis-format zo 1D histogram files % from Capella, Mkn 421, and MARX capella simulation % and scale and plot the distributions. % 1/5/00 dd % % Generalize to study focus/LRF vs time... % 1/12/00 dd plot_unit("Angstrom"); Remove_Spectrum_Gaps=1; delete_data([1:6]); % get_data variables: variable dds_obs_out_dir, data_set_name; % - - - - - % more variables ... % Simulation variables: variable simm, simp, csimm, csimp; % Variables to scale the observations, etc. variable cap_fact, sim_fact, nom_streak_fact, streak_fact; % Nominal streak factor... % Normalize the streak to the same scale as the zeroth-order % Ratio of zeroth-order exposure time to streak exposure time, % include fact that the HAK streak excludes ~240 rows: streak_fact = 3200.0 / ( (41.0/1024.0) * (1022.-240.0) ) ; nom_streak_fact = 3200.0 / ( (41.0/1024.0) * (1022.-240.0) ) ; % Modified calculation for 1.7s and streak 165 rows accepted: %%streak_fact = 1700.0 / (41.0 * ((165.0)/1022.0)) ; % variables to keep average background count in variable background; % other variables variable data_xtitle; % and pixel scales variable meg_angpix, heg_angpix, use_angpix; meg_angpix = 0.011124491; heg_angpix = 0.0055625373; % - - - - - % Default selection of counts or flux here... %% use_counts; % Load some zero-order profiles: % % For comparison: % Capella 1318 zero-order; 500 um analysis % Capella 1318 streak ()=load_data("/spectra/d2/GtoHAK/Obs_out/" + "obsid_1318_Sky_zo_1d_isis.dat"); ()=load_data("/spectra/d2/GtoHAK/Obs_out/" + "obsid_1318_Sky_strk_1d_isis.dat"); %****************************************** % Observation under measurement: % TBD test obs zero-order % TBD test obs streak % Set cap_fact so that Observation and Capella % zeroth-order wings agree %****************************************** % Data sets: % - - - %%data_set_name = "Capella (obsid 1235)"; %%()=load_data("/spectra/d2/GtoHAK/Obs_out/" + %% "obsid_1235_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/GtoHAK/Obs_out/" + %% "obsid_1235_Sky_strk_1d_isis.dat"); %%cap_fact = 0.5; % - - - %%data_set_name = "HR1099 (obsid 62538)"; %%()=load_data("/spectra/d2/GtoHAK/Obs_out/" + %% "obsid_62538_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/GtoHAK/Obs_out/" + %% "obsid_62538_Sky_strk_1d_isis.dat"); %%cap_fact = 6.0; % - - - %%data_set_name = "HR1099 (obsid 1252)"; %%()=load_data("/spectra/d2/GtoHAK/Obs_out/" + %% "obsid_1252_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/GtoHAK/Obs_out/" + %% "obsid_1252_Sky_strk_1d_isis.dat"); %%cap_fact = 0.9; % - - - %%data_set_name = "SS 433 (obsid 106)"; %%()=load_data("/spectra/d2/GtoHAK/Obs_out/" + %% "obsid_106_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/GtoHAK/Obs_out/" + %% "obsid_106_Sky_strk_1d_isis.dat"); %%cap_fact = 0.8; % - - - % 500 um CD %%data_set_name = "Capella (obsid 1103)"; %%()=load_data("/spectra/d2/GtoHAK/Obs_out/" + %% "obsid_1103_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/GtoHAK/Obs_out/" + %% "obsid_1103_Sky_strk_1d_isis.dat"); %%cap_fact = 1.4; % - - - % 500 um CD %%data_set_name = "Capella (obsid 1318)"; %%()=load_data("/spectra/d2/GtoHAK/Obs_out/" + %% "obsid_1318_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/GtoHAK/Obs_out/" + %% "obsid_1318_Sky_strk_1d_isis.dat"); %%cap_fact = 1.01; % - - - %%data_set_name = "AB Dor (obsid 16)"; %%()=load_data("/spectra/d2/Fmeo/hak_16/Obs_out/" + %% "obsid_16_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/Fmeo/hak_16/Obs_out/" + %% "obsid_16_Sky_strk_1d_isis.dat"); %%cap_fact = 0.8; % - - - %%data_set_name = "Q 0836 (obsid 1450)"; %%()=load_data("/spectra/d2/GtoHAK/Obs_out/" + %% "obsid_1450_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/GtoHAK/Obs_out/" + %% "obsid_1450_Sky_strk_1d_isis.dat"); %%cap_fact = 0.5; % - - - %%data_set_name = "II Peg (obsid 1451)"; %%()=load_data("/spectra/d2/GtoHAK/Obs_out/" + %% "obsid_1451_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/GtoHAK/Obs_out/" + %% "obsid_1451_Sky_strk_1d_isis.dat"); %%cap_fact = 1.3; % - - - %%data_set_name = "4U 1636 (obsid 105)"; %%()=load_data("/spectra/d2/GtoHAK/Obs_out/" + %% "obsid_105_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/GtoHAK/Obs_out/" + %% "obsid_105_Sky_strk_1d_isis.dat"); %%cap_fact =3.0; % - - - %%data_set_name = "Cyg X-3 (obsid 1456)"; %%()=load_data("/spectra/d2/Fmeo/hak_1456/Obs_out/" + %% "obsid_1456_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/Fmeo/hak_1456/Obs_out/" + %% "obsid_1456_Sky_strk_1d_isis.dat"); %%cap_fact = 4.0; % - - - %%data_set_name = "Mkn 421 (obsid 457)"; %%()=load_data("/spectra/d2/GtoHAK/Obs_out/" + %% "obsid_457_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/GtoHAK/Obs_out/" + %% "obsid_457_Sky_strk_1d_isis.dat"); %%cap_fact = 1.0/3.8; % - - - %%data_set_name = "PKS 2149-306 (obsid 1481)"; %%()=load_data("/spectra/d2/GtoHAK/Obs_out/" + %% "obsid_1481_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/GtoHAK/Obs_out/" + %% "obsid_1481_Sky_strk_1d_isis.dat"); %%cap_fact = 0.2; % - - - %%data_set_name = "Lambda And (obsid 609)"; %%()=load_data("/spectra/d2/Fmeo/hak_609/Obs_out/" + %% "obsid_609_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/Fmeo/hak_609/Obs_out/" + %% "obsid_609_Sky_strk_1d_isis.dat"); %%cap_fact = 1.0; % - - - %%data_set_name = "Cyg X-3,II (obsid 1456b)"; %%()=load_data("/spectra/d2/Fmeo/hak_1456p2/Obs_out/" + %% "obsid_1456_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/Fmeo/hak_1456p2/Obs_out/" + %% "obsid_1456_Sky_strk_1d_isis.dat"); %%cap_fact =2.0; % - - - %%data_set_name = "3C 273 (obsid 459)"; %%()=load_data("/spectra/d2/GtoHAK/Obs_out/" + %% "obsid_459_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/GtoHAK/Obs_out/" + %% "obsid_459_Sky_strk_1d_isis.dat"); %%cap_fact = 2.0; % - - - %%data_set_name = "Cyg X-1 (obsid 1511)"; %%()=load_data("/spectra/d2/GtoHAK/Obs_out/" + %% "obsid_1511_Sky_zo_1d_isis.dat"); % no streak (all streak?!) because CC mode %%()=load_data("/spectra/d2/GtoHAK/Obs_out/" + %% "obsid_1511_Sky_zo_1d_isis.dat"); %%cap_fact = 5.0; % - - - %%data_set_name = "UX Ari (obsid 605)"; %%()=load_data("/spectra/d2/Fmeo/hak_605/Obs_out/" + %% "obsid_605_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/Fmeo/hak_605/Obs_out/" + %% "obsid_605_Sky_strk_1d_isis.dat"); %%cap_fact = 0.8; % - - - %%data_set_name = "delta Ori (obsid 639)"; %%()=load_data("/spectra/d2/Fmeo/hak_639/Obs_out/" + %% "obsid_639_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/Fmeo/hak_639/Obs_out/" + %% "obsid_639_Sky_strk_1d_isis.dat"); %%cap_fact = 0.05; % - - - %%data_set_name = "NGC 3783 (obsid 373)"; %%()=load_data("/spectra/d2/Fmeo/hak_373/Obs_out/" + %% "obsid_373_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/Fmeo/hak_373/Obs_out/" + %% "obsid_373_Sky_strk_1d_isis.dat"); %%cap_fact = 1.0; % - - - %%data_set_name = "NGC 5548 (obsid 837)"; %%()=load_data("/spectra/d2/Fmeo/hak_837/Obs_out/" + %% "obsid_837_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/Fmeo/hak_837/Obs_out/" + %% "obsid_837_Sky_strk_1d_isis.dat"); %%cap_fact = 0.8; % - - - %%data_set_name = "LMC X-3 (obsid 92,N002)"; %%()=load_data("/spectra/d2/Fmeo/hak_92/Obs_out/" + %% "obsid_92N002_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/Fmeo/hak_92/Obs_out/" + %% "obsid_92N002_Sky_strk_1d_isis.dat"); %%cap_fact = 1.0; % - - - %%data_set_name = "Cir X-1 (obsid 706)"; % NOTE: exposures are 1.7s; streak is from 165 rows % For obsid 706 frame=1.7s; streak rows = 165: %%streak_fact = 1700.0 / ( (41.0/1024.0) * (165.0) ) ; % %%()=load_data("/spectra/d2/Fmeo/hak_706/Obs_out/" + %% "obsid_706_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/Fmeo/hak_706/Obs_out/" + %% "obsid_706_Sky_strk_1d_isis.dat"); %%cap_fact = 8.0; % - - - % 500 um CD %%data_set_name = "Capella (obsid 57)"; %%()=load_data("/spectra/d2/GtoHAK/Obs_out/" + %% "obsid_57_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/GtoHAK/Obs_out/" + %% "obsid_57_Sky_strk_1d_isis.dat"); %%cap_fact = 1.0; % - - - %%data_set_name = "Cen X-3 (obsid 705)"; %%()=load_data("/spectra/d2/Fmeo/hak_705/Obs_out/" + %% "obsid_705_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/Fmeo/hak_705/Obs_out/" + %% "obsid_705_Sky_strk_1d_isis.dat"); %%cap_fact = 1.0; % - - - %%data_set_name = "NGC 4151 (obsid 335)"; %%()=load_data("/spectra/d2/GtoHAK/Obs_out/" + %% "obsid_335_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/GtoHAK/Obs_out/" + %% "obsid_335_Sky_strk_1d_isis.dat"); %%cap_fact = 0.5; % - - - %%data_set_name = "HD 68273 (obsid 629)"; %%()=load_data("/spectra/d2/Fmeo/hak_629/Obs_out/" + %% "obsid_629_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/Fmeo/hak_629/Obs_out/" + %% "obsid_629_Sky_strk_1d_isis.dat"); %%cap_fact = 0.3; % - - - %%data_set_name = "Zeta Pup (obsid 640)"; %%()=load_data("/spectra/d2/Fmeo/hak_640/Obs_out/" + %% "obsid_640_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/Fmeo/hak_640/Obs_out/" + %% "obsid_640_Sky_strk_1d_isis.dat"); %%cap_fact = 0.3; % - - - %%data_set_name = "Algol (obsid 604)"; %%()=load_data("/spectra/d2/Fmeo/hak_604/Obs_out/" + %% "obsid_604_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/Fmeo/hak_604/Obs_out/" + %% "obsid_604_Sky_strk_1d_isis.dat"); %%cap_fact = 2.0; % - - - %%data_set_name = "Cyg X-3 TOO (obsid 425)"; %%()=load_data("/spectra/d2/Fmeo/hak_425/Obs_out/" + %% "obsid_425_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/Fmeo/hak_425/Obs_out/" + %% "obsid_425_Sky_strk_1d_isis.dat"); %%cap_fact = 3.0; % - - - %%data_set_name = "GRS 1915+105 (obsid 660)"; %%()=load_data("/spectra/d2/Fmeo/hak_660/Obs_out/" + %% "obsid_660_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/Fmeo/hak_660/Obs_out/" + %% "obsid_660_Sky_strk_1d_isis.dat"); %%cap_fact = 3.0; % - - - %%data_set_name = "EX Hydra (obsid 1706)"; %%()=load_data("/spectra/d2/Fmeo/hak_1706/Obs_out/" + %% "obsid_1706_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/Fmeo/hak_1706/Obs_out/" + %% "obsid_1706_Sky_strk_1d_isis.dat"); %%cap_fact = 2.0; % - - - %%data_set_name = "Mkn 421 (b) (obsid 1714)"; %%()=load_data("/spectra/d2/Fmeo/hak_1714/Obs_out/" + %% "obsid_1714_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/Fmeo/hak_1714/Obs_out/" + %% "obsid_1714_Sky_strk_1d_isis.dat"); %%cap_fact = 3.0; % - - - %%data_set_name = "4U 1735-44 (obsid 704)"; %%()=load_data("/spectra/d2/Fmeo/hak_704/Obs_out/" + %% "obsid_704_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/Fmeo/hak_704/Obs_out/" + %% "obsid_704_Sky_strk_1d_isis.dat"); %%cap_fact = 3.0; % - - - %%data_set_name = "3C 273 no HETG (obsid 1712)"; %%()=load_data("/spectra/d2/Fmeo/hak_1712/Obs_out/" + %% "obsid_1712_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/Fmeo/hak_1712/Obs_out/" + %% "obsid_1712_Sky_strk_1d_isis.dat"); %%cap_fact = 10.0; % - - - %%data_set_name = "TZ CRB (obsid 15)"; %%()=load_data("/spectra/d2/Fmeo/hak_15/Obs_out/" + %% "obsid_15_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/Fmeo/hak_15/Obs_out/" + %% "obsid_15_Sky_strk_1d_isis.dat"); %%cap_fact = 2.0; % - - - % So weak, use the Zero-order in place of streak... for fun only: %%data_set_name = "TW Hydra (obsid 5)"; %%()=load_data("/spectra/d2/Fmeo/hak_5/Obs_out/" + %% "obsid_5_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/Fmeo/hak_5/Obs_out/" + %% "obsid_5_Sky_zo_1d_isis.dat"); %%cap_fact = 0.1; % - - - % HAK was run by MTA folks - copied from web page %%data_set_name = "4U 1822-37 (obsid 671)"; %%()=load_data("/spectra/d2/Fmeo/hak_671/Obs_out/" + %% "obsid_0671_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/Fmeo/hak_671/Obs_out/" + %% "obsid_0671_Sky_strk_1d_isis.dat"); %%cap_fact = 4.0; % - - - %%data_set_name = "SS Cygni (obsid 646)"; %%()=load_data("/spectra/d2/Fmeo/hak_646/Obs_out/" + %% "obsid_646Elt3_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/Fmeo/hak_646/Obs_out/" + %% "obsid_646Elt3_Sky_strk_1d_isis.dat"); %%cap_fact = 1.5; % - - - %%data_set_name = "X1543-62 (obsid 702)"; %%()=load_data("/spectra/d2/Fmeo/hak_702/Obs_out/" + %% "obsid_702_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/Fmeo/hak_702/Obs_out/" + %% "obsid_702_Sky_strk_1d_isis.dat"); %%cap_fact = 5.0; % - - - %%data_set_name = "X1543-62 (obsid 702, E<3keV)"; %%()=load_data("/spectra/d2/Fmeo/hak_702/Obs_out/" + %% "obsid_702Elt3_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/Fmeo/hak_702/Obs_out/" + %% "obsid_702Elt3_Sky_strk_1d_isis.dat"); %%cap_fact = 4.0; % - - - %%data_set_name = "AU Mic (obsid 17)"; %%()=load_data("/spectra/d2/Fmeo/hak_17/Obs_out/" + %% "obsid_17_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/Fmeo/hak_17/Obs_out/" + %% "obsid_17_Sky_strk_1d_isis.dat"); %%cap_fact = 0.3; % - - - %%data_set_name = "Eta Carinae (obsid 632)"; %%()=load_data("/spectra/d2/Fmeo/hak_632/Obs_out/" + %% "obsid_632_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/Fmeo/hak_632/Obs_out/" + %% "obsid_632_Sky_strk_1d_isis.dat"); %%cap_fact = 1.0; % - - - % This is extended - here to see if pileup % is noticable in the core of the image... %%data_set_name = "NGC 1068 (obsid 332)"; %%()=load_data("/spectra/d2/Fmeo/hak_332/Obs_out/" + %% "obsid_332_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/Fmeo/hak_332/Obs_out/" + %% "obsid_332_Sky_strk_1d_isis.dat"); %%cap_fact = 1.0; % - - - % This has reduced cross-disp width... % 350 rows total and 1.3 s exposures %%streak_fact = 1300.0 / ( (41.0/1024.0) * (350.-240.0) ) ; %%data_set_name = "Cen X-3 (obsid 1943)"; %%()=load_data("/spectra/d2/Fmeo/hak_1943/Obs_out/" + %% "obsid_1943N500k_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/Fmeo/hak_1943/Obs_out/" + %% "obsid_1943N500k_Sky_strk_1d_isis.dat"); %%cap_fact = 2.0; % E < 2 keV : %%streak_fact = 1300.0 / ( (41.0/1024.0) * (350.-240.0) ) ; %%data_set_name = "Cen X-3 (obsid 1943, E < 2keV )"; %%()=load_data("/spectra/d2/Fmeo/hak_1943/Obs_out/" + %% "obsid_1943Elt2_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/Fmeo/hak_1943/Obs_out/" + %% "obsid_1943Elt2_Sky_strk_1d_isis.dat"); %%cap_fact = 4.5; % - - - % VELA X-1 % E < 2 keV : %%streak_fact = 1700.0 / ( (41.0/1024.0) * (512.-240.0) ) ; %%data_set_name = "Vela X-1 (obsid 1928, E < 2keV )"; %%()=load_data("/spectra/d2/Fmeo/hak_1928/Obs_out/" + %% "obsid_1928Elt2_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/Fmeo/hak_1928/Obs_out/" + %% "obsid_1928Elt2_Sky_strk_1d_isis.dat"); %%cap_fact = 1.5; % - - - % Capella %%data_set_name = "Capella (obsid 1010)"; %%()=load_data("/spectra/d2/GtoHAK/Obs_out/" + %% "obsid_1010_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/GtoHAK/Obs_out/" + %% "obsid_1010_Sky_strk_1d_isis.dat"); %%cap_fact = 1.0; % - - - % SS 433 %%streak_fact = 2500.0 / ( (41.0/1024.0) * (773.-240.0) ) ; %%data_set_name = "SS 433 (obsid 1019)"; %%()=load_data("/spectra/d2/Fmeo/hak_1019/Obs_out/" + %% "obsid_1019_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/Fmeo/hak_1019/Obs_out/" + %% "obsid_1019_Sky_strk_1d_isis.dat"); %%cap_fact = 0.5; % - - - % VELA X-1 % E < 2 keV : %%streak_fact = 1700.0 / ( (41.0/1024.0) * (512.-240.0) ) ; %%data_set_name = "GRS 1758-258 (obsid 2429, E < 2keV )"; %%()=load_data("/spectra/d2/Fmeo/hak_2429/Obs_out/" + %% "obsid_2429Elt2_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/Fmeo/hak_2429/Obs_out/" + %% "obsid_2429Elt2_Sky_strk_1d_isis.dat"); %%cap_fact = 0.8; % - - - %%data_set_name = "3C 273 (obsid 2463)"; %%()=load_data("/spectra/d2/GtoHAK/Obs_out/" + %% "obsid_2463_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/GtoHAK/Obs_out/" + %% "obsid_2463_Sky_strk_1d_isis.dat"); %%streak_fact = 1800.0 / (41.0 * ((541.-240.0)/1022.0)) ; %%cap_fact = 1.4; % - - - %%data_set_name = "GX340+0 (obsid 1921)"; %%()=load_data("/spectra/d2/Fmeo/hak_1921/Obs_out/" + %% "obsid_1921_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/Fmeo/hak_1921/Obs_out/" + %% "obsid_1921_Sky_strk_1d_isis.dat"); %%streak_fact = 1400.0 / (41.0 * ((399.-240.0)/1022.0)) ; %%cap_fact = 19.0; % - - - %%data_set_name = "GX340+0 (obsid 1921, E < 4 keV)"; %%()=load_data("/spectra/d2/Fmeo/hak_1921/Obs_out/" + %% "obsid_1921Elt4_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/Fmeo/hak_1921/Obs_out/" + %% "obsid_1921Elt4_Sky_strk_1d_isis.dat"); %%streak_fact = 1400.0 / (41.0 * ((399.-240.0)/1022.0)) ; %%cap_fact = 7.0; % - - - %%data_set_name = "PKS 2155-304 (obsid 3167)"; %%()=load_data("/spectra/d2/Fmeo/hak_3167/Obs_out/" + %% "obsid_3167L1_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/Fmeo/hak_3167/Obs_out/" + %% "obsid_3167L1_Sky_strk_1d_isis.dat"); %%streak_fact = nom_streak_fact; %%cap_fact = 2.0; % - - - %%data_set_name = "Cyg X-1 (obsid 2741, E < 3 keV)"; %%()=load_data("/spectra/d2/Fmeo/hak_2741/Obs_out/" + %% "obsid_2741L1_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/Fmeo/hak_2741/Obs_out/" + %% "obsid_2741L1_Sky_strk_1d_isis.dat"); % Ratio of zeroth-order exposure time to streak exposure time, % include fact that the HAK streak excludes ~240 rows: %%streak_fact = 1700.0 / ( (41.0/1024.0) * (512.-240.0) ) ; %%cap_fact = 14.0; % - - - %%data_set_name = "Capella (obsid 2583)"; %%()=load_data("/spectra/d2/GtoHAK/Obs_out/" + %% "obsid_2583_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/GtoHAK/Obs_out/" + %% "obsid_2583_Sky_strk_1d_isis.dat"); % Ratio of zeroth-order exposure time to streak exposure time, % include fact that the HAK streak excludes ~240 rows: %%streak_fact = 2500.0 / ( (41.0/1024.0) * (774.-240.0) ) ; %%cap_fact = 1.0; % - - - %%data_set_name = "GS 1826 (obsid 2739, E < 3 keV)"; %%()=load_data("/spectra/d2/Fmeo/hak_2739/Obs_out/" + %% "obsid_2739Elt3_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/Fmeo/hak_2739/Obs_out/" + %% "obsid_2739Elt3_Sky_strk_1d_isis.dat"); % Ratio of zeroth-order exposure time to streak exposure time, % include fact that the HAK streak excludes ~240 rows: %%streak_fact = 1700.0 / ( (41.0/1024.0) * (542.-240.0) ) ; %%cap_fact = 8.0; % - - - %%data_set_name = "GX 13+1 (obsid 2708, E < 3 keV)"; %%()=load_data("/spectra/d2/Fmeo/hak_2708/Obs_out/" + %% "obsid_2708Elt3_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/Fmeo/hak_2708/Obs_out/" + %% "obsid_2708Elt3_Sky_strk_1d_isis.dat"); % Ratio of zeroth-order exposure time to streak exposure time, % include fact that the HAK streak excludes ~240 rows: %%streak_fact = 1200.0 / ( (41.0/1024.0) * (350.-240.0) ) ; %%cap_fact = 8.0; % - - - %%data_set_name = "GX 339-4 (obsid 4420, E < 3 keV)"; %%()=load_data("/spectra/d2/Fmeo/hak_4420/Obs_out/" + %% "obsid_4420Elt3_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/Fmeo/hak_4420/Obs_out/" + %% "obsid_4420Elt3_Sky_strk_1d_isis.dat"); % Ratio of zeroth-order exposure time to streak exposure time, % include fact that the HAK streak excludes ~240 rows: %%streak_fact = 1400.0 / ( (41.0/1024.0) * (400.-240.0) ) ; %%cap_fact = 8.0; % - - - %%data_set_name = "Cyg X-1 (obsid 3814, E < 3 keV)"; %%()=load_data("/spectra/d2/Fmeo/hak_3814/Obs_out/" + %% "obsid_3814Elt3_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/Fmeo/hak_3814/Obs_out/" + %% "obsid_3814Elt3_Sky_strk_1d_isis.dat"); % Ratio of zeroth-order exposure time to streak exposure time, % include fact that the HAK streak excludes ~240 rows: %%streak_fact = 1700.0 / ( (41.0/1024.0) * (512.-240.0) ) ; %%cap_fact = 7.0; % - - - %%data_set_name = "Her X-1 [no HETG] (obsid 3814, E < 3 keV)"; %%()=load_data("/spectra/d2/Fmeo/hak_3662/Obs_out/" + %% "obsid_3662Elt3_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/Fmeo/hak_3662/Obs_out/" + %% "obsid_3662Elt3_Sky_strk_1d_isis.dat"); % Ratio of zeroth-order exposure time to streak exposure time, % include fact that the HAK streak excludes ~240 rows: %%streak_fact = 3100.0 / ( (41.0/1024.0) * (1024.-240.0) ) ; %%cap_fact = 150.0; % - - - %%data_set_name = "4U 1626-67 (obsid 3504, E < 3 keV)"; %%()=load_data("/spectra/d2/Fmeo/hak_3504/Obs_out/" + %% "obsid_3504Elt3_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/Fmeo/hak_3504/Obs_out/" + %% "obsid_3504Elt3_Sky_strk_1d_isis.dat"); % Ratio of zeroth-order exposure time to streak exposure time, % include fact that the HAK streak excludes ~240 rows: %%streak_fact = 1700.0 / ( (41.0/1024.0) * (512.-240.0) ) ; %%cap_fact = 1.0; % - - - %%data_set_name = "3C 273 (obsid 4430)"; %%()=load_data("/spectra/d2/Fmeo/hak_4430/Obs_out/" + %% "obsid_4430Elt3_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/Fmeo/hak_4430/Obs_out/" + %% "obsid_4430Elt3_Sky_strk_1d_isis.dat"); % Ratio of zeroth-order exposure time to streak exposure time, % include fact that the HAK streak excludes ~240 rows: %%streak_fact = 2500.0 / ( (41.0/1024.0) * (774.-240.0) ) ; %%cap_fact = 1.5; % - - - %%data_set_name = "XTE J1814 (obsid 4439, E < 3 keV)"; %%% *** Not a point source ?!? *** %%()=load_data("/spectra/d2/Fmeo/hak_4439/Obs_out/" + %% "obsid_4439Elt3_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/Fmeo/hak_4439/Obs_out/" + %% "obsid_4439Elt3_Sky_strk_1d_isis.dat"); % Ratio of zeroth-order exposure time to streak exposure time, % include fact that the HAK streak excludes ~240 rows: %%streak_fact = 3200.0 / ( (41.0/1024.0) * (1024.-240.0) ) ; %%cap_fact = 0.33; % - - - %%data_set_name = "NGC 2110 (obsid 3418)"; %%% Do NOT use for FWHM *** Not a point source ! *** %%()=load_data("/spectra/d2/Fmeo/hak_3418/Obs_out/" + %% "obsid_3418_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/Fmeo/hak_3418/Obs_out/" + %% "obsid_3418_Sky_strk_1d_isis.dat"); % Ratio of zeroth-order exposure time to streak exposure time, % include fact that the HAK streak excludes ~240 rows: %%streak_fact = 3200.0 / ( (41.0/1024.0) * (774.-240.0) ) ; %%cap_fact = 0.33; % - - - %%data_set_name = "Capella (obsid 3674)"; %%()=load_data("/spectra/d2/GtoHAK/Obs_out/" + %% "obsid_3674_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/GtoHAK/Obs_out/" + %% "obsid_3674_Sky_strk_1d_isis.dat"); % Ratio of zeroth-order exposure time to streak exposure time, % include fact that the HAK streak excludes ~240 rows: %%streak_fact = 2500.0 / ( (41.0/1024.0) * (774.-240.0) ) ; %%cap_fact = 1.0; % - - - %%data_set_name = "EXO 0748 (obsid 4574)"; %%()=load_data("/spectra/d2/GtoHAK/Obs_out/" + %% "obsid_4574Elt3_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/GtoHAK/Obs_out/" + %% "obsid_4574Elt3_Sky_strk_1d_isis.dat"); % Ratio of zeroth-order exposure time to streak exposure time, % include fact that the HAK streak excludes ~240 rows: %%streak_fact = 1700.0 / ( (41.0/1024.0) * (542.-240.0) ) ; %%cap_fact = 1.1; % - - - %%data_set_name = "4U 1254 (obsid 3823)"; %%()=load_data("/spectra/d2/GtoHAK/Obs_out/" + %% "obsid_3823Elt3_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/GtoHAK/Obs_out/" + %% "obsid_3823Elt3_Sky_strk_1d_isis.dat"); % Ratio of zeroth-order exposure time to streak exposure time, % include fact that the HAK streak excludes ~240 rows: %%streak_fact = 1700.0 / ( (41.0/1024.0) * (512.-240.0) ) ; %%cap_fact = 6.0; % - - - %%data_set_name = "AC211 (double, obsid 4572)"; %%()=load_data("/spectra/d2/Fmeo/hak_4572/Obs_out/" + %% "obsid_4572_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/Fmeo/hak_4572/Obs_out/" + %% "obsid_4572_Sky_strk_1d_isis.dat"); % Ratio of zeroth-order exposure time to streak exposure time, % include fact that the HAK streak excludes ~240 rows: %%streak_fact = 2500.0 / ( (41.0/1024.0) * (774.-240.0) ) ; %%cap_fact = 5.0; % - - - %%data_set_name = "MCG--6-30-15 (obsid 4760)"; %%()=load_data("/spectra/d2/Fmeo/hak_4760/Obs_out/" + %% "obsid_4760_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/Fmeo/hak_4760/Obs_out/" + %% "obsid_4760_Sky_strk_1d_isis.dat"); % Ratio of zeroth-order exposure time to streak exposure time, % include fact that the HAK streak excludes ~240 rows: %%streak_fact = 2500.0 / ( (41.0/1024.0) * (774.-240.0) ) ; %%cap_fact = 2.5; % - - - %%data_set_name = "4U 1624 (obsid 4559, E<3 keV)"; %%()=load_data("/spectra/d2/Fmeo/hak_4559/Obs_out/" + %% "obsid_4559Elt3_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/Fmeo/hak_4559/Obs_out/" + %% "obsid_4559Elt3_Sky_strk_1d_isis.dat"); % Ratio of zeroth-order exposure time to streak exposure time, % include fact that the HAK streak excludes ~240 rows: %%streak_fact = 1700.0 / ( (41.0/1024.0) * (512.-240.0) ) ; %%cap_fact = 1.0; % - - - data_set_name = "3C 273 (obsid 4574)"; ()=load_data("/spectra/d2/GtoHAK/Obs_out/" + "obsid_5169_Sky_zo_1d_isis.dat"); ()=load_data("/spectra/d2/GtoHAK/Obs_out/" + "obsid_5169_Sky_strk_1d_isis.dat"); % Ratio of zeroth-order exposure time to streak exposure time, % include fact that the HAK streak excludes ~240 rows: streak_fact = 2500.0 / ( (41.0/1024.0) * (774.-240.0) ) ; cap_fact = 0.8; % - - - % % Simulated data to measure... % - - - % Nominal Capella simulation %%data_set_name = "Simulation (Capella, 500 um CD)"; %%()=load_data("/spectra/d2/GtoHAK/MARX_out/" + %% "Capella_wbkg_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/GtoHAK/MARX_out/" + %% "Capella_wbkg_Sky_strk_1d_isis.dat"); %%cap_fact = 1.0; % - - - % Narrow lines only, no background % % evt1a file from MARX-->CXC_processing %%data_set_name = "Simulation (CXC-NoRand Narrow16, E<3, 0.2mm CD)"; %%()=load_data("/spectra/d2/GtoHAK/MARX_out/" + %% "Narrow16_NoRand_Elt3_evt1a_Sky_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/GtoHAK/MARX_out/" + %% "Narrow16_NoRand_Elt3_evt1a_Sky_strk_1d_isis.dat"); %%cap_fact = 3.5; % HAK analysis using the (MARX) TDET coordinates... %%data_set_name = "Simulation (HAKaTDET Narrow16RN, E<3, 0.5mm CD)"; %%()=load_data("/spectra/d2/GtoHAK/MARX_out/" + %% "Narrow16_NoRand_Elt3_evt1a_aTDET_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/GtoHAK/MARX_out/" + %% "Narrow16_NoRand_Elt3_evt1a_aTDET_strk_1d_isis.dat"); %%cap_fact = 0.66*3.5; %****************************************** % default is for cannonical set of four HETG spectra: megm=1; megp=2; hegm=3; hegp=4; % modify colors cmegm=red; cmegp=orange; chegm=yellow; chegp=green; % Also load simulation % Simulation: aTDET w/background %%()=load_data("/spectra/d2/GtoHAK/MARX_out/" + %% "sim_capella_aTDET_zo_1d_isis.dat"); %%()=load_data("/spectra/d2/GtoHAK/MARX_out/" + %% "sim_capella_aTDET_strk_1d_isis.dat"); % Set sim_fact so that sim_capella and Capella STREAKs agree % 1/0.66 factor from time apodization in sim aTDET analysis %%sim_fact = 1.0/0.66; % Simulation: Sky, no background ()=load_data("/spectra/d2/GtoHAK/MARX_out/" + "Capella_Sky0.35_Sky_zo_1d_isis.dat"); ()=load_data("/spectra/d2/GtoHAK/MARX_out/" + "Capella_Sky0.35_Sky_strk_1d_isis.dat"); % Set sim_fact so that sim_capella and Capella STREAKs agree % Factor of 1.0 should be close ... sim_fact = 1.0; % define simulation indices and colors simm=5; simp=6; csimm=purple; csimp=mediblue; %--------------- remove the average background level -------------------- % Loop over data sets! message(""); fdata=1; while ( fdata < 7 ) { rsout = region_sum(fdata, 1.0,50.5); background = rsout.sum; rsout = region_sum(fdata, 149.5,199.0); background = rsout.sum + background; % get the data set gds = get_data_counts(fdata); % the sum was over 1/2 of the bins background = background/(0.5*length(gds.value)); message("hist index "+string(fdata)+" background cnts/bin = "+ string(background)); gds.value = gds.value - background; put_data_counts(fdata, gds); fdata = fdata + 1; } message(""); %--------------- set the data colors -------------------- set_data_color(megm,cmegm); set_data_color(megp,cmegp); set_data_color(hegm,chegm); set_data_color(hegp,chegp); set_data_color(simm,csimm); set_data_color(simp,csimp); % % Process the 1D profiles (in pixels) to be offset to % realistic Angstroms... % use_angpix = meg_angpix; % Capella % Zero-order gds = get_data_counts(megm); gds.bin_lo = 15.0+use_angpix*(gds.bin_lo-100.0); gds.bin_hi = 15.0+use_angpix*(gds.bin_hi-100.0); gds.value = cap_fact * gds.value; gds.err = cap_fact * gds.err; put_data_counts(megm, gds); % Streak, scale the counts gds = get_data_counts(megp); gds.bin_lo = 15.0+use_angpix*(gds.bin_lo-100.0); gds.bin_hi = 15.0+use_angpix*(gds.bin_hi-100.0); gds.value = cap_fact * gds.value * nom_streak_fact; gds.err = cap_fact * gds.err * nom_streak_fact; put_data_counts(megp, gds); % Capella SIMULATION % Zero-order gds = get_data_counts(simm); gds.bin_lo = 15.0+use_angpix*(gds.bin_lo-100.0); gds.bin_hi = 15.0+use_angpix*(gds.bin_hi-100.0); gds.value = cap_fact * sim_fact * gds.value; gds.err = cap_fact * sim_fact * gds.err; put_data_counts(simm, gds); % Streak, scale the counts gds = get_data_counts(simp); gds.bin_lo = 15.0+use_angpix*(gds.bin_lo-100.0); gds.bin_hi = 15.0+use_angpix*(gds.bin_hi-100.0); gds.value = cap_fact * sim_fact * gds.value * nom_streak_fact; gds.err = cap_fact * sim_fact * gds.err * nom_streak_fact; put_data_counts(simp, gds); % Obs being measured % Zero-order gds = get_data_counts(hegm); gds.bin_lo = 15.0+use_angpix*(gds.bin_lo-100.0); gds.bin_hi = 15.0+use_angpix*(gds.bin_hi-100.0); gds.value = gds.value; gds.err = gds.err; put_data_counts(hegm, gds); % Streak, scale the counts gds = get_data_counts(hegp); gds.bin_lo = 15.0+use_angpix*(gds.bin_lo-100.0); gds.bin_hi = 15.0+use_angpix*(gds.bin_hi-100.0); gds.value = gds.value * streak_fact; gds.err = gds.err * streak_fact; put_data_counts(hegp, gds); % Notice only the core of the streak histograms fdata=2; while ( fdata < 7 ) { xnotice(fdata, 15.0-2.0*use_angpix, 15.0+2.0*use_angpix); fdata = fdata + 2; } % Plot the data... %%LABEL_BY_DEFAULT=0; window(p1); plot_bin_density; ylog; yrange(10.0^3.0, 10.0^8.5); %%yrange(10.0^4.0, 10.0^9.5); xrange(15.0-10.0*use_angpix, 15.0+10.0*use_angpix); color(white); % Capella erase; plot_data(megp); oplot_data(megm); color(cmegp); xylabel(14.92, 10^6.4, "Capella Streak"); color(cmegm); xylabel(14.92, 10^6.2, "Capella Zeroth-order (500 um CD)"); color(white); xylabel(15.04, 10^6.4, "streak factor = "+string(streak_fact) ); % Observation oplot_data(hegm); color(chegm); xylabel(14.92, 10^5.8, data_set_name+" Zeroth-order"); % Observation streak oplot_data(hegp); color(chegp); xylabel(14.92, 10^6.0, data_set_name+" Streak"); % color(white); xylabel(15.04, 10^6.2, "Cap factor = "+string(cap_fact) ); data_xtitle="Pseudo Wavelength (" + string(use_angpix) + " A per pixel)"; label(data_xtitle,"log Counts/Pseudo-A", data_set_name+": Zeroth Order 1D Projections"); % simulation while (0 > 1) { erase; plot_data(simm); oplot_data(simp); color(white); xylabel(15.04, 10^6.0, "simulation factor = "+string(sim_fact) ); color(csimp); xylabel(14.985, 10^3.6, "Sim_1318 Streak"); color(csimm); xylabel(14.985, 10^3.4, "Sim_1318 Zeroth-order"); color(white); xylabel(15.04, 10^6.4, "streak factor = "+string(streak_fact) ); } % % Fit the observation's zo and streak % % Do it by manually pasting these lines... % if (2 > 1) { % Zeroth order fdata = hegm; cfdata = chegm; ()=evalfile("fit_line.i"); % Streak fdata = hegp; cfdata = chegp; ()=evalfile("fit_line.i"); } %%LABEL_BY_DEFAULT=1;