Difference between revisions of "JB Absolute theta"

From New IAC Wiki
Jump to navigation Jump to search
Line 23: Line 23:
  
 
Here is an outline of the steps of analysis:
 
Here is an outline of the steps of analysis:
# To correct the beam-on distribution for dead-time,  multiply the rate in each detector by <math>\frac{1}{1-R_{G}}</math>, where <math>R_{G}</math> is the gamma rate in the detector. This assumes the gamma rate is constant, which makes for a very good approximation since it is very nearly constant.  
+
# To correct the target distribution for dead-time,  multiply the rate in each detector by <math>\frac{1}{1-R_{G}}</math>, where <math>R_{G}</math> is the gamma rate in the detector. This assumes the gamma rate is constant, which makes for a good approximation since the gamma rate fluctuations are small.
#  
+
# Subtract from the target distribution, a cosmics-subtracted Al target distribution. This serves to eliminate beam-related noise. 
 +
## scale the Al distribution individually in each detector, such that the total number of gammas is equal to that of the target distribution.
 +
## Correct Al for dead-time in the same way described in step 1.
 +
## Scale cosmics distribution with the same scaling factors used for the Al distribution, and subtract the result from the Al distribution. This distribution should now contain only beam-related noise. 
 +
#Subtract cosmics distribution from target distribution.  
 +
# Divide by Cf252 distribution to get relative rate over theta_abs, in each detector.
 +
 
 +
Note: Cosmic noise from Cf252 distribution was less than 0.1% effect.
  
 
[[File:D2OTheta absAnalysisSteps(1).png|650px]]
 
[[File:D2OTheta absAnalysisSteps(1).png|650px]]

Revision as of 16:56, 9 November 2017

go_back

Overview

Here I use neutron single events to measure the distribution of theta_abs, or the angle between an incident photon and a resultaning photo-neutron. The distribution of uncorrelated neutrons from the SF of californium 252 is used to "divide out" the effects of detector geometry, efficiency, drifts, ect. For D2O, the result is compared to an MCNP simulation which was built to model as many aspects of the experiment as possible.

Simulation

I performed an MCNP simulation with a D2O target (axis length = 2"; dia. = 0.75") subject to a 10.5 MeV end-point bremsstrahlung beam. A mock-up of the entire neutron detector array is included in the simulation. Detector physics is modeled by applying a detection threshold in terms of light output (MeVee), which is equal to the typical MeVee produced by 0.5 MeV neutrons within the scintillator.

A Cf252 source was also simulated, allowing me to apply the exact same analysis technique to simulation and experimental data .

The plot below shows the relative distribution of neutron direction cosines w.r.t. the incident photon beam. These neutrons are not effected by scattering or detector geometry, since the direction cosine are taken right as neutrons are created.

MCNPSimD20Theta abs.png

Analysis

Ten detectors are used in the measurement, 5 beam left and 5 beam right. For each detector on beam left, there is a corresponding detector on beam right which covers the same range in theta_abs. Each left-right pair of detectors produces a peak in the theta_abs distribution, seen below.

NCorrTheta abs histogram.png

During the analysis, each peak is condensed into a single point whose value on the y-axis is the mean of the peak, and is the weighted average, or "center of mass", along the x-axis.

Here is an outline of the steps of analysis:

  1. To correct the target distribution for dead-time, multiply the rate in each detector by [math]\frac{1}{1-R_{G}}[/math], where [math]R_{G}[/math] is the gamma rate in the detector. This assumes the gamma rate is constant, which makes for a good approximation since the gamma rate fluctuations are small.
  2. Subtract from the target distribution, a cosmics-subtracted Al target distribution. This serves to eliminate beam-related noise.
    1. scale the Al distribution individually in each detector, such that the total number of gammas is equal to that of the target distribution.
    2. Correct Al for dead-time in the same way described in step 1.
    3. Scale cosmics distribution with the same scaling factors used for the Al distribution, and subtract the result from the Al distribution. This distribution should now contain only beam-related noise.
  3. Subtract cosmics distribution from target distribution.
  4. Divide by Cf252 distribution to get relative rate over theta_abs, in each detector.

Note: Cosmic noise from Cf252 distribution was less than 0.1% effect.

D2OTheta absAnalysisSteps(1).png

Final results

MCNP-POLIMI

Below is an MCNP-POLIMI simulation of a cylindrical D20 target (axis length = 2"; dia. = 0.75") subject to a bremsstrahlung photon beam with an end point of 10.5MeV. The plot below shows the relative distribution of neutron direction cosines w.r.t. the incident photon beam. All neutrons are from the photodisintegration of D20 and the direction cosine is taken as neutrons exit the target geometry.

MCNPSimD20Theta abs.png