LB SageBrushWork PostDefense

From New IAC Wiki
Jump to navigation Jump to search

Mass Information

Note that only 2 samples were produced using the oven method (0.04% and 10%). The target masses are as follows

50% Se/Sage Sample

Selenium Target (50% Sagebrush Se Mix) Mass (g)
Outer Front Ni Foil 0.6948
Se Mass 0.51
Outer Se Mass 0.1
Sagebrush Leaves Mass 0.4969



10% Sample

Selenium Target (10% Sagebrush Se Mix) Mass (g)
Outer Front Ni Foil 0.2783
Se Mass 0.0523
Outer Se Mass 0.0971
Sagebrush Leaves Mass 0.5111


0.1% Se/Sage Sample

Selenium Target (0.1% Sagebrush Se Mix) Mass (g)
Outer Front Ni Foil 1.1437
Se Mass 0.0559
Outer Se Mass 0.1118
Sagebrush Leaves Mass 50.1080



0.04% Sample

Selenium Target (0.04% Sagebrush Se Mix) Mass (g)
Outer Front Ni Foil 0.5447
Se Mass 0.0020
Outer Se Mass 0.0992
Sagebrush Leaves Mass 4.9326


Time Cuts

The first step in the PAA process is to identify the time cuts used for the split run. This must be done for all samples (50%,10%,0.1%, and 0.04%).

50% Sage/Se Mixture Time Cuts

The 50% Se/Sage mixture was measure on Detector B (unshielded, so watch the SNR) on 5/24/17 for a total of 2130.540 seconds. The outer pure witness Se was measure first, and sub-run was roughly 300 seconds for the samples of interest, and 60 seconds for the Co-60 flag. The root command used to draw the timing information was

TTree* tree = MPA;

MPA->Draw("evt.Chan:evt.Sec>>hist2","evt.ADCid == 1");

which produced the histogram below

LB 50 Percent SeSage TimeCutInfo.png

This seems a little sloppy, lets try a 1D histogram instead to see if the cuts are more clear

ROOT Command: MPA->Draw("evt.Sec>>hist","ADCid==1","");

This produced the histogram below

LB 50 Percent SeSage TimeCutInfo 1DHistogram.png

This seems much cleaner and clearly shows the time cuts

The time cuts are

0-300 (Pure Se)

375-440 (Co-60)

450-700 (Mixture)

715-780 (Co-60)

800-1010 (Pure Se)

1055-1150 (Co-60)

1160-1420 (Mixture)

1425-1495 (Co-60)

1500-1800 (Pure Se)

1810-1870 (Co-60)

1880-2100 (Mixture)

10% Sage/Se Mixture Time Cuts

The 10% Se/Sage mixture was measured on Detector A (shielded) on 5/25/17 for a total of 5136.343 seconds. The mixture was measured first, and a sub-run was roughly 300 seconds for the Se samples and 60 seconds for the Co-60 flag. The root commands used were

TTree* tree = MPA;

MPA->Draw("evt.Chan:evt.Sec>>hist2","evt.ADCid == 0");

10Percent SeSage TimingCutInfo.png

10Percent SeSage TimingCutInfo 1DHistogram.png

The time cut information is as follows

0-300 (Mixture)

300-360 (Co-60)

400-640 ( Pure Se)

680-710 (Co-60)

730-1020 (Mixture)

1030-1080 (Co-60)

1100-1360 ( Pure Se)

1400-1440 (Co-60)

1480-1775 (Mixture)

1800-1840 (Co-60)

1875-2150 (Pure Se)

2190-2220 (Co-60)

2250-2550 (Mixture)

2590-2620 (Co-60)

2650-2930 (Pure Se)

2950-3000 (Co-60)

3050-3300 (Mixture)

3390-3350 (Co-60)

3400-3690 (Pure Se)

3700-3750 (Co-60)

3775-4050 (Mixture)

4060-4100 (Co-60)

4120-4400 (Pure Se)

4420-4470 (Co-60)

4480-4770 (Mixture)

4790-4820 (Co-60)

4840-5130 (Pure Se)

0.1% Sage/Se Mixture Time Cuts

The 0.1% Se/Sage mixture was measured on Detector A (shielded) on 5/24/17 for a total of 2130.540 seconds. The mixture was measured first, and a sub-run was roughly 300 seconds for the Se samples and 60 seconds for the Co-60 flag. The root commands used were

TTree* tree = MPA;

MPA->Draw("evt.Chan:evt.Sec>>hist2","evt.ADCid == 0");

which produced the histogram below

0 1 Percent SeSage TimeCutInfo.png

This seems a little sloppy, so let's repeat the procedure for the 50% sample

ROOT Command:ROOT Command: MPA->Draw("evt.Sec>>hist","ADCid==0","");

This produced the histogram below

0 1 Percent SeSage TimeCutInfo 1DHistogram.png

This is again much better than the 2D histogram.

The time cuts are

0-300 (Mixture)

306-340 (Co-60)

376-640 (Pure Se)

650-720 (Co-60)

730-1010 (Mixture)

1020-1050 (Co-60)

1100-1350 (Pure Se)

1400-1460 (Co-60)

1470-1730 (Mixture)

1735-1800 (Co-60)

1804-2140 (Pure Se)

The histogram for the time cuts of the second run is shown below

0 1 Percent SeSage TimeCutInfo 1DHistogram Run2.png

The time cuts are

0-240 (Outer Se)

245-300 (Co-60)

305-820 (Mixture)

0.04% Sample Time Cuts

The 0.04% Se/Sage mixture was measured on Detector A (shielded) on 5/25/17 for a total of 1711.891 seconds. The mixture was measured first, and a sub-run was roughly 300 seconds for the Se samples and 60 seconds for the Co-60 flag. The root commands used were

TTree* tree = MPA;

MPA->Draw("evt.Sec>>hist2","evt.ADCid == 0");


The histogram representing the time cuts for the first run are shown below

0 04Percent SeSage TimeCuts 1DHistogram.png

The time cuts for this run are

0-290 (Mixture)

295-360 (Co-60)

370-684 (Pure Se)

The histogram representing the time cuts for the second run are shown below

0 04Percent SeSage TimeCuts 1DHistogram Run2.png

The time cuts for this run are

0-300 (Mixture)

305-355 (Co-60)

360-650 (Pure Se)

660-735 (Co-60)

740-1020 (Mixture)

Gamma Spectrum Analysis

50% Sample Analysis

Below is the analysis table for the 50% Se/Sage mixture. The 103 keV line of Se-81m corresponds to a channel # of 130. Due to the sloppiness of the spectrum, I will use a linear + Gaussian fitting function and restrict my analysis to 6 channels because there are other energy lines nearby. Below is a sample spectrum for an idea of the resolution of Detector B for this measurement

LB 50Percent SeSage SampleSpectrum.png

As a comparison for the constant fit, here is a sample spectrum of the Se/Sage measured by Detector A

10Percent SeSage SampleSpectrum.png


0<t<300
Thin Window Histogram
Signal in Thin Window
Integrated Background
Signal - Background
Runtime (s)
Rate (Hz)
Integral Decay Correction (Hz)
Dead Time (%)
Dead Time Corrected Signal (Hz)
.Dat File Entry For HL Plot

10% Sample Analysis

Below is the analysis table for the 10% Se/Sage mixture. The 103 keV line of Se-81m corresponds to a channel # of 113


0<t<300 730<t<1020 1480<t<1775 2250<t<2550 3050<t<3300 3775<t<4050 4480<t<4770
Thin Window Histogram 10Percent SeSageMix 0 t 300Sec 10Chan Se81m.png 10Percent SeSageMix 730 t 1020Sec 10Chan Se81m.png 10Percent SeSageMix 1480 t 1775Sec 10Chan Se81m.png 10Percent SeSageMix 2250 t 2550Sec 10Chan Se81m.png 10Percent SeSageMix 3050 t 3300Sec 10Chan Se81m.png 10Percent SeSageMix 3775 t 4050Sec 10Chan Se81m.png 10Percent SeSageMix 4480 t 4770Sec 10Chan Se81m.png
Signal in Thin Window [math] 6.019 \times 10^5 \pm 775.82[/math] [math] 4.852 \times 10^5 \pm 696.56 [/math] [math] 4.161 \times 10^5 \pm 645.06 [/math] [math] 3.667 \times 10^5 \pm 605.56 [/math] [math] 2.645 \times 10^5 \pm 514.30 [/math] [math] 2.451 \times 10^5 \pm 495.08 [/math] [math] 2.344 \times 10^5 \pm 484.15 [/math]
Integrated Background [math] 2.956 \times 10^5 \pm 885.1 [/math] [math] 2.074 \times 10^5 \pm 728.7 [/math] [math] 1.668 \times 10^5 \pm 653.0 [/math] [math] 1.47 \times 10^5 \pm 606.0 [/math] [math] 1.115 \times 10^5 \pm 528.6 [/math] [math] 1.055 \times 10^5 \pm 418.7 [/math] [math] 9.935 \times 10^4 \pm 496 [/math]
Signal - Background [math] 3.063 \times 10^5 \pm 1177.01 [/math] [math] 2.778 \times 10^5 \pm 1008.07 [/math] [math] 2.493 \times 10^5 \pm 917.88[/math] [math] 2.197 \times 10^5 \pm 856.70 [/math] [math] 1.53 \times 10^5 \pm 737.51 [/math] [math] 1.396 \times 10^5 \pm 648.39 [/math] [math] 1.35050 \times 10^5 \pm 693.12 [/math]
Runtime (s) 300 290 295 300 250 275 290
Rate (Hz) [math] 1021 \pm 3.92 [/math] [math] 957.93 \pm 3.48 [/math] [math] 845.08 \pm 3.11 [/math] [math] 732.33 \pm 2.86 [/math] [math] 612 \pm 2.95 [/math] [math] 507.64 \pm 2.36 [/math] [math] 465.69 \pm 2.39 [/math]
Integral Decay Correction (Hz) [math] 1052.2 \pm 4.04 [/math] [math] 986.23 \pm 3.58 [/math] [math] 870.47 \pm 3.20 [/math] [math] 754.71 \pm 2.95 [/math] [math] 627.56 \pm 3.03 [/math] [math] 521.85 \pm 2.43 [/math] [math] 479.44 \pm 2.46 [/math]
Dead Time (%) 5.06 +/- 0.39 3.06 +/- 0.30 2.53 +/- 0.31 2.26 +/- 0.33 1.92 +/- 0.26 1.65 +/- 0.34 1.55 +/- 0.18
Dead Time Corrected Signal (Hz)
.Dat File Entry For HL Plot

0.1% Sample Analysis

Below is the analysis table for the 0.1% Mixture, the 103 keV line for Se-81m corresponds to a channel # of 113.


0<t<300 730<t<1010 1470<t<1730 305<t<820 (Run 2)
Thin Window Histogram 0 1PercentSeSageMix 0 t 300Sec.png 0 1PercentSeSageMix 730 t 1010Sec.png 0 1PercentSeSageMix 1470 t 1730Sec.png 0 1PercentSeSageMix 305 t 820Sec Run2.png
Signal in Thin Window No Signal Visible No Signal Visible No Signal Visible Possible Signal Visible?
Integrated Background - - - -
Signal - Background - - - -
Runtime (s) 300 280 260 515
Rate (Hz) - - - -
Integral Decay Correction (Hz) - - - -
Dead Time (%) - - - -
Dead Time Corrected Signal (Hz) - - - -
.Dat File Entry For HL Plot - - - -

0.04% Sample Analysis

0<t<290 0<t<300 (Run 2) 740<t<1020 (Run 2)
Thin Window Histogram 0 04Percent SeSage Mix 0 t 290Sec Run1.png 0 04Percent SeSage Mix 0 t 300Sec Run2.png 0 04Percent SeSage Mix 740 t 1020Sec Run2.png
Signal in Thin Window No Signal Visible No Signal Visible No Signal Visible
Integrated Background - - - -
Signal - Background - - - -
Runtime (s) 290 300 280
Rate (Hz) - - - -
Integral Decay Correction (Hz) - - - -
Dead Time (%) - - - -
Dead Time Corrected Signal (Hz) - - - -
.Dat File Entry For HL Plot - - - -

Fall 2018 Possible Sage Irradiation Activity Prediction

Isotope Search link: http://nucleardata.nuclear.lu.se/toi/

For this section I will attempt to identify the most active line within a sagebrush sample. The histograms that follow in this section will be weighted by the mass of the sage that was used for the sample. I will use the 10% Sage/Se mixture because it was processed using the oven method for dehydration. The mass was 0.5111g which gives a weighting factor of 1.9566 [math] g^{-1} [/math] Below is a full spectrum for the time cut of 0 <t< 300 seconds.

10Percent SeSageMix 0 t 300Sec FullSpectrum PreIrrAnalysis.png

10Percent SeSageMix 0 t 300Sec 0 E 600 PreIrrAnalysis.png 10Percent SeSageMix 0 t 300Sec 600 E 1500 PreIrrAnalysis.png

It seems the two most prominent lines are around 150 keV and 250 keV, so let's try to study the decay of those first.

144.5 keV Line

Note: the x-axis for the first 2 histograms should be Energy (keV). This line is NOT present in the Pure Se Sample, so it should be a line produced by the sage itself.

0<t<300 730<t<1020 1480<t<1775 2250<t<2550 3050<t<3300 3775<t<4050 4480<t<4770
Histogram 10Percent SeSageMix 0 t 300Sec 144keVLine PreIrrAnalysis.png 10Percent SeSageMix 730 t 1020Sec 144keVLine PreIrrAnalysis.png 10Percent SeSageMix 1480 t 1775Sec 144keVLine PreIrrAnalysis.png 10Percent SeSageMix 2250 t 2550Sec 144keVLine PreIrrAnalysis.png 10Percent SeSageMix 3050 t 3300Sec 144keVLine PreIrrAnalysis.png 10Percent SeSageMix 3775 t 4050Sec 144keVLine PreIrrAnalysis.png 10Percent SeSageMix 4480 t 4770Sec 144keVLine PreIrrAnalysis.png
Signal in Thin Window 3.424 [math] \times 10^4 \pm [/math] 185.04 2.535 [math] \times 10^4 \pm [/math] 159.22 2.037 [math] \times 10^4 \pm [/math] 142.72 1.695 [math] \times 10^4 \pm [/math] 130.19 1.14 [math] \times 10^4 \pm [/math] 106.77 1.154 [math] \times 10^4 \pm [/math] 107.42 1.115 [math] \times 10^4 \pm [/math] 105.59
Integrated Background 1.9628 [math] \times 10^4 \pm [/math] 167.4 1.4035 [math] \times 10^4 \pm [/math] 139.3 1.1473 [math] \times 10^4 \pm [/math] 128.80 9877 [math] \pm [/math] 118.3 7133 [math] \pm [/math] 86.1 6951.7 [math] \pm [/math] 84.0 7182 [math] \pm [/math] 85.4
Signal - Background 1.4612 [math] \times 10^4 \pm [/math] 249.52 1.1315 [math] \times 10^4 \pm [/math] 211.55 8897 [math] \pm [/math] 192.24 7073 [math] \pm [/math] 175.91 4267 [math] \pm [/math] 137.16 4588.3 [math] \pm [/math] 136.36 3968 [math] \pm [/math] 135.80
Runtime (s) 300 290 295 300 250 275 290
Rate (Hz) 48.71 [math]\pm[/math] 0.83 39.02 [math] \pm [/math] 0.73 30.16 [math] \pm [/math] 0.65 23.58 [math] \pm [/math] 0.59 17.07 [math] \pm [/math] 0.55 16.68 [math] \pm [/math] 0.49 13.68 [math] \pm [/math] 0.47

Without any corrections, the half life here is 38.61 [math] \pm [/math] 0.84 minutes. This is likely Cl-34m (32 min half life) produced by a (g,n) reaction on Cl-35 (75.76% natural abundance), which is stable, because a second line characteristic of Cl-34m was found (1176 keV). This could also be produced by a (g,pn) knockout on Ar-36 in air. This is chlorine since using the isotope search with an energy uncertainty of 5 keV and a half life uncertainty of 10 minutes, the only isotope that shared 146 and 1176 keV with branching ratios that made sense was Cl-34m.


9/10/18 Irradiation

The mass information is

Al Cylinder: 26.0249g

Sage Ash (oven method): 2.0913g

Beam off: 18:30

Calibration

Below is the runlist used for the calibration on Detector D


Source Serial No. Reference Date Activiy Half life Expected Energies (keV) Line Seen (Channel No.) Runtime (s) Position File Name Date Measured
Mn-54 J4-348 8/1/12 9.882 uCi 312.20 days 834.827 835.117 300.601 10 LB_Mn_54_9_10_18_Pos10_DetD 9/10/18
Cs-137 129793 7/1/08 1.066 uCi 30.08 years 661.657 661.887 141.356 10 LB_Cs137_9_10_18_Pos010_DetD 9/10/18
Co-60 129739 07/01/08 1.082 uCi 1925.28 days 1173.228, 1332.492 1173.338, 1332.564 232.521 10 LB_Co60_9_10_18_Pos010_DetD 9/10/18
Ba-133 129790 07/01/08 1.188 uCi 10.551 years 80.9979, 276.3989, 302.8508, 356.0129, 383.8485 81.391, 276.741, 303.166, 356.315, 384.137 250.816 10 LB_Ba133_9_10_18_Pos010_DetD 9/10/18

The calibration coefficients are (according to my fits)

Intercept = -0.370972 +/- 0.616

Slope = 0.490352 +/- 0.000534241

The calibration used at the IAC was

Intercept = 0.272366

Slope = 0.490198

Efficiency

Below is a table of the runlist used for the efficiency on Detector D.

Source Serial No. Reference Date Activiy Half life Expected Energies (keV) Start Stop Runtime (s) Position File Name Date Measured
Mn-54 J4-348 8/1/12 9.882 uCi 312.20 days 828.859 12:45 12:50 300.601 10 LB_Mn_54_9_10_18_Pos10_DetD 9/10/18
Mn-54 J4-348 8/1/12 9.882 uCi 312.20 days 828.859 12:52 12:57 305.787 20 LB_Mn_54_9_10_Pos020_DetD 9/10/18
Mn-54 J4-348 8/1/12 9.882 uCi 312.20 days 828.859 12:59:41 113:04 300.734 30 LB_Mn_54_9_10_18_Pos030_DetD 9/10/18
Co-60 129740 07/01/08 10.42 uCi 1925.28 days 1173.228, 1332.492 18:17:35 18:24 409.011 80 LB_Co60_9_10_18_Pos080_DetD 9/10/18
Co-60 129740 07/01/08 10.42 uCi 1925.28 days 1173.228, 1332.492 18:26:11 18:29 201.084 70 LB_Co60_9_10_18_Pos070_DetD 9/10/18

Runlist for 9/10/18 Irradiated Samples

Sample Date Start Stop Runtime (s) Position Detector % Dead File Name
Sage Ash 9/10/18 20:31:35 09:26 46507.875 40 Detector D < 2% LB_SageAsh_OvenMethod_9_10_18_DetD_002
Sage Ash 9/11/18 14:59:50 10:57 71854.534 40 Detector D <1% LB_SageAsh_OvenMethod_9_11_18_DetD_003
Sage Ash 9/12/18 11:01:05 20 Detector D <1% LB_SageAsh_OvenMethod_9_12_18_DetD_Pos20_001

Sage Background Analysis

Below is the first measured histogram generated from the file LB_SageAsh_OvenMethod_9_10_18_DetD_002. Note that the Al cylinder and the sage ash were not separated, so there will be some lines from Al present. The histogram shown below is not weighted because there are two different samples in front of the detector at the same time (Al cylinder and sage ash). Once the lines have been identified, the mass will be taken into account

The command lines used in the directory /data/IAC/LB_Thesis_Analysis/Sept2018_SageAshBackground/Samples

.L Eff.C

Eff t

t.Loop(RangeMin,RangeMax);


LB SageAshOven 9 10 18 FirstMeasure.png


Below is a histogram of the second measurement taken overnight on 9/11/18

LB SageAshOven 9 11 18 SecondMeasure.png

Now let's overlay these two histograms as a function of activity rather than counts to see if there are any lines that have decayed away that could possibly be identified. I did this as Hz vs. 0.5*(Chan)-0.4. The dithering proved to be quite difficult for both histograms.

LB SageAshOven Measure1and2 Overlay.png

Below is a table detailing the decaying/decayed lines. The relative rates were calculated using a linear fit plus a Gaussian fit, then the background subtracted activity was found.


Histogram Line Channel Energy (keV) Relative Rate ([math] \frac{A_2}{A_1} [/math])
863 1726.8 846 0
1846 3692 1810 0
2155 4310.8 2133.4 0

846 keV Line Analysis

Below is the analysis table. The energy window under study is [841,851]


0<t<7751.3125 7751.3125 <t< 15502.625 15502.325 <t< 23253.9375
Histogram PureSageAsh 1stMeasurement 0 t 7751Sec 846keVLine.png PureSageAsh 1stMeasurement 7751 t 15502Sec 846keVLine.png
Signal in Window 1.013 [math] \times 10^5 \pm [/math] 318.28 5.856 [math] \times 10^4 \pm [/math] 241.99
Integrated Background 6040 [math] \pm [/math] 110 4912 [math] \pm [/math] 98
Signal - Background 9.5260 [math] \times 10 ^4 \pm [/math] 336.75 5.3648 [math] \times 10^4 \pm [/math] 261.08
Runtime (s) 7751.3125 7751.3125 7751.3125
Rate (Hz) 12.29 [math] \pm [/math] 0.04 6.92 [math] \pm [/math] 0.03