macro limit col=2
 
 h/del *
 v/del *
 
 chain -lim
 chain lim 188.6+183_comb_l3_all_80.0-90.0_syst.ntp 188.6+183_comb_l3_all_90.5-100.0_syst.ntp 188.6+183_comb_l3_all_93.1-96.9_syst.ntp
 chain lim -P ~/temp/hist/
 
 cd //lim
 ntuple/loop 1 x.f(1,mh,type.eq.4)
 ntuple/loop 1 x.f(2,(min(max(prb,0.0),1.0)),type.eq.5)
 
 sigma mh=order(v1,v1)
 sigma prb=order(v2,v1)
 
* Make the plot
 exec style_1by2
 opt utit
 opt liny
 
* Luck plot
 set xlab 1.8
 set xval 0.65
 tit 'm?H! (GeV)' u
 null 80 97 0.0001 1.01
 set mscf 0.00001
 set dmod 1
 set plci [col]
 v/pl prb%mh ! sl
 atit ! 'CL?b!'
 opt liny
 null 0 1 0 1 sab
 set csiz 0.8
 key 0.15 0.2 1 'CL?b! for 189 GeV data' ! L