# general options
OUTPUT: 3

# event count
EVENTS: 1000000/8  # / num_cpus

# save results
GENERATE_RESULT_DIRECTORY: true

# set up $p$ $\bar{p}$ beams, each at 100 GeV
BEAMS: [2212, 2212]
BEAM_ENERGIES: 6500
PDF_LIBRARY: LHAPDFSherpa
PDF_SET: NNPDF31_lo_as_0118
PDF_SET_VERSIONS: [0, 0]

# SHERPA_LDADD:
# - FlatPDF
# PDF_LIBRARY: FlatPDF
# PDF_SET: FlatPDF

# matrix-element calculation
ME_GENERATORS:
- Comix
- Amegic

SCALES: 'VAR{1/2*Abs2(p[0]+p[1])}'

# the main scattering process
PROCESSES:
- 94 -94 -> 22 22:
    Order: {QCD: 0, EW: 2}
    Integration_Error: 0.001

# use pure QED
EW_SCHEME: alpha0

# disable all advanced processes
FRAGMENTATION: None
SHOWER_GENERATOR: None
MI_HANDLER: None

# cuts
SELECTORS:
- [Eta, 22, -1, 1]
- [PT, 22, 20, 200000000]

# no transverse impulses
BEAM_REMNANTS: false

ANALYSIS: Rivet
RIVET:
  IGNOREBEAMS: 1
  ANALYSES:
    - MC_DIPHOTON_PROTON