ZEUS interface

Here is a short demonstration of how to access the ZEUS transfer matrix and bin interface.

using PartonDensity, CSV, DelimitedFiles

Transfer matrix

eMPp = 1 # e+/e- switch 0/1

1

Read in an example integrated cross section

numbers_from_file = readdlm("data/HERAPDF20_NNLO_EIG_ePp.txt")

429×3 Matrix{Float64}:
   400.0    650.0  49.414
   400.0    650.0  52.149
   400.0    650.0  27.263
   400.0    650.0  33.421
   400.0    650.0  18.594
   400.0    650.0  16.348
   400.0    650.0  21.211
   400.0    650.0  12.374
   400.0    650.0  11.224
   400.0    650.0  10.247
     ⋮             
  8000.0  11000.0   0.035885
  8000.0  11000.0   0.020569
  8000.0  11000.0   0.010562
  8000.0  11000.0   0.00056915
 11000.0  20000.0   0.6079
 11000.0  20000.0   0.60542
 11000.0  20000.0   0.020995
 20000.0  30000.0   0.13832
 20000.0  30000.0   0.0061146

List of integrated cross section values in 429 bins

integ_xsec = numbers_from_file[:,3]

429-element Vector{Float64}:
 49.414
 52.149
 27.263
 33.421
 18.594
 16.348
 21.211
 12.374
 11.224
 10.247
  ⋮
  0.035885
  0.020569
  0.010562
  0.00056915
  0.6079
  0.60542
  0.020995
  0.13832
  0.0061146

Corresponding list of expected event numbers

prediction = get_pred_N(integ_xsec, eMPp);

integ_xsec[153]

prediction[151]

0.15272861711226063

Bins in detector space

The transfer matrix projects into a list of 153 bins with irregular x and Q2 spacing. The bin edges are provided in data/zeus_bin_edges.jl but can also be quickly and easily accessed with the helper function get_bin_info as shown below.

get_bin_info(10)

([800, 950], [0.059, 0.076])

This page was generated using Literate.jl.