* ***********************************
* 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
 exec allplot
* exec test
   
 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
  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 0.08 _
    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 min bias 99 -
* --------------------
* preliminary H1 min bias (Doris Eckstein)
*  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
  
* ** V.S.
  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.045 _
		       7.9E-3  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 + extended (99-00) (A. Dubak)
  ve/cr h1_x(12)  r  _
    0.00139791 0.0140 0.0098  0.509  0.80   0.80  0.8  0.5 _
    0.08   0.0155  0.004  .00199
  ve/cr h1_q2(12) r _
    89.125 891.25  891.25 46334. 46334. 1000. 400. 180. _
    180.   141.25 141.25 89.125

  
  
* H1 high Q2 FL (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.
  
* ---------------------------------------
* --- cdf/d0 inclusive jets |eta|<0.7 ---
* ---------------------------------------
*  vec/create cdfd0_x(9) R 1.E-0.2018 1.0E-0.132 1.0E-0.0099 1.0E-0.0992 1.0E-0.3506 1.0E-1.1808 _
*    1.0E-1.6505 1.0E-1.7696 1.0E-1.7166 
*  vec/create cdfd0_q2(9) R 1.0E+3.3004 1.0E+4.2790  1.0E+4.6025  1.0E+5.1996 1.0E+5.3987 _
*    1.0E+4.5818   1.0E+3.6819  1.0E+3.4123  1.0E+3.3004 
  
  vec/create cd0_x(9) R 0.6283 0.9701 0.9774 0.7958 0.4461  0.0659 0.0224 0.016998 0.01920
  
  vec/create cd0_q2(9) R 1.997E+3 19.010E+3 40.040E+3 158.343E+3 250.437E+3 38.176E+3 4.807E+3 2.584E+3 1.997E+3
  
* -------------------------------------------------------
* --- d0 inclusive jets |eta|<3 ,present measuremnts  ---
* -------------------------------------------------------
  vec/create  d0_x(12) R 0.002513 0.00219 0.00192 0.00206 0.01995 0.44607 0.76015 0.89888 0.98492 0.98492 0.94080 0.93368
  
  vec/create  d0_q2(12) R 2.487E+3 2.51E+3 2.815E+3 3.609E+3  38.176E+3 250.553E+3 188.061E+3 _
    134.617E+3 31.535E+3 17.782E+3 2.81E+3 2.487E+3
  
  
  
* ----------------
* --- 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 0
   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 bord 1
   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
* =============================
   set bord 0
   set faci 2
   set plci 2
   set fais 3
   set fasi 144
   exec plot2 h95
   exec plot2 hi97
   exec plot2 hm99
   exec plot2 hsp
   exec plot2 hfl
   exec plot2 h1
   
   
   
   
   
* all in common style
   set bord 1
   set faci 4
   set plci 1
   set fais 3
   set fasi 1
*** H1 shifted vertex 95
   set fasi 359
   exec plot2 h95
   i=[i]-1; y=[y[i]]
   oldfasi=$GRAFINFO(FASI)
   set faci 2; set plci 2; set fasi 144;       key [x] [y] $GRAFINFO(FASI) ' ' 0.8 H
   set faci 4; set plci 4; set fasi [oldfasi]; 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]]
   oldfasi=$GRAFINFO(FASI)
   set faci 2; set plci 2; set fasi 144;       key [x] [y] $GRAFINFO(FASI) ' ' 0.8 H
   set faci 4; set plci 4; set fasi [oldfasi]; 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]]
   oldfasi=$GRAFINFO(FASI)
   set faci 2; set plci 2; set fasi 144;       key [x] [y] $GRAFINFO(FASI) ' ' 0.8 H
   set faci 4; set plci 4; set fasi [oldfasi];   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]]
   oldfasi=$GRAFINFO(FASI)
   set faci 2; set plci 2; set fasi 144;       key [x] [y] $GRAFINFO(FASI) ' ' 0.8 H
   set faci 4; set plci 4; set fasi [oldfasi]; 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]]
   oldfasi=$GRAFINFO(FASI)
   set faci 2; set plci 2; set fasi 144;       key [x] [y] $GRAFINFO(FASI) ' ' 0.8 H
   set faci 4; set plci 4; set fasi [oldfasi]; 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]]
   oldfasi=$GRAFINFO(FASI)
   set faci 2; set plci 2; set fasi 144;       key [x] [y] $GRAFINFO(FASI) ' ' 0.8 H
   set faci 4; set plci 4; set fasi [oldfasi]; 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
 
* **********************
 macro allplot
   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.
   xmin=4.e-7
 xmax=1.1
* xmax=1.5
   ymin=0.035
 ymax=350200.
* ymax=4.e+8
* null 4.e-7 1.001 0.035 250200.
   null [xmin] [xmax] [ymin] [ymax]
   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=8
  
   y=[y[i]]

* =============================
* put fixed target experiments
* =============================
   set bord 1
   facizeus=5
   plcizeus=1
   set faci [facizeus]
   set plci [plcizeus]
   set fais 1
   set fasi 1
   exec plot2 z98
   exec plot2 z967
   exec plot2 zbpt 
   exec plot2 z95
* all in common style
   set bord 1
   facifixed=7
   set faci [facifixed]
   set plci 1
   set fais 3
   set fasi 244
   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
   i=[i]-1; y=[y[i]];  key [x] [y] $GRAFINFO(FASI) 'Fixed Target Experiments:' 0.8 H
   
* CCFR
   exec plot2 ccf
   
* E665
   exec plot2 e665
   
* SLAC
   exec plot2 slc
   
* NMC
   exec plot2 nmc
   
* BCDMS
   exec plot2 bcd
   
   
   
   
   set bord 1
   set faci 4
   set plci 1
   set fais 3
* set fasi 359
   set fasi 105
   exec plot2 d0
* exec plot2 d0
* i=[i]-2; y=[y[i]] ;  key [x] [y] $GRAFINFO(FASI) 'D"2#  Inclusive jets [c]"L#3' 0.8 H
   i=[i]-1; y=[y[i]] ;  key [x] [y] $GRAFINFO(FASI) 'D0  Inclusive jets [c]"L#3' 0.8 H

* CDF
   facicdf=3
   plcicdf=1
   set bord 1
   set faci [facicdf]
   set plci [plcicdf]
   set fais 3
* set fasi 359
   set fasi 159
   exec plot2 cd0
* exec plot2 cd0
* i=[i]-1; y=[y[i]] ;  key [x] [y] $GRAFINFO(FASI) 'CDF/D"2#   Inclusive jets [c]"L#0.7' 0.8 H
   i=[i]-1; y=[y[i]] ;  key [x] [y] $GRAFINFO(FASI) 'CDF/D0   Inclusive jets [c]"L#0.7' 0.8 H
  
   
   
* =============================
* H1 kinematic regions
* =============================
   set faci [facizeus]
   set plci [plcizeus]
   i=[i]-1; y=[y[i]];  key [x] [y] $GRAFINFO(FACI) 'ZEUS' 0.8 F
   
   facih1=2
   plcih1=1
   set bord 1
   set faci [facih1]
   set plci [plcih1]
   set fais 3
* set fasi 235
   set fasi 144
   exec plot2 h95
   exec plot2 hi97
   exec plot2 hm99
   exec plot2 hsp
   exec plot2 hfl
   exec plot2 h1
   i=[i]-1; y=[y[i]] ;  key [x] [y] $GRAFINFO(FASI) 'H1' 0.8 H
     
   
   null [xmin] [xmax] [ymin] [ymax] sa   
   
   pi/pri kinplane.eps
   pi/pri kinplane.ps
 Return
 
 
 macro test
   opt logx
   opt logy
* vec/create h1a_q2() R   0.3     0.3     0.7    .5100     .2830 _    
   .2830      _ 
* vec/create h1a_x()  R   5.E-6   7.E-6   2.E-5  .148E-04  .148E-04  
     _
     
     null 4.e-7 1.001 0.035 90200.
*   exec plot2 z98
*   exec plot2 z967
*   exec plot2 zbpt 
*   exec plot2 z95
   
   
* exec plot2 h95
   exec plot2 hi97
* exec plot2 hm99
* exec plot2 hsp
*   exec plot2 hfl
*   exec plot2 h1
   
 return