PYTHIA8/pythia8170/phpdoc/TimelikeShowers.php

   1 <html>
   2 <head>
   3 <title>Timelike Showers</title>
   4 <link rel="stylesheet" type="text/css" href="pythia.css"/>
   5 <link rel="shortcut icon" href="pythia32.gif"/>
   6 </head>
   7 <body>
   8
   9 <script language=javascript type=text/javascript>
  10 function stopRKey(evt) {
  11 var evt = (evt) ? evt : ((event) ? event : null);
  12 var node = (evt.target) ? evt.target :((evt.srcElement) ? evt.srcElement : null);
  13 if ((evt.keyCode == 13) && (node.type=="text"))
  14 {return false;}
  15 }
  16
  17 document.onkeypress = stopRKey;
  18 </script>
  19 <?php
  20 if($_POST['saved'] == 1) {
  21 if($_POST['filepath'] != "files/") {
  22 echo "<font color='red'>SETTINGS SAVED TO FILE</font><br/><br/>"; }
  23 else {
  24 echo "<font color='red'>NO FILE SELECTED YET.. PLEASE DO SO </font><a href='SaveSettings.php'>HERE</a><br/><br/>"; }
  25 }
  26 ?>
  27
  28 <form method='post' action='TimelikeShowers.php'>
  29
  30 <h2>Timelike Showers</h2>
  31
  32 The PYTHIA algorithm for timelike final-state showers is based on
  33 the article [<a href="Bibliography.php" target="page">Sjo05</a>], where a transverse-momentum-ordered
  34 evolution scheme is introduced, with the extension to fully interleaved
  35 evolution covered in [<a href="Bibliography.php" target="page">Cor10a</a>]. This algorithm is influenced by
  36 the previous mass-ordered algorithm in PYTHIA [<a href="Bibliography.php" target="page">Ben87</a>] and by
  37 the dipole-emission formulation in Ariadne [<a href="Bibliography.php" target="page">Gus86</a>]. From the
  38 mass-ordered algorithm it inherits a merging procedure for first-order
  39 gluon-emission matrix elements in essentially all two-body decays
  40 in the standard model and its minimal supersymmetric extension
  41 [<a href="Bibliography.php" target="page">Nor01</a>].
  42
  43 <p/>
  44 The normal user is not expected to call <code>TimeShower</code> directly,
  45 but only have it called from <code>Pythia</code>. Some of the parameters
  46 below, in particular <code>TimeShower:alphaSvalue</code>, would be of
  47 interest for a tuning exercise, however.
  48
  49 <h3>Main variables</h3>
  50
  51 Often the maximum scale of the FSR shower evolution is understood from the
  52 context. For instance, in a resonace decay half the resonance mass sets an
  53 absolute upper limit. For a hard process in a hadronic collision the choice
  54 is not as unique. Here the <?php $filepath = $_GET["filepath"];
  55 echo "<a href='CouplingsAndScales.php?filepath=".$filepath."' target='page'>";?>factorization
  56 scale</a> has been chosen as the maximum evolution scale. This would be
  57 the <i>pT</i> for a <i>2 -> 2</i> process, supplemented by mass terms
  58 for massive outgoing particles. For some special applications we do allow
  59 an alternative.
  60
  61 <br/><br/><table><tr><td><strong>TimeShower:pTmaxMatch  </td><td>  &nbsp;&nbsp;(<code>default = <strong>1</strong></code>; <code>minimum = 0</code>; <code>maximum = 2</code>)</td></tr></table>
  62 Way in which the maximum shower evolution scale is set to match the
  63 scale of the hard process itself.
  64 <br/>
  65 <input type="radio" name="1" value="0"><strong>0 </strong>: <b>(i)</b> if the final state of the hard process  (not counting subsequent resonance decays) contains at least one quark  (<ei>u, d, s, c ,b</ei>), gluon or photon then <ei>pT_max</ei>  is chosen to be the factorization scale for internal processes  and the <code>scale</code> value for Les Houches input;  <b>(ii)</b> if not, emissions are allowed to go all the way up to  the kinematical limit (i.e. to half the dipole mass).  This option agrees with the corresponding one for  <aloc href="SpacelikeShowers">spacelike showers</aloc>. There the  reasoning is that in the former set of processes the ISR emission of yet another quark, gluon or photon could lead to doublecounting, while no such danger exists in the latter case. The argument is less compelling for timelike showers, but could be a reasonable starting point. <br/>
  66 <input type="radio" name="1" value="1" checked="checked"><strong>1 </strong>: always use the factorization scale for an internal process and the <code>scale</code> value for Les Houches input,  i.e. the lower value. This should avoid doublecounting, but may leave out some emissions that ought to have been simulated. (Also known as wimpy showers.) <br/>
  67 <input type="radio" name="1" value="2"><strong>2 </strong>: always allow emissions up to the kinematical limit  (i.e. to half the dipole mass). This will simulate all possible event  topologies, but may lead to doublecounting.  (Also known as power showers.) <br/>
  68 <br/><b>Note:</b> These options only apply to the hard interaction.
  69 Emissions off subsequent multiparton interactions are always constrainted
  70 to be below the factorization scale of the process itself. They also
  71 assume you use interleaved evolution, so that FSR is in direct
  72 competition with ISR for the hardest emission. If you already
  73 generated a number of ISR partons at low <ei>pT</ei>, it would not
  74 make sense to have a later FSR shower up to the kinematical for all
  75 of them.
  76
  77 <br/><br/><table><tr><td><strong>TimeShower:pTmaxFudge </td><td></td><td> <input type="text" name="2" value="1.0" size="20"/>  &nbsp;&nbsp;(<code>default = <strong>1.0</strong></code>; <code>minimum = 0.25</code>; <code>maximum = 2.0</code>)</td></tr></table>
  78 In cases where the above <code>pTmaxMatch</code> rules would imply
  79 that <i>pT_max = pT_factorization</i>, <code>pTmaxFudge</code>
  80 introduces a multiplicative factor <i>f</i> such that instead
  81 <i>pT_max = f * pT_factorization</i>. Only applies to the hardest
  82 interaction in an event, cf. below. It is strongly suggested that
  83 <i>f = 1</i>, but variations around this default can be useful to
  84 test this assumption.
  85 <br/><b>Note:</b>Scales for resonance decays are not affected, but can
  86 be set separately by <?php $filepath = $_GET["filepath"];
  87 echo "<a href='UserHooks.php?filepath=".$filepath."' target='page'>";?>user hooks</a>.
  88
  89
  90 <br/><br/><table><tr><td><strong>TimeShower:pTmaxFudgeMPI </td><td></td><td> <input type="text" name="3" value="1.0" size="20"/>  &nbsp;&nbsp;(<code>default = <strong>1.0</strong></code>; <code>minimum = 0.25</code>; <code>maximum = 2.0</code>)</td></tr></table>
  91 A multiplicative factor <i>f</i> such that
  92 <i>pT_max = f * pT_factorization</i>, as above, but here for the
  93 non-hardest interactions (when multiparton interactions are allowed).
  94
  95
  96 <br/><br/><table><tr><td><strong>TimeShower:pTdampMatch  </td><td>  &nbsp;&nbsp;(<code>default = <strong>0</strong></code>; <code>minimum = 0</code>; <code>maximum = 2</code>)</td></tr></table>
  97 These options only take effect when a process is allowed to radiate up
  98 to the kinematical limit by the above <code>pTmaxMatch</code> choice,
  99 and no matrix-element corrections are available. Then, in many processes,
 100 the fall-off in <ei>pT</ei> will be too slow by one factor of <ei>pT^2</ei>.
 101 That is, while showers have an approximate <ei>dpT^2/pT^2</ei> shape, often
 102 it should become more like <ei>dpT^2/pT^4</ei> at <ei>pT</ei> values above
 103 the scale of the hard process. This argument is more obvious for ISR,
 104 but is taken over unchanged for FSR to have a symmetric description.
 105 <br/>
 106 <input type="radio" name="4" value="0" checked="checked"><strong>0 </strong>: emissions go up to the kinematical limit,  with no special dampening. <br/>
 107 <input type="radio" name="4" value="1"><strong>1 </strong>: emissions go up to the kinematical limit,   but dampened by a factor <ei>k^2 Q^2_fac/(pT^2 + k^2 Q^2_fac)</ei>,  where <ei>Q_fac</ei> is the factorization scale and <ei>k</ei> is a  multiplicative fudge factor stored in <code>pTdampFudge</code> below. <br/>
 108 <input type="radio" name="4" value="2"><strong>2 </strong>: emissions go up to the kinematical limit,  but dampened by a factor <ei>k^2 Q^2_ren/(pT^2 + k^2 Q^2_ren)</ei>,  where <ei>Q_ren</ei> is the renormalization scale and <ei>k</ei> is a  multiplicative fudge factor stored in <code>pTdampFudge</code> below.  <br/>
 109 <br/><b>Note:</b> These options only apply to the hard interaction.
 110 Emissions off subsequent multiparton interactions are always constrainted
 111 to be below the factorization scale of the process itself.
 112
 113 <br/><br/><table><tr><td><strong>TimeShower:pTdampFudge </td><td></td><td> <input type="text" name="5" value="1.0" size="20"/>  &nbsp;&nbsp;(<code>default = <strong>1.0</strong></code>; <code>minimum = 0.25</code>; <code>maximum = 4.0</code>)</td></tr></table>
 114 In cases 1 and 2 above, where a dampening is imposed at around the
 115 factorization or renormalization scale, respectively, this allows the
 116 <i>pT</i> scale of dampening of radiation by a half to be shifted
 117 by this factor relative to the default <i>Q_fac</i> or <i>Q_ren</i>.
 118 This number ought to be in the neighbourhood of unity, but variations
 119 away from this value could do better in some processes.
 120
 121
 122 <p/>
 123 The amount of QCD radiation in the shower is determined by
 124 <br/><br/><table><tr><td><strong>TimeShower:alphaSvalue </td><td></td><td> <input type="text" name="6" value="0.1383" size="20"/>  &nbsp;&nbsp;(<code>default = <strong>0.1383</strong></code>; <code>minimum = 0.06</code>; <code>maximum = 0.25</code>)</td></tr></table>
 125 The <i>alpha_strong</i> value at scale <i>M_Z^2</i>. The default
 126 value corresponds to a crude tuning to LEP data, to be improved.
 127
 128
 129 <p/>
 130 The actual value is then regulated by the running to the scale
 131 <i>pT^2</i>, at which the shower evaluates <i>alpha_strong</i>.
 132
 133 <br/><br/><table><tr><td><strong>TimeShower:alphaSorder  </td><td>  &nbsp;&nbsp;(<code>default = <strong>1</strong></code>; <code>minimum = 0</code>; <code>maximum = 2</code>)</td></tr></table>
 134 Order at which <ei>alpha_strong</ei> runs,
 135 <br/>
 136 <input type="radio" name="7" value="0"><strong>0 </strong>: zeroth order, i.e. <ei>alpha_strong</ei> is kept  fixed.<br/>
 137 <input type="radio" name="7" value="1" checked="checked"><strong>1 </strong>: first order, which is the normal value.<br/>
 138 <input type="radio" name="7" value="2"><strong>2 </strong>: second order. Since other parts of the code do  not go to second order there is no strong reason to use this option,  but there is also nothing wrong with it.<br/>
 139
 140 <p/>
 141 QED radiation is regulated by the <i>alpha_electromagnetic</i>
 142 value at the <i>pT^2</i> scale of a branching.
 143
 144 <br/><br/><table><tr><td><strong>TimeShower:alphaEMorder  </td><td>  &nbsp;&nbsp;(<code>default = <strong>1</strong></code>; <code>minimum = -1</code>; <code>maximum = 1</code>)</td></tr></table>
 145 The running of <ei>alpha_em</ei>.
 146 <br/>
 147 <input type="radio" name="8" value="1" checked="checked"><strong>1 </strong>: first-order running, constrained to agree with <code>StandardModel:alphaEMmZ</code> at the <ei>Z^0</ei> mass. <br/>
 148 <input type="radio" name="8" value="0"><strong>0 </strong>: zeroth order, i.e. <ei>alpha_em</ei> is kept  fixed at its value at vanishing momentum transfer.<br/>
 149 <input type="radio" name="8" value="-1"><strong>-1 </strong>: zeroth order, i.e. <ei>alpha_em</ei> is kept  fixed, but at <code>StandardModel:alphaEMmZ</code>, i.e. its value at the <ei>Z^0</ei> mass. <br/>
 150
 151 <p/>
 152 The natural scale for couplings, and PDFs for dipoles stretching out
 153 to the beam remnants, is <i>pT^2</i>. To explore uncertainties it
 154 is possibly to vary around this value, however, in analogy with what
 155 can be done for <?php $filepath = $_GET["filepath"];
 156 echo "<a href='CouplingsAndScales.php?filepath=".$filepath."' target='page'>";?>hard processes</a>.
 157
 158 <br/><br/><table><tr><td><strong>TimeShower:renormMultFac </td><td></td><td> <input type="text" name="9" value="1." size="20"/>  &nbsp;&nbsp;(<code>default = <strong>1.</strong></code>; <code>minimum = 0.1</code>; <code>maximum = 10.</code>)</td></tr></table>
 159 The default <i>pT^2</i> renormalization scale is multiplied by
 160 this prefactor. For QCD this is equivalent to a change of
 161 <i>Lambda^2</i> in the opposite direction, i.e. to a change of
 162 <i>alpha_strong(M_Z^2)</i> (except that flavour thresholds
 163 remain at fixed scales).
 164
 165
 166 <br/><br/><table><tr><td><strong>TimeShower:factorMultFac </td><td></td><td> <input type="text" name="10" value="1." size="20"/>  &nbsp;&nbsp;(<code>default = <strong>1.</strong></code>; <code>minimum = 0.1</code>; <code>maximum = 10.</code>)</td></tr></table>
 167 The default <i>pT^2</i> factorization scale is multiplied by
 168 this prefactor.
 169
 170
 171 <p/>
 172 The rate of radiation if divergent in the <i>pT -> 0</i> limit. Here,
 173 however, perturbation theory is expected to break down. Therefore an
 174 effective <i>pT_min</i> cutoff parameter is introduced, below which
 175 no emissions are allowed. The cutoff may be different for QCD and QED
 176 radiation off quarks, and is mainly a technical parameter for QED
 177 radiation off leptons.
 178
 179 <br/><br/><table><tr><td><strong>TimeShower:pTmin </td><td></td><td> <input type="text" name="11" value="0.4" size="20"/>  &nbsp;&nbsp;(<code>default = <strong>0.4</strong></code>; <code>minimum = 0.1</code>; <code>maximum = 2.0</code>)</td></tr></table>
 180 Parton shower cut-off <i>pT</i> for QCD emissions.
 181
 182
 183 <br/><br/><table><tr><td><strong>TimeShower:pTminChgQ </td><td></td><td> <input type="text" name="12" value="0.4" size="20"/>  &nbsp;&nbsp;(<code>default = <strong>0.4</strong></code>; <code>minimum = 0.1</code>; <code>maximum = 2.0</code>)</td></tr></table>
 184 Parton shower cut-off <i>pT</i> for photon coupling to coloured particle.
 185
 186
 187 <br/><br/><table><tr><td><strong>TimeShower:pTminChgL </td><td></td><td> <input type="text" name="13" value="0.0005" size="20"/>  &nbsp;&nbsp;(<code>default = <strong>0.0005</strong></code>; <code>minimum = 0.0001</code>; <code>maximum = 2.0</code>)</td></tr></table>
 188 Parton shower cut-off <i>pT</i> for pure QED branchings.
 189 Assumed smaller than (or equal to) <code>pTminChgQ</code>.
 190
 191
 192 <p/>
 193 Shower branchings <i>gamma -> f fbar</i>, where <i>f</i> is a
 194 quark or lepton, in part compete with the hard processes involving
 195 <i>gamma^*/Z^0</i> production. In order to avoid overlap it makes
 196 sense to correlate the maximum <i>gamma</i> mass allowed in showers
 197 with the minumum <i>gamma^*/Z^0</i> mass allowed in hard processes.
 198 In addition, the shower contribution only contains the pure
 199 <i>gamma^*</i> contribution, i.e. not the <i>Z^0</i> part, so
 200 the mass spectrum above 50 GeV or so would not be well described.
 201
 202 <br/><br/><table><tr><td><strong>TimeShower:mMaxGamma </td><td></td><td> <input type="text" name="14" value="10.0" size="20"/>  &nbsp;&nbsp;(<code>default = <strong>10.0</strong></code>; <code>minimum = 0.001</code>; <code>maximum = 50.0</code>)</td></tr></table>
 203 Maximum invariant mass allowed for the created fermion pair in a
 204 <i>gamma -> f fbar</i> branching in the shower.
 205
 206
 207 <h3>Interleaved evolution</h3>
 208
 209 Multiparton interactions (MPI) and initial-state showers (ISR) are
 210 always interleaved, as follows. Starting from the hard interaction,
 211 the complete event is constructed by a set of steps. In each step
 212 the <i>pT</i> scale of the previous step is used as starting scale
 213 for a downwards evolution. The MPI and ISR components each make
 214 their respective Monte Carlo choices for the next lower <i>pT</i>
 215 value. The one with larger <i>pT</i> is allowed to carry out its
 216 proposed action, thereby modifying the conditions for the next steps.
 217 This is relevant since the two components compete for the energy
 218 contained in the beam remnants: both an interaction and an emission
 219 take avay some of the energy, leaving less for the future. The end
 220 result is a combined chain of decreasing <i>pT</i> values, where
 221 ones associated with new interactions and ones with new emissions
 222 are interleaved.
 223
 224 <p/>
 225 There is no corresponding requirement for final-state radiation (FSR)
 226 to be interleaved. Such an FSR emission does not compete directly for
 227 beam energy (but see below), and also can be viewed as occuring after
 228 the other two components in some kind of time sense. Interleaving is
 229 allowed, however, since it can be argued that a high-<i>pT</i> FSR
 230 occurs on shorter time scales than a low-<i>pT</i> MPI, say.
 231 Backwards evolution of ISR is also an example that physical time
 232 is not the only possible ordering principle, but that one can work
 233 with conditional probabilities: given the partonic picture at a
 234 specific <i>pT</i> resolution scale, what possibilities are open
 235 for a modified picture at a slightly lower <i>pT</i> scale, either
 236 by MPI, ISR or FSR? Complete interleaving of the three components also
 237 offers advantages if one aims at matching to higher-order matrix
 238 elements above some given scale.
 239
 240 <br/><br/><strong>TimeShower:interleave</strong>  <input type="radio" name="15" value="on" checked="checked"><strong>On</strong>
 241 <input type="radio" name="15" value="off"><strong>Off</strong>
 242  &nbsp;&nbsp;(<code>default = <strong>on</strong></code>)<br/>
 243 If on, final-state emissions are interleaved in the same
 244 decreasing-<i>pT</i> chain as multiparton interactions and initial-state
 245 emissions. If off, final-state emissions are only addressed after the
 246 multiparton interactions and initial-state radiation have been considered.
 247
 248
 249 <p/>
 250 As an aside, it should be noted that such interleaving does not affect
 251 showering in resonance decays, such as a <i>Z^0</i>. These decays are
 252 only introduced after the production process has been considered in full,
 253 and the subsequent FSR is carried out inside the resonance, with
 254 preserved resonance mass.
 255
 256 <p/>
 257 One aspect of FSR for a hard process in hadron collisions is that often
 258 colour diples are formed between a scattered parton and a beam remnant,
 259 or rather the hole left behind by an incoming partons. If such holes
 260 are allowed as dipole ends and take the recoil when the scattered parton
 261 undergoes a branching then this translates into the need to take some
 262 amount of remnant energy also in the case of FSR, i.e. the roles of
 263 ISR and FSR are not completely decoupled. The energy taken away is
 264 bokkept by increasing the <i>x</i> value assigned to the incoming
 265 scattering parton, and a reweighting factor
 266 <i>x_new f(x_new, pT^2) / x_old f(x_old, pT^2)</i>
 267 in the emission probability ensures that not unphysically large
 268 <i>x_new</i> values are reached. Usually such <i>x</i> changes are
 269 small, and they can be viewed as a higher-order effect beyond the
 270 accuracy of the leading-log initial-state showers.
 271
 272 <p/>
 273 This choice is not unique, however. As an alternative, if nothing else
 274 useful for cross-checks, one could imagine that the FSR is completely
 275 decoupled from the ISR and beam remnants.
 276
 277 <br/><br/><strong>TimeShower:allowBeamRecoil</strong>  <input type="radio" name="16" value="on" checked="checked"><strong>On</strong>
 278 <input type="radio" name="16" value="off"><strong>Off</strong>
 279  &nbsp;&nbsp;(<code>default = <strong>on</strong></code>)<br/>
 280 If on, the final-state shower is allowed to borrow energy from
 281 the beam remnants as described above, thereby changing the mass of the
 282 scattering subsystem. If off, the partons in the scattering subsystem
 283 are constrained to borrow energy from each other, such that the total
 284 four-momentum of the system is preserved. This flag has no effect
 285 on resonance decays, where the shower always preserves the resonance
 286 mass, cf. the comment above about showers for resonances never being
 287 interleaved.
 288
 289
 290 <br/><br/><strong>TimeShower:dampenBeamRecoil</strong>  <input type="radio" name="17" value="on" checked="checked"><strong>On</strong>
 291 <input type="radio" name="17" value="off"><strong>Off</strong>
 292  &nbsp;&nbsp;(<code>default = <strong>on</strong></code>)<br/>
 293 When beam recoil is allowed there is still some ambiguity how far
 294 into the beam end of the dipole that emission should be allowed.
 295 It is dampened in the beam region, but probably not enough.
 296 When on an additional suppression factor
 297 <i>4 pT2_hard / (4 pT2_hard + m2)</i> is multiplied on to the
 298 emission probability. Here <i>pT_hard</i> is the transverse momentum
 299 of the radiating parton and <i>m</i> the off-shell mass it acquires
 300 by the branching, <i>m2 = pT2/(z(1-z))</i>. Note that
 301 <i>m2 = 4 pT2_hard</i> is the kinematical limit for a scattering
 302 at 90 degrees without beam recoil.
 303
 304
 305 <h3>Global recoil</h3>
 306
 307 The final-state algorithm is based on dipole-style recoils, where
 308 one single parton takes the full recoil of a branching. This is unlike
 309 the initial-state algorithm, where the complete already-existing
 310 final state shares the recoil of each new emission. As an alternative,
 311 also the final-state algorithm contains an option where the recoil
 312 is shared between all partons in the final state. Thus the radiation
 313 pattern is unrelated to colour correlations. This is especially
 314 convenient for some matching algorithms, like MC@NLO, where a full
 315 analytic knowledge of the shower radiation pattern is needed to avoid
 316 doublecountning. (The <i>pT</i>-ordered shower is described in
 317 [<a href="Bibliography.php" target="page">Sjo05</a>], and the corrections for massive radiator and recoiler
 318 in [<a href="Bibliography.php" target="page">Nor01</a>].)
 319
 320 <p/>
 321 Technically, the radiation pattern is most conveniently represented
 322 in the rest frame of the final state of the hard subprocess. Then, for
 323 each parton at a time, the rest of the final state can be viewed as
 324 a single effective parton. This "parton" has a fixed invariant mass
 325 during the emission process, and takes the recoil without any changed
 326 direction of motion. The momenta of the individual new recoilers are
 327 then obtained by a simple common boost of the original ones.
 328
 329 <p/>
 330 This alternative approach will miss out on the colour coherence
 331 phenomena. Specifically, with the whole subcollision mass as "dipole"
 332 mass, the phase space for subsequent emissions is larger than for
 333 the normal dipole algorithm. The phase space difference grows as
 334 more and more gluons are created, and thus leads to a way too steep
 335 multiplication of soft gluons. Therefore the main application is
 336 for the first one or few emissions of the shower, where a potential
 337 overestimate of the emission rate is to be corrected for anyway,
 338 by matching to the relevant matrix elements. Thereafter, subsequent
 339 emissions should be handled as before, i.e. with dipoles spanned
 340 between nearby partons. Furthermore, only the first (hardest)
 341 subcollision is handled with global recoils, since subsequent MPI's
 342 would not be subject to matrix element corrections anyway.
 343
 344 <p/>
 345 In order for the mid-shower switch from global to local recoils
 346 to work, colours are traced and bookkept just as for normal showers;
 347 it is only that this information is not used in those steps where
 348 a global recoil is requested. (Thus, e.g., a gluon is still bookkept
 349 as one colour and one anticolour dipole end, with half the charge
 350 each, but with global recoil those two ends radiate identically.)
 351
 352 <br/><br/><strong>TimeShower:globalRecoil</strong>  <input type="radio" name="18" value="on"><strong>On</strong>
 353 <input type="radio" name="18" value="off" checked="checked"><strong>Off</strong>
 354  &nbsp;&nbsp;(<code>default = <strong>off</strong></code>)<br/>
 355 Alternative approach as above, where all final-state particles share
 356 the recoil of an emission.
 357 <br/>If off, then use the standard dipole-recoil approach.
 358 <br/>If on, use the alternative global recoil, but only for the first
 359 interaction, and only while the number of particles in the final state
 360 is at most <code>TimeShower:nMaxGlobalRecoil</code> before the
 361 branching.
 362
 363
 364 <br/><br/><table><tr><td><strong>TimeShower:nMaxGlobalRecoil  </td><td></td><td> <input type="text" name="19" value="2" size="20"/>  &nbsp;&nbsp;(<code>default = <strong>2</strong></code>; <code>minimum = 1</code>)</td></tr></table>
 365 Represents the maximum number of particles in the final state for which
 366 the next final-state emission can be performed with the global recoil
 367 strategy. This number counts all particles, whether they are
 368 allowed to radiate or not, e.g. also <i>Z^0</i>. Also partons
 369 created by initial-state radiation emissions counts towards this sum,
 370 as part of the interleaved evolution. Without interleaved evolution
 371 this option would not make sense, since then a varying and large
 372 number of partons could already have been created by the initial-state
 373 radiation before the first final-state one, and then there is not
 374 likely to be any matrix elements available for matching.
 375
 376
 377 <p/>
 378 The global-recoil machinery does not work well with rescattering in the
 379 MPI machinery, since then the recoiling system is not uniquely defined.
 380 <code>MultipartonInteractions:allowRescatter = off</code> by default,
 381 so this is not a main issue. If both options are switched on,
 382 rescattering will only be allowed to kick in after the global recoil
 383 has ceased to be active, i.e. once the <code>nMaxGlobalRecoil</code>
 384 limit has been exceeded. This should not be a major conflict,
 385 since rescattering is mainly of interest at later stages of the
 386 downwards <i>pT</i> evolution.
 387
 388 <p/>
 389 Further, it is strongly recommended to set
 390 <code>TimeShower:MEcorrections = off</code> (not default!), i.e. not
 391 to correct the emission probability to the internal matrix elements.
 392 The internal ME options do not cover any cases relevant for a multibody
 393 recoiler anyway, so no guarantees are given what prescription would
 394 come to be used. Instead, without ME corrections,  a process-independent
 395 emission rate is obtained, and  <?php $filepath = $_GET["filepath"];
 396 echo "<a href='UserHooks.php?filepath=".$filepath."' target='page'>";?>user hooks</a>
 397 can provide the desired process-specific rejection factors.
 398
 399 <h3>Radiation off octet onium states</h3>
 400
 401 In the current implementation, charmonium and bottomonium production
 402 can proceed either through colour singlet or colour octet mechanisms,
 403 both of them implemented in terms of <i>2 -> 2</i> hard processes
 404 such as <i>g g -> (onium) g</i>.
 405 In the former case the state does not radiate and the onium therefore
 406 is produced in isolation, up to normal underlying-event activity. In
 407 the latter case the situation is not so clear, but it is sensible to
 408 assume that a shower can evolve. (Assuming, of course, that the
 409 transverse momentum of the onium state is sufficiently high that
 410 radiation is of relevance.)
 411
 412 <p/>
 413 There could be two parts to such a shower. Firstly a gluon (or even a
 414 quark, though less likely) produced in a hard <i>2 -> 2</i> process
 415 can undergo showering into many gluons, whereof one branches into the
 416 heavy-quark pair. Secondly, once the pair has been produced, each quark
 417 can radiate further gluons. This latter kind of emission could easily
 418 break up a semibound quark pair, but might also create a new semibound
 419 state where before an unbound pair existed, and to some approximation
 420 these two effects should balance in the onium production rate.
 421 The showering "off an onium state" as implemented here therefore should
 422 not be viewed as an accurate description of the emission history
 423 step by step, but rather as an effective approach to ensure that the
 424 octet onium produced "in the hard process" is embedded in a realistic
 425 amount of jet activity.
 426 Of course both the isolated singlet and embedded octet are likely to
 427 be extremes, but hopefully the mix of the two will strike a reasonable
 428 balance. However, it is possible that some part of the octet production
 429 occurs in channels where it should not be accompanied by (hard) radiation.
 430 Therefore reducing the fraction of octet onium states allowed to radiate
 431 is a valid variation to explore uncertainties.
 432
 433 <p/>
 434 If an octet onium state is chosen to radiate, the simulation of branchings
 435 is based on the assumption that the full radiation is provided by an
 436 incoherent sum of radiation off the quark and off the antiquark of the
 437 onium state. Thus the splitting kernel is taken to be the normal
 438 <i>q -> q g</i> one, multiplied by a factor of two. Obviously this is
 439 a simplification of a more complex picture, averaging over factors pulling
 440 in different directions. Firstly, radiation off a gluon ought
 441 to be enhanced by a factor 9/4 relative to a quark rather than the 2
 442 now used, but this is a minor difference. Secondly, our use of the
 443 <i>q -> q g</i> branching kernel is roughly equivalent to always
 444 following the harder gluon in a <i>g -> g g</i> branching. This could
 445 give us a bias towards producing too hard onia. A soft gluon would have
 446 little phase space to branch into a heavy-quark pair however, so the
 447 bias may not be as big as it would seem at first glance. Thirdly,
 448 once the gluon has branched into a quark pair, each quark carries roughly
 449 only half of the onium energy. The maximum energy per emitted gluon should
 450 then be roughly half the onium energy rather than the full, as it is now.
 451 Thereby the energy of radiated gluons is exaggerated, i.e. onia become too
 452 soft. So the second and the third points tend to cancel each other.
 453
 454 <p/>
 455 Finally, note that the lower cutoff scale of the shower evolution depends
 456 on the onium mass rather than on the quark mass, as it should be. Gluons
 457 below the octet-onium scale should only be part of the octet-to-singlet
 458 transition.
 459
 460 <br/><br/><table><tr><td><strong>TimeShower:octetOniumFraction </td><td></td><td> <input type="text" name="20" value="1." size="20"/>  &nbsp;&nbsp;(<code>default = <strong>1.</strong></code>; <code>minimum = 0.</code>; <code>maximum = 1.</code>)</td></tr></table>
 461 Allow colour-octet charmonium and bottomonium states to radiate gluons.
 462 0 means that no octet-onium states radiate, 1 that all do, with possibility
 463 to interpolate between these two extremes.
 464
 465
 466 <br/><br/><table><tr><td><strong>TimeShower:octetOniumColFac </td><td></td><td> <input type="text" name="21" value="2." size="20"/>  &nbsp;&nbsp;(<code>default = <strong>2.</strong></code>; <code>minimum = 0.</code>; <code>maximum = 4.</code>)</td></tr></table>
 467 The colour factor used used in the splitting kernel for those octet onium
 468 states that are allowed to radiate, normalized to the <i>q -> q g</i>
 469 splitting kernel. Thus the default corresponds to twice the radiation
 470 off a quark. The physically preferred range would be between 1 and 9/4.
 471
 472
 473 <h3>Further variables</h3>
 474
 475 There are several possibilities you can use to switch on or off selected
 476 branching types in the shower, or in other respects simplify the shower.
 477 These should normally not be touched. Their main function is for
 478 cross-checks.
 479
 480 <br/><br/><strong>TimeShower:QCDshower</strong>  <input type="radio" name="22" value="on" checked="checked"><strong>On</strong>
 481 <input type="radio" name="22" value="off"><strong>Off</strong>
 482  &nbsp;&nbsp;(<code>default = <strong>on</strong></code>)<br/>
 483 Allow a QCD shower, i.e. branchings <i>q -> q g</i>, <i>g -> g g</i>
 484 and <i>g -> q qbar</i>; on/off = true/false.
 485
 486
 487 <br/><br/><table><tr><td><strong>TimeShower:nGluonToQuark  </td><td></td><td> <input type="text" name="23" value="5" size="20"/>  &nbsp;&nbsp;(<code>default = <strong>5</strong></code>; <code>minimum = 0</code>; <code>maximum = 5</code>)</td></tr></table>
 488 Number of allowed quark flavours in <i>g -> q qbar</i> branchings
 489 (phase space permitting). A change to 4 would exclude
 490 <i>g -> b bbar</i>, etc.
 491
 492
 493 <br/><br/><strong>TimeShower:QEDshowerByQ</strong>  <input type="radio" name="24" value="on" checked="checked"><strong>On</strong>
 494 <input type="radio" name="24" value="off"><strong>Off</strong>
 495  &nbsp;&nbsp;(<code>default = <strong>on</strong></code>)<br/>
 496 Allow quarks to radiate photons, i.e. branchings <i>q -> q gamma</i>;
 497 on/off = true/false.
 498
 499
 500 <br/><br/><strong>TimeShower:QEDshowerByL</strong>  <input type="radio" name="25" value="on" checked="checked"><strong>On</strong>
 501 <input type="radio" name="25" value="off"><strong>Off</strong>
 502  &nbsp;&nbsp;(<code>default = <strong>on</strong></code>)<br/>
 503 Allow leptons to radiate photons, i.e. branchings <i>l -> l gamma</i>;
 504 on/off = true/false.
 505
 506
 507 <br/><br/><strong>TimeShower:QEDshowerByGamma</strong>  <input type="radio" name="26" value="on" checked="checked"><strong>On</strong>
 508 <input type="radio" name="26" value="off"><strong>Off</strong>
 509  &nbsp;&nbsp;(<code>default = <strong>on</strong></code>)<br/>
 510 Allow photons to branch into lepton or quark pairs, i.e. branchings
 511 <i>gamma -> l+ l-</i> and <i>gamma -> q qbar</i>;
 512 on/off = true/false.
 513
 514
 515 <br/><br/><table><tr><td><strong>TimeShower:nGammaToQuark  </td><td></td><td> <input type="text" name="27" value="5" size="20"/>  &nbsp;&nbsp;(<code>default = <strong>5</strong></code>; <code>minimum = 0</code>; <code>maximum = 5</code>)</td></tr></table>
 516 Number of allowed quark flavours in <i>gamma -> q qbar</i> branchings
 517 (phase space permitting). A change to 4 would exclude
 518 <i>g -> b bbar</i>, etc.
 519
 520
 521 <br/><br/><table><tr><td><strong>TimeShower:nGammaToLepton  </td><td></td><td> <input type="text" name="28" value="3" size="20"/>  &nbsp;&nbsp;(<code>default = <strong>3</strong></code>; <code>minimum = 0</code>; <code>maximum = 3</code>)</td></tr></table>
 522 Number of allowed lepton flavours in <i>gamma -> l+ l-</i> branchings
 523 (phase space permitting). A change to 2 would exclude
 524 <i>gamma -> tau+ tau-</i>, and a change to 1 also
 525 <i>gamma -> mu+ mu-</i>.
 526
 527
 528 <br/><br/><strong>TimeShower:MEcorrections</strong>  <input type="radio" name="29" value="on" checked="checked"><strong>On</strong>
 529 <input type="radio" name="29" value="off"><strong>Off</strong>
 530  &nbsp;&nbsp;(<code>default = <strong>on</strong></code>)<br/>
 531 Use of matrix element corrections where available; on/off = true/false.
 532
 533
 534 <br/><br/><strong>TimeShower:MEafterFirst</strong>  <input type="radio" name="30" value="on" checked="checked"><strong>On</strong>
 535 <input type="radio" name="30" value="off"><strong>Off</strong>
 536  &nbsp;&nbsp;(<code>default = <strong>on</strong></code>)<br/>
 537 Use of matrix element corrections also after the first emission,
 538 for dipole ends of the same system that did not yet radiate.
 539 Only has a meaning if <code>MEcorrections</code> above is
 540 switched on.
 541
 542
 543 <br/><br/><strong>TimeShower:phiPolAsym</strong>  <input type="radio" name="31" value="on" checked="checked"><strong>On</strong>
 544 <input type="radio" name="31" value="off"><strong>Off</strong>
 545  &nbsp;&nbsp;(<code>default = <strong>on</strong></code>)<br/>
 546 Azimuthal asymmetry induced by gluon polarization; on/off = true/false.
 547
 548
 549 <br/><br/><strong>TimeShower:recoilToColoured</strong>  <input type="radio" name="32" value="on" checked="checked"><strong>On</strong>
 550 <input type="radio" name="32" value="off"><strong>Off</strong>
 551  &nbsp;&nbsp;(<code>default = <strong>on</strong></code>)<br/>
 552 In the decays of coloured resonances, say <i>t -> b W</i>, it is not
 553 possible to set up dipoles with matched colours. Originally the
 554 <i>b</i> radiator therefore has <i>W</i> as recoiler, and that
 555 choice is unique. Once a gluon has been radiated, however, it is
 556 possible either to have the unmatched colour (inherited by the gluon)
 557 still recoiling against the <i>W</i> (<code>off</code>), or else
 558 let it recoil against the <i>b</i> also for this dipole
 559 (<code>on</code>). Before version 8.160 the former was the only
 560 possibility, which could give unphysical radiation patterns. It is
 561 kept as an option to check backwards compatibility. The same issue
 562 exists for QED radiation, but obviously is less significant. Consider
 563 the example <i>W -> e nu</i>, where originally the <i>nu</i>
 564 takes the recoil. In the old (<code>off</code>) scheme the <i>nu</i>
 565 would remain recoiler, while in the new (<code>on</code>) instead
 566 each newly emitted photon becomes the new recoiler.
 567
 568
 569 <input type="hidden" name="saved" value="1"/>
 570
 571 <?php
 572 echo "<input type='hidden' name='filepath' value='".$_GET["filepath"]."'/>"?>
 573
 574 <table width="100%"><tr><td align="right"><input type="submit" value="Save Settings" /></td></tr></table>
 575 </form>
 576
 577 <?php
 578
 579 if($_POST["saved"] == 1)
 580 {
 581 $filepath = $_POST["filepath"];
 582 $handle = fopen($filepath, 'a');
 583
 584 if($_POST["1"] != "1")
 585 {
 586 $data = "TimeShower:pTmaxMatch = ".$_POST["1"]."\n";
 587 fwrite($handle,$data);
 588 }
 589 if($_POST["2"] != "1.0")
 590 {
 591 $data = "TimeShower:pTmaxFudge = ".$_POST["2"]."\n";
 592 fwrite($handle,$data);
 593 }
 594 if($_POST["3"] != "1.0")
 595 {
 596 $data = "TimeShower:pTmaxFudgeMPI = ".$_POST["3"]."\n";
 597 fwrite($handle,$data);
 598 }
 599 if($_POST["4"] != "0")
 600 {
 601 $data = "TimeShower:pTdampMatch = ".$_POST["4"]."\n";
 602 fwrite($handle,$data);
 603 }
 604 if($_POST["5"] != "1.0")
 605 {
 606 $data = "TimeShower:pTdampFudge = ".$_POST["5"]."\n";
 607 fwrite($handle,$data);
 608 }
 609 if($_POST["6"] != "0.1383")
 610 {
 611 $data = "TimeShower:alphaSvalue = ".$_POST["6"]."\n";
 612 fwrite($handle,$data);
 613 }
 614 if($_POST["7"] != "1")
 615 {
 616 $data = "TimeShower:alphaSorder = ".$_POST["7"]."\n";
 617 fwrite($handle,$data);
 618 }
 619 if($_POST["8"] != "1")
 620 {
 621 $data = "TimeShower:alphaEMorder = ".$_POST["8"]."\n";
 622 fwrite($handle,$data);
 623 }
 624 if($_POST["9"] != "1.")
 625 {
 626 $data = "TimeShower:renormMultFac = ".$_POST["9"]."\n";
 627 fwrite($handle,$data);
 628 }
 629 if($_POST["10"] != "1.")
 630 {
 631 $data = "TimeShower:factorMultFac = ".$_POST["10"]."\n";
 632 fwrite($handle,$data);
 633 }
 634 if($_POST["11"] != "0.4")
 635 {
 636 $data = "TimeShower:pTmin = ".$_POST["11"]."\n";
 637 fwrite($handle,$data);
 638 }
 639 if($_POST["12"] != "0.4")
 640 {
 641 $data = "TimeShower:pTminChgQ = ".$_POST["12"]."\n";
 642 fwrite($handle,$data);
 643 }
 644 if($_POST["13"] != "0.0005")
 645 {
 646 $data = "TimeShower:pTminChgL = ".$_POST["13"]."\n";
 647 fwrite($handle,$data);
 648 }
 649 if($_POST["14"] != "10.0")
 650 {
 651 $data = "TimeShower:mMaxGamma = ".$_POST["14"]."\n";
 652 fwrite($handle,$data);
 653 }
 654 if($_POST["15"] != "on")
 655 {
 656 $data = "TimeShower:interleave = ".$_POST["15"]."\n";
 657 fwrite($handle,$data);
 658 }
 659 if($_POST["16"] != "on")
 660 {
 661 $data = "TimeShower:allowBeamRecoil = ".$_POST["16"]."\n";
 662 fwrite($handle,$data);
 663 }
 664 if($_POST["17"] != "on")
 665 {
 666 $data = "TimeShower:dampenBeamRecoil = ".$_POST["17"]."\n";
 667 fwrite($handle,$data);
 668 }
 669 if($_POST["18"] != "off")
 670 {
 671 $data = "TimeShower:globalRecoil = ".$_POST["18"]."\n";
 672 fwrite($handle,$data);
 673 }
 674 if($_POST["19"] != "2")
 675 {
 676 $data = "TimeShower:nMaxGlobalRecoil = ".$_POST["19"]."\n";
 677 fwrite($handle,$data);
 678 }
 679 if($_POST["20"] != "1.")
 680 {
 681 $data = "TimeShower:octetOniumFraction = ".$_POST["20"]."\n";
 682 fwrite($handle,$data);
 683 }
 684 if($_POST["21"] != "2.")
 685 {
 686 $data = "TimeShower:octetOniumColFac = ".$_POST["21"]."\n";
 687 fwrite($handle,$data);
 688 }
 689 if($_POST["22"] != "on")
 690 {
 691 $data = "TimeShower:QCDshower = ".$_POST["22"]."\n";
 692 fwrite($handle,$data);
 693 }
 694 if($_POST["23"] != "5")
 695 {
 696 $data = "TimeShower:nGluonToQuark = ".$_POST["23"]."\n";
 697 fwrite($handle,$data);
 698 }
 699 if($_POST["24"] != "on")
 700 {
 701 $data = "TimeShower:QEDshowerByQ = ".$_POST["24"]."\n";
 702 fwrite($handle,$data);
 703 }
 704 if($_POST["25"] != "on")
 705 {
 706 $data = "TimeShower:QEDshowerByL = ".$_POST["25"]."\n";
 707 fwrite($handle,$data);
 708 }
 709 if($_POST["26"] != "on")
 710 {
 711 $data = "TimeShower:QEDshowerByGamma = ".$_POST["26"]."\n";
 712 fwrite($handle,$data);
 713 }
 714 if($_POST["27"] != "5")
 715 {
 716 $data = "TimeShower:nGammaToQuark = ".$_POST["27"]."\n";
 717 fwrite($handle,$data);
 718 }
 719 if($_POST["28"] != "3")
 720 {
 721 $data = "TimeShower:nGammaToLepton = ".$_POST["28"]."\n";
 722 fwrite($handle,$data);
 723 }
 724 if($_POST["29"] != "on")
 725 {
 726 $data = "TimeShower:MEcorrections = ".$_POST["29"]."\n";
 727 fwrite($handle,$data);
 728 }
 729 if($_POST["30"] != "on")
 730 {
 731 $data = "TimeShower:MEafterFirst = ".$_POST["30"]."\n";
 732 fwrite($handle,$data);
 733 }
 734 if($_POST["31"] != "on")
 735 {
 736 $data = "TimeShower:phiPolAsym = ".$_POST["31"]."\n";
 737 fwrite($handle,$data);
 738 }
 739 if($_POST["32"] != "on")
 740 {
 741 $data = "TimeShower:recoilToColoured = ".$_POST["32"]."\n";
 742 fwrite($handle,$data);
 743 }
 744 fclose($handle);
 745 }
 746
 747 ?>
 748 </body>
 749 </html>
 750
 751 <!-- Copyright (C) 2012 Torbjorn Sjostrand -->