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 0.5 pt 7 set style line 2 lc 3 lw 1 ps 0.5 pt 7 set style line 3 lc 0 lw 2 ps 1.5 set style line 4 lc 0 lw 1 ps 1.5 lt 1 #### 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 [-100:100] #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 size 8.6cm,6cm font "Helvetica,16pt" set output '../pulses/excitation_prob_exc_pulse_and_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 x1=exp_refl_delay tau0=15.3907 eta=0.3 lambda=0.00381085*14750000 g(x)= eta*lambda1*exp((x1-x)/tau)*(exp(time_width/tau)-1) lambda1=0.0184655*2103400 overlap=0.03 p(x)=x/(overlap*eta*(time_width/tau)) # to normalize reflection counts into probabilities coeff=p(4.754E-7) set multiplot set format y "%4.0f" # bottom plot is the fluorescence data set origin -0.025,0 set size 1.05,0.55 #set title "t=15ns;n=104; Forward OD=4.3" set yrange [0:600] set ytics nomirror set y2range [0:coeff*600] set y2tics nomirror set y2label offset -1,0 "P_e(t)" set label 1 at -90,500 "(b)" set xlabel 't (ns)' # set ylabel 'N_{d} per 1ns timebin' set ylabel offset 1,0 'Events' plot exp_p_15ns_104_filename u \ (column(time)-exp_refl_delay):at_refl_count:(sqrt(column(at_refl_count)))\ w yerrorbars notitle ls 2,\ g((x+x1)) notitle ls 4 #t 'falling exponential fit' ls 4 # top plot excitation pulse data set origin -0.025,0.5 set size 1.05,0.5 set xlabel "" set format x "" set y2label " " set format y2 " " set y2tics 100 set label 1 at -90,1200 "(a)" set yrange [0:1500] plot exp_p_15ns_104_filename u \ (column(time)-exp_trigger_delay):bg_tran_count:(sqrt(column(bg_tran_count))) \ w yerrorbars ls 1 notitle,\ f((x+x0)) notitle ls 4 #t 'rising exponential fit' ls 4