From 5cad836f18eff0ce9bbc63be10a4355fbc551c73 Mon Sep 17 00:00:00 2001 From: Tobias Neumann Date: Thu, 28 Mar 2024 11:22:14 -0700 Subject: [PATCH] Demo input file for 161/164 --- Bin/input_top.ini | 265 ++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 265 insertions(+) create mode 100644 Bin/input_top.ini diff --git a/Bin/input_top.ini b/Bin/input_top.ini new file mode 100644 index 0000000..d908107 --- /dev/null +++ b/Bin/input_top.ini @@ -0,0 +1,265 @@ +mcfm_version = 10.3 + +[general] + # process number + nproc = 164 + # part: lo, nlo, nnlo, nlocoeff, nnlocoeff + part = lo + # 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 = 30 + # maximum jet rapidity; can also have etajetmin + etajetmax = 2.4 + 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 = .false. + + # minimum lepton pT; can also have ptleptmax + ptleptmin = 0.0 + # maximum lepton rapidity, absolute value; can also have etaleptmin + etaleptmax = 99.0 + # 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 = 0.0 + # maximum rapidity for additional leptons, absolute value; can also have etalept2min + etalept2max = 99.0 + # 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 = .false. + +[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 = 2.5 + ndmx = 30 + iterbatchwarmup = 15 + initcallsnloreal = 8000000 + initcallsnlovirt = 2000000 + +# 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 = 1.0 + +[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 = .false. + nnlo_enable_interf_lxd = .false. + nnlo_enable_interf_hxd = .false. + nnlo_fully_inclusive = .false.