#set multiplot layout 2,1 rootdir="./" square="sq_pulse/" dir1="/processed_histogram_with_yerror_bars/" exp="exp_pulse/vary_photon_number_perkin_elmer_" #### 5 ns pulses #### exp_p_5ns_123_filename=exp."5ns".dir1."OD_in_0.32_OD_out_2.2_2.9_after_4_histo" exp_p_5ns_5_filename=exp."5ns".dir1."OD_in_2.3_OD_out_1.4_1.1_after_histo" exp_p_5ns_562_filename=exp."5ns".dir1."OD_in_0.16_OD_out_2.2_2.9_after_histo" exp_p_5ns_1028_filename=exp."5ns".dir1."OD_in_0.16_OD_out_2.3_2.9_after_4_histo" sq_p_5ns_117_filename=square."5ns".dir1."OD_in_0.98_OD_out_1.1_2.9_after_histo" sq_p_5ns_6_filename=square."5ns".dir1."OD_in_2.3_OD_out_2.9_after_histo" sq_p_5ns_43_filename=square."5ns".dir1."OD_in_1.4_OD_out_1.1_2.9_after_histo" sq_p_5ns_523_filename=square."5ns".dir1."OD_in_0.32_OD_out_2.2_2.9_after_2_histo" sq_p_5ns_980_filename=square."5ns".dir1."OD_in_0.00_OD_out_2.2_2.9_after_histo" #### 15 ns pulses #### exp_p_15ns_6_filename=exp."15ns".dir1."OD_in_1.8_0.51_OD_out_2.9_after_histo" exp_p_15ns_104_filename=exp."15ns".dir1."OD_in_1.1_OD_out_1.4_2.9_after_histo" exp_p_15ns_241_filename=exp."15ns".dir1."OD_in_0.76_OD_out_1.4_2.9_after_histo" exp_p_15ns_637_filename=exp."15ns".dir1."OD_in_0.16_OD_out_1.8_2.9_after_histo" exp_p_15ns_342_filename=exp."15ns".dir1."OD_in_0.32_OD_out_1.8_2.9_after_histo" exp_p_15ns_1229_filename=exp."15ns".dir1."OD_in_0.00_OD_out_1.8_2.9_after_histo" sq_p_15ns_6_filename=square."15ns".dir1."OD_in_0.76_2.3_OD_out_2.9_after_histo" sq_p_15ns_110_filename=square."15ns".dir1."OD_in_1.8_OD_out_1.1_2.9_after_histo" sq_p_15ns_256_filename=square."15ns".dir1."OD_in_1.4_OD_out_1.1_2.9_after_histo" sq_p_15ns_538_filename=square."15ns".dir1."OD_in_0.98_OD_out_2.2_2.9_after_histo" sq_p_15ns_374_filename=square."15ns".dir1."OD_in_1.1_OD_out_2.2_2.9_after_histo" sq_p_15ns_1456_filename=square."15ns".dir1."OD_in_0.51_OD_out_2.2_2.9_after_histo" #### 60ns pulses #### exp_p_60ns_4_5_filename=exp."60ns".dir1."OD_in_1.1_1.5_OD_out_2.9_after_histo" exp_p_60ns_90_filename=exp."60ns".dir1."OD_in_1.4_OD_out_0.51_2.9_after_histo" exp_p_60ns_913_filename=exp."60ns".dir1."OD_in_0.51_OD_out_1.5_2.9_after_histo" sq_p_60ns_4_5_filename=square."60ns".dir1."OD_in_1.1_2.0_OD_out_2.9_after_3_histo" sq_p_60ns_80_filename=square."60ns".dir1."OD_in_0.32_2.0_OD_out_2.2_2.9_after_histo" sq_p_60ns_1260_filename=square."60ns".dir1."OD_in_0.32_1.0_OD_out_2.2_2.9_after_histo" #### Columns for pulse data file #### bin=1 time=2 at_tran_count=3 at_refl_count=4 at_tran_prob=5 at_refl_prob=6 at_tran_prob_err=7 at_refl_prob_err=8 bg_tran_count=9 bg_refl_count=10 bg_tran_prob=11 bg_refl_prob=12 bg_tran_prob_err=13 bg_refl_prob_err=14 #### Styles #### set style line 1 lc 1 lw 1 ps 1.5 set style line 2 lc 3 lw 1 ps 1.5 set style line 3 lc 4 lw 2 ps 1.5 set style line 4 lc 2 lw 3 ps 1.5 #### Some experimental parameters #### exp_trigger_delay=555 exp_refl_delay=295+exp_trigger_delay exp_refl_delay2=480+exp_refl_delay sq_trigger_delay=324 sq_refl_delay=302+sq_trigger_delay tau=26.24 set xrange [-5:5] #set yrange [1e-7:] set arrow 1 from graph 0.5,0 to graph 0.5,1 ls 0 lw 2 nohead # plot exp_p_5ns_123_filename u (column(time)-exp_trigger_delay)/tau:bg_tran_prob:bg_tran_prob_err w yerrorl,\ # exp_p_5ns_123_filename u (column(time)-exp_refl_delay)/tau:at_refl_prob:at_refl_prob_err w yerrorl,\ # exp_p_5ns_5_filename u (column(time)-exp_trigger_delay)/tau:bg_tran_prob:bg_tran_prob_err w yerrorl,\ # exp_p_5ns_5_filename u (column(time)-exp_refl_delay)/tau:at_refl_prob:at_refl_prob_err w yerrorl #unset key #set logscale y # set title "t=5ns;n=123,117" # plot exp_p_5ns_123_filename u (column(time)-exp_refl_delay)/tau:at_refl_prob:at_refl_prob_err w yerrorl t 'exp' ls 2,\ # sq_p_5ns_117_filename u (column(time)-sq_refl_delay)/tau:at_refl_prob:at_refl_prob_err w yerrorl t 'square' ls 1 #pause -1 #set title "t=5ns;n=562,523" #plot exp_p_5ns_562_filename u (column(time)-exp_refl_delay)/tau:at_refl_prob:at_refl_prob_err w yerrorl t 'exp' ls 2,\ # sq_p_5ns_523_filename u (column(time)-sq_refl_delay)/tau:at_refl_prob:at_refl_prob_err w yerrorl t 'square' ls 1 #pause -1 # set title "t=5ns;n=5,6" # plot exp_p_5ns_5_filename u (column(time)-exp_refl_delay)/tau:at_refl_prob:at_refl_prob_err w yerrorl t 'exp' ls 2,\ # sq_p_5ns_6_filename u (column(time)-sq_refl_delay)/tau:at_refl_prob:at_refl_prob_err w yerrorl t 'square' ls 1 # set title "t=5ns;n=562,523" # plot exp_p_5ns_562_filename u (column(time)-exp_refl_delay)/tau:at_refl_prob:at_refl_prob_err w yerrorl t 'exp' ls 2,\ # sq_p_5ns_523_filename u (column(time)-sq_refl_delay)/tau:at_refl_prob:at_refl_prob_err w yerrorl t 'square' ls 1 # set title "t=5ns;n=1028,980" # plot exp_p_5ns_1028_filename u (column(time)-exp_refl_delay)/tau:at_refl_prob:at_refl_prob_err w yerrorl t 'exp' ls 2,\ # sq_p_5ns_980_filename u (column(time)-sq_refl_delay)/tau:at_refl_prob:at_refl_prob_err w yerrorl t 'square' ls 1 # set title "t=15ns;n=6.3,6.3" # plot exp_p_15ns_6_filename u (column(time)-exp_refl_delay)/tau:at_refl_prob:at_refl_prob_err w yerrorl t 'exp' ls 2,\ # sq_p_15ns_6_filename u (column(time)-sq_refl_delay)/tau:at_refl_prob:at_refl_prob_err w yerrorl t 'square' ls 1 # set title "t=15ns;n=104,110" # plot exp_p_15ns_104_filename u (column(time)-exp_refl_delay)/tau:at_refl_prob:at_refl_prob_err w yerrorl t 'exp' ls 2,\ # sq_p_15ns_110_filename u (column(time)-sq_refl_delay)/tau:at_refl_prob:at_refl_prob_err w yerrorl t 'square' ls 1 # set title "t=15ns;n=241,256" # plot exp_p_15ns_241_filename u (column(time)-exp_refl_delay)/tau:at_refl_prob:at_refl_prob_err w yerrorl t 'exp' ls 2,\ # sq_p_15ns_256_filename u (column(time)-sq_refl_delay)/tau:at_refl_prob:at_refl_prob_err w yerrorl t 'square' ls 1 # set title "t=15ns;n=342,374" # plot exp_p_15ns_342_filename u (column(time)-exp_refl_delay)/tau:at_refl_prob:at_refl_prob_err w yerrorl t 'exp' ls 2,\ # sq_p_15ns_374_filename u (column(time)-sq_refl_delay)/tau:at_refl_prob:at_refl_prob_err w yerrorl t 'square' ls 1 # set title "t=15ns;n=637,538" # plot exp_p_15ns_637_filename u (column(time)-exp_refl_delay)/tau:at_refl_prob:at_refl_prob_err w yerrorl t 'exp' ls 2,\ # sq_p_15ns_538_filename u (column(time)-sq_refl_delay)/tau:at_refl_prob:at_refl_prob_err w yerrorl t 'square' ls 1 # set title "t=15ns;n=1229,1456" # plot exp_p_15ns_1229_filename u (column(time)-exp_refl_delay)/tau:at_refl_prob:at_refl_prob_err w yerrorl t 'exp' ls 2,\ # sq_p_15ns_1456_filename u (column(time)-sq_refl_delay)/tau:at_refl_prob:at_refl_prob_err w yerrorl t 'square' ls 1 #### Example pulse fit #### set terminal postscript enhanced color eps set output '../pulses/excitation_fluo_pulse_15ns_n100.eps' f(x)= eta*lambda*exp((x-x0)/tau0)*(exp(time_width/tau0)-1) x0=exp_trigger_delay; time_width=1 tau0=15.3907 eta=0.3 lambda=0.00381085*14750000 g(x)= eta*lambda1*exp((x1-x)/tau)*(exp(time_width/tau)-1) x1=exp_refl_delay lambda1=0.0184655*2103400 #set title "t=15ns;n=104; Forward OD=4.3" set yrange [0:1500] set key right set xlabel 't {/Symbol G}' font 'Helvetica, 20' offset 0,-1 set ylabel 'N_{d} per 1ns timebin' font 'Helvetica, 20' offset -2,0 set xtics font 'Helvetica, 20' offset 0,-0.5 set ytics font 'Helvetica, 20' offset -0.5,0 plot exp_p_15ns_104_filename u (column(time)-exp_refl_delay)/tau:at_refl_count:(sqrt(column(at_refl_count))) w yerrorl t 'APD 1' ls 2,\ exp_p_15ns_104_filename u (column(time)-exp_trigger_delay)/tau:bg_tran_count:(sqrt(column(bg_tran_count))) w yerrorl t 'APD 2' ls 1,\ f((x*tau+x0)) t 'rising exponential fit' ls 3,\ g((x*tau+x1)) t 'falling exponential fit' ls 4 unset output set term wxt #### Example excitation pulse #### set terminal postscript enhanced color eps set xrange [-5:2] set yrange [0:1600] unset title lambda=0.00381085*14750000 set arrow 1 from graph 0.714,0 to graph 0.714,1 ls 0 lw 2 nohead set output '../pulses/excitation_exp_pulse_15ns_n100.eps' set label 'ND 2 = 4.3 dB' at -4,1200 font 'Helvetica, 20' plot exp_p_15ns_104_filename u (column(time)-exp_trigger_delay)/tau:bg_tran_count:(sqrt(column(bg_tran_count))) w yerrorl t '' ls 1,\ f((x*tau+x0)) t '' ls 4 unset label set output '../pulses/excitation_sq_pulse_15ns_n100.eps' set xrange [-2:2] set yrange [0:*] set label 'ND 2 = 4.0 dB' at -1.5,2500 font 'Helvetica, 20' set arrow 1 from graph 0.5,0 to graph 0.5,1 ls 0 lw 2 nohead plot sq_p_15ns_110_filename u (column(time)-sq_trigger_delay)/tau:bg_tran_count:(sqrt(column(bg_tran_count))) w yerrorl t '' ls 2 set output '../pulses/excitation_exp_pulse_5ns_n100.eps' tau0=5.1938 lambda=0.000437931*6750000 set xrange [-5:2] set yrange [0:250] plot exp_p_5ns_123_filename u (column(time)-exp_trigger_delay)/tau:bg_tran_count:(sqrt(column(bg_tran_count))) w yerrorl t '' ls 1,\ f((x*tau+x0)) t '' ls 4 set output '../pulses/excitation_sq_pulse_5ns_n100.eps' set yrange [0:*] plot sq_p_5ns_117_filename u (column(time)-sq_trigger_delay)/tau:bg_tran_count:(sqrt(column(bg_tran_count))) w yerrorl t '' ls 2 unset output set output '../pulses/excitation_pulse_5ns_n100.eps' tau0=5.1938 lambda=0.000437931*6750000 set xrange [-2:2] set yrange [0:*] plot exp_p_5ns_123_filename u (column(time)-exp_trigger_delay)/tau:bg_tran_prob:(column(bg_tran_prob_err)) axes x1y2 w yerrorl t '' ls 2,\ sq_p_5ns_117_filename u (column(time)-sq_trigger_delay)/tau:bg_tran_prob:(column(bg_tran_prob_err)) axes x1y1 w yerrorl t '' ls 1 # f((x*tau+x0)) t '' ls 4,\ unset output set term wxt