Skip to content

Commit

Permalink
Merge pull request #3310 from alintulu/run3-hbb
Browse files Browse the repository at this point in the history
Add Run 3 13p6TeV Higgs cards
  • Loading branch information
menglu21 authored Dec 7, 2022
2 parents 7f850c8 + 7625d99 commit 53f69ee
Show file tree
Hide file tree
Showing 13 changed files with 857 additions and 0 deletions.
Original file line number Diff line number Diff line change
@@ -0,0 +1,60 @@
numevts NEVENTS
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

hmass 125
hwidth 0.004221d0

hdecaywidth 0 ! If equals to 1 read total decay width from HDECAY sm.br2 file
! if 0 the hwidth value is used


bwcutoff 15 ! how many W widths below and above the Higgs mass are included

! 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)

! 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 150000 ! number of calls for initializing the integration grid
itmx1 1 ! number of iterations for initializing the integration grid
ncall2 100000 ! number of calls for computing the integral and finding upper bound
itmx2 5 ! number of iterations for computing the integral and finding upper bound
fastbtlbound 1 ! (default 0) if 1 use fast btilde bound
storemintupb 1 ! store calls for grids
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 2000 ! number of bbarra calls to setup norm of upper bounding function

! OPTIONAL PARAMETERS

testplots 1 ! (default 0, do not) do NLO and PWHG distributions

withnegweights 1 ! (1 default) If 1 output negative weighted events.
! If 0 descard them

renscfact 1d0 ! (default 1d0) ren scale factor: muren = muref * renscfact
facscfact 1d0 ! (default 1d0) fac scale factor: mufact = muref * facscfact

bornktmin 200.0 ! Minimum transverse momentum of the Higgs at the underlying Born level

minlo 1 ! default 0, set to 1 to use minlo

factsc2min 2 ! at this value the factorization scale is frozen (neede with minlo)

sudscalevar 1 ! (default 1) scale variation also in Sudakov form factors in minlo

iseed SEED ! initialize random number sequence

# variables for including quark mass effects
quarkmasseffects 1 ! Include quark mass effects
topmass 172.5d0 ! top quark mass (172.5 in HNNLO-patches/mdata.f)
bottommass 3.38d0 ! bottom quark mass in MSbar at MH
bmass_in_minlo 1 ! Include quark mass effects in Sudakov exponent
Original file line number Diff line number Diff line change
@@ -0,0 +1,70 @@
vdecaymode 0 ! 1: e, 2: mu, 3: tau, 4:nu_e, 5: nu_mu, 6: nu_tau, 0: hadronic, 10: inclusive, 11: inclusive leptons, 12: inclusive neutrinos
hdecaymode -1 ! for PYTHIA or HERWIG: do not decay the Higgs boson

hmass 125d0 ! Higgs boson mass
hwidth 0.00407d0 ! Higgs boson width

min_h_mass 12d0
max_h_mass 1250d0

min_z_mass 10d0
max_z_mass 1000d0



numevts NEVENTS
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

bornktmin 0.26d0 ! (default 0d0) generation cut. Minimum kt in underlying Born
bornsuppfact 0d0 ! (default 0d0) mass param for Born suppression factor. If < 0 suppfact = 1

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

renscfact 1d0 ! (default 1d0) ren scale factor: muren = muref * renscfact
facscfact 1d0 ! (default 1d0) fac scale factor: mufact = muref * facscfact

! 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

! the typical call uses 1/1400 seconds (1400 calls per second)
ncall1 100000 ! Number of calls for the construction of the importance sampling grid
itmx1 1 ! Number of iterations for grid
ncall2 100000 ! Number of calls for the computation of the upper bounding envelope
! for the generation of radiation
itmx2 1 ! Number of iterations for the above
! Notice: the total number of calls is ncall2*itmx2*foldcsi*foldy*foldphi
foldcsi 2 ! number of folds on csi integration
foldy 5 ! number of folds on y integration
foldphi 2 ! number of folds on phi integration

nubound 100000 ! number of calls to set up the upper bounding norms for radiation
! This is performed using only the Born cross section (fast)


! OPTIONAL PARAMETERS

withnegweights 1 ! (default 0) if on (1) use negative weights
testplots 1 ! (default 0, do not) do NLO and PWHG distributions

iseed SEED

xgriditeration 1
fastbtlbound 1

minlo 1
minlo_nnll 1

# if running with minlo, set the following to 0
massivetop 0

sudscalevar 1

doublefsr 1

nohad 1

Original file line number Diff line number Diff line change
@@ -0,0 +1,70 @@
vdecaymode 11 ! 1: e, 2: mu, 3: tau, 4:nu_e, 5: nu_mu, 6: nu_tau, 0: hadronic, 10: inclusive, 11: inclusive leptons, 12: inclusive neutrinos
hdecaymode -1 ! for PYTHIA or HERWIG: do not decay the Higgs boson

hmass 125d0 ! Higgs boson mass
hwidth 0.00407d0 ! Higgs boson width

min_h_mass 12d0
max_h_mass 1250d0

min_z_mass 10d0
max_z_mass 1000d0



numevts NEVENTS
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

bornktmin 0.26d0 ! (default 0d0) generation cut. Minimum kt in underlying Born
bornsuppfact 0d0 ! (default 0d0) mass param for Born suppression factor. If < 0 suppfact = 1

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

renscfact 1d0 ! (default 1d0) ren scale factor: muren = muref * renscfact
facscfact 1d0 ! (default 1d0) fac scale factor: mufact = muref * facscfact

! 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

! the typical call uses 1/1400 seconds (1400 calls per second)
ncall1 100000 ! Number of calls for the construction of the importance sampling grid
itmx1 1 ! Number of iterations for grid
ncall2 100000 ! Number of calls for the computation of the upper bounding envelope
! for the generation of radiation
itmx2 1 ! Number of iterations for the above
! Notice: the total number of calls is ncall2*itmx2*foldcsi*foldy*foldphi
foldcsi 2 ! number of folds on csi integration
foldy 5 ! number of folds on y integration
foldphi 2 ! number of folds on phi integration

nubound 100000 ! number of calls to set up the upper bounding norms for radiation
! This is performed using only the Born cross section (fast)


! OPTIONAL PARAMETERS

withnegweights 1 ! (default 0) if on (1) use negative weights
testplots 1 ! (default 0, do not) do NLO and PWHG distributions

iseed SEED

xgriditeration 1
fastbtlbound 1

minlo 1
minlo_nnll 1

# if running with minlo, set the following to 0
massivetop 0

sudscalevar 1

doublefsr 1

nohad 1

Original file line number Diff line number Diff line change
@@ -0,0 +1,70 @@
vdecaymode 12 ! 1: e, 2: mu, 3: tau, 4:nu_e, 5: nu_mu, 6: nu_tau, 0: hadronic, 10: inclusive, 11: inclusive leptons, 12: inclusive neutrinos
hdecaymode -1 ! for PYTHIA or HERWIG: do not decay the Higgs boson

hmass 125d0 ! Higgs boson mass
hwidth 0.00407d0 ! Higgs boson width

min_h_mass 12d0
max_h_mass 1250d0

min_z_mass 10d0
max_z_mass 1000d0



numevts NEVENTS
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

bornktmin 0.26d0 ! (default 0d0) generation cut. Minimum kt in underlying Born
bornsuppfact 0d0 ! (default 0d0) mass param for Born suppression factor. If < 0 suppfact = 1

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

renscfact 1d0 ! (default 1d0) ren scale factor: muren = muref * renscfact
facscfact 1d0 ! (default 1d0) fac scale factor: mufact = muref * facscfact

! 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

! the typical call uses 1/1400 seconds (1400 calls per second)
ncall1 100000 ! Number of calls for the construction of the importance sampling grid
itmx1 1 ! Number of iterations for grid
ncall2 100000 ! Number of calls for the computation of the upper bounding envelope
! for the generation of radiation
itmx2 1 ! Number of iterations for the above
! Notice: the total number of calls is ncall2*itmx2*foldcsi*foldy*foldphi
foldcsi 2 ! number of folds on csi integration
foldy 5 ! number of folds on y integration
foldphi 2 ! number of folds on phi integration

nubound 100000 ! number of calls to set up the upper bounding norms for radiation
! This is performed using only the Born cross section (fast)


! OPTIONAL PARAMETERS

withnegweights 1 ! (default 0) if on (1) use negative weights
testplots 1 ! (default 0, do not) do NLO and PWHG distributions

iseed SEED

xgriditeration 1
fastbtlbound 1

minlo 1
minlo_nnll 1

# if running with minlo, set the following to 0
massivetop 0

sudscalevar 1

doublefsr 1

nohad 1

Original file line number Diff line number Diff line change
@@ -0,0 +1,47 @@
numevts NEVENTS
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

lhans1 325300 ! pdf set for hadron 1 (LHA numbering)
lhans2 325300 ! 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 125000 ! number of calls for initializing the integration grid
itmx1 7 ! number of iterations for initializing the integration grid
ncall2 75000 ! 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 2 ! number of folds on csi integration
foldy 5 ! number of folds on y integration
foldphi 2 ! number of folds on phi integration
nubound 175000 ! 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
withdamp 1 ! (default 0, do not use) use Born-zero damping factor
testplots 1 ! (default 0, do not) do NLO and PWHG distributions
#pdfreweight 0
#storeinfo_rwgt 0
#withnegweights 1 1 ! default 0,

iseed SEED

hmass 125d0 ! Higgs boson mass
hwidth 0.00407d0 ! Higgs boson width
hdecaymode -1 ! -1 no decay
! 0 all decay channels open
! 1-6 d dbar, u ubar,..., t tbar
! 7-9 e+ e-, mu+ mu-, tau+ tau-
! 10 W+W-
! 11 ZZ
! 12 gamma gamma

Original file line number Diff line number Diff line change
@@ -0,0 +1,70 @@
idvecbos -24 ! 24: W+ boson, -24: W- boson
vdecaymode 0 ! 1: e-nu, 2: mu-nu, 3: tau-nu, 0: hadronic, 10: inclusive, 11: inclusive leptons
hdecaymode -1 ! undecayed Higgs boson (for PYTHIA and HERWIG)

hmass 125d0 ! Higgs boson mass
hwidth 0.00407d0 ! Higgs boson width

min_h_mass 12d0
max_h_mass 1250d0

min_w_mass 10d0
max_w_mass 1000d0


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 6800d0 ! energy of beam 1
ebeam2 6800d0 ! energy of beam 2

bornktmin 0.26d0 ! (default 0d0) generation cut. Minimum kt in underlying Born
bornsuppfact 0d0 ! (default 0d0) mass param for Born suppression factor. If < 0 suppfact = 1

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

renscfact 1d0 ! (default 1d0) ren scale factor: muren = muref * renscfact
facscfact 1d0 ! (default 1d0) fac scale factor: mufact = muref * facscfact

! 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

! the typical call uses 1/1400 seconds (1400 calls per second)
ncall1 400000 ! Number of calls for the construction of the importance sampling grid
itmx1 1 ! Number of iterations for grid
ncall2 400000 ! Number of calls for the computation of the upper bounding envelope
! for the generation of radiation
itmx2 1 ! Number of iterations for the above
! Notice: the total number of calls is ncall2*itmx2*foldcsi*foldy*foldphi
foldcsi 2 ! number of folds on csi integration
foldy 5 ! number of folds on y integration
foldphi 2 ! number of folds on phi integration

nubound 600000 ! number of calls to set up the upper bounding norms for radiation
! This is performed using only the Born cross section (fast)


! OPTIONAL PARAMETERS

withnegweights 1 ! (default 0) if on (1) use negative weights
testplots 1 ! (default 0, do not) do NLO and PWHG distributions

iseed SEED ! Start the random number generator with seed iseed

fastbtlbound 1
storemintupb 1

minlo 1
minlo_nnll 1

# if running with minlo, set the following to 0
massivetop 0

sudscalevar 1

doublefsr 1

nohad 1

Loading

0 comments on commit 53f69ee

Please sign in to comment.