* KUMAC file to create loe plots from loe.prt output from nuscan.f angle 30 -120 opt eah opt date opt fit opt linx opt liny opt linz opt htit set gsiz 0.4 set asiz 0.2 set tsiz 0.3 tit 'L/E study: 1489D dat, 100yr on-mc, 0.2-10 GeV' ve/cr elbin(10) r ve/cr fdel(10) r ve/cr fdmu(10) r ve/cr fmel(10) r ve/cr fmmu(10) r ve/cr ffdel(10) r ve/cr ffdmu(10) r ve/cr ffmel(10) r ve/cr ffmmu(10) r ve/cr rel(10) r ve/cr rmu(10) r ve/cr frel(10) r ve/cr frmu(10) r ve/cr ffrel(10) r ve/cr ffrmu(10) r ve/cr R(10) r ve/cr Rmc(10) r * and their errors ve/cr elbiner(10) r ve/cr fdeler(10) r ve/cr fdmuer(10) r ve/cr fmeler(10) r ve/cr fmmuer(10) r ve/cr ffdeler(10) r ve/cr ffdmuer(10) r ve/cr ffmeler(10) r ve/cr ffmmuer(10) r ve/cr reler(10) r ve/cr rmuer(10) r ve/cr freler(10) r ve/cr frmuer(10) r ve/cr ffreler(10) r ve/cr ffrmuer(10) r ve/cr Rer(10) r ve/cr Rmcer(10) r ve/read elbin,fdel,fdmu,fmel,fmmu,ffdel,ffdmu,ffmel,ffmmu,rel,rmu,frel,frmu loe.prt 13e10.2 ve/read elbiner,fdeler,fdmuer,fmeler,fmmuer,ffdeler,ffdmuer,ffmeler,ffmmuer,reler,rmuer,freler,frmuer loer.prt 13e10.2 1dh 101 'N dat electrons vs log10(L/E/km/GeV)' 10 0.0 5.0 0. 1dh 102 ' N dat muons vs log10(L/E/km/GeV)' 10 0.0 5.0 0. 1dh 103 ' N mc electrons vs log10(L/E/km/GeV)' 10 0.0 5.0 0. 1dh 104 ' N mc muons vs log10(L/E/km/GeV)' 10 0.0 5.0 0. 1dh 105 'N mc osc-electr vs log10(L/E/km/GeV)' 10 0.0 5.0 0. 1dh 106 ' N mc osc-muons vs log10(L/E/km/GeV)' 10 0.0 5.0 0. 1dh 107 ' N mc electr vs log10(L/E/km/GeV)' 10 0.0 5.0 0. 1dh 108 ' N mc muons vs log10(L/E/km/GeV)' 10 0.0 5.0 0. 1dh 110 ' Data/MC Es vs log10(L/E/km/GeV)' 10 0.0 5.0 0. 1dh 111 ' Data/MC Mus vs log10(L/E/km/GeV)' 10 0.0 5.0 0. 1dh 112 ' Osc-MC/MC Es vs log10(L/E/km/GeV)' 10 0.0 5.0 0. 1dh 113 ' Osc-MC/MC Mus vs log10(L/E/km/GeV)' 10 0.0 5.0 0. 1dh 114 ' Dat/Osc-MC Es vs log10(L/E/km/GeV)' 10 0.0 5.0 0. 1dh 115 ' Dat/Osc-MC Mus vs log10(L/E/km/GeV)' 10 0.0 5.0 0. 1dh 116 ' R vs log10(L/E/km/GeV)' 10 0.0 5.0 0. 1dh 117 ' MC R with Osc vs log10(L/E/km/GeV)' 10 0.0 5.0 0. hi/put_vect/contents 101 fdel hi/put_vect/contents 102 fdmu hi/put_vect/contents 103 fmel hi/put_vect/contents 104 fmmu hi/put_vect/contents 105 ffdel hi/put_vect/contents 106 ffdmu hi/put_vect/contents 107 ffmel hi/put_vect/contents 108 ffmmu hi/put_vect/contents 110 rel hi/put_vect/contents 111 rmu hi/put_vect/contents 112 frel hi/put_vect/contents 113 frmu hi/put_vect/errors 101 fdeler hi/put_vect/errors 102 fdmuer hi/put_vect/errors 103 fmeler hi/put_vect/errors 104 fmmuer hi/put_vect/errors 105 ffdeler hi/put_vect/errors 106 ffdmuer hi/put_vect/errors 107 ffmeler hi/put_vect/errors 108 ffmmuer hi/put_vect/errors 110 reler hi/put_vect/errors 111 rmuer hi/put_vect/errors 112 freler hi/put_vect/errors 113 frmuer * fit for s2th: hi/fit 111(8:10) p0 * make ratio of data to oscillated signals, see nuscan for dm2 test value sigma ffrel=fdel/ffdel sigma ffrmu=fdmu/ffdmu hi/put_vect/contents 114 ffrel hi/put_vect/contents 115 ffrmu * this is a overestimate of the error because mc is scaled down sigma ffreler=ffrel*sqrt(1/fdel+1/ffdel) sigma ffrmuer=ffrmu*sqrt(1/fdmu+1/ffdmu) hi/put_vect/errors 114 ffreler hi/put_vect/errors 115 ffrmuer *hi/set/max 101 100. *hi/set/min 101 0.0 * make big R plot sigma R=(fdmu/fdel)/(fmmu/fmel) sigma Rer = R*sqrt((reler/rel)**2+(rmuer/rmu)**2) hi/put_vect/contents 116 R hi/put_vect/errors 116 Rer * make big R plot for mc oscillated data sigma Rmc=(ffdmu/ffdel)/(fmmu/fmel) sigma Rmcer = Rmc*sqrt((freler/frel)**2+(frmuer/frmu)**2) hi/put_vect/contents 117 Rmc hi/put_vect/errors 117 Rmcer zon 1 2 opt liny hi/pl 101 E1 hi/pl 103 S E1 hi/set/max 110 2.0 hi/pl 110 E1 *pict/print 'loe_electrons.ps' *pict/print 'loe_electrons.gif' hi/pl 102 E1 hi/pl 104 S E1 hi/set/max 111 2.0 hi/pl 111 E1 *pict/print 'loe_muons.ps' *pict/print 'loe_muons.gif' hi/pl 105 E1 hi/pl 107 S E1 hi/set/max 112 2.0 hi/pl 112 E1 *pict/print 'loe_mc_electrons.ps' *pict/print 'loe_mc_electrons.gif' hi/pl 106 E1 hi/pl 108 S E1 hi/set/max 113 2.0 hi/pl 113 E1 *pict/print 'loe_mc_muons.ps' *pict/print 'loe_mc_muons.gif' hi/set/max 114 2.0 hi/set/max 115 2.0 hi/pl 114 E1 hi/pl 115 E1 *pict/print 'loe_datomc.ps' *pict/print 'loe_datomc.gif' zon 1 1 hi/set/min 116 0.0 hi/set/max 116 1.2 set mtyp 20 hi/pl 116 E1 set mtype 1 lin 0 1 5 1 set mtyp 0 hi/pl 117 s' ' set mtyp 1 pict/print 'r_vs_loe.ps' *pict/print 'r_vs_loe.gif' zon 2 2 set mtyp 20 hi/set/max 101 900 hi/pl 101 E1 set mtyp 1 set htyp 344 hi/pl 103 s set htype 000 set mtyp 20 hi/set/max 102 900 hi/pl 102 E1 set mtyp 1 set htyp 344 hi/pl 104 s set htyp 000 set mtyp 20 hi/set/max 110 2.0 hi/set/min 110 0.0 hi/pl 110 E1 set mtyp 1 hi/pl 112 sc set ltyp 3 line 0 1 5 1 hi/set/max 111 2.0 hi/set/min 111 0.0 hi/pl 111 E1 set mtyp 1 hi/pl 113 sc set ltyp 3 line 0 1 5 1 pic/print loe.ps *pic/print loe.gif *pic/print loe.eps