Privacy and Security Notice

heep

Correction to my hee'p method 03/28/06 (rev. 03/29/06)


In setting all of the offsets in the analyzer to zero as described here, I set some offsets to zero that I shouldn't have touched.  These are the offsets in hmsflags.param that look like this:
; Do not to change these values, since these are the zero order
;  CMOP matrix elements. If you do change then your hms sieve
;  plots will be screwed up.
  hdelta_offset = 0.       ; (%)   hdelta_tar = hdelta_tar + hdelta_offset
  htheta_offset = 4.83e-5  ; (rad) hyp_tar = hyp_tar + htheta_offset
  hphi_offset   = -4.73e-3 ; (rad) hxp_tar = hxp_tar + hphi_offset

Then, when I fitted the offsets in my heep analysis, I put the fitted offsets back in as an h_oop_central offset.


So, I repeat the heep analysis again, starting with setting all of the adjustable and previously determined offsets to zero.

Plots showing the positions of the peaks with all offsets

Table of the offsets from the hee'p analysis
Ebeam(GeV)
 P hms(GeV)
 TH hms(deg)
 P sos(GeV)
 TH sos(deg)
 dE
 dth_ss
 dp_ss
 dth_hs
 dp_hs
 dssxptar
 dhsxptar
2.038
 -1.692
 25.005
 0.874
 54.90
 0.0
 0.3
 0.8
 1.0
 1.1
 -0.0035 -0.0011
5.012
 4.116
 18.70
 -1.7227
 50.00
 0.0
 -0.9
 10.5
 0.5
 1.8
 -0.0006 -0.0011
5.766
 4.867
 16.19
 -1.6868
 51.61
 0.0
 -0.4
 8.4
 0.0
 0.6
 -0.0039
-0.0011
4.021
 3.232
 22.15
 -1.588
 50.00
 0.0
 -0.8
 3.5
 1.2
 1.5
 -0.0023
-0.0011
The format of the offsets "d[variablename]" are the same as those in "heepcheck.f".  For example, for the electron in the sos
         p_e = ssp*(1-dp_ss/1000.0)
         th_e = mysstheta-dth_ss/1000.0  [rad]
         p_p = hsp*(1-dp_hs/1000.0) 
         th_p = myhstheta-dth_hs/1000.0   [rad]
         newE=e_beam*(1+dE/1000.0)
         newssxptar= ssxptar- dssxptar  [rad]
         newhsxptar=hsxptar-dhsxptar  [rad]


This plot below is the sos fieldcorr (black points showing dp_ss) with Tanja's results (red points).  It can be seen that Tanja's fit works very well, and no other parameterization is required.  This parameterization is -0.36-0.0632*exp(spcentral**3.004).


The work on the oop offset will now be repeated below, similar to what was done here.

Plots showing the positions of the peaks without the xptar offsets

xptar_sos -(0.43-1.14*hsp/ssp)*d + P_hms/P_sos*xptar_hms = 0

  Ebeam  Phms/Psos    Centriod of ssxptar+(hsp/ssp)*hsxptar [mrad]
 -----  ----------   -----------------------------------------
 2038      1.94             -5.64
 4021      2.03             -4.48
 5012      2.89             -3.32
 5766      2.33             -7.08

It is difficult to draw a straight line through these points, even if the 5012 MeV point is thrown out.


After parameterizing dp_ss and using a global offset for dssxptar, these are the final hee'p offsets that were inserted into the Analyzer.
Ebeam(GeV)
 P hms(GeV)
 TH hms(deg)
 P sos(GeV)
 TH sos(deg)
 dE
 dth_ss
 dp_ss
 dth_hs
 dp_hs
 dssxptar
 dhsxptar
2.038
 -1.692
 25.005
 0.874
 54.90
 0.0
 0.3
 0.48
 1.0
 1.1
 -0.0032 -0.0011
5.012
 4.116
 18.70
 -1.7227
 50.00
 0.0
 -0.9
 11.0
 0.5
 1.8
 -0.0032 -0.0011
5.766
 4.867
 16.19
 -1.6868
 51.61
 0.0
 -0.4
 8.11
 0.0
 0.6
 -0.0032
-0.0011
4.021
 3.232
 22.15
 -1.588
 50.00
 0.0
 -0.8
 3.85
 1.2
 1.5
 -0.0032
-0.0011
The position of the peaks after inserting these offsets into replay and the analyzer are here.