CH HPGe Efficiency

From New IAC Wiki
Revision as of 14:24, 31 October 2022 by Harpconn (talk | contribs) (→‎Efficiency Measurements Summer 2022)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to navigation Jump to search

Efficiency Measurements Summer 2022

This page is mostly correct, but will need to be updated to reflect the final methodologies.

Electronics Chain

Click here to view the electronics chain used for these efficiency measurements.

Experimental Setup

The Eu152 source provided by Roy Dunker was used to determine the efficiency and energy calibration for future experiments. As of 7/26/22, the Eu152 source had decayed to an activity of approximately 37.7μCi (~1.4•106 Bq).

Measurements

Energy Calibration

A linear regression was performed on the (ADC Channel Number, Peak Energy) data set. The data set is included in the table below.

Index ADC Channel Number Eu152 Peak Energy (KeV)
1 689 244.70
2 937 344.28
3 1105 411.12
4 1187 443.96
5 2025 778.90
6 2247 867.37
7 2490 964.08
8 2796 1085.9
9 2862 1112.1
10 3607 1408.0

The conversion from ADC Channel to KeV is as follows:

[math] \text{Photon Energy (KeV)}=0.398344*(\text{Channel Number})-28.5131 [/math]

The KeV per bin can be found by taking the derivative of the energy conversion, leading to a constant 0.398 KeV per bin.

Note that while the ADC is meant to be linear, there may be a nonlinear response as a function of energy, it may be wise to switch to a higher order energy calibration to account for this.

Background Spectrum

Run Parameters

Date Performed 7/27/22 - 7/28/22
File Name r7224.root
Run Time 82129 Seconds
Count Rate 24.97 Hz

On Axis

Run Parameters

Date Performed 7/26/22
File Name r7222.root
Run Time 726 Seconds
Count Rate 2442.06 Hz

90 Degrees Off Axis

Run Parameters

Date Performed 7/26/22
File Name r7223.root
Run Time 542 Seconds
Count Rate 2634.47 Hz

Overlays

Dead Time Analysis

Efficiency Calculations

Photon Emission Rate and Flux

The radioactive source used in the efficiency experiments (on and off axis as defined above) was 152Eu button source with an activity of ~37.7μCi on the day of irradiation.

The activity of the source was converted from 37.7μCi to Becquerels using the conversion of 37,000[math]\frac{\text{Bq}}{\text{μCi}}[/math], yielding a decay rate for the 152Eu source of ~1.3949•106 [math]\frac{\text{decays}}{\text{sec}}[/math].

The γ emission rate into 4π was found by totaling the intensities of photon lines produced via the EC and β- decay channels of 152Eu, for a result of 159.69%. Data sourced from IAEA Nuclear Data Section, decay radiation subsection for 152Eu.

γ-rate = 1.3949×106 [math]\frac{decays}{sec}[/math] × 1.5969 [math]\frac{photons}{decay}[/math] = 2.22752×106 [math]\frac{{photons}}{{sec}}[/math] into 4π.

The flux of photons 45cm from the source was calculated by using the γ-rate from above and the surface area of a sphere at 45cm.

flux (Φ) = [math]\frac{γ-rate}{4π(45^2)}[/math] = [math]\frac{2.22752 {×10^6}} {25446.9} \frac{photons}{sec\ cm^2}[/math] = 87.5358 [math]\frac{photons}{sec\ cm^2}[/math]

Intrinsic Efficiency

Intrinsic efficiency as a function of energy is calculated by integrating the energy peak of interest and dividing by the number of photons crossing the detector plane at that energy. I used the method outlined below to calculate the intrinsic efficiency as a function of energy for the HPGe detector in the clean room. The number of photons crossing the detector plane was found by multiplying the flux from above by the cross sectional area of the detector crystal in each orientation. The rate of the peak energy photons uses

γ-rate into detectorOn Axis = flux (Φ) × X-area = 87.5358 [math]\frac{photons}{sec\ cm^2}[/math] × π (4.9)2cm2 = 6602.8 [math]\frac{photons}{sec}[/math]

γ-rate into detectorOff Axis = flux (Φ) × X-area = 87.5358 [math]\frac{photons}{sec\ cm^2}[/math] × (9.8)2cm2 = 8406.94 [math]\frac{photons}{sec}[/math]

Peak Energy (KeV) Line Intensity (%) [math]\frac{Photons}{sec}[/math] Across Detector Plane On Axis Efficiency (%) Off Axis Efficiency (%)
344.28 26.58 1755.02 11.0864 9.0585
411.12 2.23 147.242 11.1684 9.60532
443.96 3.12 206.007 11.5711 9.45957
778.9 12.96 855.722 8.74022 7.44913
867.37 4.23 279.298 8.26167 7.12743
964.08 14.62 965.329 8.11188 6.78583

Absolute Efficiency

The absolute efficiency as a function of energy is calculated by integrating the energy peak of interest and dividing by photons of that energy emitted into 4π from the source.

γ-rateEi = 2.22752×106 [math]\frac{{photons}}{{sec}}[/math] × ProbEi


Peak Energy (KeV) Line Intensity (%) [math]\frac{Photons}{sec}[/math] Emitted On Axis Efficiency (%) Off Axis Efficiency (%)
344.28 26.58 592074 0.0328624 0.034188
411.12 2.23 49673.6 0.0331053 0.0362518
443.96 3.12 69498.5 0.0342989 0.0357017
778.9 12.96 288686 0.0259078 0.028114
867.37 4.23 94223.9 0.0244892 0.0268999
964.08 14.62 325663 0.0240452 0.0256106

Past Work

Most likely incorrect

Eu152 Calibration

  • HPGe calibrated using Eu152 source
CH Eu152 Calibration Spectrum.png

As of 9/28/21

The calibration used to map from ADC channel number to KeV is:

[math]Photon\ Energy\ =\ -36.2928\ +\ 0.449519\ x\ -\ 0.000000861277\ x^2[/math]

Where x = channel 9 of the ADC and is input as PADC.PADC785N[9].

Total Efficiency

This is total efficiency using the Gadolinium foil as the source.

[math] (line\ intensity)\ (incident\ thermal\ neutrons\ per\ second)\ (\sigma_{total})\ (MCNP\ fudge\ factor\ from\ Paul\ approx\ 4)\ =\ Theoretical\ Rate [/math]

For the ease of calculation, I will use the outgoing photon rate from the Gd foil as tabulated by MCNP, thanks to Paul Stonaha. The number given was into 2pi, towards the HPGe detector, multiplying by 2 should give the rate into 4pi.

The 897 KeV peak found in the Gd in/out measurements can be attributed to two different lines created by neutron capture on Gd157. The high purity germanium detector used in this experiment would not be able to resolve the two distinct lines as the resolution is not fine enough. However, adding the intensities together should yield the theoretical rate for the peak measured at that energy.

The first peak is 897.506 KeV with a relative intensity of 17.846%, according the capgam nndc data repository.

The second peak is 897.62 KeV with a relative intensity of 16%, according the capgam nndc data repository.

(0.17846 + 0.16) (2*1374Hz) (234000) = Rate

234000 is the total thermal neutron capture cross section for 157Gd

Note: the fudge factor comes from the nucleus being excited and then de-exciting from the neutron absorption and emits on average, 4 photons per neutron capture

Then find detected rate using HPGe

Divide (detected rate) by (theoretical rate) for total detector efficiency

An easier way to calculate the detector efficiency is

[math]Det.\ Eff.\ =\ \frac{Detected\ Rate}{(Normalized\ Line\ Intensity)(Theoretical\ Total\ Rate)}[/math]

Theoretical total rate = 2748 Hz from simulations

Line Energy Line Weighting Detected Rate Detector Efficiency
897 KeV (2 lines) 0.03967382135 0.357 Hz 0.33% [math]\pm[/math] 0.03
1107 KeV 0.02885336857 0.264 Hz 0.33% [math]\pm[/math] 0.03
~1186 KeV (2 lines) 0.05145546699 0.511 Hz 0.36% [math]\pm[/math] 0.03

Previous Page

CaGaS Phase II

PAS_Experiment

PAS HPGe Detector