]> git.uio.no Git - u/mrichter/AliRoot.git/blob - PYTHIA8/pythia8170/phpdoc/BeamParameters.php
Protection for dereferencing fTDCErrorBuffer. see ALIROOT-5749
[u/mrichter/AliRoot.git] / PYTHIA8 / pythia8170 / phpdoc / BeamParameters.php
1 <html>
2 <head>
3 <title>Beam Parameters</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='BeamParameters.php'>
29
30 <h2>Beam Parameters</h2>
31
32 The settings on this page relate to the beam identities and energies, 
33 to a beam momentum spread and to a beam interaction spot. 
34 As always, momenta and energies are to be given in units of GeV,
35 and of space and time in mm. 
36
37 <h3>Incoming beams</h3>
38
39 There are two ways to set the identities and energies of the two incoming 
40 beam particles. One is to use the <code>init()</code> method with no 
41 arguments. Then the settings variables below will be read and used. The 
42 alternative is to call <code><?php $filepath = $_GET["filepath"];
43 echo "<a href='ProgramFlow.php?filepath=".$filepath."' target='page'>";?>init(...)</a></code> 
44 with arguments that provide this information. Then you need not use the 
45 variables below (although it would still be possible). Note that, if nothing 
46 is done, you will default to LHC at 14 TeV.
47
48 <p/>
49 Currently the beam particles must be either a hadron pair or a lepton
50 pair. In the former category <i>p p</i> and <i>pbar p</i> 
51 combinations dominate, but it is also possible to combine with 
52 <i>pi^+</i>, <i>pi^-</i> and <i>pi^0</i>. In the latter 
53 <i>e^+ e^-</i> and <i>mu^+ mu^-</i> would be the most useful 
54 combinations, but also others should work if combined with an 
55 appropriate hard process.
56
57 <br/><br/><table><tr><td><strong>Beams:idA  </td><td></td><td> <input type="text" name="1" value="2212" size="20"/>  &nbsp;&nbsp;(<code>default = <strong>2212</strong></code>)</td></tr></table>
58 The PDG <code>id</code> code for the first incoming particle.
59 Allowed codes include 
60 <br/><i>2212 = p</i>, <i>-2212 = pbar</i>, 
61 <br/><i>211 = pi^+</i>, <i>-211 = pi^-</i>, <i>111 = pi^0</i>,  
62 <br/><i>990 = Pomeron</i> (used in diffractive machinery;
63 here mainly for debug purposes),
64 <br/><i>11 = e^-</i>, <i>-11 = e^+</i>, 
65 <br/><i>13 = mu^-</i>, <i>-13 = mu^+</i>, 
66 <br/>and a few more leptons/neutrinos in a few combinations.
67   
68
69 <br/><br/><table><tr><td><strong>Beams:idB  </td><td></td><td> <input type="text" name="2" value="2212" size="20"/>  &nbsp;&nbsp;(<code>default = <strong>2212</strong></code>)</td></tr></table>
70 The PDG <code>id</code> code for the second incoming particle.
71   
72  
73 <br/><br/><table><tr><td><strong>Beams:frameType  </td><td>  &nbsp;&nbsp;(<code>default = <strong>1</strong></code>; <code>minimum = 1</code>; <code>maximum = 5</code>)</td></tr></table>
74 Choice of frame for the two colliding particles. For options
75 1 - 3 the beam identities are specified above, while they are
76 obtained by the Les Houches information for options 4 and 5.
77 <br/>
78 <input type="radio" name="3" value="1" checked="checked"><strong>1 </strong>: the beams are colliding in their CM frame,  and therefore only the CM energy needs to be provided, see  <code>Beams:eCM</code> below. <br/>
79 <input type="radio" name="3" value="2"><strong>2 </strong>: the beams are back-to-back, but with different energies, see <code>Beams:eA</code> and <code>Beams:eB</code> below. This option could also be used for fixed-target configurations. <br/>
80 <input type="radio" name="3" value="3"><strong>3 </strong>: the beams are not back-to-back, and therefore the three-momentum of each incoming particle needs to be specified, see <code>Beams:pxA</code> through <code>Beams:pzB</code> below. <br/>
81 <input type="radio" name="3" value="4"><strong>4 </strong>: the beam and event information is stored in a  <aloc href="LesHouchesAccord">Les Houches Event File</aloc>,  see <code>Beams:LHEF</code> below. <br/>
82 <input type="radio" name="3" value="5"><strong>5 </strong>: the beam and event information is obtained by a pointer to an <code><aloc href="LesHouchesAccord">LHAup</aloc></code>  class instance.   <br/>
83
84 <br/><br/><table><tr><td><strong>Beams:eCM </td><td></td><td> <input type="text" name="4" value="14000." size="20"/>  &nbsp;&nbsp;(<code>default = <strong>14000.</strong></code>; <code>minimum = 10.</code>)</td></tr></table>
85 Collision CM energy, to be set if <code>Beams:frameType</code> = 1. 
86   
87
88 <br/><br/><table><tr><td><strong>Beams:eA </td><td></td><td> <input type="text" name="5" value="7000." size="20"/>  &nbsp;&nbsp;(<code>default = <strong>7000.</strong></code>; <code>minimum = 0.</code>)</td></tr></table>
89 The energy of the first incoming particle, moving in the 
90 <i>+z </i>direction, to be set if <code>Beams:frameType</code> = 2. 
91 If the particle energy is smaller than its mass
92 it is assumed to be at rest. 
93   
94
95 <br/><br/><table><tr><td><strong>Beams:eB </td><td></td><td> <input type="text" name="6" value="7000." size="20"/>  &nbsp;&nbsp;(<code>default = <strong>7000.</strong></code>; <code>minimum = 0.</code>)</td></tr></table>
96 The energy of the second incoming particle, moving in the 
97 <i>-z</i> direction, to be set if <code>Beams:frameType</code> = 2. 
98 If the particle energy is smaller than its mass
99 it is assumed to be at rest.
100   
101
102 <br/><br/><table><tr><td><strong>Beams:pxA </td><td></td><td> <input type="text" name="7" value="0." size="20"/>  &nbsp;&nbsp;(<code>default = <strong>0.</strong></code>)</td></tr></table>
103 The <i>p_x</i> component of the first incoming particle, 
104 to be set if <code>Beams:frameType</code> = 3. 
105   
106
107 <br/><br/><table><tr><td><strong>Beams:pyA </td><td></td><td> <input type="text" name="8" value="0." size="20"/>  &nbsp;&nbsp;(<code>default = <strong>0.</strong></code>)</td></tr></table>
108 The <i>p_y</i> component of the first incoming particle, 
109 to be set if <code>Beams:frameType</code> = 3. 
110   
111
112 <br/><br/><table><tr><td><strong>Beams:pzA </td><td></td><td> <input type="text" name="9" value="7000." size="20"/>  &nbsp;&nbsp;(<code>default = <strong>7000.</strong></code>)</td></tr></table>
113 The <i>p_z</i> component of the first incoming particle, 
114 to be set if <code>Beams:frameType</code> = 3. 
115   
116
117 <br/><br/><table><tr><td><strong>Beams:pxB </td><td></td><td> <input type="text" name="10" value="0." size="20"/>  &nbsp;&nbsp;(<code>default = <strong>0.</strong></code>)</td></tr></table>
118 The <i>p_x</i> component of the second incoming particle, 
119 to be set if <code>Beams:frameType</code> = 3. 
120   
121
122 <br/><br/><table><tr><td><strong>Beams:pyB </td><td></td><td> <input type="text" name="11" value="0." size="20"/>  &nbsp;&nbsp;(<code>default = <strong>0.</strong></code>)</td></tr></table>
123 The <i>p_y</i> component of the second incoming particle, 
124 to be set if <code>Beams:frameType</code> = 3. 
125   
126
127 <br/><br/><table><tr><td><strong>Beams:pzB </td><td></td><td> <input type="text" name="12" value="-7000." size="20"/>  &nbsp;&nbsp;(<code>default = <strong>-7000.</strong></code>)</td></tr></table>
128 The <i>p_z</i> component of the second incoming particle, 
129 to be set if <code>Beams:frameType</code> = 3. 
130   
131
132 <br/><br/><table><tr><td><strong>Beams:LHEF  </td><td></td><td> <input type="text" name="13" value="void" size="20"/>  &nbsp;&nbsp;(<code>default = <strong>void</strong></code>)</td></tr></table>
133 The name of a Les Houches Event File,
134 to be set if <code>Beams:frameType</code> = 4. 
135   
136
137 <br/><br/><table><tr><td><strong>Beams:LHEFheader  </td><td></td><td> <input type="text" name="14" value="void" size="20"/>  &nbsp;&nbsp;(<code>default = <strong>void</strong></code>)</td></tr></table>
138 As some information in a Les Houches Event File init block is only known
139 at the end of generation, some programs choose to output this as a
140 separate file. If <code>Beams:LHEFheader</code> is given, information up
141 till the end of the init block is read from this file, with
142 the events themselves read as usual from the file given in
143 <code>Beams:LHEF</code>.
144   
145
146 <br/><br/><strong>Beams:newLHEFsameInit</strong>  <input type="radio" name="15" value="on"><strong>On</strong>
147 <input type="radio" name="15" value="off" checked="checked"><strong>Off</strong>
148  &nbsp;&nbsp;(<code>default = <strong>off</strong></code>)<br/>
149 Allow to begin reading events from a new LHEF or or a new 
150 <code>LHAup</code> instance without a completely new initialization. 
151 Only useful when <code>Beams:frameType</code> = 4 or 5.
152   
153
154 <br/><br/><strong>Beams:readLHEFheaders</strong>  <input type="radio" name="16" value="on" checked="checked"><strong>On</strong>
155 <input type="radio" name="16" value="off"><strong>Off</strong>
156  &nbsp;&nbsp;(<code>default = <strong>on</strong></code>)<br/>
157 Read in LHEF header blocks and store them in the
158 <?php $filepath = $_GET["filepath"];
159 echo "<a href='EventInformation.php?filepath=".$filepath."' target='page'>";?>Info</a> class. See also
160 <?php $filepath = $_GET["filepath"];
161 echo "<a href='LesHouchesAccord.php?filepath=".$filepath."' target='page'>";?>LHAupLHEF</a> for more information.
162   
163
164 <p/><code>mode&nbsp; </code><strong> Beams:nSkipLHEFatInit &nbsp;</strong> 
165  (<code>default = <strong>0</strong></code>)<br/>
166 Skip the first <i>nSkip</i> events of the input stream 
167 (cf. the <code>LHAup::skipEvent(nSkip)</code> method).
168 Only used when <code>Beams:frameType</code> = 4 or 5.
169   
170
171 <h3>Beam momentum spread</h3>
172
173 This framework currently is intended for a modest beam spread, such as
174 experienced at hadron colliders. Thus it can be safely assumed that the 
175 physics does not change over the CM energy range probed, so that the 
176 parameters of the physics initialization at the nominal energy can be
177 used as is. Currently it can <b>not</b> be used for the more extensive
178 energy spread expected at linear <i>e^+ e^-</i> colliders. Also,
179 any attempt to combine it with external Les Houches input of 
180 parton-level events is at own risk. 
181      
182 <p/>
183 On this page you can set the momentum spread according to a simple 
184 Gaussian distribution. If you instead want a more sophisticated 
185 parametrization, you can write and link your own 
186 <code><?php $filepath = $_GET["filepath"];
187 echo "<a href='BeamShape.php?filepath=".$filepath."' target='page'>";?>BeamShape</a></code> class.
188
189 <br/><br/><strong>Beams:allowMomentumSpread</strong>  <input type="radio" name="17" value="on"><strong>On</strong>
190 <input type="radio" name="17" value="off" checked="checked"><strong>Off</strong>
191  &nbsp;&nbsp;(<code>default = <strong>off</strong></code>)<br/>
192 Allow the beam momenta to be smeared around their initialization
193 nominal values. 
194   
195
196 <br/><br/><table><tr><td><strong>Beams:sigmaPxA </td><td></td><td> <input type="text" name="18" value="0." size="20"/>  &nbsp;&nbsp;(<code>default = <strong>0.</strong></code>; <code>minimum = 0.</code>)</td></tr></table>
197 The width of a Gaussian distribution of the <i>p_x</i> spread of the
198 first incoming particle.
199   
200
201 <br/><br/><table><tr><td><strong>Beams:sigmaPyA </td><td></td><td> <input type="text" name="19" value="0." size="20"/>  &nbsp;&nbsp;(<code>default = <strong>0.</strong></code>; <code>minimum = 0.</code>)</td></tr></table>
202 The width of a Gaussian distribution of the <i>p_y</i> spread of the
203 first incoming particle.
204   
205
206 <br/><br/><table><tr><td><strong>Beams:sigmaPzA </td><td></td><td> <input type="text" name="20" value="0." size="20"/>  &nbsp;&nbsp;(<code>default = <strong>0.</strong></code>; <code>minimum = 0.</code>)</td></tr></table>
207 The width of a Gaussian distribution of the <i>p_z</i> spread of the
208 first incoming particle.
209   
210
211 <br/><br/><table><tr><td><strong>Beams:maxDevA </td><td></td><td> <input type="text" name="21" value="5." size="20"/>  &nbsp;&nbsp;(<code>default = <strong>5.</strong></code>; <code>minimum = 0.</code>)</td></tr></table>
212 The triply Gaussian distribution <i>(p_x, p_y, p_z)</i> is restricted to 
213 a maximal total deviation from the nominal values <i>(p_x0, p_y0, p_z0)</i>
214 for the first incoming particle, like
215 <br/><i>
216 (p_x - p_x0)^2/sigma_px^2 + (p_y - p_y0)^2/sigma_py^2 +
217 (p_z - p_z0)^2/sigma_pz^2 < maxDev^2
218 </i><br/>
219 (Note the absence of a factor 2 in the denominator, unlike the Gaussians 
220 used to pick <i>(p_x, p_y, p_z)</i>.) 
221   
222
223 <br/><br/><table><tr><td><strong>Beams:sigmaPxB </td><td></td><td> <input type="text" name="22" value="0." size="20"/>  &nbsp;&nbsp;(<code>default = <strong>0.</strong></code>; <code>minimum = 0.</code>)</td></tr></table>
224 The width of a Gaussian distribution of the <i>p_x</i> spread of the
225 second incoming particle.
226   
227
228 <br/><br/><table><tr><td><strong>Beams:sigmaPyB </td><td></td><td> <input type="text" name="23" value="0." size="20"/>  &nbsp;&nbsp;(<code>default = <strong>0.</strong></code>; <code>minimum = 0.</code>)</td></tr></table>
229 The width of a Gaussian distribution of the <i>p_y</i> spread of the
230 second incoming particle.
231   
232
233 <br/><br/><table><tr><td><strong>Beams:sigmaPzB </td><td></td><td> <input type="text" name="24" value="0." size="20"/>  &nbsp;&nbsp;(<code>default = <strong>0.</strong></code>; <code>minimum = 0.</code>)</td></tr></table>
234 The width of a Gaussian distribution of the <i>p_z</i> spread of the
235 second incoming particle.
236   
237
238 <br/><br/><table><tr><td><strong>Beams:maxDevB </td><td></td><td> <input type="text" name="25" value="5." size="20"/>  &nbsp;&nbsp;(<code>default = <strong>5.</strong></code>; <code>minimum = 0.</code>)</td></tr></table>
239 The triply Gaussian distribution <i>(p_x, p_y, p_z)</i> is restricted to 
240 a maximal total deviation from the nominal values <i>(p_x0, p_y0, p_z0)</i>, 
241 for the second incoming particle, like
242 <br/><i>
243 (p_x - p_x0)^2/sigma_px^2 + (p_y - p_y0)^2/sigma_py^2 +
244 (p_z - p_z0)^2/sigma_pz^2 < maxDev^2
245 </i><br/>
246 (Note the absence of a factor 2 in the denominator, unlike the Gaussians 
247 used to pick <i>(p_x, p_y, p_z)</i>.) 
248   
249
250 <h3>Beam interaction vertex</h3>
251
252 On this page you can set the spread of the interaction vertex according to 
253 a simple Gaussian distribution. If you instead want a more sophisticated 
254 parametrization, you can write and link your own 
255 <code><?php $filepath = $_GET["filepath"];
256 echo "<a href='BeamShape.php?filepath=".$filepath."' target='page'>";?>BeamShape</a></code> class.
257
258 <br/><br/><strong>Beams:allowVertexSpread</strong>  <input type="radio" name="26" value="on"><strong>On</strong>
259 <input type="radio" name="26" value="off" checked="checked"><strong>Off</strong>
260  &nbsp;&nbsp;(<code>default = <strong>off</strong></code>)<br/>
261 Allow the interaction vertex of the two colliding beams to be smeared.
262 If off, then the vertex is set to be the origin.
263   
264
265 <br/><br/><table><tr><td><strong>Beams:sigmaVertexX </td><td></td><td> <input type="text" name="27" value="0." size="20"/>  &nbsp;&nbsp;(<code>default = <strong>0.</strong></code>; <code>minimum = 0.</code>)</td></tr></table>
266 The width of a Gaussian distribution of the <i>x</i> location of the
267 interaction vertex. 
268   
269
270 <br/><br/><table><tr><td><strong>Beams:sigmaVertexY </td><td></td><td> <input type="text" name="28" value="0." size="20"/>  &nbsp;&nbsp;(<code>default = <strong>0.</strong></code>; <code>minimum = 0.</code>)</td></tr></table>
271 The width of a Gaussian distribution of the <i>y</i> location of the
272 interaction vertex. 
273   
274
275 <br/><br/><table><tr><td><strong>Beams:sigmaVertexZ </td><td></td><td> <input type="text" name="29" value="0." size="20"/>  &nbsp;&nbsp;(<code>default = <strong>0.</strong></code>; <code>minimum = 0.</code>)</td></tr></table>
276 The width of a Gaussian distribution of the <i>z</i> location of the
277 interaction vertex. 
278   
279
280 <br/><br/><table><tr><td><strong>Beams:maxDevVertex </td><td></td><td> <input type="text" name="30" value="5." size="20"/>  &nbsp;&nbsp;(<code>default = <strong>5.</strong></code>; <code>minimum = 0.</code>)</td></tr></table>
281 The triply Gaussian distribution of interaction vertex position 
282 <i>(x, y, z)</i> is restricted to a maximal total deviation from the 
283 origin, like
284 <br/><i>
285 x^2/sigma_x^2 + y^2/sigma_y^2 + z^2/sigma_z^2 < maxDevVertex^2
286 </i><br/>
287 (Note the absence of a factor 2 in the denominator, unlike the Gaussians 
288 used to pick <i>(x, y, z)</i>.) 
289   
290
291 <br/><br/><table><tr><td><strong>Beams:sigmaTime </td><td></td><td> <input type="text" name="31" value="0." size="20"/>  &nbsp;&nbsp;(<code>default = <strong>0.</strong></code>; <code>minimum = 0.</code>)</td></tr></table>
292 The width of a Gaussian distribution of the collision time (in units of
293 mm/c). Note that, if the above space parametrization is viewed as the
294 effect of two incoming beams along the <i>+-z</i> axis, with each beam
295 having a Gaussian spread, then the spread of the time would also become 
296 a Gaussian with the same width as the <i>z</i> one (times the 
297 velocity of the beams, which we expect is close to unity). For flexibility
298 we have not enforced any such relation, however. 
299   
300
301 <br/><br/><table><tr><td><strong>Beams:maxDevTime </td><td></td><td> <input type="text" name="32" value="5." size="20"/>  &nbsp;&nbsp;(<code>default = <strong>5.</strong></code>; <code>minimum = 0.</code>)</td></tr></table>
302 The collision time is restricted to be in the range 
303 <i>|t| &lt; sigma_t * maxDevTime</i>. 
304   
305
306 <p/>
307 The distributions above are all centered at the origin. It is also 
308 possible to shift the above distributions to be centered around another
309 nominal position. You must have <code>Beams:allowVertexSpread = on</code>
310 to use this possibility.
311
312 <br/><br/><table><tr><td><strong>Beams:offsetVertexX </td><td></td><td> <input type="text" name="33" value="0." size="20"/>  &nbsp;&nbsp;(<code>default = <strong>0.</strong></code>)</td></tr></table>
313 The <i>x</i> location of the interaction vertex is centered at this value. 
314   
315
316 <br/><br/><table><tr><td><strong>Beams:offsetVertexY </td><td></td><td> <input type="text" name="34" value="0." size="20"/>  &nbsp;&nbsp;(<code>default = <strong>0.</strong></code>)</td></tr></table>
317 The <i>y</i> location of the interaction vertex is centered at this value. 
318   
319
320 <br/><br/><table><tr><td><strong>Beams:offsetVertexZ </td><td></td><td> <input type="text" name="35" value="0." size="20"/>  &nbsp;&nbsp;(<code>default = <strong>0.</strong></code>)</td></tr></table>
321 The <i>z</i> location of the interaction vertex is centered at this value. 
322   
323
324 <br/><br/><table><tr><td><strong>Beams:offsetTime </td><td></td><td> <input type="text" name="36" value="0." size="20"/>  &nbsp;&nbsp;(<code>default = <strong>0.</strong></code>)</td></tr></table>
325 The time <i>t</i> of the interaction vertex is centered at this value. 
326   
327
328 <input type="hidden" name="saved" value="1"/>
329
330 <?php
331 echo "<input type='hidden' name='filepath' value='".$_GET["filepath"]."'/>"?>
332
333 <table width="100%"><tr><td align="right"><input type="submit" value="Save Settings" /></td></tr></table>
334 </form>
335
336 <?php
337
338 if($_POST["saved"] == 1)
339 {
340 $filepath = $_POST["filepath"];
341 $handle = fopen($filepath, 'a');
342
343 if($_POST["1"] != "2212")
344 {
345 $data = "Beams:idA = ".$_POST["1"]."\n";
346 fwrite($handle,$data);
347 }
348 if($_POST["2"] != "2212")
349 {
350 $data = "Beams:idB = ".$_POST["2"]."\n";
351 fwrite($handle,$data);
352 }
353 if($_POST["3"] != "1")
354 {
355 $data = "Beams:frameType = ".$_POST["3"]."\n";
356 fwrite($handle,$data);
357 }
358 if($_POST["4"] != "14000.")
359 {
360 $data = "Beams:eCM = ".$_POST["4"]."\n";
361 fwrite($handle,$data);
362 }
363 if($_POST["5"] != "7000.")
364 {
365 $data = "Beams:eA = ".$_POST["5"]."\n";
366 fwrite($handle,$data);
367 }
368 if($_POST["6"] != "7000.")
369 {
370 $data = "Beams:eB = ".$_POST["6"]."\n";
371 fwrite($handle,$data);
372 }
373 if($_POST["7"] != "0.")
374 {
375 $data = "Beams:pxA = ".$_POST["7"]."\n";
376 fwrite($handle,$data);
377 }
378 if($_POST["8"] != "0.")
379 {
380 $data = "Beams:pyA = ".$_POST["8"]."\n";
381 fwrite($handle,$data);
382 }
383 if($_POST["9"] != "7000.")
384 {
385 $data = "Beams:pzA = ".$_POST["9"]."\n";
386 fwrite($handle,$data);
387 }
388 if($_POST["10"] != "0.")
389 {
390 $data = "Beams:pxB = ".$_POST["10"]."\n";
391 fwrite($handle,$data);
392 }
393 if($_POST["11"] != "0.")
394 {
395 $data = "Beams:pyB = ".$_POST["11"]."\n";
396 fwrite($handle,$data);
397 }
398 if($_POST["12"] != "-7000.")
399 {
400 $data = "Beams:pzB = ".$_POST["12"]."\n";
401 fwrite($handle,$data);
402 }
403 if($_POST["13"] != "void")
404 {
405 $data = "Beams:LHEF = ".$_POST["13"]."\n";
406 fwrite($handle,$data);
407 }
408 if($_POST["14"] != "void")
409 {
410 $data = "Beams:LHEFheader = ".$_POST["14"]."\n";
411 fwrite($handle,$data);
412 }
413 if($_POST["15"] != "off")
414 {
415 $data = "Beams:newLHEFsameInit = ".$_POST["15"]."\n";
416 fwrite($handle,$data);
417 }
418 if($_POST["16"] != "on")
419 {
420 $data = "Beams:readLHEFheaders = ".$_POST["16"]."\n";
421 fwrite($handle,$data);
422 }
423 if($_POST["17"] != "off")
424 {
425 $data = "Beams:allowMomentumSpread = ".$_POST["17"]."\n";
426 fwrite($handle,$data);
427 }
428 if($_POST["18"] != "0.")
429 {
430 $data = "Beams:sigmaPxA = ".$_POST["18"]."\n";
431 fwrite($handle,$data);
432 }
433 if($_POST["19"] != "0.")
434 {
435 $data = "Beams:sigmaPyA = ".$_POST["19"]."\n";
436 fwrite($handle,$data);
437 }
438 if($_POST["20"] != "0.")
439 {
440 $data = "Beams:sigmaPzA = ".$_POST["20"]."\n";
441 fwrite($handle,$data);
442 }
443 if($_POST["21"] != "5.")
444 {
445 $data = "Beams:maxDevA = ".$_POST["21"]."\n";
446 fwrite($handle,$data);
447 }
448 if($_POST["22"] != "0.")
449 {
450 $data = "Beams:sigmaPxB = ".$_POST["22"]."\n";
451 fwrite($handle,$data);
452 }
453 if($_POST["23"] != "0.")
454 {
455 $data = "Beams:sigmaPyB = ".$_POST["23"]."\n";
456 fwrite($handle,$data);
457 }
458 if($_POST["24"] != "0.")
459 {
460 $data = "Beams:sigmaPzB = ".$_POST["24"]."\n";
461 fwrite($handle,$data);
462 }
463 if($_POST["25"] != "5.")
464 {
465 $data = "Beams:maxDevB = ".$_POST["25"]."\n";
466 fwrite($handle,$data);
467 }
468 if($_POST["26"] != "off")
469 {
470 $data = "Beams:allowVertexSpread = ".$_POST["26"]."\n";
471 fwrite($handle,$data);
472 }
473 if($_POST["27"] != "0.")
474 {
475 $data = "Beams:sigmaVertexX = ".$_POST["27"]."\n";
476 fwrite($handle,$data);
477 }
478 if($_POST["28"] != "0.")
479 {
480 $data = "Beams:sigmaVertexY = ".$_POST["28"]."\n";
481 fwrite($handle,$data);
482 }
483 if($_POST["29"] != "0.")
484 {
485 $data = "Beams:sigmaVertexZ = ".$_POST["29"]."\n";
486 fwrite($handle,$data);
487 }
488 if($_POST["30"] != "5.")
489 {
490 $data = "Beams:maxDevVertex = ".$_POST["30"]."\n";
491 fwrite($handle,$data);
492 }
493 if($_POST["31"] != "0.")
494 {
495 $data = "Beams:sigmaTime = ".$_POST["31"]."\n";
496 fwrite($handle,$data);
497 }
498 if($_POST["32"] != "5.")
499 {
500 $data = "Beams:maxDevTime = ".$_POST["32"]."\n";
501 fwrite($handle,$data);
502 }
503 if($_POST["33"] != "0.")
504 {
505 $data = "Beams:offsetVertexX = ".$_POST["33"]."\n";
506 fwrite($handle,$data);
507 }
508 if($_POST["34"] != "0.")
509 {
510 $data = "Beams:offsetVertexY = ".$_POST["34"]."\n";
511 fwrite($handle,$data);
512 }
513 if($_POST["35"] != "0.")
514 {
515 $data = "Beams:offsetVertexZ = ".$_POST["35"]."\n";
516 fwrite($handle,$data);
517 }
518 if($_POST["36"] != "0.")
519 {
520 $data = "Beams:offsetTime = ".$_POST["36"]."\n";
521 fwrite($handle,$data);
522 }
523 fclose($handle);
524 }
525
526 ?>
527 </body>
528 </html>
529
530 <!-- Copyright (C) 2012 Torbjorn Sjostrand -->