!!!!!!!!!   parameters that define the process !!!!!!!!!!!!!!!!!!
higgstype 3 ! higgs type: 1 = light CP even; 2 = heavy CP even; 3 = pseudoscalar
ih1   1           ! hadron 1 (1 for protons, -1 for antiprotons)
ih2   1           ! hadron 2 (1 for protons, -1 for antiprotons)
ebeam1 6500     ! energy of beam 1
ebeam2 6500    ! energy of beam 2
hdecaymode 8      ! PDG code for first decay product of the higgs
! 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
hwidth 10d0	
hdecaywidth 0	! 1 read Higgs decay width from HDECAY file
zerowidth 0 ! Control if the Higgs boson is to be produced on-shell or not: 1 = On-Shell; 0 = Off-shell with Breit-Wigner
bwshape 1	! BW shape: 1 running width BW, 2 fixed width BW
masswindow_high 5d0 ! number of widths around hmass in the BW for an off-shell Higgs boson  upper limit
masswindow_low 5d0  ! number of widths around hmass in the BW for an off-shell Higgs boson  lower limit
ew 0                     ! ew = 0 disable EW corrections - ew = 1 enable EW corrections
fastew 1                   ! fast ew corrections evaluations by mass sampling
nosquarks 0 ! 1 = No scalars, 0/not defined = scalars enabled
lhsquarks 0 ! 1 = Scalars contribution are taken in the infinite mass limit, 0 = full mass dependence


!!!!!!!!!!!!  parameters of the simulation   !!!!!!!!!!!!!!!!!!!!
numevts NEVENTS


!!!!!!!!!!!!  QCD parameters   !!!!!!!!!!!!!!!!!!!!
!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
lhans1 325300       ! pdf set for hadron 1 (LHA numbering)
lhans2 325300      ! 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 
! mur,muf settings
runningscale 0    ! 0 = scales equal to the Higgs pole mass; 1  = scales equal to the Higgs virtuality; 
                  ! 2 = scales equal to the Higgs pole mass for Born-like configuration and to the transverse mass for real emission contribution
! When using 2 uncomment the following option
! btlscalereal 1
renscfact  1   ! (default 1d0) ren scale factor: muren  = muref * renscfact
facscfact  1   ! (default 1d0) fac scale factor: mufact = muref * facscfact
!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
!charmthrpdf  1.5  ! (default 1.5 GeV) pdf charm treshold
!bottomthrpdf 5.0  ! (default 5.0 GeV) pdf bottom treshold
!hfact    %hfact    ! (default no dumping factor) dump factor for high-pt radiation: > 0 dumpfac=h**2/(pt2+h**2)



!!!!!!!!!!!!  generator parameters     !!!!!!!!!!!!!!!
! 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
! grids preparation
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 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

!ptsupp     0d0   ! (default 0d0)  mass param for Born suppression factor (generation cut) If < 0 suppfact = 1
!bornonly   0      ! (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
testplots  1      ! (default 0, do not) do NLO and PWHG distributions
!testsuda  1       ! (default 0, do not test) test Sudakov form factor
!radregion 1       ! (default all regions) only generate radiation in the selected singular region  
iseed SEED
!rand1    32      ! initialize random number sequence 
!rand2    0      ! 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 



!!!!!!!!!!!! physical parameters  !!!!!!!!!!!!!!!
gfermi 0.116637D-04        ! GF

!!!!!!!!!!! MSSM parameter !!!!!!!!!
scheme 2      ! Mass renormalization scheme. 2 = OS/Feynhiggs, 1 = DRBAR

!!!!!!!!!!! MSSM OS parameter  !!!!!!!!!
FHslhaflag 0  ! relevant when scheme=2 . If 0 read FeynHiggs input param from powheg-fh.in, otherwise from powheg-fh.slha
FHdecaywidth 1 ! 1 = use the decay width provided by FeynHiggs
! bottommass 4.75d0 ! If not set, the bottom OS mass is calculated with a 1-loop scheme change formula.


!!!!!!!!!!! MSSM DRbar parameter  !!!!!!!!!
mssmmbatmh 0 ! 1 = mb(mh) is used (from mssm-param-mh.slha), 0 = mb(M_susy) is used
#pdfreweight 0
#storeinfo_rwgt 0
#withnegweights 1 1 ! default 0, 


#manyseeds 1

#parallelstage 4

#xgriditeration 1 1
rwl_group_events 2000
lhapdf6maxsets 50
rwl_file 'pwg-rwl.dat'
rwl_format_rwgt 1
