File: StandardModelParameters.php

package info (click to toggle)
pythia8 8.1.65-1
  • links: PTS, VCS
  • area: main
  • in suites: wheezy
  • size: 22,660 kB
  • sloc: cpp: 59,593; xml: 30,509; php: 6,649; sh: 796; makefile: 353; ansic: 33
file content (431 lines) | stat: -rw-r--r-- 17,979 bytes parent folder | download
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
<html>
<head>
<title>Standard-Model Parameters</title>
<link rel="stylesheet" type="text/css" href="pythia.css"/>
<link rel="shortcut icon" href="pythia32.gif"/>
</head>
<body>

<script language=javascript type=text/javascript>
function stopRKey(evt) {
var evt = (evt) ? evt : ((event) ? event : null);
var node = (evt.target) ? evt.target :((evt.srcElement) ? evt.srcElement : null);
if ((evt.keyCode == 13) && (node.type=="text"))
{return false;}
}

document.onkeypress = stopRKey;
</script>
<?php
if($_POST['saved'] == 1) {
if($_POST['filepath'] != "files/") {
echo "<font color='red'>SETTINGS SAVED TO FILE</font><br/><br/>"; }
else {
echo "<font color='red'>NO FILE SELECTED YET.. PLEASE DO SO </font><a href='SaveSettings.php'>HERE</a><br/><br/>"; }
}
?>

<form method='post' action='StandardModelParameters.php'>

<h2>Standard-Model Parameters</h2>

<h3>The strong coupling</h3> 

The <code>AlphaStrong</code> class is used to provide a first- or 
second-order running <i>alpha_strong</i> (or, trivially, a 
zeroth-order fixed one). Formulae are the standard ones found in 
[<a href="Bibliography.php" target="page">Yao06</a>]. The second-order expression used, eq. (9.5),
may be somewhat different in other approaches (with differences
formally of higher order), so do not necessarily expect perfect
agreement, especially not at small <i>Q^2</i> scales. The starting 
<i>alpha_strong</i> value is defined at the <i>M_Z</i> mass scale.
The <i>Lambda</i> values are matched at the <i>b</i> and <i>c</i> 
flavour thresholds, such that <i>alpha_strong</i> is continuous.
For second-order matching an approximate iterative method is used.
 
<p/>
Since we allow <i>alpha_strong</i> to vary separately for 
hard processes, timelike showers, spacelike showers and  multiparton 
interactions, the relevant values can be set in each of these classes. 
The default behaviour is everywhere first-order running.
 
<p/>
The <i>alpha_strong</i> calculation is initialized by 
<code>init( value, order)</code>, where <code>value</code> 
is the <i>alpha_strong</i> value at <i>M_Z</i> and <code>order</code> 
is the order of the running, 0, 1 or 2.   Thereafter the value can be 
calculated by <code>alphaS(scale2)</code>, where 
<code>scale2</code> is the <i>Q^2</i> scale in GeV^2. 

<p/>
For applications inside shower programs, a second-order <code>alpha_s</code> 
value can be obtained as the product of the two functions 
<code>alphaS1Ord(scale2)</code> and <code>alphaS2OrdCorr(scale2)</code>, 
where the first gives a simple first-order running (but with the 
second-order <i>Lambda</i>) and the second the correction factor, 
below unity, for the second-order terms. This allows a compact handling 
of evolution equations.

<h3>The electromagnetic coupling</h3> 

The <code>AlphaEM</code> class is used to generate a running
<i>alpha_em</i>. The input <code>StandardModel:alphaEMmZ</code>
value at the <i>M_Z</i> mass is matched to a low-energy behaviour
with running starting at the electron mass threshold. The matching
is done by fitting an effective running coefficient in the region
betweeen the light-quark treshold and the charm/tau threshold. This
procedure is approximate, but good enough for our purposes. 

<p/>
Since we allow <i>alpha_em</i> to vary separately for 
hard processes, timelike showers, spacelike showers and  multiparton 
interactions, the choice between using a fixed or a running 
<i>alpha_em</i> can be made in each of these classes. 
The default behaviour is everywhere first-order running.
The actual values assumed at zero momentum transfer and 
at <i>M_Z</i> are only set here, however. 

<br/><br/><table><tr><td><strong>StandardModel:alphaEM0 </td><td></td><td> <input type="text" name="1" value="0.00729735" size="20"/>  &nbsp;&nbsp;(<code>default = <strong>0.00729735</strong></code>; <code>minimum = 0.0072973</code>; <code>maximum = 0.0072974</code>)</td></tr></table>
The <i>alpha_em</i> value at vanishing momentum transfer
(and also below <i>m_e</i>). 
  

<br/><br/><table><tr><td><strong>StandardModel:alphaEMmZ </td><td></td><td> <input type="text" name="2" value="0.00781751" size="20"/>  &nbsp;&nbsp;(<code>default = <strong>0.00781751</strong></code>; <code>minimum = 0.00780</code>; <code>maximum = 0.00783</code>)</td></tr></table>
The <i>alpha_em</i> value at the <i>M_Z</i> mass scale. 
Default is taken from [<a href="Bibliography.php" target="page">Yao06</a>].
  

<p/>
The <i>alpha_em</i> calculation is initialized by 
<code>init(order)</code>, where <code>order</code> is the order of 
the running, 0 or 1, with -1 a special option to use the fix value
provided at <i>M_Z</i>.   Thereafter the value can be 
calculated by <code>alphaEM(scale2)</code>, where 
<code>scale2</code> is the <i>Q^2</i> scale in GeV^2. 

<h3>The electroweak couplings</h3> 

There are two degrees of freedom that can be set, related to the 
electroweak mixing angle:

<br/><br/><table><tr><td><strong>StandardModel:sin2thetaW </td><td></td><td> <input type="text" name="3" value="0.2312" size="20"/>  &nbsp;&nbsp;(<code>default = <strong>0.2312</strong></code>; <code>minimum = 0.225</code>; <code>maximum = 0.240</code>)</td></tr></table>
The sine-squared of the weak mixing angle, as used in all <i>Z^0</i> 
and <i>W^+-</i> masses and couplings, except for the vector couplings 
of fermions to the <i>Z^0</i>, see below. Default is the MSbar value 
from [<a href="Bibliography.php" target="page">Yao06</a>].
  

<br/><br/><table><tr><td><strong>StandardModel:sin2thetaWbar </td><td></td><td> <input type="text" name="4" value="0.2315" size="20"/>  &nbsp;&nbsp;(<code>default = <strong>0.2315</strong></code>; <code>minimum = 0.225</code>; <code>maximum = 0.240</code>)</td></tr></table>
The sine-squared of the weak mixing angle, as used to derive the vector 
couplings of fermions to the <i>Z^0</i>, in the relation 
<i>v_f = a_f - 4 e_f sin^2(theta_W)bar</i>. Default is the
effective-angle value from [<a href="Bibliography.php" target="page">Yao06</a>].
  

<p/>
The Fermi constant is not much used in the currently coded matrix elements,
since it is redundant, but it is available:

<br/><br/><table><tr><td><strong>StandardModel:GF </td><td></td><td> <input type="text" name="5" value="1.16637e-5" size="20"/>  &nbsp;&nbsp;(<code>default = <strong>1.16637e-5</strong></code>; <code>minimum = 1.0e-5</code>; <code>maximum = 1.3e-5</code>)</td></tr></table>
The Fermi coupling constant, in units of GeV<i>^-2</i>. 
  

<h3>The quark weak-mixing matrix</h3>

The absolute values of the Cabibbo-Kobayashi-Maskawa matrix elements are 
set by the following nine real values taken from [<a href="Bibliography.php" target="page">Yao06</a>] - 
currently the CP-violating phase is not taken into account in this 
parametrization. It is up to the user to pick a consistent unitary 
set of new values whenever changes are made.  

<br/><br/><table><tr><td><strong>StandardModel:Vud </td><td></td><td> <input type="text" name="6" value="0.97383" size="20"/>  &nbsp;&nbsp;(<code>default = <strong>0.97383</strong></code>; <code>minimum = 0.973</code>; <code>maximum = 0.975</code>)</td></tr></table>
The <i>V_ud</i> CKM matrix element.
  

<br/><br/><table><tr><td><strong>StandardModel:Vus </td><td></td><td> <input type="text" name="7" value="0.2272" size="20"/>  &nbsp;&nbsp;(<code>default = <strong>0.2272</strong></code>; <code>minimum = 0.224</code>; <code>maximum = 0.230</code>)</td></tr></table>
The <i>V_us</i> CKM matrix element.
  

<br/><br/><table><tr><td><strong>StandardModel:Vub </td><td></td><td> <input type="text" name="8" value="0.00396" size="20"/>  &nbsp;&nbsp;(<code>default = <strong>0.00396</strong></code>; <code>minimum = 0.0037</code>; <code>maximum = 0.0042</code>)</td></tr></table>
The <i>V_ub</i> CKM matrix element.
  

<br/><br/><table><tr><td><strong>StandardModel:Vcd </td><td></td><td> <input type="text" name="9" value="0.2271" size="20"/>  &nbsp;&nbsp;(<code>default = <strong>0.2271</strong></code>; <code>minimum = 0.224</code>; <code>maximum = 0.230</code>)</td></tr></table>
The <i>V_cd</i> CKM matrix element.
  

<br/><br/><table><tr><td><strong>StandardModel:Vcs </td><td></td><td> <input type="text" name="10" value="0.97296" size="20"/>  &nbsp;&nbsp;(<code>default = <strong>0.97296</strong></code>; <code>minimum = 0.972</code>; <code>maximum = 0.974</code>)</td></tr></table>
The <i>V_cs</i> CKM matrix element.
  

<br/><br/><table><tr><td><strong>StandardModel:Vcb </td><td></td><td> <input type="text" name="11" value="0.04221" size="20"/>  &nbsp;&nbsp;(<code>default = <strong>0.04221</strong></code>; <code>minimum = 0.0418</code>; <code>maximum = 0.0426</code>)</td></tr></table>
The <i>V_cb</i> CKM matrix element.
  

<br/><br/><table><tr><td><strong>StandardModel:Vtd </td><td></td><td> <input type="text" name="12" value="0.00814" size="20"/>  &nbsp;&nbsp;(<code>default = <strong>0.00814</strong></code>; <code>minimum = 0.006</code>; <code>maximum = 0.010</code>)</td></tr></table>
The <i>V_td</i> CKM matrix element.
  

<br/><br/><table><tr><td><strong>StandardModel:Vts </td><td></td><td> <input type="text" name="13" value="0.04161" size="20"/>  &nbsp;&nbsp;(<code>default = <strong>0.04161</strong></code>; <code>minimum = 0.039</code>; <code>maximum = 0.043</code>)</td></tr></table>
The <i>V_ts</i> CKM matrix element.
  

<br/><br/><table><tr><td><strong>StandardModel:Vtb </td><td></td><td> <input type="text" name="14" value="0.9991" size="20"/>  &nbsp;&nbsp;(<code>default = <strong>0.9991</strong></code>; <code>minimum = 0.99907</code>; <code>maximum = 0.9992</code>)</td></tr></table>
The <i>V_tb</i> CKM matrix element.
  

<h3>The CoupSM class</h3> 

The <code><?php $filepath = $_GET["filepath"];
echo "<a href='ProgramFlow.php?filepath=".$filepath."' target='page'>";?>Pythia</a></code> class contains a
public instance <code>coupSM</code> of the <code>CoupSM</code> class.
This class contains one instance each of the <code>AlphaStrong</code>    
and <code>AlphaEM</code> classes, and additionally stores the weak couplings
and the quark mixing matrix mentioned above. This class is used especially
in the calculation of cross sections and resonance widths, but could also
be used elsewhere. Specifically, as already mentioned, there are separate 
<code>AlphaStrong</code> and <code>AlphaEM</code> instances for timelike 
and spacelike showers and for multiparton interactions, while weak couplings 
and the quark mixing matrix are only stored here. With the exception of the 
first two methods below, which are for internal use, the subsequent ones
could also be used externally.

<a name="method1"></a>
<p/><strong>CoupSM::CoupSM() &nbsp;</strong> <br/>
the constructor does nothing. Internal.
  

<a name="method2"></a>
<p/><strong>void CoupSM::init(Settings& settings, Rndm* rndmPtr) &nbsp;</strong> <br/>
this is where the <code>AlphaStrong</code> and <code>AlphaEM</code>
instances are initialized, and weak couplings and the quark mixing matrix
are read in and set. This is based on the values stored on this page and
among the <?php $filepath = $_GET["filepath"];
echo "<a href='CouplingsAndScales.php?filepath=".$filepath."' target='page'>";?>Couplings and Scales</a>. 
Internal.
  

<a name="method3"></a>
<p/><strong>double CoupSM::alphaS(double scale2) &nbsp;</strong> <br/>
the <i>alpha_strong</i> value at the quadratic scale <code>scale2</code>.
  

<a name="method4"></a>
<p/><strong>double CoupSM::alphaS1Ord(double scale2) &nbsp;</strong> <br/>
a first-order overestimate of the full second-order <i>alpha_strong</i> 
value at the quadratic scale <code>scale2</code>.
  

<a name="method5"></a>
<p/><strong>double CoupSM::alphaS2OrdCorr(double scale2) &nbsp;</strong> <br/>
a multiplicative correction factor, below unity, that brings the 
first-order overestimate above into agreement with the full second-order
<i>alpha_strong</i> value at the quadratic scale <code>scale2</code>.
  

<a name="method6"></a>
<p/><strong>double CoupSM::Lambda3() &nbsp;</strong> <br/>
  
<strong>double CoupSM::Lambda4() &nbsp;</strong> <br/>
  
<strong>double CoupSM::Lambda5() &nbsp;</strong> <br/>
the three-, four-, and five-flavour <i>Lambda</i> scale.
  

<a name="method7"></a>
<p/><strong>double CoupSM::alphaEM(double scale2) &nbsp;</strong> <br/>
the <i>alpha_em</i> value at the quadratic scale <code>scale2</code>.
  

<a name="method8"></a>
<p/><strong>double CoupSM::sin2thetaW() &nbsp;</strong> <br/>
  
<strong>double CoupSM::cos2thetaW() &nbsp;</strong> <br/>
the sine-squared and cosine-squared of the weak mixing angle, as used in 
the gauge-boson sector.
  

<a name="method9"></a>
<p/><strong>double CoupSM::sin2thetaWbar() &nbsp;</strong> <br/>
the sine-squared of the weak mixing angle, as used to derive the vector 
couplings of fermions to the <i>Z^0</i>.
  

<a name="method10"></a>
<p/><strong>double CoupSM::GF() &nbsp;</strong> <br/>
the Fermi constant of weak decays, in GeV<i>^-2</i>.
  

<a name="method11"></a>
<p/><strong>double CoupSM::ef(int idAbs) &nbsp;</strong> <br/>
the electrical charge of a fermion, by the absolute sign of the PDF code,
i.e. <code>idAbs</code> must be in the range between 1 and 18.
  

<a name="method12"></a>
<p/><strong>double CoupSM::vf(int idAbs) &nbsp;</strong> <br/>
  
<strong>double CoupSM::af(int idAbs) &nbsp;</strong> <br/>
the vector and axial charges of a fermion, by the absolute sign of the PDF 
code (<i>a_f = +-1, v_f = a_f - 4. * sin2thetaWbar * e_f</i>).
  

<a name="method13"></a>
<p/><strong>double CoupSM::t3f(int idAbs) &nbsp;</strong> <br/>
  
<strong>double CoupSM::lf(int idAbs) &nbsp;</strong> <br/>
  
<strong>double CoupSM::rf(int idAbs) &nbsp;</strong> <br/>
the weak isospin, left- and righthanded charges of a fermion, by the 
absolute sign of the PDF code (<i>t^3_f = a_f/2, l_f = (v_f + a_f)/2,
r_f = (v_f - a_f)/2</i>; you may find other conventions in the literature
that differ by a factor of 2).
  

<a name="method14"></a>
<p/><strong>double CoupSM::ef2(int idAbs) &nbsp;</strong> <br/>
  
<strong>double CoupSM::vf2(int idAbs) &nbsp;</strong> <br/>
  
<strong>double CoupSM::af2(int idAbs) &nbsp;</strong> <br/>
  
<strong>double CoupSM::efvf(int idAbs) &nbsp;</strong> <br/>
  
<strong>double CoupSM::vf2af2(int idAbs) &nbsp;</strong> <br/>
common quadratic combinations of the above couplings:
<i>e_f^2, v_f^2, a_f^2, e_f * v_f, v_f^2 + a_f^2</i>.
  

<a name="method15"></a>
<p/><strong>double CoupSM::VCKMgen(int genU, int genD) &nbsp;</strong> <br/>
  
<strong>double CoupSM::V2CKMgen(int genU, int genD) &nbsp;</strong> <br/>
the CKM mixing element,or the square of it, for
up-type generation index <code>genU</code> 
(<i>1 = u, 2 = c, 3 = t, 4 = t'</i>) and
down-type generation index <code>genD</code>
(<i>1 = d, 2 = s, 3 = b, 4 = b'</i>).
  

<a name="method16"></a>
<p/><strong>double CoupSM::VCKMid(int id1, int id2) &nbsp;</strong> <br/>
  
<strong>double CoupSM::V2CKMid(int id1, int id2) &nbsp;</strong> <br/>
the CKM mixing element,or the square of it, for
flavours <code>id1</code> and <code>id2</code>, both in the 
range from <i>-18</i> to <i>+18</i>. The sign is here not 
checked (so it can be used both for <i>u + dbar -> W+</i>
and <i>u -> d + W+</i>, say), but impossible flavour combinations
evaluate to zero. The neutrino sector is numbered by flavor
eigenstates, so there is no mixing in the lepton-neutrino system. 
  

<a name="method17"></a>
<p/><strong>double CoupSM::V2CKMsum(int id) &nbsp;</strong> <br/>
the sum of squared CKM mixing element that a given flavour can couple to, 
excluding the top quark and fourth generation. Is close to unity
for the first two generations. Returns unity for the lepton-neutrino
sector. 
  

<a name="method18"></a>
<p/><strong>int CoupSM::V2CKMpick(int id) &nbsp;</strong> <br/>
picks a random CKM partner quark or lepton (with the same sign as 
<code>id</code>) according to the respective squared elements, again 
excluding the top quark and fourth generation from the list of 
possibilities. Unambiguous choice for the lepton-neutrino sector. 
  

<input type="hidden" name="saved" value="1"/>

<?php
echo "<input type='hidden' name='filepath' value='".$_GET["filepath"]."'/>"?>

<table width="100%"><tr><td align="right"><input type="submit" value="Save Settings" /></td></tr></table>
</form>

<?php

if($_POST["saved"] == 1)
{
$filepath = $_POST["filepath"];
$handle = fopen($filepath, 'a');

if($_POST["1"] != "0.00729735")
{
$data = "StandardModel:alphaEM0 = ".$_POST["1"]."\n";
fwrite($handle,$data);
}
if($_POST["2"] != "0.00781751")
{
$data = "StandardModel:alphaEMmZ = ".$_POST["2"]."\n";
fwrite($handle,$data);
}
if($_POST["3"] != "0.2312")
{
$data = "StandardModel:sin2thetaW = ".$_POST["3"]."\n";
fwrite($handle,$data);
}
if($_POST["4"] != "0.2315")
{
$data = "StandardModel:sin2thetaWbar = ".$_POST["4"]."\n";
fwrite($handle,$data);
}
if($_POST["5"] != "1.16637e-5")
{
$data = "StandardModel:GF = ".$_POST["5"]."\n";
fwrite($handle,$data);
}
if($_POST["6"] != "0.97383")
{
$data = "StandardModel:Vud = ".$_POST["6"]."\n";
fwrite($handle,$data);
}
if($_POST["7"] != "0.2272")
{
$data = "StandardModel:Vus = ".$_POST["7"]."\n";
fwrite($handle,$data);
}
if($_POST["8"] != "0.00396")
{
$data = "StandardModel:Vub = ".$_POST["8"]."\n";
fwrite($handle,$data);
}
if($_POST["9"] != "0.2271")
{
$data = "StandardModel:Vcd = ".$_POST["9"]."\n";
fwrite($handle,$data);
}
if($_POST["10"] != "0.97296")
{
$data = "StandardModel:Vcs = ".$_POST["10"]."\n";
fwrite($handle,$data);
}
if($_POST["11"] != "0.04221")
{
$data = "StandardModel:Vcb = ".$_POST["11"]."\n";
fwrite($handle,$data);
}
if($_POST["12"] != "0.00814")
{
$data = "StandardModel:Vtd = ".$_POST["12"]."\n";
fwrite($handle,$data);
}
if($_POST["13"] != "0.04161")
{
$data = "StandardModel:Vts = ".$_POST["13"]."\n";
fwrite($handle,$data);
}
if($_POST["14"] != "0.9991")
{
$data = "StandardModel:Vtb = ".$_POST["14"]."\n";
fwrite($handle,$data);
}
fclose($handle);
}

?>
</body>
</html>

<!-- Copyright (C) 2012 Torbjorn Sjostrand -->