% 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 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; % % 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^2.0, 10.0^7.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, 6.4, "Capella Streak"); color(cmegm); xylabel(14.92, 6.2, "Capella Zeroth-order (500 um CD)"); color(white); xylabel(15.04, 6.4, "streak factor = "+string(streak_fact) ); % Observation oplot_data(hegm); color(chegm); xylabel(14.92, 5.8, data_set_name+" Zeroth-order"); % Observation streak oplot_data(hegp); color(chegp); xylabel(14.92, 6.0, data_set_name+" Streak"); % color(white); xylabel(15.04, 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, 6.0, "simulation factor = "+string(sim_fact) ); color(csimp); xylabel(14.985, 3.6, "Sim_1318 Streak"); color(csimm); xylabel(14.985, 3.4, "Sim_1318 Zeroth-order"); color(white); xylabel(15.04, 6.4, "streak factor = "+string(streak_fact) ); } % % Fit the observation's zo and streak % % Do it by manually pasting these lines... % while (0 > 1) { % Zeroth order fdata = hegm; cfdata = chegm; ()=evalfile("fit_line.i"); % Streak fdata = hegp; cfdata = chegp; ()=evalfile("fit_line.i"); } LABEL_BY_DEFAULT=1;