ATLAS Offline Software
Zgamj_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 
10 class Zgamj_MiNNLO(PowhegRES):
11  """! Default Powheg configuration for Z gamma 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(Zgamj_MiNNLO, self).openloops_error()
25  warnings = super(Zgamj_MiNNLO, self).hoppet_warning()
26  infos = super(Zgamj_MiNNLO, self).hoppet_info()
27  infos.append("qqvvamp: increasing precision to")
28  super(Zgamj_MiNNLO, self).__init__(base_directory, "ZgamJ",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 = ["z > e+ e-",
38  "z > mu+ mu-",
39  "z > tau+ tau-",
40  "z > e+ e- / mu+ mu-",
41  "z > l+ l-",
42  "z > vl vl~"]
43 
44  # Add all keywords for this process, overriding defaults if required
45  self.add_keyword("ih1")
46  self.add_keyword("ih2")
47  self.add_keyword("lhans1", self.default_PDFs_nnlo)
48  self.add_keyword("lhans2", self.default_PDFs_nnlo)
49  self.add_keyword("alphas_from_pdf", 1)
50  self.add_keyword("fixedscale", 0) # if 0 use dynamical scale below (set by whichscale), if 1 scale is fixed to Z-boson mass (leave this to 0 when MiN(N)LO is used)
51  self.add_keyword("whichscale", 1) # only used if fixedscale=0, 0 = sqrt(M_Z^2+pt_gamma^2), 1 = M_llgamma, (leave this to 1 when MiN(N)LO is used)
52  self.add_keyword("minlo", 1)
53  self.add_keyword("minnlo", 1)
54  self.add_keyword("modlog_p", -1)
55  self.add_keyword("Q0", 0)
56  self.add_keyword("largeptscales", 1)
57  self.add_keyword("kappaQ", 1)
58  self.add_keyword("smartMiNLO", 1)
59  self.add_keyword("run_mode")
60  self.add_keyword("rwl_group_events", 1)
61  self.add_keyword("renscfact", self.default_scales[1])
62  self.add_keyword("facscfact", self.default_scales[0])
63  self.add_keyword("storeinfo_rwgt", 1)
64  self.add_keyword("rwl_file")
65  self.add_keyword("rwl_add")
66  self.add_keyword("rwl_format_rwgt")
67  self.add_keyword("clobberlhe")
68  self.add_keyword("ewscheme")
69  self.add_keyword("gfermi")
70  self.add_keyword("zmass")
71  self.add_keyword("wmass")
72  self.add_keyword("tmass")
73  self.add_keyword("zwidth")
74  self.add_keyword("wwidth")
75  self.add_keyword("twidth")
76  self.add_keyword("vdecaymode", "w+ w- > l+ vl l'- vl'~", name="decay_mode", hidden=False)
77  self.add_keyword("e+e-")
78  self.add_keyword("mu+mu-")
79  self.add_keyword("tau+tau-")
80  self.add_keyword("leptonic_notau")
81  self.add_keyword("leptonic")
82  self.add_keyword("sum_over_families")
83  self.add_keyword("massive_leptons")
84  self.add_keyword("e_mass")
85  self.add_keyword("mu_mass")
86  self.add_keyword("tau_mass")
87  self.add_keyword("useOL")
88  self.add_keyword("OL_CMS")
89  self.add_keyword("OL_nf")
90  self.add_keyword("OL_onshellphoton")
91  self.add_keyword("anomcoup")
92  self.add_keyword("hZ1")
93  self.add_keyword("hZ2")
94  self.add_keyword("hZ3")
95  self.add_keyword("hZ4")
96  self.add_keyword("hg1")
97  self.add_keyword("hg2")
98  self.add_keyword("hg3")
99  self.add_keyword("hg4")
100  self.add_keyword("anommode")
101  self.add_keyword("pt_j1_cut")
102  self.add_keyword("pt_a_cut")
103  self.add_keyword("m_lepg_cut")
104  self.add_keyword("invmass_min")
105  self.add_keyword("smooth_dyn")
106  self.add_keyword("smooth_eps")
107  self.add_keyword("smooth_pt0")
108  self.add_keyword("smooth_R")
109  self.add_keyword("smooth_n")
110  self.add_keyword("withdamp", 1)
111  self.add_keyword("suppmodel")
112  self.add_keyword("bornsuppfact", 1)
113  self.add_keyword("remnsuppfact")
114  self.add_keyword("ptj_suppfact")
115  self.add_keyword("powj_suppfact")
116  self.add_keyword("pta_suppfact")
117  self.add_keyword("powa_suppfact")
118  self.add_keyword("ptnunu_suppfact")
119  self.add_keyword("DRal_suppfact")
120  self.add_keyword("DRj_suppfact")
121  self.add_keyword("powdr_suppfact")
122  self.add_keyword("use-old-grid")
123  self.add_keyword("ncall1", 500000)
124  self.add_keyword("itmx1", 1)
125  self.add_keyword("ncall2", 500000)
126  self.add_keyword("itmx2",1)
127  self.add_keyword("foldcsi")
128  self.add_keyword("foldy")
129  self.add_keyword("foldphi")
130  self.add_keyword("testplots")
131  self.add_keyword("use-old-ubound")
132  self.add_keyword("nubound", 50000)
133  self.add_keyword("storemintupb")
134  self.add_keyword("xupbound", 2)
135  self.add_keyword("mintupbratlim", 1000)
136  self.add_keyword("ubexcess_correct", 1)
137  self.add_keyword("iymax", 3)
138  self.add_keyword("icsimax", 3)
139  # self.add_keyword("btildeviol")
140  # self.add_keyword("corr_btilde")
141  # self.add_keyword("corr_remnant")
142  self.add_keyword("colltest", 0)
143  self.add_keyword("softtest", 0)
144  self.add_keyword("smartsig")
145  self.add_keyword("fastbtlbound")
146  self.add_keyword("check_bad_st1")
147  self.add_keyword("check_bad_st2")
148  self.add_keyword("manyseeds")
149  self.add_keyword("parallelstage")
150  self.add_keyword("xgriditeration")
151  self.add_keyword("maxseeds")
152 
153  def validate_decays(self):
154  """! Validate the various decay mode keywords."""
155  self.expose() # convenience call to simplify syntax
156  self.check_decay_mode(self.decay_mode, self.allowed_decay_modes)
157  # Enable appropriate decay mode
158  list(self.parameters_by_keyword("vdecaymode"))[0].value = 0
159  __decay_mode_lookup = { "z > e+ e-": "e+e-",
160  "z > mu+ mu-": "mu+mu-",
161  "z > tau+ tau-": "tau+tau-",
162  "z > e+ e- / mu+ mu-": "leptonic_notau",
163  "z > l+ l-": "leptonic",
164  "z > vl vl~": "vv"}
165  if __decay_mode_lookup[self.decay_mode] == "vv":
166  list(self.parameters_by_keyword("vdecaymode"))[0].value = 0
167  else:
168  list(self.parameters_by_keyword("vdecaymode"))[0].value = 1
169  list(self.parameters_by_keyword(__decay_mode_lookup[self.decay_mode]))[0].value = 1
170  list(self.parameters_by_keyword("anomcoup"))[0].value = 0
171 
172 
python.processes.configurable.Configurable.expose
def expose(self)
Add all names to the interface of this object.
Definition: configurable.py:46
python.processes.powheg_base.PowhegBase.hoppet_info
def hoppet_info(self)
Definition: powheg_base.py:23
python.processes.powheg.Zgamj_MiNNLO.Zgamj_MiNNLO.reweight_for_MiNNLO
reweight_for_MiNNLO
Definition: Zgamj_MiNNLO.py:34
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
python.processes.powheg.Zgamj_MiNNLO.Zgamj_MiNNLO.__init__
def __init__(self, base_directory, **kwargs)
Constructor: all process options are set here.
Definition: Zgamj_MiNNLO.py:18
dumpHVPathFromNtuple.append
bool append
Definition: dumpHVPathFromNtuple.py:91
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_RES.PowhegRES.default_PDFs_nnlo
def default_PDFs_nnlo(self)
Default PDFs for this process.
Definition: powheg_RES.py:43
python.processes.powheg.Zgamj_MiNNLO.Zgamj_MiNNLO.allowed_decay_modes
allowed_decay_modes
List of allowed decay modes.
Definition: Zgamj_MiNNLO.py:37
python.processes.powheg.Zgamj_MiNNLO.Zgamj_MiNNLO
Default Powheg configuration for Z gamma production plus one jet using MiNNLOPS.
Definition: Zgamj_MiNNLO.py:10
python.processes.powheg.Zgamj_MiNNLO.Zgamj_MiNNLO.validate_decays
def validate_decays(self)
Validate the various decay mode keywords.
Definition: Zgamj_MiNNLO.py:153
python.processes.powheg_RES.PowhegRES
Base class for PowhegBox RES processes.
Definition: powheg_RES.py:7
Generated on Thu Sep 25 2025 21:22:31 for ATLAS Offline Software by   doxygen 1.8.18