Proton Calorimetry/User Interface

From PBTWiki
Jump to navigation Jump to search

To do list for GUI

  • Implement start/stop button
  • Stop button functionality to pause
  • Check GUI access with NUC PC (install apache, root, FTDI libraries)
  • Fix cache issue for Chrome
  • Replay script: run replay script on GUI and see if that works
  • Trigger replay script with GUI: add dropdown menu and space to specify all the variables needed for the replay code/


Files for GUI visualisation

The experimental data acquired needs to be calibrated with the acquired shoot-through measurements (back and front) and the background needs to be subtracted. Then, each of the measurements can be fitted to the models, extracting the range and the energy of the beam. There are two different codes to perform this.

Calibrated data: Development of a python code (replay.py) that generates a file with header info followed by the 30k measurements calibrated and with background subtracted. If our initial file is called Run0xy.txt the output of this will be Run0xy_calibrated.txt

python3 replay.py t_INT framerate skipahead reverse folder run background backST frontST

Fitted data: Once the data has been calibrated, we need to perform the fittings and extract the parameters. This code (fit.cpp) utilises Run_0xy.txt and outputs Run0xy_fitted_frequency.txt, where the frequency value is 25Hz by default but can be changed by adding an extra argument to the code line. This takes about 15 second per run for Clatterbridge data.

g++ -o fit fit.cpp `root-config --cflags --glibs` -O3
./fit reverse energy facility folder run (frequency)​

​ Both, fitted and calibrated data for all the runs acquired at Clatterbridge experiment in 2023 are in https://www.hep.ucl.ac.uk/pbt/wikiData/data/GUI-Clatterbridge/


Live Update: The webpage is reading from file outputGUI.csv, located in the same directory than the GUI js and html files. To visualise it live, we create a shell script that goes through each measurement at a certain time interval, writing the output to the outputGUI.csv file. To execute this, simply we cd into ReplayFinal and then execute the script.

./shellGUI.sh

GUI version 2.1.b

This version has fixed all of the issues we had in previous versions. Design is the same but the code has been rewritten completely. The key improvements are:

  • Zooming feature: Shifting between PDdata and curves when zooming in: two different axis were used for PDdata and fitted data. This has been modified to have only one xScale, defined as a linear scale. This means that the x-axis now utilises numerical datasets instead of categorical, and therefore the zooming of both sets of data is performed equally at the same time.

var xScale = d3.scaleLinear()

.domain([0, setWET])

.range([margins.left, width - margins.right ]);

  • Bandwidth separation: Previous version assumed that the separation between PDdata was uniform, but this did not correspond to reality as the PD thicknesses varies. The bandwidth needs to be redefined based on the separation between each PD and the previous one.

bandwidths.push(xPDposition[i+1] - xPDposition[i]);


This version has been developed for Clatterbridge data (12PD) - If using for larger data sets we need to modify the numer of setPhotodiodes and setWET in GUI.js

GUI version 2.06

Located on: https://www.hep.ucl.ac.uk/pbt/PDdisplay/escribano/v2.06/

Improved version based on v 2.05. The start and stop buttons have been changed by a single button that performs both actions, keeping in screen the results from the last data plotted instead of deleting the graph and beam parameters.

Live update + fitting curves.

GUI version 2.05

This version of the GUI can be seen on: https://www.hep.ucl.ac.uk/pbt/PDdisplay/escribano/v2.05/

General information

Developed in JavaScript with d3 and papaparse libraries. This version consists of a single tab, where the user can adjust WET number (mm) and max Y (pC) and select which elements are plotted on the graph, via checkboxes.

Data used for testing GUI: 245 MeV, acquired at The Christie in October 2022

Disclaimer!: The error bars from this version do not correspond to the real error value, as the magnitude has been increased by 5pC for visualisation purposes (defined in the code as extraerror).


Live update

This version provides a live visualisation of the data file. In order to generate a dynamic set of data this two scripts have been developed in Python: GUIrandom.py and runauto.py.

- GUIrandom.py: calculates a random number and multiplies the PD data and the curves by this number. The random number is shown in the terminal.

- runauto.py: runs GUIrandom.py at the desired time interval (initially 0.02s, to have a refresh rate of 50Hz)

Uploading Data to HEP Server

Data can be uploaded to the HEP cluster in order to run a webpage with live photodiode data.

  • On your local machine:
    • Move to the directory containing photodiode data being saved by CoolTerm and the C++ script to generate the CSV file: cd /Path/To/Data
    • Generate the photodiode data with: root -l "writeCSV.cpp"
    • In a new terminal, move to the same directory and run the following command in order to copy this data, line by line, to the plus1 web area: tail -F -n 1 output.csv | ssh username@plus2.hep.ucl.ac.uk -T 'cat - > /unix/www/html/pbt/PDdisplay/photodiodeValuesAll.csv'
  • On plus2 :
    • In a new terminal, move into the directory where the photodiode data is written: cd /unix/www/html/pbt/PDdisplay/
    • Run the following script in that directory: ./outputSingleLine.sh -v photodiodeValuesAll.csv photodiodeValues.csv 0.02
      • This will output a single line CSV file - photodiodeValues.csv - every 10ms (100Hz). It's a good idea to set the write rate about 2x faster than the desired frame rate.
    • To output this any more or less often, adjust the "waitTime" value at the end (needs to be in seconds).
    • To stop displaying each line as it is written, remove the -v flag.
    • Note that this script contains an infinite loop, so to stop it running you need to kill it by hand with Ctrl-C. You may need to do this a few times to get it to quit...