4Misc_StartN o;ROGI'XX RtEP/winspoolHP LaserJet 4V/4MV PostScripthp4mvHP LaserJet 4V/4MV PostScriptG 4dXA4PRIV'''' o;ROGI'XX RtEP/winspoolHP LaserJet 4V/4MV PostScripthp4mvHP LaserJet 4V/4MV PostScriptG 4dXA4PRIV'''' o;ROGI'XX RtEP/winspoolHP LaserJet 4V/4MV PostScripthp4mvHP LaserJet 4V/4MV PostScriptG 4dXA4PRIV'''' o;ROGI'XX RtEP/winspoolHP LaserJet 4V/4MV PostScripthp4mvHP LaserJet 4V/4MV PostScriptG 4dXA4PRIV'''' dSummary_Layout_3D0ROGI'XX-K :P$winspoolcolour laser paperNe04:colour laser paper_ 4dXA4PRIV''''J ROGI'XXp>><4!*8 WINSPOOLHP LaserJet 4MLPT1:HP LaserJet 4M@g XXLetter 8 1/2 x 11 in@MSUDNHP LaserJet 4M;d  ROGI'XXp>><4!*8 WINSPOOLHP LaserJet 4MLPT1:HP LaserJet 4M@g XXLetter 8 1/2 x 11 in@MSUDNHP LaserJet 4M;d ^Graph*@@wwwwww?wwwwww?lyWDashSettings#  ! CNormal@ Arial<HHHH$$ CNormal@ Arial<HHHH$$444444 Normal@ Arial<HHHH$$4 4 4 4 4 4 ihome5-dF\D:meis:macros:D:meis:macrosBD:\meismacros: UwD:\kwX'/w`w Nw`'|EIww5 0t /wwTw$0`otemppath;3d%C:MEIS:lateral:ErSi:C:MEIS:lateral:ErSiC:\MEIS\lateralErSi UwD:\kwX'/w`w Nw`'|EIww5 0t /ww$0` 4Misc_EndN tXOPState_Start@N Data BrowserW  t"Surface Plotterr Help Files\4XOPState_EndN L q@k3V_qV_RabV_PrV_sigaV_sigbspike_flagsmooth_flag?V_minChunkLocV_maxChunkLocV_minLayerLocV_maxLayerLocV_maxRowLocV_minRowLocV_maxColLocV_minColLocskew_flag?calib_flagmulti_flag?print_flagS_filename/wPwersiadj100_111.0009tempathme/wPwrawdatame/wPw no_filenameexpdatame/wPw er100mat.datsimdatame/wPwersiadj100_111.0009simfolder/wPwC:MEIS:lateral:ErSi:str1lder/wPwS_waveNames/wPwsim0;sim1;sim2;sim3;sim4;sim5;sim6;sim7;sim8;sim9;sim10;sim11;sim12;sim13;sim14;sim15;sim16;sim17;sim18;sim19;sim20;sim21;sim22;sim23;sim24;sim25;sim26;sim27;sim28;sim29;sim30;sim31;sim32;sim33;sim34;sim35;sim36;sim37;sim38;sim39;sim40;sim41;sim42;sim43;sim44;sim45;sim46;sim47;sim48;sim49;sim50;sim51;sim52;sim53;sim54;sim55;sim56;sim57;sim58;sim59;sim60;sim61;sim62;sim63;sim64;sim65;sim66;sim67;sim68;sim69;sim70;sim71;sim72;sim73;sim74;sim75;sim76;sim77;sim78;sim79;sim80;sim81;sim82;sim83;sim84;sim85;sim86;sim87;sim88;sim89;sim90;S_pathames/wPwC:MEIS:lateral:ErSi:Msimulation_summary() KillWaves/A/Z/F path: "D:user experiments:york:ho-er-si:" the rvalue is 500.104 *********************** the rvalue is 500.104 *********************** *********************** the rvalue is 390.533 *********************** *********************** the rvalue is 255.719 *********************** *********************** the rvalue is 64.3405 *********************** *********************** the rvalue is 28.0279 *********************** *********************** the rvalue is 28.0279 *********************** % of spikes detected per file= Inf *********************** the rvalue is 28.0279 *********************** % of spikes detected per file= 28.8288 *********************** the rvalue is 18.9033 *********************** % of spikes detected per file= 34.2342 *********************** the rvalue is 18.9033 *********************** % of spikes detected per file= 39.6396 *********************** the rvalue is 18.9033 *********************** % of spikes detected per file= 45.045 ShowTools/A ShowTools/A HideTools/A HideTools/A DoWindow/C/R comparison_single ShowTools/A HideTools/A ShowTools/A comparison() ShowTools/A comparison_single() ShowTools/A Button start proc=single_start HideTools/A *********************** the rvalue is 2.89692 *********************** % of spikes detected per file= 0 *********************** the rvalue is 0.983588 *********************** % of spikes detected per file= 5.94059 *********************** the rvalue is 0.870068 *********************** % of spikes detected per file= 0 *********************** the rvalue is 0.825048 *********************** % of spikes detected per file= 0 *********************** the rvalue is 0.801602 *********************** % of spikes detected per file= 0 *********************** the rvalue is 0.788644 *********************** % of spikes detected per file= 0 *********************** the rvalue is 0.78198 *********************** % of spikes detected per file= 0 *********************** the rvalue is 0.779465 *********************** % of spikes detected per file= 0 *********************** the rvalue is 0.779783 *********************** % of spikes detected per file= 0 *********************** the rvalue is 0.782048 *********************** % of spikes detected per file= 0 *********************** the rvalue is 0.785643 *********************** % of spikes detected per file= 0 *********************** the rvalue is 0.790125 *********************** % of spikes detected per file= 0 *********************** the rvalue is 0.795183 *********************** % of spikes detected per file= 0 *********************** the rvalue is 0.800596 *********************** % of spikes detected per file= 0 *********************** the rvalue is 0.806211 *********************** % of spikes detected per file= 0 *********************** the rvalue is 0.811923 *********************** % of spikes detected per file= 0 *********************** the rvalue is 0.81766 *********************** % of spikes detected per file= 0 *********************** the rvalue is 0.823378 *********************** % of spikes detected per file= 0 *********************** the rvalue is 0.829049 *********************** % of spikes detected per file= 0 *********************** the rvalue is 0.834658 *********************** % of spikes detected per file= 0 *********************** the rvalue is 0.840198 *********************** % of spikes detected per file= 0 comparison_single() *********************** the rvalue is 1.22504 *********************** % of spikes detected per file= 5.40541 *********************** the rvalue is 1.11138 *********************** % of spikes detected per file= 0 *********************** the rvalue is 1.05764 *********************** % of spikes detected per file= 0 *********************** the rvalue is 1.02341 *********************** % of spikes detected per file= 0 *********************** the rvalue is 0.998837 *********************** % of spikes detected per file= 0 *********************** the rvalue is 0.980372 *********************** % of spikes detected per file= 0 *********************** the rvalue is 0.966279 *********************** % of spikes detected per file= 0 *********************** the rvalue is 0.955487 *********************** % of spikes detected per file= 0 *********************** the rvalue is 0.947239 *********************** % of spikes detected per file= 0 *********************** the rvalue is 0.940966 *********************** % of spikes detected per file= 0 *********************** the rvalue is 0.93624 *********************** % of spikes detected per file= 0 *********************** the rvalue is 0.93273 *********************** % of spikes detected per file= 0 *********************** the rvalue is 0.930188 *********************** % of spikes detected per file= 0 *********************** the rvalue is 0.928424 *********************** % of spikes detected per file= 0 *********************** the rvalue is 0.927295 *********************** % of spikes detected per file= 0 *********************** the rvalue is 0.926691 *********************** % of spikes detected per file= 0 *********************** the rvalue is 0.926532 *********************** % of spikes detected per file= 0 *********************** the rvalue is 0.926754 *********************** % of spikes detected per file= 0 *********************** the rvalue is 0.927309 *********************** % of spikes detected per file= 0 *********************** the rvalue is 0.928161 *********************** % of spikes detected per file= 0 *********************** the rvalue is 0.92928 *********************** % of spikes detected per file= 0 *********************** the rvalue is 1.22504 *********************** % of spikes detected per file= 5.40541 *********************** the rvalue is 0.954985 *********************** % of spikes detected per file= 5.82524 simulation_summary() ShowTools/A SetVariable setvar8 limits={1,Inf,1} HideTools/A simulation_summary() *********************** the rvalue is 1.41582 *********************** % of spikes detected per file= 5.40541 variable/g spike_flag variable/g smooth_flag variable/g skew_flag variable/g calib_flag path: "I:meis:dysi:dysisimuls:dy100t48:" variable/g multi_flag variable/g print_flag Started on Tue, May 08, 2001 at 4:47:42 PM 121 files to process Finished on Tue, May 08, 2001 at 4:48:19 PM path: "I:meis:dysi:dysisimuls:dy111t48:" Started on Tue, May 08, 2001 at 4:49:52 PM 121 files to process Finished on Tue, May 08, 2001 at 4:50:29 PM Finished on Tue, May 08, 2001 at 4:50:54 PM batch_end() Finished on Tue, May 08, 2001 at 5:09:41 PM path: "H:user_experiments:steve tear:Er_Si:ErSi1:" Started on Mon, May 21, 2001 at 8:52:35 AM 121 files to process tidynow() Started on Mon, May 21, 2001 at 8:55:03 AM 121 files to process Finished on Mon, May 21, 2001 at 8:55:35 AM Started on Mon, May 21, 2001 at 9:01:23 AM 121 files to process Finished on Mon, May 21, 2001 at 9:01:55 AM path: "H:user_experiments:steve tear:Er_Si:ErSi2:" General text load from "er111mat.dat" LoadWave is unable to find column names on line 0 Data length: 116, waves: wave0, wave1 Edit wave0,wave1 Edit wave0,wave1 General text load from "er111mat.dat" LoadWave is unable to find column names on line 0 Data length: 116, waves: wave2, wave3 Edit wave2,wave3 Save/G/M="\r\n" wave2,wave3 as "wave2++.dat" Started on Mon, May 21, 2001 at 9:06:15 AM 121 files to process Finished on Mon, May 21, 2001 at 9:06:47 AM General text load from "er111data.dat" LoadWave is unable to find column names on line 0 Data length: 115, waves: wave0, wave1 Edit wave0,wave1 Save/G/M="\r\n" wave0,wave1 as "wave0++.dat" Started on Mon, May 21, 2001 at 9:07:56 AM 121 files to process Finished on Mon, May 21, 2001 at 9:08:29 AM Finished on Mon, May 21, 2001 at 9:11:01 AM batch_end() Finished on Mon, May 21, 2001 at 9:11:06 AM Printout() Summary_layout_3D() Printout() General text load from "rgd111_1cut.dat" LoadWave is unable to find column names on line 0 General text load from "gdsi111_1cut.dat" Data length: 136, waves: angle, data Display data vs angle Edit angle,data Save/G/M="\r\n" angle,data as "angle++.dat" tidynow() batch_end() Finished on Mon, May 21, 2001 at 2:26:29 PM batch_end() Finished on Mon, May 21, 2001 at 2:26:55 PM path: "H:user_experiments:steve tear:GdSi:GdSi2:" General text load from "111_data.dat" LoadWave is unable to find column names on line 0 Data length: 127, waves: wave0, wave1 General text load from "GdSi2.0001" LoadWave is unable to find column names on line 4 Data length: 126, waves: wave2, wave3, wave4, wave5 LoadWave is unable to find column names on line 135 Edit wave0,wave1,wave2,wave3,wave4,wave5 Edit wave0,wave1 Save/G/M="\r\n" wave0,wave1 as "wave0++.dat" Started on Mon, May 21, 2001 at 2:30:10 PM 121 files to process Finished on Mon, May 21, 2001 at 2:30:30 PM path: "H:user_experiments:Paul Dastoor:investigative simulations:PtCu_05:" Started on Thu, Oct 25, 2001 at 11:07:27 AM 121 files to process Finished on Thu, Oct 25, 2001 at 11:07:50 AM Started on Thu, Oct 25, 2001 at 11:07:59 AM 121 files to process Finished on Thu, Oct 25, 2001 at 11:08:23 AM General text load from "PtCu_06.out" LoadWave is unable to find column names on line 2 LoadWave is unable to find column names on line 208 LoadWave is unable to find column names on line 414 Data length: 200, waves: wave0, wave1, wave2, wave3 LoadWave is unable to find column names on line 650 Display wave1,wave2,wave3 vs wave0 Save/G/M="\r\n" wave0,wave2 as "wave0++.dat" RemoveFromGraph wave1,wave3 wave2+=gnoise(1) wave2-=gnoise(1) General text load from "data.dat" LoadWave is unable to find column names on line 0 KillWaves/A/Z/F General text load from "data.dat" LoadWave is unable to find column names on line 0 Data length: 200, waves: wave1, wave3 KillWaves/A/Z/F General text load from "data.dat" LoadWave is unable to find column names on line 0 Data length: 200, waves: wave0, wave1 Display wave1 vs wave0 wave1+=gnoise(0.01) wave1+=gnoise(0.01) Save/G/M="\r\n" wave0,wave1 as "wave0++.dat" Started on Thu, Oct 25, 2001 at 11:17:05 AM 121 files to process Finished on Thu, Oct 25, 2001 at 11:17:28 AM Calc_RValue_Xsq() rvalue=0 Calc_RValue_Xsq() print ne1 2 print ne2 197 print exp_avg 0.739824 print sim_avg 0.73707 General text load from "data.dat" LoadWave is unable to find column names on line 0 Data length: 200, waves: wave0, wave1 Display wave1 vs wave0 KillWaves/A/Z General text load from "data+noise.dat" LoadWave is unable to find column names on line 0 Data length: 200, waves: wave0, wave1 Display wave1 vs wave0 wave1*=100 Save/G/M="\r\n" wave0,wave1 as "wave0++.dat" tidynow() Started on Thu, Oct 25, 2001 at 11:26:26 AM 121 files to process Finished on Thu, Oct 25, 2001 at 11:26:48 AM rvalue=0 Calc_RValue_Xsq() General text load from "data.dat" LoadWave is unable to find column names on line 0 Data length: 200, waves: wave0, wave1 Display wave1 vs wave0 wave1*=100 wave1+=gnoise(10) General text load from "data.dat" LoadWave is unable to find column names on line 0 Data length: 200, waves: wave2, wave3 wave1=wave3 wave1*=1000 wave1+=gnoise(31) wave1=wave3 wave1*=10000 wave1+=gnoise(100) wave1=wave3 wave1*=3000 wave1+=gnoise(55) wave1=wave3 wave1*=6000 wave1+=gnoise(77) wave1=wave3 wave1*=4000 wave1+=gnoise(63) Save/G/M="\r\n" wave0,wave1 as "wave0++.dat" tidynow() Started on Thu, Oct 25, 2001 at 11:32:22 AM 121 files to process Finished on Thu, Oct 25, 2001 at 11:32:45 AM printout() General text load from "ptcu_07.0061" LoadWave is unable to find column names on line 4 LoadWave is unable to find column names on line 209 LoadWave is unable to find column names on line 414 LoadWave is unable to find column names on line 657 General text load from "ptcu_07.0061" LoadWave is unable to find column names on line 4 Data length: 200, waves: wave0, wave1, wave2, wave3 LoadWave is unable to find column names on line 209 Display wave1 vs wave0 Display wave2 vs wave0 wave2*=4000 wave2+=gnoise(64) KillWaves/A/Z General text load from "ptcu_07.0061" LoadWave is unable to find column names on line 4 Data length: 200, waves: wave0, wave1, wave2, wave3 LoadWave is unable to find column names on line 209 Display wave2 vs wave0 wave1=wave2 wave2*=4000 WaveStats wave2 V_npnts= 200; V_numNaNs= 0; V_numINFs= 0; V_avg= 2808.06; V_sdev= 380.359; V_rms= 2833.58; V_adev= 338.527; V_skew= -0.622643; V_kurt= -1.11398; V_minloc= 52; V_min= 2064; V_maxloc= 35; V_max= 3392; wave2=wave1 wave2*=5714 wave2+=gnoise(64) Save/G/M="\r\n" wave0,wave2 as "wave0++.dat" WaveStats wave2 V_npnts= 200; V_numNaNs= 0; V_numINFs= 0; V_avg= 4004.35; V_sdev= 549.379; V_rms= 4041.67; V_adev= 487.797; V_skew= -0.597463; V_kurt= -1.13065; V_minloc= 52; V_min= 2884.28; V_maxloc= 34; V_max= 4833.85; path: "H:user_experiments:Paul Dastoor:investigative simulations:PtCu_07:" Started on Thu, Oct 25, 2001 at 2:59:40 PM 121 files to process Finished on Thu, Oct 25, 2001 at 3:00:03 PM General text load from "ptcu_07.0061" LoadWave is unable to find column names on line 4 Data length: 200, waves: wave0, wave1, wave2, wave3 LoadWave is unable to find column names on line 209 Display wave2 vs wave0 wave2/=(sin(wave0/57.3/2)^4) WaveStats wave2 V_npnts= 200; V_numNaNs= 0; V_numINFs= 0; V_avg= 4.80413; V_sdev= 2.29475; V_rms= 5.32158; V_adev= 1.8497; V_skew= 0.749233; V_kurt= -0.447258; V_minloc= 152; V_min= 2.16221; V_maxloc= 0; V_max= 10.4988; wave2*=5714 wave2+=gnoise(wave2^0.5) wave2*=(sin(wave0/57.3/2)^4) wave2/=(sin(wave0/57.3/2)^4) wave2/=2 wave2+=gnoise(wave2^0.5) wave2*=(sin(wave0/57.3/2)^4) Save/G/M="\r\n" wave0,wave2 as "wave0++.dat" Started on Thu, Oct 25, 2001 at 3:39:56 PM 121 files to process Finished on Thu, Oct 25, 2001 at 3:40:22 PM printout() path: "E:inas:1dta:" the rvalue is 2.07297 Started on 03 Jun 2003 at 15:39:39 121 files to process Finished on 03 Jun 2003 at 15:40:17 % of spikes detected per file= 0 path: "E:inas:inas1dta:" the rvalue is 10.1038 Started on 03 Jun 2003 at 15:40:56 625 files to process the rvalue is 10.1038 Started on 03 Jun 2003 at 15:42:41 625 files to process Finished on 03 Jun 2003 at 15:46:13 % of spikes detected per file= 1.70182 path: "E:inas:out1:" path: "E:inas:1dta:" the rvalue is 4.88416 Finished on 03 Jun 2003 at 15:48:04 % of spikes detected per file= 1.6733 the rvalue is 4.88416 Started on 03 Jun 2003 at 15:48:42 121 files to process Finished on 03 Jun 2003 at 15:48:56 % of spikes detected per file= 0.215221 the rvalue is 3.79363 Started on 03 Jun 2003 at 15:49:42 121 files to process Finished on 03 Jun 2003 at 15:49:56 % of spikes detected per file= 0.398388 Started on 03 Jun 2003 at 15:50:34 121 files to process Finished on 03 Jun 2003 at 15:50:47 % of spikes detected per file= 0.398388 Started on 03 Jun 2003 at 15:51:55 121 files to process Finished on 03 Jun 2003 at 15:52:09 % of spikes detected per file= 0.398388 path: "E:inas:1dtb:" the rvalue is 1.88594 Started on 03 Jun 2003 at 15:54:18 0 files to process path: "E:inas:1dta:" the rvalue is 1.88594 Started on 03 Jun 2003 at 15:54:49 121 files to process Finished on 03 Jun 2003 at 15:55:03 % of spikes detected per file= 0.398388 the rvalue is 2.48342 Started on 03 Jun 2003 at 15:55:29 121 files to process Finished on 03 Jun 2003 at 15:55:42 % of spikes detected per file= 0.398388 the rvalue is 5.13149 Started on 03 Jun 2003 at 15:56:11 121 files to process Finished on 03 Jun 2003 at 15:56:24 % of spikes detected per file= 0.398388 the rvalue is 1.44102 Started on 03 Jun 2003 at 15:56:42 121 files to process Finished on 03 Jun 2003 at 15:56:52 the rvalue is 1.50406 Started on 03 Jun 2003 at 15:57:41 121 files to process Finished on 03 Jun 2003 at 15:57:54 % of spikes detected per file= 0.398388 the rvalue is 3.92356 Started on 04 Jun 2003 at 09:44:48 121 files to process Finished on 04 Jun 2003 at 09:45:03 % of spikes detected per file= 0.398388 the rvalue is 1.88594 Started on 04 Jun 2003 at 09:49:01 121 files to process Finished on 04 Jun 2003 at 09:49:15 % of spikes detected per file= 1.13106 Started on 04 Jun 2003 at 09:49:33 121 files to process Finished on 04 Jun 2003 at 09:49:46 % of spikes detected per file= 1.20511 the rvalue is 2.26004 the rvalue is 2.26004 Started on 04 Jun 2003 at 09:51:05 121 files to process Finished on 04 Jun 2003 at 09:51:18 % of spikes detected per file= 0.398388 the rvalue is 4.54597 Started on 04 Jun 2003 at 09:51:39 121 files to process Finished on 04 Jun 2003 at 09:51:53 % of spikes detected per file= 0.398388 the rvalue is 4.54597 Started on 04 Jun 2003 at 09:52:40 121 files to process Finished on 04 Jun 2003 at 09:52:54 % of spikes detected per file= 0.398388 the rvalue is 2.48342 Started on 04 Jun 2003 at 09:53:28 121 files to process Finished on 04 Jun 2003 at 09:53:41 % of spikes detected per file= 0.398388 the rvalue is 2.48342 Started on 04 Jun 2003 at 09:54:08 121 files to process Finished on 04 Jun 2003 at 09:54:21 % of spikes detected per file= 0.398388 the rvalue is 5.13149 Started on 04 Jun 2003 at 09:55:30 121 files to process Finished on 04 Jun 2003 at 09:55:44 % of spikes detected per file= 0.398388 the rvalue is 2.48342 the rvalue is 2.48342 Started on 04 Jun 2003 at 09:56:14 121 files to process Finished on 04 Jun 2003 at 09:56:27 % of spikes detected per file= 0.398388 the rvalue is 2.48342 Started on 04 Jun 2003 at 10:33:18 121 files to process Finished on 04 Jun 2003 at 10:33:31 % of spikes detected per file= 0.398388 path: "D:meis:map13-08:test:" path: "C:MEIS:lateral:HoSi:hosi100_111adj:" filename= hosi100_111adj.0001 Started on 16/09/03 at 12:05:23 45 files to process Finished on 16/09/03 at 12:05:29 path: "D:meis:dysi:Lateral:leedmulti2:" filename= leedmulti2.0001 Started on 16/09/03 at 12:06:21 10 files to process Finished on 16/09/03 at 12:06:22 path: "C:MEIS:lateral:ErSi:" filename= ersiadj100_111.0001 Started on 16/09/03 at 12:07:00 9 files to process Finished on 16/09/03 at 12:07:01 .sim0t?eD@D@@33D@LD@`ffD@D@D@@3D@D@`fD@E@E@@33E@LE@`ffE@E@E@@3E@E@`fE@F@F@@33F@LF@`ffF@F@F@@3F@F@`fF@G@G@@33G@LG@`ffG@G@G@@3G@G@`fG@H@H@@33H@LH@`ffH@H@H@@3H@H@`fH@I@I@@33I@LI@`ffI@I@I@@3I@I@`fI@J@J@@33J@LJ@`ffJ@J@J@@3J@J@`fJ@K@K@@33K@LK@`ffK@K@K@@3K@K@`fK@L@L@@33L@LL@`ffL@L@L@@3L@L@`fL@M@M@@33M@LM@`ffM@M@M@@3M@M@`fM@N@N@@33N@LN@`ffN@N@N@@3N@N@`fN@O@O@@33O@LO@`ffO@O@.yhsim1t?e-?`"Y?(\µ?V-?-'?ʡE?S㥛? rh?S?Cl?/$?h|?U?-6?NbX9?x&1?A`"+?&1?1Z?Q?v?/d?*?MbX?@5^I?'1z?l?J +?}?5^ ?Gz?.?ʡE?ףp= g?x?"~"?fffff?Ed?v?v/?}?5^i?zG1??X9v?Gz?(\e?> ףp=?Q?Q?(\o?Dl?Q;? rha?> ףp]?p= 3?Dl?-? ףp?zG?v?ClI?+ٮ?/$?(1?V-=?'1?-w?o?jt?(\/?x?Cl7?x&1?x ?$c?d;O?tV?Gz4?Zd;O]?|?A`В?V?ƛ?A`В?X9Ȇ?y&1|?d;Ow?.sim2t?ex@`"@X9v@~jt@i|?5@!rh@> ףp=@ r@d;O@+@ˡE@Mb@"~@n@@"~j@Zd;@kt@5^I @S㥛@ʡE@㥛 @L7A`@"~@ rh@ʡE@"~j @jt @K @S@Zd;O@1Zd@^I +@/$@"~@Mb@K7A@ˡE@(\@n@M@5^I @"~@|?5^@|?5^@/$ @i|?5 @RQ @ rh@i|?5@ʡE@@(\@E @/$ @Pn @Mb@ r@x@T㥛 @-@X9v@Q@7A`@MbX@X9@ʡE6@L7A`e@v/@S㥛 @M @/$ @S@J +@Cl@/$@Zd;@l@K@Dl @Q @K7A @I + @Zd;O @㥛 @ˡE @tV @B`" @#~j @!rh @ rh @Zd;O @|?5^ @ @RQ @S㥛 @/$ @Cl @Zd; @Cl @- @"~ @"~ @n @(\ @333333 @Cl @Zd; @Zd; @- @(\ @I +@Q@i|?5@V-@Pn@.asim_smootht?``J +@Cl@m@V-@On@@'1Z@(\@~jt@RQ@J +@S@ʡE @ r @p= ף @\(\ @ rh @5^I  @Cl @q= ףp @ rh @NbX9 @'1Z@On@On@+@V-@#~j<@5^I @S㥛 @@S㥛@1Zd@I + @V- @(\ @/?@Nɳ?hA?29?p? ?rc;?+5W?ܧ:q?Q~?Dz?>({??" =?Q ?fG+$?;?4zBS?j?v+?f?nz}.?\?lhE?vB?2f? ?D.+?r z>?X}Q?gd?D %v?0ш?sx?+@@?&?t 3?/]?Fh *?Q0? f ?PXT?`r(?]6?$(D?e\Q?:1^?tpj?pbev?*8u?JG΍?XȘ?\* c?j0&?ΞP|?!t?۰?ݟ-?";? >?\?z X&?  ?$͒?+4?x ?4? 录?m7?h]?`("?p|o&?kfݟ)?0[Tq,?p;.?`=|v0?!2?4?;4?5?.exp1t?`eѦW1!?Yk?*[ɛ?\[?i fm??ff?u*x?HO ?Hm?G6?kR?!?\?*3?7>K?f D?EP Vw?n=ǵ?y9G?0??flb??9/?[Z?-?f3?$<[?+֬r+V?W=?ms?r_?̎(??}?E|??D?{2o\?þy?%Y??JV?I|5??/~?:@?jZ[Q? rc?[(}`?KrF?i}?~f?<$D?F)?X? vq? ͌׎?'Z?1/?]/?R?Gx?Рjy?eA;?T ?;vo?J8?o(s?r?q*$?c򩩵? p>d#?5N?Tb?qIF?o[+O?b?ˠa?A4?u4gB?n*?Y?Uό%?<4P?xm?~f?Y]n?zOͬ?I~?ֻ?jMZ?ᅶ\?U<7?!Qk@ l@sl@im@In@m@Gl@ M,m@m@@,m@Kl@ =l@@j@ |i@ 0i@mi@,h@Jf@@[e@lf@5e@`fe@re@`e@ 7f@e@`e@`&f@g@@Bg@`g@`~h@6i@֯i@j@Kk@ k@ Fk@`k@>l@fm@0Wl@k@`_l@ m@ Il@k@Zl@.aexp_fitt?`eBFVuWX?RJ?!k?]%Ia?ՁV?q,4L?dǝ+C?YL9?/?P {{&?>p^?As?q- ?}4?e? )??כ}?]Ԕ?_.?_?͒8? &6KR?D[¸?48e?$9:?I6A?Jz?6?@?FrAT?qV?S[ \?FI?ez}?Ͽ[Vy?Tu?;q?m?4V{j?Bg?{;d?Y_fa?@N^? S\?|;Z?i X?|s0V?oT?DS?2Q?9O yP?4HUO?du7O?[׽N?x{vN?baN?s$~N?N?OOO?V/P?]:VP??R?P8LS?d-T?X4/yV?>@ZX? nnZ?U#)s\?E,_?S4Ca?T_ d?j8g?Ukk?Gwn?3Mr?lu?Cz?<~?"x?`L? &B ?[%?]?vQZǛ?R#nd?3?^ 4?$/Ag?j̹?Vrud?j.?+*?#-Y?$*?\L?G9?\YO ?L03?oRValue ??de%N@)(@F9F@}+|@aQH@ʂlP@ x*@iTx"@8p3(@rbfRvalue0y?%. rBBRB}7B#BcB,;B|B:YBBBbBwB?B*BBBQBPǪByBB.BBB!BDBVBwhBŨB0BBBB@B[BJqBMBB9BbBاBHըB KBTBFB*B~}BB^B;B BABB#!BRBBBNѹB=B%BtBB>B?hBrBo1BhB:eB\hB_B7BWBwԭB7B{BB ѱB'GBChBOɿBݳB&BB5IBb*BBBWBLB`B$BBBZB?B1B8B8BhBrBm"BOBNByBֿB[BBfBBBBFCBqBzBBcBQ3B]BB*B BOBVFRvalue_origy??$.$.k,@@b-@Xt ./@K10@=1@o`2@=ȭ3@ڋf$5@m6@7@OE~%9@R4J)@YH+@b,@b.@ѿ30@-a1@0R?2@PhhU3@d%y4@u牿5@[ш6@K\'@Ky))@ p*@V+@}J-@B t/@R0@p1@92@U3@v*5@~e`k&@ov'@и*(@b*@ s+@~@;-@:ݡ.@D0_0@#Zmj1@5Szz2@u.3@3_M%@BZ&@)ԑ'@(@I*@(|zw+@d-@]O.@3F0@;1@S6*2@QR2%@T#.&@^h'@93n(@4)@һA+o*@}չ+@)|;P-@.,.@ aq`0@O 8<.1@“}%@1~4&@. ~6 '@_S'@K((@$Jڞ*@M0~>+@/͸:,@m ~.@o^/@0@D:m&@6L|"'@Ʈ`|j'@txl(@J:)@6G?)*@MuB+@P5}yg,@ |ѓ-@̲.@֒10@+E4(@jx(@9Z(@a*)@$͐Օ9*@n+@@dw@@@[@z@j@J@@]@D@@4@H@C@+@UX@ @{@{s@)@^@F@T@ @̾@*@J@ @ͽ@@V@)@@@@x@D@}@n@α@@V@+@@ٺ@@i@i@?G@m&@ @Y@͹@f@=@@k@W@,D@2@#@@X@a@@F@&@@޸@ݸ@ݸ@ @@'@5@@U@k@2@,@=@VFfit_P2_RValue29/? AW A Ah A+C AyA}A`AA3JAAyAdAaAAAOAuANAMALAAA2AAA*AAA4A@5@@ @2@9@n@@c@@D^@ @zY@@P@[@>@x@q@@@^Nfit_P3_RValuek?EsEs`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@`@_@_@_@_@_@_@_@_@_@_@_@_@_@_@_@_@_@_@_@_@_@_@_@_@_@_@_@_@_@_@_@_@_@_@_@_@_@_@_@_@_@_@_@_@_@_@_@^Nnexpt_bin1t?PVBBjBΧBBB:BBHҴB./BtB- B:ggC fCGkCuC~C{CqC|^CMC۸C_KCTCYCW`C[dCLnCOrChemC=eCibiC6oCbYlC\eCaCRUCLCɄICEnKCfAC>P2CR/Cd3C.CS7+C+CT/Cѽ1C>,Cs/C35Cm?C8C#>CKGCLC~MCt6VCXz\C_P^C!6ZC^CgC5kCbC^C+bCMhCHbC_CrgC4`* -matrix12 ????| AYAoUA@@Ԁ A3 AcA@T@ AyA"Au@ԝ@ A-AlAA#d@r_ AA0AA@8 AS* A?AAzAJ A AA\AAtA[. AU ALAA޻A A; A AK A.:sim3t?q{eV-?Q?V-?"~?(\?MbX9?M?MbX9?Fx?S㥛?oʡ?y&1?Zd;O?\(\?Cl?ʡE?"~?E?V-?Cl?i|?5?p= ף?sh|??㥛 ?B`"?V-?jt?V-?/$?Q?zG? +?ʡE?I +?333333?Fx?~jt??\(\?/$?7A`?d;O?"~j?M?T㥛 ?V-?ʡE?S?x?/$?Zd;O?}?5^I?^I +?n?{Gz? +?/$?x&1?Cl?n?X9v?n?K?x&1?HzG?d;O?V-?uV?&1? +?-??+?V-?1Zd?l?M?V-?ˡE?RQ?X9v?S?x?㥛 ?"~?i|?5?(\?㥛 ?l?$C?-?K7A?#~j?K7A`?@5^I ?e;O?K?-?'1Z?NbX9?K7A`?(\?kt?M?n?x?^I +?333333?Cl?Cl?rh|?Cl?\(\?K?$C?A`"?^Nexpt_bin0t?emP B B!Bff"B33#B$B$B%Bff&B33'B(B(B)Bff*B33+B,B,B-Bff.B33/B0B0B1Bff2B333B4B4B5Bff6B337B8B8B9Bff:B33;BB33?B@B@BABffBB33CBDBDBEBffFB33GBHBHBIBffJB33KBLBLBMBffNB33OBPBPBQBffRB33SBTBTBUBffVB33WBXBXBYBffZB33[B\B\B]Bff^B33_B`B`BaBffbB33cBdBdBeBfffB33gBhBhBiBffjB33kBlBlBmBffnB33oBpBpBqBffrB33sBtBtBuBffvB33wBxBxByBffzB33{B|B.c~sim4t?deD@D@333333D@LD@gfffffD@D@D@33333D@D@gffffD@E@E@333333E@LE@gfffffE@E@E@33333E@E@gffffE@F@F@333333F@LF@gfffffF@F@F@33333F@F@gffffF@G@G@333333G@LG@gfffffG@G@G@33333G@G@gffffG@H@H@333333H@LH@gfffffH@H@H@33333H@H@gffffH@I@I@333333I@LI@gfffffI@I@I@33333I@I@gffffI@J@J@333333J@LJ@gfffffJ@J@J@33333J@J@gffffJ@K@K@333333K@LK@gfffffK@K@K@33333K@K@gffffK@L@L@333333L@LL@gfffffL@L@L@33333L@L@gffffL@M@M@333333M@LM@gfffffM@M@M@33333M@M@gffffM@N@N@333333N@LN@gfffffN@N@N@33333N@N@gffffN@O@O@333333O@LO@gfffffO@O@.sim5t?deCl@K7@"~j@S@S@~jt@I +@/$@(\@)\( @V-@+@q= ףp @ @- @Cl @B`" @MbX9 @ ףp=@9v@n@p= ף@On@ףp= @> ףp= @~jt @|?5^ @K7@X9v@I +@ ףp= @T㥛 @ˡE@-@ʡE@`"@I +@{Gz@n@n@@Zd;@-F@Fx@Qk@T㥛 @Cl@C@)\(@Gz@(\@+w@"~@Dl)@m@X9v@1Zd@kt@w/@K7A@"~@x&@ ףp=@Q@/$@Cl@.6sim6t?dep= ף?%C?"~?{Gz?RQ?'1Z?ʡE?> ףp=?x&1?i|?5?V-?MbX9?/$?T㥛 ?)\(?T㥛 ?/$?&1?T㥛 ?jt?Mb?@5^I ?x&1?kt????+??V-?+????+??+???+?+?+??+???+?+?+?+??+??+?+?+?+??+??+?+?+?+????+?+?+???+??+??+??+?+???+???+?+??+??+?+?+??+?+???+?+???+??+?+?+?+?+?+??+?+?+?+?+?.Hsim7t?de-?(\?"~?;On?|?5^?RQ?ʡE?MbX9?i|?5?x&1?/$?V-?)\(?)\(?)\(?V-?/$?&1?T㥛 ?jt?jt?Mb?@5^I ?@5^I ?kt?kt???kt?kt????V-??kt?+??+?+?V-?+??+?kt?+??+??+?+?+??+??kt?+??+???+?kt????kt??????+?+??+?+????+??+??+???+?+???+?+???+????+?+????+??+??+?+?+?+?+?+??+?."Osim8t?dem?(\? ףp=?> ףp=?> ףp=?> ףp=?MbX9?i|?5?i|?5?x&1?x&1?i|?5?)\(?x&1?i|?5?x&1?x&1?V-?x&1?x&1?V-?x&1?V-?x&1?i|?5?V-?x&1?x&1?x&1?x&1?V-?V-?V-?.sim9t?de/$?B`"?tV?T㥛 ?-?/$?K7A?-?RQ?/$?K7A`? rh?q= ףp?333333? ףp= ?S㥛?"~j?V-?+?\(\?-?M? rh?Mb??333333?X9v?@5^I ?HzG?(\?RQ?jt?Zd;?V-?ʡE?ʡE?/$?jt?RQ?NbX9?Q?n?|?5^? +?ʡE?K7A?w/?q= ףp?&1?䥛 ?On?Cl?Cl?MbX9?RQ?;On?%C?9v?9v?Q?E?S㥛??n?Q?NbX9?#~j?#~j?-?-?p= ף?~jt?S㥛?"~?Pn?"~? r? r?;On?"~j?tV?M?}?5^I?M?M?}?5^I?ʡE?MbX9?i|?5?MbX9?MbX9?MbX9?i|?5?x&1?i|?5?i|?5?V-?x&1?x&1?i|?5?i|?5?V-?i|?5?i|?5?MbX9?x&1?)\(?i|?5?i|?5?x&1?x&1?x&1?V-?i|?5?x&1?x&1?sim10?de?@ sim11?de,@Q@sim12?deףp= <@Qd@sim13?de?@sim14?de??Osim15?deQ?Esim16?de @Rsim17?deMsim18?deQ?Q?hsim19?de@@sim20?de??Nsim21?deQ?Asim22?de#~j@^I +@psim23?deq= ףp?q= ףp?{8sim24?de??Xsim25?de@@sim26?de@@Nsim27?deQ?lsim28?deZd@|?5^?sim29?deQ@Q@ sim30?de??Lsim31?de@ @sim32?de??Msim33?deQ?Csim34?de @Asim35?de)\(@)\(@%sim36?de(\µ?(\µ?<sim37?de"@$@sim38?de??Msim39?deQ?tsim40?deZd@|?5^?$Fsim41?deV-@V-@)sim42?de(\µ?(\µ?$sim43?de*@&@(@,@sim44?de????rOsim45?deQ?Q?Qsim46?de#~j@Zd@|?5^?^I +@sim47?de^I k!@33333@33333@^I k!@sim48?de(\µ?(\µ?(\µ?(\µ?sim49?de1@.@2@0@sim50?de????vNsim51?deQ?Q?sim52?de#~j@ @^I +@ǽsim53?de㥛 '@Q&@㥛 '@Q&@sim54?de(\µ?(\µ?(\µ?(\µ? sim55?de3@5@6@4@sim56?de????pNsim57?deQ?Q?Lsim58?deZd@|?5^? @sim59?dep= c,@i|?-@i|?-@p= c,@sim60?de(\µ?(\µ?(\µ?(\µ?sim61?de9@7@8@:@sim62?de????tMsim63?deQ?Q?Ssim64?de#~j@Zd@|?5^?^I +@sim65?dew/2@GzT1@GzT1@w/2@sim66?de(\µ?(\µ?(\µ?(\µ?sim67?de=@;@>@<@sim68?de????nMsim69?deQ?Q?sim70?de#~j@ @^I +@sim71?de9v?5@ףp= w4@9v?5@ףp= w4@sim72?de(\µ?(\µ?(\µ?(\µ?Ksim73?de?@@@A@@@sim74?de????rLsim75?deQ?Q?Jsim76?deZd@|?5^? @sim77?de7@Mb8@Mb8@7@sim78?de(\µ?(\µ?(\µ?(\µ?sim79?deB@A@B@C@sim80?de????vKsim81?deQ?Q?Usim82?de#~j@Zd@|?5^?^I +@6sim83?de/݄;@\(:@\(:@/݄;@sim84?de(\µ?(\µ?(\µ?(\µ?sim85?deD@C@E@D@sim86?de????pKsim87?deQ?Q?sim88?de#~j@ @^I +@Msim89?deCl>@Q=@Cl>@Q=@sim90?de(\µ?(\µ?(\µ?(\µ?.JxVeg1t?dd}?5^I?ʡ%?A`b?/?433333?&C?Zd;?Cl?(\?Cl{?:v?:v?䥛 p?n@?FԈ??S?(\?ףp= ?S ?ҽ??:v_?(\?Gz?Q~?I +'?C?S?o1?NbX94?x&Q? ףp=? r?\(4??Gz>?V-U?X9v>?~jt ?"?Q?Gz\?ʡE+?-?\(?ƣ?'1j?5^I "??䥛 ?ˡE?Zd;_?S㥛?w/?S3?X9?On? ?Zd?)\?K7!?䥛 ?/ݴ?v?QH?l?ʡE?X9f?1Z?S?X9v?*g? rhm?`"9?NbX9?gffffF? Zd?Q? ףp=?43333S?rh?:v?Q(?v/?n*?rh?*·?Onc?v_?S㥛t?(\?43333?43333C?MbX?V=??Qe?S?MbX?"~J?ףp= ?zG? +?x?M"?Mb?"~?+?(\-?+F?nZ? rhi?(\r?I +w?jtx?W_sigma?eeXz9b?Z 8q?:Id?>?9@Hk>.Veg2t?ddT㥛`?K7)? +9?x&?On?Q?-??(\µ?NbX9t?ʡE?jtD?\(\?@5^IL?.&?ObX9?(1?~jt? ?#~j?Pn3?+?ʡE?x/?x?d;O?+ί? ףp=r?K79?|?5?x&1?h|??"~?GzT??x&? rh?S%?n?h|??nʑ?rh,?ʡE? ףp=j?y&1?jt? Zd?n ?^I K?vO?-?M?ˡE?e;O?_I [?S?QK?V?(\?Hz'?Qu?q= ף?+>? ?Zd;_???$C??C,?S㥛t?6A`?p= ף@?Gz?KW?x&?V?Mb(?x&?)\X?S?Mr?V-?Pnc?V-?S+?Et?"?Zd;O? r?#~j ?Pn?V-?h|??$C?V-?.YVeg6t?ddL?(\?|GZ?"~j?w/?ʡE?vO?B`"?~jt?Mb?2Z$?Q?X9f?/$v?K?\(\?I +'?8A`?33333S?T㥛?tf?(\?/$&?/$&? ףp=?MbX?,?&1?䥛 ?8A`p?ClW?MbX9?(\?zG1?fffff??K7A?m?bX9?䥛 p?Pn3?p= ?Zd?Cl? rhu?y&1L? ?X9? ףp?kt?ˡE?(\o?E?+َ?x&a?kt?oʱ?MbXI?S??d;O?~jT?Zd;?V-O?|?5?"~j?x?ʡ?S[?V-?NbX9?x&?SE?7A`?p= ף?䥛 ?Q[?(\?ClG?rh? ףp=?S;?n?"~j ? ?Q.?MbX9?"~?ʡE? +?ʡ?"~:?J +?E?jt?m?y&1?bX9H?S㥛?x&1|? ?Q?o!??NbX9?Fh?ʡE?S㥛?/$Q?(\?~jt?V?'1?S?/$?ʡE?fffff?33333?.oVeg8t?ee-?`"Y?(\µ?V-?-'?ʡE?S㥛? rh?S?Cl?/$?h|?U?-6?NbX9?x&1?A`"+?&1?1Z?Q?v?/d?*?MbX?@5^I?'1z?l?J +?}?5^ ?Gz?.?ʡE?ףp= g?x?"~"?fffff?Ed?v?v/?}?5^i?zG1??X9v?Gz?(\e?> ףp=?Q?Q?(\o?Dl?Q;? rha?> ףp]?p= 3?Dl?-? ףp?zG?v?ClI?+ٮ?/$?(1?V-=?'1?-w?o?jt?(\/?x?Cl7?x&1?x ?$c?d;O?tV?Gz4?Zd;O]?|?A`В?V?ƛ?A`В?X9Ȇ?y&1|?d;Ow?./Veg9t?eeA`Т?^I  ?QX?-f?d;O?.&?X9v>?SC?M2?Zd;?V-/?'1Z?Gz?(\?p= ׃?a"ۉ?"?.?䥛 P?+N?S?K?ʡ%?Gz?{G?:v_?x?J +v?&1?$?F?xy?L7A`?sh|?M?V-J? rp?rhl?/$I?~jt?Gz?v?lq? +I?(1"?fffff?p= ף?Zd;??n?%C ?zG1?Zd;/?K7 ?T㥛?+^?/?U<-<5<:<matrix12=7.9.x? W$@"z@@2ՏE @`y3 @@Vmatrix12=7.9? ? ?@@ n_ @?matrix12=8.x? @@ @@ݴ4'@J@B@matrix12=8? -h\@@@c@@@?matrix12=8.1.x? d_ˏ?@jr2,@aJ[Q@gQ/ @@n @@^matrix12=8.1? hsX??@@4\lm @@@)x@ematrix12=8.2.x ? اڍ?XXig?tVh?@;$@:Xݺ@@ d @@;matrix12=8.2 ? ?@PQ@@@TJ\o@@2ء@vmatrix12=8.3.x ? @\ý @@㮔@Fqw @@@M+#??Fgu ?۟%?matrix12=8.3 ? ;&gcr@@Jٵ@@@ |wx@@@?71matrix12=8.4.x ? Xc ?j??~mO?ھQ?TeL ?6o -?@ O@@@)matrix12=8.4 ? ?Ll7?@@@@u#V@@ ʿ@P]c@آmatrix12=8.5.x ? @z=} @@Nl"@@JH/e?zǰsF ??P3C?zDW)?x_?t& .?matrix12=8.5 ? P@@C|@@h@@@x@@@?Uimatrix12=8.6.x ? @Fiq;??&gv?Ҷ.?N$?r??"i_?@ SH@@@ matrix12=8.6 ? ?@@@@.@@Nݳys@@ i@t@matrix12=8.7.x? b @@&,U??|#<3?ԏuJ_a?iJ?6matrix12=8.7? @Z"(b@@Ru@@@@,Pֶ @matrix12=8.8.x? YY?j} L??lOR??matrix12=8.8? +@@@B!S@ @matrix12=8.9.x? J?@vM$?%~matrix12=8.9? @@Am.+@$20 -matrix12=xymarkers.x-?????@@@?@@@?@@@?@@@?@@@?@@@?@@@?@@@?@@@^=matrix12=xymarkers-? ?????@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @ @ @ @ @ $| Platform=WindowsNT, IGORVersion=3.150 NewPath/Z temppath "C:MEIS:lateral:ErSi:" DefaultFont "Arial" MoveWindow/C 13.5,268.25,493.5,578.75 dummy_graph() Contour_12() P1_RValue_Graph() r_value() comparison() simulation_summary() MoveWindow/P 261,203,762,498.5 Window simulation_summary() : Panel PauseUpdate; Silent 1 // building window... NewPanel /W=(0,43.25,641.25,484.25) as "Simulation analysis" SetDrawLayer UserBack SetDrawEnv fsize= 14,fstyle= 1 DrawText 33,31,"This program produces R-values from Vegas simulations" SetDrawEnv fstyle= 1 DrawText 30,80,"1." SetDrawEnv fstyle= 1 DrawText 30,110,"2. " SetDrawEnv fstyle= 1 DrawText 30,140,"3. " SetDrawEnv fstyle= 1 DrawText 30,170,"4. " SetDrawEnv fstyle= 1 DrawText 30,200,"5. " SetDrawEnv fstyle= 1 DrawText 30,230,"6. " SetDrawEnv fstyle= 1 DrawText 30,260,"7." SetDrawEnv fstyle= 1 DrawText 30,290,"8." SetDrawEnv fstyle= 1 DrawText 30,320,"9." SetDrawEnv fstyle= 1 DrawText 25,350,"10." SetDrawEnv fstyle= 1 DrawText 25,380,"11." SetDrawEnv fstyle= 1 DrawText 25,410,"12." SetDrawEnv fstyle= 1 DrawText 315,380,"13." SetDrawEnv fstyle= 1 DrawText 316,411,"14." SetVariable setvar0,pos={100,66},size={288,18},title="Experimental data; " SetVariable setvar0,frame=0,limits={-Inf,Inf,1},value= expdata SetVariable setvar1,pos={102,95},size={474,18},title="Simulation folder; " SetVariable setvar1,frame=0,limits={-Inf,Inf,1},value= simfolder Button button0,pos={497,291},size={91,53},proc=run,title="Run" CheckBox check0l,pos={50,123},size={185,20},proc=multi,title="Multiple simulation files",value=1 CheckBox check1,pos={50,306},size={147,20},proc=printbest,title="Print best result?",value=0 CheckBox check2,pos={50,243},size={147,20},proc=smoothset,title="Smooth simulation",value=1 CheckBox check3,pos={50,273},size={147,20},proc=spikeset,title="De-spike simulation",value=1 CheckBox check4,pos={50,212},size={147,20},proc=skewset,title="Skew data",value=1 CheckBox check4_1,pos={50,152},size={135,20},proc=element_1set,title="Use first element",value=0 CheckBox check4_2,pos={185,152},size={147,20},proc=element_2set,title="Use second element",value=1 CheckBox calibflag,pos={50,182},size={122,20},proc=calibset,title="Use calibration",value=0 SetVariable setvar4,pos={50,333},size={225,18},title="Multi parameter number 1 ;" SetVariable setvar4,limits={1,100,1},value= K11 SetVariable setvar5,pos={50,363},size={225,18},title="Multi parameter number 2 ;" SetVariable setvar5,limits={1,100,1},value= K12 SetVariable setvar6,pos={50,393},size={225,18},title="Multi parameter number 3 ;" SetVariable setvar6,limits={1,100,1},value= K13 Button find_expt,pos={45,61},size={50,20},proc=find_expt,title="find" Button find_sim,pos={46,94},size={50,20},proc=find_sim,title="find" SetVariable setcalib,pos={175,182},size={215,18},proc=SetCalibration,title="Value: counts per layer ;" SetVariable setcalib,limits={-Inf,Inf,1},value= calibration SetVariable set_calib_range,pos={399,182},size={94,18},title="% range" SetVariable set_calib_range,limits={-Inf,Inf,1},value= calibration_range SetVariable setvar7,pos={199,273},size={207,18},title="set spike tolerance level" SetVariable setvar7,limits={0.01,2,0.05},value= spikelevel SetVariable setvar8,pos={199,244},size={127,18},title="set amount " SetVariable setvar8,limits={1,Inf,2},value= smooth_amount SetVariable setvar2,pos={337,394},size={238,18},title="Parameter polynomial fit order" SetVariable setvar2,limits={3,9,2},value= poly_amount SetVariable setvar3,pos={337,364},size={238,18},title="Skew range overlap " SetVariable setvar3,limits={0,25,1},value= Ne_extra EndMacro Window comparison() : Graph PauseUpdate; Silent 1 // building window... Display /W=(0,37.25,459.75,237.5) sim1 vs sim0 as "Comparison of experiment and simulation" AppendToGraph exp1 vs sim0 ModifyGraph mode(exp1)=4 ModifyGraph marker(exp1)=19 ModifyGraph lSize(sim1)=2 ModifyGraph rgb(exp1)=(0,0,0) ModifyGraph nticks(bottom)=10 Label left "scattering yield [layers]" Label bottom "scattering angle [degrees]" Legend/J/N=text0/F=0/A=MC/X=35.00/Y=-48.00 "\\JL\\Z09\\s(sim1) simulation \\s(exp1) experiment" Cursor A exp1 11;Cursor B exp1 113 Button Restart,pos={69,30},size={55,17},proc=Vegas,title="Restart" SetVariable Browse,pos={129,7},size={110,18},proc=Select_current_sim,title="Get sim." SetVariable Browse,limits={1,Inf,1},value= sim_num Button start,pos={69,51},size={55,17},proc=batch_start,title="Run" Button reset,pos={69,8},size={55,17},proc=reset,title="Reset" EndMacro Window r_value() : Graph PauseUpdate; Silent 1 // building window... Display /W=(201.75,250.25,459.75,370.25) RValue as "progress chart" ModifyGraph zero(bottom)=1 ModifyGraph nticks(bottom)=10 ModifyGraph minor(bottom)=1 Label left "R-value" Label bottom "number of files" Textbox/N=text0/F=0/A=MC/X=-6.27/Y=39.39 "min. R-value = \\{Rv}" Cursor A RValue 8 ShowInfo ControlBar 17 Button button0,pos={5,1},size={91,16},proc=save_rvalue,title="save R-value" Button button1,pos={122,1},size={90,16},proc=load_Rvalue,title="load R-value" Button button2,pos={244,1},size={52,16},proc=revert_Rvalue,title="revert" EndMacro Window P1_RValue_Graph() : Graph PauseUpdate; Silent 1 // building window... Display /W=(470.25,37.25,680.25,167.75) P1_RValue,fit_P1_RValue ModifyGraph mode(P1_RValue)=3 ModifyGraph marker(P1_RValue)=19 ModifyGraph lSize(fit_P1_RValue)=2 ModifyGraph rgb(fit_P1_RValue)=(0,0,0) Label left "Parameter one R-value" Label bottom "Parameter one" Textbox/N=text0/F=0/A=MC/X=3.00/Y=40.00 "\\JCmin. value@ \\{p1_min}\rsigma = \\{p1_sigma}" Tag/N=text1/F=0/X=-6.70/Y=50.00 fit_P1_RValue, 0, "- sigma" Tag/N=text2/F=0/X=1.30/Y=50.00 fit_P1_RValue, 8, "+ sigma" Cursor A fit_P1_RValue -1.97753;Cursor B fit_P1_RValue 16 EndMacro Window Contour_12() : Graph PauseUpdate; Silent 1 // building window... Display /W=(0,250.25,200.25,370.25) as "Contour 1 and 2" AppendMatrixContour matrix12 ModifyContour matrix12 ctabLines={*,*,YellowHot,1}, xymarkers=1, labels=0, tracesDF=root:WM_CTraces_matrix12: AppendImage matrix12 ModifyImage matrix12 ctab= {*,*,Grays,0} ModifyGraph mode('matrix12=xymarkers')=3 Label left "parameter one" Label bottom "parameter two" EndMacro Window dummy_graph() : Graph PauseUpdate; Silent 1 // building window... Display /W=(160.5,48.5,556.5,255.5) expt_bin1 as "dummy" EndMacro o |version 2 - 6/3/96 - expt and sim matched to average between 60 and 135 degrees |version 3 - 12/3/96 - low angle intensity fall-off corrected for |version 4 - 15/3/96 - new exptl. data |version 5 - 18/3/96 - overall curvature messed about with |version 6 - 19/3/96 - smoothing added to simulation |version 7 - 19/3/96 - IGOR pro |version 8 - 21/3/96 - batch processing and smooth curve correction for each file |version 9 - 8/9/96 - working towards all singing, all dancing versuion |version 10 - 12/1/97 - ??? |version 11 - 15/1/97 - Contours added |version 12 - 21/1/97 - R-value error fixed |version 59 - 28/1/98 - files loaded by name (for Dave Spence, York) |version 60 - 11/2/98 - oops! test Xsq expression was being used (only minor effect) |version 61 - 19/3/98 - Rvalue parameter fit now uses 5th order polygon |version 62 - 20/3/98 - Rvalue scaling error fixed (sim_avg) |version 63 - 20/3/98 - first attempt at including absolute calibration |version 68 - 23/3/98 - refined |version 70 - 25/3/98 - 2D output changed |version 71 - 1/6/98 - matching section removed |version 72 - 3/7/98 - calibration bug in batch_end fixed |version 73 - 3/8/98 - modified for the new element specific output from vegas |version 74 - 5/8/98 - check included for matching file length between data and simulation |version 75 - 27/8/98 - bug fix in 2D layout macro ("error cursor not on top graph") - not true!!!!! |version 76 - 10/9/98 - order of parameters changed to match the order in Vegas |version 77 - 24/9/98 - moved from mac to pc platform |version 78 - 23/10/98 - fixed (hopefully) "error cursor not on top graph" |version 79 - 9/12/98 - contour axes corrected + few minor mods |version 80 - 23/3/99 - "error cursor not on top graph" fixed for 3D simulations (still don't understand it) |version 81 - 23/3/99 - de-spiking option included |version 82 - 23/3/99 - smoothing of simulation can be user-defined |version 83 - 14/4/99 - 3D summary prinout changed |version 84 - 15/4/99 - 3D contours tweaked |version 85 - 21/4/99 - improved de-spiking |version 86 - 22/4/99 - few minor changes |version 87 - 30/4/99 - changed contour levels to auto |version 88 - 19/5/99 - moved skew and calib into subroutines |version 89 - 20/5/99 - STILL found things to change! |version 90 - 2/6/99 - Element specific options introduced. |version 91 - 4/6/99 - Bugs introduced in versions 88 onwards fixed! |version 92 - 7/6/99 - Minor bugs fixed. |version 93 - 8/6/99 - Control buttons moved to LHS of comparison graph. |version 94 - 10/6/99 - Very subtle bug destroyed! |version 95 - 22/7/99 - Modified to accept just one multimodel parameter. |version 96 - 12/8/99 - Modified to cope with eclipse. #pragma rtglobals=0 //sod this silly new global parameter stuff!! Menu "Tidy" "-" "tidy_now/0" end Macro batch_set_up() silent 1; pauseupdate variable M string pathname, filename if (k12==0) k12=1 endif if (k13==0) k13=1 endif pathname="temppath" filename=indexedfile($pathname, 0, "????") print "filename=" print filename Loadwave/Q/G/D/O/W/P=temppath/N=sim/L={0,0,0,0,6} filename sim0=expt_bin0 duplicate/O sim1 exp1, exp_fit, sim_fit, sim_smooth dowindow/k comparison2 dowindow/k comparison dowindow/k comparison2_1 dowindow/k comparison_1 execute "comparison()" if (numpnts(expt_bin1)!=numpnts(sim1)) tidy_now() beep abort "the data and simulation files are not the same length" endif if ((element_1flag==0) %& (element_2flag==0)) beep abort "you have to select a least one element!" endif if ((element_1flag==1) %& (element_2flag==1)) sim1=sim1+sim2 //using both elements endif if ((element_1flag==0) %& (element_2flag==1)) sim1=sim3 //using second element only endif exp1=expt_bin1 duplicate/O exp1, exp_orig Ne1=2 Ne2=numpnts(exp1)-3 if (spikeflag>0) spike() endif if (smoothflag>0) Smooth smooth_amount, sim1 endif WaveStats/Q/R=[Ne1, Ne2] exp1 exp_avg=V_avg WaveStats/Q/R=[Ne1, Ne2] sim1 Sim_Avg=V_avg exp1*=Sim_avg/exp_avg WaveStats/Q/R=[Ne1, Ne2] sim0 Ang_Avg=V_avg Legend/K/N=text0 Legend/J/N=text0/F=0/A=MC/X=18.13/Y=-43.81 "\\JL\\Z09\\s(sim1) simulation \\s(exp1) experiment" if (calibflag>0) AppendText "scaling factor = \\{scale/calibration}, calibration range = \\{calibration_range}%" endif cursor/C=(0,0,65000)/P A, exp1, Ne1 cursor/C=(0,65000,0)/P B, exp1, Ne2 N_over=ne2-ne1+1 if (skewflag>0) skew() endif if (calibflag>0) calib() endif DoWindow/B simulation_summary DoWindow/B match_data DoWindow/F R_Value DoWindow/F Comparison KillWaves/A/Z/F end macro batch_start(ctrlName) : ButtonControl String ctrlName, filename, pathname, str1 pauseupdate silent 1 spikedetect=0 Ne1=pcsr(A) Ne2=pcsr(B) N_over=Ne2-Ne1+1 WaveStats/Q/R=[Ne1, Ne2] sim0 Ang_Avg=V_avg DoWindow/K R_Value DoWindow/K Comparison R_Value2() comparison2() DoWindow/B contour12 DoWindow/B contour23 DoWindow/B contour13 print "Started on "+date()+" at "+time() batchprocess() end Function Calc_RValue_is() variable N, R R=0 N=Ne1-1 do N+=1 R+=((sim1[N]-exp1[N])/exp1[N])^2 while (Nmax) max=sim1[N] endif while (Nmax) max=exp1[N] endif while (N0) spike() endif if (smoothflag>0) Smooth smooth_amount, sim1 endif sprintf sim_name, "Veg" +"%g" index duplicate/O sim1 $sim_name //creates wave from modified simulation WaveStats/Q/R=[Ne1, Ne2] sim1 sim_avg=v_avg WaveStats/Q/R=[Ne1, Ne2] exp_orig exp_avg=V_avg scale=exp_avg/sim_avg //scale is the emperical calibration exp1=exp_orig/scale //scales experimental data - simulation is left unchanged if (skewflag>0) skew() endif if (calibflag>0) calib() endif Calc_RValue_Xsq() index+=1 doupdate while (index<(filenum+1)) simdata=S_Filename Execute "batch_end()" End Macro batch_end() print "Finished on "+date()+" at "+time() if (spikeflag>0) print "% of spikes detected per file= ", (spikedetect/index/(Ne2-Ne1)*100) endif tidy_now() silent 1 variable N pauseupdate wavestats/Q RValue sim_num=V_minloc //sim_num points to simulation with lowest R Value String sim_name sprintf sim_name, "Veg"+"%g" sim_num sim1=$sim_name WaveStats/Q/R=[Ne1, Ne2] sim1 sim_avg=v_avg WaveStats/Q/R=[NE1,Ne2] exp_orig exp_avg=v_avg scale=exp_avg/sim_avg //scale is the emperical calibration exp1=exp_orig/scale //scales experimental data - simulation is left unchanged DoWindow/F Comparison if (skewflag>0) skew() endif if (calibflag>0) calib() endif DoWindow/F R_Value if ((K13>=2) %& (K12>=2)) contour_3D() RValue_Graphs_3D() endif if ((K13==1) %& (K12>=2)) contour_2D() RValue_Graphs_2D() endif if ((K13==1) %& (K12==1)) RValue_Graphs_1D() endif Rv=RValue[sim_num-1] if (printflag==1) printout() endif tidy_now() end Macro single() silent 1 pauseupdate string pathname, filename pathname="temppath" Loadwave/Q/G/D/O/W/P=temppath/N=sim/L={0,0,0,0,6} filename duplicate/O sim1 exp1, exp_fit, sim_fit exp1=expt_bin1 duplicate/O exp1, exp_orig Ne1=2 Ne2=numpnts(exp1)-3 WaveStats/Q/R=[Ne1, Ne2] exp1 exp_avg=V_avg WaveStats/Q/R=[Ne1, Ne2] sim1 Sim_Avg=V_avg exp1*=Sim_avg/exp_avg WaveStats/Q/R=[Ne1, Ne2] sim0 Ang_Avg=V_avg cursor/C=(0,0,65000)/P A, exp1, Ne1 cursor/C=(0,65000,0)/P B, exp1, Ne2 N_over=ne2-ne1+1 DoWindow/B simulation_summary DoWindow/B match_data DoWindow/F R_Value DoWindow/F Comparison KillWaves/A/Z/F if ((element_1flag==0) %& (element_2flag==0)) beep abort "you have to select a least one element!" endif if ((element_1flag==1) %& (element_2flag==1)) sim1=sim1+sim2 //using both elements endif if ((element_1flag==0) %& (element_2flag==1)) sim1=sim2 //using second element only endif if (smoothflag>0) Smooth smooth_amount, sim1 endif WaveStats/Q/R=[Ne1, Ne2] sim1 exp1=exp1*V_Avg/exp_avg //scales experimental data - simulation is left unchanged if (skewflag>0) skew() endif if (calibflag>0) calib() endif Calc_RValue_Xsq() doupdate PRINT "the rvalue is", Rv end Macro Printout() silent 1 Summary_Layout_3D() printlayout Summary_Layout_3D DoWindow /K Summary_Layout_3D Return Window r_value() : Graph PauseUpdate; Silent 1 Display /W=(202,250,460,370) RValue as "progress chart" ModifyGraph zero(bottom)=1 ModifyGraph nticks(bottom)=10 ModifyGraph minor(bottom)=1 Label left "R-value" Label bottom "number of files" Wavestats/Q RValue Textbox/N=text0/F=0/A=MC/X=-6.27/Y=39.39 "min. R-value = \\{Rv}" ShowInfo ControlBar 17 Button button0,pos={5,1},size={91,16},proc=save_rvalue,title="save R-value" Button button1,pos={122,1},size={90,16},proc=load_Rvalue,title="load R-value" Button button2,pos={244,1},size={52,16},proc=revert_Rvalue,title="revert" EndMacro Window comparison() : Graph PauseUpdate; Silent 1 Display /W=(0,30,460,230) sim1 vs sim0 as "Comparison of experiment and simulation" Append exp1 vs sim0 ModifyGraph mode(exp1)=4 ModifyGraph marker(exp1)=19 ModifyGraph lSize(sim1)=2 ModifyGraph rgb(exp1)=(0,0,0) ModifyGraph nticks(bottom)=10 Label left "scattering yield [layers]" Label bottom "scattering angle [degrees]" Legend/J/N=text0/F=0/A=MC/X=35/Y=-48 "\\JL\\Z09\\s(sim1) simulation \\s(exp1) experiment" if (calibflag>0) AppendText "scaling factor = \\{scale/calibration}, calibration range = \\{calibration_range}%" endif cursor/C=(0,0,65000)/P A, exp1, Ne1 cursor/C=(0,65000,0)/P B, exp1, Ne2 Button Restart,pos={69,30},size={55,17},proc=Vegas,title="Restart" SetVariable Browse,pos={129,7},size={110,18},proc=Select_current_sim,title="Get sim." SetVariable Browse,limits={1,Inf,1},value= sim_num Button start,pos={69,51},size={55,17},proc=batch_start,title="Run" Button reset,pos={69,8},size={55,17},proc=reset,title="Reset" EndMacro Function reset(ctrlName) : ButtonControl String ctrlName execute "batch_end()" End Window simulation_summary() : Panel PauseUpdate; Silent 1 // building window... NewPanel /W=(0,42.8,640.2,483.2) as "Simulation analysis" SetDrawLayer UserBack SetDrawEnv fsize= 14,fstyle= 1 DrawText 33,31,"This program produces R-values from Vegas simulations" SetDrawEnv fstyle= 1 DrawText 30,80,"1." SetDrawEnv fstyle= 1 DrawText 30,110,"2. " SetDrawEnv fstyle= 1 DrawText 30,140,"3. " SetDrawEnv fstyle= 1 DrawText 30,170,"4. " SetDrawEnv fstyle= 1 DrawText 30,200,"5. " SetDrawEnv fstyle= 1 DrawText 30,230,"6. " SetDrawEnv fstyle= 1 DrawText 30,260,"7." SetDrawEnv fstyle= 1 DrawText 30,290,"8." SetDrawEnv fstyle= 1 DrawText 30,320,"9." SetDrawEnv fstyle= 1 DrawText 25,350,"10." SetDrawEnv fstyle= 1 DrawText 25,380,"11." SetDrawEnv fstyle= 1 DrawText 25,410,"12." SetDrawEnv fstyle= 1 DrawText 315,380,"13." SetDrawEnv fstyle= 1 DrawText 316,411,"14." SetVariable setvar0,pos={100,66},size={288,18},title="Experimental data; " SetVariable setvar0,frame=0,limits={-Inf,Inf,1},value= expdata SetVariable setvar1,pos={102,95},size={474,18},title="Simulation folder; " SetVariable setvar1,frame=0,limits={-Inf,Inf,1},value= simfolder Button button0,pos={497,291},size={91,53},proc=run,title="Run" CheckBox check0l,pos={50,123},size={185,20},proc=multi,title="Multiple simulation files",value=1 CheckBox check1,pos={50,306},size={147,20},proc=printbest,title="Print best result?",value=0 CheckBox check2,pos={50,243},size={147,20},proc=smoothset,title="Smooth simulation",value=1 CheckBox check3,pos={50,273},size={147,20},proc=spikeset,title="De-spike simulation",value=1 CheckBox check4,pos={50,212},size={147,20},proc=skewset,title="Skew data",value=1 CheckBox check4_1,pos={50,152},size={135,20},proc=element_1set,title="Use first element",value=0 CheckBox check4_2,pos={185,152},size={147,20},proc=element_2set,title="Use second element",value=1 CheckBox calibflag,pos={50,182},size={122,20},proc=calibset,title="Use calibration",value=0 SetVariable setvar4,pos={50,333},size={225,18},title="Multi parameter number 1 ;" SetVariable setvar4,limits={1,100,1},value= K11 SetVariable setvar5,pos={50,363},size={225,18},title="Multi parameter number 2 ;" SetVariable setvar5,limits={1,100,1},value= K12 SetVariable setvar6,pos={50,393},size={225,18},title="Multi parameter number 3 ;" SetVariable setvar6,limits={1,100,1},value= K13 Button find_expt,pos={45,61},size={50,20},proc=find_expt,title="find" Button find_sim,pos={46,94},size={50,20},proc=find_sim,title="find" SetVariable setcalib,pos={175,182},size={215,18},proc=SetCalibration,title="Value: counts per layer ;" SetVariable setcalib,limits={-Inf,Inf,1},value= calibration SetVariable set_calib_range,pos={399,182},size={94,18},title="% range" SetVariable set_calib_range,limits={-Inf,Inf,1},value= calibration_range SetVariable setvar7,pos={199,273},size={207,18},title="set spike tolerance level" SetVariable setvar7,limits={0.01,2,0.05},value= spikelevel SetVariable setvar8,pos={199,244},size={127,18},title="set amount " SetVariable setvar8,limits={1,Inf,2},value= smooth_amount SetVariable setvar2,pos={337,394},size={238,18},title="Parameter polynomial fit order" SetVariable setvar2,limits={3,9,2},value= poly_amount SetVariable setvar3,pos={337,364},size={238,18},title="Skew range overlap " SetVariable setvar3,limits={0,25,1},value= Ne_extra EndMacro Macro element_2set(ctrlName, flag) : checkboxControl String ctrlName variable flag element_2flag=flag EndMacro Macro element_1set(ctrlName, flag) : checkboxControl String ctrlName variable flag element_1flag=flag EndMacro macro find_expt(ctrlName) : ButtonControl String ctrlName loadwave /G/O/A/Q/N=expt_bin expdata=S_filename exp_show() End macro find_sim(ctrlName) : ButtonControl String ctrlName newpath/O temppath pathinfo temppath simfolder=S_path End Macro multi(ctrlName, flag) : checkboxControl String ctrlName variable flag multiflag=flag EndMacro Macro printbest(ctrlName, flag) : checkboxControl String ctrlName variable flag printflag=flag EndMacro Macro smoothset(ctrlName, flag) : checkboxControl String ctrlName variable flag smoothflag=flag EndMacro Macro spikeset(ctrlName, flag) : checkboxControl String ctrlName variable flag spikeflag=flag EndMacro Macro calibset(ctrlName, flag) : checkboxControl String ctrlName variable flag calibflag=flag EndMacro Macro skewset(ctrlName, flag) : checkboxControl String ctrlName variable flag skewflag=flag EndMacro Macro run(ctrlName) : ButtonControl String ctrlName DoWindow/B Simulation_summary DoWindow/B Match_data if (multiflag==1) batch_set_up() else single() endif EndMacro Function ButtonRestart(ctrlName) : ButtonControl String ctrlName dowindow/f Match_data End //________________________________________________________________ // START OF CONTOUR STUFF //________________________________________________________________ //K11 is number of values of multi-model paramter 1 //K12 is number of values of multi-model paramter 2 //K13 is number of values of multi-model paramter 3 macro test() K16=trunc(k19/K11/k12) K15=trunc((k19-K16*K11*K12)/K11) K14=k19-K15*K11-K16*K12*K11 print k19 print k14 print k15 print k16 end function contour_2D() variable N, M wavestats/q rvalue V_minloc-=1 K14=trunc(V_minloc/K12) K15=V_minloc-K14*K12 Make/N=(K12, K11)/O matrix12 N=0 M=0 do | picks K11 and K12 do matrix12[N] [M]=rvalue[(N+M*K12)] N+=1 while (N0) skew() endif if (calibflag>0) calib() endif End Macro Vegas(ctrlName) : ButtonControl String ctrlName DoWindow/K exp_show DoWindow/F Simulation_summary End Window comparison2() : Graph PauseUpdate; Silent 1 Display /W=(0,30,680,320) sim1 vs sim0 as "Comparison of experiment and simulation" Append exp1 vs sim0 ModifyGraph mode(exp1)=4 ModifyGraph marker(exp1)=19 ModifyGraph lSize(sim1)=2 ModifyGraph rgb(exp1)=(0,0,0) ModifyGraph nticks(bottom)=10 ModifyGraph lblMargin(bottom)=3 ModifyGraph lblLatPos(bottom)=-1 Textbox/N=text0/S=3/A=MC/X=32.14/Y=-44.25 "simulation number = \\{index}" cursor/C=(0,0,65000)/P A, exp1, Ne1 cursor/C=(0,65000,0)/P B, exp1, Ne2 EndMacro Window R_Value2() : Graph PauseUpdate; Silent 1 Display /W=(170,330,510,480) RValue as "progress chart" ModifyGraph zero(bottom)=1 ModifyGraph nticks(bottom)=10 ModifyGraph minor(bottom)=1 SetAxis/A/E=1 bottom EndMacro function tidy_now() PauseUpdate; Silent 1 dowindow/k comparison2 dowindow/k simulation_summary dowindow/k r_value2 dowindow/k comparison dowindow/k r_value dowindow/k comparison2_1 dowindow/k r_value2_1 dowindow/k comparison_1 dowindow/k r_value_1 dowindow/k Graph0 dowindow/k Graph1 dowindow/k table0 dowindow/k table1 execute "simulation_summary()" wavestats/Q RValue execute "r_value()" cursor/p A, RValue, (sim_num-1) execute "comparison()" end function chi_square_precision_3D() variable N, sigma execute "CurveFit/Q/L=(numpnts(p3_rvalue)*10) poly poly_amount, P3_RValue /D " wavestats/Q fit_p3_rvalue p3_min=V_minloc Duplicate/O fit_p3_rvalue diff_p3_rvalue Differentiate diff_p3_rvalue Differentiate diff_p3_rvalue sigma=(abs(2/ diff_p3_rvalue[9+10*trunc(V_minloc)]))^0.5 p3_sigma=sigma dowindow/F P3_Rvalue_graph cursor A, fit_p3_rvalue, V_minloc-sigma cursor B, fit_p3_rvalue, V_minloc+sigma Textbox/N=text0/F=0/A=MC/X=10/Y=40 "\\JCmin. value@ \\{p3_min}\rsigma = \\{p3_sigma}" Tag/N=text1/F=0/X=(xcsr(A)-6.7)/Y=50 fit_P3_RValue, xcsr(A), "- sigma" Tag/N=text2/F=0/X=(xcsr(B)-6.7)/Y=50 fit_P3_RValue, xcsr(B), "+ sigma" execute "CurveFit/Q/L=(numpnts(p1_rvalue)*10) poly poly_amount, P1_RValue /D " wavestats/Q fit_p1_rvalue p1_min=V_minloc Duplicate/O fit_p1_rvalue diff_p1_rvalue Differentiate diff_p1_rvalue Differentiate diff_p1_rvalue sigma=(abs(2/ diff_p1_rvalue[9+10*trunc(V_minloc)]))^0.5 p1_sigma=sigma dowindow/F P1_Rvalue_graph cursor A, fit_p1_rvalue, V_minloc-sigma cursor B, fit_p1_rvalue, V_minloc+sigma Textbox/N=text0/F=0/A=MC/X=3/Y=40 "\\JCmin. value@ \\{p1_min}\rsigma = \\{p1_sigma}" Tag/N=text1/F=0/X=(xcsr(A)-6.7)/Y=50 fit_P1_RValue, xcsr(A), "- sigma" Tag/N=text2/F=0/X=(xcsr(B)-6.7)/Y=50 fit_P1_RValue, xcsr(B), "+ sigma" execute "CurveFit/Q/L=(numpnts(p2_rvalue)*10) poly poly_amount, P2_RValue /D " wavestats/Q fit_p2_rvalue p2_min=V_minloc Duplicate/O fit_p2_rvalue diff_p2_rvalue Differentiate diff_p2_rvalue Differentiate diff_p2_rvalue sigma=(abs(2/ diff_p2_rvalue[9+10*trunc(V_minloc)]))^0.5 p2_sigma=sigma dowindow/F P2_Rvalue_graph cursor A, fit_p2_rvalue, V_minloc-sigma cursor B, fit_p2_rvalue, V_minloc+sigma Textbox/N=text0/F=0/A=MC/X=10/Y=40 "\\JCmin. value@ \\{p2_min}\rsigma = \\{p2_sigma}" Tag/N=text1/F=0/X=(xcsr(A)-6.7)/Y=50 fit_P2_RValue, xcsr(A), "- sigma" Tag/N=text2/F=0/X=(xcsr(B)-6.7)/Y=50 fit_P2_RValue, xcsr(B), "+ sigma" end function chi_square_precision_2D() variable N, sigma execute "CurveFit/Q/L=(numpnts(p2_rvalue)*10) poly poly_amount, P2_RValue /D " wavestats/Q fit_p2_rvalue p2_min=V_minloc Duplicate/O fit_p2_rvalue diff_p2_rvalue Differentiate diff_p2_rvalue Differentiate diff_p2_rvalue sigma=(abs(2/ diff_p2_rvalue[9+10*trunc(V_minloc)]))^0.5 p2_sigma=sigma dowindow/F P2_Rvalue_graph cursor A, fit_p2_rvalue, V_minloc-sigma cursor B, fit_p2_rvalue, V_minloc+sigma Textbox/N=text0/F=0/A=MC/X=3/Y=40 "\\JCmin. value@ \\{p2_min}\rsigma = \\{p2_sigma}" Tag/N=text1/F=0/X=(xcsr(A)-6.7)/Y=50 fit_P2_RValue, xcsr(A), "- sigma" Tag/N=text2/F=0/X=(xcsr(B)-6.7)/Y=50 fit_P2_RValue, xcsr(B), "+ sigma" execute "CurveFit/Q/L=(numpnts(p1_rvalue)*10) poly poly_amount, P1_RValue /D " wavestats/Q fit_p1_rvalue p1_min=V_minloc Duplicate/O fit_p1_rvalue diff_p1_rvalue Differentiate diff_p1_rvalue Differentiate diff_p1_rvalue sigma=(abs(2/ diff_p1_rvalue[9+10*trunc(V_minloc)]))^0.5 p1_sigma=sigma dowindow/F P1_Rvalue_graph cursor A, fit_p1_rvalue, V_minloc-sigma cursor B, fit_p1_rvalue, V_minloc+sigma Textbox/N=text0/F=0/A=MC/X=3/Y=40 "\\JCmin. value@ \\{p1_min}\rsigma = \\{p1_sigma}" Tag/N=text1/F=0/X=(xcsr(A)-6.7)/Y=50 fit_P1_RValue, xcsr(A), "- sigma" Tag/N=text2/F=0/X=(xcsr(B)-6.7)/Y=50 fit_P1_RValue, xcsr(B), "+ sigma" end function chi_square_precision_1D() variable sigma p1_rvalue=rvalue execute "CurveFit/Q/L=(numpnts(p1_rvalue)*10) poly poly_amount, P1_RValue /D " wavestats/Q fit_p1_rvalue p1_min=V_minloc Duplicate/O fit_p1_rvalue diff_p1_rvalue Differentiate diff_p1_rvalue Differentiate diff_p1_rvalue sigma=(abs(2/ diff_p1_rvalue[9+10*trunc(V_minloc)]))^0.5 p1_sigma=sigma dowindow/F P1_Rvalue_graph cursor A, fit_p1_rvalue, V_minloc-sigma cursor B, fit_p1_rvalue, V_minloc+sigma Textbox/N=text0/F=0/A=MC/X=3/Y=40 "\\JCmin. value@ \\{p1_min}\rsigma = \\{p1_sigma}" Tag/N=text1/F=0/X=(xcsr(A)-6.7)/Y=50 fit_P1_RValue, xcsr(A), "- sigma" Tag/N=text2/F=0/X=(xcsr(B)-6.7)/Y=50 fit_P1_RValue, xcsr(B), "+ sigma" end Function save_rvalue(ctrlName) : ButtonControl String ctrlName Save/G RValue as "RValue" End Function load_Rvalue(ctrlName) : ButtonControl String ctrlName ctrlName = S_Filename Duplicate/O RValue Rvalue_orig loadwave /G/O/A/Q/N=Rvalue Rvalue+=Rvalue0 Rvalue/=2 S_Filename = ctrlName execute "batch_end()" End Function revert_Rvalue(ctrlName) : ButtonControl String ctrlName Rvalue=Rvalue_orig execute "batch_end()" End Function calc_rvalue(ctrlName) : ButtonControl String ctrlName Ne1=pcsr(A) Ne2=pcsr(B) N_over=Ne2-Ne1+1 Calc_RValue_Xsq() execute "print rvalue(index)" End Function SetCalibration(ctrlName,varNum,varStr,varName) : SetVariableControl String ctrlName Variable varNum String varStr String varName End Window dummy_graph() : Graph PauseUpdate; Silent 1 Display /W=(5,42,400,250) expt_bin1 as "dummy" EndMacro Window Summary_layout_3D() : Layout PauseUpdate; Silent 1 // building window... Layout/C=1/W=(4.8,42.2,573.6,445.4) Contour_12(601.8,226.2,768.6,369.6)/O=1,comparison(72.6,72.6,414.6,216)/O=1 Append r_value(72.6,226.2,415.8,369.6)/O=1,P1_RValue_Graph(425.4,72.6,592.2,216)/O=1 Append P2_RValue_Graph(425.4,226.2,592.2,369.6)/O=1,P3_RValue_Graph(425.4,378.6,592.2,522)/O=1 Append Contour_31(601.8,72.6,768.6,216)/O=1,Contour_23(601.8,378.6,768.6,522)/O=1 Textbox/N=text0/F=0/A=LB/X=0.60/Y=0.8 "\\Z12Printed on "+ date() + "@" + time() AppendText "simulation file; "+simdata+".\{sim_num}" AppendText "exptl. data file; "+expdata sprintf str1 "The R-value is %g", Rv AppendText str1 sprintf str1 "For %g data points", (ne2-ne1) //sprintf str1 "The R-value-2 is %g", ((RValue[k8-1]-1)/avg_y*1e3) AppendText str1 sprintf str1 "from \{hcsr(A)} to \{hcsr(B)} degrees" AppendText str1 if (calibflag>0) AppendText "Using absolute calibration" else AppendText "Not using absolute calibration" endif if (spikeflag>0) sprintf str1 "percentage of spikes in file=\{spikedetect/index/(Ne2-Ne1)*100}" AppendText str1 sprintf str1 "spike tolerance range=\{spikelevel}" AppendText str1 else AppendText "Simulation not de-spiked" endif if (smoothflag>0) sprintf str1 "Simulation smoothed, value = \{smooth_amount}" AppendText str1 else AppendText "Simulation not smoothed" endif if (skewflag>0) sprintf str1 "Simulation skewed" AppendText str1 else AppendText "Simulation not skewed" endif EndMacro Function spike() variable M=0 do sim_smooth=sim1 smooth 15, sim_smooth if (sim1[M+1]>(sim_smooth[M]+spikelevel)) // checks for a greater than 'spikelevel' increase.... sim1[M+1]=sim_smooth[M+1] // between adjacent points spikedetect+=1 endif M+=1 while ((M+1)(1+calibration_range/100))) |sets the upper limit to the range of calibration allowed exp1=exp1*scale/(calibration*(1+calibration_range/100)) endif if (scale/calibration<(1-calibration_range/100))) |sets the lower limit to the range of calibration allowed exp1=exp1*scale/(calibration*(1-calibration_range/100)) endif end Function skew() CurveFit/Q poly 3, exp1 ((xcsr(A)-Ne_extra),(xcsr(B)+Ne_extra)) /X=sim0 /D=exp_fit exp_fit= W_coef[0]+ W_coef[1]*sim0+ W_coef[2]*(sim0^2) CurveFit/Q poly 3, sim1 ((xcsr(A)-Ne_extra),(xcsr(B)+Ne_extra)) /X=sim0 /D=sim_fit sim_fit= W_coef[0]+ W_coef[1]*sim0+ W_coef[2]*(sim0^2) exp1*=sim_fit/exp_fit end Window exp_show() : Graph PauseUpdate; Silent 1 Display /W=(0,40,640,480) expt_bin1 vs expt_bin0 as "Experimental data" ModifyGraph mode=4 ModifyGraph marker=19 ModifyGraph lStyle=1 ModifyGraph nticks(bottom)=10 Label left "scattering yield corrected for sin\\S4\\M\\F'Symbol'(q/2)" Label bottom "scattering angle [degrees]" Textbox/N=text0/F=0/S=2/A=MC/X=-7.93/Y=47.96 "Data from file - " + expdata Button exp_OK,pos={523,209},size={70,20},proc=Vegas,title="Continue" EndMacro