]>
Commit | Line | Data |
---|---|---|
c6b60c38 | 1 | <chapter name="Multiparton Interactions"> |
2 | ||
3 | <h2>Multiparton Interactions</h2> | |
4 | ||
5 | The starting point for the multiparton interactions physics scenario in | |
6 | PYTHIA is provided by <ref>Sjo87</ref>. Recent developments have | |
7 | included a more careful study of flavour and colour correlations, | |
8 | junction topologies and the relationship to beam remnants | |
9 | <ref>Sjo04</ref>, interleaving with initial-state radiation | |
10 | <ref>Sjo05</ref>, making use of transverse-momentum-ordered | |
11 | initial- and final-state showers, with the extension to fully | |
12 | interleaved evolution covered in <ref>Cor10a</ref>. A framework to | |
13 | handle rescattering is described in <ref>Cor09</ref>. | |
14 | ||
15 | <p/> | |
16 | A big unsolved issue is how the colour of all these subsystems is | |
17 | correlated. For sure there is a correlation coming from the colour | |
18 | singlet nature of the incoming beams, but in addition final-state | |
19 | colour rearrangements may change the picture. Indeed such extra | |
20 | effects appear necessary to describe data, e.g. on | |
21 | <ei><pT>(n_ch)</ei>. A simple implementation of colour | |
22 | rearrangement is found as part of the | |
23 | <aloc href="BeamRemnants">beam remnants</aloc> description. | |
24 | ||
25 | <h3>Main variables</h3> | |
26 | ||
27 | <h4>Matching to hard process</h4> | |
28 | ||
29 | The maximum <ei>pT</ei> to be allowed for multiparton interactions is | |
30 | related to the nature of the hard process itself. It involves a | |
31 | delicate balance between not double-counting and not leaving any | |
32 | gaps in the coverage. The best procedure may depend on information | |
33 | only the user has: how the events were generated and mixed (e.g. with | |
34 | Les Houches Accord external input), and how they are intended to be | |
35 | used. Therefore a few options are available, with a sensible default | |
36 | behaviour. | |
37 | <modepick name="MultipartonInteractions:pTmaxMatch" default="0" min="0" | |
38 | max="2"> | |
39 | Way in which the maximum scale for multiparton interactions is set | |
40 | to match the scale of the hard process itself. | |
41 | <option value="0"><b>(i)</b> if the final state of the hard process | |
42 | (not counting subsequent resonance decays) contains only quarks | |
43 | (<ei>u, d, s, c, b</ei>), gluons and photons then <ei>pT_max</ei> | |
44 | is chosen to be the factorization scale for internal processes | |
45 | and the <code>scale</code> value for Les Houches input; | |
46 | <b>(ii)</b> if not, interactions are allowed to go all the way up | |
47 | to the kinematical limit. | |
48 | The reasoning is that the former kind of processes are generated by | |
49 | the multiparton-interactions machinery and so would double-count hard | |
50 | processes if allowed to overlap the same <ei>pT</ei> range, | |
51 | while no such danger exists in the latter case. | |
52 | </option> | |
53 | <option value="1">always use the factorization scale for an internal | |
54 | process and the <code>scale</code> value for Les Houches input, | |
55 | i.e. the lower value. This should avoid double-counting, but | |
56 | may leave out some interactions that ought to have been simulated. | |
57 | </option> | |
58 | <option value="2">always allow multiparton interactions up to the | |
59 | kinematical limit. This will simulate all possible event topologies, | |
60 | but may lead to double-counting. | |
61 | </option> | |
62 | <note>Note:</note> If a "second hard" process is present, the two | |
63 | are analyzed separately for the default 0 option. It is enough that | |
64 | one of them only consists of quarks, gluons and photons to restrict | |
65 | the <ei>pT</ei> range. The maximum for MPI is then set by the hard | |
66 | interaction with lowest scale. | |
67 | </modepick> | |
68 | ||
69 | <h4>Cross-section parameters</h4> | |
70 | ||
71 | The rate of interactions is determined by | |
72 | <parm name="MultipartonInteractions:alphaSvalue" default="0.127" | |
73 | min="0.06" max="0.25"> | |
74 | The value of <ei>alpha_strong</ei> at <ei>m_Z</ei>. Default value is | |
75 | picked equal to the one used in CTEQ 5L. | |
76 | </parm> | |
77 | ||
78 | <p/> | |
79 | The actual value is then regulated by the running to the scale | |
80 | <ei>pT^2</ei>, at which it is evaluated | |
81 | <modepick name="MultipartonInteractions:alphaSorder" default="1" | |
82 | min="0" max="2"> | |
83 | The order at which <ei>alpha_strong</ei> runs at scales away from | |
84 | <ei>m_Z</ei>. | |
85 | <option value="0">zeroth order, i.e. <ei>alpha_strong</ei> is kept | |
86 | fixed.</option> | |
87 | <option value="1">first order, which is the normal value.</option> | |
88 | <option value="2">second order. Since other parts of the code do | |
89 | not go to second order there is no strong reason to use this option, | |
90 | but there is also nothing wrong with it.</option> | |
91 | </modepick> | |
92 | ||
93 | <p/> | |
94 | QED interactions are regulated by the <ei>alpha_electromagnetic</ei> | |
95 | value at the <ei>pT^2</ei> scale of an interaction. | |
96 | ||
97 | <modepick name="MultipartonInteractions:alphaEMorder" default="1" | |
98 | min="-1" max="1"> | |
99 | The running of <ei>alpha_em</ei> used in hard processes. | |
100 | <option value="1">first-order running, constrained to agree with | |
101 | <code>StandardModel:alphaEMmZ</code> at the <ei>Z^0</ei> mass. | |
102 | </option> | |
103 | <option value="0">zeroth order, i.e. <ei>alpha_em</ei> is kept | |
104 | fixed at its value at vanishing momentum transfer.</option> | |
105 | <option value="-1">zeroth order, i.e. <ei>alpha_em</ei> is kept | |
106 | fixed, but at <code>StandardModel:alphaEMmZ</code>, i.e. its value | |
107 | at the <ei>Z^0</ei> mass. | |
108 | </option> | |
109 | </modepick> | |
110 | ||
111 | <p/> | |
112 | Note that the choices of <ei>alpha_strong</ei> and <ei>alpha_em</ei> | |
113 | made here override the ones implemented in the normal process machinery, | |
114 | but only for the interactions generated by the | |
115 | <code>MultipartonInteractions</code> class. | |
116 | ||
117 | <p/> | |
118 | In addition there is the possibility of a global rescaling of | |
119 | cross sections (which could not easily be accommodated by a | |
120 | changed <ei>alpha_strong</ei>, since <ei>alpha_strong</ei> runs) | |
121 | <parm name="MultipartonInteractions:Kfactor" default="1.0" min="0.5" | |
734883cb | 122 | max="32.0"> |
c6b60c38 | 123 | Multiply all cross sections by this fix factor. |
124 | </parm> | |
125 | ||
126 | <p/> | |
127 | The processes used to generate multiparton interactions form a subset | |
128 | of the standard library of hard processes. The input is slightly | |
129 | different from the standard hard-process machinery, however, | |
130 | since incoming flavours, the <ei>alpha_strong</ei> value and most | |
131 | of the kinematics are already fixed when the process is called. | |
132 | It is possible to regulate the set of processes that are included in the | |
133 | multiparton-interactions framework. | |
134 | ||
135 | <modepick name="MultipartonInteractions:processLevel" default="3" | |
136 | min="0" max="3"> | |
137 | Set of processes included in the machinery. | |
138 | <option value="0">only the simplest <ei>2 -> 2</ei> QCD processes | |
139 | between quarks and gluons, giving no new flavours, i.e. dominated by | |
140 | <ei>t</ei>-channel gluon exchange.</option> | |
141 | <option value="1">also <ei>2 -> 2</ei> QCD processes giving new flavours | |
142 | (including charm and bottom), i.e. proceeding through <ei>s</ei>-channel | |
143 | gluon exchange.</option> | |
144 | <option value="2">also <ei>2 -> 2</ei> processes involving one or two | |
145 | photons in the final state, <ei>s</ei>-channel <ei>gamma</ei> | |
146 | boson exchange and <ei>t</ei>-channel <ei>gamma/Z^0/W^+-</ei> | |
147 | boson exchange.</option> | |
148 | <option value="3">also charmonium and bottomonium production, via | |
149 | colour singlet and colour octet channels.</option> | |
150 | </modepick> | |
151 | ||
152 | <h4>Cross-section regularization</h4> | |
153 | ||
154 | There are two complementary ways of regularizing the small-<ei>pT</ei> | |
155 | divergence, a sharp cutoff and a smooth dampening. These can be | |
156 | combined as desired, but it makes sense to coordinate with how the | |
157 | same issue is handled in <aloc href="SpacelikeShowers">spacelike | |
158 | showers</aloc>. Actually, by default, the parameters defined here are | |
159 | used also for the spacelike showers, but this can be overridden. | |
160 | ||
161 | <p/> | |
162 | Regularization of the divergence of the QCD cross section for | |
163 | <ei>pT -> 0</ei> is obtained by a factor <ei>pT^4 / (pT0^2 + pT^2)^2</ei>, | |
164 | and by using an <ei>alpha_s(pT0^2 + pT^2)</ei>. An energy dependence | |
165 | of the <ei>pT0</ei> choice is introduced by two further parameters, | |
166 | so that <ei>pT0Ref</ei> is the <ei>pT0</ei> value for the reference | |
167 | CM energy, <ei>pT0Ref = pT0(ecmRef)</ei>. | |
168 | <note>Warning:</note> if a large <ei>pT0</ei> is picked for multiparton | |
169 | interactions, such that the integrated interaction cross section is | |
170 | below the nondiffractive inelastic one, this <ei>pT0</ei> will | |
171 | automatically be scaled down to cope. | |
172 | ||
173 | <p/> | |
174 | The actual <ei>pT0</ei> parameter used at a given CM energy scale, | |
175 | <ei>ecmNow</ei>, is obtained as | |
176 | <eq> | |
177 | pT0 = pT0(ecmNow) = pT0Ref * (ecmNow / ecmRef)^ecmPow | |
178 | </eq> | |
179 | where <ei>pT0Ref</ei>, <ei>ecmRef</ei> and <ei>ecmPow</ei> are the | |
180 | three parameters below. | |
181 | ||
182 | <parm name="MultipartonInteractions:pT0Ref" default="2.15" min="0.5" | |
183 | max="10.0"> | |
184 | The <ei>pT0Ref</ei> scale in the above formula. | |
185 | <note>Note:</note> <ei>pT0Ref</ei> is one of the key parameters in a | |
186 | complete PYTHIA tune. Its value is intimately tied to a number of other | |
187 | choices, such as that of colour flow description, so unfortunately it is | |
188 | difficult to give an independent meaning to <ei>pT0Ref</ei>. | |
189 | </parm> | |
190 | ||
191 | <parm name="MultipartonInteractions:ecmRef" default="1800.0" min="1."> | |
192 | The <ei>ecmRef</ei> reference energy scale introduced above. | |
193 | </parm> | |
194 | ||
195 | <parm name="MultipartonInteractions:ecmPow" default="0.24" min="0.0" | |
196 | max="0.5"> | |
197 | The <ei>ecmPow</ei> energy rescaling pace introduced above. | |
198 | </parm> | |
199 | ||
200 | <p/> | |
201 | Alternatively, or in combination, a sharp cut can be used. | |
202 | <parm name="MultipartonInteractions:pTmin" default="0.2" min="0.1" | |
203 | max="10.0"> | |
204 | Lower cutoff in <ei>pT</ei>, below which no further interactions | |
205 | are allowed. Normally <ei>pT0</ei> above would be used to provide | |
206 | the main regularization of the cross section for <ei>pT -> 0</ei>, | |
207 | in which case <ei>pTmin</ei> is used mainly for technical reasons. | |
208 | It is possible, however, to set <ei>pT0Ref = 0</ei> and use | |
209 | <ei>pTmin</ei> to provide a step-function regularization, or to | |
210 | combine them in intermediate approaches. Currently <ei>pTmin</ei> | |
211 | is taken to be energy-independent. | |
212 | </parm> | |
213 | ||
214 | <p/> | |
215 | Gösta Gustafson has proposed (private communication, unpublished) | |
216 | that the amount of screening, as encapsulated in the <ei>pT0</ei> | |
217 | parameter, fluctuates from one event to the next. Specifically, | |
218 | high-activity event are more likely to lead to interactions at large | |
219 | <ei>pT</ei> scales, but the high activity simultaneously leads to a | |
220 | larger screening of interactions at smaller <ei>pT</ei>. Such a scenario | |
221 | can approximately be simulated by scaling up the <ei>pT0</ei> by a | |
222 | factor <ei>sqrt(n)</ei>, where <ei>n</ei> is the number of interactions | |
223 | considered so far, including the current one. That is, for the first | |
224 | interaction the dampening factor is <ei>pT^4 / (pT0^2 + pT^2)^2</ei>, | |
225 | for the second <ei>pT^4 / (2 pT0^2 + pT^2)^2</ei>, for the third | |
226 | <ei>pT^4 / (3 pT0^2 + pT^2)^2</ei>, and so on. Optionally the scheme | |
227 | may also be applied to ISR emissions. For simplicity the same | |
228 | <ei>alpha_s(pT0^2 + pT^2)</ei> is used throughout. Note that, in this | |
229 | scenario the <ei>pT0</ei> scale must be lower than in the normal case | |
230 | to begin with, since it later is increased back up. Also note that the | |
231 | idea with this scenario is to propose an alternative to colour | |
232 | reconnection to understand the rise of <ei><pT>(n_ch)</ei>, | |
233 | so that the amount of colour reconnection should be reduced. | |
234 | <modepick name="MultipartonInteractions:enhanceScreening" default="0" | |
235 | min="0" max="2"> | |
236 | Choice to activate the above screening scenario, i.e. an increasing | |
237 | effective <ei>pT0</ei> for consecutive interactions. | |
238 | <option value="0">No activity-dependent screening, i.e. <ei>pT0</ei> | |
239 | is fixed.</option> | |
240 | <option value="1">The <ei>pT0</ei> scale is increased as a function | |
241 | of the number of MPI's, as explained above. ISR is not affected, | |
242 | but note that, if <code>SpaceShower:samePTasMPI</code> is on, | |
243 | then <code>MultipartonInteractions:pT0Ref</code> is used also for ISR, | |
244 | which may or may not be desirable. | |
245 | </option> | |
246 | <option value="2">Both MPI and ISR influence and are influenced by the | |
247 | screening. That is, the dampening is reduced based on the total number | |
248 | of MPI and ISR steps considered so far, including the current one. | |
249 | This dampening is implemented both for MPI and for ISR emissions, | |
250 | for the latter provided that <code>SpaceShower:samePTasMPI</code> is on | |
251 | (default). | |
252 | </option> | |
253 | </modepick> | |
254 | ||
255 | <h4>Impact-parameter dependence</h4> | |
256 | ||
257 | The choice of impact-parameter dependence is regulated by several | |
258 | parameters. The ones listed here refer to nondiffractive topologies | |
259 | only, while their equivalents for diffractive events are put in the | |
260 | <aloc href="Diffraction">Diffraction</aloc> description. Note that | |
261 | there is currently no <code>bProfile = 4</code> option for diffraction. | |
262 | Other parameters are assumed to agree between diffractive and | |
263 | nondiffractive topologies. | |
264 | ||
265 | <modepick name="MultipartonInteractions:bProfile" default="1" | |
266 | min="0" max="4"> | |
267 | Choice of impact parameter profile for the incoming hadron beams. | |
268 | <option value="0">no impact parameter dependence at all.</option> | |
269 | <option value="1">a simple Gaussian matter distribution; | |
270 | no free parameters.</option> | |
271 | <option value="2">a double Gaussian matter distribution, | |
272 | with the two free parameters <ei>coreRadius</ei> and | |
273 | <ei>coreFraction</ei>.</option> | |
274 | <option value="3">an overlap function, i.e. the convolution of | |
275 | the matter distributions of the two incoming hadrons, of the form | |
276 | <ei>exp(- b^expPow)</ei>, where <ei>expPow</ei> is a free | |
277 | parameter.</option> | |
278 | <option value="4">a Gaussian matter distribution with a width | |
279 | that varies according to the selected <ei>x</ei> value of an interaction, | |
280 | <ei>1. + a1 log (1 / x)</ei>, where <ei>a1</ei> is a free parameter. | |
281 | Note that once <ei>b</ei> has been selected for the hard process, | |
282 | it remains fixed for the remainder of the evolution. | |
283 | </option> | |
284 | </modepick> | |
285 | ||
286 | <parm name="MultipartonInteractions:coreRadius" default="0.4" | |
287 | min="0.1" max="1."> | |
288 | When assuming a double Gaussian matter profile, <ei>bProfile = 2</ei>, | |
289 | the inner core is assumed to have a radius that is a factor | |
290 | <ei>coreRadius</ei> smaller than the rest. | |
291 | </parm> | |
292 | ||
293 | <parm name="MultipartonInteractions:coreFraction" default="0.5" | |
294 | min="0." max="1."> | |
295 | When assuming a double Gaussian matter profile, <ei>bProfile = 2</ei>, | |
296 | the inner core is assumed to have a fraction <ei>coreFraction</ei> | |
297 | of the matter content of the hadron. | |
298 | </parm> | |
299 | ||
300 | <parm name="MultipartonInteractions:expPow" default="1." min="0.4" max="10."> | |
301 | When <ei>bProfile = 3</ei> it gives the power of the assumed overlap | |
302 | shape <ei>exp(- b^expPow)</ei>. Default corresponds to a simple | |
303 | exponential drop, which is not too dissimilar from the overlap | |
304 | obtained with the standard double Gaussian parameters. For | |
305 | <ei>expPow = 2</ei> we reduce to the simple Gaussian, <ei>bProfile = 1</ei>, | |
306 | and for <ei>expPow -> infinity</ei> to no impact parameter dependence | |
307 | at all, <ei>bProfile = 0</ei>. For small <ei>expPow</ei> the program | |
308 | becomes slow and unstable, so the min limit must be respected. | |
309 | </parm> | |
310 | ||
311 | <parm name="MultipartonInteractions:a1" default="0.15" min="0." max="2."> | |
312 | When <ei>bProfile = 4</ei>, this gives the <ei>a1</ei> constant in the | |
313 | Gaussian width. When <ei>a1 = 0.</ei>, this reduces back to the single | |
314 | Gaussian case. | |
315 | </parm> | |
316 | ||
317 | <modepick name="MultipartonInteractions:bSelScale" default="1" | |
318 | min="1" max="3"> | |
319 | The selection of impact parameter is related to the scale of the | |
320 | hard process: the harder this scale is, the more central the collision. | |
321 | In practice this centrality saturates quickly, however, and beyond | |
322 | a scale of roughly 20 GeV very little changes. (The relevant | |
323 | quantity is that the QCD jet cross section above the scale should be | |
324 | a tiny fraction of the total cross section.) In <ei>2 -> 1</ei> and | |
325 | <ei>2 -> 2</ei> processes traditional scale choices work fine, but | |
326 | ambiguities arise for higher multiplicities, in particular when the | |
327 | scale is used for matching between the multiparton matrix elements | |
328 | and parton showers. Then the event scale may | |
329 | be chosen as that of a very low-<ei>pT</ei> parton, i.e. suggesting a | |
330 | peripheral collision, while the much harder other partons instead | |
331 | would favour a central collision. Therefore the default here is to | |
332 | override whatever scale value have been read in from an LHEF, say. | |
333 | Notice that the scale used here is decoupled from the maximum scale | |
334 | for MPIs (<code>MultipartonInteractions:pTmaxMatch</code>). | |
335 | <option value="1"> | |
336 | Use the mass for a <ei>2 -> 1</ei> process. | |
337 | For <ei>2 -> n, n > 1</ei> processes order the particles in | |
338 | falling <ei>mmT = m + mT</ei> and then let the scale be | |
339 | <ei>(mmT_1 + mmT_2)/2 + mmT_3/3 + mmT_4/4 + ... + mmT_n/n</ei>. | |
340 | This is constructed always to be above <ei>m1</ei>, and to assign | |
341 | decreasing importance to softer particles that are less likely to | |
342 | be associated with the hard process.</option> | |
343 | <option value="2">Use the <code>scale</code> parameter of the event. | |
344 | </option> | |
345 | <option value="3">use the same scale as chosen by the rules for | |
346 | <code>MultipartonInteractions:pTmaxMatch</code>.</option> | |
347 | </modepick> | |
348 | ||
349 | <h4>Rescattering</h4> | |
350 | ||
351 | It is possible that a parton may rescatter, i.e. undergo a further | |
352 | interaction subsequent to the first one. The machinery to model this | |
353 | kind of physics has only recently become fully operational | |
354 | <ref>Cor09</ref>, and is therefore not yet so well explored. | |
355 | ||
356 | <p/> | |
357 | The rescattering framework has ties with other parts of the program, | |
358 | notably with the <aloc href="BeamRemnants">beam remnants</aloc>. | |
359 | ||
360 | <flag name="MultipartonInteractions:allowRescatter" default="off"> | |
361 | Switch to allow rescattering of partons; on/off = true/false.<br/> | |
362 | <b>Note:</b> the rescattering framework has not yet been implemented | |
363 | for the <code>MultipartonInteractions:bProfile = 4</code> option, | |
364 | and can therefore not be switched on in that case. | |
365 | <b>Warning:</b> use with caution since machinery is still not | |
366 | so well tested. | |
367 | </flag> | |
368 | ||
369 | <flag name="MultipartonInteractions:allowDoubleRescatter" default="off"> | |
370 | Switch to allow rescattering of partons, where both incoming partons | |
371 | have already rescattered; on/off = true/false. Is only used if | |
372 | <code>MultipartonInteractions:allowRescatter</code> is switched on.<br/> | |
373 | <b>Warning:</b> currently there is no complete implementation that | |
374 | combines it with shower evolution, so you must use | |
375 | <code>PartonLevel:ISR = off</code> and <code>PartonLevel:FSR = off</code>. | |
376 | If not, a warning will be issued and double rescattering will not be | |
377 | simulated. The rate also comes out to be much lower than for single | |
378 | rescattering, so to first approximation it can be neglected. | |
379 | </flag> | |
380 | ||
381 | <modepick name="MultipartonInteractions:rescatterMode" default="0" | |
382 | min="0" max="4"> | |
383 | Selection of which partons rescatter against unscattered partons | |
384 | from the incoming beams A and B, based on their rapidity value | |
385 | <ei>y</ei> in the collision rest frame. Here <ei>ySep</ei> is | |
386 | shorthand for <code>MultipartonInteractions:ySepRescatter</code> and | |
387 | <ei>deltaY</ei> for <code>MultipartonInteractions:deltaYRescatter</code>, | |
388 | defined below. The description is symmetric between the two beams, | |
389 | so only one case is described below. | |
390 | <option value="0">only scattered partons with <ei>y > 0</ei> | |
391 | can collide with unscattered partons from beam B.</option> | |
392 | <option value="1">only scattered partons with <ei>y > ySep</ei> | |
393 | can collide with unscattered partons from beam B.</option> | |
394 | <option value="2">the probability for a scattered parton to be considered | |
395 | as a potential rescatterer against unscattered partons in beam B increases | |
396 | linearly from zero at <ei>y = ySep - deltaY</ei> to unity at | |
397 | <ei>y = ySep + deltaY</ei>.</option> | |
398 | <option value="3">the probability for a scattered parton to be considered | |
399 | as a potential rescatterer against unscattered partons in beam B increases | |
400 | with <ei>y</ei> according to | |
401 | <ei>(1/2) * (1 + tanh( (y - ySep) / deltaY))</ei>.</option> | |
402 | <option value="4">all partons are potential rescatterers against both | |
403 | beams.</option> | |
404 | </modepick> | |
405 | ||
406 | <parm name="MultipartonInteractions:ySepRescatter" default="0."> | |
407 | used for some of the <code>MultipartonInteractions:rescatterMode</code> | |
408 | options above, as the rapidity for which a scattered parton has a 50% | |
409 | probability to be considered as a potential rescatterer. | |
410 | A <ei>ySep > 0</ei> generally implies that some central partons cannot | |
411 | rescatter at all, while a <ei>ySep < 0</ei> instead allows central | |
412 | partons to scatter against either beam. | |
413 | </parm> | |
414 | ||
415 | <parm name="MultipartonInteractions:deltaYRescatter" default="1." min="0.1"> | |
416 | used for some of the <code>MultipartonInteractions:rescatterMode</code> | |
417 | options above, as the width of the rapidity transition region, where the | |
418 | probability rises from zero to unity that a scattered parton is considered | |
419 | as a potential rescatterer. | |
420 | </parm> | |
421 | ||
422 | ||
423 | <h3>Further variables</h3> | |
424 | ||
425 | These should normally not be touched. Their only function is for | |
426 | cross-checks. | |
427 | ||
428 | <modeopen name="MultipartonInteractions:nQuarkIn" default="5" min="0" | |
429 | max="5"> | |
430 | Number of allowed incoming quark flavours in the beams; a change | |
431 | to 4 would thus exclude <ei>b</ei> and <ei>bbar</ei> as incoming | |
432 | partons, etc. | |
433 | </modeopen> | |
434 | ||
435 | <modeopen name="MultipartonInteractions:nSample" default="1000" min="100"> | |
436 | The allowed <ei>pT</ei> range is split (unevenly) into 100 bins, | |
437 | and in each of these the interaction cross section is evaluated in | |
438 | <ei>nSample</ei> random phase space points. The full integral is used | |
439 | at initialization, and the differential one during the run as a | |
440 | "Sudakov form factor" for the choice of the hardest interaction. | |
441 | A larger number implies increased accuracy of the calculations. | |
442 | </modeopen> | |
443 | ||
444 | <h3>Technical notes</h3> | |
445 | ||
446 | Relative to the articles mentioned above, not much has happened. | |
447 | The main news is a technical one, that the phase space of the | |
448 | <ei>2 -> 2</ei> (massless) QCD processes is now sampled in | |
449 | <ei>dy_3 dy_4 dpT^2</ei>, where <ei>y_3</ei> and <ei>y_4</ei> are | |
450 | the rapidities of the two produced partons. One can show that | |
451 | <eq> | |
452 | (dx_1 / x_1) * (dx_2 / x_2) * d(tHat) = dy_3 * dy_4 * dpT^2 | |
453 | </eq> | |
454 | Furthermore, since cross sections are dominated by the "Rutherford" | |
455 | one of <ei>t</ei>-channel gluon exchange, which is enhanced by a | |
456 | factor of 9/4 for each incoming gluon, effective structure functions | |
457 | are defined as | |
458 | <eq> | |
459 | F(x, pT2) = (9/4) * xg(x, pT2) + sum_i xq_i(x, pT2) | |
460 | </eq> | |
461 | With this technical shift of factors 9/4 from cross sections to parton | |
462 | densities, a common upper estimate of | |
463 | <eq> | |
464 | d(sigmaHat)/d(pT2) < pi * alpha_strong^2 / pT^4 | |
465 | </eq> | |
466 | is obtained. | |
467 | ||
468 | <p/> | |
469 | In fact this estimate can be reduced by a factor of 1/2 for the | |
470 | following reason: for any configuration <ei>(y_3, y_4, pT2)</ei> also | |
471 | one with <ei>(y_4, y_3, pT2)</ei> lies in the phase space. Not both | |
472 | of those can enjoy being enhanced by the <ei>tHat -> 0</ei> | |
473 | singularity of | |
474 | <eq> | |
475 | d(sigmaHat) propto 1/tHat^2. | |
476 | </eq> | |
477 | Or if they are, which is possible with identical partons like | |
478 | <ei>q q -> q q</ei> and <ei>g g -> g g</ei>, each singularity comes | |
479 | with half the strength. So, when integrating/averaging over the two | |
480 | configurations, the estimated <ei>d(sigmaHat)/d(pT2)</ei> drops. | |
481 | Actually, it drops even further, since the naive estimate above is | |
482 | based on | |
483 | <eq> | |
484 | (4 /9) * (1 + (uHat/sHat)^2) < 8/9 < 1 | |
485 | </eq> | |
486 | The 8/9 value would be approached for <ei>tHat -> 0</ei>, which | |
487 | implies <ei>sHat >> pT2</ei> and thus a heavy parton-distribution | |
488 | penalty, while parton distributions are largest for | |
489 | <ei>tHat = uHat = -sHat/2</ei>, where the above expression | |
490 | evaluates to 5/9. A fudge factor is therefore introduced to go the | |
491 | final step, so it can easily be modified when further non-Rutherford | |
492 | processes are added, or should parton distributions change significantly. | |
493 | ||
494 | <p/> | |
495 | At initialization, it is assumed that | |
496 | <eq> | |
497 | d(sigma)/d(pT2) < d(sigmaHat)/d(pT2) * F(x_T, pT2) * F(x_T, pT2) | |
498 | * (2 y_max(pT))^2 | |
499 | </eq> | |
500 | where the first factor is the upper estimate as above, the second two | |
501 | the parton density sum evaluated at <ei>y_3 = y_ 4 = 0</ei> so that | |
502 | <ei>x_1 = x_2 = x_T = 2 pT / E_cm</ei>, where the product is expected | |
503 | to be maximal, and the final is the phase space for | |
504 | <ei>-y_max < y_{3,4} < y_max</ei>. | |
505 | The right-hand side expression is scanned logarithmically in <ei>y</ei>, | |
506 | and a <ei>N</ei> is determined such that it always is below | |
507 | <ei>N/pT^4</ei>. | |
508 | ||
509 | <p/> | |
510 | To describe the dampening of the cross section at <ei>pT -> 0</ei> by | |
511 | colour screening, the actual cross section is multiplied by a | |
512 | regularization factor <ei>(pT^2 / (pT^2 + pT0^2))^2</ei>, and the | |
513 | <ei>alpha_s</ei> is evaluated at a scale <ei>pT^2 + pT0^2</ei>, | |
514 | where <ei>pT0</ei> is a free parameter of the order of 2 - 4 GeV. | |
515 | Since <ei>pT0</ei> can be energy-dependent, an ansatz | |
516 | <eq> | |
517 | pT0(ecm) = pT0Ref * (ecm/ecmRef)^ecmPow | |
518 | </eq> | |
519 | is used, where <ei>ecm</ei> is the current CM frame energy, | |
520 | <ei>ecmRef</ei> is an arbitrary reference energy where <ei>pT0Ref</ei> | |
521 | is defined, and <ei>ecmPow</ei> gives the energy rescaling pace. For | |
522 | technical reasons, also an absolute lower <ei>pT</ei> scale <ei>pTmin</ei>, | |
523 | by default 0.2 GeV, is introduced. In principle, it is possible to | |
524 | recover older scenarios with a sharp <ei>pT</ei> cutoff by setting | |
525 | <ei>pT0 = 0</ei> and letting <ei>pTmin</ei> be a larger number. | |
526 | ||
527 | <p/> | |
528 | The above scanning strategy is then slightly modified: instead of | |
529 | an upper estimate <ei>N/pT^4</ei> one of the form | |
530 | <ei>N/(pT^2 + r * pT0^2)^2</ei> is used. At first glance, <ei>r = 1</ei> | |
531 | would seem to be fixed by the form of the regularization procedure, | |
532 | but this does not take into account the nontrivial dependence on | |
533 | <ei>alpha_s</ei>, parton distributions and phase space. A better | |
534 | Monte Carlo efficiency is obtained for <ei>r</ei> somewhat below unity, | |
535 | and currently <ei>r = 0.25</ei> is hardcoded. | |
536 | ||
537 | In the generation a trial <ei>pT2</ei> is then selected according to | |
538 | <eq> | |
539 | d(Prob)/d(pT2) = (1/sigma_ND) * N/(pT^2 + r * pT0^2)^2 * ("Sudakov") | |
540 | </eq> | |
541 | For the trial <ei>pT2</ei>, a <ei>y_3</ei> and a <ei>y_4</ei> are then | |
542 | selected, and incoming flavours according to the respective | |
543 | <ei>F(x_i, pT2)</ei>, and then the cross section is evaluated for this | |
544 | flavour combination. The ratio of trial/upper estimate gives the | |
545 | probability of survival. | |
546 | ||
547 | <p/> | |
548 | Actually, to profit from the factor 1/2 mentioned above, the cross | |
549 | section for the combination with <ei>y_3</ei> and <ei>y_4</ei> | |
550 | interchanged is also tried, which corresponds to exchanging <ei>tHat</ei> | |
551 | and <ei>uHat</ei>, and the average formed, while the final kinematics | |
552 | is given by the relative importance of the two. | |
553 | ||
554 | <p/> | |
555 | Furthermore, since large <ei>y</ei> values are disfavoured by dropping | |
556 | PDF's, a factor | |
557 | <eq> | |
558 | WT_y = (1 - (y_3/y_max)^2) * (1 - (y_4/y_max)^2) | |
559 | </eq> | |
560 | is evaluated, and used as a survival probability before the more | |
561 | time-consuming PDF+ME evaluation, with surviving events given a | |
562 | compensating weight <ei>1/WT_y</ei>. | |
563 | ||
564 | <p/> | |
565 | An impact-parameter dependence is also allowed. Based on the hard | |
566 | <ei>pT</ei> scale of the first interaction, and enhancement/depletion | |
567 | factor is picked, which multiplies the rate of subsequent interactions. | |
568 | ||
569 | <p/> | |
570 | Parton densities are rescaled and modified to take into account the | |
571 | energy-momentum and flavours kicked out by already-considered | |
572 | interactions. | |
573 | ||
574 | </chapter> | |
575 | ||
576 | <!-- Copyright (C) 2013 Torbjorn Sjostrand --> |