Difference between revisions of "AW Singles Prob Y88 2015"

From New IAC Wiki
Jump to navigation Jump to search
 
(80 intermediate revisions by 2 users not shown)
Line 6: Line 6:
  
  
 +
<math>N\left (\frac{1836}{898}\right ) =\frac{\frac{ UN(1836)}{.31}}{ \frac{UN(898)}{.78}}=\frac{5}{2} \frac{UN(1836)}{ UN(898)}=1 =</math> Ratio
  
  
 +
UCR<math>\equiv \frac{UN(1836)}{ UN(898)}=\frac{2}{5} = 0.4=</math>ratio before efficiency correction
  
 
{| border="1"  |cellpadding="20" cellspacing="0  
 
{| border="1"  |cellpadding="20" cellspacing="0  
 
|-
 
|-
| Days ||<math>T_{1/2}</math> ||  (898 keV) Integral Raw||  (1836 keV) Integral Raw||  Ratio || (898 keV) Integral Eff Cor||  (1836 keV) Integral Eff Cor || Ratio
+
| Days ||<math>T_{1/2}</math> ||  (898 keV) Sub||  (1836 keV)Sub||  UCR || (898 keV) Integral Eff Cor||  (1836 keV) Integral Eff Cor || Ratio
 
|-
 
|-
| 108 ||  1.01 ||1.364 <math>\pm 0.062 [[File:Run7022_898_raw.png]]</math>
+
| 108 ||  1.01 ||1.364 <math>\pm 0.062</math>  [[Media:Run7022_898_raw.png]] ||  <math>\pm</math> 0.048 [[Media:Run7022_1836_raw1.png]]  || 0.62 || 174.2 <math>\pm</math> 13.24 || 269.8 <math>\pm</math> 30.91 || 1.55
 
|-
 
|-
| 159 ||  1.49 ||  0.8697  <math>\pm 0.055  </math>
+
| 159 ||  1.49 ||  0.8697  <math>\pm 0.055  </math> [[Media:Run7107_898_raw.png]]  ||  <math>\pm</math> 0.042 [[Media:Run7107_1836_raw.png]]  || 0.58 || 111.07 <math>\pm</math> 9.753 || 160.8 <math>\pm</math> 20.9 || 1.45
 
|-
 
|-
|  229 || 2.15 ||  0.4733<math>\pm 0.058  </math>
+
|  229 || 2.15 ||  0.4733<math>\pm 0.058  </math> [[Media:Run7203_898_raw_ratio3.png]]  || <math>\pm</math> 0.059 [[Media:Run7203_1836_raw_nocut.png]] || 0.60 || 60.45 <math>\pm</math> 8.272 || 90.48 <math>\pm</math> 20.95 || 1.50
 
|-
 
|-
|  245 ||  2.30 || 0.5807 <math>\pm 0.049 </math>
+
|  245 ||  2.30 || 0.5807 <math>\pm 0.049 </math> [[Media:Run7235_898_raw.png]]  ||  <math>\pm</math> 0.032 [[Media:Run7235_1836_raw.png]] || 0.42 || 74.16 <math>\pm</math> 7.718 || 78.52 <math>\pm</math> 12.89 || 1.06
 
|-
 
|-
 
|}
 
|}
  
 +
 +
What is the background contribution to the above rates ( Background is all channels with rate < 0.1 Hz)
 +
 +
(Set energy window to <math> \mu \pm 3 \sigma  </math>)
 +
 +
 +
 +
{| border="1"  |cellpadding="20" cellspacing="0
 +
|-
 +
| Run # ||Days ||<math>T_{1/2}</math> || Energy (keV) || Raw Integral || Background Integral || Background (Left) || Background (Right) || Ratio (Raw/Background) || Sub Integral || Sub Back (Left) || Sub Back (Left)  || Ratio
 +
|-
 +
|  rowspan="2" | 7022 || rowspan="2" | 108 || rowspan="2" | 1.01 || 898 || 1.34  || 0.0891  || 0.2145 || 0.3148 || 15.04 || 1.252|| 0.04951 || 0.08752 
 +
|-
 +
| 1836 || 0.8524  || 0.0200 || 0.02507 || 0.01393 || 42.62  || 0.8325 ||  ||  || 0.67
 +
|-
 +
|  rowspan="2" | 7107 || rowspan="2" | 159 || rowspan="2" | 1.49 || 898 || 0.9576  || 0.1068 || 0.2332 || 0.1767 || 8.966  || 0.8501|| 0.03603 || 0.02687 
 +
|-
 +
| 1836 || 0.5300  || 0.0299  || 0.01413 || 0.02827 || 17.73  || 0.5007||  ||  || .58
 +
|-
 +
|  rowspan="2" | 7203 || rowspan="2" | 229 || rowspan="2" | 2.15 || 898 ||  0.6682 || 0.2094 || 0.3165 || 0.8652 || 3.191 ||0.4589 ||  || 
 +
|-
 +
| 1836 || 0.3539  || 0.0799 || 0.02119 || 0.7627 || 4.429  || 0.2740 ||  ||  || 0.59
 +
|-
 +
|  rowspan="2" | 7235 ||rowspan="2" | 245 || rowspan="2" | 2.30 || 898 ||  0.541  ||  0.0803|| 0.2951 || 0.2459  ||  6.737 || 0.4595|| 0.07145 || 0.0648 
 +
|-
 +
| 1836 || 0.3033  || 0.0636 || 0.01639 || 0.02869  || 4.769  || 0.2414||  ||  || 0.52
 +
|-
 +
|}
 +
 +
[[File:7022_7107_overlay.png | 200 px]]
 +
[[File:7022_7107_overlay_left.png | 200 px]]
 +
[[File:7022_7107_overlay_right.png | 200 px]]
 +
 +
 +
Weighted Averaged Background:
 +
 +
898: 0.0255 +/- 0.0206
 +
 +
1836: 0.0072 +/- 0.0120
 +
 +
{| border="1"  |cellpadding="20" cellspacing="0
 +
|-
 +
| Run # ||Days ||<math>T_{1/2}</math> || Energy (keV) || Raw Integral || Sub Integral || UCR || Amp (raw) || Amp-UCR (raw) || Amp (sub) || Amp-UCR (sub)
 +
|-
 +
|  rowspan="2" | 7022 || rowspan="2" | 108 || rowspan="2" | 1.01 || 898 || 1.34  || 1.238  || rowspan="2" | 0.65 || <math>0.7035 \pm 0.9116</math> || rowspan="2" | 0.44 || <math>0.6780 \pm 1.527</math> || rowspan="2" | <math>0.45 \pm 1.44</math>
 +
|-
 +
| 1836 || 0.8524  || 0.8092 || <math>0.3115 \pm 0.6806</math> || <math>0.3043 \pm 0.6895</math>
 +
|-
 +
|  rowspan="2" | 7107 || rowspan="2" | 159 || rowspan="2" | 1.49 || 898 || 0.9576  || 0.8301 || rowspan="2" | 0.59 ||<math> 0.3923 \pm 0.2435</math> || rowspan="2" | 0.45 || <math>0.3668 \pm 0.3737</math> || rowspan="2" | <math>0.46 \pm 2.17</math>
 +
|-
 +
| 1836 || 0.5300  || 0.4868 || <math>0.1757 \pm 0.7528</math> || <math>0.1685 \pm 0.7747</math>
 +
|-
 +
|  rowspan="2" | 7235 ||rowspan="2" | 245 || rowspan="2" | 2.30 || 898 ||  0.541  || 0.4645 || rowspan="2" | 0.53 ||<math> 0.2930 \pm 0.6434</math> || rowspan="2" | 0.25 || <math>0.2675 \pm 0.6259</math> || rowspan="2" | <math>0.25 \pm 0.46</math>
 +
|-
 +
| 1836 || 0.3033  || 0.2457 || <math>0.0740 \pm 0.0371</math> || <math>0.06677 \pm 0.1162</math>
 +
|}
 +
 +
 +
 +
 +
{| border="1"  |cellpadding="20" cellspacing="0
 +
|-
 +
| Run # ||Days ||<math>T_{1/2}</math> || Energy (keV) || Amp || Amp-UCR
 +
|-
 +
|  rowspan="2" | 7022 || rowspan="2" | 108 || rowspan="2" | 1.01 || 898 || <math>1.11421 \pm 0.057009</math> || rowspan="2" | <math>0.4425 \pm 0.0431</math>
 +
|-
 +
| 1836 || <math>0.493043 \pm 0.040830</math>
 +
|-
 +
|  rowspan="2" | 7107 || rowspan="2" | 159 || rowspan="2" | 1.49 || 898 ||<math> 0.558352 \pm 0.044758</math>  || rowspan="2" | <math>0.4746 \pm 0.0613</math>
 +
|-
 +
| 1836 || <math>0.265019 \pm 0.026870</math> 
 +
|-
 +
|  rowspan="2" | 7235 ||rowspan="2" | 245 || rowspan="2" | 2.30 || 898 || <math>0.434429 \pm 0.042097</math>  || rowspan="2" |<math> 0.3396 \pm 0.0631</math>
 +
|-
 +
| 1836 || <math>0.147541 \pm 0.023407</math>
 +
|}
 +
 +
 +
Efficiency corrected plot of the two activities as a function of days
 +
 +
[[File:integral_eff_boosted.png | 400 px]]
 +
 +
 +
==Ntuple commands==
 +
 +
;r7022 (0.002786)
 +
:TH1F *r7022=new TH1F("r7022","r7022",2200,0,4400)
 +
: ntuple->Draw("ADC7*0.599803-50.9271 >>r7022");
 +
:TH1F *N7022=new TH1F("N7022","N7022",2200,0,4400)
 +
:N7022->Add(r7022,0.002786);
 +
 +
;r7107 (0.003534)
 +
:TH1F *r7107=new TH1F("r7107","r7107",2200,0,4400)
 +
:ntuple->Draw("ADC7*0.601126-51.3834 >>r7107");
 +
:TH1F *N7107=new TH1F("N7107","N7107",2200,0,4400)
 +
:N7107->Add(r7107,0.002786);
 +
 +
;r7235 ()
 +
:TH1F *r7235=new TH1F("r7235","r7235",2200,0,4400)
 +
:ntuple->Draw("ADC7*0.61389-50.2085 >>r7235");
 +
:TH1F *N7235=new TH1F("N7235","N7235",2200,0,4400)
 +
:N7235->Add(r7235,0.004098);
  
  
 
[[A_W_thesis#Final_problems_to_sort_out]]
 
[[A_W_thesis#Final_problems_to_sort_out]]

Latest revision as of 14:51, 31 March 2015

Problem description:

The singles spectra production for the initial Y-88 activity is off by a factor of 2.
Right now the high energy line is a factor of 2 higher.


[math]N\left (\frac{1836}{898}\right ) =\frac{\frac{ UN(1836)}{.31}}{ \frac{UN(898)}{.78}}=\frac{5}{2} \frac{UN(1836)}{ UN(898)}=1 =[/math] Ratio


UCR[math]\equiv \frac{UN(1836)}{ UN(898)}=\frac{2}{5} = 0.4=[/math]ratio before efficiency correction

Days [math]T_{1/2}[/math] (898 keV) Sub (1836 keV)Sub UCR (898 keV) Integral Eff Cor (1836 keV) Integral Eff Cor Ratio
108 1.01 1.364 [math]\pm 0.062[/math] Media:Run7022_898_raw.png [math]\pm[/math] 0.048 Media:Run7022_1836_raw1.png 0.62 174.2 [math]\pm[/math] 13.24 269.8 [math]\pm[/math] 30.91 1.55
159 1.49 0.8697 [math]\pm 0.055 [/math] Media:Run7107_898_raw.png [math]\pm[/math] 0.042 Media:Run7107_1836_raw.png 0.58 111.07 [math]\pm[/math] 9.753 160.8 [math]\pm[/math] 20.9 1.45
229 2.15 0.4733[math]\pm 0.058 [/math] Media:Run7203_898_raw_ratio3.png [math]\pm[/math] 0.059 Media:Run7203_1836_raw_nocut.png 0.60 60.45 [math]\pm[/math] 8.272 90.48 [math]\pm[/math] 20.95 1.50
245 2.30 0.5807 [math]\pm 0.049 [/math] Media:Run7235_898_raw.png [math]\pm[/math] 0.032 Media:Run7235_1836_raw.png 0.42 74.16 [math]\pm[/math] 7.718 78.52 [math]\pm[/math] 12.89 1.06


What is the background contribution to the above rates ( Background is all channels with rate < 0.1 Hz)

(Set energy window to [math] \mu \pm 3 \sigma [/math])


Run # Days [math]T_{1/2}[/math] Energy (keV) Raw Integral Background Integral Background (Left) Background (Right) Ratio (Raw/Background) Sub Integral Sub Back (Left) Sub Back (Left) Ratio
7022 108 1.01 898 1.34 0.0891 0.2145 0.3148 15.04 1.252 0.04951 0.08752
1836 0.8524 0.0200 0.02507 0.01393 42.62 0.8325 0.67
7107 159 1.49 898 0.9576 0.1068 0.2332 0.1767 8.966 0.8501 0.03603 0.02687
1836 0.5300 0.0299 0.01413 0.02827 17.73 0.5007 .58
7203 229 2.15 898 0.6682 0.2094 0.3165 0.8652 3.191 0.4589
1836 0.3539 0.0799 0.02119 0.7627 4.429 0.2740 0.59
7235 245 2.30 898 0.541 0.0803 0.2951 0.2459 6.737 0.4595 0.07145 0.0648
1836 0.3033 0.0636 0.01639 0.02869 4.769 0.2414 0.52

7022 7107 overlay.png 7022 7107 overlay left.png 7022 7107 overlay right.png


Weighted Averaged Background:

898: 0.0255 +/- 0.0206

1836: 0.0072 +/- 0.0120

Run # Days [math]T_{1/2}[/math] Energy (keV) Raw Integral Sub Integral UCR Amp (raw) Amp-UCR (raw) Amp (sub) Amp-UCR (sub)
7022 108 1.01 898 1.34 1.238 0.65 [math]0.7035 \pm 0.9116[/math] 0.44 [math]0.6780 \pm 1.527[/math] [math]0.45 \pm 1.44[/math]
1836 0.8524 0.8092 [math]0.3115 \pm 0.6806[/math] [math]0.3043 \pm 0.6895[/math]
7107 159 1.49 898 0.9576 0.8301 0.59 [math] 0.3923 \pm 0.2435[/math] 0.45 [math]0.3668 \pm 0.3737[/math] [math]0.46 \pm 2.17[/math]
1836 0.5300 0.4868 [math]0.1757 \pm 0.7528[/math] [math]0.1685 \pm 0.7747[/math]
7235 245 2.30 898 0.541 0.4645 0.53 [math] 0.2930 \pm 0.6434[/math] 0.25 [math]0.2675 \pm 0.6259[/math] [math]0.25 \pm 0.46[/math]
1836 0.3033 0.2457 [math]0.0740 \pm 0.0371[/math] [math]0.06677 \pm 0.1162[/math]



Run # Days [math]T_{1/2}[/math] Energy (keV) Amp Amp-UCR
7022 108 1.01 898 [math]1.11421 \pm 0.057009[/math] [math]0.4425 \pm 0.0431[/math]
1836 [math]0.493043 \pm 0.040830[/math]
7107 159 1.49 898 [math] 0.558352 \pm 0.044758[/math] [math]0.4746 \pm 0.0613[/math]
1836 [math]0.265019 \pm 0.026870[/math]
7235 245 2.30 898 [math]0.434429 \pm 0.042097[/math] [math] 0.3396 \pm 0.0631[/math]
1836 [math]0.147541 \pm 0.023407[/math]


Efficiency corrected plot of the two activities as a function of days

Integral eff boosted.png


Ntuple commands

r7022 (0.002786)
TH1F *r7022=new TH1F("r7022","r7022",2200,0,4400)
ntuple->Draw("ADC7*0.599803-50.9271 >>r7022");
TH1F *N7022=new TH1F("N7022","N7022",2200,0,4400)
N7022->Add(r7022,0.002786);
r7107 (0.003534)
TH1F *r7107=new TH1F("r7107","r7107",2200,0,4400)
ntuple->Draw("ADC7*0.601126-51.3834 >>r7107");
TH1F *N7107=new TH1F("N7107","N7107",2200,0,4400)
N7107->Add(r7107,0.002786);
r7235 ()
TH1F *r7235=new TH1F("r7235","r7235",2200,0,4400)
ntuple->Draw("ADC7*0.61389-50.2085 >>r7235");
TH1F *N7235=new TH1F("N7235","N7235",2200,0,4400)
N7235->Add(r7235,0.004098);


A_W_thesis#Final_problems_to_sort_out