1
|
Hi again,
|
2
|
|
3
|
PeterB agreed with the idea of repeating the fit with a constraint on the beam
|
4
|
energy offset based on the arc energy measurement uncertainty. The new results
|
5
|
are below.
|
6
|
|
7
|
According to the log, each measurement has an uncertainty of about 2.6 MeV,
|
8
|
independent of energy at least for the range 2.7-4.9 GeV. That corresponds to
|
9
|
a fractional uncertainty of 0.59-0.87E-3, so I took +/-0.7E-3 as the constraint
|
10
|
on the beam energy offset dE.
|
11
|
|
12
|
However, one of the measurements, had a significant issue in that the arc
|
13
|
energy tool gave an energy of 2.7832 GeV, which was surprisingly far from the
|
14
|
expected 2.75 GeV. After some discussion in June-2019, MarkJ came up with a
|
15
|
revised estimate of 2.74996 GeV, which I think Vijay is using his analysis. So
|
16
|
I'll assume the uncertainty on this measurement is the same as the others for
|
17
|
now, but this might need revisting.
|
18
|
|
19
|
I also fixed an error, which was that the Heepcheck derivatives were calculated
|
20
|
for a wrong setting for one of the energies.
|
21
|
|
22
|
With this beam energy offset constraint, the global fit offset results look
|
23
|
quite a bit more reasonable. As expected, the dpp(SHMS) and dpe(HMS) offsets
|
24
|
are also now quite a bit smaller. The dthp(SHMS) angle offset is also quite a
|
25
|
bit smaller, while the dthe(HMS) angle offset has grown slightly.
|
26
|
|
27
|
Finally, the chi-square per energy is also more uniform. The 3.6 GeV Heep
|
28
|
setting still has the largest chi-square contribution to the global fit, but by
|
29
|
a lesser amount than before. Of course, the overall chi-square is slightly
|
30
|
larger now, as the constrained fit no longer corresponds to the absolutely best
|
31
|
fit, but the difference is not huge.
|
32
|
|
33
|
If you have questions or comments, please let me know. We can discuss these
|
34
|
further at Thursday's meeting.
|
35
|
|
36
|
Thanks, Garth
|
37
|
|
38
|
|
39
|
*** RESULTS ***
|
40
|
! units are 0.1% for momenta/energy, 1 mrad for angles
|
41
|
|
42
|
Global Offsets:
|
43
|
dthe 1.3000 dpe 0.0000 dthp 0.7000 dpp 1.8000
|
44
|
dE: -0.7000 0.7000 0.7000 0.0000 -0.2000
|
45
|
chisq: 48.357 5.608 17.334 9.608 5.458 10.350
|
46
|
|
47
|
Offsets for Energy: 2750.0
|
48
|
dthe 2.3000 dpe 0.0000 dthp 0.3000 dpp 4.8000 dE 0.7000 chisq 3.201
|
49
|
|
50
|
Offsets for Energy: 3660.2
|
51
|
dthe 0.4000 dpe 0.1000 dthp 0.1000 dpp 1.4000 dE 0.6000 chisq 0.082
|
52
|
|
53
|
Offsets for Energy: 3834.9
|
54
|
dthe 0.1000 dpe 0.6000 dthp 0.7000 dpp 0.0000 dE -0.1000 chisq 0.603
|
55
|
|
56
|
Offsets for Energy: 4559.7
|
57
|
dthe 0.4000 dpe 1.6000 dthp 1.8000 dpp 0.0000 dE -0.1000 chisq 2.795
|
58
|
|
59
|
Offsets for Energy: 4932.0
|
60
|
dthe 0.9000 dpe 0.8000 dthp 2.2000 dpp 0.2000 dE -0.5000 chisq 4.194
|
61
|
|
62
|
*** Expected shifts ***
|
63
|
units are 0.1% for momenta/energy, 1 mrad for angles
|
64
|
|
65
|
--- Original --- Global --- Individual Fit
|
66
|
|
67
|
Energy 2750.0
|
68
|
W 6.3945 1.2125 0.5605
|
69
|
EM 4.0667 -0.5763 -1.2563
|
70
|
PMZ 4.0486 1.5776 1.1436
|
71
|
PMX -0.5784 -1.1474 0.0186
|
72
|
|
73
|
Energy 3660.2
|
74
|
W 1.0657 -3.7233 0.0367
|
75
|
EM 0.8051 -0.4669 -0.1919
|
76
|
PMZ 0.8696 1.4406 0.2086
|
77
|
PMX -0.8696 1.0854 0.0244
|
78
|
|
79
|
Energy 3834.9
|
80
|
W 3.1973 -2.1227 0.1733
|
81
|
EM 1.0916 -0.6014 -0.5094
|
82
|
PMZ 1.1614 1.2994 0.5594
|
83
|
PMX -0.3871 1.7469 0.0269
|
84
|
|
85
|
Energy 4559.7
|
86
|
W 10.6576 1.8176 0.3706
|
87
|
EM 3.3951 -1.3749 -1.1409
|
88
|
PMZ 2.5072 0.4362 1.1642
|
89
|
PMX -0.7163 -0.2713 -0.0173
|
90
|
|
91
|
Energy 4932.0
|
92
|
W 12.7891 2.3761 0.5001
|
93
|
EM 4.0027 -1.3573 -1.4443
|
94
|
PMZ 2.7100 0.7800 1.3630
|
95
|
PMX -0.7743 -1.5013 -0.0113
|
96
|
|
97
|
Mean 2.7718 0.0315 0.0789
|
98
|
|