/usr/share/doc/pythia8-doc/html/RHadrons.html is in pythia8-doc-html 8.1.86-1.2.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 | <html>
<head>
<title>R-hadrons</title>
<link rel="stylesheet" type="text/css" href="pythia.css"/>
<link rel="shortcut icon" href="pythia32.gif"/>
</head>
<body>
<h2>R-hadrons</h2>
When a coloured SUSY particle is longer-lived than typical
hadronization scales, i.e. around c*tau > 1 fm, or equivalently
width Gamma < 0.2 GeV, it will have time to hadronize into a colour
singlet hadronic state, a R-hadron. Currently a set of such
R-hadrons have been implemented for the case of a long-lived
gluino, stop or sbottom. Needless to say, the normal case would be
that only one of them will be long-lived enough to form R-hadrons.
<p/>
For simplicity all gluino-mesons are assumed to have light-flavour
spin 1, since those are the lightest and favoured by spin-state
counting. Further, all gluino-baryons are bookkept as having
light-flavour spin 3/2, and flavours are listed in descending order.
This is more for convenience of notation, however, since the normal
baryon octet e.g. has no uuu = "p++" state. When a diquark is
extracted, a mixture of spin 0 and spin 1 is allowed. Names and codes
are essentially in agreement with the PDG conventions, e.g.
<br/>1000993 <code>R0(~g g)</code> (or gluinoball)
<br/>1009213 <code>R+(~g u dbar)</code> (or gluino-rho+)
<br/>1092214 <code>R+(~g uud)</code> (or gluino-Delta+)
<br/>For internal bookkeeping of momenta, the code 1009002,
<code>Rtemp(~g q)</code>, is used to denote the intermediate
state formed when only one of the two string pieces attached to
the gluino has broken.
<p/>
For the stop- and sbottom-hadrons the spin counting is simpler,
since it is entirely given by the constituent quark or diquark spin.
Again names and codes follow PDG conventions, e.g.
<br/>1000612 <code>R+(~t dbar)</code>
<br/>1006211 <code>R+(~t ud0)</code>
<p/>
The spin and electromagnetic charge of the new particle plays only
a minor role in the hadronization process, that can be neglected
to first approximation. Therefore it is possible to use the same
R-hadrons framework instead for other BSM scenarios with long-lived
coloured particles, e.g. with massive extra-dimensions copies
of gluons and quarks, or with leptoquarks. This can be regulated by
the switches below. Note that the codes and names of the R-hadrons
is not changed when the heavy particle involved is switched, for
reasons of administrative simplicity. R-hadron mass spectra and
other relevant particle data is automatically updated to reflect
the change, however.
<p/><code>flag </code><strong> RHadrons:allow </strong>
(<code>default = <strong>off</strong></code>)<br/>
Allows the gluino, stop and sbottom to hadronize if their respective
widths are below the limit <code>RHadrons:maxWidth</code>.
<p/><code>parm </code><strong> RHadrons:maxWidth </strong>
(<code>default = <strong>0.2</strong></code>; <code>minimum = 0.0</code>; <code>maximum = 1.0</code>)<br/>
The maximum width of the gluino for which it is possible to form
R-hadrons, provided that <code>RHadrons:allow</code> is on.
<p/><code>mode </code><strong> RHadrons:idGluino </strong>
(<code>default = <strong>1000021</strong></code>)<br/>
The gluino identity code. For other scenarios than SUSY this code
could be changed to represent another long-lived uncharged colour
octet particle, that then would be treated in the same spirit.
Could be set to 0 to forbid any gluino R-hadron formation even when
the above two criteria, <code>RHadrons:allow</code>
and <code>RHadrons:maxWidth</code>, are met.
<p/><code>mode </code><strong> RHadrons:idStop </strong>
(<code>default = <strong>1000006</strong></code>)<br/>
The lightest stop identity code. For other scenarios than SUSY this
code could be changed to represent another long-lived charge 2/3
colour triplet particle, that then would be treated in the same
spirit. As above it could be set to 0 to forbid any stop R-hadron
formation.
<p/><code>mode </code><strong> RHadrons:idSbottom </strong>
(<code>default = <strong>1000005</strong></code>)<br/>
The lightest sbottom identity code. For other scenarios than SUSY this
code could be changed to represent another long-lived charge -1/3
colour triplet particle, that then would be treated in the same
spirit. As above it could be set to 0 to forbid any sbottom R-hadron
formation.
<p/><code>flag </code><strong> RHadrons:allowDecay </strong>
(<code>default = <strong>on</strong></code>)<br/>
Allows the R-hadrons to decay or not. If the gluino/stop/sbottom is
stable or too long-lived to decay inside the detector this switch
has no real function, since then no decays will be performed anyway.
If the sparticle is so short-lived that it decays before reaching
the beam pipe then having the decay on is the logical choice.
So the interesting region is when the decays happens after the
R-hadron has passed through part of the detector, and changed its
momentum and quite possibly its flavour content before it is to
decay. Then normal decays should be switched off, and the R-hadron
tracked through matter by a program like GEANT
[<a href="Bibliography.html" target="page">Kra04,Mac07</a>]. After that, the new R-hadron info can be
overwritten into the event record and the
<code>Pythia::forceRHadronDecay()</code> method can be called
to force this modified R-hadron to decay.
<p/><code>flag </code><strong> RHadrons:setMasses </strong>
(<code>default = <strong>on</strong></code>)<br/>
Use simple mass formulae to construct all available R-hadron masses
based on the currently initialized gluino/squark masses and the
constituent masses of the other partons in the hadron. If you switch
this off, it is your responsibility to set each of the R-hadron masses
on your own, and set them in an internally consistent way. If you
mess up on this you may generate accordingly crazy results.
Specifically, it is to be assumed that none of the R-hadrons has a
mass below its constituent sparticle, i.e. that the light degrees
of freedom and the additional confinement gluon field gives a net
positive contribution to the R-hadron mass.
<p/><code>parm </code><strong> RHadrons:probGluinoball </strong>
(<code>default = <strong>0.1</strong></code>; <code>minimum = 0.0</code>; <code>maximum = 1.0</code>)<br/>
The fraction of produced gluino R-hadrons that are contain a "valence"
gluon, with the rest containing a meson or baryon quark flavour content.
<p/><code>parm </code><strong> RHadrons:mOffsetCloud </strong>
(<code>default = <strong>0.2</strong></code>; <code>minimum = 0.0</code>)<br/>
Extra mass (in GeV) added to each of the one or two extra constituent
masses in an R-hadron, to calculate the mass of a R-hadron. The same
offset is also used when the R-hadron momentum and mass is split
between the squark or gluino and the one or two light (di)quarks,
one for a squark and two for a gluino. Thus once or twice this amount
represents a part of the nominal squark or gluino mass that will not
decay weakly, since it is taken to correspond to the cloud of gluons
that surround the squark or gluino.
<p/><code>parm </code><strong> RHadrons:mCollapse </strong>
(<code>default = <strong>1.0</strong></code>; <code>minimum = 0.0</code>)<br/>
A colour singlet system with an invariant mass less than this amount,
above the R-hadron mass with the given flavour content, is assumed to
collapse to this single R-hadron, whereas a full fragmentation handling
is applied above this mass.
<p/><code>parm </code><strong> RHadrons:diquarkSpin1 </strong>
(<code>default = <strong>0.5</strong></code>; <code>minimum = 0.0</code>; <code>maximum = 1.0</code>)<br/>
Probability that a diquark extracted from the flavour code of a gluino
R-hadron should be assigned spin 1, with the rest being spin 0. Does
not apply for two identical quarks, where spin 1 is only possibility.
Note that gluino R-hadron codes for simplicity are assigned as if spin
is 1 always, and so give no guidance. For stop and sbottom the diquark
spin is preserved in the particle code, so there is no corresponding
issue.
</body>
</html>
<!-- Copyright (C) 2014 Torbjorn Sjostrand -->
|