https://docs.google.com/spreadsheets/d/1W1T2rC10UtOtrAaT8-ZAKFvX70v2PwDk_4-goYL-brI/edit?usp=sharing

Both Au samples (20mm and 22mm) are single mylar.

Both Si30 (22mm) and natSi (20mm) are double mylar.

Empty is single mylar ONLY

Dummy is double mylar + glue 

Isoflex certificates for enriched silicon samples

New 22mm sample compared to old 20mm sample used at the beginning of the campaign + picture of dummy target in the setup.

Au and Si29 runs normalised by proton number show consistency (dedicated pulses, example detector 4). Plots for tof spectra zooming on one resonance and histogram of counts in this resonance divided by protons with statistical errors.

Plots were produced with the macros attached.

Attached are the pictures of enriched Si-30 metal pieces.

The EAR2 ratios of BCT/PKUP, BCT/SILI and PKUP/SILI plotted for each run for samples Si, Sinat, Au20, Au22 and Dummy. 

Also included is the normalised counts (counts/BCT) with 7 detectors on one plot (STED8 was ommited).  'counts' for Si was integreated over the Si resonance, for Au it was integrated over the largest resonance. For Dummy and Sinat the counts were integrated over a large range of 1e+4 to 1e+6ns. The Sinat resonance was too small for adequate statistics to only integrate over the resonance hence why a large range was used. The C6D6 normalised counts were not plotted as there seems to be issues with this detector in EAR2. 

Calibration parameters for each STED detector: With function f=p0+p1*x+p2*x*x

A quadratic function was fitted to the reference point channel (from experimental histogram) and reference point energy (from simulation histogram). The reference point 'r' was r=mean+(FWHM/2) for a gaussian fitted to each compton edge. A quadratic function was used as it fitted the higher energy AmBe reference point better than a linear function. This is highlighted by the comparisson of "STED 2 Calibration" and "STED 2 Example Linear Fit Calibration" plots below. The p0 error is quite large and p2 error very large (larger than the parameter value in STED1) - however the parameters are highly correlated which may have inflated the error. 

CmC histograms were not used as the compton edge could not be located. 

For reference the STED 1 correlation matrix:

Attatched are the overlayed compton edges for the calibrated experimental and the GEANT4 simulations. Note p1 and p2 denotes the lower and higher energy compton edges for Yttrium. 

The experimental and simulated histograms are shifted in quite a few of these plots suggesting the calibration was not fully successful. The compton edge for AmBe was difficult to locate and may have contributed to the poor calibration. 

There seems to be issues with the EAR2 C6D6 T histograms. For Au the earlier run seems to be shifted for detector 2 compared to the rest. 

For Si plotted is 3 different runs 216123, 216128 and 216135 for both detectors. The histogram is messy but highlights that there is a problem. The single Si resonace should be around 20 000ns and this is only the case for run216128. 

Apologies for no legend on plots. 

Au22
det1 run216158 black
det1 run216109 red
det2 run216158 green
det2 run216109 blue

Au22
det1 run216157 black
det1 run216108 red
det2 run216157 green
det2 run216108 blue
 

Below are the histograms of parasitic vs dedicated runs for EAR1. Also included are histograms with ambient background removed which closes the gap between parasitic and dedicated histograms. This was preformed by scaling the ambient by counts by a multiplication of (dedicated or parasitic) bunch number/ambient bunch number. There are no legends but the dedicated counts are in black and parasitic counts in red. Only included are C6D6 1 detectors but the other three detectors were equivelant. 

Here are all the files needed for calculation of the yield. There is some information that still needs to be added to the energyandyield.c file, for example the atomic mass of your target in amu, the neutron separation energy for the compound nucleus in MeV (this is needed due to using the Weighting technique). Also, you need to add some lines  to read the tof histogram, which will then be converted in a neutron energy hisogram and then divided by the flux.

For the conversion to neutron energy we assume an approximate flight path for now. We will determine this more accurately once we have the yields and fit resonances. Emma will need to use a different flux than Annie, because we have a different spallation target. Emma's flux is preliminary, so we also should look at the SILI detectors at some point.

There are some more instances where names for rootfiles and ascii files need to be added. In the first instance, please try and add the missing information and try running the code. I haven't tested it so there may be still problem.  To run it you will need to download the 3 rootfiles into the same folder (or move somewhere else and adjust the path name in the .c file.

PLEASE set variable THRESHSCALER to 1 for now. We dont know yet what it is :-)

Campaign started yesterday. C6D6s placed at 9 cm distance, angle 125 degrees w.r.t beamline.

Measured calibration sources: Cs137, Y88, AmBe, CmC. Measuring Si-nat over the weekend.

Detector          Voltage

Distance Sample to Detector: Det1: ~10.5 cm, Det2: 11 cm, Det3: 11 cm, Det4: 11 cm

My TTOFSort macro is located in

/afs/cern.ch/work/n/nsosnin/public/C6D6_Calibration/Calibration_Data/Calibrator.C

Updated plots for Zn and Au