H1 data for fits:



QCD Parameterizations
from the H1 QCD fits (DESY 99-107)

You may need the CERNLIB routine dchsum.f
PDF's are available in MS-bar and DIS scheme via subroutine calls to:

call pdf_xxxx_ms(x,q2,uv,dv,us,ds,ss,cs,bs,gl) (MS-bar scheme)
call pdf_xxxx_ds(x,q2,uv,dv,us,ds,ss,cs,bs,gl) ( DIS scheme)

For MRSH they are also available in leading order:
call pdf_xxxx_lo(x,q2,uv,dv,us,ds,ss,cs,bs,gl)

NC structure functions

(including EW terms) are available via function calls to:

FL=flncxxxx00(x,q2)
F2=f2ncxxxx00(x,q2)
xF3=f3ncxxxx00(x,q2)

CC structure functions for electron scattering:

WL=wlcexxxx00(x,q2)
W2=w2cexxxx00(x,q2)
xW3=w3cexxxx00(x,q2)

CC structure functions for positron scattering:

WL=wlcpxxxx00(x,q2)
W2=w2cpxxxx00(x,q2)
xW3=w3cpxxxx00(x,q2)

where xxxx determines the fit type:

The parametrisations are valid for 4 < Q2 < 47000 GeV2 and 10-5 < x < 0.97

Filename Description Filename Description
pdfh1_hiq2.f PDF's from the high Q2 fit (xxxx = hiq2) stfh1_hiq2.f SF's from the high Q2 fit (xxxx = hiq2)
pdfh1_loq2.f PDF's from the low Q2 fit (xxxx = loq2) stfh1_loq2.f SF's from the low Q2 fit (xxxx = loq2)
pdfh1_loal.f PDF's from the low Q2 fit with alphas=0.115 (xxxx = loal) stfh1_loal.f SF's from the low Q2 fit with alphas=0.115 (xxxx = loal)
pdfh1_lal1.f PDF's from the low Q2 fit with alphas=0.123 (xxxx = lal1) stfh1_lal1.f SF's from the low Q2 fit with alphas=0.123 (xxxx = lal1)
pdfh1_lal2.f PDF's from the low Q2 fit with alphas=0.113 (xxxx = lal2) stfh1_lal2.f SF's from the low Q2 fit with alphas=0.113 (xxxx = lal2)
pdfh1_lal3.f PDF's from the low Q2 fit with alphas=0.108 (xxxx = lal3) stfh1_lal3.f SF's from the low Q2 fit with alphas=0.108 (xxxx = lal3)
pdfh1_lal4.f PDF's from the low Q2 fit with alphas=0.128 (xxxx = lal4) stfh1_lal4.f SF's from the low Q2 fit with alphas=0.128 (xxxx = lal4)
pdfh1_mrsh.f PDF's from MRSH in Leading Order (xxxx = mrsh) stfh1_mrsh.f SF's from MRSH (xxxx = mrsh)
pdfh1_ct5d.f PDF's from CTEQ-5D (xxxx = ct5d) stfh1_ct5d.f SF's from CTEQ-5D (xxxx = ct5d)


You should be able to plot the functions in paw, depending on your set-up, by doing the following.
Note the second statement may take some seconds to complete:

call dchsum.f77
call stfh1_hiq2.f77
fun1 100 f2nchiq200(10**x,100.0) 50 -3 0

qcdfunc.f is a collection of functions and subroutines to help plot cross-sections for e+ or e- scattering at any beam energy. You will also find functions to plot integrated cross-sections ds/dq2 and ds/dx as well as total cross-sections. The functions are designed to be used in paw and require 3 vector inputs:

elepin(1) --> input lepton beam energy
epin(1) --> input proton beam energy
qin(1) --> input lepton charge

Some functions overide these vectors - please look at the code. EW parameters have been fixed. There are also functions to allow you to reweight one cross-section to another e.g. nc-->cc, or e+ --> e- etc

F2 parameterizations

Following are a number of routines that calculate the proton structure function F2 according to the cited parameterisations. The corresponding source and data files can be downloaded for implementation in private code.



Filename Description Reference
f2allm.f ALLM F2 parameterisation for low Q2 HEP-PH/9712415
f2pqcd96.f H1 F2 as presented in Jerusalem 1997  
f2pqcd96.data   f2jeru.ps
f2p96ext.f H1 F2 as presented in Jerusalem 1997 extended to high Q2 f2jeru.ps
f2p96ext.data    



Parameterizations of the gluon density in the proton

This routine returns the gluon density in the proton.

Filename Description Reference
glujeru.f H1 gluon density as presented in Jerusalem 1997 f2jeru.ps
glujeru.data    



Parameterizations of parton density functions (pdfs)

A number of parameterizations of parton density functions for implementation in private code.



Filename Description Reference
f2d3_94_qcd.f Parton distributions from 1994 f2d3 DESY-97-158 / Z.Phys. C76 (1997) 613






 Last updated  12.10.99   by E. Rizvi