! gg_H production parameters:
hmass  1000.       ! Higgs boson mass
hwidth 0.00407     ! Higgs boson width
masswindow 30d0
 
! GGF_H production:
! **** Mandatory parameters for ALL models ****
massren 0           ! Mass renormalization scheme. 0 = OS, 1 = MSBAR , 2 = DRBAR
zerowidth 0         ! Control if the Higgs boson is to be produced on-shell or not: 1 = On-Shell; 0
model 0             ! model: 0 = SM
ew 1                ! ew = 0 disable EW corrections - ew = 1 enable EW corrections
!passarino 1 (by "ew")       ! Use Passarino's grid instead of ABDV functions.
!fastew 1 (by "ew")          ! fast ew corrections evaluations by mass sampling
bottommass 4.95d0   ! bottom quark mass - if defined it enables the bottom quark
charmmass 1.55d0    ! char quark mass - if defined it enables the charm quark
! Optional
hdecaywidth 0        ! If equals to 1 read total decay width from HDECAY sm.br2 file
bwshape 1
 
Zmass 91.1876d0
Zwidth 2.4952d0
Wmass 80.399d0
Wwidth 2.085d0
alphaem 0.78165303985d-02
sthw2 0.23113d0
gfermi 0.116637D-04  ! GF
topmass 172.5        ! top quark mass
hdecaymode 0         ! PDG code for Higgs boson decay products (it affects only the SMC)
!                       allowed values are:
!                       0 all decay channels open
!                       1-6 d dbar, u ubar,..., t tbar (as in HERWIG)
!                       7-9 e+ e-, mu+ mu-, tau+ tau-
!                       10  W+W-
!                       11  ZZ
!                       12  gamma gamma
!                       -1  all decay channels closed
 
# Not used in ggH
#CKM_Vud 0.97428
#CKM_Vus 0.2253
#CKM_Vub 0.00347
#CKM_Vcd 0.2252
#CKM_Vcs 0.97345
#CKM_Vcb 0.0410
#CKM_Vtd 0.00862
#CKM_Vts 0.0403
#CKM_Vtb 0.999152
 
!
pdfreweight 1
 
! general parameters:
 
numevts 1000000  ! 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)
#ndns1 131         ! pdf set for hadron 1 (mlm numbering)
#ndns2 131         ! pdf set for hadron 2 (mlm numbering)
ebeam1 4000d0     ! energy of beam 1
ebeam2 4000d0     ! energy of beam 2
 
! To be set only if using LHA pdfs
lhans1 21100 ! pdf set for hadron 1 (LHA numbering)
lhans2 21100 ! pdf set for hadron 2 (LHA numbering)
! To be set only if using different pdf sets for the two incoming hadrons
! QCDLambda5  0.25 ! for not equal pdf sets
 
! 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; <>
 
ncall1 120000  ! number of calls for initializing the integration grid
itmx1     5    ! number of iterations for initializing the integration grid
ncall2 250000  ! number of calls for computing the integral and finding upper bound
itmx2     5    ! number of iterations for computing the integral and finding upper bound
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 50000  ! number of bbarra calls to setup norm of upper bounding function
icsimax   1    ! <= 100, number of csi subdivision when computing the upper bounds
iymax     1    ! <= 100, number of y subdivision when computing the upper bounds
xupbound 2d0   ! increase upper bound for radiation generation
 
! OPTIONAL PARAMETERS
 
#renscfact  1d0    ! (default 1d0) ren scale factor: muren  = muref * renscfact
#facscfact  1d0    ! (default 1d0) fac scale factor: mufact = muref * facscfact
#ptsupp     0d0    ! (default 0d0)  mass param for Born suppression factor (generation cut) If < 0 s
#bornonly   1      ! (default 0) if 1 do Born only
#smartsig   0      ! (default 1) remember equal amplitudes (0 do not remember)
#withsubtr  0      ! (default 1) subtract real counterterms (0 do not subtract)
#withdamp    1     ! (default 0, do not use) use Born-zero damping factor
#ptsqmin    0.8    ! (default 0.8 GeV) minimum pt for generation of radiation
#charmthr   1.5    ! (default 1.5 GeV) charm treshold for gluon splitting
#bottomthr  5.0    ! (default 5.0 GeV) bottom treshold for gluon splitting
#testplots  1      ! (default 0, do not) do NLO and PWHG distributions
hfact    833.3333333333d0    ! (default no dumping factor) dump factor for high-pt radiation: > 0 du
#testsuda  1       ! (default 0, do not test) test Sudakov form factor
#radregion 1       ! (default all regions) only generate radiation in the selected singular region
#charmthrpdf  1.5  ! (default 1.5 GeV) pdf charm treshold
#bottomthrpdf 5.0  ! (default 5.0 GeV) pdf bottom treshold
 
iseed    3  ! initialize random number sequence
rand1    -1     ! initialize random number sequence
rand2    -1     ! initialize random number sequence
 
#iupperisr 1 ! (default 1) choice of ISR upper bounding functional form
#iupperfsr 2 ! (default 2) choice of FSR upper bounding functional form