1 | <html> |
---|
2 | <head> |
---|
3 | <title>Spacelike 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='SpacelikeShowers.php'> |
---|
29 | |
---|
30 | <h2>Spacelike Showers</h2> |
---|
31 | |
---|
32 | The PYTHIA algorithm for spacelike initial-state showers is |
---|
33 | based on the article [<a href="Bibliography.php" target="page">Sjo05</a>], where a |
---|
34 | transverse-momentum-ordered backwards evolution scheme is introduced, |
---|
35 | with the extension to fully interleaved evolution covered in |
---|
36 | [<a href="Bibliography.php" target="page">Cor10a</a>]. |
---|
37 | This algorithm is a further development of the virtuality-ordered one |
---|
38 | presented in [<a href="Bibliography.php" target="page">Sj085</a>], with matching to first-order matrix |
---|
39 | element for <i>Z^0</i>, <i>W^+-</i> and Higgs (in the |
---|
40 | <i>m_t -> infinity</i> limit) production as introduced in |
---|
41 | [<a href="Bibliography.php" target="page">Miu99</a>]. |
---|
42 | |
---|
43 | <p/> |
---|
44 | The normal user is not expected to call <code>SpaceShower</code> |
---|
45 | directly, but only have it called from <code>Pythia</code>, |
---|
46 | via <code>PartonLevel</code>. Some of the parameters below, |
---|
47 | in particular <code>SpaceShower:alphaSvalue</code>, |
---|
48 | would be of interest for a tuning exercise, however. |
---|
49 | |
---|
50 | <h3>Main variables</h3> |
---|
51 | |
---|
52 | The maximum <i>pT</i> to be allowed in the shower evolution is |
---|
53 | related to the nature of the hard process itself. It involves a |
---|
54 | delicate balance between not doublecounting and not leaving any |
---|
55 | gaps in the coverage. The best procedure may depend on information |
---|
56 | only the user has: how the events were generated and mixed (e.g. with |
---|
57 | Les Houches Accord external input), and how they are intended to be |
---|
58 | used. Therefore a few options are available, with a sensible default |
---|
59 | behaviour. |
---|
60 | |
---|
61 | <br/><br/><table><tr><td><strong>SpaceShower:pTmaxMatch </td><td> (<code>default = <strong>0</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" checked="checked"><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. 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. <br/> |
---|
66 | <input type="radio" name="1" value="1"><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. This will simulate all possible event topologies, but may lead to doublecounting. (Also known as power showers.) <br/> |
---|
68 | <br/><b>Note 1:</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. |
---|
71 | <br/><b>Note 2:</b> Some processes contain matrix-element matching |
---|
72 | to the first emission; this is the case notably for single |
---|
73 | <ei>gamma^*/Z^0, W^+-</ei> and <ei>H^0</ei> production. Then default |
---|
74 | and option 2 give the correct result, while option 1 should never |
---|
75 | be used. |
---|
76 | |
---|
77 | <br/><br/><table><tr><td><strong>SpaceShower:pTmaxFudge </td><td></td><td> <input type="text" name="2" value="1.0" size="20"/> (<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 | |
---|
86 | |
---|
87 | <br/><br/><table><tr><td><strong>SpaceShower:pTmaxFudgeMPI </td><td></td><td> <input type="text" name="3" value="1.0" size="20"/> (<code>default = <strong>1.0</strong></code>; <code>minimum = 0.25</code>; <code>maximum = 2.0</code>)</td></tr></table> |
---|
88 | A multiplicative factor <i>f</i> such that |
---|
89 | <i>pT_max = f * pT_factorization</i>, as above, but here for the |
---|
90 | non-hardest interactions (when multiparton interactions are allowed). |
---|
91 | |
---|
92 | |
---|
93 | <br/><br/><table><tr><td><strong>SpaceShower:pTdampMatch </td><td> (<code>default = <strong>0</strong></code>; <code>minimum = 0</code>; <code>maximum = 2</code>)</td></tr></table> |
---|
94 | These options only take effect when a process is allowed to radiate up |
---|
95 | to the kinematical limit by the above <code>pTmaxMatch</code> choice, |
---|
96 | and no matrix-element corrections are available. Then, in many processes, |
---|
97 | the fall-off in <ei>pT</ei> will be too slow by one factor of <ei>pT^2</ei>. |
---|
98 | That is, while showers have an approximate <ei>dpT^2/pT^2</ei> shape, often |
---|
99 | it should become more like <ei>dpT^2/pT^4</ei> at <ei>pT</ei> values above |
---|
100 | the scale of the hard process. Whether this actually is the case |
---|
101 | depends on the particular process studied, e.g. if <ei>t</ei>-channel |
---|
102 | gluon exchange is likely to dominate. If so, the options below could |
---|
103 | provide a reasonable high-<ei>pT</ei> behaviour without requiring |
---|
104 | higher-order calculations. |
---|
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>SpaceShower:pTdampFudge </td><td></td><td> <input type="text" name="5" value="1.0" size="20"/> (<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>SpaceShower:alphaSvalue </td><td></td><td> <input type="text" name="6" value="0.137" size="20"/> (<code>default = <strong>0.137</strong></code>; <code>minimum = 0.06</code>; <code>maximum = 0.25</code>)</td></tr></table> |
---|
125 | The <i>alpha_strong</i> value at scale <code>M_Z^2</code>. |
---|
126 | Default value is picked equal to the one used in CTEQ 5L. |
---|
127 | |
---|
128 | |
---|
129 | <p/> |
---|
130 | The actual value is then regulated by the running to the scale |
---|
131 | <i>pT^2</i>, at which it is evaluated |
---|
132 | <br/><br/><table><tr><td><strong>SpaceShower:alphaSorder </td><td> (<code>default = <strong>1</strong></code>; <code>minimum = 0</code>; <code>maximum = 2</code>)</td></tr></table> |
---|
133 | Order at which <ei>alpha_strong</ei> runs, |
---|
134 | <br/> |
---|
135 | <input type="radio" name="7" value="0"><strong>0 </strong>: zeroth order, i.e. <ei>alpha_strong</ei> is kept fixed.<br/> |
---|
136 | <input type="radio" name="7" value="1" checked="checked"><strong>1 </strong>: first order, which is the normal value.<br/> |
---|
137 | <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/> |
---|
138 | |
---|
139 | <p/> |
---|
140 | QED radiation is regulated by the <i>alpha_electromagnetic</i> |
---|
141 | value at the <i>pT^2</i> scale of a branching. |
---|
142 | |
---|
143 | <br/><br/><table><tr><td><strong>SpaceShower:alphaEMorder </td><td> (<code>default = <strong>1</strong></code>; <code>minimum = -1</code>; <code>maximum = 1</code>)</td></tr></table> |
---|
144 | The running of <ei>alpha_em</ei>. |
---|
145 | <br/> |
---|
146 | <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/> |
---|
147 | <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/> |
---|
148 | <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/> |
---|
149 | |
---|
150 | <p/> |
---|
151 | The natural scale for couplings and PDFs is <i>pT^2</i>. To explore |
---|
152 | uncertainties it is possibly to vary around this value, however, in |
---|
153 | analogy with what can be done for |
---|
154 | <?php $filepath = $_GET["filepath"]; |
---|
155 | echo "<a href='CouplingsAndScales.php?filepath=".$filepath."' target='page'>";?>hard processes</a>. |
---|
156 | |
---|
157 | <br/><br/><table><tr><td><strong>SpaceShower:renormMultFac </td><td></td><td> <input type="text" name="9" value="1." size="20"/> (<code>default = <strong>1.</strong></code>; <code>minimum = 0.1</code>; <code>maximum = 10.</code>)</td></tr></table> |
---|
158 | The default <i>pT^2</i> renormalization scale is multiplied by |
---|
159 | this prefactor. For QCD this is equivalent to a change of |
---|
160 | <i>Lambda^2</i> in the opposite direction, i.e. to a change of |
---|
161 | <i>alpha_strong(M_Z^2)</i> (except that flavour thresholds |
---|
162 | remain at fixed scales). Below, when <i>pT^2 + pT_0^2</i> is used |
---|
163 | as scale, it is this whole expression that is multiplied by the prefactor. |
---|
164 | |
---|
165 | |
---|
166 | <br/><br/><table><tr><td><strong>SpaceShower:factorMultFac </td><td></td><td> <input type="text" name="10" value="1." size="20"/> (<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 | There are two complementary ways of regularizing the small-<i>pT</i> |
---|
173 | divergence, a sharp cutoff and a smooth dampening. These can be |
---|
174 | combined as desired but it makes sense to coordinate with how the |
---|
175 | same issue is handled in multiparton interactions. |
---|
176 | |
---|
177 | <br/><br/><strong>SpaceShower:samePTasMPI</strong> <input type="radio" name="11" value="on"><strong>On</strong> |
---|
178 | <input type="radio" name="11" value="off" checked="checked"><strong>Off</strong> |
---|
179 | (<code>default = <strong>off</strong></code>)<br/> |
---|
180 | Regularize the <i>pT -> 0</i> divergence using the same sharp cutoff |
---|
181 | and smooth dampening parameters as used to describe multiparton interactions. |
---|
182 | That is, the <code>MultipartonInteractions:pT0Ref</code>, |
---|
183 | <code>MultipartonInteractions:ecmRef</code>, |
---|
184 | <code>MultipartonInteractions:ecmPow</code> and |
---|
185 | <code>MultipartonInteractions:pTmin</code> parameters are used to regularize |
---|
186 | all ISR QCD radiation, rather than the corresponding parameters below. |
---|
187 | This is a sensible physics ansatz, based on the assumption that colour |
---|
188 | screening effects influence both MPI and ISR in the same way. Photon |
---|
189 | radiation is regularized separately in either case. |
---|
190 | <br/><b>Warning:</b> if a large <code>pT0</code> is picked for multiparton |
---|
191 | interactions, such that the integrated interaction cross section is |
---|
192 | below the nondiffractive inelastic one, this <code>pT0</code> will |
---|
193 | automatically be scaled down to cope. Information on such a rescaling |
---|
194 | does NOT propagate to <code>SpaceShower</code>, however. |
---|
195 | |
---|
196 | |
---|
197 | <p/> |
---|
198 | The actual <code>pT0</code> parameter used at a given CM energy scale, |
---|
199 | <i>ecmNow</i>, is obtained as |
---|
200 | <br/><i> |
---|
201 | pT0 = pT0(ecmNow) = pT0Ref * (ecmNow / ecmRef)^ecmPow |
---|
202 | </i><br/> |
---|
203 | where <i>pT0Ref</i>, <i>ecmRef</i> and <i>ecmPow</i> are the |
---|
204 | three parameters below. |
---|
205 | |
---|
206 | <br/><br/><table><tr><td><strong>SpaceShower:pT0Ref </td><td></td><td> <input type="text" name="12" value="2.0" size="20"/> (<code>default = <strong>2.0</strong></code>; <code>minimum = 0.5</code>; <code>maximum = 10.0</code>)</td></tr></table> |
---|
207 | Regularization of the divergence of the QCD emission probability for |
---|
208 | <i>pT -> 0</i> is obtained by a factor <i>pT^2 / (pT0^2 + pT^2)</i>, |
---|
209 | and by using an <i>alpha_s(pT0^2 + pT^2)</i>. An energy dependence |
---|
210 | of the <i>pT0</i> choice is introduced by the next two parameters, |
---|
211 | so that <i>pT0Ref</i> is the <i>pT0</i> value for the reference |
---|
212 | cm energy, <i>pT0Ref = pT0(ecmRef)</i>. |
---|
213 | |
---|
214 | |
---|
215 | <br/><br/><table><tr><td><strong>SpaceShower:ecmRef </td><td></td><td> <input type="text" name="13" value="1800.0" size="20"/> (<code>default = <strong>1800.0</strong></code>; <code>minimum = 1.</code>)</td></tr></table> |
---|
216 | The <i>ecmRef</i> reference energy scale introduced above. |
---|
217 | |
---|
218 | |
---|
219 | <br/><br/><table><tr><td><strong>SpaceShower:ecmPow </td><td></td><td> <input type="text" name="14" value="0.0" size="20"/> (<code>default = <strong>0.0</strong></code>; <code>minimum = 0.</code>; <code>maximum = 0.5</code>)</td></tr></table> |
---|
220 | The <i>ecmPow</i> energy rescaling pace introduced above. |
---|
221 | |
---|
222 | |
---|
223 | <br/><br/><table><tr><td><strong>SpaceShower:pTmin </td><td></td><td> <input type="text" name="15" value="0.2" size="20"/> (<code>default = <strong>0.2</strong></code>; <code>minimum = 0.1</code>; <code>maximum = 10.0</code>)</td></tr></table> |
---|
224 | Lower cutoff in <i>pT</i>, below which no further ISR branchings |
---|
225 | are allowed. Normally the <i>pT0</i> above would be used to |
---|
226 | provide the main regularization of the branching rate for |
---|
227 | <i>pT -> 0</i>, in which case <i>pTmin</i> is used mainly for |
---|
228 | technical reasons. It is possible, however, to set <i>pT0Ref = 0</i> |
---|
229 | and use <i>pTmin</i> to provide a step-function regularization, |
---|
230 | or to combine them in intermediate approaches. Currently <i>pTmin</i> |
---|
231 | is taken to be energy-independent. |
---|
232 | |
---|
233 | |
---|
234 | <br/><br/><table><tr><td><strong>SpaceShower:pTminChgQ </td><td></td><td> <input type="text" name="16" value="0.5" size="20"/> (<code>default = <strong>0.5</strong></code>; <code>minimum = 0.01</code>)</td></tr></table> |
---|
235 | Parton shower cut-off <i>pT</i> for photon coupling to a coloured |
---|
236 | particle. |
---|
237 | |
---|
238 | |
---|
239 | <br/><br/><table><tr><td><strong>SpaceShower:pTminChgL </td><td></td><td> <input type="text" name="17" value="0.0005" size="20"/> (<code>default = <strong>0.0005</strong></code>; <code>minimum = 0.0001</code>)</td></tr></table> |
---|
240 | Parton shower cut-off mass for pure QED branchings. |
---|
241 | Assumed smaller than (or equal to) <i>pTminChgQ</i>. |
---|
242 | |
---|
243 | |
---|
244 | <br/><br/><strong>SpaceShower:rapidityOrder</strong> <input type="radio" name="18" value="on"><strong>On</strong> |
---|
245 | <input type="radio" name="18" value="off" checked="checked"><strong>Off</strong> |
---|
246 | (<code>default = <strong>off</strong></code>)<br/> |
---|
247 | Force emissions, after the first, to be ordered in rapidity, |
---|
248 | i.e. in terms of decreasing angles in a backwards-evolution sense. |
---|
249 | Could be used to probe sensitivity to unordered emissions. |
---|
250 | Only affects QCD emissions. |
---|
251 | |
---|
252 | |
---|
253 | <h3>Further variables</h3> |
---|
254 | |
---|
255 | These should normally not be touched. Their only function is for |
---|
256 | cross-checks. |
---|
257 | |
---|
258 | <p/> |
---|
259 | There are three flags you can use to switch on or off selected |
---|
260 | branchings in the shower: |
---|
261 | |
---|
262 | <br/><br/><strong>SpaceShower:QCDshower</strong> <input type="radio" name="19" value="on" checked="checked"><strong>On</strong> |
---|
263 | <input type="radio" name="19" value="off"><strong>Off</strong> |
---|
264 | (<code>default = <strong>on</strong></code>)<br/> |
---|
265 | Allow a QCD shower; on/off = true/false. |
---|
266 | |
---|
267 | |
---|
268 | <br/><br/><strong>SpaceShower:QEDshowerByQ</strong> <input type="radio" name="20" value="on" checked="checked"><strong>On</strong> |
---|
269 | <input type="radio" name="20" value="off"><strong>Off</strong> |
---|
270 | (<code>default = <strong>on</strong></code>)<br/> |
---|
271 | Allow quarks to radiate photons; on/off = true/false. |
---|
272 | |
---|
273 | |
---|
274 | <br/><br/><strong>SpaceShower:QEDshowerByL</strong> <input type="radio" name="21" value="on" checked="checked"><strong>On</strong> |
---|
275 | <input type="radio" name="21" value="off"><strong>Off</strong> |
---|
276 | (<code>default = <strong>on</strong></code>)<br/> |
---|
277 | Allow leptons to radiate photons; on/off = true/false. |
---|
278 | |
---|
279 | |
---|
280 | <p/> |
---|
281 | There are some further possibilities to modify the shower: |
---|
282 | |
---|
283 | <br/><br/><strong>SpaceShower:MEcorrections</strong> <input type="radio" name="22" value="on" checked="checked"><strong>On</strong> |
---|
284 | <input type="radio" name="22" value="off"><strong>Off</strong> |
---|
285 | (<code>default = <strong>on</strong></code>)<br/> |
---|
286 | Use of matrix element corrections; on/off = true/false. |
---|
287 | |
---|
288 | |
---|
289 | <br/><br/><strong>SpaceShower:MEafterFirst</strong> <input type="radio" name="23" value="on" checked="checked"><strong>On</strong> |
---|
290 | <input type="radio" name="23" value="off"><strong>Off</strong> |
---|
291 | (<code>default = <strong>on</strong></code>)<br/> |
---|
292 | Use of matrix element corrections also after the first emission, |
---|
293 | for dipole ends of the same system that did not yet radiate. |
---|
294 | Only has a meaning if <code>MEcorrections</code> above is |
---|
295 | switched on. |
---|
296 | |
---|
297 | |
---|
298 | <br/><br/><strong>SpaceShower:phiPolAsym</strong> <input type="radio" name="24" value="on" checked="checked"><strong>On</strong> |
---|
299 | <input type="radio" name="24" value="off"><strong>Off</strong> |
---|
300 | (<code>default = <strong>on</strong></code>)<br/> |
---|
301 | Azimuthal asymmetry induced by gluon polarization; on/off = true/false. |
---|
302 | |
---|
303 | |
---|
304 | <br/><br/><strong>SpaceShower:phiIntAsym</strong> <input type="radio" name="25" value="on" checked="checked"><strong>On</strong> |
---|
305 | <input type="radio" name="25" value="off"><strong>Off</strong> |
---|
306 | (<code>default = <strong>on</strong></code>)<br/> |
---|
307 | Azimuthal asymmetry induced by interference; on/off = true/false. |
---|
308 | |
---|
309 | |
---|
310 | <br/><br/><table><tr><td><strong>SpaceShower:strengthIntAsym </td><td></td><td> <input type="text" name="26" value="0.7" size="20"/> (<code>default = <strong>0.7</strong></code>; <code>minimum = 0.</code>; <code>maximum = 0.9</code>)</td></tr></table> |
---|
311 | Size of asymmetry induced by interference. Natural value of order 0.5; |
---|
312 | expression would blow up for a value of 1. |
---|
313 | |
---|
314 | |
---|
315 | <br/><br/><table><tr><td><strong>SpaceShower:nQuarkIn </td><td></td><td> <input type="text" name="27" value="5" size="20"/> (<code>default = <strong>5</strong></code>; <code>minimum = 0</code>; <code>maximum = 5</code>)</td></tr></table> |
---|
316 | Number of allowed quark flavours in <i>g -> q qbar</i> branchings, |
---|
317 | when kinematically allowed, and thereby also in incoming beams. |
---|
318 | Changing it to 4 would forbid <i>g -> b bbar</i>, etc. |
---|
319 | |
---|
320 | |
---|
321 | <h3>Technical notes</h3> |
---|
322 | |
---|
323 | Almost everything is equivalent to the algorithm in [1]. Minor changes |
---|
324 | are as follows. |
---|
325 | <ul> |
---|
326 | <li> |
---|
327 | It is now possible to have a second-order running <i>alpha_s</i>, |
---|
328 | in addition to fixed or first-order running. |
---|
329 | </li> |
---|
330 | <li> |
---|
331 | The description of heavy flavour production in the threshold region |
---|
332 | has been modified, so as to be more forgiving about mismatches |
---|
333 | between the <i>c/b</i> masses used in Pythia relative to those |
---|
334 | used in a respective PDF parametrization. The basic idea is that, |
---|
335 | in the threshold region of a heavy quark <i>Q</i>, <i>Q = c/b</i>, |
---|
336 | the effect of subsequent <i>Q -> Q g</i> branchings is negligible. |
---|
337 | If so, then |
---|
338 | <br/><i> |
---|
339 | f_Q(x, pT2) = integral_mQ2^pT2 dpT'2/pT'2 * alpha_s(pT'2)/2pi |
---|
340 | * integral P(z) g(x', pT'2) delta(x - z x') |
---|
341 | </i><br/> |
---|
342 | so use this to select the <i>pT2</i> of the <i>g -> Q Qbar</i> |
---|
343 | branching. In the old formalism the same kind of behaviour should |
---|
344 | be obtained, but by a cancellation of a <i>1/f_Q</i> that diverges |
---|
345 | at the theshold and a Sudakov that vanishes. |
---|
346 | <br/> |
---|
347 | The strategy therefore is that, once <i>pT2 < f * mQ2</i>, with |
---|
348 | <i>f</i> a parameter of the order of 2, a <i>pT2</i> is chosen |
---|
349 | like <i>dpT2/pT2</i> between <i>mQ2</i> and <i>f * mQ2</i>, a |
---|
350 | nd a <i>z</i> flat in the allowed range. Thereafter acceptance |
---|
351 | is based on the product of three factors, representing the running |
---|
352 | of <i>alpha_strong</i>, the splitting kernel (including the mass term) |
---|
353 | and the gluon density weight. At failure, a new <i>pT2</i> is chosen |
---|
354 | in the same range, i.e. is not required to be lower since no Sudakov |
---|
355 | is involved. |
---|
356 | </li> |
---|
357 | <li> |
---|
358 | The QED algorithm now allows for hadron beams with non-zero photon |
---|
359 | content. The backwards-evolution of a photon in a hadron is identical |
---|
360 | to that of a gluon, with <i>CF -> eq^2</i> and <i>CA -> 0</i>. |
---|
361 | Note that this will only work in conjunction with |
---|
362 | parton distribution that explicitly include photons as part of the |
---|
363 | hadron structure (such as the MRST2004qed set). Since Pythia's |
---|
364 | internal sets do not allow for photon content in hadrons, it is thus |
---|
365 | necessary to use the LHAPDF interface to make use of this feature. The |
---|
366 | possibility of a fermion backwards-evolving to a photon has not yet |
---|
367 | been included, nor has photon backwards-evolution in lepton beams. |
---|
368 | </li> |
---|
369 | </ul> |
---|
370 | |
---|
371 | <input type="hidden" name="saved" value="1"/> |
---|
372 | |
---|
373 | <?php |
---|
374 | echo "<input type='hidden' name='filepath' value='".$_GET["filepath"]."'/>"?> |
---|
375 | |
---|
376 | <table width="100%"><tr><td align="right"><input type="submit" value="Save Settings" /></td></tr></table> |
---|
377 | </form> |
---|
378 | |
---|
379 | <?php |
---|
380 | |
---|
381 | if($_POST["saved"] == 1) |
---|
382 | { |
---|
383 | $filepath = $_POST["filepath"]; |
---|
384 | $handle = fopen($filepath, 'a'); |
---|
385 | |
---|
386 | if($_POST["1"] != "0") |
---|
387 | { |
---|
388 | $data = "SpaceShower:pTmaxMatch = ".$_POST["1"]."\n"; |
---|
389 | fwrite($handle,$data); |
---|
390 | } |
---|
391 | if($_POST["2"] != "1.0") |
---|
392 | { |
---|
393 | $data = "SpaceShower:pTmaxFudge = ".$_POST["2"]."\n"; |
---|
394 | fwrite($handle,$data); |
---|
395 | } |
---|
396 | if($_POST["3"] != "1.0") |
---|
397 | { |
---|
398 | $data = "SpaceShower:pTmaxFudgeMPI = ".$_POST["3"]."\n"; |
---|
399 | fwrite($handle,$data); |
---|
400 | } |
---|
401 | if($_POST["4"] != "0") |
---|
402 | { |
---|
403 | $data = "SpaceShower:pTdampMatch = ".$_POST["4"]."\n"; |
---|
404 | fwrite($handle,$data); |
---|
405 | } |
---|
406 | if($_POST["5"] != "1.0") |
---|
407 | { |
---|
408 | $data = "SpaceShower:pTdampFudge = ".$_POST["5"]."\n"; |
---|
409 | fwrite($handle,$data); |
---|
410 | } |
---|
411 | if($_POST["6"] != "0.137") |
---|
412 | { |
---|
413 | $data = "SpaceShower:alphaSvalue = ".$_POST["6"]."\n"; |
---|
414 | fwrite($handle,$data); |
---|
415 | } |
---|
416 | if($_POST["7"] != "1") |
---|
417 | { |
---|
418 | $data = "SpaceShower:alphaSorder = ".$_POST["7"]."\n"; |
---|
419 | fwrite($handle,$data); |
---|
420 | } |
---|
421 | if($_POST["8"] != "1") |
---|
422 | { |
---|
423 | $data = "SpaceShower:alphaEMorder = ".$_POST["8"]."\n"; |
---|
424 | fwrite($handle,$data); |
---|
425 | } |
---|
426 | if($_POST["9"] != "1.") |
---|
427 | { |
---|
428 | $data = "SpaceShower:renormMultFac = ".$_POST["9"]."\n"; |
---|
429 | fwrite($handle,$data); |
---|
430 | } |
---|
431 | if($_POST["10"] != "1.") |
---|
432 | { |
---|
433 | $data = "SpaceShower:factorMultFac = ".$_POST["10"]."\n"; |
---|
434 | fwrite($handle,$data); |
---|
435 | } |
---|
436 | if($_POST["11"] != "off") |
---|
437 | { |
---|
438 | $data = "SpaceShower:samePTasMPI = ".$_POST["11"]."\n"; |
---|
439 | fwrite($handle,$data); |
---|
440 | } |
---|
441 | if($_POST["12"] != "2.0") |
---|
442 | { |
---|
443 | $data = "SpaceShower:pT0Ref = ".$_POST["12"]."\n"; |
---|
444 | fwrite($handle,$data); |
---|
445 | } |
---|
446 | if($_POST["13"] != "1800.0") |
---|
447 | { |
---|
448 | $data = "SpaceShower:ecmRef = ".$_POST["13"]."\n"; |
---|
449 | fwrite($handle,$data); |
---|
450 | } |
---|
451 | if($_POST["14"] != "0.0") |
---|
452 | { |
---|
453 | $data = "SpaceShower:ecmPow = ".$_POST["14"]."\n"; |
---|
454 | fwrite($handle,$data); |
---|
455 | } |
---|
456 | if($_POST["15"] != "0.2") |
---|
457 | { |
---|
458 | $data = "SpaceShower:pTmin = ".$_POST["15"]."\n"; |
---|
459 | fwrite($handle,$data); |
---|
460 | } |
---|
461 | if($_POST["16"] != "0.5") |
---|
462 | { |
---|
463 | $data = "SpaceShower:pTminChgQ = ".$_POST["16"]."\n"; |
---|
464 | fwrite($handle,$data); |
---|
465 | } |
---|
466 | if($_POST["17"] != "0.0005") |
---|
467 | { |
---|
468 | $data = "SpaceShower:pTminChgL = ".$_POST["17"]."\n"; |
---|
469 | fwrite($handle,$data); |
---|
470 | } |
---|
471 | if($_POST["18"] != "off") |
---|
472 | { |
---|
473 | $data = "SpaceShower:rapidityOrder = ".$_POST["18"]."\n"; |
---|
474 | fwrite($handle,$data); |
---|
475 | } |
---|
476 | if($_POST["19"] != "on") |
---|
477 | { |
---|
478 | $data = "SpaceShower:QCDshower = ".$_POST["19"]."\n"; |
---|
479 | fwrite($handle,$data); |
---|
480 | } |
---|
481 | if($_POST["20"] != "on") |
---|
482 | { |
---|
483 | $data = "SpaceShower:QEDshowerByQ = ".$_POST["20"]."\n"; |
---|
484 | fwrite($handle,$data); |
---|
485 | } |
---|
486 | if($_POST["21"] != "on") |
---|
487 | { |
---|
488 | $data = "SpaceShower:QEDshowerByL = ".$_POST["21"]."\n"; |
---|
489 | fwrite($handle,$data); |
---|
490 | } |
---|
491 | if($_POST["22"] != "on") |
---|
492 | { |
---|
493 | $data = "SpaceShower:MEcorrections = ".$_POST["22"]."\n"; |
---|
494 | fwrite($handle,$data); |
---|
495 | } |
---|
496 | if($_POST["23"] != "on") |
---|
497 | { |
---|
498 | $data = "SpaceShower:MEafterFirst = ".$_POST["23"]."\n"; |
---|
499 | fwrite($handle,$data); |
---|
500 | } |
---|
501 | if($_POST["24"] != "on") |
---|
502 | { |
---|
503 | $data = "SpaceShower:phiPolAsym = ".$_POST["24"]."\n"; |
---|
504 | fwrite($handle,$data); |
---|
505 | } |
---|
506 | if($_POST["25"] != "on") |
---|
507 | { |
---|
508 | $data = "SpaceShower:phiIntAsym = ".$_POST["25"]."\n"; |
---|
509 | fwrite($handle,$data); |
---|
510 | } |
---|
511 | if($_POST["26"] != "0.7") |
---|
512 | { |
---|
513 | $data = "SpaceShower:strengthIntAsym = ".$_POST["26"]."\n"; |
---|
514 | fwrite($handle,$data); |
---|
515 | } |
---|
516 | if($_POST["27"] != "5") |
---|
517 | { |
---|
518 | $data = "SpaceShower:nQuarkIn = ".$_POST["27"]."\n"; |
---|
519 | fwrite($handle,$data); |
---|
520 | } |
---|
521 | fclose($handle); |
---|
522 | } |
---|
523 | |
---|
524 | ?> |
---|
525 | </body> |
---|
526 | </html> |
---|
527 | |
---|
528 | <!-- Copyright (C) 2012 Torbjorn Sjostrand --> |
---|
529 | |
---|