INCNLO-direct photon and inclusive hadron production code

INCNLO version 1.4

Aims of the program

This program can compute
  1. the NLO hadro-production of high transverse momentum photon combining
    1. the code from Aurenche-Baier-Fontannaz-Schiff for NLO direct (abfs/pholib) and
    2. the code from Aversa-Chiappetta-Greco-Guillet for NLO fragmentation (pioincl/hadlib)
  2. which should reproduce results presented in hep-ph/9811382 where all references can be found.
  3. the NLO hadro-production of high transverse momentum hadron using the code from Aversa-Chiappetta-Greco-Guillet which should reproduce results presented in hep-ph/9910252 where all references can be found.

Structure of the program

INCNLO program directories
alfas merged frag hadlib pholib bases
Directory with two loop's strong coupling evaluation solving implicit equation (author, L.Bourhis) Directory with the merged photon/hadron program Directory with the fragmentation functions packages bfg_photon (photon fragmentation of Bourhis-Fontannaz-Guillet), and hadron (hadron fragmentation of Binnewies-Kniehl-Kramer, of Kniehl-Kramer-Poetter and of Bourhis-Fontannaz-Guillet-Werlen) Directory with the hadronic NLO package Directory with the photon NLO package Directory with the BASES/SPRING package from S.Kawabata
INCNLO user directories
working merged doc
Directory which contains the files: "param.dat", the executable "inclnll", some examples Directory containing the merged photon/hadron program Directory with some documentation

Download the program inc.tar.gz (version 1.4, december 2002)

To run the program

  • Some of the program features have been tested on DEC Alpha (Makefile_lapp), Linux PC (Makefile_home) and RSPLUS IAX (Makefile_cern). Select one of these in the Makefile driver. Some modifications may be needed to adapt the chosen Makefile to your computer system.
  • Verify that the CERN library (see CERNLIB ) and PDFLIB (see PDFLIB in CERN/IT/DOC ) are installed on your computer and properly linked in the chosen Makefile stored in merged/src

  • Run the Makefile driver (with "make") to create the libraries in pholib/lib, hadlib/lib, frag/bfg_photon/lib, frag/hadron/lib, alfas/lib and bases/lib. Then, with "make working", produce the executable file which will be called "inclnll" and stored in the directory working.

  • Select your input via the data file param.dat (see commented example in directory working). You want to calculate the cross section for
    IH1 + IH2 -> IH3 + X.
    First select the Path for the output files.
    Then select the structure functions from pdflib with standard ntype, ngroup and nset

    Make your choice of projectile: IH1=0 (proton), 1 (anti-proton), 3 (nucleus), 4 (pi+), 5 (pi-)
    and target: IH2=0 (proton), 3 (nucleus).
    If IH1 and/or IH2=3, define the projectile and/or target nucleus with Z and A .

    You have the choice of the following produced particles
  • IH3= 0 photon
  • IH3= 1 (pi^+ + pi^-)/2
  • IH3= 2 (K^+ + K^-)/2
  • IH3= 3 (K^0 + K^0_bar)/2
  • IH3= 5 pi^0
  • IH3= 7 h^+ + h^-

  • Then LHIOR selects the precision of the calculation
  • 0 BORN only (for photon production)
  • 1 BORN + ANOMA (for photon production)
  • 2 NLO abfs (for photon production)
  • 3 NLO abfs + NLO fragmentation (for photon production)
  • 4 NLO Hadron production

  • IHIOR : selects lowest order (if 1) or includes higher order (if 2) in Hadron production (should be always 2)
    IANORD: fragmentation set choice
  • selects Bourhis-Fontannaz-Guillet NLO photon fragmentation function set I (if 1) or set II (if 2);
  • selects NLO hadron fragmentation functions of BKK (if 1), of KKP (if 2), of BFGW (if 3).

  • ILOOP : number of loops in alphas calculation. Note that alfas is calculated from lambda used in the structure functions (which may not be the one used in fragmentation functions).

    The others inputs shoud be straigtforward.

  • The photon production program works only with N_F=4. It is not the case for the hadron production program.
  • There is a limitation on the number of P_t points computed by the program, this number must be less or equal to 20.
  • This program does not give the errors on the computed cross sections. They can be estimated by varying the number of integration points (for the direct part for the photon) or looking at test.out for the fragmentation part or hadron production.

  • The output files will be collected in the directory defined by the path for outputs files in param.dat (directory which you have to create). The files are:
    "error.out" error messages
    "test.out" parameters read and some intermediate results
    "result.out" results of integration of cross sections
    Note that the default output file, used by bases and pdflib for example, is not saved automatically


    Back to Homepage