1 | ! #184 |
---|
2 | |
---|
3 | process W1 = u, D -> E1, n1 |
---|
4 | process W2 = D, u -> E1, n1 |
---|
5 | |
---|
6 | compile |
---|
7 | |
---|
8 | sqrts = 1000 |
---|
9 | |
---|
10 | !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
---|
11 | iterations = 1:40000 |
---|
12 | ! setting calls to >=40000 => integrate gives the output |
---|
13 | ! | Contributing Helicity Combinations: 2 of 16 |
---|
14 | ! | Threshold: amp / avg > 0.22E-05 = 0.10E+11 * epsilon() |
---|
15 | ! | 1: -1 1 -1 -1 |
---|
16 | ! | 2: -1 1 1 -1 |
---|
17 | ! | Contributing Helicity Combinations: 2 of 16 |
---|
18 | ! | Threshold: amp / avg > 0.22E-05 = 0.10E+11 * epsilon() |
---|
19 | ! | 1: 1 -1 -1 -1 |
---|
20 | ! | 2: 1 -1 1 -1 |
---|
21 | !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
---|
22 | !iterations = 1:39999 |
---|
23 | ! setting calls to <40000 => integrate gives the output |
---|
24 | ! | Contributing Helicity Combinations: 2 of 16 |
---|
25 | ! | Threshold: amp / avg > 0.22E-05 = 0.10E+11 * epsilon() |
---|
26 | ! | 1: -1 1 -1 -1 |
---|
27 | ! | 2: -1 1 1 -1 |
---|
28 | ! | Contributing Helicity Combinations: 3 of 16 |
---|
29 | ! | Threshold: amp / avg > 0.22E-05 = 0.10E+11 * epsilon() |
---|
30 | ! | 1: 1 -1 -1 -1 |
---|
31 | ! | 2: 1 -1 1 -1 |
---|
32 | ! | 3: 1 -1 1 1 |
---|
33 | !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
---|
34 | ! changing number of iterations has no effect on this |
---|
35 | ! tested with r1671 compiled with nagfor and gfotran |
---|
36 | !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
---|
37 | |
---|
38 | beams = p, p -> lhapdf |
---|
39 | |
---|
40 | integrate(W1) |
---|
41 | |
---|
42 | integrate(W2) |
---|