# minimum suggested values for ncall and itmx: 
# (ncall1,itmx1,ncall2,itmx2) = (250K,5,1M,5)

idDM  1000021   ! id of DM particles 
                ! (+- 19 is the default, but might not always work.
                ! Using the id of neutralinos is recommended)


vdecaymode  -1 ! 1 for V-mediated fermionic DM
                 ! 2 for A-mediated fermionic DM
                 ! -1, -2 for full theory (mediator mass and width needed)

DMmass      350 ! DM-candidate mass

##############################
# EFT
DMLambda     500 ! UV scale
##############################
# FULL
DMVmass    1750   ! mediator mass
DMVwidth   98.1796574297 ! mediator width
DMgSM      0.25   ! q-q-med coupling
DMgDM      1   ! X-Xbar-med coupling
#runningwidth 0  ! (default 0) if 1, use running width in propagator
###############################

#testplots  0      ! (default 0, do not) do NLO and PWHG distributions
#bornonly 0        ! (default 0) if 1 do Born only in fixed-order
#LOevents 0        ! (default 0) if 1, produce LOPS events (scalup=ptj)

runningscale 3     ! (default 3) 0, 1, 2, 3 correspond to 2mX (fixed), ptj, XXbar invmass, Ht/2
                   ! Details in Born_phsp.f
renscfact  1       ! (default 1d0) ren scale factor: muren  = muref * renscfact 
facscfact  1       ! (default 1d0) fac scale factor: mufact = muref * facscfact 

bornktmin  1       ! (default 0d0) generation cut. Minimum kt in underlying Born
bornsuppfact 400   ! (default 0d0) mass param for Born suppression factor. If < 0 suppfact = 1
withnegweights 0   ! (default 0) if 1, use negative weights

smartsig 1
doublefsr 1

numevts NEVENTS   ! number of events to be generated
ih1   1           ! hadron 1 (1 for protons, -1 for antiprotons)
ih2   1           ! hadron 2 (1 for protons, -1 for antiprotons)
ebeam1 6500d0     ! energy of beam 1
ebeam2 6500d0     ! energy of beam 2

! To be set only if using internal (mlm) pdfs
! ndns1 131         ! pdf set for hadron 1 (mlm numbering)
! ndns2 131         ! pdf set for hadron 2 (mlm numbering)

! To be set only if using LHA pdfs
! 21000 MSTW2008 (LO central)
! 21100 MSTW2008 (NLO central)

lhans1  260000        ! pdf set for hadron 1 (LHA numbering)
lhans2  260000        ! pdf set for hadron 2 (LHA numbering)	

! Parameters to allow or not the use of stored data
use-old-grid    1 ! if 1 use old grid if file pwggrids.dat is present (<> 1 regenerate)
use-old-ubound  1 ! if 1 use norm of upper bounding function stored in pwgubound.dat, if present; <> 1 regenerate

ncall1 50000    ! number of calls for initializing the integration grid
itmx1    5      ! number of iterations for initializing the integration grid
ncall2 50000    ! number of calls for computing the integral and finding upper bounding envelope
itmx2    5      ! number of iterations for computing the integral and finding upper bouding envelope
foldcsi   1     ! number of folds on csi integration
foldy     1     ! number of folds on  y  integration
foldphi   1     ! number of folds on phi integration
nubound   100000  ! number of calls to set up the upper bounding norms for radiation

! OPTIONAL PARAMETERS
xupbound 2d0      ! increase upper bound for radiation generation by this factor

#iseed    SEED     !  Start the random number generator with seed iseed
#rand1     0       !  skipping  rand2*100000000+rand1 numbers.         
#rand2     0       !  (see RM48 short writeup in CERNLIB)              

#manyseeds 1       ! Used to perform multiple runs with different random
                   ! seeds in the same directory.
                   ! If set to 1, the program asks for an integer j.
                   ! The file pwgseeds.dat at line j is read, and the
                   ! integer at line j is used to initialize the random
                   ! sequence for the generation of the events.
                   ! The event file is called pwgevents-'j'.lhe

storeinfo_rwgt 1
pdfreweight 1