Difference between revisions of "2n Position resolution"

From New IAC Wiki
Jump to navigation Jump to search
Line 122: Line 122:
 
| [[File:resp_resol_DetK.png|400px|Region 3]]
 
| [[File:resp_resol_DetK.png|400px|Region 3]]
 
|- align="center"
 
|- align="center"
| Experiment|| Simulation w/ y-resolution; Solid angle = 0.093 sr||Simulation w/o y-resolution; Solid angle = 0.092 sr
+
| Fit:  cm/ns; Exp:  cm/ns|| Fit: cm/ns; Exp:  cm/ns||Fit: cm/ns; Exp: cm/ns
 
|}
 
|}
  
  
 
[[File:resp_resol_DetI.png | 400 px]]
 
[[File:resp_resol_DetI.png | 400 px]]

Revision as of 15:08, 10 May 2013

Discrete FFT algorithm test

Convolution deconvolution stepFcn.png

So as can be seen, deconvolution of the final curve (blue) with the resolution function (red) will give the response step function (black) which has the same shape as the initial step function.

If we take narrower response function we will get

Convolution deconvolution stepFcn2.png

FFT application

In this section it was used run 4133 data with DU target.

Det G

Experimental distribution of TDC difference time spectrum (black dots) and gaussian fit of the data (red line) are shown below:

Experimental data DetG gauss.png

The fit of the experimental data with a function which is the result of convolution of the detector resolution function and detector response function is shown below

Experimental data DetG fftFit.png

Deconv data fit.png

As can be seen the gaussian fit parameter "w" is equal to the parameter "Sigma3" of my_fcn fit of the experimental data and also "xc" is equal to "M3". Hence fitting the experimental data with my_fcn one can get the position resolution of the detector.

Comparison of the response function obtained as a result of the deconvolution of my_fit function and gaussian with the parameters from the fit and response function drawn using the parameters from the fit ("A3" and "B3").

Response fcn comparison.png

The width of the ideal response function is 11 ns which corresponds to the result of the calibration with Co-60 source: for Det G (current case) it was obtained the conversion factor 7.24 cm/ns and, as the detector active area width is 75 cm, we obtain 75/7.24 = 10.36 ns width of the ideal response function.

Hence, any coordinate of the neutron hit has the following uncertainty in: sigma = 2.129 ns and its value converted into cm is 2.129 ns x 7.24 cm/ns = 15.4 cm. The total length of the active area of the neutron detector is 75 cm.

If we cut off all the photons from the distribution we will get:

Neutron data DetG fftFit2.png

In this case the position resolution is 1.59 ns x 7.24 cm/ns = 11.5 cm

Det F

Fit of the experimental data in the case when photons and neutrons are combined:

Experimental data DetF fftFit.png

Fit of the experimental data in the case when only neutrons are included:

Neutrons data DetF fftFit.png

In this case the position resolution is 1.47 ns x 7.37 cm/ns = 10.8 cm

Det H

Fit of the experimental data in the case when only neutrons are included:

Neutrons data DetH fftFit.png

In this case the position resolution is 1.4 ns x 7.75 cm/ns = 10.85 cm

Det M

Fit of the experimental data in the case when only neutrons are included:

Neutrons data DetM fftFit.png

In this case the position resolution is 1.54 ns x 6.82 cm/ns = 10.5 cm

Det E

Fit of the experimental data in the case when only neutrons are included:

Neutrons data DetE fftFit2.png

In this case the position resolution is 0.8 ns x 7.41 cm/ns = 5.9 cm

Det K

Fit of the experimental data in the case when only neutrons are included:

Neutrons data DetK fftFit2.png

In this case the position resolution is 1.4 ns x 6.71 cm/ns = 9.4 cm

Det I

Fit of the experimental data in the case when only neutrons are included:

Neutrons data DetI fftFit2.png

D2O data analysis

Position resolution1.png

Position resolution2.png

Response fcns and resolution fcns reconstructed

Reconstructed Reconstructed Reconstructed
Region 1 Region 2 Region 3
Fit: 6.7 cm/ns; Exp: 7.41 cm/ns Fit: 5.7 cm/ns; Exp: 6.82 cm/ns Fit: 5.8 cm/ns; Exp: 7.37 cm/ns


Reconstructed Reconstructed Reconstructed
Region 1 Region 2 Region 3
Fit: cm/ns; Exp: cm/ns Fit: cm/ns; Exp: cm/ns Fit: cm/ns; Exp: cm/ns


Resp resol DetI.png