* *********************************** * KInematic Plane kumac * * Inherited from Christoph Amelung * * Author: N. Tuning * * Date: 17/04/01 * * * * Usage: * * exec kinplane#read * * exec kinplane#plot * * *********************************** * =========== Macro kinplane * =========== exec read exec h1plot RETURN * =========== Macro read * =========== opt * set * ve/del * opt zfl1 * -------------------- * --- ZEUS 1998/99 --- * -------------------- * -- A.Kappes:-- v/cre z98_x(17) R_ 0.0036 0.0036 0.01 0.03 0.1 0.56 0.99 0.99 0.99 0.99 _ 0.56 0.56 0.32 0.32 0.125 0.12 0.12 v/cre z98_q2(17) R _ 185. 260. 902 2706 9020 50000 50000 2500 2500 2500 _ 2500 900. 900. 520. 400.889 400. 185 * -------------------- * --- ZEUS 1996/97 --- * -------------------- * -- N.Tuning:-- v/cre z967_x(22) R _ 3.1E-05 0.0001 0.0003 0.001 0.003 0.01 0.03 0.1 0.56 _ 0.99 0.926 0.862 0.714 0.556 0.385 0.273 0.2 0.111 _ 0.03 0.01 0.00547 3.1E-05 v/cre z967_q2(22) R _ 2. 6.41513 19.3612 65.9258 210.731 902 2706 9020 50000 _ 50000 5000. 2500. 1000. 500. 250. 150. 100. 50. _ 11.5039 3.68126 2. 2. * -------------------- * --- ZEUS BPT 97 --- * -------------------- * -- Ch.Amelung:-- v/cre zbpt_x(35) R _ 5.275E-07 7E-07 1E-06 2E-06 4E-06 6E-06 8E-06 1E-05 2E-05 4E-05 _ 6E-05 8E-05 0.0001 0.0002 0.0004 0.0006 0.0008 0.001 0.00165 _ 0.00048 0.0002 0.0001 8E-05 6E-05 4E-05 2E-05 1E-05 8E-06 6E-06 _ 4E-06 2E-06 1E-06 8E-07 6E-07 5.275E-07 v/cre zbpt_q2(35) R _ 0.04 0.0530387 0.0757703 0.148248 0.251646 0.327872 0.386394 0.432737 _ 0.569299 0.675956 0.720981 0.74 0.74 0.74 0.74 0.74 0.74 0.74 _ 0.74 0.215703 0.214209 0.211695 0.21046 0.208434 0.204496 0.193527 _ 0.174778 0.166702 0.154783 0.135418 0.0984626 0.0636968 0.0532685 _ 0.04 0.04 * ---------------------- * --- ZEUS SVTX 1995 --- * ---------------------- vec/create z95_q2(9) R 0.45 20.0 20.0 1.08 1.08 0.50 0.50 0.45 0.45 vec/create z95_y(9) R 0.61 0.61 0.05 0.05 0.23 0.23 0.45 0.45 0.61 sigma z95_x=z95_q2/(4*27.52*820*z95_y) * -------------------- * --- H1 Spacal 96-97 --- * -------------------- * -- Ch.Amelung: -- * ve/cr hsp_x(9) r 0.12 2E-3 1E-4 1.5E-5 2.7E-5 4E-4 5E-4 1.1E-3 7E-3 * ve/cr hsp_q2(9) r 120 1.8 1.8 1.0 2. 31 31 70 120 * -- N.Tuning:-- * ve/cr hsp_x(14) r _ * 1.E-4 1.3E-4 4.E-4 0.001 0.003 0.005 0.02 _ * 0.3 0.3 0.01 6.E-3 1.5E-4 0.8E-4 2.3E-5 * ve/cr hsp_q2(14) r _ * 6. 10. 30. 55. 100. 130. 180. _ * 180. 100. 2.5 1.8 1.8 1.5 1.5 * replaced by * published low Q2 96/97 (Doris Eckstein) ve/cr hsp_x(14) r _ 1.E-4 1.23E-4 3.9E-4 1.E-3 1.4E-3 4.E-3 1.6E-2 1.6E-2 _ 0.368 0.368 .00655 1.E-4 1.E-4 2.5E-5 ve/cr hsp_q2(14) r _ 6. 10. 31.6 56. 75. 135. 135. 180. _ 180. 100. 1.78 1.78 1.33 1.33 * -------------------- * --- H1 Spacal 96-97 --- * -------------------- v/crea hm99_x(12) r 1.6E-5 2.2E-5 1.97E-5 6.2E-5 _ 7.E-5 9.3E-5 1.1E-4 1.95E-4 0.13 _ 2.3E-2 8.E-5 4.5E-5 v/crea hm99_q2(12) r 1.33 1.78 1.78 5.6 _ 5.6 7.5 7.5 13.3 13.3 _ 2.33 2.33 1.33 * -------------------- * --- H1 94-97 --- * -------------------- * -- N.Tuning:-- * ve/cr h1_x(12) r _ * 0.0025 0.003 0.007 0.018 0.56 0.99 0.95 0.8 0.5 _ * 0.08 0.015 0.005 * ve/cr h1_q2(12) r _ * 130. 180. 500. 1500. 50000. 50000. 1000. 400. 180. _ * 180. 180. 130. ve/cr h1_x(11) r _ 0.0022 0.0140 0.0098 0.509 0.80 0.80 0.8 0.5 _ 0.08 0.0155 0.005 ve/cr h1_q2(11) r _ 141.25 891.25 891.25 46334. 46334. 1000. 400. 180. _ 180. 141.25 141.25 * H1 high Q2 (Ana Dubak & Burkard Reisert) ve/cr hfl_x(4) r _ 0.000978535 0.00139791 0.0139791 0.00978535 ve/cr hfl_q2(4) r _ 89.125 89.125 891.25 891.25 * ------------------- * --- H1 ISR 1995 --- * ------------------- * H1 ISR v/cre his_x(7) r 0.000022 0.000047 0.00047 0.0032 0.0001 _ 0.000047 0.000047 v/cre his_q2(7) r 0.2 0.2 0.4 5.0 5.0 _ 2.5 0.4 * ------------------- * --- H1 ISR 1997 --- * ------------------- * H1 ISR 97 * updated analysis (Cigdem Issever) v/cre hi97_x(17) r .148E-04 .322E-03 .780E-03 .204E-02 .204E-02 _ .839E-02 .839E-02 .469E-01 .469E-01 .198E-02 .322E-03 .148E-03 _ .686E-04 .686E-04 .322E-04 .148E-04 .148E-04 v/cre hi97_q2(17) r .2830 .2830 .5100 .9870 1.9360 3.5400 6.5200 _ 13.0400 34.6400 21.2500 13.0400 6.5200 3.5400 1.9360 .9870 _ .5100 .2830 * -------------------- * --- H1 SVTX 1995 --- * -------------------- * --Ch.Amelung:-- * vec/create h95_q2(19) R _ * 0.3162 4.217 4.217 3.162 3.162 2.371 2.371 _ * 1.334 1.334 0.99999 0.99999 0.7499 0.7499 0.5623 _ * 0.5623 0.4217 0.4217 0.3162 0.3162 * vec/create h95_x(19) R _ * 4.6E-6 0.0000631 0.00158 0.00158 0.001 0.001 0.000631 _ * 0.000631 0.000398 0.000398 0.0001 0.0001 2.07E-5 1.56E-5 _ * 1.13E-5 0.85E-5 0.85E-5 0.637E-5 4.6E-6 * --N.Tuning:-- vec/create h95_q2(6) R _ 0.3 4. 4. 1. 0.7 0.3 vec/create h95_x(6) R _ 5.E-6 7.E-5 0.0015 2.5E-4 2.E-5 7.E-6 * ------------ * --- SLAC --- * ------------ vec/create SLC_x(8) R _ 7.1969E-2 2.2230E-1 2.4421E-1 4.4984E-1 8.2864E-1 8.2864E-1 _ 7.1969E-2 7.1969E-2 vec/create SLC_q2(8) R _ 6.0209E-1 7.0381E-1 9.1059E-1 9.8451E-1 9.2491E+0 2.9825E+1 _ 2.0991E+0 6.0209E-1 * ------------ * --- E665 --- * ------------ vec/create e665_x(7) R _ 7.1969E-4 1.6768E-3 5.1795E-1 5.1795E-1 1.2649E-1 7.1969E-4 7.1969E-4 vec/create e665_q2(7) R _ 2.1826E-1 2.1826E-1 1.0000E+1 5.3557E+1 5.3557E+1 6.9291E-1 2.1826E-1 * ----------- * --- NMC --- * ----------- vec/create nmc_x(9) R _ 5.6899E-3 1.6000E-2 1.6000E-2 9.5409E-2 5.5000E-1 5.5000E-1 3.0000E-1 _ 9.0000E-3 5.6899E-3 vec/create nmc_q2(9) R _ 4.6538E-1 4.6538E-1 1.0000E-0 1.0000E-0 7.3182E+0 7.0000E+1 7.0000E+1 _ 5.5000E+0 5.5000E+0 * ------------- * --- BCDMS --- * ------------- vec/create bcd_x(16) R _ 7.1968E-2 1.0481E-1 1.5264E-1 1.8421E-1 2.0235E-1 2.6826E-1 _ 3.2374E-1 5.1794E-1 6.2505E-1 7.5431E-1 7.5431E-1 5.1794E-1 _ 3.2374E-1 1.6768E-1 7.1968E-2 7.1968E-2 vec/create bcd_q2(16) R _ 7.0381E+0 7.0381E+0 9.2491E+0 9.2491E+0 1.0000E+1 1.0000E+1 _ 1.0981E+1 1.0981E+1 1.1873E+1 2.3598E+1 2.4537E+2 2.4537E+2 _ 1.1242E+2 5.1506E+1 2.0187E+1 7.0381E+0 * ---------------- * --- JLAB --- * ---------------- vec/create jlab_x(4) R 0.72 0.42 0.09 0.2 vec/create jlab_q2(4) R 3.3 3.3 0.45 0.45 * ------------ * --- jinr/ihep --- * ------------ vec/create jinr_x(6) R 8.E-2 2.5E-1 6.5E-1 6.5E-1 4.5E-1 8.0E-2 vec/create jinr_q2(6) R 2.E-1 2.0E-1 7.0E-1 1.0E+1 1.0E+1 3.0E+0 * -------------- * --- CHORUS --- * -------------- vec/create chor_x(6) R 1.E-2 1.E-2 5.2E-2 7.E-1 7.E-1 5.E-1 vec/create chor_q2(6) R 2.E+0 4.E-1 4.0E-1 6.E+0 1.E+2 1.E+2 * ------------ * --- ccfr --- * ------------ vec/create ccf_x(6) R 5.E-3 5.E-3 5.2E-2 7.E-1 7.E-1 3.8E-1 vec/create ccf_q2(6) R 3.E+0 1.E+0 1.E+0 1.7E+1 2.7E+2 2.7E+2 * ---------------------------- * --- FIXED TRAGET COMBINED: --- * ---------------------------- vec/create ft_x(14) R 8.0E-2 1.5E-2 1.5E-2 5.0E-3 1.5E-3 7.0E-4 _ 7.0E-4 1.0E-2 4.0E-1 8.0E-1 6.5E-1 2.5E-1 8.0E-2 8.0E-2 vec/create ft_q2(14) R 1.0E+0 1.0E+0 5.0E-1 5.0E-1 2.1E-1 2.1E-1 _ 7.0E-1 6.0E+0 2.7E+2 2.7E+2 6.0E-1 2.0E-1 2.0E-1 1.0E+0 * ---------------- * --- jets cdf --- * ---------------- * From Christoph Amelung: vec/create cdf_x(5) R 9.E-2 2.7E-1 7.E-1 7.5E-1 5.7E-1 vec/create cdf_q2(5) R 3.E+3 1.8E+3 4.E+4 1.0E+5 2.3E+5 * From CTEQ5D (hep-ph/9903282) Jets DY DY Dir.ph Dir.ph W.as Jets vec/create cdf2_x(7) R 0.5 0.08 0.008 0.0133 0.05 0.3 0.5 vec/create cdf2_q2(7) R 122500. 22500. 150. 150. 200. 6400. 122500. * ---------------- * --- jets LHC --- (see MRST, hep-ph 9907231) * ---------------- vec/create lhc_x(5) R 5.E-7 4.E-1 0.99E+0 0.99E+0 5.E-7 vec/create lhc_q2(5) R 1.E+2 1.E+8 1.E+8 1.E+2 1.E+2 * ---------------- * --- THERA --- (see http://www-zeuthen.desy.de/thera/) * ---------------- vec/create th_x(5) R 1.E-6 1.E+0 1.E+0 1.E-5 1.E-6 vec/create th_q2(5) R 1. 1.E+6 1.E+5 1. 1. Return * ================================================================= *================================================= Macro plot z98=0 z967=1 zbpt=1 zsvt=1 _ h1=1 nmc=1 bcdms=1 ccfr=0 e665=1 slac=1 other=0 *================================================= set xval 0.45 set xlab 1.3 set fais 0 set faci 1 set fasi 000 set txfp -20 set *fon -20 set asiz 0.5 set vsiz 0.4 * --- opt ndate opt logx opt logy * --- KN. set xmgl 2.3 opt nbox set xval 0.55 set ylab 1.25 set asiz 0.55 set vsiz 0.45 set xlab 1.5 * ========================== null 4.e-7 1.001 0.035 90200. atit 'x' 'Q^2! (GeV^2!)' igset chhe 0.35 igset tang 45. itx 0.00003 3. 'Kinematic limit y=1' itx 0.00023 0.045 'y=0.004' *========================== * ------- H1: --------- If ([h1]=1) Then * set fais 3 * set fasi 145 set fais 1 set faci 5 set plci 5 set lwid 5. exec plot2 hsp exec plot2 h1 set plci 3 set faci 3 exec plot2 h95 set faci 1 set lwid 1. Endif * ------- ZEUS: --------- If ([z967]=1) Then * set fais 1 * set faci 2 * exec plot2 z967 set fais 3 set fasi 145 set faci 2 set plci 2 set lwid 5. exec plot2 z967 set plci 1 set lwid 1. Endif If ([z98]=1) Then set fais 3 set fasi 200 set faci 4 set plci 4 set lwid 5. exec plot2 z98 set plci 1 set lwid 1. Endif If ([zbpt]=1) Then set fais 1 set fais 3 set fasi 154 set faci 7 set plci 7 set lwid 5. exec plot2 zbpt set plci 1 set lwid 1. Endif If ([zsvt]=1) Then set fais 1 set fais 3 set fasi 144 set faci 4 set plci 4 set lwid 5. exec plot2 z95 set plci 1 set lwid 1. Endif * ============== set fais 0 set faci 1 set fasi 000 * ============== * ----- H1:-------- If ([h1]=1) Then set fais 2 set lwid 5 set plci 5 * exec plot2 hsp * exec plot2 h1 set plci 3 * exec plot2 h95 set lwid 1 set faci 1 set plci 1 Endif * ------------ set faci 3 * ------------ * --- CCFR: ---- If ([ccfr]=1) Then set fais 3 set fasi 350 exec plot2 ccf Endif * ---E665:--- If ([e665]=1) Then set fais 3 set fasi 359 exec plot2 e665 Endif * ---SLAC:---- If ([slac]=1) Then set fais 3 set fasi 159 exec plot2 slc Endif If ([other]=1) Then set faci 1 set faci 0 set fais 0 * ---CHORUS:--- exec plot2 chor * --- JLAB & JINR: ---- exec plot2 jlab exec plot2 jinr Endif * ------------- set faci 1 set lwid 5 * ------------- * ---NMC:--- If ([nmc]=1) Then set fais 3 set fasi 345 exec plot2 nmc Endif * --- BCDMS:----- If ([bcdms]=1) Then set fais 3 set fasi 354 exec plot2 bcd Endif set lwid 1 * ====================================================== set dmod 1 line -7. -2.0448 0. 4.9552 | y=1 line -5. -2.4427335 0. 2.5572665 | y=0.004 * ====================================================== Exec Label 0.1 0.93 [z98] [z967] [zbpt] [zsvt] [h1] _ [nmc] [bcdms] [ccfr] [e665] [slac] [other] * ====================================================== pi/pri kinplane.ps Return * ===================== Macro plot2 wat=z96 * ===================== set bord 1 sigma vecx=[wat]_x sigma vecq=[wat]_q2 graph $vdim(vecx) vecx vecq f Return * ----------------------------------------------------------------- *=============================================================== Macro label xs=0.15 ys=0.93 _ z98=1 z967=1 zbpt=1 zsvt=0 _ h1=1 nmc=1 bcdms=1 ccfr=1 e665=1 slac=1 other=0 * =============================================================== set cfon -20 set csiz 0.3 set dmod 1 * ------------------------------------ xmin=$grafinfo('wnxmin') xmax=$grafinfo('wnxmax') x=[xmin]+[xs]*([xmax]-[xmin]) maxw=$grafinfo('wnymax') y1=$sigma(10.**([ys]*[maxw])) y2=$sigma(10.**(([ys]-0.07)*[maxw])) y3=$sigma(10.**(([ys]-0.14)*[maxw])) y4=$sigma(10.**(([ys]-0.21)*[maxw])) y5=$sigma(10.**(([ys]-0.28)*[maxw])) y6=$sigma(10.**(([ys]-0.36)*[maxw])) y7=$sigma(10.**(([ys]-0.43)*[maxw])) y8=$sigma(10.**(([ys]-0.50)*[maxw])) y9=$sigma(10.**(([ys]-0.57)*[maxw])) y10=$sigma(10.**(([ys]-0.64)*[maxw])) y11=$sigma(10.**(([ys]-0.71)*[maxw])) x=$sigma(10.**([x])) * ------------------------------------------- i=0 y=[y[i]] * ---ZEUS:--- If ([z967]=1) Then i=[i]+1 y=[y[i]] set fais 3 set fasi 154 set faci 2 key [x] [y] 145 'ZEUS 1996-97' 0.8 H EndIf If ([z98]=1) Then i=[i]+1 y=[y[i]] set fais 3 set fasi 200 set faci 4 key [x] [y] 200 'ZEUS 1998-99 (Preliminary)' 0.8 H EndIf If ([zbpt]=1) Then i=[i]+1 y=[y[i]] set fais 3 set fasi 154 set faci 7 key [x] [y] 154 'ZEUS BPT 1997' 0.8 H EndIf If ([zsvt]=1) Then i=[i]+1 y=[y[i]] set fais 3 set fasi 144 set faci 4 key [x] [y] 144 'ZEUS SVX 1995' 0.8 H EndIf * ============== set fais 0 set faci 1 set fasi 000 * ============== * ---H1:--- If ([h1]=1) Then i=[i]+1 y=[y[i]] * set fais 3 * set fasi 145 set fais 1 set faci 5 * key [x] [y] 145 'H1 1994-97' 0.8 H key [x] [y] 5 'H1 1994-97' 0.8 F i=[i]+1 y=[y[i]] * set faci 3 set faci 3 * key [x] [y] 145 'H1 SVX 1995' 0.8 H key [x] [y] 3 'H1 SVX 1995' 0.8 F set faci 1 EndIf * ------------ If ([nmc]=1) Then set fais 3 set fasi 345 i=[i]+1 y=[y[i]] key [x] [y] 345 'NMC' 0.8 H EndIf * ------------ If ([bcdms]=1) Then set fais 3 set fasi 354 i=[i]+1 y=[y[i]] key [x] [y] 354 'BCDMS' 0.8 H Endif * ------------ If ([ccfr]=1) Then set faci 3 set fais 3 set fasi 350 i=[i]+1 y=[y[i]] key [x] [y] 350 'CCFR' 0.8 H EndIf * ------------ If ([e665]=1) Then set faci 3 set fais 3 set fasi 359 i=[i]+1 y=[y[i]] key [x] [y] 359 'E665' 0.8 H EndIf * ------------ If ([slac]=1) Then set faci 3 i=[i]+1 y=[y[i]] key [x] [y] 159 'SLAC' 0.8 H EndIf * ----Other Fixed Target (at high x):----- If ([other]=1) Then set faci 0 set fais 0 i=[i]+1 y=[y[i]] key [x] [y] 1 'CHORUS+JLAB+JINR' 0.8 H EndIf Return macro h1plot set xval 0.45 set xlab 1.3 set fais 0 set faci 1 set fasi 000 set txfp -20 set *fon -20 set asiz 0.5 set vsiz 0.4 set xmgr 0.2 set ymgu 0.2 * --- opt ndate opt logx opt logy * --- KN. set xmgl 2.3 opt nbox set xval 0.55 set ylab 1.25 set asiz 0.55 set vsiz 0.45 set xlab 1.5 * ========================== * Frame and axis labels * ========================== null 4.e-7 1.001 0.035 90200. atit 'x' 'Q^2! (GeV^2!)' * =========================== * prepare labeling * ========================== set cfon -20 set csiz 0.3 set dmod 1 * ------------------------------------ xs=0.1 ys=0.93 xmin=$grafinfo('wnxmin') xmax=$grafinfo('wnxmax') x=[xmin]+[xs]*([xmax]-[xmin]) maxw=$grafinfo('wnymax') y1=$sigma(10.**([ys]*[maxw])) y2=$sigma(10.**(([ys]-0.07)*[maxw])) y3=$sigma(10.**(([ys]-0.14)*[maxw])) y4=$sigma(10.**(([ys]-0.21)*[maxw])) y5=$sigma(10.**(([ys]-0.28)*[maxw])) y6=$sigma(10.**(([ys]-0.36)*[maxw])) y7=$sigma(10.**(([ys]-0.43)*[maxw])) y8=$sigma(10.**(([ys]-0.50)*[maxw])) y9=$sigma(10.**(([ys]-0.57)*[maxw])) y10=$sigma(10.**(([ys]-0.64)*[maxw])) y11=$sigma(10.**(([ys]-0.71)*[maxw])) x=$sigma(10.**([x])) * i=number of lines for labels i=11 y=[y[i]] * ============================= * put fixed target experiments * ============================= * all in common style set faci 5 set plci 5 set fais 1 set fasi 1 i=[i]-1; y=[y[i]]; key [x] [y] 1 'E665, SLAC ' 0 L i=[i]-1; y=[y[i]]; key [x] [y] 1 'CCFR, NMC, BCDMS,' 0 L i=[i]-1; y=[y[i]]; key [x] [y] $GRAFINFO(FACI) 'Fixed Target Experiments:' 0.8 F * CCFR exec plot2 ccf * E665 exec plot2 e665 * SLAC exec plot2 slc * NMC exec plot2 nmc * BCDMS exec plot2 bcd * ============================= * Zeus kinematic regions * ============================= * all in common style set faci 2 set plci 2 set fais 3 set fasi 144 i=[i]-1; y=[y[i]]; key [x] [y] $GRAFINFO(FASI) 'ZEUS' 0.8 H * standard analysis 96 97 exec plot2 z967 * beam pipe trackers exec plot2 zbpt * shifted vertex exec plot2 z95 * ============================= * H1 kinematic regions * ============================= * all in common style set faci 4 set plci 4 set fais 3 set fasi 1 *** H1 shifted vertex 95 set fasi 359 exec plot2 h95 i=[i]-1; y=[y[i]]; key [x] [y] $GRAFINFO(FASI) 'H1 shifted vertex' 0.8 H ** H1 ISR (Cigdem Issever) set fasi 354 exec plot2 hi97 i=[i]-1; y=[y[i]]; key [x] [y] $GRAFINFO(FASI) 'H1 ISR' 0.8 H *** H1 minimum bias 99 set fasi 345 exec plot2 hm99 i=[i]-1; y=[y[i]]; key [x] [y] $GRAFINFO(FASI) 'H1 min. bias' 0.8 H *** H1 low Q2 (Spacal) set fasi 350 exec plot2 hsp i=[i]-1; y=[y[i]]; key [x] [y] $GRAFINFO(FASI) 'H1 low Q^2!' 0.8 H * ** H1 hiq2 fl set fasi 244 exec plot2 hfl i=[i]-1; y=[y[i]]; key [x] [y] $GRAFINFO(FASI) 'H1 high Q^2! F?L!' 0.8 H * ** H1 hiq2 set fasi 344 exec plot2 h1 i=[i]-1; y=[y[i]]; key [x] [y] $GRAFINFO(FASI) 'H1 high Q^2!' 0.8 H * ===================================== * lines of constant y (sqrt(s)=320 GeV set dmod 1 set plci 1 line -7. -1.99482 0. 5.00518 | y=1 line -5. -2.39276 0. 2.60724 | y=0.004 igset chhe 0.35 igset tang 45. itx 0.000025 3. 'Kinematic limit y=1' itx 0.00023 0.045 'y=0.004' * ======================================== pi/pri h1kinplane.eps Return