Difference between revisions of "JB Absolute theta"

From New IAC Wiki
Jump to navigation Jump to search
 
(10 intermediate revisions by the same user not shown)
Line 1: Line 1:
 
[[Production Analysis|go_back]]
 
[[Production Analysis|go_back]]
  
= Overview =
+
[[JB Absolute theta(oldmethod)]]
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 =
+
=Analysis=
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 .
+
After my first attempt at reconstructing  the theta_abs distribution gave mediocre results, I decided to try again, but this time without integrating over experimental observables throught analysis, except for at the final step. To illustrate what I mean, see the histogram below:
  
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.  
+
[[File:D2O3DHistogram.png|700px|thumb|center|alt=Large | Every singles event lies in a three dimensional space consisting of a PMT top and bottom time, and a specific detector. These observables can be recast as energy, vertical z position, and detector angle (54,78,102, ect). ]]
  
[[File:MCNPSimD20Theta abs.png|650px]]
+
The benefit of conducting the analysis in this 3D space, is that the neutrons correlated with the beam can be normalized to uncorrelated Cf252 neutrons detected in the same detector with similar positions and energies.
 
 
= 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.
 
 
 
[[File:NCorrTheta abs histogram.png||350px]]
 
 
 
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.
 
 
 
Below is a step by step outline of the of analysis procedure. All rates are in terms of event per pulse.
 
# 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 reasonable approximation since the gamma rate fluctuations are small. This assumption was tested by calculating the average (weighted by counts) of dead-time corrections measured over 10 minute intervals. The result differed by <1%.
 
# Subtract from the target distribution, a cosmics-subtracted, dead-time corrected Al distribution. This will eliminate noise that results strictly from the beam
 
## scale the Al distribution in each detector individually so that its gamma rate is equal to that of the target distribution.
 
## Scale the cosmics distribution with the same scaling factors.
 
## Correct Al for dead-time in the same manor described in step 1.
 
## Subtract the scaled cosmics distribution from the scaled and dead-time corrected Al distribution. The resulting distribution should now contain only beam-related noise. 
 
#Subtract cosmics distribution from deat-time corrected target distribution.
 
# Divide by Cf252 distribution to get the relative rate over theta_abs, in each detector.
 
 
 
Note: Cosmic noise from Cf252 distribution was less than 0.1% effect.
 
 
 
==non-normalized D2O and Cf252 distributions==
 
 
 
[[File:D2OTheta absAnalysisSteps(1).png|650px|thumb|center|alt=Large | Thet_abs distribution of neutron singles from D2O. No corrections applied. ]]
 
 
 
[[File:D2OTheta absAnalysisSteps(2).png|650px|thumb|center|alt=Large | Cf252 neutron distribution. ]]
 
 
 
==Normalized distributions==
 
[[File:D2OTheta absAnalysisSteps(3).png |650px|thumb|center|alt=Large |  D2O dist. normalized to Cf252. No corrections applied to D2O distribution.]]
 
 
 
[[File:D2OTheta absAnalysisSteps(4).png |650px|thumb|center|alt=Large | Normalized D2O distribution after applying dead-time-correction. Correction for beam right:  1.816, 1.068, 1.011, 1.01, 1.007. Correction for beam left: 1.592, 1.079, 1.015, 1.01, 1.012]]
 
 
 
[[File:D2OTheta absAnalysisSteps(5).png|650px|thumb|center|alt=Large | Dead-time corrected cosmic/noise subtracted D2O. Cosmic and Al subtractions produced a change of ~3%]]
 
 
 
==Final results ==
 
[[File:D2OTheta absAnalysisSteps(6).png|650px]]
 
 
 
[[File:ThNeutronsNormalized2Cf252.png|650px]] [[File:ThNeutronsNormalized2Cf252LRAverage.png|650px]]
 
 
 
[[File:DUNeutronsNormalized2Cf252.png|650px]]  [[File:DUNeutronsNormalized2Cf252LRAverage.png|650px]]
 
 
 
= Idiot checks =
 
[[File:D2OGammasNormalized2Cf252.png |650px|thumb|center|alt=Large | Same analysis applied to Gammas. ]]
 
 
 
[[File:D2ONoiseNormalized2Cf252.png |650px|thumb|center|alt=Large |Same analysis applied to noise. ]]
 
 
 
[[File:D2O Cf252NeutronsLeftRightRatio.png |650px|thumb|center|alt=Large|  ]]
 

Latest revision as of 21:28, 22 November 2017

go_back

JB Absolute theta(oldmethod)

Analysis

After my first attempt at reconstructing the theta_abs distribution gave mediocre results, I decided to try again, but this time without integrating over experimental observables throught analysis, except for at the final step. To illustrate what I mean, see the histogram below:

Large
Every singles event lies in a three dimensional space consisting of a PMT top and bottom time, and a specific detector. These observables can be recast as energy, vertical z position, and detector angle (54,78,102, ect).

The benefit of conducting the analysis in this 3D space, is that the neutrons correlated with the beam can be normalized to uncorrelated Cf252 neutrons detected in the same detector with similar positions and energies.