ATLAS Offline Software
Z_SMEFT.py
Go to the documentation of this file.
1 # Copyright (C) 2002-2021 CERN for the benefit of the ATLAS collaboration
2 
3 from AthenaCommon import Logging
4 from ..powheg_V2 import PowhegV2
5 import itertools
6 
7 
8 logger = Logging.logging.getLogger("PowhegControl")
9 
10 
12  """!
13  Powheg configuration for Z boson production in Standard Model Effective Field Theory
14 
15  This process is an implementation of the Drell-Yan NLO cross sections pp -> Z/y* -> l+l- by SMEFT operators
16 
17  It includes 196 SMEFT coefficients that the user can set. For symmetry, the authors included coefficients
18  for operators involving top-quarks. These should not be used as they are only relevant when initial-state
19  tops are considered, which is never the case.
20 
21  Publication (please cite if using this process): https://arxiv.org/abs/1804.07407
22  Brief manual (highly recommended reading if using this process): accessible on LXPlus etc. at $POWHEGPATH/POWHEG-BOX-V2/Z_smeft/Docs/manual-BOX-Z_smeft.pdf
23 
24  @author Stefan Richter <stefan.richter@cern.ch>
25  """
26 
27  def __init__(self, base_directory, **kwargs):
28  """!
29  Constructor: all process options are set here.
30 
31  @param base_directory: path to PowhegBox code.
32  @param kwargs dictionary of arguments from Generate_tf.
33  """
34  super(Z_SMEFT, self).__init__(base_directory, "Z_smeft", **kwargs)
35 
36  # Add parameter validation functions
37  self.validation_functions.append("validate_decays")
38 
39 
40  self.allowed_decay_modes = ["z > e+ e-", "z > mu+ mu-", "z > tau+ tau-", "z > ve ve~", "z > vm vm~", "z > vt vt~"]
41 
42  # Add all keywords for this process, overriding defaults if required
43  self.add_keyword("alphaem")
44  self.add_keyword("bmass_lhe")
45  self.add_keyword("bornktmin")
46  self.add_keyword("bornonly")
47  self.add_keyword("bornsuppfact")
48  self.add_keyword("bornzerodamp")
49  self.add_keyword("bottomthr")
50  self.add_keyword("bottomthrpdf")
51  self.add_keyword("btlscalect")
52  self.add_keyword("btlscalereal")
53  self.add_keyword("charmthr")
54  self.add_keyword("charmthrpdf")
55  self.add_keyword("check_bad_st1")
56  self.add_keyword("check_bad_st2")
57  # clobberlhe: needed to make reweighting work --- otherwise Powheg crashes because it
58  # cannot overwrite the LHE file with events before reweighting.
59  # The value given to the parameter here doesn't matter, it just needs to be present in
60  # the powheg.input file. The reweighter postprocessor will replace the value with 1
61  # anyway.
62  # self.add_keyword("clobberlhe")
63  self.add_keyword("cmass_lhe")
64  self.add_keyword("colltest")
65  self.add_keyword("compute_rwgt")
66  self.add_keyword("dim6", 1)
67  self.add_keyword("facscfact", self.default_scales[0])
68  self.add_keyword("fcnc", 0)
69  self.add_keyword("flg_debug")
70  self.add_keyword("foldcsi", 5)
71  self.add_keyword("foldphi", 2)
72  self.add_keyword("foldy", 2)
73  self.add_keyword("fullrwgt")
74  self.add_keyword("fullrwgtmode")
75  self.add_keyword("hdamp")
76  self.add_keyword("hfact")
77  self.add_keyword("icsimax")
78  self.add_keyword("ih1")
79  self.add_keyword("ih2")
80  self.add_keyword("itmx1", 5)
81  self.add_keyword("itmx2", 2)
82  self.add_keyword("iupperfsr")
83  self.add_keyword("iupperisr")
84  self.add_keyword("iymax")
85  self.add_keyword("LambdaNP", 1000.0)
86  self.add_keyword("lhans1", self.default_PDFs)
87  self.add_keyword("lhans2", self.default_PDFs)
88  self.add_keyword("lhrwgt_descr")
89  self.add_keyword("lhrwgt_group_combine")
90  self.add_keyword("lhrwgt_group_name")
91  self.add_keyword("lhrwgt_id")
92  self.add_keyword("manyseeds")
93  self.add_keyword("mass_high", 2.0 * list(self.parameters_by_name("beam_energy"))[0].value)
94  self.add_keyword("mass_low", 60.0)
95  self.add_keyword("masswindow_high")
96  self.add_keyword("masswindow_low")
97  self.add_keyword("ncall1", 40000)
98  self.add_keyword("ncall2", 70000)
99  self.add_keyword("nubound", 10000)
100  self.add_keyword("par_2gsupp")
101  self.add_keyword("par_diexp")
102  self.add_keyword("par_dijexp")
103  self.add_keyword("parallelstage")
104  self.add_keyword("pdfreweight")
105  self.add_keyword("ptsqmin")
106  self.add_keyword("ptsupp")
107  self.add_keyword("radregion")
108  self.add_keyword("rand1")
109  self.add_keyword("rand2")
110  self.add_keyword("renscfact", self.default_scales[1])
111  self.add_keyword("running_width")
112  self.add_keyword("runningscale", 1) # if 1, central scale equal to Z/gamma virtuality. If 0, central scale equal to M_Z
113  self.add_keyword("rwl_add")
114  self.add_keyword("rwl_file")
115  self.add_keyword("rwl_format_rwgt")
116  self.add_keyword("rwl_group_events")
117  self.add_keyword("smartsig")
118  self.add_keyword("softtest")
119  self.add_keyword("sthw2")
120  self.add_keyword("storeinfo_rwgt")
121  self.add_keyword("testplots")
122  self.add_keyword("testsuda")
123  self.add_keyword("ubsigmadetails")
124  self.add_keyword("use-old-grid")
125  self.add_keyword("use-old-ubound")
126  self.add_keyword("vdecaymode", self.allowed_decay_modes[0], name="decay_mode")
127  self.add_keyword("withdamp")
128  self.add_keyword("withnegweights")
129  self.add_keyword("withsubtr")
130  self.add_keyword("Wmass")
131  self.add_keyword("Wwidth")
132  self.add_keyword("xupbound")
133  self.add_keyword("Zmass")
134  self.add_keyword("Zwidth")
135  # SMEFT operator coefficients
136  for coefficient_name in ["ReGEw", "ImGEw", "ReGEe", "ImGEe"]:
137  self.add_keyword(coefficient_name)
138  for coefficient in ["ReGUw_{flavours}", "ImGUw_{flavours}", "ReGUe_{flavours}", "ImGUe_{flavours}", "ReLeQu3_{flavours}", "ImLeQu3_{flavours}", "ReLeQu_{flavours}", "ImLeQu_{flavours}"]:
139  for flavours in ["{}{}".format(*x) for x in itertools.product(["u", "c", "t"], repeat=2)]:
140  coefficient_name = coefficient.format(flavours=flavours)
141  self.add_keyword(coefficient_name)
142  for coefficient in ["ReGDw_{flavours}", "ImGDw_{flavours}", "ReGDe_{flavours}", "ImGDe_{flavours}", "ReLedQ_{flavours}", "ImLedQ_{flavours}"]:
143  for flavours in ["{}{}".format(*x) for x in itertools.product(["d", "s", "b"], repeat=2)]:
144  coefficient_name = coefficient.format(flavours=flavours)
145  self.add_keyword(coefficient_name)
146  for coefficient in ["QphiU_{flavours}", "Uphi_{flavours}", "QLu_{flavours}", "Ceu_{flavours}", "CLu_{flavours}"]:
147  for flavours in ["uu", "uc", "ut", "cc", "ct", "tt"]:
148  coefficient_name = coefficient.format(flavours=flavours)
149  self.add_keyword(coefficient_name)
150  for coefficient in ["QphiD_{flavours}", "Dphi_{flavours}", "QLd_{flavours}", "Ced_{flavours}", "CLd_{flavours}", "Qe_{flavours}"]:
151  for flavours in ["dd", "ds", "db", "ss", "sb", "bb"]:
152  coefficient_name = coefficient.format(flavours=flavours)
153  self.add_keyword(coefficient_name)
154 
155 
156 
157  def validate_decays(self):
158  """! Validate vdecaymode keyword."""
159  self.expose() # convenience call to simplify syntax
160  self.check_decay_mode(self.decay_mode, self.allowed_decay_modes)
161  # Calculate appropriate decay mode numbers
162  __decay_mode_lookup = {"e+ e-": 1, "mu+ mu-": 2, "tau+ tau-": 3, "ve ve~": 4, "vm vm~": 5, "vt vt~": 6}
163  list(self.parameters_by_keyword("vdecaymode"))[0].value = __decay_mode_lookup[self.decay_mode.split("> ")[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.Z_SMEFT.Z_SMEFT.__init__
def __init__(self, base_directory, **kwargs)
Constructor: all process options are set here.
Definition: Z_SMEFT.py:27
vtune_athena.format
format
Definition: vtune_athena.py:14
python.processes.configurable.Configurable.parameters_by_name
def parameters_by_name(self, name)
Retrieve all parameters that use a given name.
Definition: configurable.py:64
python.processes.powheg.Z_SMEFT.Z_SMEFT.allowed_decay_modes
allowed_decay_modes
List of allowed decay modes.
Definition: Z_SMEFT.py:40
python.processes.powheg_base.PowhegBase.default_scales
def default_scales(self)
Default scale variations for this process.
Definition: powheg_base.py:246
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.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:179
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:283
histSizes.list
def list(name, path='/')
Definition: histSizes.py:38
python.processes.powheg.Z_SMEFT.Z_SMEFT
Powheg configuration for Z boson production in Standard Model Effective Field Theory.
Definition: Z_SMEFT.py:11
python.processes.powheg_base.PowhegBase.default_PDFs
def default_PDFs(self)
Default PDFs for this process.
Definition: powheg_base.py:241
python.processes.powheg_V2.PowhegV2
Base class for PowhegBox V2 processes.
Definition: powheg_V2.py:6
python.processes.powheg.Z_SMEFT.Z_SMEFT.validate_decays
def validate_decays(self)
Validate vdecaymode keyword.
Definition: Z_SMEFT.py:157
Trk::split
@ split
Definition: LayerMaterialProperties.h:38