Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

ggH with SMEFTatNLO and get_scaling refactor #26

Merged
merged 3 commits into from
Jun 2, 2022
Merged

ggH with SMEFTatNLO and get_scaling refactor #26

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
bbd0d07
Adding 1jet cards for ggH and some basic docs
ajgilbert Oct 20, 2021
e1618f4
Adding EFT2ObsHist class that takes over functionality previously in …
ajgilbert Feb 15, 2022
101830d
Finish refactor: EFT2ObsHist as a bare representation of the set of r…
ajgilbert Jun 2, 2022
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
121 changes: 121 additions & 0 deletions cards/ggF-SMEFTatNLO-1jet-loop/param_card.dat
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,121 @@
######################################################################
## PARAM_CARD AUTOMATICALY GENERATED BY MG5 ####
######################################################################
###################################
## INFORMATION FOR DIM6
###################################
BLOCK DIM6 #
1 1.000000e+03 # lambda
8 1.000000e-08 # cpg
###################################
## INFORMATION FOR DIM62F
###################################
BLOCK DIM62F #
19 1.000000e-08 # ctp
24 1.000000e-08 # ctg
###################################
## INFORMATION FOR LOOP
###################################
BLOCK LOOP #
1 1.250000e+02 # mu_r
###################################
## INFORMATION FOR MASS
###################################
BLOCK MASS #
6 1.720000e+02 # mt
23 9.118760e+01 # mz
24 7.982440e+01 # mw
25 1.250000e+02 # mh
1 0.000000e+00 # d : 0.0
2 0.000000e+00 # u : 0.0
3 0.000000e+00 # s : 0.0
4 0.000000e+00 # c : 0.0
5 0.000000e+00 # b : 0.0
11 0.000000e+00 # e- : 0.0
12 0.000000e+00 # ve : 0.0
13 0.000000e+00 # mu- : 0.0
14 0.000000e+00 # vm : 0.0
15 0.000000e+00 # ta- : 0.0
16 0.000000e+00 # vt : 0.0
21 0.000000e+00 # g : 0.0
22 0.000000e+00 # a : 0.0
9000002 9.118760e+01 # ghz : mz
9000003 7.982440e+01 # ghwp : mw
9000004 7.982440e+01 # ghwm : mw
###################################
## INFORMATION FOR RENOR
###################################
BLOCK RENOR #
1 9.118760e+01 # mueft
###################################
## INFORMATION FOR SMINPUTS
###################################
BLOCK SMINPUTS #
2 1.166370e-05 # gf
3 1.184000e-01 # as
###################################
## INFORMATION FOR YUKAWA
###################################
BLOCK YUKAWA #
6 1.720000e+02 # ymt
###################################
## INFORMATION FOR DECAY
###################################
DECAY 1 0.000000e+00 # d : 0.0
DECAY 2 0.000000e+00 # u : 0.0
DECAY 3 0.000000e+00 # s : 0.0
DECAY 4 0.000000e+00 # c : 0.0
DECAY 5 0.000000e+00 # b : 0.0
DECAY 6 0.000000e+00 # t : 0.0
DECAY 11 0.000000e+00 # e- : 0.0
DECAY 12 0.000000e+00 # ve : 0.0
DECAY 13 0.000000e+00 # mu- : 0.0
DECAY 14 0.000000e+00 # vm : 0.0
DECAY 15 0.000000e+00 # ta- : 0.0
DECAY 16 0.000000e+00 # vt : 0.0
DECAY 21 0.000000e+00 # g : 0.0
DECAY 22 0.000000e+00 # a : 0.0
DECAY 23 0.000000e+00 # z : 0.0
DECAY 24 0.000000e+00 # w+ : 0.0
DECAY 25 0.000000e+00 # h : 0.0
###################################
## INFORMATION FOR QNUMBERS 9000001
###################################
BLOCK QNUMBERS 9000001 # gha
1 0 # 3 times electric charge
2 1 # number of spin states (2s+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # particle/antiparticle distinction (0=own anti)
###################################
## INFORMATION FOR QNUMBERS 9000002
###################################
BLOCK QNUMBERS 9000002 # ghz
1 0 # 3 times electric charge
2 1 # number of spin states (2s+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # particle/antiparticle distinction (0=own anti)
###################################
## INFORMATION FOR QNUMBERS 9000003
###################################
BLOCK QNUMBERS 9000003 # ghwp
1 3 # 3 times electric charge
2 1 # number of spin states (2s+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # particle/antiparticle distinction (0=own anti)
###################################
## INFORMATION FOR QNUMBERS 9000004
###################################
BLOCK QNUMBERS 9000004 # ghwm
1 -3 # 3 times electric charge
2 1 # number of spin states (2s+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # particle/antiparticle distinction (0=own anti)
###################################
## INFORMATION FOR QNUMBERS 9000005
###################################
BLOCK QNUMBERS 9000005 # ghg
1 0 # 3 times electric charge
2 1 # number of spin states (2s+1)
3 8 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # particle/antiparticle distinction (0=own anti)

6 changes: 6 additions & 0 deletions cards/ggF-SMEFTatNLO-1jet-loop/proc_card.dat
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,6 @@
import model SMEFTatNLO-ggFshape

generate p p > H QED=1 NP=0 [noborn=QCD]
add process p p > H j QED=1 NP=0 [noborn=QCD]

output ggF-SMEFTatNLO-1jet-loop
87 changes: 87 additions & 0 deletions cards/ggF-SMEFTatNLO-1jet-loop/pythia8_card.dat
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,87 @@
!
! Pythia8 cmd card automatically generated by MadGraph5_aMC@NLO
! For more information on the use of the MG5aMC / Pythia8 interface, visit
! https://cp3.irmp.ucl.ac.be/projects/madgraph/wiki/LOPY8Merging
!
! ==================
! General parameters
! ==================
!
Main:numberOfEvents = -1
!
! -------------------------------------------------------------------
! Specify the HEPMC output of the Pythia8 shower. You can set it to:
! auto : MG5aMC will automatically place it the run_<i> directory
! /dev/null : to turn off the HEPMC output.
! <path> : to select where the HEPMC file must written. It will
! therefore not be placed in the run_<i> directory. The
! specified path, if not absolute, will be relative to
! the Event/run_<i> directory of the process output.
! fifo : to have MG5aMC setup the piping of the PY8 output to
! analysis tools such as MadAnalysis5.
! fifo@<fifo_path> :
! Same as 'fifo', but selecting a custom path to create the
! fifo pipe. (useful to select a mounted drive that supports
! fifo). Note that the fifo file extension *must* be '.hepmc.fifo'.
! -------------------------------------------------------------------
!
HEPMCoutput:file = auto
!
! --------------------------------------------------------------------
! Parameters relevant only when performing MLM merging, which can be
! turned on by setting ickkw to '1' in the run_card and chosing a
! positive value for the parameter xqcut.
! For details, see section 'Jet Matching' on the left-hand menu of
! http://home.thep.lu.se/~torbjorn/pythia81html/Welcome.html
! --------------------------------------------------------------------
! If equal to -1.0, MadGraph5_aMC@NLO will set it automatically based
! on the parameter 'xqcut' of the run_card.dat
JetMatching:qCut = -1.0
! Use default kt-MLM to match parton level jets to those produced by the
! shower. But the other Shower-kt scheme is available too with this option.
JetMatching:doShowerKt = off
! A value of -1 means that it is automatically guessed by MadGraph.
! It is however always safer to explicitly set it.
JetMatching:nJetMax = -1
!
! --------------------------------------------------------------------
! Parameters relevant only when performing CKKW-L merging, which can
! be turned on by setting the parameter 'ptlund' *or* 'ktdurham' to
! a positive value.
! For details, see section 'CKKW-L Merging' on the left-hand menu of
! http://home.thep.lu.se/~torbjorn/pythia81html/Welcome.html
! --------------------------------------------------------------------
! Central merging scale values you want to be used.
! If equal to -1.0, then MadGraph5_aMC@NLO will set this automatically
! based on the parameter 'ktdurham' of the run_card.dat
Merging:TMS = -1.0
! This must be set manually, according to Pythia8 directives.
! An example of possible value is 'pp>LEPTONS,NEUTRINOS'
! Alternatively, from Pythia v8.223 onwards, the value 'guess' can be
! used to instruct Pythia to guess the hard process. The guess would mean
! that all particles apart from light partons will be considered as a part
! of the hard process. This guess is prone to errors if the desired hard
! process is complicated (i.e. contains light partons). The user should
! then be wary of suspicious error messages in the Pythia log file.
Merging:Process = <set_by_user>
! A value of -1 means that it is automatically guessed by MadGraph.
! It is however always safer to explicitly set it.
Merging:nJetMax = -1
!
! For all merging schemes, decide whehter you want the merging scale
! variation computed for only the central weights or all other
! PDF and scale variation weights as well
SysCalc:fullCutVariation = off
!
! ==========================
! User customized parameters
! ==========================
!
! By default, Pythia8 generates multi-parton interaction events. This is
! often irrelevant for phenomenology and very slow. You can turn this
! feature off by uncommenting the line below if so desired.
partonlevel:mpi = off
SLHA:keepSM = on
SLHA:minMassSM = 1000.
25:onMode = off
25:onIfMatch = 22 22
43 changes: 43 additions & 0 deletions cards/ggF-SMEFTatNLO-1jet-loop/reweight_card.dat
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,43 @@
change rwgt_dir rwgt
change process p p > h NP=2 QED=1 [noborn=QCD]
change process p p > h j NP=2 QED=1 [noborn=QCD] --add
launch --rwgt_name=rw0000
set dim6 8 0
set dim62f 19 0
set dim62f 24 0
launch --rwgt_name=rw0001
set dim6 8 0.005
set dim62f 19 0
set dim62f 24 0
launch --rwgt_name=rw0002
set dim6 8 0.01
set dim62f 19 0
set dim62f 24 0
launch --rwgt_name=rw0003
set dim6 8 0
set dim62f 19 0.005
set dim62f 24 0
launch --rwgt_name=rw0004
set dim6 8 0
set dim62f 19 0.01
set dim62f 24 0
launch --rwgt_name=rw0005
set dim6 8 0
set dim62f 19 0
set dim62f 24 0.005
launch --rwgt_name=rw0006
set dim6 8 0
set dim62f 19 0
set dim62f 24 0.01
launch --rwgt_name=rw0007
set dim6 8 0.01
set dim62f 19 0.01
set dim62f 24 0
launch --rwgt_name=rw0008
set dim6 8 0.01
set dim62f 19 0
set dim62f 24 0.01
launch --rwgt_name=rw0009
set dim6 8 0
set dim62f 19 0.01
set dim62f 24 0.01
Loading