Difference between revisions of "A W thesis old"

From New IAC Wiki
Jump to navigation Jump to search
Line 190: Line 190:
 
|-
 
|-
 
|}
 
|}
 +
=Fit method=
 +
 +
 +
Do the fit:
 +
 +
get parameters for line
 +
 +
Then fill 1 histogram with line
 +
 +
Then subtract
 +
 +
<pre>
 +
 +
TH1F *coin1=new TH1F("coin1","coin1",30,1800,1860);
 +
 +
ntuple->Draw("ADC7*0.604963-49.7001 >>coin1")
 +
 +
TH1F *lin1=new TH1F("lin1","line1",30,1800,1860);
 +
 +
for(int i=1800;i<1861;i++){
 +
lin1->Fill(i,-2028+1.12*i)
 +
}
 +
 +
TH1F *sub1=new TH1F("sub1","sub1",30,1800,1860);
 +
 +
sub1->Add(coin1,1);
 +
sub1->Add(lin1,-1);
 +
 +
sub1->Draw();
 +
 +
</pre>
 +
 +
 +
A pdf of the Mathematica notebook used to calculate background area, gaussian area, and plot signal/noise vs. activity.
 +
 +
[[File:AW_Background_Noise_custom4.pdf]]
 +
 +
All the ROOT fit parameters used to find the background and the resulting peaks.
 +
 +
[[File:Y-88_Fit_Log_Scaled&Cut.pdf]]
 +
 +
Integrating the gaussian of the HpGe detector signal.
 +
 +
[[File:AW_Gaussian_Integral2.pdf]]
 +
 +
 +
 +
 +
 +
:<math> f(x)=A \textstyle \int_{\mu-2\sigma}^{\mu+2\sigma} e^{-\frac{1}{2}\left(\frac{x-\mu}{\sigma}\right)^2} \operatorname{d}\!x +B x +C</math>
 +
  
  
 
[[A_W_thesis]]
 
[[A_W_thesis]]

Revision as of 20:10, 29 August 2014

8 keV signal table

[math]\sigma_Cut[/math]

[math]T_{1/2}[/math] Trig Signal BackG Subtracted Fit Parameters Signal Area Noise Area SNR
0.44 coin Run6980 898 raw.png Run6980 898 sub.png [math]\mu=[/math] 895.744 +/- 5.774e-4 [math]\sigma=[/math] 0.8189 +/- 0.0034 A= 0.0097 +/- 0.0013 B= 4.891e-9 +/- 6.707e-9 C= 2.543e-5 +/- 1.053e-4 0.0217 +/- 6.383e-4 4.643e-4 +/- 5.804e-5 46.678 +/- 0.0124
0.98 sing Run7022 898 raw.png Run7022 898 sub.png [math]\mu=[/math] 897.255 +/- 0.0047 [math]\sigma=[/math] 0.7277 +/- 0.0280 A=0.7027 +/- 0.0129 B= 8.448e-4+/- 4.457e-5 C=-0.7355 +/- 0.0408 1.2938 +/- 0.0381 0.3123 +/- 0.0390 4.143 +/- 0.6977
coin Run7023 898 raw1.png Run7023 898 sub1.png [math]\mu=[/math] 896.818 +/- 0.0075 [math]\sigma=[/math] 0.9856 +/- 0.0276 A= 0.0020 +/- 2.036e-5 B= 1.263e-8 +/- 7.797e-10 C= 1.095e-5 +/- 5.103e-7 0.0049 +/- 1.446e-4 9.518e-5 +/- 1.190e-5 51.591 +/- 0.0036
1.46 sing Run7107 898 raw.png Run7107 898 sub.png [math]\mu=[/math] 897.663 +/- 0.00987 [math]\sigma=[/math] 0.8348 +/- 0.0580 A=0.3925 +/- 0.1107 B=-2.696e-4 +/- 2.725e-7 C= 0.2633 +/- 1.939e-4 1.075 +/- 0.0316 0.0777 +/- 0.0097 13.827 +/- 0.3763
coin Run7108 898 raw.png Run7108 898 sub.png [math]\mu=[/math] 897.780 +/- 0.0025 [math]\sigma=[/math] 0.9737 +/- 0.0186 A= 1.679e-3 +/- 5.996e-7 B= 1.217e-8 +/- 7.211e-11 C= 1.098e-5 +/- 2.962e-8 3.987e-3 +/- 1.174e-4 8.732e-5 +/- 1.092e-5 45.661 +/- 0.0030
2.12 sing Run7203 898 raw.png Run7203 898 sub.png [math]\mu=[/math] 894.701 +/- 0.0127 [math]\sigma=[/math] 0.6371 +/- 0.0569 A=0.3740 +/- 6.200e-5 B= 8.285e-5+/- 4.503e-8 C=-0.0482+/- 5.103e-7 0.6212 +/- 0.0183 0.0723 +/- 0.0090 8.588 +/- 0.3858
coin Run7204 898 raw.png Run7204 898 sub.png [math]\mu=[/math] 890.896 +/- 0.0144 [math]\sigma=[/math] 1.137 +/- 0.0056 A= 7.293e-4 +/- 8.083e-9 B= 8.441e-7 +/- 6.284e-11 C= -7.342e-4 +/- 2.106e-7 0.00199 +/- 5.872e-5 8.296e-5 +/- 1.037e-5 24.032 +/- 0.0015
2.27 sing Run7235 898 raw2.png Run7235 898 sub2.png [math]\mu=[/math] 903.478 +/- 0.0064 [math]\sigma=[/math] 0.7121 +/-0.0877 A=0.2929 +/- 2.373e-4 B= -7.136e-4+/- 1.414e-9 C= 0.6715 +/- 1.777e-4 0.5609+/- 0.0165 0.0862 +/- 0.0108 6.505 +/- 0.3063
coin Run7236 898 raw.png Run7236 898 sub.png [math]\mu=[/math] 905.932 +/- 0.0050 [math]\sigma=[/math] 0.8066 +/- 0.0236 A= 5.888e-4 +/- 2.063e-7 B= 1.238e-6 +/- 4.497e-10 C= -1.106e-3 +/- 4.701e-7 1.170e-3 +/- 3.446e-5 5.310e-5 +/- 6.638e-6 22.033 +/- 8.137e-4


Re-Fit of Run 7203 without Event Cuts

[math]T_{1/2}[/math] Trig Signal BackG Subtracted Fit Parameters Signal Area Noise Area SNR
2.12 sing Run7203 898 raw nocut.png Run7203 898 sub nocut.png [math]\mu=[/math] 894.859 +/- 0.0377 [math]\sigma=[/math] 0.8294 +/- 0.0389 A= 0.1064 +/- 1.281e-4 B=4.131e-4 +/- 3e-9 C= -0.3439 +/- 1.047e-4 0.2214 +/- 0.0065 0.0900 +/- 0.0113 2.459 +/- 0.0082


Re-Fit of Run 7203 without Event Cuts and Scaled by 1.66

[math]T_{1/2}[/math] Trig Signal BackG Subtracted Fit Parameters Signal Area Noise Area SNR
2.12 sing Run7203 898 raw ratio3.png Run7203 898 sub ratio3.png [math]\mu=[/math] 894.901 +/- 0.0878 [math]\sigma=[/math] 0.8380 +/- 0.0389 A= 0.1774 +/- 4.505e-4 B=4.931e-5 +/- 9.842e-7 C= -2.231e-3 +/- 9.027e-4 0.3731 +/- 0.0110 0.1532 +/- 0.0192 2.435 +/- 0.1543



Signal, Noise, and Signal to Noise Ratio plots for 898 keV singles runs.

S N SNR 898 singles2.png SNR 898 singles2.png

Averaging over the last three noise values:

S N SNR 898 singles avg2.png SNR 898 singles avg2.png

Signal, Noise, and Signal to Noise Ratio plots for 898 keV coincidence runs.


S N SNR 898 coin.png SNR 898 coin.png

Averaging over the noise values:

S N SNR 898 coin avg.png SNR 898 coin avg.png

Signal to Noise Ratio (Singles and Coincidence): SNR 898 both.png

SNR 898 both2.png


=other

Run7204 Y88 898 signalnoise.png Run7204 Y88 898 signal.png

Run7235 Y88 898 signalnoise.png Run7235 Y88 898 signal.png

Run7236 Y88 898 signalnoise.png Run7236 Y88 898 signal.png

898 keV Signal and Background Noise Table

Run # Area of Signal Area of Background Signal to Noise Ratio
7203 130.17 120.59 1.08
7204 825.34 558.66 1.48
7235 194.38 60.23 3.23
7236 1429.75 436.93 3.27

A_W_thesis

Run7204 Y88 1836 signalnoise.png Run7204 Y88 1836 signal.png

Run7235 Y88 1836 signalnoise.png Run7235 Y88 1836 signal.png

Run7236 Y88 1836 signalnoise.png Run7236 Y88 1836 signal.png

1836.1 keV Signal and Background Noise Table

Run # Area of Signal Area of Background Signal to Noise Ratio
7203 83.29 24.73 3.37
7204 944.58 101.68 9.29
7235 82.25 20.90 3.94
7236 1849.86 97.64 18.95


Efficiency on 2/28/2014

Run Source Energy Expected rate (Hz) HpGe Det B Rate (ADC 7) (Hz) HpGe Eff (%)
7028 Ba-133 80.9971 1067.468 +/- 1.131 14.483 +/- 1.442 1.357
7028 Ba-133 276.3997 221.297 +/- 0.2345 4.239 +/- 0.9671 1.916
7028 Ba-133 302.8510 574.310 +/- 0.6085 9.523 +/- 0.5539 1.658
7028 Ba-133 356.0134 1941.418 +/- 2.057 27.130 +/- 1.316 1.397
7028 Ba-133 383.8480 278.416 +/- 0.2950 3.982 +/- 0.2440 1.430
7027 Na-22 511 1589.218 +/- 1.842 17.545 +/- 0.4560 1.104
7025 Cs-137 661.657 2548.022 +/- 2.222 25.622+/- 0.3286 1.006
7029 Mn-54 834.848 4032.782 +/- 0.0114 0.2911 +/- 0.1667 0.0072
7026 Co-60 1173.228 1876.619 +/- 0.2687 11.062 +/- 0.1215 0.589
7027 Na-22 1274.537 888.264 +/- 1.029 4.685 +/- 0.0825 0.527
7026 Co-60 1332.492 1878.167 +/- 0.2690 9.616 +/- 0.1181 0.512


Signal, Noise, and Signal to Noise Ratio plots for 898 keV singles runs.

SNR 898 singles plot.png Ratio 898 singles plot.png

Signal, Noise, and Signal to Noise Ratio plots for 898 keV coincidence runs.

SNR 898 coin plot.png Ratio 898 coin plot.png

Signal Diagram.png Signal Diagram 2.png


[math]\sigma_Fit[/math]

[math]T_{1/2}[/math] Trig Fit Parameters Signal Area Noise Area SNR
0.98 sing [math]\mu=[/math] 897.260 +/- 0.0047 [math]\sigma=[/math] 0.7074 +/- 0.0280 A=0.725 +/- 0.0129 B= 8.448e-4+/- 4.457e-5 C=-0.7355 +/- 0.0408 0.9938 +/- 0.0054 0.3039 +/- 0.0380 3.270 +/- 0.0243
coin [math]\mu=[/math] +/- [math]\sigma=[/math] +/- A= +/- B= +/- C= +/- +/- +/- +/-
2.27 sing [math]\mu=[/math] 903.482 +/- 0.0064 [math]\sigma=[/math] 0.6595 +/- 0.0877 A= 0.3199 +/- 0.0156 B= -7.248e-4+/- 9.968e-7 C= 6.816e-1 +/- 8.138e-4 0.5146 +/- 0.0073 0.0852 +/- 0.0107 6.037 +/- 0.0082
coin [math]\mu=[/math] +/- [math]\sigma=[/math] +/- A= +/- B= +/- C= +/- +/- +/- +/-

Fit method

Do the fit:

get parameters for line

Then fill 1 histogram with line

Then subtract 


TH1F *coin1=new TH1F("coin1","coin1",30,1800,1860);

ntuple->Draw("ADC7*0.604963-49.7001 >>coin1")

TH1F *lin1=new TH1F("lin1","line1",30,1800,1860);

for(int i=1800;i<1861;i++){
lin1->Fill(i,-2028+1.12*i)
}

TH1F *sub1=new TH1F("sub1","sub1",30,1800,1860);

sub1->Add(coin1,1);
sub1->Add(lin1,-1);

sub1->Draw();


A pdf of the Mathematica notebook used to calculate background area, gaussian area, and plot signal/noise vs. activity.

File:AW Background Noise custom4.pdf

All the ROOT fit parameters used to find the background and the resulting peaks.

File:Y-88 Fit Log Scaled&Cut.pdf

Integrating the gaussian of the HpGe detector signal.

File:AW Gaussian Integral2.pdf



[math] f(x)=A \textstyle \int_{\mu-2\sigma}^{\mu+2\sigma} e^{-\frac{1}{2}\left(\frac{x-\mu}{\sigma}\right)^2} \operatorname{d}\!x +B x +C[/math]


A_W_thesis

Retrieved from "https://wiki.iac.isu.edu/index.php?title=A_W_thesis_old&oldid=94110"