Analysis of Compton Scattering data by MCC.  July/August 2017.

Setup Parameters


{FIXME ${/download/attachments/164080117/Table_layout.png?version=1&modificationDate=1501004539000&api=v2}$ {FIXME ${/download/attachments/164080117/Compton_pig.png?version=1&modificationDate=1502292371000&api=v2}$ {FIXME ${/download/attachments/164080117/Compton_geometry.png?version=2&modificationDate=1502292505000&api=v2}$ Set PMT to 1500v

PHA preamp gain.  Coarse = 8  Fine = 1.68

Detector set to 90deg.

PMT remained at HV = 1500v throughout all measurements.

PHA Channel measurements made by creating a region of interest in USX software and reading the centroid.  ROI was chosen to encompass symmetric portion of the peak as shown for the 356keV peak of Ba133 in the figure below.

{FIXME ${/download/attachments/164080117/Ba133_roi_selection.PNG?version=1&modificationDate=1501090462000&api=v2}$

Pulse Height to Energy Calibration


Calibration data using Cs137, Na22 and Ba133 sources were taken over multiple days.

Calibration Runs

PMT Energy Calibration 072017

PHA Ch E(keV)
841 662
658 511
466.4 356
110.1 81
396.9 302

PMT Energy Calibration  072417

PHA Ch E(keV)
824 662
643 511
456 356
387 302
106 81

PMT Energy Calibration  072517

PHA Ch E(keV)
831 662
646 511
459 356
388 302
107 81

PMT Energy Calibration  080717.  PMT has been on continuously since 072017.

PHA Ch E(keV)
809 662
629 511
446 356
379 302
104 81

Fit all three calibration data sets to a straight line.  

PHA Ch = m * Energy + b

Best fit values from the three sets of calibration data.

Date b m
072017 13.44 1.26
072417 10.99 1.23
072517 10.26 1.24
080717 10.09* 1.21*

*Averages do not include last data pair which was taken over a week after the scattering data.

Y-intercept b: Avg = 11.6  STD = 1.4.

Slope m: Avg = 1.24  STD = 0.01.

So to get energy from PHA channel,

E = ( P - b ) / m

dE = E * ( db / ( P-b ) + dm / m )  Assuming negligible uncertainty in peak location.

{FIXME ${/download/attachments/164080117/PHA_calib.png?version=1&modificationDate=1501101041000&api=v2}$

Gain Dependence on Angle

Does pmt gain change with angle (i.e. orientation in earths magnetic field).

PHA Ch E(keV) Angle (Deg)
832 662 90
832 662 120
833.8 662 60
836 662 30
838 662 0
838 662 140

Measure again on 080917.  Use 662keV line from Cs137 calibration source.  Pig closed and scatterer removed.

Angle PHA Ch Gross Counts in peak in 20s Counts in full spectrum
0 810 38301 89618
30 809 38183 89805
60 809 38169 89826
90 810 38069 89529
120 810 37830 89205
140 810 37851 89084

Uncertainty in peak centroid measurement.


Estimate PHA uncertainty by selecting ROI multiple times.

Include making some obviously non-symmetric selections.  The following data table contains 16 separate measurements of the centroid of the 356keV peak of Ba133.  Note that spread in measurements is exaggerated as I intentionally made my ROI selections sloppy. 

466 466 465 466
467 466 467 467
466 466 467 466
466 466 466 467

Average = 466.18

STD = 0.58

Uncertainty in the selection of the peak location is 0.1%.  

Test effect of run length on centroid estimation.

To test how sensitive the centroid estimation is to number of counts in the peak two Cs137 spectra were taken.  One where the full energy peak amplitude was about 50 ( Lo counts ) and one where the full energy peak reached about 1000 counts (Hi counts).  Visually the Hi count full energy peak appears smoother and would seem to provide a better estimate of the centroid location than the Lo count spectra.  However when each spectra is analyzed by reselecting the ROI 10 times, there is stastically no difference in the value of the centroid.

Lo Count Spectra

{FIXME ${/download/attachments/164080117/Lo_counts.PNG?version=1&modificationDate=1502300929000&api=v2}$ Ten ROI centroid measurements:  824, 823, 824, 825, 823, 824, 825, 824, 824, 824.

Avg = 824.0  Std = 0.6

Hi Count Spectra

{FIXME ${/download/attachments/164080117/Hi_counts.PNG?version=1&modificationDate=1502300929000&api=v2}$ Ten ROI centroid measurements: 825, 824, 825, 825, 824, 823, 824, 825, 825, 824.

Avg = 824.4  Std = 0.7

Conclusion

For a user with an experienced eye in selecting ROI's, using the USX programs centroid calculation provides very repeatable results.  Determination of full energy peak centroid seems to be fairly insensitive to the precise boundaries of the ROI as well as the visual smoothness of the peaks.  We have probably been having students collect data for much longer than is necessary.

Scattering Data


Angle PHA Ch E (keV) dE (keV)
90 377.9 295.5 3.5
140 263.2 202.7 2.8
120 294.4 227.7 3.0
105 336.0 261.6 3.2
75 429.0 336.6 3.9
60 506.0 398.7 4.4
45 603.0 476.9 5.0
30 708.0 561.6 5.7
15 796 632.6 6.3
145 259 199.5 2.7
130 272 210.0 2.8
110 313 243.1 3.1
120 290 224.5 3.0
40 647
90 375
100 343

{FIXME ${/download/attachments/164080117/Energy_vs_angle.png?version=1&modificationDate=1501101093000&api=v2}$ 080717 - Take another data set over about 1.5 hours after the pmt has been on for a couple weeks straight.  5 minute runs for all angles.  Calibration data was taken on the same day as the scattering data.  This coupled with the fact the pmt has been held at HV since 072017 should minimize the effect of gain drift.

Angle PHA Ch E (keV) dE (keV)
15 787 642.1 1.5
30 701 571.0 1.4
45 596 484.2 1.3
60 502 406.5 1.2
75 427 344.6 1.1
90 365 285.9 1.1
105 317 253.6 1.1
120 284 226.4 1.0
135 259 205.7 1.0
140 253 200.8 1.0

Prior to the fit the measured angles were corrected for the known 0.7 degree offset in the table scale.  The only free parameter of the fit was the electron rest mass.  The fit returned a rest mass for the electron of 513.3keV ± 1.3keV.  The reduced Chi Sq was 6.2.  However the uncertainties on the full energy peak centroids are not due to statistically random fluctuations, so a reduced Chi Sq of 1 is not necessarily expected.  

Note that this result was obtained from spectra whose run times were considerably shorter than what we have considered necessary.

{FIXME ${/download/attachments/164080117/Compton_scatter_080717.png?version=1&modificationDate=1502303790000&api=v2}$

Angular Acceptance Effects


Can the non-zero angular acceptance of the detector produce a systematic shift in the measured energies as a function of angle.  We test for this possible systematic effect by recording and comparing spectra of the 511keV line from Na22 with that of 511keV compton scattered gammas.  We also do the comparison at 244keV using Eu152 in place of Na22.  Comparison of 511keV spectra from Na22 source and compton scattered photons at 40deg which corresponds to a scattered energy of ~511keV.

511keV Comparison

{FIXME ${/download/attachments/164080117/511kev_spectra.png?version=1&modificationDate=1501101135000&api=v2}$ Zooming in on the full energy peaks and fitting both to gaussian functions gives.

{FIXME ${/download/attachments/164080117/Gaussian_fits.png?version=1&modificationDate=1501533950000&api=v2}$ The table below gives the best fit parameters for a gaussian function of the form,

{FIXME ${/download/attachments/164080117/image2017-7-31%2015%3A56%3A11.png?version=1&modificationDate=1501534571000&api=v2}$.

Best fit values @ 511keV
P1 P2 P3 P4
40deg 15 9031 29 648
Na22 8 9177 23 647

Both fits give a reduces chisq of about 1.5.  While the compton scattered data show a broader full energy peak, there is no meaningful shift in the location of the centroid.

244keV Comparison

Use Eu152 to get a line at 244keV which corresponds to a compton scattered angle of 108.8.

{FIXME ${/download/attachments/164080117/244keV_comparison.png?version=1&modificationDate=1501607529000&api=v2}$ The best fit values for the gaussians are given below.  As with the 511keV energy the compton scattered data have a wider sigma, but the centroids are essentially the same.

Best fit values @ 244keV
P1 P2 P3 P4
109deg 541 29854 17 324
Eu152 71 17914 14 323

Conclusion

It does not appear that there is any meaningful shift in measured energies as a function of scattering angle.  The full energy peak of the compton scattered gammas is slightly wider which may impact measurements of the differential scattering cross section.  But there should be no effect on measurements of the compton scattered photon energy.

Test for Offset of Table Angle Scale


Scan through the beam to check the 0 deg mark.  All runs at 100s live time.  Used lead collimator in detector.  ROI encompasses entire spectrum, not just full energy peak.  ROI remained unchanged for all measurements.  Dead time never exceeded ~5%, so there should be no saturation effects.  No scatterer.

Angle Gross
-6 11491
-4 56757
-3 86190
-2 136445
0 198120
1 201009
2 191332
3 154442
4 102471
6.2 36946
7 15431

Best fit parameters of fit to Gaussian.

 {FIXME ${/download/attachments/164080117/image2017-7-31%2015%3A56%3A11.png?version=1&modificationDate=1501534571000&api=v2}$ P1 = -5092

P2 = 1.56e6

P3 = 3.0

P4 = 0.7

{FIXME ${/download/attachments/164080117/detector_scan.png?version=1&modificationDate=1501620859000&api=v2}$ The data show a 0.7° offset in the orientation of the angle scale with respect to the beam axis.  This should be added to the angle readings taken from the table scale.

Uncertainty in Angle measurement.  

Uncertainty in reading the angle markings on the table is less than ±0.5deg.  Accuracy of the scale however is unknown.