* projdata3.kumac * plot projections of bands of DE vs Mes plots for JpL and Jpp * exec pub-plot projdata3 1 1 * clean up everything first close 0 h/del 0 vec/del * * need to adjust some parameters set ylab 1.1 | distance from axis to title set xlab 1.2 set asiz 0.4 set chhe 0.35 set mtyp 20 lwidth = 4 ltype = 2 set mscf 0.7 maxy = 3.4 laby = 2.1 * Start with JpL * 2D data for proper scatter plots appl data JpL-mes-DE.dat 5.27595 -0.01113 1.00000 -1 5.28308 0.12180 1.00000 -1 5.27967 -0.00832 1.00000 1 5.23511 0.15526 1.00000 1 5.22422 0.10544 1.00000 1 5.24035 0.22405 1.00000 1 5.23631 0.20489 1.00000 1 5.28168 -0.00340 1.00000 1 5.25914 0.22836 1.00000 -1 5.23965 -0.01312 1.00000 -1 5.28176 0.02333 1.00000 1 5.23941 0.22772 1.00000 -1 5.25992 0.15003 1.00000 -1 5.26841 0.15068 1.00000 1 5.28606 0.17209 1.00000 1 5.24736 0.12454 1.00000 -1 5.22898 0.21024 1.00000 1 5.26806 0.21374 1.00000 1 5.28502 -0.00485 1.00000 1 5.25122 0.18924 1.00000 1 5.26146 0.24032 1.00000 1 5.21919 0.21911 1.00000 -1 5.21812 0.19503 1.00000 1 5.24124 0.11708 1.00000 -1 5.24261 0.19053 1.00000 1 5.21533 0.21389 1.00000 -1 5.21611 0.18286 1.00000 1 5.25107 0.10642 1.00000 1 5.20558 0.06228 1.00000 -1 5.26451 0.18244 1.00000 -1 5.22778 0.04929 1.00000 1 5.20567 0.19737 1.00000 1 5.25792 0.14076 1.00000 -1 5.23025 0.20760 1.00000 -1 5.28862 0.13164 1.00000 -1 5.21962 0.15566 1.00000 -1 5.24414 0.20279 1.00000 -1 5.25925 0.16536 1.00000 -1 5.25951 0.08006 1.00000 1 5.20172 0.12150 1.00000 -1 5.22145 0.13269 1.00000 -1 5.21023 0.13667 1.00000 -1 5.26804 0.14518 1.00000 -1 JpL-mes-DE.dat nt/create 1 'jpl' 4 ! ! jplmes jplde jplwt jplcharge nt/read 1 JpL-mes-DE.dat id = 100 2d [id] 'JpL De vs Mes' 50 5.2 5.3 70 -0.10 0.25 * uncomment to get points with errors on projections * H/SET/IDOPT [id] error * size of ellipse mNomChB = 5.279 resMChB = 0.0031 resDEChB = 0.0065 SigEllChB = 2.4 minDEChB = 0 - [SigEllChB]*[resDEChB] maxDEChB = 0 + [SigEllChB]*[resDEChB] minMChB = [mNomChB] - [SigEllChB]*[resMChB] maxMChB = [mNomChB] + [SigEllChB]*[resMChB] * Scatter plot zone 2 2 set ndvy 404 nt/plot 1.jplde%jplmes IDH=[id] atitle 'M?ES! (GeV/c^2!)' '[D]E (GeV)' itx 5.202 -0.08 (a) itx 5.202 0.26 B^+! "5# J/[y]p[L]&^-! set lwid [lwidth]; set ltyp [ltype] dline [minMChB] [minMChB] -0.1 0.25 dline [maxMChB] [maxMChB] -0.1 0.25 dline 5.2 5.3 [minDEChB] [minDEChB] dline 5.2 5.3 [maxDEChB] [maxDEChB] rx = [SigEllChB]*[resMChB] ry = [SigEllChB]*[resDEChB] ellipse [mNomChB] 0. [rx] [ry] h/create/banx [id] [minDEChB] [maxDEChB] h/create/bany [id] [minMChB] [maxMChB] h/project [id] * Mes projection zone 2 4 5 s max [id] [maxy] set ndvy 204 h/plot [id].banx.1 atit 'M?ES! (GeV/c^2!)' 'Entries' itx 5.202 [laby] (b) set lwid [lwidth]; set ltyp [ltype] dline [minMChB] [minMChB] 0. [maxy] dline [maxMChB] [maxMChB] 0. [maxy] zone 2 4 7 s set ndvy 204 h/plot [id].bany.1 atit '[D]E (GeV)' 'Entries' itx -0.095 [laby] (c) set lwid [lwidth]; set ltyp [ltype] dline [minDEChB] [minDEChB] 0. [maxy] dline [maxDEChB] [maxDEChB] 0. [maxy] * Now Jpp appl data Jpp-mes-DE.dat 5.27864 0.14614 1.00000 5.21567 -0.11749 1.00000 5.24644 0.17027 1.00000 5.28286 -0.18816 1.00000 5.21736 0.09066 1.00000 5.25116 0.14129 1.00000 5.24421 0.09860 1.00000 5.26241 0.17507 1.00000 5.25294 0.14676 1.00000 5.25220 -0.02851 1.00000 5.20750 -0.05198 1.00000 5.21918 0.23580 1.00000 5.24213 -0.07994 1.00000 5.21099 0.17125 1.00000 5.25764 0.16213 1.00000 5.28135 -0.16236 1.00000 5.25764 0.08894 1.00000 5.21539 -0.03269 1.00000 5.22415 -0.06578 1.00000 5.22338 0.18278 1.00000 5.21092 0.17114 1.00000 5.23332 -0.16550 1.00000 5.27726 -0.15465 1.00000 5.21748 -0.01192 1.00000 5.24267 0.24815 1.00000 5.20018 -0.02546 1.00000 5.21929 0.17788 1.00000 5.26376 0.14289 1.00000 5.24065 0.01973 1.00000 5.21605 -0.03481 1.00000 5.22530 0.05302 1.00000 5.20919 -0.08260 1.00000 5.26617 -0.01656 1.00000 5.27019 -0.09655 1.00000 5.21690 0.06283 1.00000 5.25879 0.07412 1.00000 5.25050 0.17992 1.00000 5.28276 -0.00022 1.00000 5.27157 0.07621 1.00000 5.26836 -0.11628 1.00000 5.23278 0.14616 1.00000 5.26551 -0.02032 1.00000 5.22000 -0.04623 1.00000 5.20774 0.11266 1.00000 5.25644 -0.13987 1.00000 5.27900 -0.03252 1.00000 5.23998 -0.05414 1.00000 5.20675 -0.14672 1.00000 5.27266 -0.12808 1.00000 5.25119 0.17851 1.00000 5.23961 0.17021 1.00000 5.21345 0.23810 1.00000 5.27987 0.18247 1.00000 5.26104 0.19500 1.00000 5.27856 0.11421 1.00000 5.22259 0.11555 1.00000 5.24438 0.14795 1.00000 5.28090 0.10547 1.00000 5.25358 0.12119 1.00000 5.21355 -0.07831 1.00000 5.24212 0.04498 1.00000 5.20537 0.10839 1.00000 5.25383 0.24718 1.00000 5.22532 -0.01647 1.00000 5.27864 -0.02632 1.00000 5.20058 0.11215 1.00000 5.21267 -0.00018 1.00000 5.27121 0.08367 1.00000 5.22937 -0.11924 1.00000 5.23996 -0.00839 1.00000 5.25507 -0.08667 1.00000 5.25007 -0.02519 1.00000 5.23814 0.17878 1.00000 5.25599 -0.10151 1.00000 5.23527 0.10584 1.00000 5.20317 0.07283 1.00000 5.22976 -0.12875 1.00000 5.27086 0.14588 1.00000 5.23584 -0.06446 1.00000 5.25819 -0.09205 1.00000 5.23574 0.11132 1.00000 5.23275 0.07683 1.00000 5.21554 0.07653 1.00000 5.28753 0.01353 1.00000 5.20593 0.00423 1.00000 5.27353 0.07799 1.00000 5.24978 -0.18366 1.00000 5.21003 -0.11930 1.00000 5.20197 -0.02106 1.00000 5.20593 0.06307 1.00000 5.20036 -0.03915 1.00000 5.20239 0.18693 1.00000 5.26180 0.07595 1.00000 5.25454 0.00795 1.00000 5.25904 -0.06072 1.00000 5.23797 0.14815 1.00000 5.20009 -0.03960 1.00000 5.20945 -0.01432 1.00000 5.21056 0.14672 1.00000 5.27529 0.10689 1.00000 5.27991 -0.06214 1.00000 5.26362 0.05441 1.00000 5.23118 -0.01877 1.00000 5.23895 0.15394 1.00000 5.26550 -0.04378 1.00000 5.20777 0.10267 1.00000 5.21110 0.15937 1.00000 5.27765 -0.11391 1.00000 5.23499 0.07516 1.00000 5.28503 -0.16189 1.00000 5.28467 0.09286 1.00000 5.22252 0.17977 1.00000 5.25921 0.13369 1.00000 5.20247 0.21435 1.00000 5.25572 -0.09366 1.00000 5.27710 -0.04939 1.00000 5.27431 0.19721 1.00000 5.21837 0.22309 1.00000 5.20427 0.07476 1.00000 5.26196 0.04970 1.00000 5.20105 0.09329 1.00000 5.24752 -0.15122 1.00000 5.24982 0.00323 1.00000 5.23725 0.14668 1.00000 5.26533 -0.03735 1.00000 5.23553 0.07501 1.00000 5.20361 0.12104 1.00000 Jpp-mes-DE.dat nt/create 2 'jpl' 3 ! ! jppmes jppde jppwt nt/read 2 Jpp-mes-DE.dat id = 101 2d [id] 'Jpp De vs Mes' 50 5.2 5.3 50 -0.25 0.25 * uncomment to get points with errors on projections * H/SET/IDOPT [id] error * size of ellipse mNomB0 = 5.279 resMB0 = 0.0027 resDEB0 = 0.0055 SigEllB0 = 2.2 minDEB0 = 0 - [SigEllB0]*[resDEB0] maxDEB0 = 0 + [SigEllB0]*[resDEB0] minMB0 = [mNomB0] - [SigEllB0]*[resMB0] maxMB0 = [mNomB0] + [SigEllB0]*[resMB0] * Scatter plot zone 2 2 2 s set ndvy 404 nt/plot 2.jppde%jppmes IDH=[id] atitle 'M?ES! (GeV/c^2!)' '[D]E (GeV)' itx 5.204 -0.2214 (d) itx 5.204 0.2643 B^0! "5# J/[y]pp&^_! set lwid [lwidth]; set ltyp [ltype] dline [minMB0] [minMB0] -0.25 0.25 dline [maxMB0] [maxMB0] -0.25 0.25 dline 5.2 5.3 [minDEB0] [minDEB0] dline 5.2 5.3 [maxDEB0] [maxDEB0] rx = [SigEllB0]*[resMB0] ry = [SigEllB0]*[resDEB0] ellipse [mNomB0] 0. [rx] [ry] h/create/banx [id] [minDEB0] [maxDEB0] h/create/bany [id] [minMB0] [maxMB0] h/project [id] * Mes projection zone 2 4 6 s max [id] [maxy] set ndvy 204 h/plot [id].banx.1 atit 'M?ES! (GeV/c^2!)' 'Entries' itx 5.204 [laby] (e) set lwid [lwidth]; set ltyp [ltype] dline [minMB0] [minMB0] 0. [maxy] dline [maxMB0] [maxMB0] 0. [maxy] zone 2 4 8 s set ndvy 204 h/plot [id].bany.1 atit '[D]E (GeV)' 'Entries' itx -.23 [laby] (f) set lwid [lwidth]; set ltyp [ltype] dline [minDEB0] [minDEB0] 0. [maxy] dline [maxDEB0] [maxDEB0] 0. [maxy] zone exec makeps projdata3.ps sh ps2epsi projdata3.ps sh eps2eps projdata3.epsi projdata3.eps sh /bin/rm projdata3.ps sh /bin/rm projdata3.epsi