!!!!!!!!!!!!!!!!!!!!!!!
! NECESSARY POWHEG PARAMETERS
! Basic collider parameters
ih1     1               ! hadron 1 (1 for protons, -1 for antiprotons)
ih2     1               ! hadron 2 (1 for protons, -1 for antiprotons)
ebeam1  6800d0          ! energy of beam 1
ebeam2  6800d0          ! energy of beam 2
!
! Parameters corresponding to the PDF sets
! ndns is 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)
! lhans is to be set only if using LHA pdfs
lhans1  93100          ! pdf set for hadron 1 (LHA numbering)
lhans2  93100          ! pdf set for hadron 2 (LHA numbering)	
!
! Parameters to control the statistics, i.e. number of calls and iterations
ncall1      40000       ! number of calls for initializing the integration grid
itmx1       3          ! number of iterations for initializing the integration grid
ncall2      40000       ! number of calls for computing the integral and finding upper bounding envelope
itmx2       3           ! number of iterations for computing the integral and finding upper bouding envelope
foldcsi     2           ! number of folds on csi integration
foldy       2           ! number of folds on  y  integration
foldphi     2           ! number of folds on phi integration
nubound     20000       ! number of calls to set up the upper bounding norms for radiation
numevts 10000            ! number of events to be generated
!
! 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
!!!!!!!!!!!!!!!!!!!!!!!
! OPTIONAL POWHEG PARAMETERS
withnegweights      1       ! (default 0) if 1, use negative weights
!xupbound           2d0     ! increase upper bound for radiation generation
alphas_from_lhapdf  0       ! obtain the value of alpha_s from the PDF set
!renscfact          1d0     ! (default 1d0) ren scale factor: muren  = muref * renscfact 
!facscfact          1d0     ! (default 1d0) fac scale factor: mufact = muref * facscfact
fastbtlbound        1
storemintupb        0
check_bad_st2       1       ! Check for bad sees after stage 2 and discard
bornonly            0       ! (default 0) if 1 do Born only
#LOevents            1
nores               1       ! DO NOT use resonance histories
withdamp            1       ! (default 0, do not) use Born-zero damping factor
testplots           1       ! (default 0, do not) do NLO and PWHG distributions
allrad 0
!
! RNG input
iseed   4337            !  Start the random number generator with seed iseed
!rand1  0               !  skipping  rand2*100000000+rand1 numbers.         
!rand2  0               !  (see RM48 short writeup in CERNLIB)              
!
! For running on multiple cores in parallel
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
#parallelstage 2
#xgriditeration 1
maxseeds 1000           ! max. number of seeds
!!!!!!!!!!!!!!!!!!!!!!!
! OPTIONAL PROCESS PARAMETERS
! Physical parameters
hmass      125d0           ! Higgs boson mass
hwidth	    0d0             ! Higgs boson width   
wmass      80.379d0        ! W boson mass
wwidth     2.085d0         ! W boson width
zmass      91.1876d0        ! Z boson mass
zwidth     2.4952d0         ! Z boson width
tmass      172.5         ! top quark mass
gfermi     1.16638d-05   ! Fermi constant
!
! Parameters for NLO calculations
runningscales       3       ! 0: use fixed scale, 1: use sum_{i=1}^{npartfin} (pt_i+sqrt(MH^2+ptH^2))/2, 2: use (P_T(p1)+P_T(p2))/2
fakevirt            0       ! Approximate virtuals by alpha*Born or alpha_s*Born
qed_qcd             0       ! 0: NLO-QCD only, 1: NLO-EW only

!!!!!!!!!!!!!!! CUTS USED FOR ANALYSIS and PHASE SPACE: !!!!!!!!!!!!!!!!!!!!!!!
ptalljetmin    20d0 ! Minimum pt for all jets
yjetmax        4.7d0 ! rapidity acceptance for jets
Rsep_jjmin     0.4d0 ! R to be used in jet algorithm
ptjetmin       20d0 ! Minimum pt for two tagging jets
mjjmin         300d0 ! Minimum invariant mass of dijet system
deltay_jjmin   0d0 ! rapidity separation between tagging jets
jet_opphem     0 ! 1: require the tagging jets in opposite detector hemispheres.
!!!!!!!!!!!!!!! END USER PARAMETERS !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
storeinfo_rwgt 1

#rwl_add 1
rwl_file '-'
<initrwgt>
<weightgroup name='nominal' combine='None'>
<weight id='0'>default</weight>
<weight id='1'> renscfact=2d0 facscfact=2d0 </weight>
<weight id='2'> renscfact=0.5d0 facscfact=0.5d0 </weight>
<weight id='3'> renscfact=1d0 facscfact=2d0 </weight>
<weight id='4'> renscfact=1d0 facscfact=0.5d0 </weight>
<weight id='5'> renscfact=2d0 facscfact=1d0 </weight>
<weight id='6'> renscfact=0.5d0 facscfact=1d0 </weight>
</weightgroup>
</initrwgt>

rwl_group_events 20