Difference between revisions of "Theoretical analysis of 2n accidentals rates"

From New IAC Wiki
Jump to navigation Jump to search
Line 6: Line 6:
  
 
==Probability of detecting a given pair of neutrons in a single pulse==
 
==Probability of detecting a given pair of neutrons in a single pulse==
 
 
=Old, needed to be simplified=
 
Consider a pulse that causes three neutron producing reactions, two 1n-knockouts and a photofission event with multiplicity of 4. In terms of the notation, <math>N=3</math>, <math>V_{1}=1</math>, <math>V_{2}=1</math>, and <math>V_{3}=4</math>. Now let's find the probability, <math>p_{acc}</math>, of detecting only and exactly the following pair of neutrons:
 
# the neutron from the 1n-knockout, <math>V_{2}</math>, and
 
# A given neutron from the photofission event, <math>V_{3}</math>, of which there are four neutrons to choose from, but here we consider the detection of a single particular neutron out of the four.
 
This example is an accidental, since each neutron is created in a separate reaction. Let <math>D1</math> represent the event where the neutron from the 1n-knockout in <math>V_{2}</math> is detected, and <math>D2</math> for the detection of the photofission neutron. Define <math>\overline{D_{others}}</math> as the event that all other neutrons are ''not'' detected. Using standard notation, the probability can be written:
 
 
::<math>p_{acc} = P\left(D_1,D_2,\overline{D_{others}},N=3,V_1=1,V_2=1,V_3=4\right)</math>
 
 
,which is interpreted as the Boolean product (AND) between all events separated by the commas, which are used to enhance readability.
 
 
By an application of the general product rule of probability (also called the chain rule), <math>p_{acc}</math> can be factored as follows:
 
 
::<math>p_{acc} = P\left(D_1,D_2,\overline{D_{others}}\big| N=3,V_1=1,V_2=1,V_3=4\right)\times P(N=3,V_1=1,V_2=1,V_3=4)</math>
 
 
where the factor on the left is a conditional probability which may be read:
 
:"The probability of the detection of both and only the two neutrons under consideration, given that two 1n-knochouts and a single photofission event, with a multiplicity of four, occurred during this pulse."
 
The factor on the right can be simplified into four factors by recognizing that the events are all mutually independent, since each photo-nuclear interaction is assumed to be uncorrelated. This gives
 
 
::<math>p_{acc} = P\left(D_1,D_2,\overline{D_{others}}\big| N=3,V_1=1,V_2=1,V_3=4\right)\times P(N=3)P(V_1=1)P(V_2=1)P(V_3=4)</math>
 
 
In order to save space, the event <math>(N=3,V_1=1,V_2=1,V_3=4)</math> will hereafter be denoted as <math>E^{3}_{1,1,4}</math>.
 
 
Now to tackle the factor on the left, <math>P(D_1,D_2,\overline{D_{others}}|E^{3}_{1,1,4})</math>. Application of the product rule gives:
 
 
::<math>P(D_1,D_2,\overline{D_{others}}\big|E^{3}_{1,1,4}) =P(D_1,D_2 \big|,\overline{D_{others}},E^{3}_{1,1,4})P(\overline{D_{others}}) </math>
 
 
With ten detectors in total, each covering a solid angle of 0.5% of <math>4\pi</math>, and having an intrinsic efficiency of less than 25%, <math>P(\overline{D_{others}})</math> may be approximated by 1, or quantitatively: 
 
::<math>P(\overline{D_{others}})\approx (1-8*0.005*0.20)^{n_{total}}>0.96 \qquad  \text{for } n_{total}<7 </math>
 
::thus,
 
::<math>P(\overline{D_{others}})\approx 1</math>
 
:where,
 
:<math>8</math> is the number of detectors available.
 
:<math>n_{total}</math> is the total number neutrons emitted.
 
 
Seven neutrons emitted in a single pulse would be an extreme case under the conditions of the experiment, and so  <math>P(\overline{D_{others}})</math> may be set to unity as a approximation.
 
 
The expression for the probability of detecting our two neutrons becomes:
 
 
::<math>p_{acc} = P(D_1,D_2,|E^{3}_{1,1,4})\times P(E^{3}_{1,1,4})</math>
 
 
Since we are looking at accidentals only, <math>D_1</math> and <math>D_2</math> are uncorrelated, and thus their only dependence on <math>E^{3}_{1,1,4}</math> is the existence of the neutrons which <math>D_1</math> and <math>D_2</math> are referring to. There's no use in picking either <math>D_1</math> or <math>D_2</math> which refer to neutrons that were never created, so the product, <math>D_1D_2</math>, may be considered independent of <math>E^{3}_{1,1,4}</math>.
 
 
This brings us to an important point in the analysis, where <math>p_{acc}</math> is evidently a product of mutually independent variables:
 
 
::<math>p_{acc} = P(E^{3}_{1,1,4})\times P(D_1D_2) = P(N=3)P(V_1=1)P(V_2=1)P(V_3=4)P(D_1D_2)</math>
 
 
==probability of detecting any pair of accidentals in a given pulse==
 
 
The above result for <math>p_{acc}</math> would be the same had <math>D_1</math> and <math>D_2</math> had been chosen to refer to different neutrons. In the specific example above, there are 9 possible pairs of accidentals, the four fission neutrons paired with each 1n-knockout neutron, plus the case of the two knockout neutrons pairs with each other. Thus, the probability of detecting ''any'' accidentall:
 
 
::<math>P(A) = 9*P(N=3)P(V_1=1)P(V_2=1)P(V_3=4)P(D_1D_2)</math>
 
 
The number of pairs of accidentals is, in general:
 
 
::<math>\sum V_iV_j \qquad \text{for every} \{V_i,V_j\} \text{that is a subset of } \{V_1,V_2,...V_N\}</math>
 

Revision as of 03:06, 20 January 2018

go_back

Introduction

A given photon pulse may cause multiple neutron-producing reactions, ranging from zero to "infinity" reactions. The number of neutron-producing reactions in a pulse is hereafter denoted by [math]N[/math]. Being the number of neutron-producing reactions actually occurring per pulse, [math]N[/math] is assumed to follow the Poissonian distribution as a limiting case of the binomial distribution. Each neutron-producing interaction is said to produce [math]V_{i}[/math] correlated neutrons, where the random variable [math]V_{i}[/math] is the distribution of the number of neutrons produced in a single neutron-producing reaction. Each of the [math]V_i\text{'}s[/math] are independent and identically distributed random variables, so the purpose of the subscript is to distinguish between several distinct neutron-producing interactions which may occur in a single pulse.

The beam has a Bremsstrahlung end point of 10.5 MeV, which energetically allows for only two possible neutron-producing interactions, 1n-knochout and photofission. Thus, [math]V_{i}[/math] is equal to the photofission neutron multiplicity plus a contribution at [math]V_{i}=1[/math] from 1n-knockout events. The analysis that follows does not need to distinguish between 1n-knockout events and photofission events that emit a single neutron. In both cases, a single neutron is emitted and is uncorrelated with all and any other neutrons.

Probability of detecting a given pair of neutrons in a single pulse