Context Navigation

← Previous Revision
Latest Revision
Next Revision →
Blame
Revision Log

ParticleDecays.php @ 1

Last change on this file since 1 was 1, checked in by zerwas, 11 years ago
first import of structure, PYTHIA8 and DELPHES
File size: 23.0 KB

Line
1	<html>
2	<head>
3	<title>Particle Decays</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='ParticleDecays.php'>
29
30	<h2>Particle Decays</h2>
31
32	The <code>ParticleDecays</code> class performs the sequential decays of
33	all unstable hadrons produced in the string fragmentation stage,
34	i.e. up to and including <i>b</i> hadrons and their decay products,
35	such as the <i>tau</i> lepton. It is not to be used for the decay of
36	more massive <?php $filepath = $_GET["filepath"];
37	echo "<a href='ResonanceDecays.php?filepath=".$filepath."' target='page'>";?>resonances</a>, such as top,
38	<i>Z^0</i> or SUSY, where decays must be performed already at the
39	<code>ProcessLevel</code> of the event generation.
40
41	<p/>
42	The decay description essentially copies the one present in
43	PYTHIA since many years, but with some improvements, e.g. in the decay
44	tables and the number of decay models available. Recently a more
45	sophisticated handling of <i>tau</i> decays has also been introduced.
46	Some issues may need further polishing.
47
48	<h3>Variables determining whether a particle decays</h3>
49
50	Before a particle is actually decayed, a number of checks are made.
51
52	<p/>
53	(i) Decay modes must have been defined for the particle kind;
54	tested by the <code>canDecay()</code> method of <code>Event</code>
55	(and <code>ParticleData</code>).
56
57	<p/>
58	(ii) The main switch for allowing this particle kind to decay must
59	be on; tested by the <code>mayDecay()</code> method of <code>Event</code>
60	(and <code>ParticleData</code>).
61
62	<p/>
63	(iii) Particles may be requested to have a nominal proper lifetime
64	<i>tau0</i> below a threshold.
65
66	<br/><br/><strong>ParticleDecays:limitTau0</strong> <input type="radio" name="1" value="on"><strong>On</strong>
67	<input type="radio" name="1" value="off" checked="checked"><strong>Off</strong>
68	(<code>default = <strong>off</strong></code>)<br/>
69	When on, only particles with <i>tau0 < tau0Max</i> are decayed.
70
71
72	<br/><br/><table><tr><td><strong>ParticleDecays:tau0Max </td><td></td><td> <input type="text" name="2" value="10." size="20"/>   (<code>default = <strong>10.</strong></code>; <code>minimum = 0.</code>)</td></tr></table>
73	The above <i>tau0Max</i>, expressed in mm/c.
74
75
76	<p/>
77	(iv) Particles may be requested to have an actual proper lifetime
78	<i>tau</i> below a threshold.
79
80	<br/><br/><strong>ParticleDecays:limitTau</strong> <input type="radio" name="3" value="on"><strong>On</strong>
81	<input type="radio" name="3" value="off" checked="checked"><strong>Off</strong>
82	(<code>default = <strong>off</strong></code>)<br/>
83	When on, only particles with <i>tau < tauMax</i> are decayed.
84
85
86	<br/><br/><table><tr><td><strong>ParticleDecays:tauMax </td><td></td><td> <input type="text" name="4" value="10." size="20"/>   (<code>default = <strong>10.</strong></code>; <code>minimum = 0.</code>)</td></tr></table>
87	The above <i>tauMax</i>, expressed in mm/c.<br/>
88	In order for this and the subsequent tests to work, a <i>tau</i>
89	is selected and stored for each particle, whether in the end it
90	decays or not. (If each test would use a different temporary
91	<i>tau</i> it would lead to inconsistencies.)
92
93
94	<p/>
95	(v) Particles may be requested to decay within a given distance
96	of the origin.
97
98	<br/><br/><strong>ParticleDecays:limitRadius</strong> <input type="radio" name="5" value="on"><strong>On</strong>
99	<input type="radio" name="5" value="off" checked="checked"><strong>Off</strong>
100	(<code>default = <strong>off</strong></code>)<br/>
101	When on, only particles with a decay within a radius <i>r < rMax</i>
102	are decayed. There is assumed to be no magnetic field or other
103	detector effects.
104
105
106	<br/><br/><table><tr><td><strong>ParticleDecays:rMax </td><td></td><td> <input type="text" name="6" value="10." size="20"/>   (<code>default = <strong>10.</strong></code>; <code>minimum = 0.</code>)</td></tr></table>
107	The above <i>rMax</i>, expressed in mm.
108
109
110	<p/>
111	(vi) Particles may be requested to decay within a given cylidrical
112	volume around the origin.
113
114	<br/><br/><strong>ParticleDecays:limitCylinder</strong> <input type="radio" name="7" value="on"><strong>On</strong>
115	<input type="radio" name="7" value="off" checked="checked"><strong>Off</strong>
116	(<code>default = <strong>off</strong></code>)<br/>
117	When on, only particles with a decay within a volume limited by
118	<i>rho = sqrt(x^2 + y^2) < xyMax</i> and <i>\|z\| < zMax</i>
119	are decayed. There is assumed to be no magnetic field or other
120	detector effects.
121
122
123	<br/><br/><table><tr><td><strong>ParticleDecays:xyMax </td><td></td><td> <input type="text" name="8" value="10." size="20"/>   (<code>default = <strong>10.</strong></code>; <code>minimum = 0.</code>)</td></tr></table>
124	The above <i>xyMax</i>, expressed in mm.
125
126
127	<br/><br/><table><tr><td><strong>ParticleDecays:zMax </td><td></td><td> <input type="text" name="9" value="10." size="20"/>   (<code>default = <strong>10.</strong></code>; <code>minimum = 0.</code>)</td></tr></table>
128	The above <i>zMax</i>, expressed in mm.
129
130
131	<h3>Mixing</h3>
132
133	<br/><br/><strong>ParticleDecays:mixB</strong> <input type="radio" name="10" value="on" checked="checked"><strong>On</strong>
134	<input type="radio" name="10" value="off"><strong>Off</strong>
135	(<code>default = <strong>on</strong></code>)<br/>
136	Allow or not <i>B^0 - B^0bar</i> and <i>B_s^0 - B_s^0bar</i> mixing.
137
138
139	<br/><br/><table><tr><td><strong>ParticleDecays:xBdMix </td><td></td><td> <input type="text" name="11" value="0.776" size="20"/>   (<code>default = <strong>0.776</strong></code>; <code>minimum = 0.74</code>; <code>maximum = 0.81</code>)</td></tr></table>
140	The mixing parameter <i>x_d = Delta(m_B^0)/Gamma_B^0</i> in the
141	<i>B^0 - B^0bar</i> system. (Default from RPP2006.)
142
143
144	<br/><br/><table><tr><td><strong>ParticleDecays:xBsMix </td><td></td><td> <input type="text" name="12" value="26.05" size="20"/>   (<code>default = <strong>26.05</strong></code>; <code>minimum = 22.0</code>; <code>maximum = 30.0</code>)</td></tr></table>
145	The mixing parameter <i>x_s = Delta(m_B_s^0)/Gamma_B_s^0</i> in the
146	<i>B_s^0 - B_s^0bar</i> system. (Delta-m from CDF hep-ex-0609040,
147	Gamma from RPP2006.)
148
149
150	<h3>Tau decays</h3>
151
152	A new machinery has been introduced to handle <i>tau</i> lepton decays,
153	with helicity information related to the production process and with
154	the form of the hadronic current fitted to data. It is largely based
155	on the corresponding Herwig++ implementation [<a href="Bibliography.php" target="page">Gre07</a>], with
156	some input from Tauola [<a href="Bibliography.php" target="page">Jad90</a>]. A complete writeup is
157	in preparation [<a href="Bibliography.php" target="page">Ilt11</a>].
158
159	<p/>
160	For <i>tau</i>s in external processes, interfaced with Les Houches
161	Acccord information available, e.g. via Les Houches Event Files (LHEF),
162	the new machinery interprets the SPINUP number for <i>tau</i> leptons
163	as giving their helicity, and decays them accordingly. The only exceptions
164	are when a specific polarization is forced by the user (see below),
165	which then overrides the SPINUP value, or when SPINUP has the special
166	value 9 (unpolarized). In the latter case, PYTHIA defaults back to
167	attempting to determine the helicity structure from the production
168	process, in the same way as for internal processes.
169
170	<p/>
171	This new machinery is on by default, but it is possible to revert to
172	the simpler old decay handling, e.g. to study differences. Furthermore
173	the spin tracing framework does not yet cover all possibilities; notably
174	it cannot handle taus coming from SUSY decay chains
175	(except via LHEF), so it makes sense
176	to switch off the new machinery in such instances, for speed reasons if
177	nothing else. In case only one tau mother species is undefined, the
178	polarization involved can be set by hand.
179
180	<br/><br/><table><tr><td><strong>ParticleDecays:sophisticatedTau </td><td>   (<code>default = <strong>1</strong></code>; <code>minimum = 0</code>; <code>maximum = 3</code>)</td></tr></table>
181	Choice of <ei>tau</ei> decay model.
182	<br/>
183	<input type="radio" name="13" value="0"><strong>0 </strong>: old decay model, with isotropic decays. When reading LHEF files, the SPINUP digit will be ignored.<br/>
184	<input type="radio" name="13" value="1" checked="checked"><strong>1 </strong>: sophisticated decays where <ei>tau</ei> polarization is calculated from the <ei>tau</ei> production mechanism. When reading LHEF files, the SPINUP digit will be used. <br/>
185	<input type="radio" name="13" value="2"><strong>2 </strong>: sophisticated decays as above, but additionally <ei>tau</ei> polarization is set to <code>ParticleDecaus:tauPolarization</code> for <ei>tau</ei>s produced from <code>ParticleDecays:tauMother</code>. When reading LHEF files, this overrides the SPINUP digit. <br/>
186	<input type="radio" name="13" value="3"><strong>3 </strong>: sophisticated decays where <ei>tau</ei> polarization is set to <code>ParticleDecaus:tauPolarization</code> for all <ei>tau</ei> decays. When reading LHEF files, this overrides the SPINUP digit. <br/>
187	<br/><b>Note</b>: options <code>2</code> and <code>3</code>,
188	to force a specific <ei>tau</ei> polarization, only affect the decay
189	of the <ei>tau</ei>. The angular distribution of the <ei>tau</ei> itself,
190	given by its production, is not modified by these options. If you want, e.g.,
191	a righthanded <ei>W</ei>, or a SUSY decay chain, the kinematics should
192	be handled by the corresponding cross section class(es), supplemented by
193	the resonance decay one(s). The options here could then still be used
194	to ensure the correct polarization at the <ei>tau</ei> decay stage.
195
196	<br/><br/><table><tr><td><strong>ParticleDecays:tauPolarization </td><td></td><td> <input type="text" name="14" value="0" size="20"/>   (<code>default = <strong>0</strong></code>; <code>minimum = -1.</code>; <code>maximum = 1.</code>)</td></tr></table>
197	Polarization of the <i>tau</i> when mode <i>2</i> or <i>3</i> of
198	<code>ParticleDecays:sophisticatedTau</code> is selected.
199
200
201	<br/><br/><table><tr><td><strong>ParticleDecays:tauMother </td><td></td><td> <input type="text" name="15" value="0" size="20"/>   (<code>default = <strong>0</strong></code>; <code>minimum = 0</code>)</td></tr></table>
202	Mother of the <i>tau</i> for forced polarization when mode <i>2</i> of
203	<code>ParticleDecays:sophisticatedTau</code> is selected. You should give the
204	positive identity code; to the extent an antiparticle exists it will
205	automatically obtain the inverse polarization.
206
207
208	<h3>Other variables</h3>
209
210	<br/><br/><table><tr><td><strong>ParticleDecays:mSafety </td><td></td><td> <input type="text" name="16" value="0.0005" size="20"/>   (<code>default = <strong>0.0005</strong></code>; <code>minimum = 0.</code>; <code>maximum = 0.01</code>)</td></tr></table>
211	Minimum mass difference required between the decaying mother mass
212	and the sum of the daughter masses, kept as a safety margin to avoid
213	numerical problems in the decay generation.
214
215
216	<br/><br/><table><tr><td><strong>ParticleDecays:sigmaSoft </td><td></td><td> <input type="text" name="17" value="0.5" size="20"/>   (<code>default = <strong>0.5</strong></code>; <code>minimum = 0.2</code>; <code>maximum = 2.</code>)</td></tr></table>
217	In semileptonic decays to more than one hadron, such as
218	<i>B -> nu l D pi</i>, decay products after the first three are
219	dampened in momentum by an explicit weight factor
220	<i>exp(-p^2/sigmaSoft^2)</i>, where <i>p</i> is the
221	three-momentum in the rest frame of the decaying particle.
222	This takes into account that such further particles come from the
223	fragmentation of the spectator parton and thus should be soft.
224
225
226	<p/>
227	When a decay mode is defined in terms of a partonic content, a random
228	multiplicity (and a random flavour set) of hadrons is to be picked,
229	especially for some charm and bottom decays. This is done according to
230	a Poissonian distribution, for <i>n_p</i> normal particles and
231	<i>n_q</i> quarks the average value is chosen as
232	<br/><i>
233	n_p/ 2 + n_q/4 + multIncrease * ln ( mDiff / multRefMass)
234	</i><br/>
235	with <i>mDiff</i> the difference between the decaying particle mass
236	and the sum of the normal-particle masses and the constituent quark masses.
237	For gluonic systems <i>multGoffset</i> offers and optional additional
238	term to the multiplicity. The lowest possible multiplicity is
239	<i>n_p + n_q/2</i> (but at least 2) and the highest possible 10.
240	If the picked hadrons have a summed mass above that of the mother a
241	new try is made, including a new multiplicity. These constraints
242	imply that the actual average multiplicity does not quite agree with
243	the formula above.
244
245	<br/><br/><table><tr><td><strong>ParticleDecays:multIncrease </td><td></td><td> <input type="text" name="18" value="4." size="20"/>   (<code>default = <strong>4.</strong></code>; <code>minimum = 2.</code>; <code>maximum = 6.</code>)</td></tr></table>
246	The above <i>multIncrease</i> parameter, except for
247	<code>meMode = 23</code>.
248
249
250	<br/><br/><table><tr><td><strong>ParticleDecays:multIncreaseWeak </td><td></td><td> <input type="text" name="19" value="2.5" size="20"/>   (<code>default = <strong>2.5</strong></code>; <code>minimum = 1.</code>; <code>maximum = 4.</code>)</td></tr></table>
251	The above <i>multIncrease</i> parameter, specifically for
252	<code>meMode = 23</code>. Here the weak decay implies that only the
253	virtual W mass should contribute to the production of new particles,
254	rather than the full meson mass.
255
256
257	<br/><br/><table><tr><td><strong>ParticleDecays:multRefMass </td><td></td><td> <input type="text" name="20" value="0.7" size="20"/>   (<code>default = <strong>0.7</strong></code>; <code>minimum = 0.2</code>; <code>maximum = 2.0</code>)</td></tr></table>
258	The above <i>multRefMass</i> parameter.
259
260
261	<br/><br/><table><tr><td><strong>ParticleDecays:multGoffset </td><td></td><td> <input type="text" name="21" value="0.5" size="20"/>   (<code>default = <strong>0.5</strong></code>; <code>minimum = 0.0</code>; <code>maximum = 2.0</code>)</td></tr></table>
262	The above <i>multGoffset</i> parameter.
263
264
265	<br/><br/><table><tr><td><strong>ParticleDecays:colRearrange </td><td></td><td> <input type="text" name="22" value="0.5" size="20"/>   (<code>default = <strong>0.5</strong></code>; <code>minimum = 0.</code>; <code>maximum = 1.0</code>)</td></tr></table>
266	When a decay is given as a list of four partons to be turned into
267	hadrons (primarily for modes 41 - 80) it is assumed that they are
268	listed in pairs, as a first and a second colour singlet, which could
269	give rise to separate sets of hadrons. Here <i>colRearrange</i> is
270	the probability that this original assignment is not respected, and
271	default corresponds to no memory of this original colour topology.
272
273
274	<br/><br/><strong>ParticleDecays:FSRinDecays</strong> <input type="radio" name="23" value="on"><strong>On</strong>
275	<input type="radio" name="23" value="off"><strong>Off</strong>
276	(<code>default = <strong>true</strong></code>)<br/>
277	When a particle decays to <i>q qbar</i>, <i>g g</i>, <i>g g g</i>
278	or <i>gamma g g</i>, with <code>meMode > 90</code>, allow or not a
279	shower to develop from it, before the partonic system is hadronized.
280	(The typical example is <i>Upsilon</i> decay.)
281
282
283	In addition, some variables defined for string fragmentation and for
284	flavour production are used also here.
285
286	<h3>Modes for Matrix Element Processing</h3>
287
288	Some decays can be treated better than what pure phase space allows,
289	by reweighting with appropriate matrix elements. In others a partonic
290	content has to be converted to a set of hadrons. The presence of such
291	corrections is signalled by a nonvanishing <code>meMode()</code> value
292	for a decay mode in the <?php $filepath = $_GET["filepath"];
293	echo "<a href='ParticleDataScheme.php?filepath=".$filepath."' target='page'>";?>particle
294	data table</a>. The list of allowed possibilities almost agrees with the
295	PYTHIA 6 ones, but several obsolete choices have been removed,
296	a few new introduced, and most have been moved for better consistency.
297	Here is the list of currently allowed <code>meMode()</code> codes:
298	<ul>
299	<li> 0 : pure phase space of produced particles ("default");
300	input of partons is allowed and then the partonic content is
301	converted into the minimal number of hadrons (i.e. one per
302	parton pair, but at least two particles in total)</li>
303	<li> 1 : <i>omega</i> and <i>phi -> pi+ pi- pi0</i></li>
304	<li> 2 : polarization in <i>V -> PS + PS</i> (<i>V</i> = vector,
305	<i>PS</i> = pseudoscalar), when <i>V</i> is produced by
306	<i>PS -> PS + V</i> or <i>PS -> gamma + V</i></li>
307	<li> 11 : Dalitz decay into one particle, in addition to the
308	lepton pair (also allowed to specify a quark-antiquark pair that
309	should collapse to a single hadron)</li>
310	<li> 12 : Dalitz decay into two or more particles in addition
311	to the lepton pair</li>
312	<li> 13 : double Dalitz decay into two lepton pairs</li>
313	<li> 21 : decay to phase space, but weight up <i>neutrino_tau</i> spectrum
314	in <i>tau</i> decay</li>
315	<li> 22 : weak decay; if there is a quark spectator system it collapses to
316	one hadron; for leptonic/semileptonic decays the <i>V-A</i> matrix element
317	is used, for hadronic decays simple phase space</li>
318	<li> 23 : as 22, but require at least three particles in decay</li>
319	<li> 31 : decays of type B -> gamma X, very primitive simulation where
320	X is given in terms of its flavour content, the X multiplicity is picked
321	according to a geometrical distribution with average number 2, and
322	the photon energy spectrum is weighted up relative to pure phase space</li>
323	<li> 42 - 50 : turn partons into a random number of hadrons, picked according
324	to a Poissonian with average value as described above, but at least
325	<code>code</code> - 40 and at most 10, and then distribute then in pure
326	phase space; make a new try with another multiplicity if the sum of daughter
327	masses exceed the mother one </li>
328	<li> 52 - 60 : as 42 - 50, with multiplicity between <code>code</code> - 50
329	and 10, but avoid already explicitly listed non-partonic channels</li>
330	<li> 62 - 70 : as 42 - 50, but fixed multiplicity <code>code</code> - 60</li>
331	<li> 72 - 80 : as 42 - 50, but fixed multiplicity <code>code</code> - 70,
332	and avoid already explicitly listed non-partonic channels</li>
333	<li> 91 : decay to <i>q qbar</i> or <i>g g</i>, which should shower
334	and hadronize</li>
335	<li> 92 : decay onium to <i>g g g</i> or <i>g g gamma</i>
336	(with matrix element), which should shower and hadronize</li>
337	<li> 100 - : reserved for the description of partial widths of
338	<?php $filepath = $_GET["filepath"];
339	echo "<a href='ResonanceDecays.php?filepath=".$filepath."' target='page'>";?>resonances</a></li>
340	</ul>
341
342	Three special decay product identity codes are defined.
343	<ul>
344	<li>81: remnant flavour. Used for weak decays of c and b hadrons, where the
345	c or b quark decays and the other quarks are considered as a spectator
346	remnant in this decay. In practice only used for baryons with multiple
347	c and b quarks, which presumably would never be used, but have simple
348	(copied) just-in-case decay tables. Assumed to be last decay product.</li>
349	<li>82: random flavour, picked by the standard fragmentation flavour
350	machinery, used to start a sequence of hadrons, for matrix element
351	codes in 41 - 80. Assumed to be first decay product, with -82 as second
352	and last. Where multiplicity is free to be picked it is selected as for
353	normal quarkonic systems. Currently unused.</li>
354	<li>83: as for 82, with matched pair 83, -83 of decay products. The
355	difference is that here the pair is supposed to come from a closed gluon
356	loop (e.g. <i>eta_c -> g g</i>) and so have a somewhat higher average
357	multiplicity than the simple string assumed for 82, see the
358	<code>ParticleDecays:multGoffset</code> parameter above.</li>
359	</ul>
360
361	<input type="hidden" name="saved" value="1"/>
362
363	<?php
364	echo "<input type='hidden' name='filepath' value='".$_GET["filepath"]."'/>"?>
365
366	<table width="100%"><tr><td align="right"><input type="submit" value="Save Settings" /></td></tr></table>
367	</form>
368
369	<?php
370
371	if($_POST["saved"] == 1)
372	{
373	$filepath = $_POST["filepath"];
374	$handle = fopen($filepath, 'a');
375
376	if($_POST["1"] != "off")
377	{
378	$data = "ParticleDecays:limitTau0 = ".$_POST["1"]."\n";
379	fwrite($handle,$data);
380	}
381	if($_POST["2"] != "10.")
382	{
383	$data = "ParticleDecays:tau0Max = ".$_POST["2"]."\n";
384	fwrite($handle,$data);
385	}
386	if($_POST["3"] != "off")
387	{
388	$data = "ParticleDecays:limitTau = ".$_POST["3"]."\n";
389	fwrite($handle,$data);
390	}
391	if($_POST["4"] != "10.")
392	{
393	$data = "ParticleDecays:tauMax = ".$_POST["4"]."\n";
394	fwrite($handle,$data);
395	}
396	if($_POST["5"] != "off")
397	{
398	$data = "ParticleDecays:limitRadius = ".$_POST["5"]."\n";
399	fwrite($handle,$data);
400	}
401	if($_POST["6"] != "10.")
402	{
403	$data = "ParticleDecays:rMax = ".$_POST["6"]."\n";
404	fwrite($handle,$data);
405	}
406	if($_POST["7"] != "off")
407	{
408	$data = "ParticleDecays:limitCylinder = ".$_POST["7"]."\n";
409	fwrite($handle,$data);
410	}
411	if($_POST["8"] != "10.")
412	{
413	$data = "ParticleDecays:xyMax = ".$_POST["8"]."\n";
414	fwrite($handle,$data);
415	}
416	if($_POST["9"] != "10.")
417	{
418	$data = "ParticleDecays:zMax = ".$_POST["9"]."\n";
419	fwrite($handle,$data);
420	}
421	if($_POST["10"] != "on")
422	{
423	$data = "ParticleDecays:mixB = ".$_POST["10"]."\n";
424	fwrite($handle,$data);
425	}
426	if($_POST["11"] != "0.776")
427	{
428	$data = "ParticleDecays:xBdMix = ".$_POST["11"]."\n";
429	fwrite($handle,$data);
430	}
431	if($_POST["12"] != "26.05")
432	{
433	$data = "ParticleDecays:xBsMix = ".$_POST["12"]."\n";
434	fwrite($handle,$data);
435	}
436	if($_POST["13"] != "1")
437	{
438	$data = "ParticleDecays:sophisticatedTau = ".$_POST["13"]."\n";
439	fwrite($handle,$data);
440	}
441	if($_POST["14"] != "0")
442	{
443	$data = "ParticleDecays:tauPolarization = ".$_POST["14"]."\n";
444	fwrite($handle,$data);
445	}
446	if($_POST["15"] != "0")
447	{
448	$data = "ParticleDecays:tauMother = ".$_POST["15"]."\n";
449	fwrite($handle,$data);
450	}
451	if($_POST["16"] != "0.0005")
452	{
453	$data = "ParticleDecays:mSafety = ".$_POST["16"]."\n";
454	fwrite($handle,$data);
455	}
456	if($_POST["17"] != "0.5")
457	{
458	$data = "ParticleDecays:sigmaSoft = ".$_POST["17"]."\n";
459	fwrite($handle,$data);
460	}
461	if($_POST["18"] != "4.")
462	{
463	$data = "ParticleDecays:multIncrease = ".$_POST["18"]."\n";
464	fwrite($handle,$data);
465	}
466	if($_POST["19"] != "2.5")
467	{
468	$data = "ParticleDecays:multIncreaseWeak = ".$_POST["19"]."\n";
469	fwrite($handle,$data);
470	}
471	if($_POST["20"] != "0.7")
472	{
473	$data = "ParticleDecays:multRefMass = ".$_POST["20"]."\n";
474	fwrite($handle,$data);
475	}
476	if($_POST["21"] != "0.5")
477	{
478	$data = "ParticleDecays:multGoffset = ".$_POST["21"]."\n";
479	fwrite($handle,$data);
480	}
481	if($_POST["22"] != "0.5")
482	{
483	$data = "ParticleDecays:colRearrange = ".$_POST["22"]."\n";
484	fwrite($handle,$data);
485	}
486	if($_POST["23"] != "true")
487	{
488	$data = "ParticleDecays:FSRinDecays = ".$_POST["23"]."\n";
489	fwrite($handle,$data);
490	}
491	fclose($handle);
492	}
493
494	?>
495	</body>
496	</html>
497
498	<!-- Copyright (C) 2012 Torbjorn Sjostrand -->
499

Note: See TracBrowser for help on using the repository browser.

Context Navigation

source: HiSusy/trunk/Pythia8/pythia8170/phpdoc/ParticleDecays.php @ 1

Download in other formats: