Proton Calorimetry/Detector Analysis - Revision history

SaadShaikh: /* Photodiode Live Fit */

2021-05-12T16:46:25Z

Photodiode Live Fit

← Older revision		Revision as of 16:46, 12 May 2021
Line 22:		Line 22:
	==Photodiode Live Fit==		==Photodiode Live Fit==
	*See below for details on the fitting procedure. The script <code>liveFit.cpp</code> can be used to fit photodiode data as it is being acquired.		*See below for details on the fitting procedure. The script <code>liveFit.cpp</code> can be used to fit photodiode data as it is being acquired.
			*In addition to background and shoot-through measurements, a configuration file "sheetConfig.txt" should be placed in the data-folder containing the lines:
			**Physical stack thickness in mm
			**Comma-separated values listing the sheet order
			**Stack WET (leave as 0 if not known)
	*To compile the script: <code>g++ -o liveFit liveFit.cpp `root-config --cflags --glibs`</code>		*To compile the script: <code>g++ -o liveFit liveFit.cpp `root-config --cflags --glibs`</code>
	*To execute the script: <code>./liveFit 2 12.5 folder background.txt backST.txt frontST.txt true</code>		*To execute the script: <code>./liveFit 2 12.5 folder background.txt backST.txt frontST.txt true</code>

SaadShaikh: /* Fitting QuARC Data */

2021-05-12T16:45:37Z

Fitting QuARC Data

← Older revision		Revision as of 16:45, 12 May 2021
Line 64:		Line 64:
	==Fitting QuARC Data==		==Fitting QuARC Data==
	*To fit previously acquired photodiode data, the scripts <code>calibrate.cpp</code> and <code>fit.cpp</code> are used.		*To fit previously acquired photodiode data, the scripts <code>calibrate.cpp</code> and <code>fit.cpp</code> are used.
			*In addition to the photodiode data, a configuration file "sheetConfig.txt" should be placed in the data-folder containing the lines:
			**Physical stack thickness in mm
			**Comma-separated values listing the sheet order
			**Stack WET (leave as 0 if not known)
	*To compile the calibration script: <code>g++ -o calibrate calibrate.cpp `root-config --cflags --glibs`</code>		*To compile the calibration script: <code>g++ -o calibrate calibrate.cpp `root-config --cflags --glibs`</code>
	*To execute the calibration script: <code>./calibrate 2 12.5 folder background.txt frontST.txt backST.txt data.txt</code>		*To execute the calibration script: <code>./calibrate 2 12.5 folder background.txt frontST.txt backST.txt data.txt</code>

SaadShaikh: /* Fitting QuARC Data */

2021-05-12T16:43:39Z

Fitting QuARC Data

← Older revision		Revision as of 16:43, 12 May 2021
Line 76:		Line 76:
	**The second argument is the folder containing the measurements.		**The second argument is the folder containing the measurements.
	**The third argument is the name of the calibrated data (which should reside in a directory called "calibrated" in "folder").		**The third argument is the name of the calibrated data (which should reside in a directory called "calibrated" in "folder").
	**The fourth argument is the beam energy in MeV, which is ~~just~~ used as a label for the plot.		**The fourth argument is the facility beam energy in MeV (MeV/u for HIT), which is used as a label for the plot and to load the reference curve.
	**The fifth argument is the facility, used to load the reference curve. Choose either "UCLH" or "HIT".		**The fifth argument is the facility, used to load the reference curve. Choose either "UCLH" or "HIT".
	*The calibration and fitting scripts can be executed consecutively using the bash script <code>fit.sh</code>.		*The calibration and fitting scripts can be executed consecutively using the bash script <code>fit.sh</code>.

SaadShaikh: /* Fitting QuARC Data */

2021-05-12T16:42:53Z

Fitting QuARC Data

← Older revision		Revision as of 16:42, 12 May 2021
Line 72:		Line 72:
	*The output of this script is a file containing two lines: calibrated average photodiode values and standard deviations.		*The output of this script is a file containing two lines: calibrated average photodiode values and standard deviations.
	*To compile the fitting script: <code>g++ -o fit fit.cpp `root-config --cflags --glibs`</code>		*To compile the fitting script: <code>g++ -o fit fit.cpp `root-config --cflags --glibs`</code>
	*To execute the fitting script: <code>./fit 2 ~~12.5~~ folder data.txt energy</code>		*To execute the fitting script: <code>./fit 2 folder data.txt energy facility</code>
	**The first argument is the number of DDC232 boards.		**The first argument is the number of DDC232 boards.
	**The second ~~argument is the DDC232 FSR.~~		**The second argument is the folder containing the measurements.
	**The third argument is the folder containing the measurements.		**The third argument is the name of the calibrated data (which should reside in a directory called "calibrated" in "folder").
	**The ~~fourth~~ argument is the name of the calibrated data (which should reside in a directory called "calibrated" in "folder").		**The fourth argument is the beam energy in MeV, which is just used as a label for the plot.
	**The ~~fifth~~ argument is the beam energy in MeV, which is just used as a label for the plot.		**The fifth argument is the facility, used to load the reference curve. Choose either "UCLH" or "HIT".
	*The calibration and fitting scripts can be executed consecutively using the bash script <code>fit.sh</code>.		*The calibration and fitting scripts can be executed consecutively using the bash script <code>fit.sh</code>.

SaadShaikh: /* Photodiode Replay */

2021-05-11T11:24:14Z

Photodiode Replay

← Older revision		Revision as of 11:24, 11 May 2021
Line 49:		Line 49:
	*This script "replays" a previously acquired measurement, as if it were being captured at that moment.		*This script "replays" a previously acquired measurement, as if it were being captured at that moment.
	*To compile the script: <code>g++ -o replay replay.cpp `root-config --cflags --glibs`</code>		*To compile the script: <code>g++ -o replay replay.cpp `root-config --cflags --glibs`</code>
	*To execute the script: ./replay 2 350 0 170 50 capture.txt step/avg 50 100 (clean folder background.txt backST.txt frontST.txt)		*To execute the script: <code>./replay 2 350 0 170 50 capture.txt step/avg 50 100 (clean folder background.txt backST.txt frontST.txt)</code>
	**The first argument is the number of DDC232 boards.		**The first argument is the number of DDC232 boards.
	**The second argument is the DDC232 FSR.		**The second argument is the DDC232 FSR.

SaadShaikh at 11:23, 11 May 2021

2021-05-11T11:23:35Z

SaadShaikh: /* Photodiode Data Processing */

2021-05-11T11:20:21Z

Photodiode Data Processing

@@ Line 8: / Line 8: @@
 **As each photodiode is split across two inputs (as of rev. B), the 32 total measurements become 16 measurements, where inputs of the same photodiode are summed.
 **Practically, this makes the full-scale range twice that of what is set on the FPGA.
 ==Photodiode Noise Analysis==

SaadShaikh at 11:06, 11 May 2021

2021-05-11T11:06:29Z

← Older revision		Revision as of 11:06, 11 May 2021
Line 36:		Line 36:
	**The tenth argument "clean" is optional and can be used to replay photodiode data with previously measured background and calibration corrections.		**The tenth argument "clean" is optional and can be used to replay photodiode data with previously measured background and calibration corrections.
	***The eleventh, argument is the folder the background and front and back shoot-through measurements are in. The remaining arguments are the filenames.		***The eleventh, argument is the folder the background and front and back shoot-through measurements are in. The remaining arguments are the filenames.


	==Fitting QuARC Data==		==Fitting QuARC Data==
Line 55:		Line 54:
	**The fifth argument is the beam energy in MeV, which is just used as a label for the plot.		**The fifth argument is the beam energy in MeV, which is just used as a label for the plot.
	*The calibration and fitting scripts can be executed consecutively using the bash script <code>fit.sh</code>.		*The calibration and fitting scripts can be executed consecutively using the bash script <code>fit.sh</code>.


	===Useful literature===		===Useful literature===

SaadShaikh: /* Fitting QuARC Data */

2021-05-11T11:06:11Z

Fitting QuARC Data

@@ Line 39: / Line 39: @@
 ==Fitting QuARC Data==
 ===Useful literature===
 *Laurent's paper on the [https://doi.org/10.1002/mp.14099 quenched Bragg model].
@@ Line 50: / Line 68: @@
 ===Notes on Fitting===
 *Sigma is estimated with Bortfeld's approximation, which is only truly valid for protons.
 **While sigma can be estimated from a Gaussian fit of the Bragg peak, the poor spatial resolution of the detector means that Bortfeld's estimation tends to work better for all particle species.

SaadShaikh: /* Photodiode Noise Analysis */

2021-05-11T10:44:30Z

Photodiode Noise Analysis

← Older revision		Revision as of 10:44, 11 May 2021
Line 12:		Line 12:
	*The script <code>analysis.cpp</code> analyses DDC232 data to report the amount of external electronic noise present in the data. This is useful in determining the quality of the data.		*The script <code>analysis.cpp</code> analyses DDC232 data to report the amount of external electronic noise present in the data. This is useful in determining the quality of the data.
	*To compile the script: <code>g++ -o analysis analysis.cpp `root-config --cflags --glibs`</code>		*To compile the script: <code>g++ -o analysis analysis.cpp `root-config --cflags --glibs`</code>
	*To execute the script: <code>./analysis 4 folder capture.txt</code		*To execute the script: <code>./analysis 4 folder capture.txt</code>
	**The first argument is the number of DDC232 boards.		**The first argument is the number of DDC232 boards.
	**The second argument is the folder containing the captured data.		**The second argument is the folder containing the captured data.

← Older revision		Revision as of 11:23, 11 May 2021
Line 11:		Line 11:
	==Photodiode Live Plot==		==Photodiode Live Plot==
	*The script <code>livePlot.cpp</code> can be used to plot photodiode data live while being acquired by <code>FTDI.cpp</code>.		*The script <code>livePlot.cpp</code> can be used to plot photodiode data live while being acquired by <code>FTDI.cpp</code>.
	*To compile the script: g++ -o livePlot livePlot.cpp `root-config --cflags --glibs`		*To compile the script: <code>g++ -o livePlot livePlot.cpp `root-config --cflags --glibs`</code>
	*To execute the script: ./livePlot 1 350 0 (clean folder background.txt backST.txt frontST.txt)		*To execute the script: <code>./livePlot 1 350 0 (clean folder background.txt backST.txt frontST.txt)</code>
	**The first argument is the number of DDC232 boards.		**The first argument is the number of DDC232 boards.
	**The second argument is the DDC232 FSR.		**The second argument is the DDC232 FSR.
Line 22:		Line 22:
	==Photodiode Live Fit==		==Photodiode Live Fit==
	*See below for details on the fitting procedure. The script <code>liveFit.cpp</code> can be used to fit photodiode data as it is being acquired.		*See below for details on the fitting procedure. The script <code>liveFit.cpp</code> can be used to fit photodiode data as it is being acquired.
	*To compile the script: g++ -o liveFit liveFit.cpp `root-config --cflags --glibs`		*To compile the script: <code>g++ -o liveFit liveFit.cpp `root-config --cflags --glibs`</code>
	*To execute the script: ./liveFit 2 12.5 folder background.txt backST.txt frontST.txt true		*To execute the script: <code>./liveFit 2 12.5 folder background.txt backST.txt frontST.txt true</code>
	**The first argument is the number of DDC232 boards.		**The first argument is the number of DDC232 boards.
	**The second argument is the DDC232 FSR.		**The second argument is the DDC232 FSR.