Difference between revisions of "Integrated asymmetry"

From New IAC Wiki
Jump to navigation Jump to search
Line 11: Line 11:
 
<math>Asymm^{detC} = \frac
 
<math>Asymm^{detC} = \frac
 
  {\left(N_{up}^{D_2O}/W_{up}^{D_2O}-N_{up}^{H_2O}/W_{up}^{H_2O}\right)
 
  {\left(N_{up}^{D_2O}/W_{up}^{D_2O}-N_{up}^{H_2O}/W_{up}^{H_2O}\right)
- \left({N_{sd^{D_2O}/W_{sd}^{D_2O}-N_{sd}^{H_2O}/W_{sd}^{H_2O}\right)
+
- \left(N_{sd}^{D_2O}/W_{sd}^{D_2O}-N_{sd}^{H_2O}/W_{sd}^{H_2O}\right)}
 
  {\left(\frac{N_{up}^{D_2O}}{W_{up}^{D_2O}}-\frac{N_{up}^{H_2O}}{W_{up}^{H_2O}}\right)
 
  {\left(\frac{N_{up}^{D_2O}}{W_{up}^{D_2O}}-\frac{N_{up}^{H_2O}}{W_{up}^{H_2O}}\right)
 
+ \left(\frac{N_{sd^{D_2O}}}{W_{sd}^{D_2O}}-\frac{N_{sd}^{H_2O}}{W_{sd}^{H_2O}}\right)}
 
+ \left(\frac{N_{sd^{D_2O}}}{W_{sd}^{D_2O}}-\frac{N_{sd}^{H_2O}}{W_{sd}^{H_2O}}\right)}

Revision as of 21:44, 8 June 2009

Go Back

Integrated asymmetry calculation

[math]Asymm^{detA} = \frac {\left(\frac{N_{up}^{D_2O}}{W_{up}^{D_2O}}-\frac{N_{up}^{H_2O}}{W_{up}^{H_2O}}\right) - \left(\frac{N_{sd^{D_2O}}}{W_{sd}^{D_2O}}-\frac{N_{sd}^{H_2O}}{W_{sd}^{H_2O}}\right)} {\left(\frac{N_{up}^{D_2O}}{W_{up}^{D_2O}}-\frac{N_{up}^{H_2O}}{W_{up}^{H_2O}}\right) + \left(\frac{N_{sd^{D_2O}}}{W_{sd}^{D_2O}}-\frac{N_{sd}^{H_2O}}{W_{sd}^{H_2O}}\right)} [/math]

[math]Asymm^{detC} = \frac {\left(N_{up}^{D_2O}/W_{up}^{D_2O}-N_{up}^{H_2O}/W_{up}^{H_2O}\right) - \left(N_{sd}^{D_2O}/W_{sd}^{D_2O}-N_{sd}^{H_2O}/W_{sd}^{H_2O}\right)} {\left(\frac{N_{up}^{D_2O}}{W_{up}^{D_2O}}-\frac{N_{up}^{H_2O}}{W_{up}^{H_2O}}\right) + \left(\frac{N_{sd^{D_2O}}}{W_{sd}^{D_2O}}-\frac{N_{sd}^{H_2O}}{W_{sd}^{H_2O}}\right)} [/math]

where

[math]N_{up}^{D_2O}=\sum_{ii=min}^{max}{count(ii)};\ \ \ \  integrated\ number\ of\ counts\ for\ detector\ Up,\ D_2O[/math]
[math]N_{up}^{H_2O}=\sum_{ii=min}^{max}{count(ii)};\ \ \ \  integrated\ number\ of\ counts\ for\ detector\ Up,\ H_2O[/math]
[math]N_{sd}^{D_2O}=\sum_{ii=min}^{max}{count(ii)};\ \ \ \  integrated\ number\ of\ counts\ for\ detector\ Side,\ D_2O[/math]
[math]N_{sd}^{H_2O}=\sum_{ii=min}^{max}{count(ii)};\ \ \ \  integrated\ number\ of\ counts\ for\ detector\ Side,\ H_2O[/math]
[math]W_{up}^{D_2O};\ \ \ \  weighted\ coefficient(NaI,\ Ref)\ for\ detector\ Up,\ D_2O[/math]


  • For detector A summation is over [1000:1600] bin numbers
  • For detector C summation is over [900:1600] bin numbers

Error calculation

[math]\delta (Asymm) = \sqrt{\left(\frac{\boldsymbol\partial A}{\boldsymbol\partial up}\right)^2(\boldsymbol\delta up)^2 + \left(\frac{\boldsymbol\partial A}{\boldsymbol\partial Wup}\right)^2(\boldsymbol\delta Wup)^2 + ...}[/math]

If assume:
[math]\boldsymbol\delta up = \sqrt{up};\quad\boldsymbol\delta Wup = \sqrt{Wup};\quad...[/math]

what is <math>\sqrt{up}</math> ?  
Does <math>\sqrt{up}</math> = number of counts in in the detector 
up spectrum as suggested above?  Why not use <math>N_{up}</math> as variable.  
If distribution is Binomial (detector yes/no) then width of distribution is 
<math>1.265\sqrt{N}</math> = 2 times e error.  
Relative error = <math>\frac{error}{Ave} = \frac{0.6325 \sqrt{N}}{N/2}</math>

If you are measuring total number of neutrons detected then you will have binomial if you 
break up the time spectrum into bins and don't integrate then probability distribution
is combination of detection probability and time measurement probability.

Then
[math]\delta (Asymm) = \sqrt{\left(\frac{(+)-(-)}{(+)^2}\right)^2\left[\frac{up}{Wup^2}+\frac{up^2}{Wup}+\frac{upbg}{Wupbg^2}+\frac{upbg^2}{Wupbg}\right] + \left(\frac{(+)+(-)}{(+)^2}\right)^2\left[\frac{sd}{Wsd^2}+\frac{sd^2}{Wsd}+\frac{sdbg}{Wsdbg^2}+\frac{sdbg^2}{Wsdbg}\right]}[/math]
where
[math](-)=\left[\left(\frac{up}{Wup}-\frac{upbg}{Wupbg}\right) - \left(\frac{sd}{Wsd}-\frac{sdbg}{Wsdbg}\right)\right][/math]

[math](+)=\left[\left(\frac{up}{Wup}-\frac{upbg}{Wupbg}\right) + \left(\frac{sd}{Wsd}-\frac{sdbg}{Wsdbg}\right)\right][/math]

Cases was analysed

Det A (was analized all possible combination):

D2O Up,   files# [40,56,102,108,134,205,210,230];
H2O Up, files# [44];
D2O Side, files# [48,74,78,82,86,90,94,146,180,190,225,235];
H2O Side, files# [52];

Det C (was analized all possible combination):

D2O Up,   files# [49,75,79,83,87,91,95,147,181,191,226,236];
H2O Up, files# [53];
D2O Side, files# [41,57,103,107,135,206,211,231];
H2O Side, files# [45];

Weighted coefficients was used

(File:Weight coeff.pdf)

Results

Table 1: Det A, weighted with [math]{\color{Red}NaI \ detector}[/math]
Table 2: Det C, weighted with [math]{\color{Red}NaI \ detector}[/math]
Table 3: Det A, weighted with [math]{\color{Red}Ref \ detector}[/math]
Table 4: Det C, weighted with [math]{\color{Red}Ref \ detector}[/math]


(File:Asymm table.pdf)

Change the X-axis to nanosecond or neutron energy (TF). A:Will do it.

Example of bin by bin asymmetry

Asymm RefDet DetA 40,44,48(74,78),52.jpg


Go Back