You can not select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
263 lines
6.5 KiB
263 lines
6.5 KiB
mcfm_version = 10.3
|
|
|
|
[general]
|
|
# process number
|
|
nproc = 1610
|
|
# part: lo, nlo, nnlo, nlocoeff, nnlocoeff
|
|
part = nnlo
|
|
# string to identify run
|
|
runstring = st_fiducial
|
|
rundir = st_fiducial
|
|
sqrts = 13000
|
|
# ih1, ih2: +1 for proton, -1 for antiproton
|
|
ih1 = +1
|
|
ih2 = +1
|
|
zerowidth = .false.
|
|
removebr = .false.
|
|
# electroweak corrections: none, sudakov or exact
|
|
ewcorr = none
|
|
|
|
[nnlo]
|
|
# optional: tau cutoff for NNLO processes
|
|
taucut = 0.001
|
|
|
|
[pdf]
|
|
# PDF label for internal routines
|
|
pdlabel = 'CT14.NN'
|
|
|
|
[lhapdf]
|
|
# LHAPDF PDF label, multiple sets separated by spaces possible
|
|
lhapdfset = CT14nnlo
|
|
# LHAPDF PDF member number, specify for each set, separated by spaces
|
|
lhapdfmember = 0
|
|
# calculate PDF uncertainties.
|
|
# note: can only do PDF uncertainties or scale variation at the same time.
|
|
dopdferrors = .false.
|
|
|
|
[scales]
|
|
;renscale = 173.3
|
|
;facscale = 173.3
|
|
;dynamicscale = none
|
|
renscale = 1.0
|
|
facscale = 1.0
|
|
dynamicscale = DDIS
|
|
# perform scale variation
|
|
doscalevar = .false.
|
|
# can be 2 or 6 for 2-point or 6-point scale-variation
|
|
maxscalevar = 6
|
|
|
|
[masses]
|
|
wmass = 80.385
|
|
# Higgs mass
|
|
hmass = 125
|
|
# Top-quark mass
|
|
mt = 173.3
|
|
# Bottom-quark mass
|
|
mb = 0.0
|
|
# Charm-quark mass
|
|
mc = 0.0
|
|
|
|
[basicjets]
|
|
# Jet-inclusive cross section
|
|
inclusive = .true.
|
|
# Jet-clustering algorithm: ankt, ktal, cone, hqrk, none
|
|
algorithm = ankt
|
|
# minimum jet pT; can also have ptjetmax
|
|
ptjetmin = 40
|
|
# maximum jet rapidity; can also have etajetmin
|
|
etajetmax = 5.0
|
|
Rcutjet = 0.5
|
|
|
|
[masscuts]
|
|
# minimum mass of 3-4 system
|
|
m34min = 0
|
|
# optional, maximum mass of 3-4 system, otherwise sqrts
|
|
# m34max = 14000
|
|
# minimum mass of 5-6 system
|
|
m56min = 0
|
|
# optional, maximum mass of 5-6 system, otherwise sqrts
|
|
# m56max = 14000
|
|
# minimum mass of 3-4-5-6 system
|
|
m3456min = 0
|
|
# optional, maximum mass of 3-4-5-6 system, otherwise sqrts
|
|
# m3456max = 14000
|
|
|
|
[cuts]
|
|
# make additional cuts in this section
|
|
# make sure to modify src/User/gencuts_user.f90 for additional b-tagging requirements
|
|
makecuts = .true.
|
|
|
|
# minimum lepton pT; can also have ptleptmax
|
|
ptleptmin = 30
|
|
# maximum lepton rapidity, absolute value; can also have etaleptmin
|
|
etaleptmax = 2.4
|
|
# rapidity veto on leptons
|
|
etaleptveto = 0.0 0.0
|
|
# minimum missing pT
|
|
ptminmiss = 0.0
|
|
|
|
# minimum pT for additional leptons; can also have ptlept2max
|
|
ptlept2min = 30.0
|
|
# maximum rapidity for additional leptons, absolute value; can also have etalept2min
|
|
etalept2max = 2.4
|
|
# rapidity veto on additional leptons
|
|
etalept2veto = 0.0 0.0
|
|
|
|
# minimum (3,4) transverse mass
|
|
m34transmin = 0.0
|
|
# minimum R(jet,lept)
|
|
Rjlmin = 0.0
|
|
# minimum R(lept,lept)
|
|
Rllmin = 0.0
|
|
# minimum rapidity separation between jets
|
|
delyjjmin = 0.0
|
|
# enforce jets to be in opposite hemispheres
|
|
jetsopphem = .false.
|
|
lbjscheme = 0
|
|
# b-jet minimum pT; can also have ptbjetmax
|
|
ptbjetmin = 0.0
|
|
# b-jet maximum rapidity; can also have etabjetmin
|
|
etabjetmax = 99.0
|
|
|
|
[photon]
|
|
# include fragmentation
|
|
fragmentation = .false.
|
|
# fragmentation set
|
|
fragmentation_set = GdRG__LO
|
|
# fragmentation scale
|
|
fragmentation_scale = 1.0
|
|
# minimum photon pT; can also have gammptmax
|
|
gammptmin = 40
|
|
# maximum photon rapidity; can also have gammrapmin
|
|
gammrapmax = 2.5
|
|
# second photon minimum pT
|
|
gammpt2 = 25
|
|
# third photon minimum pT
|
|
gammpt3 = 25
|
|
# minimum R(photon,lepton)
|
|
Rgalmin = 0
|
|
# minimum R(photon,photon)
|
|
Rgagamin = 0.4
|
|
# minimum R(photon,jet)
|
|
Rgajetmin = 0
|
|
# cone size for isolation
|
|
cone_ang = 0.4
|
|
# energy fraction for isolation
|
|
epsilon_h = 0.5
|
|
# exponent for smooth-cone isolation
|
|
n_pow = 1
|
|
|
|
[histogram]
|
|
# write top-drawer histograms
|
|
writetop = .false.
|
|
# write raw table file for each histogram
|
|
writetxt = .true.
|
|
newstyle = .true.
|
|
|
|
[integration]
|
|
# use sobol low discrepancy sequence
|
|
usesobol = .true.
|
|
# random number seed, when not using sobol, 0 means random
|
|
seed = 0
|
|
# relative precision goal
|
|
precisiongoal = 0.001
|
|
absprecisiongoal = 1.0
|
|
# resume from previous integration snapshot
|
|
readin = .false.
|
|
# write histograms and results after each vegas iteration
|
|
writeintermediate = .true.
|
|
# warmup precision goal in percent (default: 25%)
|
|
# unless this precision is reached, the number of calls will be
|
|
# increased and the warmup runs again
|
|
warmupprecisiongoal = 0.75
|
|
# warmup chisq goal (default: 2.5)
|
|
# unless this value of chisq/it is reached, the number of calls will be
|
|
# increased and the warmup runs again
|
|
warmupchisqgoal = 25.0
|
|
ndmx = 25
|
|
iterbatchwarmup = 10
|
|
|
|
# process specific initializations
|
|
|
|
# Single top SMEFT, nproc=164, 169
|
|
[singletop]
|
|
# C_phiq (O1), real-valued
|
|
c_phiq = 0
|
|
# C_phiphi (O2), real and imaginary part
|
|
c_phiphi = 0.0 0.0
|
|
# C_tW (O3), real and imaginary part
|
|
c_tw = 0.0 0.0
|
|
# C_bW (O4), real and imaginary part
|
|
c_bw = 0.0 0.0
|
|
# C_tG (O6), real and imaginary part
|
|
c_tg = 0.0 0.0
|
|
# C_bG (O7), real and imaginary part
|
|
c_bg = 0.0 0.0
|
|
# Lambda, scale of EFT breakdown in GeV
|
|
lambda = 1000
|
|
# enable 1/Lambda^4 contributions
|
|
enable_lambda4 = .false.
|
|
# disable Standard Model contributions
|
|
disable_sm = .false.
|
|
# anomalous couplings mode (only LO)
|
|
mode_anomcoup = .false.
|
|
|
|
# Anomalous couplings of the W and Z
|
|
[anom_wz]
|
|
# enable anomalous W/Z couplings
|
|
enable = .false.
|
|
# Delta g1(Z)
|
|
delg1_z = 0
|
|
# Delta K(Z)
|
|
delk_z = 0
|
|
# Delta K(gamma)
|
|
delk_g = 0
|
|
# Lambda(Z)
|
|
lambda_z = 0
|
|
# Lambda(gamma)
|
|
lambda_g = 0
|
|
# h1(Z)
|
|
h1Z = 0
|
|
# h1(gamma)
|
|
h1gam = 0
|
|
# h2(Z)
|
|
h2Z = 0
|
|
# h2(gamma)
|
|
h2gam = 0
|
|
# h3(Z)
|
|
h3Z = 0
|
|
# h3(gamma)
|
|
h3gam = 0
|
|
# h4(Z)
|
|
h4Z = 0
|
|
# h4(gamma)
|
|
h4gam = 0
|
|
# Form-factor scale, in TeV
|
|
tevscale = 2.0
|
|
|
|
# Higgs+jet with mass corrections, process 200
|
|
[hjetmass]
|
|
mtex = 0
|
|
|
|
[anom_higgs]
|
|
# Gamma_H / Gamma_H(SM)
|
|
hwidth_ratio = 1.0
|
|
cttH = 1.0
|
|
cWWH = 1.0
|
|
|
|
[dipoles]
|
|
alpha = 0.1
|
|
|
|
[extra]
|
|
verbose = .true.
|
|
nohistograms = .false.
|
|
|
|
# for NNLO
|
|
[singletop]
|
|
nnlo_enable_light = .true.
|
|
nnlo_enable_heavy_prod = .true.
|
|
nnlo_enable_heavy_decay = .true.
|
|
nnlo_enable_interf_lxh = .true.
|
|
nnlo_enable_interf_lxd = .true.
|
|
nnlo_enable_interf_hxd = .true.
|
|
nnlo_fully_inclusive = .false.
|