Program Manual and input-file guide for HiggsPt_R version 1.1

This manual describes how to use the program HiggsPt_R, which computes the Higgs boson Transverse Momentum distribution at NLO in the limit of large top quark mass, with input grid and output files for use with the resummation program RESBOS.
The code produces two output files. The "xxx_Pert.out" output file displays the order alpha_s^3 Pert_1 and order alpha_s^4 Pert_2 contributions to the perturbative cross sections and the K factor (K=(Pert_1+Pert_2)/Pert_1) as a function of the input grid variables. The "xxx_Y.out" output file displays the order alpha_s^3 Y_1 and order alpha_s^4 Y_2 contributions to Y=sigma(exact)-sigma(small-pt) and the associated K factor (K=(Y_1+Y_2)/Y_1) as a function of the input grid variables. The manual gives line-by-line instructions for the input file "HiggsPt_R.indat".


General Formatting Notes:

It should be fairly easy to use the input file, "HiggsPt_R.indat", just by modifying the example that came with the source distribution. However, there are a few general formatting requirements to be aware of.
The input file must have exactly 10 lines.
In the following explanations of the inputs, a list {a, b, c} indicates a choice of string inputs on a line, while an ordered list [x y z] indicates multiple numerical parameters input on a single line. Possible inputs are colored red.
If the input to a line is a string, then it MUST have trailing white spaces after the string (NOT TABS!). It is easiest to do this by keeping the explanatory comment after the input string (i.e, white spaces followed by #x Comment) from the example file.
Any input which is a string must be spelled exactly as given below, including capitalization and without quotation marks, or else the program will complain.
The program ignores anything in a line that comes after the expected input.


Line-by-line explanation of input file:

#1 Name of the output files.
Input a string, such as output/HiggsPt_R (without quotes) to name the output files. The program generates two output files by appending the strings "_Pert.out" and "_Y.out" to the input string. (In the example, the program would produce the files "output/HiggsPt_R_Pert.out" and "output/HiggsPt_R_Y.out". The empty directory "output/" may be used for output files. Or you may organize as you desire.
#2 Type of collider.
Choose from {pp, ppbar}.
#3 Collider COM energy
Input in GeV, such as 14000.0
#4 Type of PDF set.
Choose from {cteq5, cteq6, mrst99, mrst2001lo, mrst2002, mrst2004}.
#5 PDF individual set.
Input an integer, which specifies the specific PDF set to use.
For Cteq distributions this is the variable iset, which chooses among the various distributions in the set.
For MRST distributions this is the variable mode.
Relevant choices for this program are:
cteq6 :     iset=1 : CTEQ6M Standard MSbar scheme at NLO
  iset=3 : CTEQ6L Leading Order with NLO running
  iset=4 : CTEQ6L1 Leading Order with LO running
  iset=200 : CTEQ6.1M (updated CTEQM)
  iset=400 : CTEQ6.6M (most recently updated CTEQM)
cteq5 : iset=3 : CTEQ5L Leading Order with LO running
  iset=8 :    CTEQ5M1 Standard MSbar scheme at NLO
mrst2004 : mode=1 : Best fit MSBAR at NLO
mrst2002 : mode=1 : Best fit MSBAR at NLO
mrst2001lo : mode=1 : Leading Order best fit
mrst99 : mode=1 : COR01 central gluon, a_s (MSBAR at NLO)
  mode=2 :    COR02 higher gluon (MSBAR at NLO)

Other choices of distributions in the cteq5, cteq6, and mrst99 series may also be called, but first a symbolic link must be made in the directory "main" to the relevant ".tbl", ".pds", or ".dat" file in the PDF directories. In particular, symbolic links to the 40 Cteq6.6 Error Matrix PDFs can be made by typing "make cteq6_error_PDFs" on the command line. For more info on PDF distributions, look at the comments in the file "Cincludes/pdfqcd07.h" or in the specific PDF code files in the directories "PDFs/Cteq6-2008", "PDFs/Cteq5" "PDFs/MRST99", "PDFs/MRST2001lo", "PDFs/MRST2002", and "PDFs/MRST2004".

#6 Choice of Ren/Fact Scale Variable.
Choose from {MH, MT}.
MH :       muScale = M_H   (Higgs Mass)
MT : muScale = M_T = sqrt(M_H^2+p_T^2)  (Higgs Transverse Mass)
#7 Input scale choices.
The input format is [chiR chiF]. For example, 1.0 1.0.
The scale factors are defined by:
#8 Name of the Q grid file.
Input a string, such as q_grid_h0.inp (without quotes) to identify the input grid file in the variable Q (the Higgs Boson Mass). This file should be organized with one Q value on each line and nothing else.
#9 Name of the y grid file.
Input a string, such as y_grid_h0.inp (without quotes) to identify the input grid file in the variable y (the Higgs Boson rapidity). This file should be organized with one y value on each line and nothing else.
#10 Name of the qt grid file.
Input a string, such as qt_grid_h0.inp (without quotes) to identify the input grid file in the variable qt (the Higgs Boson transverse momentum). This file should be organized with one qt value on each line and nothing else.

Renormalization :      mu_R=chiR*muScale.
Factorization : mu_F=chiF*muScale.