* KUMAC file to create exclusion plots from excl.prt output from nuscan1.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 'inclusion plot: 1489D dat, 100yr onsite mc, 0.2-10 GeV' ve/cr chisq(20,100) r ve/read chisq excl.prt (20e11.2) * change range to agree with others going from 1e-5 <-> 1.0 *2dh 5000 'Chisq for dm2 vs S2theta' 20 0.0 1.0 100 -4.0 -1.0 0.0 2dh 5000 'Chisq for dm2 vs S2theta' 20 0.0 1.0 100 -5.0 0.0 0.0 hi/put_vect/contents 5000 chisq zon 1 1 opt liny opt linx max 5000 40.00 hi/pl 5000 surf3 opt liny pict/print 'excl.ps' *pict/print 'excl.gif' *pict/print 'excl.eps' * now make Mark Messier style contour, using adaptation of his nmtkamsk.kumac zon 1 1 *hi/file 1 /home/atmpd/tmp/hist.nmt.remap *hi/file 2 /home/atmpd/tmp/hist.kam set *; set txfp 40; set hwid 6; set lwid 6; set bwid 6; set chhe 0.5 title 'SuperK 713 Days L/E Scan' null 0. 1. 0.00001 1.0 a *vec/cre f1(1) R 2.279; vec/cre f2(1) R 4.605; vec/cre f3(1) R 9.210 * our minimum is at 7.09, so add to contours vec/cre f1(1) R 10.12; vec/cre f2(1) R 12.45; vec/cre f3(1) R 18.05 set ltyp 2; contour 5000 1 2S f1 set ltyp 1; contour 5000 1 2S f2 set ltyp 3; contour 5000 1 2S f3 *set ltyp 1; contour //lun2/5004 1 2S f2 null 0 1 0 1 abs set laof 0.021; set lasi 0.0275 axis 0 1 0 0 0 1 205 dnh set laof 0.035 axis 0 0 0 1 0.00001 1 205 dgnh axis 1 1 0 1 0.00001 1 205 +udgnh set txal 20 set chhe 0.6 itx 0.2 0.9 '[n]?[m]!-[n]?[t]!' set chhe 0.5 *itx 0.4 0.70 'Kamiokande'; arrow 0.43 0.63 0.69 0.65 itx 0.4 0.55 'Super-Kamiokande'; arrow 0.5 0.7 0.54 0.5 set csiz 0.6 key 0.15 0.25 2 '' ! l; key 0.15 0.20 1 '' ! l; key 0.15 0.15 3 '' ! l itx 0.22 0.24 '68"Y#'; itx 0.22 0.19 '90"Y#'; itx 0.22 0.14 '99"Y#' itx 0.5 -0.09 'sin^2!2[q]' set tang 90; itx -0.0925 0.5 '[D]m^2! (eV^2!)' pic/prin allowed_loe.ps