ATLAS Offline Software
WWj_MiNNLO.py
Go to the documentation of this file.
1 # Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
2 
3 from AthenaCommon import Logging
4 from ..powheg_RES import PowhegRES
5 
6 
7 logger = Logging.logging.getLogger("PowhegControl")
8 
9 
11  """! Default Powheg configuration for W-boson pair production plus one jet using MiNNLOPS.
12 
13  Create a configurable object with all applicable Powheg options.
14 
15  @author Aonan Wang <aonan.wang@cern.ch>
16  """
17 
18  def __init__(self, base_directory, **kwargs):
19  """! Constructor: all process options are set here.
20 
21  @param base_directory: path to PowhegBox code.
22  @param kwargs dictionary of arguments from Generate_tf.
23  """
24  errors = super(WWj_MiNNLO, self).openloops_error()
25  warnings = super(WWj_MiNNLO, self).hoppet_warning()
26  infos = super(WWj_MiNNLO, self).hoppet_info()
27  infos.append("qqvvamp: increasing precision to")
28  super(WWj_MiNNLO, self).__init__(base_directory, "WWJ",warning_output=warnings, info_output=infos, error_output=errors, **kwargs)
29 
30  # Add parameter validation functions
31  self.validation_functions.append("validate_decays")
32 
33  # Add flag for the MiNNLO reweight
34  self.reweight_for_MiNNLO = True
35 
36 
37  self.allowed_decay_modes = ["w+ w- > e+ ve e- ve~",
38  "w+ w- > mu+ vm mu- vm~",
39  "w+ w- > tau+ vt tau- vt~",
40  "w+ w- > e+ ve mu- vm~ / mu+ vm e- ve~",
41  "w+ w- > l+ vl l'- vl'~",
42  "w+ w- > l+ vl j j / j j l- vl~",
43  "w+ w- > e+ ve j j / j j e- ve~ / mu+ vmu j j / j j mu- vmu~",
44  "w+ w- > j j j j",
45  "w+ w- > e+ ve mu- vm~",
46  "w+ w- > mu+ vm e- ve~"]
47 
48  # Add all keywords for this process, overriding defaults if required
49  self.add_keyword("ih1")
50  self.add_keyword("ih2")
51  self.add_keyword("lhans1", self.default_PDFs_nnlo_nf_4)
52  self.add_keyword("lhans2", self.default_PDFs_nnlo_nf_4)
53  self.add_keyword("alphas_from_pdf", 1)
54  self.add_keyword("runningscales") # 0 = fixed scale 2m(W), 1=m(WW), 2=mT(W+) + mT(W-)
55  self.add_keyword("minlo", 1)
56  self.add_keyword("minnlo", 1)
57  self.add_keyword("modlog_p", 6)
58  self.add_keyword("Q0", 0)
59  self.add_keyword("largeptscales", 1)
60  self.add_keyword("smartMiNLO", 1)
61  self.add_keyword("use_interpolator", 1)
62  self.add_keyword("run_mode")
63  self.add_keyword("rwl_group_events", 1)
64  self.add_keyword("renscfact", self.default_scales[1])
65  self.add_keyword("facscfact", self.default_scales[0])
66  self.add_keyword("storeinfo_rwgt", 1)
67  self.add_keyword("rwl_file")
68  self.add_keyword("rwl_add")
69  self.add_keyword("rwl_format_rwgt")
70  self.add_keyword("clobberlhe")
71  self.add_keyword("ewscheme")
72  self.add_keyword("gfermi")
73  self.add_keyword("hmass")
74  self.add_keyword("zmass")
75  self.add_keyword("wmass")
76  self.add_keyword("tmass")
77  self.add_keyword("bmass")
78  self.add_keyword("zwidth")
79  self.add_keyword("wwidth")
80  self.add_keyword("twidth")
81  self.add_keyword("hwidth")
82  self.add_keyword("e+e-")
83  self.add_keyword("mu+mu-")
84  self.add_keyword("tau+tau-")
85  self.add_keyword("e+mu-")
86  self.add_keyword("mu+e-")
87  self.add_keyword("leptonic")
88  self.add_keyword("leptonic_notau")
89  self.add_keyword("hadronic")
90  self.add_keyword("semileptonic", "w+ w- > l+ vl l'- vl'~", name="decay_mode", hidden=False)
91  self.add_keyword("semileptonic_notau")
92  self.add_keyword("bornktmin")
93  self.add_keyword("use-old-grid")
94  self.add_keyword("ncall1", 30000)
95  self.add_keyword("itmx1", 1)
96  self.add_keyword("ncall2", 30000)
97  self.add_keyword("itmx2",1)
98  self.add_keyword("foldcsi")
99  self.add_keyword("foldy")
100  self.add_keyword("foldphi")
101  self.add_keyword("testplots")
102  self.add_keyword("use-old-ubound")
103  self.add_keyword("nubound", 50000)
104  self.add_keyword("storemintupb")
105  self.add_keyword("xupbound")
106  self.add_keyword("mintupbratlim")
107  self.add_keyword("ubexcess_correct")
108  self.add_keyword("ptsqmin")
109  # self.add_keyword("btildeviol")
110  # self.add_keyword("corr_btilde")
111  # self.add_keyword("corr_remnant")
112  self.add_keyword("colltest")
113  self.add_keyword("softtest")
114  self.add_keyword("withdamp")
115  self.add_keyword("smartsig")
116  self.add_keyword("fastbtlbound")
117  self.add_keyword("check_bad_st1")
118  self.add_keyword("check_bad_st2")
119  self.add_keyword("manyseeds")
120  self.add_keyword("parallelstage")
121  self.add_keyword("maxseeds")
122  self.add_keyword("xgriditeration")
123 
124  def validate_decays(self):
125  """! Validate the various decay mode keywords."""
126  self.expose() # convenience call to simplify syntax
127  self.check_decay_mode(self.decay_mode, self.allowed_decay_modes)
128  # Enable appropriate decay mode
129  list(self.parameters_by_keyword("semileptonic"))[0].value = 0 # disable the one used as a proxy for decay mode
130  __decay_mode_lookup = { "w+ w- > e+ ve e- ve~": "e+e-",
131  "w+ w- > mu+ vm mu- vm~": "mu+mu-",
132  "w+ w- > tau+ vt tau- vt~": "tau+tau-",
133  "w+ w- > e+ ve mu- vm~": "e+mu-",
134  "w+ w- > mu+ vm e- ve~": "mu+e-",
135  "w+ w- > l+ vl l'- vl'~": "leptonic",
136  "w+ w- > e+ ve mu- vm~ / mu+ vm e- ve~": "leptonic_notau",
137  "w+ w- > l+ vl j j / j j l- vl~": "semileptonic",
138  "w+ w- > e+ ve j j / j j e- ve~ / mu+ vmu j j / j j mu- vmu~": "semileptonic_notau",
139  "w+ w- > j j j j": "hadronic"}
140 
141  list(self.parameters_by_keyword(__decay_mode_lookup[self.decay_mode]))[0].value = 1
python.processes.configurable.Configurable.expose
def expose(self)
Add all names to the interface of this object.
Definition: configurable.py:46
python.processes.powheg.WWj_MiNNLO.WWj_MiNNLO.__init__
def __init__(self, base_directory, **kwargs)
Constructor: all process options are set here.
Definition: WWj_MiNNLO.py:18
python.processes.powheg_base.PowhegBase.hoppet_info
def hoppet_info(self)
Definition: powheg_base.py:23
python.processes.powheg_base.PowhegBase.hoppet_warning
def hoppet_warning(self)
Definition: powheg_base.py:44
python.processes.powheg_base.PowhegBase.default_scales
def default_scales(self)
Default scale variations for this process.
Definition: powheg_base.py:194
python.processes.configurable.Configurable.add_keyword
def add_keyword(self, keyword, value=None, name=None, frozen=None, hidden=None, description=None, **kwargs)
Register configurable parameter that is exposed to the user.
Definition: configurable.py:21
dumpHVPathFromNtuple.append
bool append
Definition: dumpHVPathFromNtuple.py:91
python.processes.powheg.WWj_MiNNLO.WWj_MiNNLO
Default Powheg configuration for W-boson pair production plus one jet using MiNNLOPS.
Definition: WWj_MiNNLO.py:10
python.processes.powheg_base.PowhegBase.openloops_error
def openloops_error(self)
Definition: powheg_base.py:51
python.processes.configurable.Configurable.parameters_by_keyword
def parameters_by_keyword(self, keyword)
Retrieve all parameters that use a given keyword.
Definition: configurable.py:57
python.processes.powheg_base.PowhegBase.validation_functions
validation_functions
List of validation functions to run before preparing runcard.
Definition: powheg_base.py:127
python.processes.powheg_base.PowhegBase.check_decay_mode
def check_decay_mode(self, decay_mode, allowed_decay_modes=None)
Check whether a decay mode is allowed an raise an exception if it is not.
Definition: powheg_base.py:231
histSizes.list
def list(name, path='/')
Definition: histSizes.py:38
python.processes.powheg.WWj_MiNNLO.WWj_MiNNLO.validate_decays
def validate_decays(self)
Validate the various decay mode keywords.
Definition: WWj_MiNNLO.py:124
python.processes.powheg_RES.PowhegRES.default_PDFs_nnlo_nf_4
def default_PDFs_nnlo_nf_4(self)
Default PDFs for this process.
Definition: powheg_RES.py:58
python.processes.powheg.WWj_MiNNLO.WWj_MiNNLO.reweight_for_MiNNLO
reweight_for_MiNNLO
Definition: WWj_MiNNLO.py:34
python.processes.powheg_RES.PowhegRES
Base class for PowhegBox RES processes.
Definition: powheg_RES.py:7
python.processes.powheg.WWj_MiNNLO.WWj_MiNNLO.allowed_decay_modes
allowed_decay_modes
List of allowed decay modes.
Definition: WWj_MiNNLO.py:37