DPE: Difference between revisions
(51 intermediate revisions by the same user not shown) | |||
Line 7: | Line 7: | ||
* '''Processing''': particle recognition and radiation field recognition | * '''Processing''': particle recognition and radiation field recognition | ||
* '''Post-processing''': directional analysis, coincidence analysis, Compton camera, generation of physical products and their time evolution with respect to the radiation field classification (fluxes, dose rates etc.) | * '''Post-processing''': directional analysis, coincidence analysis, Compton camera, generation of physical products and their time evolution with respect to the radiation field classification (fluxes, dose rates etc.) | ||
The following information describes its '''command line platform''' and it is '''valid for version''' '''1.1.0'''. | |||
<gallery mode="slideshow" showthumbnails=""> | <gallery mode="slideshow" showthumbnails=""> | ||
File:Pid.png|Example of particle classification for data acquired onn board of ISS with TPX3 500um. | File:Pid.png|Example of particle classification for data acquired onn board of ISS with TPX3 500um. | ||
Line 17: | Line 17: | ||
</gallery> | </gallery> | ||
== Download == | == Download == | ||
The data processing engine can be downloaded from the following sources | |||
<!-- | |||
<span style="color: red;">'''Apology for Temporary Download Unavailability Due to Maintenance'''</span> | |||
--> | |||
<span style="color: black;">'''In the case of interest, please contact: support@advacam.cz or lukas.marek@advacam.cz'''</span> | |||
<!-- | |||
The data processing engine can be downloaded from the following sources based on the needed computer configuration: | |||
{| class="wikitable" | {| class="wikitable" | ||
Line 46: | Line 54: | ||
|} | |} | ||
Examples of DPE usage for sevetal cases can be downloaded from the following link: | |||
{| class="wikitable" | {| class="wikitable" | ||
|+ | |+ | ||
Line 55: | Line 63: | ||
|[[Media:DPE_Example.zip]] | |[[Media:DPE_Example.zip]] | ||
|} | |} | ||
--> | |||
== Prerequisites, Install and Run== | == Prerequisites, Install and Run== | ||
=== Prerequisites === | |||
To create graphics, '''python3''' and several python packages are needed: | |||
To create graphics, '''python3''' and several packages are needed: | |||
* matplotlib (version >= 3.4) | * matplotlib (version >= 3.4) | ||
* numpy | * numpy | ||
* multiprocessing - only if | * multiprocessing - only if creation of graphics should be done on several threads/cores of CPU | ||
The program, especially in the case of graphic export, needs several hundreds of MB on RAM (app 500 MB but it is based on the size of exported histograms). If the multiprocessing is used then the RAM usage can be up several GB and all threads of the CPU | The program, especially in the case of graphic export, needs several '''hundreds of MB on RAM''' (app 500 MB but it is based on the size of exported histograms). If the multiprocessing is used then the '''RAM usage can be up several GB''' and all threads of the CPU might be used. | ||
=== Install === | |||
In the directory where the program should be installed, extract all the files from downloaded archive/zip. | |||
=== Run the program on linux === | === Run the program on linux === | ||
The DPE is started with the following program in the command line from the installation directory: <syntaxhighlight lang="sh"> ./dpe.sh PARAMETERS_FILE_PATH/PARATERS_FILE_NAME</syntaxhighlight> | |||
<syntaxhighlight lang="sh"> | <code>PARAMETERS_FILE_PATH</code> is path to the parameters file and PARATERS_FILE_NAME is its name. It is possible that dpe.sh and clusterer have to be allowed as an executable: <code>chmod +x clusterer</code>. | ||
./dpe.sh PARAMETERS_FILE_PATH/PARATERS_FILE_NAME</syntaxhighlight> | |||
It is possible that | |||
=== Run the program on windows === | === Run the program on windows === | ||
Line 171: | Line 176: | ||
== Issues and Help== | == Issues and Help== | ||
Any ideas or issues can be reported on github forum for issues: https://github.com/lmareksla/DPE_Issues/ or send directly on the first mail below | Any ideas or issues can be reported on github forum for issues: | ||
https://github.com/lmareksla/DPE_Issues/ | |||
or send directly on the first mail below. | |||
For more information or help: lukas.marek@advacam.cz, carlos.granja@advacam.cz. | For more information or help: '''lukas.marek@advacam.cz''', '''carlos.granja@advacam.cz'''.<gallery mode="packed-overlay" heights="200"> | ||
File:Dpe git issue 1.png | |||
File:Dpe git issue 2.png | |||
File:Dpe git issue 3.png | |||
</gallery> | |||
== Main Configuration of DPE== | == Main Configuration of DPE== | ||
Line 194: | Line 207: | ||
=== List of Parameters === | === List of Parameters === | ||
{| class="wikitable sortable" | |||
!Name | |||
!Type | |||
!Description | |||
!Example | |||
!Default | |||
|- | |||
|'''CalMat''' | |||
|STRING | |||
|Path to a directory with calibration matrices if data should be calibrated during pre-processing. | |||
|<code>CalMat="/Path/To/Cal/Matrices/"</code><br> | |||
|empty string | |||
|- | |||
|'''ClogOutName''' | |||
|STRING | |||
|Name of the cluster log file for export without suffix/ending. If it is stated in the parameters file then the clog is created (alias DoCreateClog = true). | |||
|<code>ClogOutName="ClusterLog"</code><br> | |||
|ClusterLog | |||
|- | |||
|'''ClusterCount''' | |||
|INT | |||
|Count of clusters which should be processed and evaluated. | |||
|<code>ClusterCount=1</code><br><code>ClusterCount=10000</code><br><code>ClusterCount=5</code><br> | |||
|-1 (unused) | |||
|- | |||
|'''ClusterCountInSensImage''' | |||
|INT | |||
|Count of clusters which will be exported in the integrated sensor plots (also per class). It can be less, if there is no such a count of needed clusters set by this value. | |||
|<code>ClusterCountInSensImage=10</code><br><code>ClusterCountInSensImage=1000</code><br> | |||
|100 | |||
|- | |||
|'''ClustererName''' | |||
|STRING | |||
|Name o the clusterer binary. Default, it is decided based on the operation system. | |||
|<code>ClustererName="clusterer"</code><br><code>ClustererName="clusterer.exe"</code><br> | |||
|clusterer (linux) or clusterer.exe (windows) | |||
|- | |||
|'''ClustererPath''' | |||
|STRING | |||
|Path to the clusterer program (part of the download package). | |||
|<code>ClustererPath="/Path/To/Clusterer/ Program/"</code><br> | |||
|./Clusterer/ | |||
|- | |||
|'''ClusterImageCount''' | |||
|INT | |||
|Count of clusters which will be exported as individual cluster plots. | |||
|<code>ClusterImageCount=20</code><br><code>ClusterImageCount=1</code><br> | |||
|15 | |||
|- | |||
|'''CoincEventListBaseName''' | |||
|STRING | |||
|Base name of the coinc_event list which is used for txt and also for the json format. It is combined as this name and adding .txt and .json. | |||
|<code>CoincEventListBaseName="CoincEventList"</code><br><code>CoincEventListBaseName="List"</code><br><code>CoincEventListBaseName="Name"</code><br> | |||
|CoincEventList | |||
|- | |||
|'''CoincEventListJsonName''' | |||
|STRING | |||
|Name of output coinc_event list in json format. | |||
|<code>CoincEventListJsonName="SampligList.json"</code><br> | |||
|CoincEventList.json | |||
|- | |||
|'''CoincEventListName''' | |||
|STRING | |||
|Name of the output coinc_event list. | |||
|<code>CoincEventListName="CoincEventList.txt"</code><br> | |||
|CoincEventList.txt | |||
|- | |||
|'''ComptonListBaseName''' | |||
|STRING | |||
|Base name of the compton list which is used for txt and also for the json format. It is combined as this name and adding .txt and .json. | |||
|<code>ComptonListBaseName="ComptonList"</code><br><code>ComptonListBaseName="List"</code><br><code>ComptonListBaseName="Name"</code><br> | |||
|ComptonList | |||
|- | |||
|'''ComptonListJsonName''' | |||
|STRING | |||
|Name of output compton list in json format. | |||
|<code>ComptonListJsonName="SampligList.json"</code><br> | |||
|ComptonList.json | |||
|- | |||
|'''ComptonListName''' | |||
|STRING | |||
|Name of the output compton list. | |||
|<code>ComptonListName="ComptonList.txt"</code><br> | |||
|ComptonList.txt | |||
|- | |||
|'''ComptonName''' | |||
|STRING | |||
|Name of the file with Compton direction reconstruction configuration. | |||
|<code>ComptonName="Compton.ini"</code><br> | |||
|empty string | |||
|- | |||
|'''ComptonPath''' | |||
|STRING | |||
|Path to the file with Compton direction reconstruction configuration. | |||
|<code>ComptonPath="/Path/To/Compton/Config/"</code><br> | |||
|empty string | |||
|- | |||
|'''DetGeoName''' | |||
|STRING | |||
|Name of the file with detector geometry configuration. | |||
|<code>DetGeoName="DetGeo.ini"</code><br> | |||
|empty string | |||
|- | |||
|'''DetGeoPath''' | |||
|STRING | |||
|Path to the file with detector geometry configuration. | |||
|<code>DetGeoPath="/Path/To/DetGeo/Config/"</code><br> | |||
|empty string | |||
|- | |||
|'''DirectionListBaseName''' | |||
|STRING | |||
|Base name of the direction list which is used for txt and also for the json format. It is combined as this name and adding .txt and .json. | |||
|<code>DirectionListBaseName="DirectionList"</code><br><code>DirectionListBaseName="List"</code><br><code>DirectionListBaseName="Name"</code><br> | |||
|DirectionList | |||
|- | |||
|'''DirectionListJsonName''' | |||
|STRING | |||
|Name of output direction list in json format. | |||
|<code>DirectionListJsonName="SampligList.json"</code><br> | |||
|DirectionList.json | |||
|- | |||
|'''DirectionListName''' | |||
|STRING | |||
|Name of the output direction list. | |||
|<code>DirectionListName="DirectionList.txt"</code><br> | |||
|NClusters | |||
|- | |||
|'''DirectionName''' | |||
|STRING | |||
|Name of the file with directional analysis configuration. | |||
|<code>DirectionName="Direction.ini"</code><br> | |||
|empty string | |||
|- | |||
|'''DirectionPath''' | |||
|STRING | |||
|Path to the file with directional analysis configuration. | |||
|<code>DirectionPath="/Path/To/Direction/Config/"</code><br> | |||
|empty string | |||
|- | |||
|'''DoClustererLog''' | |||
|BOOL | |||
|Export clusterer log. Only valid if clusterer is used during processing (raw data). | |||
|<code>DoClustererLog=true</code><br><code>DoClustererLog=false</code><br> | |||
|true (used) | |||
|- | |||
|'''DoCompton''' | |||
|BOOL | |||
|Control whether Compton directional reconstruction should be done. There is limited list of cases when this analysis is usable based on the detector settings (see section below). | |||
|<code>DoCompton=true</code><br><code>DoCompton=false</code><br> | |||
|true (if configuration allows) | |||
|- | |||
|'''DoCreateClog''' | |||
|BOOL | |||
|Control for a generation of the clog. | |||
|<code>DoCreateClog=true</code><br><code>DoCreateClog=false</code><br> | |||
|true (used) | |||
|- | |||
|'''DoCreateElistExt''' | |||
|BOOL | |||
|Switch to produce elist with additional columns from processing output: filter pass or not passed (1,0), PID class (e.g. from -1 to 8). It is placed as the last columns of the elist with name `FilterPass, PIDClass`. | |||
|<code>DoCreateElistExt=FileIn_Count</code><br> | |||
|true (used) | |||
|- | |||
|'''DoDirection''' | |||
|BOOL | |||
|Control whether Compton directional reconstruction should be done. There is limited list of cases when this analysis is usable based on the detector settings (see section below). | |||
|<code>DoDirection=true</code><br><code>DoDirection=false</code><br> | |||
|true | |||
|- | |||
|'''DoDirectionTrackCond''' | |||
|BOOL | |||
|Control to use tracking condition during directional analysis (additional constrain on features of clusters). | |||
|<code>DoDirectionTrackCond=true</code><br><code>DoDirectionTrackCond=false</code><br> | |||
|false | |||
|- | |||
|'''DoDoseUnitRadiology''' | |||
|BOOL | |||
|Switch to change from the dose rate in uGy/h to Gy/s. | |||
|<code>DoDoseUnitRadiology=true</code><br><code>DoDoseUnitRadiology=false</code><br> | |||
|false | |||
|- | |||
|'''DoExportElistExtFilter''' | |||
|BOOL | |||
|Control whether clusters which were filtered out should not be exported into extended elist. | |||
|<code>DoExportElistExtFilter=true</code><br><code>DoExportElistExtFilter=false</code><br> | |||
|false (exported) | |||
|- | |||
|'''DoExportElistExtSensEdge''' | |||
|BOOL | |||
|Control whether clusters which were at sensor edge should not be exported into extended elist. | |||
|<code>DoExportElistExtSensEdge=true</code><br><code>DoExportElistExtSensEdge=false</code><br> | |||
|false (exported) | |||
|- | |||
|'''DoExportGraphics''' | |||
|BOOL | |||
|Control for an export of graphical output. | |||
|<code>DoExportGraphics=true</code><br><code>DoExportGraphics=false</code><br> | |||
|true (used) | |||
|- | |||
|'''DoExportText''' | |||
|BOOL | |||
|Control of an export of the output into files. If turned off, then no export into files is done. | |||
|<code>DoExportText=true</code><br><code>DoExportText=false</code><br> | |||
|true (used) | |||
|- | |||
|'''DoFilter''' | |||
|BOOL | |||
|Control whether filtering should be done during processing. | |||
|<code>DoFilter=true</code><br><code>DoFilter=false</code><br> | |||
|true (used) | |||
|- | |||
|'''DoFrame''' | |||
|BOOL | |||
|Control to do frame analysis during processing. | |||
|<code>DoFrame=true</code><br><code>DoFrame=false</code><br> | |||
|true (used) | |||
|- | |||
|'''DoGraphicsMultiThread''' | |||
|BOOL | |||
|Switch to use multiprocessing/multitreading during graphics export. | |||
|<code>DoGraphicsMultiThread=true</code><br><code>DoGraphicsMultiThread=false</code><br> | |||
|true (used) | |||
|- | |||
|'''DoCheckPrereq''' | |||
|BOOL | |||
|Control to do or not do check of prerequisites as clusterer or python modules. | |||
|<code>DoCheckPrereq=true</code><br><code>DoCheckPrereq=false</code><br> | |||
|true | |||
|- | |||
|'''DoIgnoreClusterSensEdge''' | |||
|BOOL | |||
|Switch to ignore cluster at sensor edge. They are not used into histograms, physical products etc. (evaluation) but they remain in `Files` from clusterer. | |||
|<code>DoIgnoreClusterSensEdge=true</code><br><code>DoIgnoreClusterSensEdge=false</code><br> | |||
|false | |||
|- | |||
|'''DoLog''' | |||
|BOOL | |||
|Control for creating log into file. | |||
|<code>DoLog=true</code><br><code>DoLog=false</code><br> | |||
|true (used) | |||
|- | |||
|'''DoLogDateTime''' | |||
|BOOL | |||
|Control whether date and time should be in the log file. | |||
|<code>DoLogDateTime=true</code><br><code>DoLogDateTime=false</code><br> | |||
|true (used) | |||
|- | |||
|'''DoMaskFullPrint''' | |||
|BOOL | |||
|Switch to fully print the mask content into terminal and log file. | |||
|<code>DoMaskFullPrint=true</code><br><code>DoMaskFullPrint=false</code><br> | |||
|false | |||
|- | |||
|'''DoMultiFile''' | |||
|BOOL | |||
|Check whether the multi file processing should be used (can be turned off/false and on/true). | |||
|<code>DoMultiFile=true</code><br><code>DoMultiFile=false</code><br> | |||
|false | |||
|- | |||
|'''DoPID''' | |||
|BOOL | |||
|Switch to turn off the particle identification | |||
|<code>DoPID=true</code><br><code>DoPID=false</code><br> | |||
|true (used) | |||
|- | |||
|'''DoPrint''' | |||
|BOOL | |||
|Control for print progress into terminal. | |||
|<code>DoPrint=true</code><br><code>DoPrint=false</code><br> | |||
|true (used) | |||
|- | |||
|'''DoRadFieldRecog''' | |||
|BOOL | |||
|Control whether radiation filed recognition should be done. There is limited list of cases when this analysis is usable based on the detector settings (see section below). | |||
|<code>DoRadFieldRecog=true</code><br><code>DoRadFieldRecog=false</code><br> | |||
|true (if configuration allows) | |||
|- | |||
|'''DoRemoveElist''' | |||
|BOOL | |||
|Switch to delete elist created by the clusterer but the extended elist is not effected by this switch. If set to true, the elist is removed. It set to true because the information is doubles and extended in the elist extended. | |||
|<code>DoRemoveElist=true</code><br><code>DoRemoveElist=false</code><br> | |||
|true (used) | |||
|- | |||
|'''DoRemoveOldElist''' | |||
|BOOL | |||
|Switch to explicitly delete or not an already existing elist from a previous processing. If set to true, the elist is removed and a new one is created. If false, the DPE checks elist existence and if there is a elist, the pre-processing with clusterer is skipped. The default value is true because this could cause an error in processing if mask is used with false option -> it wont be processed again and the mask will be ignored for elist data. | |||
|<code>DoRemoveOldElist=true</code><br><code>DoRemoveOldElist=false</code><br> | |||
|true (used) | |||
|- | |||
|'''DoRunClusterer''' | |||
|BOOL | |||
|Switch to skip processing of raw data with clusterer. | |||
|<code>DoRunClusterer=true</code><br><code>DoRunClusterer=false</code><br> | |||
|true (used) | |||
|- | |||
|'''DoSensMatrixCount''' | |||
|BOOL | |||
|Switch to use counts in sensor plots instead of energy. | |||
|<code>DoSensMatrixCount=true</code><br><code>DoSensMatrixCount=false</code><br> | |||
|false | |||
|- | |||
|'''DoSigVec''' | |||
|BOOL | |||
|Control whether significant vector creation should be done. | |||
|<code>DoSigVec=true</code><br><code>DoSigVec=false</code><br> | |||
|true | |||
|- | |||
|'''DoTpxHighEnergyCorr''' | |||
|BOOL | |||
|Switch to use TPX high energy correction. | |||
|<code>DoTpxHighEnergyCorr=true</code><br><code>DoTpxHighEnergyCorr=false</code><br> | |||
|false | |||
|- | |||
|'''ElistExtOutName''' | |||
|STRING | |||
|Output name of the elist extended which also includes information as PID class etc. Ending/suffix can be stated and is used in this case (in contrast to ElistOutName). | |||
|<code>ElistExtOutName="ElistExt.advelist"</code><br><code>ElistExtOutName="EExt.advelist"</code><br> | |||
|EventListExt.advelist | |||
|- | |||
|'''ElistOutName''' | |||
|STRING | |||
|Name of the elist file for export without the suffix. | |||
|<code>ElistOutName="EventList"</code><br> | |||
|EventList | |||
|- | |||
|'''FileInCount''' | |||
|INT | |||
|Count of input files which should be processed if multi file processing is used. | |||
|<code>FileInCount=10</code><br> | |||
|-1 (unused) | |||
|- | |||
|'''FileInName''' | |||
|VSTRING | |||
|Name of the input file. Spaces can be part of the name. | |||
|<code>FileInName=tot toa.t3pa</code><br> | |||
|empty string | |||
|- | |||
|'''FileInNameEnd''' | |||
|STRING | |||
|Name end/extension of the input file. If the file name is `FileIn.txt` then name end is `.txt`. It is used for a several file processing. | |||
|<code>FileInNameEnd=".txt"</code><br> | |||
|empty string | |||
|- | |||
|'''FileInPath''' | |||
|STRING | |||
|Path of the input file(s). Directory can be also given, then processing of files in this dir is done based on valid type of files. | |||
|<code>FileInPath="/Path/ To/ File/"</code><br> | |||
|empty string | |||
|- | |||
|'''FileOutName''' | |||
|STRING | |||
|Name of the output sampling list. | |||
|<code>FileOutName="SamplingList.txt"</code><br> | |||
|SamplingList.txt | |||
|- | |||
|'''FileOutPath''' | |||
|STRING | |||
|Path for the results. All results will be exported into this directory and if the path/dir exists and it can be created. | |||
|<code>FileOutPath="/Path/File /Out/"</code><br> | |||
|current working directory | |||
|- | |||
|'''FilterName''' | |||
|STRING | |||
|Name of the file with filter configuration. | |||
|<code>FilterName="Filter.ini"</code><br> | |||
|empty string | |||
|- | |||
|'''FilterPath''' | |||
|STRING | |||
|Path to the file with filter configuration. | |||
|<code>FilterPath="/Path/To/Filter/Config/"</code><br> | |||
|empty string | |||
|- | |||
|'''FrameListBaseName''' | |||
|STRING | |||
|Base name of the frame list which is used for txt and also for the json format. It is combined as this name and adding .txt and .json. | |||
|<code>FrameListBaseName="FrameList"</code><br><code>FrameListBaseName="List"</code><br><code>FrameListBaseName="Name"</code><br> | |||
|FrameList | |||
|- | |||
|'''FrameListJsonName''' | |||
|STRING | |||
|Name of output frame list in json format. | |||
|<code>FrameListJsonName="SampligList.json"</code><br> | |||
|FrameList.json | |||
|- | |||
|'''FrameListName''' | |||
|STRING | |||
|Name of the output frame list. | |||
|<code>FrameListName="FrameList.txt"</code><br> | |||
|FrameList.txt | |||
|- | |||
|'''GraphicsMultiThreadFileSize''' | |||
|FLOAT | |||
|File size limit in B on single multi processing batch. It should be decreased if the RAM is overwhelmed during processing. | |||
|<code>GraphicsMultiThreadFileSize=2e6</code><br><code>GraphicsMultiThreadFileSize=1e6</code><br> | |||
|3e6 (30 MB) | |||
|- | |||
|'''Hist1DGraphicsRebin''' | |||
|INT | |||
|Rebin value for 1D histograms before plotting. It will rebin all 1D histograms to make them more clear and exporting graphics faster if value around 1e2 is used. If value 0 is set then this option is skipped. This option overwrite a possible user settings in the configuration files for histograms, but it is valid only for exported graphics, not for text files. | |||
|<code>Hist1DGraphicsRebin=100</code><br><code>Hist1DGraphicsRebin=2</code><br> | |||
|0 | |||
|- | |||
|'''HistName''' | |||
|STRING | |||
|Name of the file with significant vector configuration. | |||
|<code>HistName="Hist.ini"</code><br> | |||
|empty string | |||
|- | |||
|'''HistPath''' | |||
|STRING | |||
|Path to the file with significant vector configuration. | |||
|<code>HistPath="/Path/To/Hist/Config/"</code><br> | |||
|empty string | |||
|- | |||
|'''ChipType''' | |||
|STRING | |||
|Type of the chip (TPX, TPX2, TPX3). Example: `ChipType = "TPX3"` to specify that TPX3 was used chip. | |||
|<code>ChipType="TPX"</code><br><code>ChipType="TPX2"</code><br><code>ChipType="TPX3"</code><br><code>ChipType="TPX4"</code><br> | |||
|TPX | |||
|- | |||
|'''IsEnergyCalibrated''' | |||
|BOOL | |||
|Control to inform DPE that data are already energy calibrated. | |||
|<code>IsEnergyCalibrated=true</code><br><code>IsEnergyCalibrated=false</code><br> | |||
|false (uncalibrated) | |||
|- | |||
|'''LengthCorr''' | |||
|INT | |||
|Num switch to choose formula for calculation of cluster length in sensor plane. Possible values: | |||
* 1 - based on weighted standard deviation subtraction from projected length. | |||
* 2 - based on projected width subtraction from projected length. | |||
* 3 - based on model of Nabha. | |||
|<code>LengthCorr=1</code><br><code>LengthCorr=2</code><br><code>LengthCorr=3</code><br> | |||
|1 | |||
|- | |||
|'''LogName''' | |||
|STRING | |||
|Name of the log file. | |||
|<code>LogName="Log.txt"</code><br><code>LogName="file.hop"</code><br> | |||
|Log.txt | |||
|- | |||
|'''LogPath''' | |||
|STRING | |||
|Path to the log file. | |||
|<code>LogPath="/Path/To/Log/ File/"</code><br> | |||
|current working directory | |||
|- | |||
|'''MaskName''' | |||
|STRING | |||
|Name of the file with a mask configuration. | |||
|<code>MaskName="Mask.txt"</code><br> | |||
|empty string | |||
|- | |||
|'''MaskPath''' | |||
|STRING | |||
|Path to the file with a mask configuration. | |||
|<code>MaskPath="/Path/To/Mask/Config/"</code><br> | |||
|empty string | |||
|- | |||
|'''MeasMode''' | |||
|STRING | |||
|Measurement mode of detector = itot+count (itot+count). In the current version, it is used for itot+count mode recognition. | |||
|<code>MeasMode="itot+count"</code><br> | |||
|empty string | |||
|- | |||
|'''ModelPath''' | |||
|STRING | |||
|Path to models. | |||
|<code>ModelPath="/Path/To/Models/"</code><br> | |||
|./Models/ | |||
|- | |||
|'''PIDAlg''' | |||
|INT | |||
|Switch to change between different PID algorithms. All possibilities (101,102,103,201,202): | |||
* 101 - 4 class heuristic decision tree, TPX 300 um Si | |||
* 102 - 9 class heuristic decision tree, TPX 300 um Si | |||
* 103 - 16 class heuristic decision tree, TPX3 500 um Si | |||
* ...171-173 new filters for neutrons | |||
* 201 - 3 class DNN, TPX 300 um Si 30 V | |||
* 202 - 6 class DNN, TPX 300 um Si 30 V | |||
* 251 - 3 class DNN, TPX3 500 um Si 80 V | |||
* 252 - 6 class DNN, TPX3 500 um Si 80 V | |||
|<code>PIDAlg=252</code><br> | |||
|-1 (choosen based on detector settings) | |||
|- | |||
|'''RadFieldRecogDatabasePath''' | |||
|STRING | |||
|Path to database for radiation filed recognition via comparator. This database should include significant vectors which are used for comparison. | |||
|<code>RadFieldRecogDatabasePath="/Path/To/Comp/DB/"</code><br> | |||
|empty string (default database) | |||
|- | |||
|'''RadFieldRecogListBaseName''' | |||
|STRING | |||
|Base name of the rfr list which is used for txt and also for the json format. It is combined as this name and adding .txt and .json. | |||
|<code>RadFieldRecogListBaseName="RadFieldRecogList"</code><br><code>RadFieldRecogListBaseName="List"</code><br><code>RadFieldRecogListBaseName="Name"</code><br> | |||
|RadFieldRecogList | |||
|- | |||
|'''RadFieldRecogListJsonName''' | |||
|STRING | |||
|Name of output rfr list in json format. | |||
|<code>RadFieldRecogListJsonName="SampligList.json"</code><br> | |||
|RadFieldRecogList.json | |||
|- | |||
|'''RadFieldRecogListName''' | |||
|STRING | |||
|Name of the output rfr list. | |||
|<code>RadFieldRecogListName="RadFieldRecogList.txt"</code><br> | |||
|RadFieldRecogList.txt | |||
|- | |||
|'''RadFiledRecogSigVecBaseName''' | |||
|STRING | |||
|Base name (beginning) of the significant vectors in the user database for RFR via comparator. DPE only searches such a files which includes this string and loads them into database. Only valid if RadFieldRecogDatabasePath is used. | |||
|<code>RadFiledRecogSigVecBaseName="vec"</code><br><code>RadFiledRecogSigVecBaseName="sig_vec"</code><br> | |||
|SigVec | |||
|- | |||
|'''RadFiledRecogSigVecEndData''' | |||
|STRING | |||
|End/extension/suffix of the data files of significant vectors in the user database for RFR via comparator. DPE only searches such a files which includes this string and loads them into database. Only valid if RadFieldRecogDatabasePath is used. | |||
|<code>RadFiledRecogSigVecEndData=".vec"</code><br><code>RadFiledRecogSigVecEndData=".data"</code><br><code>RadFiledRecogSigVecEndData=".txt"</code><br> | |||
|.vec | |||
|- | |||
|'''RadFiledRecogSigVecEndInfo''' | |||
|STRING | |||
|End/extension/suffix of the info files of significant vectors in the user database for RFR via comparator. DPE only searches such a files which includes this string and loads them into database. Only valid if RadFieldRecogDatabasePath is used. | |||
|<code>RadFiledRecogSigVecEndInfo=".info"</code><br><code>RadFiledRecogSigVecEndInfo=".vec_info"</code><br> | |||
|.vec_info | |||
|- | |||
|'''SamplingListBaseName''' | |||
|STRING | |||
|Base name of the sampling list which is used for txt and also for the json format. It is combined as this name and adding .txt and .json. | |||
|<code>SamplingListBaseName="SamplingList"</code><br><code>SamplingListBaseName="List"</code><br><code>SamplingListBaseName="Name"</code><br> | |||
|SamplingList | |||
|- | |||
|'''SamplingListJsonName''' | |||
|STRING | |||
|Name of output sampling list in json format. | |||
|<code>SamplingListJsonName="SampligList.json"</code><br> | |||
|SamplingList.json | |||
|- | |||
|'''SamplingListJsonPath''' | |||
|STRING | |||
|Path for output of sampling list in json format. | |||
|<code>SamplingListJsonPath="/Path/File JSON/Out/"</code><br> | |||
|current working directory | |||
|- | |||
|'''SamplingListName''' | |||
|STRING | |||
|Name of the output sampling list. | |||
|<code>SamplingListName="SamplingList.txt"</code><br> | |||
|SamplingList.txt | |||
|- | |||
|'''SensBias''' | |||
|FLOAT | |||
|Sensor Bias in volts. | |||
|<code>SensBias=-50</code><br><code>SensBias=200</code><br><code>SensBias=45.234</code><br> | |||
|50 | |||
|- | |||
|'''SensDens''' | |||
|FLOAT | |||
|Density of the sensor in g/cm3 for temperature at 20 deg. If not set then its values is deduced based on the type of material (only for those mentioned in SensMat). | |||
|<code>SensDens=2.93</code><br> | |||
|-1 (deduced based on detector settings) | |||
|- | |||
|'''SensHeight''' | |||
|INT | |||
|Height of the sensor in count of pixels. | |||
|<code>SensHeight=256</code><br><code>SensHeight=512</code><br><code>SensHeight=3</code><br> | |||
|256 | |||
|- | |||
|'''SensMat''' | |||
|STRING | |||
|Sensor material (Si, CdTe, GaAs). Example: | |||
|<code>SensMat="Si"</code><br> | |||
|Si | |||
|- | |||
|'''SensThick''' | |||
|FLOAT | |||
|Sensor thickness in micrometers. | |||
|<code>SensThick=300</code><br><code>SensThick=500</code><br> | |||
|300 | |||
|- | |||
|'''SensWidth''' | |||
|INT | |||
|Width of the sensor in count of pixels. | |||
|<code>SensWidth=256</code><br><code>SensWidth=512</code><br><code>SensWidth=3</code><br> | |||
|256 | |||
|- | |||
|'''SigVecName''' | |||
|STRING | |||
|Name of the file with significant vector configuration. | |||
|<code>SigVecName="SigVec.ini"</code><br> | |||
|empty string | |||
|- | |||
|'''SigVecPath''' | |||
|STRING | |||
|Path to the file with significant vector configuration. | |||
|<code>SigVecPath="/Path/To/SigVec/Config/"</code><br> | |||
|empty string | |||
|- | |||
|'''TimeAcq''' | |||
|FLOAT | |||
|Measurement/data acquisition time which can be set to use it instead of elapsed time evaluated in engine run. It is given in seconds. | |||
|<code>TimeAcq=3</code><br><code>TimeAcq=1e2</code><br> | |||
|-1 (unused) | |||
|- | |||
|'''TimeCoinc''' | |||
|FLOAT | |||
|Time in nanoseconds for evaluations coincidence events -> clusters whose times are in between this this interval are grouped as coincidence group (first cluster time is down edge and plus coince group time is upper edge). | |||
|<code>TimeCoinc=1e2</code><br><code>TimeCoinc=10</code><br><code>TimeCoinc=1.2554</code><br> | |||
|100 | |||
|- | |||
|'''TimeFrameAcq''' | |||
|FLOAT | |||
|Acquisition time of frame (if frame are processed). DPE also tries to make an estimation if the information is included in the input data (clog etc.). | |||
|<code>TimeFrameAcq=1e-3</code><br><code>TimeFrameAcq=1</code><br> | |||
|0 (deduce from data) | |||
|- | |||
|'''TimeFrameDead''' | |||
|FLOAT | |||
|Dead time of frame (if frame are processed). DPE also tries to make an estimation if the information is included in the input data (clog etc.). | |||
|<code>TimeFrameDead=1e-3</code><br><code>TimeFrameDead=2</code><br> | |||
|0 (deduce from data) | |||
|- | |||
|'''TimeSampling''' | |||
|FLOAT | |||
|Time in seconds for individual evaluations of the data sample -> the whole data sample and its time of collection is divided into sampling intervals with duration of the TimeSampling. | |||
|<code>TimeSampling=1e-6</code><br><code>TimeSampling=0.00001</code><br><code>TimeSampling=233.2</code><br> | |||
|1 | |||
|- | |||
|'''DoMask''' | |||
|BOOL | |||
|Control to turn off masking if mask file is given. | |||
|<code>DoMask=true</code><br><code>DoMask=false</code><br> | |||
|false | |||
|- | |||
|} | |||
== Detector Settings== | |||
It is possible to set the detector configuration with following list of parameters in the main parameters file: | |||
< | * '''SensMat''' - sensor material (Si, CdTe, GaAs etc). Example: <code>SensMat = "Si"</code> to specify that the sensor material is silicon. Possible values: | ||
** <code>Si</code> - for silicon sensor. Used density: 2.329 g/cm-3. | |||
** <code>CdTe</code> - for cadmium teluride sensor. Used density: 5.85 g/cm-3. | |||
** <code>GaAs</code> - for galium arsenide. Used density: 5.318 g/cm-3. | |||
* '''SensThick''' - Sensor thickness in micrometers. Example: <code>SensThick = 500</code> to specify the thickness of the detector as 500 micrometers. This is continues variable. | |||
* '''SensBias''' - sensor bias in volts. Example: <code>SensBias = 100</code> to specify that applied was 100 V. This is a continues variable. | |||
* '''ChipType''' - type of the chip (TPX or TPX3). Example: <code>ChipType = "TPX3"</code> to specify that TPX3 was used chip. Possible values: | |||
** <code>TPX</code> - for timepix chip | |||
** <code>TPX3</code> - for timepix 3 chip | |||
** <code>TPX2</code> - for timepix 2 chip | |||
* '''SensDens''' - set density of used material in g/cm-3. DPE has default values for material specified in <code>SensMat</code>. | |||
These values are used, for example, in the calculation of dose rate because there is a dependence on the sensor material and thickness. It also determines which PID algorithm should be used. | |||
== Input== | |||
<span id="single-file-processing"></span> | |||
=== Single File Processing === | |||
To process specific file with DPE, its name <code>FILE_IN_NAME</code> and path <code>/PATH/TO /FILE/</code> must be specified via following parameters in parameters file: | |||
== | <pre> FileInName = "FILE_IN_NAME" | ||
FIleInPath = "/PATH/TO /FILE/"</pre> | |||
This will cause that DPE will try to find this file and process it. Error is produced if the file can not be found and the processing is aborted. | |||
<span id=" | <span id="supported-file-formats"></span> | ||
=== | === Supported File Formats === | ||
The current version of the DPE engine is capable to process following files: | |||
* Data driven files - t3pa, t3p, t3r | |||
* Cluster log files- calibrated/uncalibrated | |||
* Elist - only so-called full elist which includes also the header with cluster variables name. | |||
* Frame formats | |||
If a clog from itot+count measurement is used then it is needed to specify this in the main config file with following option: | |||
<pre> MeasMode = "itot+count"</pre> | |||
Supported suffixes are following: | |||
<pre> ".t3pa", ".t3p", ".t3r", ".txt", ".pmf", ".plog", ".bmf", ".clog", ".elist", ".advelist", ".extelist"</pre> | |||
Currently, the only toa formats are unsupported with undefined results deom DPE (time is handled as energy information).<br /> | |||
Examples of most of the suported files can be found in the <code>Examples/Test</code> directory (list of examples is in the section <code>First Example of Run</code>) | |||
<span id=" | <span id="multi-filebatch-processing"></span> | ||
=== | === Multi File/Batch Processing === | ||
It is possible to process several files in one run of DPE. There are several ways how to specify which files should be processed: * Process whole directory specifying path to the directory. * Process given files specifying their names in <code>FileInName</code> as an array. * Process only those files in directory which includes given strings in their names. Examples of such a processings can be found below. | |||
* | |||
* | |||
* | |||
Lets assume that our directory at path: <code>/PATH/TO/DIR/</code> includes the following files: * File_001.t3pa * File_002.t3pa * File_003.t3pa * File.t3pa * File_2.clog | |||
In the first case, when a whole directory should be processed, only the path to the directory is specified in the parameters file: | |||
The | <pre> FileInPath = "`/PATH/TO/DIR/"</pre> | ||
DPE checks the existence of the directory and finds all files. It processes only those which have the same suffix/ending, and DPE found them first.<br /> | |||
For this specific directory, let’s assume that <code>File_001.t3pa</code> was the first found by DPE, then all files ending with <code>.t3pa</code> are processed.<br /> | |||
The file <code>File_2.clog</code> is not processed, because of the different suffix.<br /> | |||
In the second case, when names of files are specified, it is needed to specify the names into parameter <code>FileInName</code>: | |||
<pre> FileInName = "File_001.t3pa", "File_003.t3pa", "File.t3pa"</pre> | |||
DPE will process files: <code>File_001.t3pa</code>, <code>File_003.t3pa</code> and <code>File.t3pa</code>. '''Files with different suffixes should not be mixed!''' (e.g. clog and t3pa, DPE handles these files differently, which might cause unexpected behavior of the engine.)<br /> | |||
In the last case, when user wants to specify only parts of the names, it is possible to specify beginning and ending of the name: | |||
<pre> FileInName = "File_" | |||
FileInNameEnd = "t3pa"</pre> | |||
This will cause that only the following files are processed: <code>File_001.t3pa</code>, <code>File_002.t3pa</code> and <code>File_003.t3pa</code>. If only the first part of the name should be specified, there is need for one additional parameter <code>DoMultiFile</code>: | |||
<pre> | <pre> FileInName = "File_0" | ||
DoMultiFile = true</pre> | |||
DPE will again process the same list of files (additional 0 is needed due to the <code>File_2.clog</code>). In the opposite case when only ending of the name is specified, there is '''no''' need to add these parameters: | |||
<pre> FileInNameEnd = ".t3pa"</pre> | |||
Files <code>File.t3pa</code>, <code>File_001.t3pa</code>, <code>File_002.t3pa</code> and <code>File_003.t3pa</code> will be processed. | |||
The switch <code>DoMultiFile</code> can also be used to suppress multi file processing, setting it to false. | |||
The | |||
== Output== | |||
== | <span id="log-and-print"></span> | ||
=== Log and Print === | |||
During a run of the program, there is an output into the terminal. This include some basic information about the DPE settings and its run.<br /> | |||
It can be turned off with <code>DoPrint = false</code> (default is true). Copy of this output is also given into log file, usually (default) called <code>LogFile.txt</code>. It can be turned off with <code>DoLog = false</code> (default is true). The preprocessing run of clusterer also created its own log file called <code>log.txt</code> where minimal log of clusterer run is stored. | |||
<span id="directory-tree"></span> | |||
=== Directory Tree === | |||
Results of DPE are exported into several directories separated in most cases based on the type of the analysis. There are some general directories for overall results. The creation of directories depends on used analysis and some of them might not be created if given processing is disabled or unsupported for given settings. | |||
* '''File''' - main data files: elist, cluster log and sampling list. It is always created. | |||
* '''Graph''' - graphical representation of the sampling list. Only if <code>DoExportGraphics = true</code>. | |||
* '''Hist''' - histograms of cluster variables. | |||
* '''SpatialMap''' - spatial maps of cluster variables also with respect to the PID classes. | |||
* '''SigVec''' - significant vector which uniquely specifies the given radiation field. | |||
* '''EventVisual''' - individual cluster plots, integrated sensor plots and plots of clusters for individual PID classes. | |||
* '''Direction''' - results of directional analysis: estimation of radiation direction. | |||
* '''CoincEvent''' - results of coincidence analysis: evaluates time correlation between clusters. | |||
* '''Compton''' - results of coincidence analysis: evaluates time correlation between clusters. | |||
* '''RadFieldRecog''' - results of coincidence analysis: evaluates time correlation between clusters. | |||
These directories can have inner structure based on the specific exports of given analysis (export of graphs in directional analysis is done into directory <code>Direction/Graph</code>). More details can be found in the sections dedicated to these specifics analysis: directional, Compton, coincidence and radiation field recognition.<br /> | |||
The whole export/output can be turned off if <code>DoExportGraphics = false</code> and <code>DoExportText = false</code>. | |||
<span id="text-basic-output"></span> | |||
=== Text basic output === | |||
* '''Elist''' - file with cluster variables/parameters (this file is in default suppressed with <code>DoRemoveElist = true</code> and replaced with extended elist below). | |||
* '''Extended Elist''' - elist extended with additional columns of PID class (-1 to N-1 of classes where -1 is for others, <code>PIDClass</code>) and filter pass (1 = passed, 0 = not passed, <code>FilterPass</code>) if the filter is used. | |||
* '''Cluster log/clog''' - detailed list of clusters containing pixelated information. Each line which is starting with <code>[</code> includes pixels of a cluster ([X,Y,E,T] = [X position,Y position,energy,time if tpx3 is sued] of a pixel). This file includes '''all clusters''' and the mask or filters are not accounted for. | |||
* '''Sampling list''' - it includes basic information about the radiation field with respect to the observables and their time dependencies. | |||
* '''Histograms''' - histograms of given variables - total and also for each PID class. | |||
* '''Spatial maps''' - matrices of the sensor filled with integrated cluster variables (energy, size, LET, E/S). | |||
* '''Significant vector''' - unique vectors characterizing processed radiation field computed based on histograms. | |||
* '''Event and sensor visualizations''' - data files with matrices of exported clusters and sensor. | |||
Each individual analysis also produces some parts of these results and more detailed descriptions are offered in dedicated sections below. | |||
<span id="graphical-output"></span> | |||
=== Graphical output === | |||
The current version uses a python scripts to convert the text output into graphics: | |||
* Histograms | |||
* Spatial maps | |||
* Individual cluster plots (with values of the cluster variables) and integrated sensor plots | |||
* Graphs of time evolution of physical products | |||
* Directional maps | |||
* Compton directional reconstruction | |||
* … | |||
For this purposes, the matplolib library is used. It can be slow for some configuration (for example high number of bins). To speed up the export of graphics, it is possible to set shown number of bins with following option: | |||
<pre> Hist1DGraphicsRebin = 100</pre> | |||
which will cause that each histogram will be rebinned to show maximally 100 bins but the exported text files are still with the original number of bins. This option is only functional for 1D histograms. | |||
Another possibility to speed up the plotting is the option <code>DoGraphicsMultiThread</code> which exploits several CPU threads for graphics creation. This option is by default on/true. The program then uses all threads on the given PC. To disable this feature use: | |||
The | <pre> DoGraphicsMultiThread = false</pre> | ||
The multi processing can use up to several GB of ram memory which can be tuned with parameter <code>GraphicsMultiThreadFileSize</code>. Its default value is 3e6 and lower this number if the ram usage should be decreased (minimal ram usage is defined with needed ram for the biggest plot/histogram which can be of order of GBs). | |||
To produce given plots, a directory <code>.temp</code> is created in the working directory. All python scripts are saved there before processing.<br /> | |||
If the program is terminated just before or in the middle of graphics creation (python run) then a next run will firstly evaluate unfinished plotting and then proceeds with new plots therefore the whole plotting is delayed. In the case of sudden program abortion, it is recommended to check for existence of the <code>.temp</code> directory and delete it before the next run of DPE. | |||
== Clusterer and Preprocessing== | |||
The main purpose of this stage is to provide conversion of '''pixelated data''' into groups of correlated pixels, '''clusters'''. The level of correlation depends on the detector used for data collection. In the case of TPX detector, only spatial correlation is given and it is afterwards used for '''clusterization process'''. On contrary for the TPX3, time information is provided along the spatial one which gives opportunity to proceed with more precise clusterization process and eventually avoid additional unwanted effects as pile-up. At this stage main '''energy per pixel''' '''calibration''' is applied to convert digital ToT information to energy in keV. This process is accompanied with additional calibrations and corrections which aim to maximally supress detector unwanted or anomalous behaviour. The clusterization stage is followed with '''cluster analysis''' which evaluates clusters '''morphological and spectral features'''. | |||
To summerize the pre-processeing stage: | |||
'' | * '''Clustering''' - grouping pixels based on coordinate and also on time information if given | ||
* '''Calibration and corrections''' - energy calibration, cluster size correction, XRF peak suppression, TPX high energy correction etc. | |||
* '''Cluster parameters''' - calculation of overall cluster parameters (total energy, roundness etc.) | |||
[[File:DPE PreProc.png|frameless|1157x1157px]] | |||
<span id="clustering"></span> | |||
=== Clustering === | |||
The | The clusterisation converts a frame or a continual stream of received pixels to sets of pixels which belonging to individual particles, clusters. A distinction is needed between two read out modes. In the case of the Timepix, the detector is only capable to measure in the frame mode and obtained data format is in a form of frames. The information is encoded into a sensor matrix of only hit pixels where each of them contain energy deposited during the acquisition time (if a measurement in ToT is performed[[DPE#sdfootnote1sym|<sup>1</sup>]]). The same mode can be also used with the Timepix3 detector. The advantage with respect to the Timepix is the possibility to measure in the data driven mode. The final data output in this mode resembles the continual stream of pixels where time of arrival is stored together with deposited energy[[DPE#sdfootnote2sym|<sup>2</sup>]]. This significantly improves the ability for the restoration of the individual particle information. | ||
The algorithm used in the case of the '''frame mode''' relies only on the coordinate information. Pixels in a frame which are coordinate neighbours (lies within 8-pixel surroundings) are grouped together and called clusters: | |||
It is also possible to measure in ToA mode, time of arrival or in a counting event mode. | |||
In the case of combined ToT and ToA measurement which is another advantage of the Timepix3 with respect to the Timepix detector. | |||
''eq with clustering condition'' | |||
During this clusterization process, an additional information about the readout can be used. It is known that the chip matrix of pixels is read out column by column which means that the stream of pixels can be efficiently parsed based on the column order. There are several disadvantages and crucial points concerning frame read-out: | |||
# '''Pile-up:''' All particles within one frame whose pixels are connecting cannot be separated within the clusterization process and they create the pile-up effect. | |||
# '''Loss of uncollected charge after the end of frame:''' It is possible that not all charge has been collected till the end of a frame (mostly some pixels with large deposited energy). | |||
# '''Dead time for matrix read out:''' This approximately several of millisecond (more than 10 ms). | |||
These issues and disadvantages motivated development of the '''data-driven''' '''read-out mode'''. The coordinate information is accompanied by time of arrival which allows to reconstruct events also based on the time information. In this case, the stream of pixels is no longer separated into frames but it is continuously read by the read-out hardware and software. | |||
In the first step, the stream has to time ordered because from the nature of the read-out, it is possible that the pixel stream violates time order. In the next stage, a block of pixels is taken from the full stream/list of pixels based on a time condition: | |||
''eq'' | |||
The parameter should reflect expected time windows in the pixel stream and therefore is disproportional to particle fluence. The value is of order of . In those cases when the condition cannot be fulfilled, a fixed block of of pixels is taken. This number is approximately 100 000 pixels which should statistically minimized the error originating from possible separation of pixels belonging to one particle. | |||
The clusterization process itself separates pixels into clusters based on two main conditions: | |||
''eq'' | |||
These conditions can be translated as a need for pixels to be coordinate neighbours and to have maximal difference in time equal or less than time . The value of this time parameter is of order of tens of nanoseconds and it is dependent on the charge sharing effect and its magnitude. | |||
=== Calibration and corrections === | |||
==== Energy calibration ==== | |||
[[File:Energy calib.png|Energy calibration of Time-Over-Threshold. Reprinted from Jakubek 2011.|376x376px|alt=Energy calibration of Time-Over-Threshold. Reprinted from jakubek 2011.|thumb]]After receiving the information from the detector, the possible energy and time of the detector are in corresponding ToT and ToA counts. These variables have to be calibrated and converted to obtain results in energy and time. Time conversion is done with following formula: | |||
''eq of toa conversion'' | |||
The first part of the calibration is the same for the Timepix and Timepix3 ASICs. This work was published by Jakubek<ref>J. Jakubek, ‘Precise energy calibration of pixel detector working in time-over-threshold mode’, ''Nucl. Instrum. Methods Phys. Res. Sect. Accel. Spectrometers Detect. Assoc. Equip.'', vol. 633, pp. S262–S266, May 2011, doi: 10.1016/j.nima.2010.06.183.</ref> and it is based on a combination of a linear and rational calibration function: | |||
''eq of calibration''where ''a,b,c'' and ''t'' are parameters obtained from the calibration procedure. | |||
=== | ==== High Energy TPX Calibration ==== | ||
[[File:EnergyHist1d OFF ON HECorr.png|Histogram of measured energy before (light) and after (dark) application of the high energy calibration. Different biases were applied during the measurement to observe increase in the mean per pixel energy which has to be mainly corrected for (higher bias -> higher per pixel energy).|360x360px|alt=High energy calibration for TPX applied on measured data from Am241 aplha particles. Shaded lines are without HE calibration, solid lines are with HE calibration.|thumb]] | |||
The calibration function has maximum of its validity up to approximately 700 keV for TPX. An additional per pixel calibration process has to be utilized for measurements in which higher per pixels energies are obtained. The additional calibration was intruced as a '''global model''' based on the results from the article of Sommer et al <ref>M. Sommer, C. Granja, S. Kodaira, and O. Ploc, ‘High-energy per-pixel calibration of timepix pixel detector with laboratory alpha source’, ''Nucl. Instrum. Methods Phys. Res. Sect. Accel. Spectrometers Detect. Assoc. Equip.'', vol. 1022, p. 165957, Jan. 2022, doi: 10.1016/j.nima.2021.165957.</ref>. Having the coefficients of the standard calibration and correlation matrix between these standards and the extended ones, it is possible to introduce an estimation of their values. | |||
< | To use this additional calibration:<syntaxhighlight lang="cpp"> | ||
DoTpxHighEnergyCorr = true | |||
</syntaxhighlight>The resulted histogram of the measured energy can be found in the fugere on the left. Two cases are shown in this image, before (light) and after (dark) an application of the high energy correction. Change in the sensor bias was also introduce during the measurement to demonstrate a different/increasing value of mean per pixel energy. Higher bias in the sensor decreases the effect of charge sharing and increases the per pixel energy and mainly the maximal per pixel energy of cluster. Therefore, it can be seen that in the most disturbed case of 200V the applied correction restores the original information to be in accordance with low bias measurement and their corrected versions. | |||
In the case of TPX3, the energy range is limited with size of the ToT counter which has nominal value 1022. This means that the pixels cannot count more than 1022 counts in contrast to the TPX with counter range of 11810. From the practical point of view, the maximal value which can be measured in one pixel is between 450-500 keV (based on the settings of DACs). The standard energy calibration covers this region and a distortion appears for energy above 2 MeV as the volcano effect. | |||
==== Cluster Energy Size Based Correction ==== | |||
... | |||
==== Halo Effect ==== | |||
... | |||
==== Correction for XRF peak suppression ==== | |||
... | |||
'''Pile-up Recognition''' | |||
... | |||
==== Volcano Effect ==== | |||
[ | [[File:Volcano effect.png|thumb|''Heavy ion cluster in energy plot with a direction almost perpendicular to the sensor plane. Two additional sensor effects can be observed: the volcano and the halo effect.'']] | ||
... | |||
Source<ref>B. Hartmann ''et al.'', ‘Distortion of the per-pixel signal in the Timepix detector observed in high energy carbon ion beams’, ''J. Instrum.'', vol. 9, no. 09, pp. P09006–P09006, Sep. 2014, doi: 10.1088/1748-0221/9/09/P09006.</ref> | |||
== | ==== Overshoot ==== | ||
... | |||
==== Time-walk Correction ==== | |||
< | ... | ||
Source<ref>D. Turecek, J. Jakubek, and P. Soukup, ‘USB 3.0 readout and time-walk correction method for Timepix3 detector’, ''J. Instrum.'', vol. 11, no. 12, pp. C12065–C12065, Dec. 2016, doi: 10.1088/1748-0221/11/12/C12065.</ref> | |||
=== Cluster Variables and Elist === | |||
<span id="SecClusterVarElist"></span> | |||
* The file includes clusters/events represented as a set of cluster variables. | |||
* The significant feature of the Elist is a presence of coincidence group expressed as coincidence number and it is adjustable based on the coincidence time interval. | |||
{| class="wikitable sortable" | |||
!LengthCorrStd | |||
!Unit | |||
!Column | |||
!Description | |||
< | |- | ||
|<code>DetectorID</code> | |||
| - | |||
|1 | |||
|It serves to as a unique ID in the cases when multi detector setup is processed. | |||
|- | |||
< | |<code>EventID</code> | ||
| - | |||
|2 | |||
|ID of events. In the case that several clusters are in the coincidence, then they have the same EventID. It starts from 0 and clusters which are not in coincidences are also accounted in this ID. | |||
|- | |||
|<code>X</code> | |||
|px | |||
|3 | |||
|X coordinate of cluster based on the weighted mean value of X coordinate of pixels, weighted with the energy of pixels. Value ha range from 0.5 to 255.5 (minimal and maximal value of X of pixels). | |||
|- | |||
|<code>Y</code> | |||
|px | |||
|4 | |||
|Y coordinate of cluster based on the weighted mean value of Y coordinate of pixels, weighted with the energy of pixels. Value ha range from 0.5 to 255.5 (minimal and maximal value of Y of pixels). | |||
|- | |||
|<code>E</code> | |||
|keV | |||
|5 | |||
|Energy of cluster, also called volume. It is sum of energies of all cluster pixels. | |||
|- | |||
|<code>T</code> | |||
|ns | |||
|6 | |||
|Time. It is the minimal time of cluster pixels. | |||
|- | |||
|<code>Flags</code> | |||
| - | |||
|7 | |||
|Free column for user needs. In the case of frame data, it includes frame number. | |||
|- | |||
|<code>Size</code> | |||
|px | |||
|8 | |||
|Count of cluster pixels. | |||
|- | |||
|<code>Height</code> | |||
|keV | |||
|9 | |||
|Maximal value of energy of cluster pixels. | |||
|- | |||
|<code>BorderPixCount</code> | |||
|px | |||
|10 | |||
|Count of border pixels in cluster. | |||
|- | |||
|<code>Roundness</code> | |||
| - | |||
|11 | |||
|Morphological feature expressing similarity of cluster to round shape. | |||
|- | |||
|<code>AngleAzim</code> | |||
|deg | |||
|12 | |||
|Estimation of | |||
|- | |||
|<code>Linearity</code> | |||
| - | |||
|13 | |||
|Morphological feature expressing similarity of cluster to linear shape. | |||
|- | |||
|<code>LengthProj</code> | |||
|px | |||
|14 | |||
|Length of cluster based on maximal distance between pixels after projection to the cluster axis. | |||
|- | |||
|<code>WidthProj</code> | |||
|px | |||
|15 | |||
|Length of cluster based on maximal distance between pixels after projection to the axis perpendicular to the cluster axis. | |||
|- | |||
|<code>IsSensEdge</code> | |||
| - | |||
|16 | |||
|Information whether cluster is at sensor edge. 0 for flase, 1 for true (is at sensor edge). | |||
|- | |||
|<code>StdAlong</code> | |||
|px | |||
|17 | |||
|Weighted standard deviation of pixels with respect to cluster axis weighted with pixels energy. | |||
|- | |||
|<code>StdPerp</code> | |||
|px | |||
|18 | |||
|Weighted standard deviation of pixels with respect to axis perpendicular to the cluster axis and weighted with pixels energy. | |||
|- | |||
|<code>Thin</code> | |||
| - | |||
|19 | |||
|Morphological feature expressing thinness of cluster. | |||
|- | |||
|<code>Thick</code> | |||
| - | |||
|20 | |||
|Morphological feature expressing thickness of cluster. | |||
|- | |||
|<code>CurlyThin</code> | |||
| - | |||
|21 | |||
|Morphological feature expressing combination of cluster thinness and inverse of linearity. | |||
|- | |||
|<code>EpixMean</code> | |||
|keV | |||
|22 | |||
|Mean value of pixels energy in cluster. | |||
|- | |||
|<code>EpixStd</code> | |||
|keV | |||
|23 | |||
|Standard deviation of pixels energy in cluster. | |||
|- | |||
|<code>LengthCorrStd</code> | |||
|px | |||
|24 | |||
|Cluster length in sensor plane corrected for charge sharing based on weighted standard deviation. | |||
|- | |||
|<code>Length3DCorrStd</code> | |||
|um | |||
|25 | |||
|Cluster length in sensor volume corrected for charge sharing based on weighted standard deviation. | |||
|- | |||
|<code>LengthCorrWidth</code> | |||
|px | |||
|24 | |||
|Cluster length in sensor plane corrected for charge sharing based on cluster width. | |||
|- | |||
|<code>Length3DCorrWidth</code> | |||
|um | |||
|25 | |||
|Cluster length in sensor volume corrected for charge sharing based on cluster width. | |||
|- | |||
|<code>LengthCorrNabha</code> | |||
|px | |||
|24 | |||
|Cluster length in sensor plane corrected for charge sharing based on Nabha model. | |||
|- | |||
|<code>Length3DCorrNabha</code> | |||
|um | |||
|25 | |||
|Cluster length in sensor volume corrected for charge sharing based on Nabha model. | |||
|- | |||
|<code>AngleElev</code> | |||
|deg | |||
|26 | |||
|Elevation angle of cluster. | |||
|- | |||
|<code>LET</code> | |||
|keV/um | |||
|27 | |||
|Linear energy transfer based on the corrected length and energy (E). | |||
|- | |||
|<code>Diameter</code> | |||
|px | |||
|28 | |||
|... | |||
|- | |||
|<code>PIDClass</code> | |||
| - | |||
|29 | |||
|If PID is used, this number expressed class into which the cluster was classified. It goes from 0 to count_class - 1. | |||
|- | |||
|<code>FilterPass</code> | |||
| - | |||
|30 | |||
|If filter is used, then it includes information whether cluster passed given filters or did not. 1 for passed, 0 for did NOT pass. | |||
|} | |||
Example of created elist: | |||
[[File:Elist text.png|1213x1213px]] | |||
=== Cluster log === | |||
Clusters can be exported as individuals pixels. This format is called cluster log or clog and it is highly dependent on the input data (detector, measured mode, etc.). In the case of '''TPX''' detector the format is following: | |||
<pre>Frame 1 (UNIX_TIME, ACQ_TIME s) | |||
[X_1, Y_1, iToT_1/E_1/] [X_2, Y_2, iToT_2/E_2/]... | |||
.... | |||
Frame 2 (UNIX_TIME, ACQ_TIME s) | |||
[X_1, Y_1, iToT_1/E_1] [X_2, Y_2, iToT_2/E_2]... | |||
....</pre> | |||
where <code>ACQ_TIME</code> is acquisition time of frames and <code>UNIX_TIME</code> is current UNIX time stamp in seconds (both). Every line starting with <code>[</code> is one cluster and in each square brackets is one pixel and its information <code>[x coordinate,y coordinate, iToT/Count/Energy]</code>. '''The ToA and Count modes are not correctly implemented and they are treated as the ToT mode (especially of importance for the cluster parameters in the elist).''' | |||
The format is different for the '''TPX3''' detector: | |||
<pre>Frame 1 (FIRST_TIME_COINC_GROUP, COINC_TIME_WINDOW ns) | |||
[X_1 ,Y_1, ToT_1/E_1, T_Diff_1] [X_2, Y_2, ToT_2/E_2, T_Diff_2]... | |||
[X_1 ,Y_1, ToT_1/E_1, T_Diff_1] [X_2, Y_2, ToT_2/E_2, T_Diff_2]... | |||
.... | |||
Frame 2 (FIRST_TIME_COINC_GROUP, COINC_TIME_WINDOW ns) | |||
[X_1 ,Y_1, ToT_1/E_1, T_Diff_1] [X_2, Y_2, ToT_2/E_2, T_Diff_2]... | |||
[X_1 ,Y_1, ToT_1/E_1, T_Diff_1] [X_2, Y_2, ToT_2/E_2, T_Diff_2]... | |||
....</pre> | |||
where each frame in this case coincidence group = clusters whose min times are within the <code>COINC_TIME_WINDOW</code> starting with the first one at <code>FIRST_TIME_COINC_GROUP</code>. The pixels time is written as difference with the <code>FIRST_TIME_COINC_GROUP</code>: <code>T_Diff = T_Pix - FIRST_TIME_COINC_GROUP</code> where <code>T_Pix</code> is original time of pixels in ns. An exception from this rule is measurement in the frame mode '''iToT+Count''' when data are already saved in cluster log within the Pixet. To properly evaluate this data with the clusterer, it is needed to add additional option into the command line: | |||
<pre>--measmode "itot+count"</pre> | |||
This data does not include additional ToA information therefore the format has the same meaning as in the first case with TPX: | |||
<pre>Frame 1 (UNIX_TIME, ACQ_TIME s) | |||
[X_1 ,Y_1, iToT_1/E_1, Count_1] [X_2, Y_2, iToT_2/E_2, Count_2]... | |||
[X_1 ,Y_1, iToT_1/E_1, Count_1] [X_2, Y_2, iToT_2/E_2, Count_2]... | |||
.... | |||
Frame 2 (UNIX_TIME, ACQ_TIME s) | |||
[X_1 ,Y_1, iToT_1/E_1, Count_1] [X_2, Y_2, iToT_2/E_2, Count_2]... | |||
[X_1 ,Y_1, iToT_1/E_1, Count_1] [X_2, Y_2, iToT_2/E_2, Count_2]... | |||
....</pre> | |||
Summary of supported input data formats: * TPX3, data driven mode, tot+toa * TPX3, frame mode, itot+count * TPX, frame mode, tot Every other format will be processed but the results might be uncorrected evaluated, for example in the cluster parameters. | |||
Reprocessing of cluster logs can create cluster logs which could have incorrect coincidence group (in the case of TPX3 with ToA). If the coincidence window <code>el-coinctoadiff</code> is '''increased''' with respect to the original one then the resulted group will not integrate several groups together which would fulfill this new time condition. | |||
== | == Masking== | ||
There is a possibility to mask a part of the sensor or individual pixels in the case of t3pa and clog files.<br /> | |||
These pixels are omitted in the pre-processing.<br /> | |||
This can be used to avoid a signal from noisy pixels or to focus on some more interesting part of the sensor.<br /> | |||
The configuration is done through a configuration file. An example can be found in the program directory.<br /> | |||
If the masking is used the during the pre-processing an additional t3pa/clog file is created with masked pixels according to the user configuration in the export directory (starts with <code>MASK_+ FileInName</code>). | |||
<span id="inclusion-of- | <span id="inclusion-of-mask.ini-into-parametrs-file"></span> | ||
=== Inclusion of | === Inclusion of Mask.ini into Parametrs File === | ||
The configuration file of the | The configuration file of the mask can be included with following options in the ParametersFile: | ||
<pre> | <pre> MaskName = "Mask.txt" // Name of the INI configuration file of the mask. | ||
MaskPath = "PATP/TO/MASK/" // Path of the INI configuration file of the mask.</pre> | |||
<span id="syntax-of-configuration-file | <span id="syntax-of-configuration-file"></span> | ||
=== Syntax of Configuration File === | === Syntax of Configuration File === | ||
The | * The pixel which should be masked is written in format [X,Y]. These are X and Y coordinates of the pixel where both of them are integers from 0 to 255. The sensor orientation is usual and X coordinate is for rows and Y is for lines. They can be written in lines or rows in the mask file but always the [X,Y] in one line (can not be separated over several lines) | ||
* If a larger part is of interest, following notation is used: [X_min - X_max, Y_min - Y_max] where X_min is the minimum coordinates of a pixel, X_max is the maximum coordinate of a pixel etc. For example, [0-255,0-120] masks almost half of the sensor - on the X axis from 0 to 255 and on the Y axis from 0 to 120. A simpler demonstration can be masking of one line of pixels (11th) : [0-255, 10]. | |||
Example of all possibilities for masking: | |||
<pre> [0,1] [2,5] | |||
[2,3] | |||
[2-43,5] | |||
[76,8-90] | |||
[0-50,50-60]</pre> | |||
<span id="masking-with-pixet-mask"></span> | |||
=== Masking with Pixet Mask === | |||
< | It is also possible to exploit the masking created in Pixet, which takes the form of a text file where the lines and rows correspond to those of the detector. The orientation in this case is the same for x axis, but it is reversed for y axis where bottom of the detector is first line and top part of the detector is last line. The delimiter in this case is simple white space <code></code>. The masked pixels are marked as <code>0</code> and unmasked are <code>1</code>. See following example (<code>...</code> used as abbreviations): | ||
<pre> | <pre> 1 1 1 1 0 1 ... 1 | ||
1 1 1 1 1 1 ... 1 | |||
... | |||
0 1 1 1 1 1 ... 1</pre> | |||
<span id=" | <span id="statistical-information"></span> | ||
=== | === Statistical Information === | ||
Statistical information can be found in the general sampling list and it is also printed into log.<br /> | |||
Example for the sampling list in json format (explanations are given in the comments after <code>//</code>): | |||
<syntaxhighlight lang=" | <syntaxhighlight lang="cpp"> "CountPixHit_Sum_MaskOk_cnt":242461, // Count of pixels which passed the mask and used for processing. | ||
"CountPixHit_Sum_MaskOk_perc":99.23, // The same as above but as part/percentage to total amount of pixels. | |||
"CountPixHit_Sum_MaskOmit_cnt":1872, // Count of pixels which did NOT pass the mask and they are used for processing. | |||
"CountPixHit_Sum_MaskOmit_perc":0.766, // The same as above but as part/percentage to total amount of pixels. | |||
"CountPixHit_Sum_All_cnt":41935, // Total count of loaded pixels for raw data.</syntaxhighlight> | |||
Example for log file: | |||
< | <pre> Mask conf file - name: Mask.txt | ||
Mask conf file - path: ./Test/data/test_029/ | |||
... | |||
-------------------------------------------------- | |||
SENSOR MASK | |||
...has been loaded. | |||
Count of masked pixels: 901 (1.37%) | |||
-------------------------------------------------- | |||
... | |||
Processing RawData T3PA: | |||
Creating masked raw file: MASK_tot_toa.t3pa | |||
... | ... | ||
-------------------------------------------------- | -------------------------------------------------- | ||
MASK | |||
-------------------------------------------------- | |||
CountPixHit_Sum_MaskOk [-]: 242461 | |||
CountPixHit_Sum_MaskOk [%]: 99.23 | |||
CountPixHit_Sum_MaskOmit [-]: 1872 | |||
CountPixHit_Sum_MaskOmit [%]: 0.766 | |||
CountPixHit_Sum_All [-]: 41935</pre> | |||
First part shows from which destination is the mask taken and what is the name of the file.<br /> | |||
The second part informs about successful loading of the file and how many pixels are masked.<br /> | |||
The third part appears in the preprocessing sections and informs about actual masking of the raw data.<br /> | |||
The last part is in the sections with results and includes the same information as it is in the sampling file. | |||
== Filtering== | |||
During the processing, filters can be used to obtain only information about particles of interest.<br /> | |||
The filters are applied on cluster variables/parameters level (e.g. energy, height etc.). It is specified trough a configuration file in the INI format.<br /> | |||
The cluster variables which should be used for filtering are specified with their unique name which is included<br /> | |||
in the header of the created elist (if elist is input then it has to be already part of the file).<br /> | |||
All results produced by DPE are only for filtered particles/for those which passed filter (histograms, graphs, spatial maps etc.).<br /> | |||
An example can be found in the program directory. | |||
<span id="inclusion-of-filter.ini-into-parameters-file"></span> | |||
=== Inclusion of Filter.ini into Parameters File === | |||
The configuration file of the filter can be included with following options in the ParametersFile: | |||
<pre> FilterName = "Filter.ini" // Name of the INI configuration file of the filter. | |||
FilterPath = "/PATH/TO/FILTER/" // Path of the INI configuration file of the filter.</pre> | |||
<span id="syntax-of-configuration-file-1"></span> | |||
=== Syntax of Configuration File === | |||
The configuration file in INI format is a set of sections where each section is dedicated to one cluster parameter for which filtering should be done.<br /> | |||
Individual parts of these sections are specific values of the filters for the given cluster parameter. Example with filter conditions for energy of clusters alias row <code>E</code> in elist: | |||
* One section of the configuration file is named based on the key <code>E</code> in the elist: <code>[E]</code>. | |||
* Condition is then written as <code>Range=100,200</code>. The name <code>Range</code> has to be first part of the condition name. This settings will produce a filter on cluster energy which should be only from 100 to 200 keV (edges are included). | |||
* More ranges can be specified for one parameter. Individual conditions/ranges have to always include string <code>Range</code> in names, but they have to differ, therefore suffixes/endings have to be introduced: <code>Range_1</code>, <code>Range_2</code> etc. (, see example below). | |||
Filter with single condition in energy <code>E</code> from 100 to 200 keV: | |||
<pre> [E] | |||
Range = 100, 200 </pre> | |||
Filter with multiple conditions on energy <code>E</code>: | |||
<pre>[ | <pre> [E] | ||
Range_1=100,200 | |||
Range_2=500,1000 | |||
Range_asdasd=300,2000 </pre> | |||
<span id="names-of-cluster-parametersvariables"></span> | |||
=== Names of Cluster Parameters/Variables === | |||
The names can be found in the table about cluster variables in the section [[#SecClusterVarElist|Cluster Variables and Elist]]. | |||
=== Output into Extended Elist === | |||
< | The output of the filtering can be found in the extended elist (in directory <code>File</code>). New column is created with name <code>FilterPass</code>.<br /> | ||
The values of this new clomun/parameter are <code>1</code> for passing the filter and <code>0</code> for NOT passing the filter.<br /> | |||
All particles are stored in the extended elist, but it is possible to suppress the export of those particles which did not pass the filter conditions: | |||
< | <pre>DoExportElistExtFilter = false</pre> | ||
with this switch in the parameters file, only those particles which passed the filter are exported into extended elist (column <code>FilterPass</code> is still present). | |||
<span id="statistical-and-log-information"></span> | |||
=== Statistical and Log Information === | |||
< | Statistical information can be found in the general sampling list and it is also printed into log.<br /> | ||
Example for the sampling list in json format (explanations are given in the comments after <code>//</code>): | |||
<syntaxhighlight lang="cpp"> "CountParticle_Sum_FilterOk_cnt":4386, // Count of particles which passed the filter. | |||
"CountParticle_Sum_FilterOk_perc":58.93, // Percentage of the above to the total count of particles. | |||
"CountParticle_Sum_FilterOmit_cnt":3057, // Count of particles which did NOT pass the filter. | |||
"CountParticle_Sum_FilterOmit_perc":41.07, // Percentage of the above to the total count of particles. | |||
"CountParticle_Sum_All_cnt":7443, // Count of all particles which were evaluated.</syntaxhighlight> | |||
<pre> Filter conf file - name: Filter.ini | |||
Filter conf file - path: ./Test/data/test_007/ | |||
... | |||
-------------------------------------------------- | |||
FILTER | |||
Conditions are logically connected as AND for different variables (within always as AND) | |||
INDEX|NAME CONDITIONS [min|max] | |||
10|Roundness | |||
|-----------[0|2] | |||
'-----------[0.4|0.45] | |||
9|BorderPixCount | |||
|-----------[1|10] | |||
'-----------[20|1e+200] | |||
4|E | |||
'-----------[100|1e+300] | |||
-------------------------------------------------- | |||
... | |||
-------------------------------------------------- | |||
FILTER | |||
-------------------------------------------------- | |||
CountParticle_Sum_FilterOk [-]: 4386 | |||
CountParticle_Sum_FilterOk [%]: 58.93 | |||
CountParticle_Sum_FilterOmit [-]: 3057 | |||
CountParticle_Sum_FilterOmit [%]: 41.07 | |||
CountParticle_Sum_All [-]: 7443</pre> | |||
The first part informing about name and path to the config file of filter can be seen in the first stage of DPE processing.<br /> | |||
The second part is in the same processing/nationalization stage of DPE and it shows which filter were recognized and what ranges are about to be used.<br /> | |||
The last part informs about statistical information and it includes the same information which are also given in the json sampling list. | |||
== Histograms== | |||
< | One of the DPE outputs are histograms of cluster variables/parameters.<br /> | ||
The DPE in default exports 1D histograms of all cluster variables in extended elist and also several their combinations as 2D histograms.<br /> | |||
It is possible to export user defined 1D and 2D histograms which can be configured with a configuration file in the INI format (an example can be found in the program directory).<br /> | |||
The DPE allows to also create histograms of algebraic combinations of the cluster variables (multiplication, division, subtraction, addition).<br /> | |||
It is important that the general label (key name of a section in a INI file) of the histogram used in the INI file is unique to each histogram.<br /> | |||
If there are more histograms with one common name then the program only updates information about the first one in the INI file.<br /> | |||
The label itself in not used in the program itself and title and name of the histogram are set separately.<br /> | |||
The histograms are used in other analysis, therefore their changes/using user configuration might disturb e.g. creation significant vectors. This might produce following error for significant vectors: | |||
<pre>[ERROR] -1041 : Error occurred during significant vector module initialization. Module will not be used in processing.</pre> | |||
It just informs that the settings of histograms is not compatible with settings of significant vectors. | |||
<span id="inclusion-of-hist.ini-into-parameters-file"></span> | |||
=== Inclusion of Hist.ini into Parameters File === | |||
To use user configuration file for histograms, it is needed to add two parameters into paramters file, name (e.g. <code>Hist.ini</code>) and path (e.g. <code>PATH/TO/HIST/CONFIG/</code>) to the configuration file: | |||
<pre> | <pre>HistName - "Hist.ini" // Name of the configuration file for histograms. | ||
HistPath = "PATH/TO/HIST/CONFIG/" // Path of the configuration file for histograms.</pre> | |||
If at least the name is set and the file is found on the given location then the DPE uses this settings of histograms otherwise default configuration is used. | |||
< | <span id="histogram-1d"></span> | ||
=== Histogram 1D === | |||
The 1D histograms can be created as fixed bin width histograms or with variable binning (different width of bins).<br /> | |||
Lets assume that histograms of size should be created <code>Hist1D_Size</code>.<br /> | |||
The fixed bin width has following settings/needed parameters (example from configuration file with explanations = everything after #): | |||
<pre> | <pre> [Hist1D_Size] | ||
VarName="Size" # Name of column with given variable (see the first line in Elist.txt - it has to be the same) | |||
The same | Title="S" # Title of histogram (can be arbitrary - X if not given) | ||
NBin=9 # Number of bins | |||
Xmin=100 # Minimum value (if X = Xmin -> it is NOT included to the first bin - it has to be > Xmin) | |||
Xmax=1000 # Maximum value (if X = Xmax -> it is included to last bin NBin)</pre> | |||
The same as above but with different choice of variable, it is based on the column position in the elist instead of the cluster variable name: | |||
<pre> | <pre> [Hist1D_Size] | ||
ColIndex=7 # Position of the column - starts from 0 (it is Size in the example) | |||
Title="S" # -||- | |||
NBin=9 # -||- | |||
Xmin=100 # -||- | |||
Xmax=1000 # -||-</pre> | |||
These cases create histogram from cluster size with 9 bins from 100 to 1000 where one bin has width 100. Variable size bin width and its needed parameters (same histogram as the one above): | |||
<pre> [Hist1D_Size] | |||
VarName="Size" # -||- | |||
Title="S" # -||- | |||
BinLowEdge=100,200,300,400,500,600,700,800,900,1000 # Low edges of bins with Xmax (Xmin,...,Xmax -> size NBin+1) </pre> | |||
<span id="histogram-2d"></span> | |||
=== Histogram 2D === | |||
The 2D histograms are constructed in similar manner as 1D histograms. The fixed bin width histograms and their needed parameters: | |||
The | |||
<pre> | <pre> [Hist2D_SE] | ||
VarName="E","Size" # Name of columns with given variable - 1st is X, 2nd is Y (see the first line in Elist.txt - it has to be the same or see special variables) | |||
Title="E,S" # Title of histogram (can be arbitrary - X if not given) | |||
. | NBinX=9 # Number of X bins where X is in this case energy,E | ||
Xmin=100 # Minimum value of X (if X = Xmin -> it is NOT included to first bin - it has to be > Xmin) | |||
Xmax=1000 # Maximum value of X (if X = Xmax -> it is included to last bin NBin) | |||
NBinY=1000 # Number of Y bins | |||
Ymin=0 # Minimum value of X (if X = Xmin -> it is NOT included to first bin - it has to be > Xmin) | |||
Ymax=1000 # Maximum value of X (if X = Xmax -> it is included to last bin NBin)</pre> | |||
An example for cluster variables based on column positions in the elist: | |||
<pre> [Hist2D_SE] | |||
ColIndex=4,7 # Position of the columns which should be processed - 1st is X, 2nd is Y | |||
#...(the same as above) </pre> | |||
These cases create 2D histogram of cluster energy and size where energy is from 100 to 1000 with bin width of 100 and size is from 0 to 1000 with bin width of 1. Variable size bin width and its needed parameters: | |||
<pre> | <pre> [Hist2D_SE] | ||
VarName="E","Size" # -||- | |||
Title="E,S" # -||- | |||
Xmin | BinLowEdgeX=100,200,300,400,500,600,700,800,900,1000 # Low edges of X bins with Xmax (Xmin,...,Xmax -> size NBinX+1) | ||
BinLowEdgeY=100,200,300,400,500,600,700,800,900,1000 # Low edges of Y bins with Ymax (Ymin,...,Ymax -> size NBinY+1) </pre> | |||
It is also possible do just variable binning in one variable the second one can be with fixed bin size: | |||
... | |||
<pre> | <pre> [Hist2D_SE] | ||
VarName="E","Size" | |||
Title="E,S" | |||
NBinX=9 | |||
Xmin | Xmin=100 | ||
Xmax=1000 | |||
BinLowEdgeY=100,200,300,400,500,600,700,800,900,1000 </pre> | |||
<span id="additional-algebraic-operations"></span> | |||
=== Additional Algebraic Operations === | |||
Both 1D and 2D histograms allow additional operations of addition, subtraction, multiplication and division on extracted variables from EList. | |||
'''Histogram 1D''': | |||
<pre> ColIndex_Div=4,7 # Division as 5th/7th column (in this case E/S = Epix) - 1st argument is divided by 2nd argument | |||
ColIndex_Mult=4,7 # Multiplication as 5th*7th column - 1st argument is multiplied with 2nd argument | |||
ColIndex_Add=4,7 # Additions 5th+7th column- 1st argument is added to 2nd argument | |||
ColIndex_Subtr=4,7 # Subtraction as 5th-7th column - 2nd argument is subtracted from 1st argument</pre> | |||
An example based on cluster variable names: | |||
<pre> VarName_Div="E","Size" # Same thing as above but with names of columns | |||
VarName_Mult="E","Size" # Same thing as above but with names of columns | |||
VarName_Add="E","Size" # Same thing as above but with names of columns | |||
VarName_Subtr="E","Size"# Same thing as above but with names of columns </pre> | |||
'''Histogram 2D''': | |||
Very similar to 1D histograms but always the first given is X and the second given is Y. | |||
<pre> ColIndex=4 # X is set to 4th column - energy | |||
ColIndex_Div=4,7 # Y is set to division of 4th/7th</pre> | |||
Operations used for both variables: | |||
<pre> ColIndex_Div=8,7,4,7 # X is set to division of 8th/7th and Y is set to division of 4th/7th</pre> | |||
The same thing can be done with names of columns/variables VarName: | |||
<pre> VarName_Div="E","Size","Height","Size" # X is E/Size and Y is Height/Size</pre> | |||
<span id="text-export"></span> | |||
=== Text export === | |||
There are two kinds of exported files: histogram data and histogram information/info file.<br /> | |||
The histogram data file includes content of bins and bins low edges. The histogram info file comprehends features of the histogram: title, name, count of bins etc. (see more details below).<br /> | |||
The data file has a suffix <code>.hist</code> and the info file <code>.hist_info</code>. The data file has following formatting for 1D histogram: | |||
< | <pre> Xmin BinCont_1 | ||
BinLowEdge_2 BinCont_2 | |||
... .... | |||
BinLowEdge_NBin BinCont_NBin | |||
Xmax Overflow</pre> | |||
The Xmax is included to allow an user easier read of this file without direct need to also read the info file (all needed information is in the data file in the base case).<br /> | |||
It is also possible to export only those bins which have non zero content with fixed binning. This can be done with following option in the Hist.ini file: | |||
<pre>DoSparseExport=1 # Check whether only nonzero bins should be exported (1 for true and 0 for false).</pre> | |||
It has to be used for each histogram separately and it is used as default settings. The exported data file has following format: | |||
<pre> | <pre> Xmin BinCont_1 | ||
BinLowEdge_2 BinCont_2 | |||
... ... | |||
BinLowEdge_i BinCont_i != 0 | |||
... ... | |||
Xmax Overflow </pre> | |||
The <code>Xmin,BinLowEdge_2,Xmax</code> are always exported even if their bin content is 0 for further reading and reconstruction of histograms. This option is not functional for variable binning to avoid a lack of information in the data file for histogram complete reconstruction in post-processing. To reconstruct the histogram, it can be done just based on the data file even in the case of sparse export because the missing bins<br /> | |||
and bin width can be calculated based on the first two bins in the data file (BinWidth = BinLowEdge_2 - Xmin).<br /> | |||
Anyway, it is recommended to read the info file for example in the case of constant bin width combined with the sparse export to ensure that given<br /> | |||
histogram is truly with constant binning and not variable binning (written in the parameter: <code>BinEquiDist=1</code> = it is const binning and 0 for variable binning). | |||
Similar approach is utilized for the 2D histograms: | |||
<pre> | <pre> Xmin Ymin BinCont_1 # First bin | ||
XBinLowEdge_2 Ymin BinCont_2 | |||
... .... .... | |||
XBinLowEdge_NBinX Ymin BinCont_NBinX | |||
Xmin YBinLowEdge_2 BinCont_NBinX+1 | |||
... .... .... | |||
XBinLowEdge_NBinX YBinLowEdge_2 BinCont_NBinX+NBinY | |||
XBinLowEdge_2 YBinLowEdge_3 BinCont_NBinX+NBinY+1 | |||
... .... .... | |||
XBinLowEdge_NBinX YBinLowEdge_NBinY BinCont_NBinX*NBinY # Last bin | |||
Xmax Ymax Overflow</pre> | |||
The sparse export has following form. | |||
<pre> | <pre> Xmin Ymin BinCont_1 | ||
XBinLowEdge_2 Ymin BinCont_2 | |||
... .... .... | |||
XBinLowEdge_i Ymin BinCont_i != 0 | |||
... .... .... | |||
Xmin YBinLowEdge_2 BinCont_NBinX+1 | |||
... .... .... | |||
XBinLowEdge_m YBinLowEdge_n BinCont_m+n*NBinX != 0 | |||
... .... .... | |||
Xmax Ymax Overflow</pre> | |||
It is similar to 1D histogram with the exception that there is also the next Y low bin edge to also estimate the width of binning for Y axis. | |||
The info file is formatted as an INI file and it includes the same information as the Hist.ini file for each histogram individually. There two additional features compared with the Hist.ini: | |||
* Statistical information (mean, err of mean, std, err of std) | |||
* Overflow and underflow information | |||
Explanation of individual parameters are written as comment in example below: | |||
<syntaxhighlight lang="ini"> [HistPar] | |||
- | Title="Hist1D_Epix_1" # Title of the histograms used for unique recognition | ||
Name="Mean energy per pixel of cluster" # General name | |||
AxisTitles="Epix [keV/px]","N [-]" # Names of axis titles/lables. | |||
NBin=300 # Count of bins | |||
Xmin=0 # HIstogram minimum | |||
Xmax=300 # Histogram maximum | |||
BinEquidist=1 # Check whether histogram has same width bins | |||
[HistCont] | |||
Underflow=0 # Count all events which are below Xmin | |||
Overflow=1 # Count all events which are above Xmax | |||
[Statistics] | |||
Mean=29.6869 # Estimation of weighted mean value of histogram (weighted with bin contents) | |||
Mean_Err=0.20951 # Error of the mean estimation | |||
Std=13.9273 # Estimation of weighted standard deviation of histogram (weighted with bin contents) | |||
Std_Err=1.40088 # Error of the std estimation</syntaxhighlight> | |||
Equivalent information is also given for histogram 2D with appropriate extension of additional dimension. | |||
<span id="graphical-export"></span> | |||
=== Graphical export === | |||
[[File:Hist1d energy.png|thumb|414x414px|Derived deposited energy histograms with equidistant binning. Data displayed in lin scale (top plot) and log scale (bottom plot).]] | |||
There is possibility to create plot of Hist1D and Hist2D via python matplotlib, it has to be preinstalled.<br /> | |||
The current version will automatically create these plots if <code>DoExportGraphics=true</code> but it can be turned off - see below. | |||
<pre> SOME_PREVIOUS_CODE # This is some code to call histogram 2D | |||
Name="Histogram2D title in plot" # Name o histogram in plot | |||
AxisTitle="X","Y","N" # Name of axis in the plot</pre> | |||
The plots can be with logarithmic scales on all axis: | |||
<pre> SOME_PREVIOUS_CODE # This is some code to call histogram 2D | |||
DoLogX=1 # 1 for true=do logarithmic X axis and 0 for false | |||
DoLogY=1 # 1 for true=do logarithmic Y axis and 0 for false</pre> | |||
To turn off plotting: | |||
[[File:2dhist es.png|thumb|418x418px|2D histogram of cluster size and energy. The figure is accompanied with statistical information about mean values, standard deviations, sum/integral, maximum and minimum value.]] | |||
<pre> SOME_PREVIOUS_CODE # This is some code to call histogram 2D | |||
DoPlot=0 # Default is 1 - do plot and 0 means not do plot</pre> | |||
It is possible to adjust the number binning of the histogram only for the graphics. This can be done with <code>Hist1DGraphicsRebin</code> in the '''main configuration file''': | |||
= | <pre> Hist1DGraphicsRebin = 100</pre> | ||
This option with value 100 will change the shown number of bins in the graphics to 100 but the exported files includes original binning.<br /> | |||
This is allowed for faster exports of histograms with more than 10000 bins which can time challenging for a PC. | |||
<span id="log-information"></span> | |||
=== Log Information === | |||
<pre> | <pre> Hist conf file - name: Hist.ini | ||
Hist conf file - path: ./Test/data/test_006/ | |||
... | |||
-------------------------------------------------- | |||
HISTOGRAMS | |||
1D: | |||
Hist1D_E_Total from column 4 | |||
Hist1D_S_Total from column 7 | |||
Hist1D_LET_Total from column 26 | |||
Hist1D_EpixMean_Total from column 21 | |||
Hist1D_Epix_Total from column-division: 4/7 | |||
Hist1D_H_Total from column 8 | |||
Hist1D_LET_VarBinWidth_Total from column 26 | |||
2D: | |||
Hist2D_ES_Total from column 4 7 | |||
Hist2D_ES_VarBinWidth_Total from column 4 7 | |||
-------------------------------------------------- | |||
... | |||
Hist are be shown with max N bins: 100 bins | |||
... | |||
Exporting histograms to: ./Test/out/test_006/./Hist/ | |||
Exporting histogram 1D: Hist1D_E_1 | |||
Exporting histogram 1D: Hist1D_S_1 | |||
Exporting histogram 1D: Hist1D_LET_1 | |||
Exporting histogram 1D: Hist1D_EpixMean_1 | |||
Exporting histogram 1D: Hist1D_Epix_1 | |||
Exporting histogram 1D: Hist1D_H_1 | |||
Exporting histogram 1D: Hist1D_LET_VarBinWidth_1 | |||
...</pre> | |||
The first part informing about name and path to the config file of histograms can be seen in the first stage of DPE processing.<br /> | |||
The second part is in the same processing/nationalization stage of DPE and it shows which histograms will be produced.<br /> | |||
The last tow parts informs exporting of histograms (first one only if change of hist binning is used in graphical plots). | |||
== Significant Vectors== | |||
The significant vectors is a unique set/vector of numbers which describes given data sample/radiation field.<br /> | |||
It can be used for radiation field recognition as it is done in the DPE.<br /> | |||
There is a possibility to create own configuration file for the generation of the significant vector.<br /> | |||
The significant vectors are exported into directory <code>SigVec</code> in path <code>FileOutPath</code>.<br /> | |||
The analysis/generation of SigVec can be suppressed with the following option in Parameters file: | |||
<pre> | <pre>DoSigVec = false</pre> | ||
The results are saved into directory <code>SigVec</code> into files <code>SigVec.vec</code> and <code>SigVec.vec_info</code>. First file includes the vector itself and second file includes information about the vector (used histograms, intervals etc.). | |||
==== Theoretical Introduction ==== | |||
=== | The significant vectors are sets of numbers which should describe '''uniquely''' a radiation field. They are created form 1D and 2D histograms based on the following rules (see the illustrations in the ): | ||
# '''1D histograms''': Each number of the significant vector is produced as a sum of bin contents within some specified ranges. Assuming that some elements of the vector should be done based on the energy distribution/histogram. Ranges of interests have to be specified where an integration is executed. This specification of the ranges is a task for an optimization method and differentiates based on the set of radiotin fields which should be recognized. One number is created for each range. These are normalized with a total integral over whole histogram to achieve comparability with different fields where more or less particles could be measured. A weighting can be also introduced as a last stage of the process to highlight some elements of the vector. This procedure can be done for several other cluster variables resulting in the vector of significant variables. | |||
# '''2D histograms''': Their conversion into significant variables obey similar rules as the 1D histograms. Except for the ranges, masks can be also specified which determines the area of interest. The integration is based on this mask where the mask has the same features as the histogram itself with values which can be used as weights. The mask values are used as a multiplicative factor in the integration therefore bins with 0 are ignored and bins with value higher than 1 are highlighted. | |||
[[File:DPE RFR SigVec 1 b.png|748x748px]][[File:DPE RFR SigVec 2 b.png|814x814px]] | |||
<span id=" | <span id="inclusion-of-sigvec.ini-into-parameters-file"></span> | ||
=== | === Inclusion of SigVec.ini into Parameters File === | ||
< | The the DPE will include a user configuration file only if it is written in the PararamtersFile.<br /> | ||
Two parameters are used for this purpose (assume that name of config file is <code>SigVec.ini</code> in path <code>/PATH/TO/SIGVEC/CONFIG/</code>): | |||
<pre> SigVecName = "SigVec.ini" // Name of the configuration file for significant vectors. | |||
SigVecPath = "/PATH/TO/SIGVEC/CONFIG/" // Path to the configuration file for significant vectors.</pre> | |||
If at least the name is set and the file is found on the given location then the DPE uses this settings of SigVec otherwise default configuration is used. | |||
<span id="histogram-1d---intervals-of-interest"></span> | |||
=== Histogram 1D - Intervals of Interest === | |||
A histogram is examined and all bins (bin contents) with bin center between up and down edge/limit of given interval is summed, number R_1.<br /> | |||
This number is used as one element of SigVec = {R_1, R_2, … , R_N} for N of intervals.<br /> | |||
At least two numbers has to be given - down and up edge of the interval where interval condition for bins are following: DOWN_EDGE < BinCenter <= UP_EDGE.<br /> | |||
An example of configuration file for 1D histogram: | |||
<syntaxhighlight lang="ini"> [Var_Epix] # Section name as unique key of the given SigVec element. It has to include the string "Var_". | |||
Title=" | Title="Epix" # Title of the variable. | ||
Intervals= | HistName="Hist_Epix" # Name of histogram which should be processed - same as name of histogram file without suffix/ending in directory Hist (data suffix) | ||
Intervals=0,200,100,500,500,5000 # Limits of intervals - 1st interval = from 0 to 200, 2nd interval = from 100 to 500 etc. (N in this case is 3)</syntaxhighlight> | |||
The intervals can be overlapping. '''The section name has to include string <code>Var</code> in its name and they have to differ for different elements of the SigVec: <code>Var_E, Var_S, ...</code>'''.<span id="histogram-2d---regions-of-interest"></span> | |||
=== Histogram 2D - Regions of Interest === | |||
< | A histogram is multiplied with mask (from 0 val to inf - can serve as weights) with the same number of bins.<br /> | ||
Then the histogram is summed and this number is as an input to the SigVec.<br /> | |||
For each given mask one number is produced. An example of configuration file: | |||
< | <syntaxhighlight lang="ini"> [Var_E_S] # Section name as unique key of the given SigVec element. It has to include the string "Var_". | ||
Title="E_S" # Title of the variables | |||
HistName="Hist_E_S_2" # Name of histogram which should be processed - same as name if histogram file without suffix/ending in directory Hist (data suffix) | |||
Mask_1="\PATH\TO\MASK\MASK1_FILE_NAME" # Path to masks which should be applied - same dimensions as histogram (see examples) | |||
Mask_2="\PATH\TO\MASK\MASK2_FILE_NAME" # Path to masks which should be applied - same dimensions as histogram (see examples)</syntaxhighlight> | |||
Lets assume that the <code>Hist_E_S_2</code> is following histogram: NBinX = 5, Xmin = 0, Xmax = 10 & NBinY = 4, Xmin = 0, Xmax = 4 then the mask should have following form: | |||
<pre>0 1 0 0 0 | |||
0 0 0 0 0 | |||
0 2 0 0 3 | |||
0 0 0 4 10</pre> | |||
A single gap is used as number separator.This will use bins [X,Y,Weight] = [2,1,1], [2,3,2], [5,2,3], [4,4,4], [4,5,10] and the matrix is read from first line to the last one (1st line is Y = 1).<span id="normalization"></span> | |||
=== Normalization === | |||
Both 1D and 2D histograms can be normalized before they are used calculation of the SigVec: | |||
<pre> Normalize=1 # 1 to DO normalization - 0 to NOT do norm. - default is 0</pre> | |||
<span id="default-configuration"></span> | |||
=== Default Configuration === | |||
<syntaxhighlight lang="ini"> [Var_S] | |||
Title="Var_S" | |||
Intervals=-1,5,5,20,20,1e+200 | |||
Normalize=1 | |||
[Var_BorderPixN] | |||
Title="Var_BorderPixN" | |||
Intervals=-1,2,2,5,5,15,15,1e+200 | |||
Normalize=1 | |||
[Var_Lin] | |||
-- | Title="Var_Lin" | ||
Intervals=-1,0.5,0.5,0.95,0.95,1e+200 | |||
Normalize=1</syntaxhighlight> | |||
Minimum values -1 and maximum 1e+200 are safety precautions to be sure that all relevant bins are taken into account. | |||
<span id="log-information-1"></span> | |||
< | === Log Information === | ||
<pre>.. | <pre> SigVec conf file - name: SigVec.ini | ||
SigVec conf file - path: ./Test/data/test_008/ | |||
... | |||
-------------------------------------------------- | |||
SIGNIFICANT VECTOR | |||
Features: | |||
Title = Epix | |||
HistName = Hist1D_EpixMean_Total | |||
Intervals | |||
Normalize= True | |||
Title = H | |||
HistName = Hist1D_H_Total | |||
Intervals | |||
Normalize= True | |||
Title = LET | |||
HistName = Hist1D_LET_Total | |||
Intervals | |||
Normalize= True | |||
Title = E | |||
HistName = Hist1D_E_Total | |||
Intervals | |||
Normalize= True | |||
Title = S | |||
HistName = Hist1D_S_Total | |||
Intervals | |||
Normalize= True | |||
-------------------------------------------------- | |||
... | |||
</pre> | |||
The first part informing about name and path to the config file of significant vector can be seen in the first stage of DPE processing (if it is not shown then default option is used).<br /> | |||
The second part is in the same processing/nationalization stage of DPE and it shows which histograms and intervals/masks will be used to creation of the vector. | |||
<pre>... | |||
======================================================================= | |||
SIGNIFICANT VECTOR PROCESSING | |||
======================================================================= | |||
* Processing histogram: Hist1D_EpixMean_Total | |||
Histogram intervaling: [0,50][50,1e+200] | |||
Intervaling results: 0.967578 0.032422 | |||
* Processing histogram: Hist1D_H_Total | |||
Histogram intervaling: [0,50][50,200][200,1e+200] | |||
Intervaling results: 0.596236 0.402650 0.001114 | |||
* | * Processing histogram: Hist1D_LET_Total | ||
Histogram intervaling: [0,1][1,3][3,1e+200] | |||
Intervaling results: 0.915691 0.084309 0.000000 | |||
* Processing histogram: Hist1D_E_Total | |||
Histogram intervaling: [0,200][200,1000][1000,1e+200] | |||
Intervaling results: 0.686559 0.313417 0.000024 | |||
* Processing histogram: Hist1D_S_Total | |||
Histogram intervaling: [0,5][5,20][20,1e+200] | |||
Intervaling results: 0.595407 0.387884 0.016709 | |||
< | </pre> | ||
This last parts includes information about created elements of the vector, from which histograms were created and what are the results. | |||
== Particle Identification== | |||
The DPE engine is also capable of basic particle identification/classification.<br /> | |||
This output can be exported into extended elist where new column/PIDClass will be created with information about recognized class. | |||
The choice of the PID algorithm is based on the switch <code>PIDAlg</code> whose numerical value are listed below. The settings can be done in the main configuration file: | |||
<pre> PIDAlg = 201</pre> | |||
the current possible values are following: | |||
<ol start=" | * 101 - Heuristic simplified DT for TPX 300 um Si | ||
<li>Dense neural network | * 102 - Heuristic decision tree with 8 classes for TPX 300 um Si | ||
* 103 - Heuristic decision tree with 16 classes for TPX3 500 um Si | |||
* 201 - Dense neural network with 3 classes for TPX 300 um Si | |||
* 202 - Dense neural network with 6 classes for TPX 300 um Si | |||
* 251 - Dense neural network with 3 classes for TPX3 500 um Si | |||
* 252 - Dense neural network with 6 classes for TPX3 500 um Si | |||
If not value is given to the program then DT is chosen based on the detector configuration. | |||
<span id="algorithms-based-on-neural-networks"></span> | |||
=== Algorithms based on Neural Networks === | |||
'''Dense neural networks for TPX 300 um Si''' | |||
<ol start="201" style="list-style-type: decimal;"> | |||
<li>Dense neural network TPX with 3 classes</li></ol> | |||
# Protons | # Protons | ||
Line 1,314: | Line 1,982: | ||
# Others | # Others | ||
<ol start=" | <ol start="202" style="list-style-type: decimal;"> | ||
<li>Dense neural network | <li>Dense neural network TPX with 6 classes<br /> | ||
</li></ol> | </li></ol> | ||
Line 1,326: | Line 1,994: | ||
# Others | # Others | ||
'''Dense neural networks for TPX3 500 um Si''' | |||
<ol start="251" style="list-style-type: decimal;"> | |||
<li>Dense neural network TPX3 with 3 classes</li></ol> | |||
<ol start="101" style="list-style-type: decimal;"> | # Protons | ||
<li>Heuristic simplified DT with 3 classes:<br /> | # Photons && Electrons | ||
</li></ol> | # Ions | ||
# Others | |||
<ol start="252" style="list-style-type: decimal;"> | |||
<li>Dense neural network TPX3 with 6 classes<br /> | |||
</li></ol> | |||
# Protons LE (<30 MeV) | |||
# Protons ME (>30 & <100 MeV) | |||
# Protons HE (>100 MeV) | |||
# Photons && Electrons | |||
# Helium ions | |||
# Ions (except He) | |||
# Others | |||
<span id="algorithms-based-on-decision-trees"></span> | |||
=== Algorithms based on Decision Trees === | |||
'''Decision tree for TPX 300 um Si''' | |||
<ol start="101" style="list-style-type: decimal;"> | |||
<li>Heuristic simplified DT with 3 classes:<br /> | |||
</li></ol> | |||
# Low LET: electrons, X-rays, gamma | # Low LET: electrons, X-rays, gamma | ||
Line 1,467: | Line 2,157: | ||
To use the RFR, it is needed to use standard/default settings of the DPE with respect to the histograms and significant vectors. | To use the RFR, it is needed to use standard/default settings of the DPE with respect to the histograms and significant vectors. | ||
[[File:Rfr sigvec distance.png|1084x1084px]] | |||
=== Distance Comparator === | === Distance Comparator === | ||
Line 1,640: | Line 2,330: | ||
"Dose_Frame_Gy+1":[ 0.000000000047616, 0.000000000052289, 0, 0.000000000011844, 0.000000000012407, 0, 0 ], // Dose in individual frames | "Dose_Frame_Gy+1":[ 0.000000000047616, 0.000000000052289, 0, 0.000000000011844, 0.000000000012407, 0, 0 ], // Dose in individual frames | ||
"DoseRate_Frame_Gy+1s-1":[ 0.0000000095232, 0.000000010458, 0, 0.0000000023689, 0.0000000024814, 0, 0 ] // Dose rate in individual frames </syntaxhighlight> | "DoseRate_Frame_Gy+1s-1":[ 0.0000000095232, 0.000000010458, 0, 0.0000000023689, 0.0000000024814, 0, 0 ] // Dose rate in individual frames </syntaxhighlight> | ||
=== Graphical Output === | |||
<gallery mode="packed" heights="300"> | |||
File:Frame analysis occupancy.png|Graph of frame occupancy in dependence on frame start time. | |||
File:Frame analysis visual.png|Visualisation of tracks for one frame (frame with highest occupancy). The plot also includes additional information about the features of the seen tracks (energy sum, occupancy etc.). | |||
</gallery> | |||
== Compton Directional Analysis== | == Compton Directional Analysis== | ||
<span id="notes-1"></span> | |||
[[File:Cocam intro.png|818x818px]] | |||
<span id="notes-1"></span> | <span id="notes-1"></span> | ||
Line 1,710: | Line 2,412: | ||
"CountCompGroup_Sum_Sample_cnt":[ 4, 8 ] // Count of all Compton groups in individual time samples</syntaxhighlight> | "CountCompGroup_Sum_Sample_cnt":[ 4, 8 ] // Count of all Compton groups in individual time samples</syntaxhighlight> | ||
== Directional Analysis== | === Graphical Output === | ||
[[File:Cocam projection.png|thumb|469x469px|Visualisation of the reconstruction and its projection. The source was placed in the position 10 of X and 20 on Y.]] | |||
To simplify the created projection alias estimation of the radiation source position, an additional image post-processing was introduced: | |||
* '''Contrast''' – removes all pixels whose values is less than fraction of the maximal value in pixels. Default is set to 70 percent alias 0.7 in the fraction terms. | |||
* '''Blurring''' – exploits Gaussian filter to blur the image and highlight the estimation. | |||
[[File:Cocam_proj_image_adjust.png|865x865px]] | |||
== Directional Analysis== | |||
DPE is capable to estimate several directional features of detected particles. | DPE is capable to estimate several directional features of detected particles. | ||
Line 1,739: | Line 2,469: | ||
"ElevationAngle_Mean_Sample_deg+1":[ 72.8143, 0, 72.7298, 0, 0, 72.6255, 0, 0, 72.5953, 0, 72.6597, 0, 0, 72.9396, 0, 0, 72.8254, 0, 72.4583, 0, 0, 72.7155, 0, 0, 72.7478, 0, 72.8675, 0, 0, 72.7720 ] </syntaxhighlight> | "ElevationAngle_Mean_Sample_deg+1":[ 72.8143, 0, 72.7298, 0, 0, 72.6255, 0, 0, 72.5953, 0, 72.6597, 0, 0, 72.9396, 0, 0, 72.8254, 0, 72.4583, 0, 0, 72.7155, 0, 0, 72.7478, 0, 72.8675, 0, 0, 72.7720 ] </syntaxhighlight> | ||
== Coincidence and Event Analysis== | === Graphical Output === | ||
* '''TimeCoinc''' - [FLOAT] Time in nanoseconds for evaluations coincidence events -> clusters whose times are in between this this interval are grouped as coincidence group (first cluster time is down edge and plus coince group time is upper edge). Example: <code>TimeCoinc = 1e2</code>, <code>TimeCoinc = 10</code>, <code>TimeCoinc = 1.2554</code>. Default: 100. | |||
<gallery mode="packed" heights="250"> | |||
<span id="information-3"></span> | File:Direction 2dhist.png|''2D histogram of cluster elevation angle and polar angle. Given information form this plot reveals a directional orientation of the source which was in this case beam of high energetic protons impacting under elevation angle of approximately 60 degrees.'' | ||
=== Information === | File:Dir analysis 2dhist sphere.png|''2D histogram shown on a sphere to more emphasize the original direction of detected particles. Data used in this plot are of same origin as in the previous image.'' | ||
</gallery> | |||
<syntaxhighlight lang="cpp"> "TimeCoincWindow_ns+1":32.8120, // Used coincidence window | |||
"CountEvent_Sum_cnt":9, // Total count of all events (not particles -> event is one particle or group of coincidence particles) | |||
"CountCoincGroup_Sum_cnt":2, // | |||
"CountCoincGroup_Sum_perc":22.2222, // | |||
"CountParticle_Mean_Event_cnt":1.4444, // | |||
"CountParticle_Sum_NoCoinc_cnt":7, // | |||
"CountParticle_Sum_NoCoinc_perc":53.8462, // | |||
"CountParticle_Sum_CoincGroup_cnt":6, // | |||
"CountParticle_Sum_CoincGroup_perc":46.1538, // | |||
"CountParticle_Mean_CoincGroup_cnt":3, // | |||
"CountParticle_Max_CoincGroup_cnt":4, // | == Coincidence and Event Analysis== | ||
"TimeLive_Int_Sample_s+1":[ 1, 2, 3, 3.0000 ], // | * '''TimeCoinc''' - [FLOAT] Time in nanoseconds for evaluations coincidence events -> clusters whose times are in between this this interval are grouped as coincidence group (first cluster time is down edge and plus coince group time is upper edge). Example: <code>TimeCoinc = 1e2</code>, <code>TimeCoinc = 10</code>, <code>TimeCoinc = 1.2554</code>. Default: 100. | ||
"CountParticle_Mean_Sample_Event_cnt":[ 1, 0, 1, 2 ],// | |||
"CountCoincGroup_Sample_cnt":[ 0, 0, 0, 2 ] //</syntaxhighlight> | <span id="information-3"></span> | ||
=== Information === | |||
<syntaxhighlight lang="cpp"> "TimeCoincWindow_ns+1":32.8120, // Used coincidence window | |||
"CountEvent_Sum_cnt":9, // Total count of all events (not particles -> event is one particle or group of coincidence particles) | |||
"CountCoincGroup_Sum_cnt":2, // | |||
"CountCoincGroup_Sum_perc":22.2222, // | |||
"CountParticle_Mean_Event_cnt":1.4444, // | |||
"CountParticle_Sum_NoCoinc_cnt":7, // | |||
"CountParticle_Sum_NoCoinc_perc":53.8462, // | |||
"CountParticle_Sum_CoincGroup_cnt":6, // | |||
"CountParticle_Sum_CoincGroup_perc":46.1538, // | |||
"CountParticle_Mean_CoincGroup_cnt":3, // | |||
"CountParticle_Max_CoincGroup_cnt":4, // | |||
"TimeLive_Int_Sample_s+1":[ 1, 2, 3, 3.0000 ], // | |||
"CountParticle_Mean_Sample_Event_cnt":[ 1, 0, 1, 2 ],// | |||
"CountCoincGroup_Sample_cnt":[ 0, 0, 0, 2 ] //</syntaxhighlight> | |||
=== Graphical Output === | |||
<gallery mode="packed" heights="250"> | |||
File:Coinc visual e.png | |||
File:Coinc visual t.png | |||
File:Coinc 1dhist t diff.png | |||
File:Coinc graph count events.png | |||
</gallery> | |||
== Spatial Maps== | |||
…. | |||
== Event Visualization== | |||
Another graphical output produced in the engine are visualisations of particle tracks in sensor plane. These plots can be only obtained for raw data containing pixelated information. Example of this output can be seen in the where deposited energy in each pixel is used for this purpose. | |||
The individual events visualisations are completed with so called integrated visualisations. In these plots several particles are shown together to obtain more statistical based information about detected tracks. Two versions of these plots are available: integration of N particles (N is usually around 100) and integration of all detected particles. | |||
<gallery mode="packed" heights="250"> | |||
File:Event visual cluster.png|''Visualisation of energy of particle track in sensor plane for alpha particle created in decay of Am241. It is accompanied with information about track/cluster parameters/variables.'' | |||
File:Event visual sum.png|''Integrated event visualisation of all detected particles in one plot. Figure is accompanied with statistical information.'' | |||
File:Event viasual partial.png|''Integrated event visualisation of 100 detected particles in one plot. Figure is accompanied with statistical information.'' | |||
</gallery> | |||
Event visualisation in graphical form is completed with text alternative in matrix form. | |||
[[File:Event visual text.png|668x668px]] | |||
== Error Codes== | == Error Codes== | ||
The engine run can produce an error code into standard error stream (<code>stderr</code>) and into the standard output stream (<code>stdout</code>). Successful engine run produces/returns <code>0.</code> Any other value signals error in the run. Possible values are listed below (negative numbers with explanation after <code>///<</code>): | The engine run can produce an error code into standard error stream (<code>stderr</code>) and into the standard output stream (<code>stdout</code>). Successful engine run produces/returns <code>0.</code> Any other value signals error in the run. Possible values are listed below (negative numbers with explanation after <code>///<</code>): | ||
Line 1,865: | Line 2,630: | ||
</syntaxhighlight> | </syntaxhighlight> | ||
== Application Programing Interface - | == Application Programing Interface - C++== | ||
… | … |
Latest revision as of 09:38, 30 August 2024
The DPE (Data Processing Engine) is a software tool which serves for a processing of data acquired with the Timepix detectors (TPX, TPX2, TPX3). The following list includes the main processing parts of the DPE:
- Pre-processing: clustering of data with application of calibration and corrections, calculation of cluster variables/parameters, creation of histograms
- Processing: particle recognition and radiation field recognition
- Post-processing: directional analysis, coincidence analysis, Compton camera, generation of physical products and their time evolution with respect to the radiation field classification (fluxes, dose rates etc.)
The following information describes its command line platform and it is valid for version 1.1.0.
Download
In the case of interest, please contact: support@advacam.cz or lukas.marek@advacam.cz
Prerequisites, Install and Run
Prerequisites
To create graphics, python3 and several python packages are needed:
- matplotlib (version >= 3.4)
- numpy
- multiprocessing - only if creation of graphics should be done on several threads/cores of CPU
The program, especially in the case of graphic export, needs several hundreds of MB on RAM (app 500 MB but it is based on the size of exported histograms). If the multiprocessing is used then the RAM usage can be up several GB and all threads of the CPU might be used.
Install
In the directory where the program should be installed, extract all the files from downloaded archive/zip.
Run the program on linux
The DPE is started with the following program in the command line from the installation directory:
./dpe.sh PARAMETERS_FILE_PATH/PARATERS_FILE_NAME
PARAMETERS_FILE_PATH
is path to the parameters file and PARATERS_FILE_NAME is its name. It is possible that dpe.sh and clusterer have to be allowed as an executable: chmod +x clusterer
.
Run the program on windows
DPE.exe PARAMETERS_FILE_PATH\PARATERS_FILE_NAME
If your path or name of the parameters includes white spaces enclose it in quotes: “PARAMETERS_F ILE_ PATH_F ILE_NAME”. The input parameter into the DPE containing the path and name of the parameter file can be omitted if the file is in the same directory as the DPE and its name is ParametersFile.txt.
First Example of Run
This sections includes information about the first run of the DPE based on the ParamteresFile given with the program example. Example directory can be found and downloaded for the same download link. It includes needed inputs for the example to run on linux and windows systems. In the case of running example on linux system then it is needed to download following directories: Linux, Clusterer, Example. The example can be run with following command from the Linux directory:
chmod +x dpe.sh
./dpe.sh ./ParametersFile.txt
The output should be exported into the Example/Output directory. The given data is a mixed sample of proton, electron and alpha ion data measured with TPX3 500 um. Run of the DPE should read the main parameters file ParametersFile.txt and create several export directories with files and graphics:
- EventVisual - individual cluster plots, integrated sensor plots and plots of clusters for individual PID classes.
- File - main data files: elist, cluster log and sampling list.
- Graph - graphical representation of the sampling list.
- Hist - histograms of cluster variables.
- SpatialMap - spatial maps of cluster variables also with respect to the PID classes.
- SigVec - significant vector which uniquely specifies the given radiation field.
- Direction - results of directional analysis: estimation of radiation direction.
- CoincEvent - results of coincidence analysis: evaluates time correlation between clusters.
The main configuration is done in the ParametersFile.txt and meaning of individual parameters can be found directly there or in this document. Other files in the Example directory determine the creation of histograms (Hist.ini), masking of the raw data (Mask.txt), filtering of cluster (Filter.ini) and creation of the SigVec (SigVec.ini). Pre-generated referential results can be found in the directory Ref. More information about their format can be found in the next sections in this document.
The similar example can be also made on Windows with alternative command from the directory Windows:
DPE.exe .\ParametersFile.txt
The output and its description is the same as the for the linux system above.
Additional Examples
Additional examples can be found in the directory Example/Test which includes several examples from testing. The input data are in directory data and its output can be found in the directory ref (referential output). The examples comprehend following cases:
- T3PA input - t3pa input with default settings (Am241 source).
- Several T3PAs input - 20x t3pa files as input with default settings (Am241 source).
- CLOG tot+toa TPX3 input - clog tot+toa from TPX3 with default settings (Am241 source).
- CLOG tot TPX input - clog tot from TPX with default settings (protons HE).
- CLOG itot+count TPX3 input - clog itot+count from TPX3 with default settings (Am241 source).
- Hist with T3PA - Hist settings with T3PA file as input.
- Filter with T3PA - Filter settings with T3PA file as input (Am241 source).
- SigVec with T3PA - SigVec settings with T3PA file as input
- Mask with T3PA - Mask settings (general file) with T3PA file as input
- Comparator with T3PA and CdTe - Comparator with T3PA file as input
- Elist old input - elist as input with default settings
- Generated t3pa input - pregenerated t3pa for physical variables etc.
- SigVec and Hist with T3PA - SigVec and Hist settings with T3PA file as input to produce SigVec from 2D hist.
- CLOG event input - clog event with default settings.
- CLOG toa input - clog toa with default settings.
- CLOG tot+time proc input - clog tot+time after processing with default settings.
- T3P input - t3p input with default setting (ba133)
- T3R input - t3r input with default settings (upper part of detector was masked)
- Matrix event input - matrix event input with default settings
- Matrix tot input - matrix tot input with default settings
- Elist extended input - elist extended input with default settings
- Graphics output CLOG tot+toa TPX3 input - graphics output for CLOG tot+toa TPX3 input with default settings
- Several T3PAs input, just dir name - several t3pa files, just specifying dir name as input with default settings
- Comparator with own DB - input is own DB for camparator and settings of SigVec and Hist for t3pa input
- Print program version - print program version with -v
- Print program help - print program help with -h
- Ignore clusters at sensor edge - remove clusters at sensor edge with default settings with T3PA file as input
- Pixet mask with T3PA - Mask settings (pixet mask) with T3PA file as input
- Acceptance of all param - param file with all params and their acceptanc.
- CLOG tot+toa meas input - clog tot+toa from measurement with default settings
- Generated t3pa 2x - pregenerated t3pa two times to check connection of two separate measurement files
- Coincidence analysis - coincidence and matching analysis (time coincidence window and products)
- Directional analysis - directional analysis with proton 200 MeV 75 deg data.
- Compton analysis single forward - Compton analysis for single layer forward (other analysis are off due to pregenerated data).
- Compton analysis single forward - Compton analysis for single layer forward (other analysis are off due to pregenerated data).
- RadFiledRecog analysis - radiation field recognition for 2mm CdTe TPX3 with Ba133 data.
- RadFiledRecog analysis - radiation field recognition for 2mm CdTe TPX3 with Cs137 data.
- Frame analysis - frame analysis.
- Detector geometry - setting of geometry via geo det config file.
- Adv elist - advanced elist.
- Adv elist extended - advanced elist extended in dpe.
- All analysis and their switches - using all analysis and turning them off compatibility.
- Length correction option - setting different length correction.
- High energy correction TPX - using high energy correction for TPX alpha Am241 (200V, peak around 3 MeV).
- Err - missing param file - error behavior without parameters file
- Err - empty file - error behavior with empty file file
- Err - simple elist - error behavior with simple elist
- Err - nonsense file - error behavior with nonsense file
- Err - empty clog - error behavior with empty clog
- Err - incorrect clog - error behavior with incorrect clog
- Err - partial elist - error behavior with partial elist
- Err - dsc file - error behavior with dsc file
- Err - matrix toa - error behavior with matrix toa
- Err - module configs without path - error behavior of module configs wrong path
All examples were produced on linux system but they should be also valid for windows system.
Note that some of the examples were written with older syntax from previous versions.
This syntax is still valid, but not stated in the list below and it should be easy to follow which names were replaced with the new ones (e.g. CalMatrix_Path to CalMat).
Issues and Help
Any ideas or issues can be reported on github forum for issues:
https://github.com/lmareksla/DPE_Issues/
or send directly on the first mail below.
For more information or help: lukas.marek@advacam.cz, carlos.granja@advacam.cz.
Main Configuration of DPE
The main settings of the DPE is done through the so-called parameters file which is a text file with all paths, names and switches to configure the DPE. It has similar notation as the c++ code and example can be found in the directory with the program.
There are several other features of DPE which are handled with different configuration files e.g. filtering, histograms etc. Their settings is mentioned in the next sections.
Syntax of the Parameters File
- Everything which is after
//
is taken as a comment and the program ignores it (gives opportunity to create comments or easily test different settings). - Every string should be enclosed in quotes
"
to correctly account for white spaces in the string. Example:FileName = "file name.txt"
. - Number can be written in the common notations, e.g. 2e5, 50, -666.12, .1 (same as 0.1) . Example:
Number = 1e7
. - If an array of values is needed then it is separated with comma symbol:
,
(white spaces are unimportant). Example:ArrayNumbers = 12,13, 14, 15
,ArrayString = "hop", "skok"
- The syntax is partially based on the c++ so it can be used in some highlighting text editors for better clarity.
- All used paths are checked for proper slashes of paths which means that same configuration files can be used on windows and linux. This is not valid for the input argument into DPE = path and name of the configuration file.
- It is also possible to convert all parameters to snake case:
CalMat
tocal_mat
(only those which are listed below, it is not compatible with older versions of DPE).
List of Parameters
Name | Type | Description | Example | Default |
---|---|---|---|---|
CalMat | STRING | Path to a directory with calibration matrices if data should be calibrated during pre-processing. | CalMat="/Path/To/Cal/Matrices/" |
empty string |
ClogOutName | STRING | Name of the cluster log file for export without suffix/ending. If it is stated in the parameters file then the clog is created (alias DoCreateClog = true). | ClogOutName="ClusterLog" |
ClusterLog |
ClusterCount | INT | Count of clusters which should be processed and evaluated. | ClusterCount=1 ClusterCount=10000 ClusterCount=5 | |
ClusterCountInSensImage | INT | Count of clusters which will be exported in the integrated sensor plots (also per class). It can be less, if there is no such a count of needed clusters set by this value. | ClusterCountInSensImage=10 ClusterCountInSensImage=1000 |
100 |
ClustererName | STRING | Name o the clusterer binary. Default, it is decided based on the operation system. | ClustererName="clusterer" ClustererName="clusterer.exe" |
clusterer (linux) or clusterer.exe (windows) |
ClustererPath | STRING | Path to the clusterer program (part of the download package). | ClustererPath="/Path/To/Clusterer/ Program/" |
./Clusterer/ |
ClusterImageCount | INT | Count of clusters which will be exported as individual cluster plots. | ClusterImageCount=20 ClusterImageCount=1 |
15 |
CoincEventListBaseName | STRING | Base name of the coinc_event list which is used for txt and also for the json format. It is combined as this name and adding .txt and .json. | CoincEventListBaseName="CoincEventList" CoincEventListBaseName="List" CoincEventListBaseName="Name" |
CoincEventList |
CoincEventListJsonName | STRING | Name of output coinc_event list in json format. | CoincEventListJsonName="SampligList.json" |
CoincEventList.json |
CoincEventListName | STRING | Name of the output coinc_event list. | CoincEventListName="CoincEventList.txt" |
CoincEventList.txt |
ComptonListBaseName | STRING | Base name of the compton list which is used for txt and also for the json format. It is combined as this name and adding .txt and .json. | ComptonListBaseName="ComptonList" ComptonListBaseName="List" ComptonListBaseName="Name" |
ComptonList |
ComptonListJsonName | STRING | Name of output compton list in json format. | ComptonListJsonName="SampligList.json" |
ComptonList.json |
ComptonListName | STRING | Name of the output compton list. | ComptonListName="ComptonList.txt" |
ComptonList.txt |
ComptonName | STRING | Name of the file with Compton direction reconstruction configuration. | ComptonName="Compton.ini" |
empty string |
ComptonPath | STRING | Path to the file with Compton direction reconstruction configuration. | ComptonPath="/Path/To/Compton/Config/" |
empty string |
DetGeoName | STRING | Name of the file with detector geometry configuration. | DetGeoName="DetGeo.ini" |
empty string |
DetGeoPath | STRING | Path to the file with detector geometry configuration. | DetGeoPath="/Path/To/DetGeo/Config/" |
empty string |
DirectionListBaseName | STRING | Base name of the direction list which is used for txt and also for the json format. It is combined as this name and adding .txt and .json. | DirectionListBaseName="DirectionList" DirectionListBaseName="List" DirectionListBaseName="Name" |
DirectionList |
DirectionListJsonName | STRING | Name of output direction list in json format. | DirectionListJsonName="SampligList.json" |
DirectionList.json |
DirectionListName | STRING | Name of the output direction list. | DirectionListName="DirectionList.txt" |
NClusters |
DirectionName | STRING | Name of the file with directional analysis configuration. | DirectionName="Direction.ini" |
empty string |
DirectionPath | STRING | Path to the file with directional analysis configuration. | DirectionPath="/Path/To/Direction/Config/" |
empty string |
DoClustererLog | BOOL | Export clusterer log. Only valid if clusterer is used during processing (raw data). | DoClustererLog=true DoClustererLog=false |
true (used) |
DoCompton | BOOL | Control whether Compton directional reconstruction should be done. There is limited list of cases when this analysis is usable based on the detector settings (see section below). | DoCompton=true DoCompton=false |
true (if configuration allows) |
DoCreateClog | BOOL | Control for a generation of the clog. | DoCreateClog=true DoCreateClog=false |
true (used) |
DoCreateElistExt | BOOL | Switch to produce elist with additional columns from processing output: filter pass or not passed (1,0), PID class (e.g. from -1 to 8). It is placed as the last columns of the elist with name `FilterPass, PIDClass`. | DoCreateElistExt=FileIn_Count |
true (used) |
DoDirection | BOOL | Control whether Compton directional reconstruction should be done. There is limited list of cases when this analysis is usable based on the detector settings (see section below). | DoDirection=true DoDirection=false |
true |
DoDirectionTrackCond | BOOL | Control to use tracking condition during directional analysis (additional constrain on features of clusters). | DoDirectionTrackCond=true DoDirectionTrackCond=false |
false |
DoDoseUnitRadiology | BOOL | Switch to change from the dose rate in uGy/h to Gy/s. | DoDoseUnitRadiology=true DoDoseUnitRadiology=false |
false |
DoExportElistExtFilter | BOOL | Control whether clusters which were filtered out should not be exported into extended elist. | DoExportElistExtFilter=true DoExportElistExtFilter=false |
false (exported) |
DoExportElistExtSensEdge | BOOL | Control whether clusters which were at sensor edge should not be exported into extended elist. | DoExportElistExtSensEdge=true DoExportElistExtSensEdge=false |
false (exported) |
DoExportGraphics | BOOL | Control for an export of graphical output. | DoExportGraphics=true DoExportGraphics=false |
true (used) |
DoExportText | BOOL | Control of an export of the output into files. If turned off, then no export into files is done. | DoExportText=true DoExportText=false |
true (used) |
DoFilter | BOOL | Control whether filtering should be done during processing. | DoFilter=true DoFilter=false |
true (used) |
DoFrame | BOOL | Control to do frame analysis during processing. | DoFrame=true DoFrame=false |
true (used) |
DoGraphicsMultiThread | BOOL | Switch to use multiprocessing/multitreading during graphics export. | DoGraphicsMultiThread=true DoGraphicsMultiThread=false |
true (used) |
DoCheckPrereq | BOOL | Control to do or not do check of prerequisites as clusterer or python modules. | DoCheckPrereq=true DoCheckPrereq=false |
true |
DoIgnoreClusterSensEdge | BOOL | Switch to ignore cluster at sensor edge. They are not used into histograms, physical products etc. (evaluation) but they remain in `Files` from clusterer. | DoIgnoreClusterSensEdge=true DoIgnoreClusterSensEdge=false |
false |
DoLog | BOOL | Control for creating log into file. | DoLog=true DoLog=false |
true (used) |
DoLogDateTime | BOOL | Control whether date and time should be in the log file. | DoLogDateTime=true DoLogDateTime=false |
true (used) |
DoMaskFullPrint | BOOL | Switch to fully print the mask content into terminal and log file. | DoMaskFullPrint=true DoMaskFullPrint=false |
false |
DoMultiFile | BOOL | Check whether the multi file processing should be used (can be turned off/false and on/true). | DoMultiFile=true DoMultiFile=false |
false |
DoPID | BOOL | Switch to turn off the particle identification | DoPID=true DoPID=false |
true (used) |
DoPrint | BOOL | Control for print progress into terminal. | DoPrint=true DoPrint=false |
true (used) |
DoRadFieldRecog | BOOL | Control whether radiation filed recognition should be done. There is limited list of cases when this analysis is usable based on the detector settings (see section below). | DoRadFieldRecog=true DoRadFieldRecog=false |
true (if configuration allows) |
DoRemoveElist | BOOL | Switch to delete elist created by the clusterer but the extended elist is not effected by this switch. If set to true, the elist is removed. It set to true because the information is doubles and extended in the elist extended. | DoRemoveElist=true DoRemoveElist=false |
true (used) |
DoRemoveOldElist | BOOL | Switch to explicitly delete or not an already existing elist from a previous processing. If set to true, the elist is removed and a new one is created. If false, the DPE checks elist existence and if there is a elist, the pre-processing with clusterer is skipped. The default value is true because this could cause an error in processing if mask is used with false option -> it wont be processed again and the mask will be ignored for elist data. | DoRemoveOldElist=true DoRemoveOldElist=false |
true (used) |
DoRunClusterer | BOOL | Switch to skip processing of raw data with clusterer. | DoRunClusterer=true DoRunClusterer=false |
true (used) |
DoSensMatrixCount | BOOL | Switch to use counts in sensor plots instead of energy. | DoSensMatrixCount=true DoSensMatrixCount=false |
false |
DoSigVec | BOOL | Control whether significant vector creation should be done. | DoSigVec=true DoSigVec=false |
true |
DoTpxHighEnergyCorr | BOOL | Switch to use TPX high energy correction. | DoTpxHighEnergyCorr=true DoTpxHighEnergyCorr=false |
false |
ElistExtOutName | STRING | Output name of the elist extended which also includes information as PID class etc. Ending/suffix can be stated and is used in this case (in contrast to ElistOutName). | ElistExtOutName="ElistExt.advelist" ElistExtOutName="EExt.advelist" |
EventListExt.advelist |
ElistOutName | STRING | Name of the elist file for export without the suffix. | ElistOutName="EventList" |
EventList |
FileInCount | INT | Count of input files which should be processed if multi file processing is used. | FileInCount=10 | |
FileInName | VSTRING | Name of the input file. Spaces can be part of the name. | FileInName=tot toa.t3pa |
empty string |
FileInNameEnd | STRING | Name end/extension of the input file. If the file name is `FileIn.txt` then name end is `.txt`. It is used for a several file processing. | FileInNameEnd=".txt" |
empty string |
FileInPath | STRING | Path of the input file(s). Directory can be also given, then processing of files in this dir is done based on valid type of files. | FileInPath="/Path/ To/ File/" |
empty string |
FileOutName | STRING | Name of the output sampling list. | FileOutName="SamplingList.txt" |
SamplingList.txt |
FileOutPath | STRING | Path for the results. All results will be exported into this directory and if the path/dir exists and it can be created. | FileOutPath="/Path/File /Out/" |
current working directory |
FilterName | STRING | Name of the file with filter configuration. | FilterName="Filter.ini" |
empty string |
FilterPath | STRING | Path to the file with filter configuration. | FilterPath="/Path/To/Filter/Config/" |
empty string |
FrameListBaseName | STRING | Base name of the frame list which is used for txt and also for the json format. It is combined as this name and adding .txt and .json. | FrameListBaseName="FrameList" FrameListBaseName="List" FrameListBaseName="Name" |
FrameList |
FrameListJsonName | STRING | Name of output frame list in json format. | FrameListJsonName="SampligList.json" |
FrameList.json |
FrameListName | STRING | Name of the output frame list. | FrameListName="FrameList.txt" |
FrameList.txt |
GraphicsMultiThreadFileSize | FLOAT | File size limit in B on single multi processing batch. It should be decreased if the RAM is overwhelmed during processing. | GraphicsMultiThreadFileSize=2e6 GraphicsMultiThreadFileSize=1e6 |
3e6 (30 MB) |
Hist1DGraphicsRebin | INT | Rebin value for 1D histograms before plotting. It will rebin all 1D histograms to make them more clear and exporting graphics faster if value around 1e2 is used. If value 0 is set then this option is skipped. This option overwrite a possible user settings in the configuration files for histograms, but it is valid only for exported graphics, not for text files. | Hist1DGraphicsRebin=100 Hist1DGraphicsRebin=2 |
0 |
HistName | STRING | Name of the file with significant vector configuration. | HistName="Hist.ini" |
empty string |
HistPath | STRING | Path to the file with significant vector configuration. | HistPath="/Path/To/Hist/Config/" |
empty string |
ChipType | STRING | Type of the chip (TPX, TPX2, TPX3). Example: `ChipType = "TPX3"` to specify that TPX3 was used chip. | ChipType="TPX" ChipType="TPX2" ChipType="TPX3" ChipType="TPX4" |
TPX |
IsEnergyCalibrated | BOOL | Control to inform DPE that data are already energy calibrated. | IsEnergyCalibrated=true IsEnergyCalibrated=false |
false (uncalibrated) |
LengthCorr | INT | Num switch to choose formula for calculation of cluster length in sensor plane. Possible values:
* 1 - based on weighted standard deviation subtraction from projected length. * 2 - based on projected width subtraction from projected length. * 3 - based on model of Nabha. |
LengthCorr=1 LengthCorr=2 LengthCorr=3 |
1 |
LogName | STRING | Name of the log file. | LogName="Log.txt" LogName="file.hop" |
Log.txt |
LogPath | STRING | Path to the log file. | LogPath="/Path/To/Log/ File/" |
current working directory |
MaskName | STRING | Name of the file with a mask configuration. | MaskName="Mask.txt" |
empty string |
MaskPath | STRING | Path to the file with a mask configuration. | MaskPath="/Path/To/Mask/Config/" |
empty string |
MeasMode | STRING | Measurement mode of detector = itot+count (itot+count). In the current version, it is used for itot+count mode recognition. | MeasMode="itot+count" |
empty string |
ModelPath | STRING | Path to models. | ModelPath="/Path/To/Models/" |
./Models/ |
PIDAlg | INT | Switch to change between different PID algorithms. All possibilities (101,102,103,201,202):
* 101 - 4 class heuristic decision tree, TPX 300 um Si * 102 - 9 class heuristic decision tree, TPX 300 um Si * 103 - 16 class heuristic decision tree, TPX3 500 um Si * ...171-173 new filters for neutrons * 201 - 3 class DNN, TPX 300 um Si 30 V * 202 - 6 class DNN, TPX 300 um Si 30 V * 251 - 3 class DNN, TPX3 500 um Si 80 V * 252 - 6 class DNN, TPX3 500 um Si 80 V |
PIDAlg=252 | |
RadFieldRecogDatabasePath | STRING | Path to database for radiation filed recognition via comparator. This database should include significant vectors which are used for comparison. | RadFieldRecogDatabasePath="/Path/To/Comp/DB/" |
empty string (default database) |
RadFieldRecogListBaseName | STRING | Base name of the rfr list which is used for txt and also for the json format. It is combined as this name and adding .txt and .json. | RadFieldRecogListBaseName="RadFieldRecogList" RadFieldRecogListBaseName="List" RadFieldRecogListBaseName="Name" |
RadFieldRecogList |
RadFieldRecogListJsonName | STRING | Name of output rfr list in json format. | RadFieldRecogListJsonName="SampligList.json" |
RadFieldRecogList.json |
RadFieldRecogListName | STRING | Name of the output rfr list. | RadFieldRecogListName="RadFieldRecogList.txt" |
RadFieldRecogList.txt |
RadFiledRecogSigVecBaseName | STRING | Base name (beginning) of the significant vectors in the user database for RFR via comparator. DPE only searches such a files which includes this string and loads them into database. Only valid if RadFieldRecogDatabasePath is used. | RadFiledRecogSigVecBaseName="vec" RadFiledRecogSigVecBaseName="sig_vec" |
SigVec |
RadFiledRecogSigVecEndData | STRING | End/extension/suffix of the data files of significant vectors in the user database for RFR via comparator. DPE only searches such a files which includes this string and loads them into database. Only valid if RadFieldRecogDatabasePath is used. | RadFiledRecogSigVecEndData=".vec" RadFiledRecogSigVecEndData=".data" RadFiledRecogSigVecEndData=".txt" |
.vec |
RadFiledRecogSigVecEndInfo | STRING | End/extension/suffix of the info files of significant vectors in the user database for RFR via comparator. DPE only searches such a files which includes this string and loads them into database. Only valid if RadFieldRecogDatabasePath is used. | RadFiledRecogSigVecEndInfo=".info" RadFiledRecogSigVecEndInfo=".vec_info" |
.vec_info |
SamplingListBaseName | STRING | Base name of the sampling list which is used for txt and also for the json format. It is combined as this name and adding .txt and .json. | SamplingListBaseName="SamplingList" SamplingListBaseName="List" SamplingListBaseName="Name" |
SamplingList |
SamplingListJsonName | STRING | Name of output sampling list in json format. | SamplingListJsonName="SampligList.json" |
SamplingList.json |
SamplingListJsonPath | STRING | Path for output of sampling list in json format. | SamplingListJsonPath="/Path/File JSON/Out/" |
current working directory |
SamplingListName | STRING | Name of the output sampling list. | SamplingListName="SamplingList.txt" |
SamplingList.txt |
SensBias | FLOAT | Sensor Bias in volts. | SensBias=-50 SensBias=200 SensBias=45.234 |
50 |
SensDens | FLOAT | Density of the sensor in g/cm3 for temperature at 20 deg. If not set then its values is deduced based on the type of material (only for those mentioned in SensMat). | SensDens=2.93 | |
SensHeight | INT | Height of the sensor in count of pixels. | SensHeight=256 SensHeight=512 SensHeight=3 |
256 |
SensMat | STRING | Sensor material (Si, CdTe, GaAs). Example: | SensMat="Si" |
Si |
SensThick | FLOAT | Sensor thickness in micrometers. | SensThick=300 SensThick=500 |
300 |
SensWidth | INT | Width of the sensor in count of pixels. | SensWidth=256 SensWidth=512 SensWidth=3 |
256 |
SigVecName | STRING | Name of the file with significant vector configuration. | SigVecName="SigVec.ini" |
empty string |
SigVecPath | STRING | Path to the file with significant vector configuration. | SigVecPath="/Path/To/SigVec/Config/" |
empty string |
TimeAcq | FLOAT | Measurement/data acquisition time which can be set to use it instead of elapsed time evaluated in engine run. It is given in seconds. | TimeAcq=3 TimeAcq=1e2 | |
TimeCoinc | FLOAT | Time in nanoseconds for evaluations coincidence events -> clusters whose times are in between this this interval are grouped as coincidence group (first cluster time is down edge and plus coince group time is upper edge). | TimeCoinc=1e2 TimeCoinc=10 TimeCoinc=1.2554 |
100 |
TimeFrameAcq | FLOAT | Acquisition time of frame (if frame are processed). DPE also tries to make an estimation if the information is included in the input data (clog etc.). | TimeFrameAcq=1e-3 TimeFrameAcq=1 |
0 (deduce from data) |
TimeFrameDead | FLOAT | Dead time of frame (if frame are processed). DPE also tries to make an estimation if the information is included in the input data (clog etc.). | TimeFrameDead=1e-3 TimeFrameDead=2 |
0 (deduce from data) |
TimeSampling | FLOAT | Time in seconds for individual evaluations of the data sample -> the whole data sample and its time of collection is divided into sampling intervals with duration of the TimeSampling. | TimeSampling=1e-6 TimeSampling=0.00001 TimeSampling=233.2 |
1 |
DoMask | BOOL | Control to turn off masking if mask file is given. | DoMask=true DoMask=false |
false |
Detector Settings
It is possible to set the detector configuration with following list of parameters in the main parameters file:
- SensMat - sensor material (Si, CdTe, GaAs etc). Example:
SensMat = "Si"
to specify that the sensor material is silicon. Possible values:Si
- for silicon sensor. Used density: 2.329 g/cm-3.CdTe
- for cadmium teluride sensor. Used density: 5.85 g/cm-3.GaAs
- for galium arsenide. Used density: 5.318 g/cm-3.
- SensThick - Sensor thickness in micrometers. Example:
SensThick = 500
to specify the thickness of the detector as 500 micrometers. This is continues variable. - SensBias - sensor bias in volts. Example:
SensBias = 100
to specify that applied was 100 V. This is a continues variable. - ChipType - type of the chip (TPX or TPX3). Example:
ChipType = "TPX3"
to specify that TPX3 was used chip. Possible values:TPX
- for timepix chipTPX3
- for timepix 3 chipTPX2
- for timepix 2 chip
- SensDens - set density of used material in g/cm-3. DPE has default values for material specified in
SensMat
.
These values are used, for example, in the calculation of dose rate because there is a dependence on the sensor material and thickness. It also determines which PID algorithm should be used.
Input
Single File Processing
To process specific file with DPE, its name FILE_IN_NAME
and path /PATH/TO /FILE/
must be specified via following parameters in parameters file:
FileInName = "FILE_IN_NAME" FIleInPath = "/PATH/TO /FILE/"
This will cause that DPE will try to find this file and process it. Error is produced if the file can not be found and the processing is aborted.
Supported File Formats
The current version of the DPE engine is capable to process following files:
- Data driven files - t3pa, t3p, t3r
- Cluster log files- calibrated/uncalibrated
- Elist - only so-called full elist which includes also the header with cluster variables name.
- Frame formats
If a clog from itot+count measurement is used then it is needed to specify this in the main config file with following option:
MeasMode = "itot+count"
Supported suffixes are following:
".t3pa", ".t3p", ".t3r", ".txt", ".pmf", ".plog", ".bmf", ".clog", ".elist", ".advelist", ".extelist"
Currently, the only toa formats are unsupported with undefined results deom DPE (time is handled as energy information).
Examples of most of the suported files can be found in the Examples/Test
directory (list of examples is in the section First Example of Run
)
Multi File/Batch Processing
It is possible to process several files in one run of DPE. There are several ways how to specify which files should be processed: * Process whole directory specifying path to the directory. * Process given files specifying their names in FileInName
as an array. * Process only those files in directory which includes given strings in their names. Examples of such a processings can be found below.
Lets assume that our directory at path: /PATH/TO/DIR/
includes the following files: * File_001.t3pa * File_002.t3pa * File_003.t3pa * File.t3pa * File_2.clog
In the first case, when a whole directory should be processed, only the path to the directory is specified in the parameters file:
FileInPath = "`/PATH/TO/DIR/"
DPE checks the existence of the directory and finds all files. It processes only those which have the same suffix/ending, and DPE found them first.
For this specific directory, let’s assume that File_001.t3pa
was the first found by DPE, then all files ending with .t3pa
are processed.
The file File_2.clog
is not processed, because of the different suffix.
In the second case, when names of files are specified, it is needed to specify the names into parameter FileInName
:
FileInName = "File_001.t3pa", "File_003.t3pa", "File.t3pa"
DPE will process files: File_001.t3pa
, File_003.t3pa
and File.t3pa
. Files with different suffixes should not be mixed! (e.g. clog and t3pa, DPE handles these files differently, which might cause unexpected behavior of the engine.)
In the last case, when user wants to specify only parts of the names, it is possible to specify beginning and ending of the name:
FileInName = "File_" FileInNameEnd = "t3pa"
This will cause that only the following files are processed: File_001.t3pa
, File_002.t3pa
and File_003.t3pa
. If only the first part of the name should be specified, there is need for one additional parameter DoMultiFile
:
FileInName = "File_0" DoMultiFile = true
DPE will again process the same list of files (additional 0 is needed due to the File_2.clog
). In the opposite case when only ending of the name is specified, there is no need to add these parameters:
FileInNameEnd = ".t3pa"
Files File.t3pa
, File_001.t3pa
, File_002.t3pa
and File_003.t3pa
will be processed.
The switch DoMultiFile
can also be used to suppress multi file processing, setting it to false.
Output
Log and Print
During a run of the program, there is an output into the terminal. This include some basic information about the DPE settings and its run.
It can be turned off with DoPrint = false
(default is true). Copy of this output is also given into log file, usually (default) called LogFile.txt
. It can be turned off with DoLog = false
(default is true). The preprocessing run of clusterer also created its own log file called log.txt
where minimal log of clusterer run is stored.
Directory Tree
Results of DPE are exported into several directories separated in most cases based on the type of the analysis. There are some general directories for overall results. The creation of directories depends on used analysis and some of them might not be created if given processing is disabled or unsupported for given settings.
- File - main data files: elist, cluster log and sampling list. It is always created.
- Graph - graphical representation of the sampling list. Only if
DoExportGraphics = true
. - Hist - histograms of cluster variables.
- SpatialMap - spatial maps of cluster variables also with respect to the PID classes.
- SigVec - significant vector which uniquely specifies the given radiation field.
- EventVisual - individual cluster plots, integrated sensor plots and plots of clusters for individual PID classes.
- Direction - results of directional analysis: estimation of radiation direction.
- CoincEvent - results of coincidence analysis: evaluates time correlation between clusters.
- Compton - results of coincidence analysis: evaluates time correlation between clusters.
- RadFieldRecog - results of coincidence analysis: evaluates time correlation between clusters.
These directories can have inner structure based on the specific exports of given analysis (export of graphs in directional analysis is done into directory Direction/Graph
). More details can be found in the sections dedicated to these specifics analysis: directional, Compton, coincidence and radiation field recognition.
The whole export/output can be turned off if DoExportGraphics = false
and DoExportText = false
.
Text basic output
- Elist - file with cluster variables/parameters (this file is in default suppressed with
DoRemoveElist = true
and replaced with extended elist below). - Extended Elist - elist extended with additional columns of PID class (-1 to N-1 of classes where -1 is for others,
PIDClass
) and filter pass (1 = passed, 0 = not passed,FilterPass
) if the filter is used. - Cluster log/clog - detailed list of clusters containing pixelated information. Each line which is starting with
[
includes pixels of a cluster ([X,Y,E,T] = [X position,Y position,energy,time if tpx3 is sued] of a pixel). This file includes all clusters and the mask or filters are not accounted for. - Sampling list - it includes basic information about the radiation field with respect to the observables and their time dependencies.
- Histograms - histograms of given variables - total and also for each PID class.
- Spatial maps - matrices of the sensor filled with integrated cluster variables (energy, size, LET, E/S).
- Significant vector - unique vectors characterizing processed radiation field computed based on histograms.
- Event and sensor visualizations - data files with matrices of exported clusters and sensor.
Each individual analysis also produces some parts of these results and more detailed descriptions are offered in dedicated sections below.
Graphical output
The current version uses a python scripts to convert the text output into graphics:
- Histograms
- Spatial maps
- Individual cluster plots (with values of the cluster variables) and integrated sensor plots
- Graphs of time evolution of physical products
- Directional maps
- Compton directional reconstruction
- …
For this purposes, the matplolib library is used. It can be slow for some configuration (for example high number of bins). To speed up the export of graphics, it is possible to set shown number of bins with following option:
Hist1DGraphicsRebin = 100
which will cause that each histogram will be rebinned to show maximally 100 bins but the exported text files are still with the original number of bins. This option is only functional for 1D histograms.
Another possibility to speed up the plotting is the option DoGraphicsMultiThread
which exploits several CPU threads for graphics creation. This option is by default on/true. The program then uses all threads on the given PC. To disable this feature use:
DoGraphicsMultiThread = false
The multi processing can use up to several GB of ram memory which can be tuned with parameter GraphicsMultiThreadFileSize
. Its default value is 3e6 and lower this number if the ram usage should be decreased (minimal ram usage is defined with needed ram for the biggest plot/histogram which can be of order of GBs).
To produce given plots, a directory .temp
is created in the working directory. All python scripts are saved there before processing.
If the program is terminated just before or in the middle of graphics creation (python run) then a next run will firstly evaluate unfinished plotting and then proceeds with new plots therefore the whole plotting is delayed. In the case of sudden program abortion, it is recommended to check for existence of the .temp
directory and delete it before the next run of DPE.
Clusterer and Preprocessing
The main purpose of this stage is to provide conversion of pixelated data into groups of correlated pixels, clusters. The level of correlation depends on the detector used for data collection. In the case of TPX detector, only spatial correlation is given and it is afterwards used for clusterization process. On contrary for the TPX3, time information is provided along the spatial one which gives opportunity to proceed with more precise clusterization process and eventually avoid additional unwanted effects as pile-up. At this stage main energy per pixel calibration is applied to convert digital ToT information to energy in keV. This process is accompanied with additional calibrations and corrections which aim to maximally supress detector unwanted or anomalous behaviour. The clusterization stage is followed with cluster analysis which evaluates clusters morphological and spectral features.
To summerize the pre-processeing stage:
- Clustering - grouping pixels based on coordinate and also on time information if given
- Calibration and corrections - energy calibration, cluster size correction, XRF peak suppression, TPX high energy correction etc.
- Cluster parameters - calculation of overall cluster parameters (total energy, roundness etc.)
Clustering
The clusterisation converts a frame or a continual stream of received pixels to sets of pixels which belonging to individual particles, clusters. A distinction is needed between two read out modes. In the case of the Timepix, the detector is only capable to measure in the frame mode and obtained data format is in a form of frames. The information is encoded into a sensor matrix of only hit pixels where each of them contain energy deposited during the acquisition time (if a measurement in ToT is performed1). The same mode can be also used with the Timepix3 detector. The advantage with respect to the Timepix is the possibility to measure in the data driven mode. The final data output in this mode resembles the continual stream of pixels where time of arrival is stored together with deposited energy2. This significantly improves the ability for the restoration of the individual particle information.
The algorithm used in the case of the frame mode relies only on the coordinate information. Pixels in a frame which are coordinate neighbours (lies within 8-pixel surroundings) are grouped together and called clusters:
It is also possible to measure in ToA mode, time of arrival or in a counting event mode.
In the case of combined ToT and ToA measurement which is another advantage of the Timepix3 with respect to the Timepix detector.
eq with clustering condition
During this clusterization process, an additional information about the readout can be used. It is known that the chip matrix of pixels is read out column by column which means that the stream of pixels can be efficiently parsed based on the column order. There are several disadvantages and crucial points concerning frame read-out:
- Pile-up: All particles within one frame whose pixels are connecting cannot be separated within the clusterization process and they create the pile-up effect.
- Loss of uncollected charge after the end of frame: It is possible that not all charge has been collected till the end of a frame (mostly some pixels with large deposited energy).
- Dead time for matrix read out: This approximately several of millisecond (more than 10 ms).
These issues and disadvantages motivated development of the data-driven read-out mode. The coordinate information is accompanied by time of arrival which allows to reconstruct events also based on the time information. In this case, the stream of pixels is no longer separated into frames but it is continuously read by the read-out hardware and software.
In the first step, the stream has to time ordered because from the nature of the read-out, it is possible that the pixel stream violates time order. In the next stage, a block of pixels is taken from the full stream/list of pixels based on a time condition:
eq
The parameter should reflect expected time windows in the pixel stream and therefore is disproportional to particle fluence. The value is of order of . In those cases when the condition cannot be fulfilled, a fixed block of of pixels is taken. This number is approximately 100 000 pixels which should statistically minimized the error originating from possible separation of pixels belonging to one particle.
The clusterization process itself separates pixels into clusters based on two main conditions:
eq
These conditions can be translated as a need for pixels to be coordinate neighbours and to have maximal difference in time equal or less than time . The value of this time parameter is of order of tens of nanoseconds and it is dependent on the charge sharing effect and its magnitude.
Calibration and corrections
Energy calibration
After receiving the information from the detector, the possible energy and time of the detector are in corresponding ToT and ToA counts. These variables have to be calibrated and converted to obtain results in energy and time. Time conversion is done with following formula:
eq of toa conversion
The first part of the calibration is the same for the Timepix and Timepix3 ASICs. This work was published by Jakubek[1] and it is based on a combination of a linear and rational calibration function:
eq of calibrationwhere a,b,c and t are parameters obtained from the calibration procedure.
High Energy TPX Calibration
The calibration function has maximum of its validity up to approximately 700 keV for TPX. An additional per pixel calibration process has to be utilized for measurements in which higher per pixels energies are obtained. The additional calibration was intruced as a global model based on the results from the article of Sommer et al [2]. Having the coefficients of the standard calibration and correlation matrix between these standards and the extended ones, it is possible to introduce an estimation of their values.
To use this additional calibration:
DoTpxHighEnergyCorr = true
The resulted histogram of the measured energy can be found in the fugere on the left. Two cases are shown in this image, before (light) and after (dark) an application of the high energy correction. Change in the sensor bias was also introduce during the measurement to demonstrate a different/increasing value of mean per pixel energy. Higher bias in the sensor decreases the effect of charge sharing and increases the per pixel energy and mainly the maximal per pixel energy of cluster. Therefore, it can be seen that in the most disturbed case of 200V the applied correction restores the original information to be in accordance with low bias measurement and their corrected versions.
In the case of TPX3, the energy range is limited with size of the ToT counter which has nominal value 1022. This means that the pixels cannot count more than 1022 counts in contrast to the TPX with counter range of 11810. From the practical point of view, the maximal value which can be measured in one pixel is between 450-500 keV (based on the settings of DACs). The standard energy calibration covers this region and a distortion appears for energy above 2 MeV as the volcano effect.
Cluster Energy Size Based Correction
...
Halo Effect
...
Correction for XRF peak suppression
...
Pile-up Recognition
...
Volcano Effect
... Source[3]
Overshoot
...
Time-walk Correction
... Source[4]
Cluster Variables and Elist
- The file includes clusters/events represented as a set of cluster variables.
- The significant feature of the Elist is a presence of coincidence group expressed as coincidence number and it is adjustable based on the coincidence time interval.
LengthCorrStd | Unit | Column | Description |
---|---|---|---|
DetectorID
|
- | 1 | It serves to as a unique ID in the cases when multi detector setup is processed. |
EventID
|
- | 2 | ID of events. In the case that several clusters are in the coincidence, then they have the same EventID. It starts from 0 and clusters which are not in coincidences are also accounted in this ID. |
X
|
px | 3 | X coordinate of cluster based on the weighted mean value of X coordinate of pixels, weighted with the energy of pixels. Value ha range from 0.5 to 255.5 (minimal and maximal value of X of pixels). |
Y
|
px | 4 | Y coordinate of cluster based on the weighted mean value of Y coordinate of pixels, weighted with the energy of pixels. Value ha range from 0.5 to 255.5 (minimal and maximal value of Y of pixels). |
E
|
keV | 5 | Energy of cluster, also called volume. It is sum of energies of all cluster pixels. |
T
|
ns | 6 | Time. It is the minimal time of cluster pixels. |
Flags
|
- | 7 | Free column for user needs. In the case of frame data, it includes frame number. |
Size
|
px | 8 | Count of cluster pixels. |
Height
|
keV | 9 | Maximal value of energy of cluster pixels. |
BorderPixCount
|
px | 10 | Count of border pixels in cluster. |
Roundness
|
- | 11 | Morphological feature expressing similarity of cluster to round shape. |
AngleAzim
|
deg | 12 | Estimation of |
Linearity
|
- | 13 | Morphological feature expressing similarity of cluster to linear shape. |
LengthProj
|
px | 14 | Length of cluster based on maximal distance between pixels after projection to the cluster axis. |
WidthProj
|
px | 15 | Length of cluster based on maximal distance between pixels after projection to the axis perpendicular to the cluster axis. |
IsSensEdge
|
- | 16 | Information whether cluster is at sensor edge. 0 for flase, 1 for true (is at sensor edge). |
StdAlong
|
px | 17 | Weighted standard deviation of pixels with respect to cluster axis weighted with pixels energy. |
StdPerp
|
px | 18 | Weighted standard deviation of pixels with respect to axis perpendicular to the cluster axis and weighted with pixels energy. |
Thin
|
- | 19 | Morphological feature expressing thinness of cluster. |
Thick
|
- | 20 | Morphological feature expressing thickness of cluster. |
CurlyThin
|
- | 21 | Morphological feature expressing combination of cluster thinness and inverse of linearity. |
EpixMean
|
keV | 22 | Mean value of pixels energy in cluster. |
EpixStd
|
keV | 23 | Standard deviation of pixels energy in cluster. |
LengthCorrStd
|
px | 24 | Cluster length in sensor plane corrected for charge sharing based on weighted standard deviation. |
Length3DCorrStd
|
um | 25 | Cluster length in sensor volume corrected for charge sharing based on weighted standard deviation. |
LengthCorrWidth
|
px | 24 | Cluster length in sensor plane corrected for charge sharing based on cluster width. |
Length3DCorrWidth
|
um | 25 | Cluster length in sensor volume corrected for charge sharing based on cluster width. |
LengthCorrNabha
|
px | 24 | Cluster length in sensor plane corrected for charge sharing based on Nabha model. |
Length3DCorrNabha
|
um | 25 | Cluster length in sensor volume corrected for charge sharing based on Nabha model. |
AngleElev
|
deg | 26 | Elevation angle of cluster. |
LET
|
keV/um | 27 | Linear energy transfer based on the corrected length and energy (E). |
Diameter
|
px | 28 | ... |
PIDClass
|
- | 29 | If PID is used, this number expressed class into which the cluster was classified. It goes from 0 to count_class - 1. |
FilterPass
|
- | 30 | If filter is used, then it includes information whether cluster passed given filters or did not. 1 for passed, 0 for did NOT pass. |
Example of created elist:
Cluster log
Clusters can be exported as individuals pixels. This format is called cluster log or clog and it is highly dependent on the input data (detector, measured mode, etc.). In the case of TPX detector the format is following:
Frame 1 (UNIX_TIME, ACQ_TIME s) [X_1, Y_1, iToT_1/E_1/] [X_2, Y_2, iToT_2/E_2/]... .... Frame 2 (UNIX_TIME, ACQ_TIME s) [X_1, Y_1, iToT_1/E_1] [X_2, Y_2, iToT_2/E_2]... ....
where ACQ_TIME
is acquisition time of frames and UNIX_TIME
is current UNIX time stamp in seconds (both). Every line starting with [
is one cluster and in each square brackets is one pixel and its information [x coordinate,y coordinate, iToT/Count/Energy]
. The ToA and Count modes are not correctly implemented and they are treated as the ToT mode (especially of importance for the cluster parameters in the elist).
The format is different for the TPX3 detector:
Frame 1 (FIRST_TIME_COINC_GROUP, COINC_TIME_WINDOW ns) [X_1 ,Y_1, ToT_1/E_1, T_Diff_1] [X_2, Y_2, ToT_2/E_2, T_Diff_2]... [X_1 ,Y_1, ToT_1/E_1, T_Diff_1] [X_2, Y_2, ToT_2/E_2, T_Diff_2]... .... Frame 2 (FIRST_TIME_COINC_GROUP, COINC_TIME_WINDOW ns) [X_1 ,Y_1, ToT_1/E_1, T_Diff_1] [X_2, Y_2, ToT_2/E_2, T_Diff_2]... [X_1 ,Y_1, ToT_1/E_1, T_Diff_1] [X_2, Y_2, ToT_2/E_2, T_Diff_2]... ....
where each frame in this case coincidence group = clusters whose min times are within the COINC_TIME_WINDOW
starting with the first one at FIRST_TIME_COINC_GROUP
. The pixels time is written as difference with the FIRST_TIME_COINC_GROUP
: T_Diff = T_Pix - FIRST_TIME_COINC_GROUP
where T_Pix
is original time of pixels in ns. An exception from this rule is measurement in the frame mode iToT+Count when data are already saved in cluster log within the Pixet. To properly evaluate this data with the clusterer, it is needed to add additional option into the command line:
--measmode "itot+count"
This data does not include additional ToA information therefore the format has the same meaning as in the first case with TPX:
Frame 1 (UNIX_TIME, ACQ_TIME s) [X_1 ,Y_1, iToT_1/E_1, Count_1] [X_2, Y_2, iToT_2/E_2, Count_2]... [X_1 ,Y_1, iToT_1/E_1, Count_1] [X_2, Y_2, iToT_2/E_2, Count_2]... .... Frame 2 (UNIX_TIME, ACQ_TIME s) [X_1 ,Y_1, iToT_1/E_1, Count_1] [X_2, Y_2, iToT_2/E_2, Count_2]... [X_1 ,Y_1, iToT_1/E_1, Count_1] [X_2, Y_2, iToT_2/E_2, Count_2]... ....
Summary of supported input data formats: * TPX3, data driven mode, tot+toa * TPX3, frame mode, itot+count * TPX, frame mode, tot Every other format will be processed but the results might be uncorrected evaluated, for example in the cluster parameters.
Reprocessing of cluster logs can create cluster logs which could have incorrect coincidence group (in the case of TPX3 with ToA). If the coincidence window el-coinctoadiff
is increased with respect to the original one then the resulted group will not integrate several groups together which would fulfill this new time condition.
Masking
There is a possibility to mask a part of the sensor or individual pixels in the case of t3pa and clog files.
These pixels are omitted in the pre-processing.
This can be used to avoid a signal from noisy pixels or to focus on some more interesting part of the sensor.
The configuration is done through a configuration file. An example can be found in the program directory.
If the masking is used the during the pre-processing an additional t3pa/clog file is created with masked pixels according to the user configuration in the export directory (starts with MASK_+ FileInName
).
Inclusion of Mask.ini into Parametrs File
The configuration file of the mask can be included with following options in the ParametersFile:
MaskName = "Mask.txt" // Name of the INI configuration file of the mask. MaskPath = "PATP/TO/MASK/" // Path of the INI configuration file of the mask.
Syntax of Configuration File
- The pixel which should be masked is written in format [X,Y]. These are X and Y coordinates of the pixel where both of them are integers from 0 to 255. The sensor orientation is usual and X coordinate is for rows and Y is for lines. They can be written in lines or rows in the mask file but always the [X,Y] in one line (can not be separated over several lines)
- If a larger part is of interest, following notation is used: [X_min - X_max, Y_min - Y_max] where X_min is the minimum coordinates of a pixel, X_max is the maximum coordinate of a pixel etc. For example, [0-255,0-120] masks almost half of the sensor - on the X axis from 0 to 255 and on the Y axis from 0 to 120. A simpler demonstration can be masking of one line of pixels (11th) : [0-255, 10].
Example of all possibilities for masking:
[0,1] [2,5] [2,3] [2-43,5] [76,8-90] [0-50,50-60]
Masking with Pixet Mask
It is also possible to exploit the masking created in Pixet, which takes the form of a text file where the lines and rows correspond to those of the detector. The orientation in this case is the same for x axis, but it is reversed for y axis where bottom of the detector is first line and top part of the detector is last line. The delimiter in this case is simple white space . The masked pixels are marked as
0
and unmasked are 1
. See following example (...
used as abbreviations):
1 1 1 1 0 1 ... 1 1 1 1 1 1 1 ... 1 ... 0 1 1 1 1 1 ... 1
Statistical Information
Statistical information can be found in the general sampling list and it is also printed into log.
Example for the sampling list in json format (explanations are given in the comments after //
):
"CountPixHit_Sum_MaskOk_cnt":242461, // Count of pixels which passed the mask and used for processing.
"CountPixHit_Sum_MaskOk_perc":99.23, // The same as above but as part/percentage to total amount of pixels.
"CountPixHit_Sum_MaskOmit_cnt":1872, // Count of pixels which did NOT pass the mask and they are used for processing.
"CountPixHit_Sum_MaskOmit_perc":0.766, // The same as above but as part/percentage to total amount of pixels.
"CountPixHit_Sum_All_cnt":41935, // Total count of loaded pixels for raw data.
Example for log file:
Mask conf file - name: Mask.txt Mask conf file - path: ./Test/data/test_029/ ... -------------------------------------------------- SENSOR MASK ...has been loaded. Count of masked pixels: 901 (1.37%) -------------------------------------------------- ... Processing RawData T3PA: Creating masked raw file: MASK_tot_toa.t3pa ... -------------------------------------------------- MASK -------------------------------------------------- CountPixHit_Sum_MaskOk [-]: 242461 CountPixHit_Sum_MaskOk [%]: 99.23 CountPixHit_Sum_MaskOmit [-]: 1872 CountPixHit_Sum_MaskOmit [%]: 0.766 CountPixHit_Sum_All [-]: 41935
First part shows from which destination is the mask taken and what is the name of the file.
The second part informs about successful loading of the file and how many pixels are masked.
The third part appears in the preprocessing sections and informs about actual masking of the raw data.
The last part is in the sections with results and includes the same information as it is in the sampling file.
Filtering
During the processing, filters can be used to obtain only information about particles of interest.
The filters are applied on cluster variables/parameters level (e.g. energy, height etc.). It is specified trough a configuration file in the INI format.
The cluster variables which should be used for filtering are specified with their unique name which is included
in the header of the created elist (if elist is input then it has to be already part of the file).
All results produced by DPE are only for filtered particles/for those which passed filter (histograms, graphs, spatial maps etc.).
An example can be found in the program directory.
Inclusion of Filter.ini into Parameters File
The configuration file of the filter can be included with following options in the ParametersFile:
FilterName = "Filter.ini" // Name of the INI configuration file of the filter. FilterPath = "/PATH/TO/FILTER/" // Path of the INI configuration file of the filter.
Syntax of Configuration File
The configuration file in INI format is a set of sections where each section is dedicated to one cluster parameter for which filtering should be done.
Individual parts of these sections are specific values of the filters for the given cluster parameter. Example with filter conditions for energy of clusters alias row E
in elist:
- One section of the configuration file is named based on the key
E
in the elist:[E]
. - Condition is then written as
Range=100,200
. The nameRange
has to be first part of the condition name. This settings will produce a filter on cluster energy which should be only from 100 to 200 keV (edges are included). - More ranges can be specified for one parameter. Individual conditions/ranges have to always include string
Range
in names, but they have to differ, therefore suffixes/endings have to be introduced:Range_1
,Range_2
etc. (, see example below).
Filter with single condition in energy E
from 100 to 200 keV:
[E] Range = 100, 200
Filter with multiple conditions on energy E
:
[E] Range_1=100,200 Range_2=500,1000 Range_asdasd=300,2000
Names of Cluster Parameters/Variables
The names can be found in the table about cluster variables in the section Cluster Variables and Elist.
Output into Extended Elist
The output of the filtering can be found in the extended elist (in directory File
). New column is created with name FilterPass
.
The values of this new clomun/parameter are 1
for passing the filter and 0
for NOT passing the filter.
All particles are stored in the extended elist, but it is possible to suppress the export of those particles which did not pass the filter conditions:
DoExportElistExtFilter = false
with this switch in the parameters file, only those particles which passed the filter are exported into extended elist (column FilterPass
is still present).
Statistical and Log Information
Statistical information can be found in the general sampling list and it is also printed into log.
Example for the sampling list in json format (explanations are given in the comments after //
):
"CountParticle_Sum_FilterOk_cnt":4386, // Count of particles which passed the filter.
"CountParticle_Sum_FilterOk_perc":58.93, // Percentage of the above to the total count of particles.
"CountParticle_Sum_FilterOmit_cnt":3057, // Count of particles which did NOT pass the filter.
"CountParticle_Sum_FilterOmit_perc":41.07, // Percentage of the above to the total count of particles.
"CountParticle_Sum_All_cnt":7443, // Count of all particles which were evaluated.
Filter conf file - name: Filter.ini Filter conf file - path: ./Test/data/test_007/ ... -------------------------------------------------- FILTER Conditions are logically connected as AND for different variables (within always as AND) INDEX|NAME CONDITIONS [min|max] 10|Roundness |-----------[0|2] '-----------[0.4|0.45] 9|BorderPixCount |-----------[1|10] '-----------[20|1e+200] 4|E '-----------[100|1e+300] -------------------------------------------------- ... -------------------------------------------------- FILTER -------------------------------------------------- CountParticle_Sum_FilterOk [-]: 4386 CountParticle_Sum_FilterOk [%]: 58.93 CountParticle_Sum_FilterOmit [-]: 3057 CountParticle_Sum_FilterOmit [%]: 41.07 CountParticle_Sum_All [-]: 7443
The first part informing about name and path to the config file of filter can be seen in the first stage of DPE processing.
The second part is in the same processing/nationalization stage of DPE and it shows which filter were recognized and what ranges are about to be used.
The last part informs about statistical information and it includes the same information which are also given in the json sampling list.
Histograms
One of the DPE outputs are histograms of cluster variables/parameters.
The DPE in default exports 1D histograms of all cluster variables in extended elist and also several their combinations as 2D histograms.
It is possible to export user defined 1D and 2D histograms which can be configured with a configuration file in the INI format (an example can be found in the program directory).
The DPE allows to also create histograms of algebraic combinations of the cluster variables (multiplication, division, subtraction, addition).
It is important that the general label (key name of a section in a INI file) of the histogram used in the INI file is unique to each histogram.
If there are more histograms with one common name then the program only updates information about the first one in the INI file.
The label itself in not used in the program itself and title and name of the histogram are set separately.
The histograms are used in other analysis, therefore their changes/using user configuration might disturb e.g. creation significant vectors. This might produce following error for significant vectors:
[ERROR] -1041 : Error occurred during significant vector module initialization. Module will not be used in processing.
It just informs that the settings of histograms is not compatible with settings of significant vectors.
Inclusion of Hist.ini into Parameters File
To use user configuration file for histograms, it is needed to add two parameters into paramters file, name (e.g. Hist.ini
) and path (e.g. PATH/TO/HIST/CONFIG/
) to the configuration file:
HistName - "Hist.ini" // Name of the configuration file for histograms. HistPath = "PATH/TO/HIST/CONFIG/" // Path of the configuration file for histograms.
If at least the name is set and the file is found on the given location then the DPE uses this settings of histograms otherwise default configuration is used.
Histogram 1D
The 1D histograms can be created as fixed bin width histograms or with variable binning (different width of bins).
Lets assume that histograms of size should be created Hist1D_Size
.
The fixed bin width has following settings/needed parameters (example from configuration file with explanations = everything after #):
[Hist1D_Size] VarName="Size" # Name of column with given variable (see the first line in Elist.txt - it has to be the same) Title="S" # Title of histogram (can be arbitrary - X if not given) NBin=9 # Number of bins Xmin=100 # Minimum value (if X = Xmin -> it is NOT included to the first bin - it has to be > Xmin) Xmax=1000 # Maximum value (if X = Xmax -> it is included to last bin NBin)
The same as above but with different choice of variable, it is based on the column position in the elist instead of the cluster variable name:
[Hist1D_Size] ColIndex=7 # Position of the column - starts from 0 (it is Size in the example) Title="S" # -||- NBin=9 # -||- Xmin=100 # -||- Xmax=1000 # -||-
These cases create histogram from cluster size with 9 bins from 100 to 1000 where one bin has width 100. Variable size bin width and its needed parameters (same histogram as the one above):
[Hist1D_Size] VarName="Size" # -||- Title="S" # -||- BinLowEdge=100,200,300,400,500,600,700,800,900,1000 # Low edges of bins with Xmax (Xmin,...,Xmax -> size NBin+1)
Histogram 2D
The 2D histograms are constructed in similar manner as 1D histograms. The fixed bin width histograms and their needed parameters:
[Hist2D_SE] VarName="E","Size" # Name of columns with given variable - 1st is X, 2nd is Y (see the first line in Elist.txt - it has to be the same or see special variables) Title="E,S" # Title of histogram (can be arbitrary - X if not given) NBinX=9 # Number of X bins where X is in this case energy,E Xmin=100 # Minimum value of X (if X = Xmin -> it is NOT included to first bin - it has to be > Xmin) Xmax=1000 # Maximum value of X (if X = Xmax -> it is included to last bin NBin) NBinY=1000 # Number of Y bins Ymin=0 # Minimum value of X (if X = Xmin -> it is NOT included to first bin - it has to be > Xmin) Ymax=1000 # Maximum value of X (if X = Xmax -> it is included to last bin NBin)
An example for cluster variables based on column positions in the elist:
[Hist2D_SE] ColIndex=4,7 # Position of the columns which should be processed - 1st is X, 2nd is Y #...(the same as above)
These cases create 2D histogram of cluster energy and size where energy is from 100 to 1000 with bin width of 100 and size is from 0 to 1000 with bin width of 1. Variable size bin width and its needed parameters:
[Hist2D_SE] VarName="E","Size" # -||- Title="E,S" # -||- BinLowEdgeX=100,200,300,400,500,600,700,800,900,1000 # Low edges of X bins with Xmax (Xmin,...,Xmax -> size NBinX+1) BinLowEdgeY=100,200,300,400,500,600,700,800,900,1000 # Low edges of Y bins with Ymax (Ymin,...,Ymax -> size NBinY+1)
It is also possible do just variable binning in one variable the second one can be with fixed bin size:
[Hist2D_SE] VarName="E","Size" Title="E,S" NBinX=9 Xmin=100 Xmax=1000 BinLowEdgeY=100,200,300,400,500,600,700,800,900,1000
Additional Algebraic Operations
Both 1D and 2D histograms allow additional operations of addition, subtraction, multiplication and division on extracted variables from EList.
Histogram 1D:
ColIndex_Div=4,7 # Division as 5th/7th column (in this case E/S = Epix) - 1st argument is divided by 2nd argument ColIndex_Mult=4,7 # Multiplication as 5th*7th column - 1st argument is multiplied with 2nd argument ColIndex_Add=4,7 # Additions 5th+7th column- 1st argument is added to 2nd argument ColIndex_Subtr=4,7 # Subtraction as 5th-7th column - 2nd argument is subtracted from 1st argument
An example based on cluster variable names:
VarName_Div="E","Size" # Same thing as above but with names of columns VarName_Mult="E","Size" # Same thing as above but with names of columns VarName_Add="E","Size" # Same thing as above but with names of columns VarName_Subtr="E","Size"# Same thing as above but with names of columns
Histogram 2D:
Very similar to 1D histograms but always the first given is X and the second given is Y.
ColIndex=4 # X is set to 4th column - energy ColIndex_Div=4,7 # Y is set to division of 4th/7th
Operations used for both variables:
ColIndex_Div=8,7,4,7 # X is set to division of 8th/7th and Y is set to division of 4th/7th
The same thing can be done with names of columns/variables VarName:
VarName_Div="E","Size","Height","Size" # X is E/Size and Y is Height/Size
Text export
There are two kinds of exported files: histogram data and histogram information/info file.
The histogram data file includes content of bins and bins low edges. The histogram info file comprehends features of the histogram: title, name, count of bins etc. (see more details below).
The data file has a suffix .hist
and the info file .hist_info
. The data file has following formatting for 1D histogram:
Xmin BinCont_1 BinLowEdge_2 BinCont_2 ... .... BinLowEdge_NBin BinCont_NBin Xmax Overflow
The Xmax is included to allow an user easier read of this file without direct need to also read the info file (all needed information is in the data file in the base case).
It is also possible to export only those bins which have non zero content with fixed binning. This can be done with following option in the Hist.ini file:
DoSparseExport=1 # Check whether only nonzero bins should be exported (1 for true and 0 for false).
It has to be used for each histogram separately and it is used as default settings. The exported data file has following format:
Xmin BinCont_1 BinLowEdge_2 BinCont_2 ... ... BinLowEdge_i BinCont_i != 0 ... ... Xmax Overflow
The Xmin,BinLowEdge_2,Xmax
are always exported even if their bin content is 0 for further reading and reconstruction of histograms. This option is not functional for variable binning to avoid a lack of information in the data file for histogram complete reconstruction in post-processing. To reconstruct the histogram, it can be done just based on the data file even in the case of sparse export because the missing bins
and bin width can be calculated based on the first two bins in the data file (BinWidth = BinLowEdge_2 - Xmin).
Anyway, it is recommended to read the info file for example in the case of constant bin width combined with the sparse export to ensure that given
histogram is truly with constant binning and not variable binning (written in the parameter: BinEquiDist=1
= it is const binning and 0 for variable binning).
Similar approach is utilized for the 2D histograms:
Xmin Ymin BinCont_1 # First bin XBinLowEdge_2 Ymin BinCont_2 ... .... .... XBinLowEdge_NBinX Ymin BinCont_NBinX Xmin YBinLowEdge_2 BinCont_NBinX+1 ... .... .... XBinLowEdge_NBinX YBinLowEdge_2 BinCont_NBinX+NBinY XBinLowEdge_2 YBinLowEdge_3 BinCont_NBinX+NBinY+1 ... .... .... XBinLowEdge_NBinX YBinLowEdge_NBinY BinCont_NBinX*NBinY # Last bin Xmax Ymax Overflow
The sparse export has following form.
Xmin Ymin BinCont_1 XBinLowEdge_2 Ymin BinCont_2 ... .... .... XBinLowEdge_i Ymin BinCont_i != 0 ... .... .... Xmin YBinLowEdge_2 BinCont_NBinX+1 ... .... .... XBinLowEdge_m YBinLowEdge_n BinCont_m+n*NBinX != 0 ... .... .... Xmax Ymax Overflow
It is similar to 1D histogram with the exception that there is also the next Y low bin edge to also estimate the width of binning for Y axis.
The info file is formatted as an INI file and it includes the same information as the Hist.ini file for each histogram individually. There two additional features compared with the Hist.ini:
- Statistical information (mean, err of mean, std, err of std)
- Overflow and underflow information
Explanation of individual parameters are written as comment in example below:
[HistPar]
Title="Hist1D_Epix_1" # Title of the histograms used for unique recognition
Name="Mean energy per pixel of cluster" # General name
AxisTitles="Epix [keV/px]","N [-]" # Names of axis titles/lables.
NBin=300 # Count of bins
Xmin=0 # HIstogram minimum
Xmax=300 # Histogram maximum
BinEquidist=1 # Check whether histogram has same width bins
[HistCont]
Underflow=0 # Count all events which are below Xmin
Overflow=1 # Count all events which are above Xmax
[Statistics]
Mean=29.6869 # Estimation of weighted mean value of histogram (weighted with bin contents)
Mean_Err=0.20951 # Error of the mean estimation
Std=13.9273 # Estimation of weighted standard deviation of histogram (weighted with bin contents)
Std_Err=1.40088 # Error of the std estimation
Equivalent information is also given for histogram 2D with appropriate extension of additional dimension.
Graphical export
There is possibility to create plot of Hist1D and Hist2D via python matplotlib, it has to be preinstalled.
The current version will automatically create these plots if DoExportGraphics=true
but it can be turned off - see below.
SOME_PREVIOUS_CODE # This is some code to call histogram 2D Name="Histogram2D title in plot" # Name o histogram in plot AxisTitle="X","Y","N" # Name of axis in the plot
The plots can be with logarithmic scales on all axis:
SOME_PREVIOUS_CODE # This is some code to call histogram 2D DoLogX=1 # 1 for true=do logarithmic X axis and 0 for false DoLogY=1 # 1 for true=do logarithmic Y axis and 0 for false
To turn off plotting:
SOME_PREVIOUS_CODE # This is some code to call histogram 2D DoPlot=0 # Default is 1 - do plot and 0 means not do plot
It is possible to adjust the number binning of the histogram only for the graphics. This can be done with Hist1DGraphicsRebin
in the main configuration file:
Hist1DGraphicsRebin = 100
This option with value 100 will change the shown number of bins in the graphics to 100 but the exported files includes original binning.
This is allowed for faster exports of histograms with more than 10000 bins which can time challenging for a PC.
Log Information
Hist conf file - name: Hist.ini Hist conf file - path: ./Test/data/test_006/ ... -------------------------------------------------- HISTOGRAMS 1D: Hist1D_E_Total from column 4 Hist1D_S_Total from column 7 Hist1D_LET_Total from column 26 Hist1D_EpixMean_Total from column 21 Hist1D_Epix_Total from column-division: 4/7 Hist1D_H_Total from column 8 Hist1D_LET_VarBinWidth_Total from column 26 2D: Hist2D_ES_Total from column 4 7 Hist2D_ES_VarBinWidth_Total from column 4 7 -------------------------------------------------- ... Hist are be shown with max N bins: 100 bins ... Exporting histograms to: ./Test/out/test_006/./Hist/ Exporting histogram 1D: Hist1D_E_1 Exporting histogram 1D: Hist1D_S_1 Exporting histogram 1D: Hist1D_LET_1 Exporting histogram 1D: Hist1D_EpixMean_1 Exporting histogram 1D: Hist1D_Epix_1 Exporting histogram 1D: Hist1D_H_1 Exporting histogram 1D: Hist1D_LET_VarBinWidth_1 ...
The first part informing about name and path to the config file of histograms can be seen in the first stage of DPE processing.
The second part is in the same processing/nationalization stage of DPE and it shows which histograms will be produced.
The last tow parts informs exporting of histograms (first one only if change of hist binning is used in graphical plots).
Significant Vectors
The significant vectors is a unique set/vector of numbers which describes given data sample/radiation field.
It can be used for radiation field recognition as it is done in the DPE.
There is a possibility to create own configuration file for the generation of the significant vector.
The significant vectors are exported into directory SigVec
in path FileOutPath
.
The analysis/generation of SigVec can be suppressed with the following option in Parameters file:
DoSigVec = false
The results are saved into directory SigVec
into files SigVec.vec
and SigVec.vec_info
. First file includes the vector itself and second file includes information about the vector (used histograms, intervals etc.).
Theoretical Introduction
The significant vectors are sets of numbers which should describe uniquely a radiation field. They are created form 1D and 2D histograms based on the following rules (see the illustrations in the ):
- 1D histograms: Each number of the significant vector is produced as a sum of bin contents within some specified ranges. Assuming that some elements of the vector should be done based on the energy distribution/histogram. Ranges of interests have to be specified where an integration is executed. This specification of the ranges is a task for an optimization method and differentiates based on the set of radiotin fields which should be recognized. One number is created for each range. These are normalized with a total integral over whole histogram to achieve comparability with different fields where more or less particles could be measured. A weighting can be also introduced as a last stage of the process to highlight some elements of the vector. This procedure can be done for several other cluster variables resulting in the vector of significant variables.
- 2D histograms: Their conversion into significant variables obey similar rules as the 1D histograms. Except for the ranges, masks can be also specified which determines the area of interest. The integration is based on this mask where the mask has the same features as the histogram itself with values which can be used as weights. The mask values are used as a multiplicative factor in the integration therefore bins with 0 are ignored and bins with value higher than 1 are highlighted.
Inclusion of SigVec.ini into Parameters File
The the DPE will include a user configuration file only if it is written in the PararamtersFile.
Two parameters are used for this purpose (assume that name of config file is SigVec.ini
in path /PATH/TO/SIGVEC/CONFIG/
):
SigVecName = "SigVec.ini" // Name of the configuration file for significant vectors. SigVecPath = "/PATH/TO/SIGVEC/CONFIG/" // Path to the configuration file for significant vectors.
If at least the name is set and the file is found on the given location then the DPE uses this settings of SigVec otherwise default configuration is used.
Histogram 1D - Intervals of Interest
A histogram is examined and all bins (bin contents) with bin center between up and down edge/limit of given interval is summed, number R_1.
This number is used as one element of SigVec = {R_1, R_2, … , R_N} for N of intervals.
At least two numbers has to be given - down and up edge of the interval where interval condition for bins are following: DOWN_EDGE < BinCenter <= UP_EDGE.
An example of configuration file for 1D histogram:
[Var_Epix] # Section name as unique key of the given SigVec element. It has to include the string "Var_".
Title="Epix" # Title of the variable.
HistName="Hist_Epix" # Name of histogram which should be processed - same as name of histogram file without suffix/ending in directory Hist (data suffix)
Intervals=0,200,100,500,500,5000 # Limits of intervals - 1st interval = from 0 to 200, 2nd interval = from 100 to 500 etc. (N in this case is 3)
The intervals can be overlapping. The section name has to include string Var
in its name and they have to differ for different elements of the SigVec: Var_E, Var_S, ...
.
Histogram 2D - Regions of Interest
A histogram is multiplied with mask (from 0 val to inf - can serve as weights) with the same number of bins.
Then the histogram is summed and this number is as an input to the SigVec.
For each given mask one number is produced. An example of configuration file:
[Var_E_S] # Section name as unique key of the given SigVec element. It has to include the string "Var_".
Title="E_S" # Title of the variables
HistName="Hist_E_S_2" # Name of histogram which should be processed - same as name if histogram file without suffix/ending in directory Hist (data suffix)
Mask_1="\PATH\TO\MASK\MASK1_FILE_NAME" # Path to masks which should be applied - same dimensions as histogram (see examples)
Mask_2="\PATH\TO\MASK\MASK2_FILE_NAME" # Path to masks which should be applied - same dimensions as histogram (see examples)
Lets assume that the Hist_E_S_2
is following histogram: NBinX = 5, Xmin = 0, Xmax = 10 & NBinY = 4, Xmin = 0, Xmax = 4 then the mask should have following form:
0 1 0 0 0 0 0 0 0 0 0 2 0 0 3 0 0 0 4 10
A single gap is used as number separator.This will use bins [X,Y,Weight] = [2,1,1], [2,3,2], [5,2,3], [4,4,4], [4,5,10] and the matrix is read from first line to the last one (1st line is Y = 1).
Normalization
Both 1D and 2D histograms can be normalized before they are used calculation of the SigVec:
Normalize=1 # 1 to DO normalization - 0 to NOT do norm. - default is 0
Default Configuration
[Var_S]
Title="Var_S"
Intervals=-1,5,5,20,20,1e+200
Normalize=1
[Var_BorderPixN]
Title="Var_BorderPixN"
Intervals=-1,2,2,5,5,15,15,1e+200
Normalize=1
[Var_Lin]
Title="Var_Lin"
Intervals=-1,0.5,0.5,0.95,0.95,1e+200
Normalize=1
Minimum values -1 and maximum 1e+200 are safety precautions to be sure that all relevant bins are taken into account.
Log Information
SigVec conf file - name: SigVec.ini SigVec conf file - path: ./Test/data/test_008/ ... -------------------------------------------------- SIGNIFICANT VECTOR Features: Title = Epix HistName = Hist1D_EpixMean_Total Intervals Normalize= True Title = H HistName = Hist1D_H_Total Intervals Normalize= True Title = LET HistName = Hist1D_LET_Total Intervals Normalize= True Title = E HistName = Hist1D_E_Total Intervals Normalize= True Title = S HistName = Hist1D_S_Total Intervals Normalize= True -------------------------------------------------- ...
The first part informing about name and path to the config file of significant vector can be seen in the first stage of DPE processing (if it is not shown then default option is used).
The second part is in the same processing/nationalization stage of DPE and it shows which histograms and intervals/masks will be used to creation of the vector.
... ======================================================================= SIGNIFICANT VECTOR PROCESSING ======================================================================= * Processing histogram: Hist1D_EpixMean_Total Histogram intervaling: [0,50][50,1e+200] Intervaling results: 0.967578 0.032422 * Processing histogram: Hist1D_H_Total Histogram intervaling: [0,50][50,200][200,1e+200] Intervaling results: 0.596236 0.402650 0.001114 * Processing histogram: Hist1D_LET_Total Histogram intervaling: [0,1][1,3][3,1e+200] Intervaling results: 0.915691 0.084309 0.000000 * Processing histogram: Hist1D_E_Total Histogram intervaling: [0,200][200,1000][1000,1e+200] Intervaling results: 0.686559 0.313417 0.000024 * Processing histogram: Hist1D_S_Total Histogram intervaling: [0,5][5,20][20,1e+200] Intervaling results: 0.595407 0.387884 0.016709
This last parts includes information about created elements of the vector, from which histograms were created and what are the results.
Particle Identification
The DPE engine is also capable of basic particle identification/classification.
This output can be exported into extended elist where new column/PIDClass will be created with information about recognized class.
The choice of the PID algorithm is based on the switch PIDAlg
whose numerical value are listed below. The settings can be done in the main configuration file:
PIDAlg = 201
the current possible values are following:
- 101 - Heuristic simplified DT for TPX 300 um Si
- 102 - Heuristic decision tree with 8 classes for TPX 300 um Si
- 103 - Heuristic decision tree with 16 classes for TPX3 500 um Si
- 201 - Dense neural network with 3 classes for TPX 300 um Si
- 202 - Dense neural network with 6 classes for TPX 300 um Si
- 251 - Dense neural network with 3 classes for TPX3 500 um Si
- 252 - Dense neural network with 6 classes for TPX3 500 um Si
If not value is given to the program then DT is chosen based on the detector configuration.
Algorithms based on Neural Networks
Dense neural networks for TPX 300 um Si
- Dense neural network TPX with 3 classes
- Protons
- Photons && Electrons
- Ions
- Others
- Dense neural network TPX with 6 classes
- Protons LE (<30 MeV)
- Protons ME (>30 & <100 MeV)
- Protons HE (>100 MeV)
- Photons && Electrons
- Helium ions
- Ions (except He)
- Others
Dense neural networks for TPX3 500 um Si
- Dense neural network TPX3 with 3 classes
- Protons
- Photons && Electrons
- Ions
- Others
- Dense neural network TPX3 with 6 classes
- Protons LE (<30 MeV)
- Protons ME (>30 & <100 MeV)
- Protons HE (>100 MeV)
- Photons && Electrons
- Helium ions
- Ions (except He)
- Others
Algorithms based on Decision Trees
Decision tree for TPX 300 um Si
- Heuristic simplified DT with 3 classes:
- Low LET: electrons, X-rays, gamma
- Mid LET: protons
- High LET: ions
- Others
- Heuristic decision tree with 8 classes:
- X-rays; LE electrons OD; HE electrons OD; muons PP
- LE protons OD; HE protons PP
- LE alphas OD; HE alphas PP
- LE ions OD; HE ions PP
- HE electrons OD; muons nPP
- HE protons nPP; UHE protons nPP; UHE alphas nPP
- He alphas nPP; UHE light ions nPP
- He ions nPP
- Others
Decision tree for TPX3 500 um Si
- Heuristic decision tree with 16 classes:
- X-rays; HE electrons PP; HE protons PP
- Gamma rays LE (e.g. 137 Cs)
- Gamma rays HE (e.g. 60 Co)
- LE electrons OD (< 10 MeV)
- HE electrons nPP (> 10 MeV)
- LE and HE electrons PP
- LE protons (< 3MeV)
- ME protons (3-10 MeV)
- HE protons (>10 MeV)
- LE alphas (<10 MeV)
- HE alphas (>10 MeV)
- LE ions (<10 MeV/u)
- HE ions (>10 MeV/u)
- LE, thermal and slow neutrons ( < 0.5 eV)
- HE fast neutrons ( > 1 MeV)
- Others
Explanation of used abbreviations: LE - Low Energy, HE - High Energy, UHE - Ultra High Energy ,OD - Omni Directional, PP - Perpendicular to the sensor, nPP - non Perpendicular. Their values differentiate for individual decision tree.
Log and Statistical Information
-------------------------------------------------- PARTICLE CLASSIFICATION Algorithm was chosen based on detector settings. Alg description: Dense neural network Alg switch: 251 Optimized for detector Chip type: TPX3 Sensor material: Si Sensor thickness: 500 um Sensor bias: 80 V Classes: 1) Protons 2) Photons && Electrons 3) Ions 4) Others -------------------------------------------------- ... -------------------------------------------------- PARTICLE CLASSIFICATION -------------------------------------------------- ClassNames: 1) Protons 2) Photons && Electrons 3) Ions 4) Others AlgorithmSwitch [-]: 251 CountParticle_Sum_Class [-]: 701, 41492, 0, 0 CountParticle_Sum_Class [%]: 1.66, 98.34, 0, 0 CountRate_Mean_Class [s-1]: 70, 4150, 0, 0 Fluence_Sum_Class [cm-2]: 353.6004, 20929.5099, 0, 0 Flux_Sum_Class [cm-2 s-1]: 35.3660, 2093.3022, 0, 0 EnergyDep_Sum_Class [keV+1]: 434964.7790, 6105693.4420, 0, 0 Dose_Sum_Class [uGy+1]: 0.30187, 4.2374, 0, 0 DoseRate_Mean_Class [uGy+1 h-1]:108.6909, 1525.7176, 0, 0 TotalValues_Class: ========================================================================================================== | | N_Class | Fluence_Class| Flux_Class | Edep_Class | Dose_Class | DR_Class | | ClassNum |-------------- -------------- -------------- -------------- -------------- --------------| | 1 | 701 | 70.1118 | 35.3660 | 434964.7790 | 0.30187 | 108.6909 | | 2 | 41492 | 4149.8963 | 2093.3022 | 6105693.4420 | 4.2374 | 1525.7176 | | 3 | 0 | 0 | 0 | 0 | 0 | 0 | | Others | 0 | 0 | 0 | 0 | 0 | 0 | ==========================================================================================================
Information about used algorithm, for which detector the given algorithm was optimized and into which classes the cluster will be separated. The second part includes the same information as those included in the sampling list, see more information in the text below and the example in json format.
Statistical information is given in the general sampling list in the directory File
.
"ClassNames":[ "Protons", "Photons && Electrons", "Ions", "Others"], // Name of particle classes.
"AlgorithmSwitch_cnt":251, // Switch of algorithm used for classification
"CountParticle_Sum_Class_cnt":[ 701, 41492, 0, 0 ], // Particle sum for given classes (numbers follows the list of classes).
"CountParticle_Sum_Class_perc":[ 1.66, 98.34, 0, 0 ], // Particle sum in percentages for given classes (numbers follows the list of classes).
"CountRate_Mean_Class_s-1":[ 70, 4150, 0, 0 ], // Sum of mean values of count rates for given classes (numbers follows the list of classes).
"Fluence_Sum_Class_cm-2":[ 353.6004, 20929.5099, 0, 0 ], // Sum of fluences for given classes (numbers follows the list of classes).
"Flux_Sum_Class_cm-2s-1":[ 35.3660, 2093.3022, 0, 0 ], // Sum of mean values of particle flux for given classes (numbers follows the list of classes).
"EnergyDep_Sum_Class_keV+1":[ 434964.7790, 6105693.4420, 0, 0 ], // Sum of deposited energies for given classes (numbers follows the list of classes).
"Dose_Sum_Class_uGy+1":[ 0.30187, 4.2374, 0, 0 ], // Sum of doses for given classes (numbers follows the list of classes).
"DoseRate_Mean_Class_uGy+1h-1":[ 108.6909, 1525.7176, 0, 0 ], // Sum of mean values oi dose rates for given classes (numbers follows the list of classes).
...
"CountParticle_Sample_Class1_cnt":[ 2, 2, 5, 4, 8, 3, 1, 3, 2, 1, 4, 2, 3, 3, 3, 3, 3, 1, 2, 1, 5, 3, 1, 3, 3, 5, 2, 3, 5, 4, 2, 3, 2, 8, 3, 0, 6, 2, 3, 2, 1, 3, 1, 2, 7, 5, 3, 1, 5, 5, 2, 4, 2, 4, 5, 1, 3, 7, 3, 4, 1, 4, 7, 2, 4, 2, 5, 3, 3, 3, 2, 6, 6, 1, 5, 4, 2, 4, 1, 6, 5, 0, 4, 4, 7, 2, 3, 3, 6, 1, 2, 5, 5, 4, 4, 3, 5, 0, 1, 5, 4, 5, 4, 3, 4, 5, 3, 4, 2, 4, 5, 3, 1, 4, 3, 2, 9, 2, 4, 6, 6, 6, 4, 3, 1, 5, 4, 3, 4, 2, 5, 3, 6, 2, 3, 5, 3, 5, 2, 2, 5, 3, 4, 0, 1, 3, 4, 2, 6, 1, 5, 5, 6, 4, 3, 2, 4, 2, 6, 7, 3, 2, 4, 6, 2, 3, 4, 0, 2, 3, 4, 3, 4, 4, 5, 4, 2, 1, 8, 7, 2, 5, 6, 3, 7, 3, 5, 1, 2, 4, 4, 5, 4, 4, 8, 3, 4, 4, 2, 1 ],
"CountParticle_Sample_Class2_cnt":[ 207, 206, 192, 205, 222, 223, 211, 225, 205, 201, 184, 231, 209, 189, 195, 225, 192, 216, 206, 184, 212, 211, 209, 188, 212, 212, 247, 203, 185, 208, 215, 224, 207, 206, 202, 195, 196, 185, 186, 181, 228, 209, 203, 193, 189, 164, 196, 224, 208, 205, 196, 188, 210, 207, 223, 219, 207, 232, 195, 207, 184, 198, 214, 232, 194, 215, 210, 214, 221, 212, 227, 221, 220, 222, 203, 202, 238, 208, 209, 222, 211, 228, 203, 206, 229, 195, 217, 202, 189, 206, 214, 227, 202, 215, 249, 223, 218, 240, 214, 220, 218, 225, 197, 178, 194, 212, 184, 204, 214, 196, 177, 210, 210, 193, 188, 193, 201, 149, 208, 183, 193, 212, 203, 180, 205, 227, 209, 170, 215, 193, 210, 196, 191, 220, 185, 199, 192, 175, 235, 228, 214, 205, 231, 191, 198, 192, 214, 220, 219, 257, 218, 191, 236, 213, 234, 207, 214, 240, 240, 223, 209, 232, 211, 217, 235, 209, 192, 191, 199, 225, 214, 169, 226, 234, 202, 197, 205, 205, 184, 208, 218, 198, 217, 190, 207, 221, 217, 172, 248, 207, 176, 213, 201, 212, 189, 202, 207, 227, 173, 215 ],
"CountParticle_Sample_Class3_cnt":[ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ],
"CountParticle_Sample_ClassOthers_cnt":[ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 ],
...
The second part of the output into sampling list contains information about individual physical products and its dependence on sampling time and used particles class. Based on the example above, one array includes values of count rate for given class in given time (sampled based on sampling time).
Radiation Field Recognition
The DPE is capable of basic radiation filed recognition (RFR). The current version supports following algorithm:
- Distance Comparator
To use the RFR, it is needed to use standard/default settings of the DPE with respect to the histograms and significant vectors.
Distance Comparator
The distance comparator uses the significant vectors of premesured known radiation fields (database of vectors) and compare them with the given data for processing.
There has to be match between the database detector configuration and the data configuration.
The comparison is based on distance evaluation between the known and unknown vector. The final probability that given data is one of the source is based on in-proportional relation between probability and the distance.
You can find below currently included databases for the detector configuration.
CdTe 2mm TPX3 and efficiency results:
IN/OUT Ba133 Cs137 Eu152 Co60 Am241 Na22 Ba133 97.40 0.40 1.31 0.26 0.33 0.30 Cs137 3.22 77.73 4.51 4.56 2.40 7.57 Eu152 2.70 1.08 94.16 0.62 0.69 0.74 Co60 1.06 2.53 1.29 89.48 1.22 4.41 Am241 1.47 1.37 1.56 1.31 93.04 1.24 Na22 0.87 3.13 1.10 3.28 0.81 90.81
…
Post-processing and Physical Products
The main physical products which can be found in the exported files are following:
- Flux
- Count of particles and pixels
- Dose rate
- Deposited energy
- Distributions of cluster variables
- Spatial event maps of cluster variables
- Spatial map of clusters and individual clusters
The flux and dose rate account for possible masking. In the output sampling list, all time coupled variables are in the cases of classes calculated with respect to the sampling time. The total variables are calculated with respect to the elapsed time which is usually shorter. The elapsed time is defined as time of the last processed event/particle. Therefore these values are bigger than the sum/mean of variables of classes.
Time Sampling of Data
Most of the physical products are calculated based on a sampling time. This means that all events/particles which are within a time interval of sampling time from some starting point contribute to the physical products.
Example, lets assume that measurement of 10 s were done and 20000 particles were registered in first 5 seconds and 30000 particles in another 5 seconds. If the sampling time is chosen as 5 s than 2 samples are created in the output file with two values for each sampled physical product. The values of flux, if no mask is used, are then: 2000 and 3000 particles/s cm-2. The total flux is calculated based on the elapsed time and if the last event was detected at time of 9 s then the total flux is: (20000 + 30000 particles)/(2 cm2 * 9 s) = 2777,8 particles/s cm-2 which is more than the mean value of the class fluxes 2500 particles/s cm-2.
In the case of the frame, all the physical products are calculated and visualized with respect tot he detector's live time. In other words, the dead time is not accounted for (it is estimation is nevertheless done and can be found in the frame list).
Sampling List
The sampling list includes general information about all physical products and their dependence on sampling time and particle class. It also includes general information about PID.
There two formats for export: plain text (copy of terminal print) and json. In the json file, the key is created from part with name and Explanation of individual parameters is given in the example below after //
:
"TimeLive_Sum_s+1":9.998322064, // Total measured time in s
"TimeSampling_s+1":0.05000000000, // Sampling time in s
"CountSample_cnt":200, // Count of samples
"Time_First_ns+1":0, // Time used as start reference in ns
"Time_Last_ns+1":9998322064.06, // Time of last event in ns
"Time_Last_s+1":9.998322064, // Time of last event in s
"CountPixHit_Sum_cnt":244333, // Total sum of all hit pixels
"CountParticle_Sum_cnt":42193, // Total sum of all particles/clusters
"CountRate_Mean_s-1":4220.0081, // Total particle rate in s-1
"CountRatePixHit_Mean_s-1":24437.4004, // Total pixel rate in s-1 (can be used to estimate readout overwhelming)
"Fluence_Sum_cm-2":21283.1103, // Total fluence of particle in cm-2
"Flux_Sum_cm-2s-1":2128.6682, // Total flux of particles in cm-2 s-1
"EnergyDep_Sum_keV":6540658.2210, // Total deposited energy in keV
"Dose_Sum_uGy+1":4.5393, // Total dose of particles in uGy or Gy (DoDoseUnitRadiology)
"DoseRate_Mean_uGy+1h-1":1634.4086, // Total dose rate of particles in uGy h-1 or Gy s-1 (DoDoseUnitRadiology)
// Described in the section about PID
"ClassNames":[ "Protons", "Photons && Electrons", "Ions", "Others"],
"AlgorithmSwitch_cnt":251,
"CountParticle_Sum_Class_cnt":[ 701, 41492, 0, 0 ],
"CountParticle_Sum_Class_perc":[ 1.66, 98.34, 0, 0 ],
"CountRate_Mean_Class_s-1":[ 70, 4150, 0, 0 ],
"Fluence_Sum_Class_cm-2":[ 353.6004, 20929.5099, 0, 0 ],
"Flux_Sum_Class_cm-2s-1":[ 35.3660, 2093.3022, 0, 0 ],
"EnergyDep_Sum_Class_keV+1":[ 434964.7790, 6105693.4420, 0, 0 ],
"Dose_Sum_Class_uGy+1":[ 0.30187, 4.2374, 0, 0 ],
"DoseRate_Mean_Class_uGy+1h-1":[ 108.6909, 1525.7176, 0, 0 ],
// Integrated sampling time for individual samples (last reflects the last time so it does not have to be whole multiple of sampling time)
"TimeLive_Int_Sample_s+1":[ 0.050000, 0.10000, 0.15000, 0.20000],
// Sampling time in each sample
"TimeSampling_Int_Sample_s+1":[ 0.050000, 0.050000, 0.050000, 0.050000],
// All physical products as in the tables above but for individual time samples
"CountPixHit_Sum_Sample_cnt":[ 1239, 1225, 1204, 11406 ],
"CountParticle_Sum_Sample_cnt":[ 209, 208, 197, 209],
"CountRate_Mean_Sample_s-1":[ 4180, 4160, 3940, 4180],
"Fluence_Sum_Sample_cm-2":[ 105.4244, 104.9199, 99.3713, 105.4244 ],
"Flux_Sum_Sample_cm-2s-1":[ 2108.4872, 2098.3988, 1987.4257, 2108.4872],
"EnergyDep_Sum_Sample_keV+1":[ 34408.3804, 31737.0849, 32648.1435, 30391.8794 ],
"Dose_Sum_Sample_uGy+1":[ 85.9668, 79.2927, 81.5689, 75.9318 ],
"DoseRate_Mean_Sample_uGy+1h-1":[ 1719.3351, 1585.8544, 1631.3787, 1518.6365 ],
// described in the section about PID
"CountParticle_Sample_Class1_cnt":[ 2, 2, 5, 4],
"CountParticle_Sample_Class2_cnt":[ 207, 206, 192, 205],
"CountParticle_Sample_Class3_cnt":[ 0, 0, 0, 0 ],
"CountParticle_Sample_ClassOthers_cnt":[ 0, 0, 0, 0],
// same as above but for count rate
"CountRate_Sample_Class1_s-1":[ 40, 40, 100, 80],
"CountRate_Sample_Class2_s-1":[ 4140, 4120, 3840, 4100 ],
"CountRate_Sample_Class3_s-1":[ 0, 0, 0, 0 ],
"CountRate_Sample_ClassOthers_s-1":[ 0, 0, 0, 0],
// same as above but for count rate
"Fluence_Sample_Class1_cm-2":[ 1.0088, 1.0088, 2.5221, 2.0177],
"Fluence_Sample_Class2_cm-2":[ 104.4155, 103.9111, 96.8492, 103.4067 ],
"Fluence_Sample_Class3_cm-2":[ 0, 0, 0, 0 ],
"Fluence_Sample_ClassOthers_cm-2":[ 0, 0, 0, 0],
// same as above but for count rate
"Flux_Sample_Class1_cm-2s-1":[ 20.1769, 20.1769, 50.4423, 40.3538 ],
"Flux_Sample_Class2_cm-2s-1":[ 2088.3103, 2078.2218, 1936.9835, 2068.1334],
"Flux_Sample_Class3_cm-2s-1":[ 0, 0, 0, 0],
"Flux_Sample_ClassOthers_cm-2s-1":[ 0, 0, 0, 0 ],
// same as above but for count rate
"EnergyDep_Sample_Class1_keV+1":[ 1377.3510, 1318.7130, 3022.6320, 2560.5120 ],
"EnergyDep_Sample_Class2_keV+1":[ 33031.0294, 30418.3719, 29625.5115, 27831.3674],
"EnergyDep_Sample_Class3_keV+1":[ 0, 0, 0, 0 ],
"EnergyDep_Sample_ClassOthers_keV+1":[ 0, 0, 0, 0],
// same as above but for count rate
"Dose_Sample_Class1_uGy+1":[ 0.00095589, 0.00091520, 0.0020977, 0.0017770 ],
"Dose_Sample_Class2_uGy+1":[ 0.022924, 0.021111, 0.020560, 0.019315],
"Dose_Sample_Class3_uGy+1":[ 0, 0, 0, 0 ],
"Dose_Sample_ClassOthers_uGy+1":[ 0, 0, 0, 0],
// same as above but for count rate
"DoseRate_Sample_Class1_uGy+1h-1":[ 68.8242, 65.8941, 151.0364, 127.9449 ],
"DoseRate_Sample_Class2_uGy+1h-1":[ 1650.5109, 1519.9603, 1480.3423, 1390.6916 ],
"DoseRate_Sample_Class3_uGy+1h-1":[ 0, 0, 0, 0],
"DoseRate_Sample_ClassOthers_uGy+1h-1":[ 0, 0, 0, 0, 0 ]
Frame Analysis
…
Notes
- acq time can be estimated only for clog input files (for matrix files it is unknown)
Information
"CountFrame_cnt":7, // Count of frames
"CountFrame_NonEmpty_cnt":4, // Count of non empty frames
"CountFrame_Empty_cnt":3, // Count of empty frames
"TimeAcq_s+1":0.005000000000, // Estimation of acquisition time of frames in s (only for clog data)
"TimeDead_Mean_s+1":1.495000000, // Mean value of dead time estimation in s -> time needed for readout of frame
"TimeDead_Std_s+1":0.5477221222, // Std of above in s
"FramesPerSecond_Mean_s-1":0.6667, // Frame rate in s-1 -> 1/(TimeAcq_s+1 + TimeDead_Mean_s+1)
"EnergyDep_MeanStdMinMax_keV+1":[ 25.5569, 32.7182, 75.3443, 0 ], // Mean, std, minimum and maximum of deposited energy in one frame in keV
"Dose_MeanStdMinMax_Gy+1":[ 0.000000000017737, 0.000000000022707, 0.000000000052289, 0 ], // Mean, std, minimum and maximum of dose in one frame Gy or uGy (DoDoseUnitRadiology)
"DoseRate_MeanStdMinMax_Gy+1s-1":[ 0.0000000035473, 0.0000000045413, 0.000000010458, 0 ], // Mean, std, minimum and maximum of dose rate in one frame in Gy s-1 or uGy h-1 (DoDoseUnitRadiology)
"EnergyDepPix_MeanStdMinMax_keV+1":[ 6.2823, 6.2487, 15.0689, 0 ], // Mean, std, minimum and maximum of per pixel deposited energy in one frame in keV
"CountParticles_MeanStdMinMax_cnt":[ 1.4286, 1.3973, 3, 0 ], // Mean, std, minimum and maximum of particle count in one frame
"CountRate_MeanStdMinMax_s-1":[ 285.7143, 279.4553, 600, 0 ], // Mean, std, minimum and maximum of particle count rate in one frame
"Fluence_MeanStdMinMax_cm-2":[ 0.72060, 0.70482, 1.5133, 0 ], // Mean, std, minimum and maximum of flunce in one frame in cm-2
"Flux_MeanStdMinMax_cm-2s-1":[ 144.1208, 140.9636, 302.6537, 0 ], // Mean, std, minimum and maximum of flux in one frame in cm-2 s-1
"CountPixHit_MeanStdMinMax_cnt":[ 2.1429, 2.4785, 6, 0 ], // Mean, std, minimum and maximum of count of hit pixels in one frame
"OccupancyPix_MeanStdMinMax_perc":[ 0.0032697, 0.0037819, 0.0091553, 0 ], // Mean, std, minimum and maximum of occupancy in one frame in percents
"Time_Frame_s+1":[ 0.0050000, 1.5050, 3.0050, 4.5050, 6.0050, 7.5050, 9.0050 ], // Time of individual frames
"CountHitPix_Frame_cnt":[ 6, 5, 0, 2, 2, 0, 0 ], // Count of hit pixels in individual frames
"CountParticle_Frame_cnt":[ 3, 3, 0, 2, 2, 0, 0 ], // Count of particles in individual frames
"CountRate_Frame_s-1":[ 600, 600, 0, 400, 400, 0, 0 ], // Count rate of particles in individual frames
"Fluence_Frame_cm-2":[ 1.5133, 1.5133, 0, 1.0088, 1.0088, 0, 0 ], // Fluence in individual frames
"Flux_Frame_cm-2s-1":[ 302.6537, 302.6537, 0, 201.7691, 201.7691, 0, 0 ], // Flux in individual frames
"OccupancyPix_Frame_perc":[ 0.0091553, 0.0076294, 0, 0.0030518, 0.0030518, 0, 0 ], // Occupancy in individual frames
"EnergyDep_Frame_keV+1":[ 68.6103, 75.3443, 0, 17.0665, 17.8773, 0, 0 ], // Deposited energy in individual frames
"EnergyDepPix_Frame_keV+1":[ 11.4350, 15.0689, 0, 8.5333, 8.9387, 0, 0 ], // Per pixel energy in individual frames
"Dose_Frame_Gy+1":[ 0.000000000047616, 0.000000000052289, 0, 0.000000000011844, 0.000000000012407, 0, 0 ], // Dose in individual frames
"DoseRate_Frame_Gy+1s-1":[ 0.0000000095232, 0.000000010458, 0, 0.0000000023689, 0.0000000024814, 0, 0 ] // Dose rate in individual frames
Graphical Output
Compton Directional Analysis
Notes
DoCompton
- can be used only with configuration file
Configuration File
[Settings]
EnergyComptonSumMin=350 # minimum sum of energy for compton event
EnergyComptonSumMax=380 # maximum sum of energy for compton event
DoSkipSelfDetectorCoinc=0 # 1 = true, 0 = false: if true, the program will skip events where more clusters are from one detector
DoOnlyPairCoinc=0 # 1 = true, 0 = false: if true, the program only process events of two clusters (more than 2 are ignored)
TimeChargeCollCdTe=86 # time for collection of charge from CdTe of given thickness and bias in the assambly settings (86 ns for 2 mm -500 V)
ProjDistance= 35 #
ProjNBinX= 100 #
ProjNBinY= 100 #
ProjBinWidthX = 1 #
ProjBinWidthY = 1 #
ForceBackward = 1 #
ForceForward = 0 #
[Assembly0]
NumID=0 #
Chip_Type="TPX3" #
Sensor_Material="CdTe" #
Sensor_Thickness=2000 #
PositionA=-7.0125, -7.0125, 0 #
PositionB=7.0675, 7.0675, 2.0 #
Vector=1, 0, 0 #
[DetectorA]
Key="J08" #
Readout="Minipix" #
Assemblies=0 #
Information
"CountPossCompEvents_cnt":12, // Count of possible Compton events of clusters passing basic energy condition (several combination in one coincidence group is taken into account creating one Compton group of combinations)
"CountAmbigCompGroup_cnt":0, // Count of possible Compton groups which can NOT be simplified to one unique Compton group
"CountUniqCompGroup_cnt":12, // Count of unique Compton groups which could be created from several combinations of clusters from one coincidence event, but only this one satisfied all criteria.
"CountCompGroup_cnt":12, // Total count of all Compton groups.
"CountDirForwCompGroup_cnt":12, // Count of Compton groups which are only forward
"CountDirBackCompGroup_cnt":0, // Count of Compton groups which are only backward
"CountDirUniqueCompGroup_cnt":10, // Count of Compton groups which are unique from directional point of view (only forward or backward)
"CountDirAmbigCompGroup_cnt":0, // Count of Compton groups which can be both, forward or backward
"CountSingleDetCompGroup_cnt":12, // Count of Compton groups from one detector
"CountMultiDetCompGroup_cnt":0, // Count of Compton groups from several detectors
"CountPassInProcCompGroup_cnt":12, // Count of Compton groups passing into projection
"CountFailInProcCompGroup_cnt":0, // Count of Compton groups failing before projection
"CountInProjCompGroup_cnt":12, // Count of Compton groups which are IN projection plane
"CountOutProjCompGroup_cnt":0, // Count of Compton groups which are OUT projection plane
"ConeHalfAngle_MeanStdMinMax_rad+1":[ 0.92899, 0.64062, 0.14930, 1.8912 ], // Compton cone half angle mean, std, maximum and minimum value in radians
"TimeLive_Int_Sample_s+1":[ 0.5, 0.98812 ], // Integrated sampling time for individual samples (last reflects the last time so it does not have to be whole multiple of sampling time)
"ConeHalfAngle_Mean_Sample_rad":[ 0.875, 0.92899 ], // Mean value of Compton cone half angle in each time sample
"CountCompGroup_Sum_Sample_cnt":[ 4, 8 ] // Count of all Compton groups in individual time samples
Graphical Output
To simplify the created projection alias estimation of the radiation source position, an additional image post-processing was introduced:
- Contrast – removes all pixels whose values is less than fraction of the maximal value in pixels. Default is set to 70 percent alias 0.7 in the fraction terms.
- Blurring – exploits Gaussian filter to blur the image and highlight the estimation.
Directional Analysis
DPE is capable to estimate several directional features of detected particles.
It is possible to use several different algorithms for correcting the 2D length of clusters for charge sharing effect. The control LengthCorr
in main configuration file can be set to following values: 1) weighted standard deviation subtraction from projected length. 2) projected width subtraction from projected length. 3) model of Nabha.
The 3D length is always calculated based on the 2D corrected length and sensor thickness assuming that the particle fully crossed the sensor.
DoDirection
- Control whether Compton directional reconstruction should be done. There is limited list of cases when this analysis is usable based on the detector settings (see section below). Example: DoDirection = true
, DoDirection = false
. Default: true. DoDirectionTrackCond
- Control to use tracking condition during directional analysis (additional constrain on features of clusters). Example: DoDirectionTrackCond = true
, DoDirectionTrackCond = false
. Default: false.
Information
"CountParticle_Sum_cnt":208, // Total count of analyzed particles
"Length2DGeoProj_MeanStdMinMax_px+1":[ 30.3412, 0.81544, 26.2376, 33.9519 ], // Mean, std, maximum and minimum value of 2D length geometrical projection of clusters in px
"Length2DCorrected_MeanStdMinMax_px+1":[ 29.2538, 0.81617, 25.3078, 32.6405 ], // Mean, std, maximum and minimum value of corrected 2D length in px
"Length3D_MeanStdMinMax_um+1":[ 1684.9090, 42.8482, 1479.0100, 1863.5600 ], // Mean, std, maximum and minimum value of 3D length (based on corrected 2D and sensor thickness) in px
"AzimuthAngle_MeanStdMinMax_deg+1":[ 15.8250, 0.57266, 12.7564, 17.7723 ], // Mean, std, maximum and minimum value of 2D length geometrical projection of clusters in px
"ElevationAngle_MeanStdMinMax_deg+1":[ 72.7252, 0.45788, 70.2410, 74.4366 ], // Mean, std, maximum and minimum value of 2D length geometrical projection of clusters in px
// The same as above but for individual time samples (first line is integrated time alive)
"TimeLive_Int_Sample_s+1":[ 10, 20, 30, 40, 50, 60, 70, 80, 90, 100, 110, 120, 130, 140, 150, 160, 170, 180, 190, 200, 210, 220, 230, 240, 250, 260, 270, 280, 290, 295.2798 ],
"Length2DGeoProj_Mean_Sample_px+1":[ 30.6116, 0, 30.2988, 0, 0, 30.1250, 0, 0, 30.1191, 0, 30.2373, 0, 0, 30.6139, 0, 0, 30.4191, 0, 30.0041, 0, 0, 30.2577, 0, 0, 30.3248, 0, 30.5943, 0, 0, 30.5333 ],
"Length2DCorrected_Mean_Sample_px+1":[ 29.4190, 0, 29.3034, 0, 0, 29.0758, 0, 0, 29.0180, 0, 29.1283, 0, 0, 29.6295, 0, 0, 29.4238, 0, 28.7745, 0, 0, 29.2299, 0, 0, 29.2946, 0, 29.5171, 0, 0, 29.3397 ],
"Length3D_Mean_Sample_um+1":[ 1693.6027, 0, 1687.6217, 0, 0, 1675.5644, 0, 0, 1672.5227, 0, 1678.3014, 0, 0, 1704.6150, 0, 0, 1693.8128, 0, 1659.7406, 0, 0, 1683.6386, 0, 0, 1687.0558, 0, 1698.7600, 0, 0, 1689.4247 ],
"AzimuthAngle_Mean_Sample_deg+1":[ 15.9028, 0, 15.8854, 0, 0, 15.9136, 0, 0, 15.7613, 0, 15.6879, 0, 0, 15.9037, 0, 0, 15.8522, 0, 15.8963, 0, 0, 15.8301, 0, 0, 16.0281, 0, 15.8330, 0, 0, 15.6919 ],
"ElevationAngle_Mean_Sample_deg+1":[ 72.8143, 0, 72.7298, 0, 0, 72.6255, 0, 0, 72.5953, 0, 72.6597, 0, 0, 72.9396, 0, 0, 72.8254, 0, 72.4583, 0, 0, 72.7155, 0, 0, 72.7478, 0, 72.8675, 0, 0, 72.7720 ]
Graphical Output
Coincidence and Event Analysis
- TimeCoinc - [FLOAT] Time in nanoseconds for evaluations coincidence events -> clusters whose times are in between this this interval are grouped as coincidence group (first cluster time is down edge and plus coince group time is upper edge). Example:
TimeCoinc = 1e2
,TimeCoinc = 10
,TimeCoinc = 1.2554
. Default: 100.
Information
"TimeCoincWindow_ns+1":32.8120, // Used coincidence window
"CountEvent_Sum_cnt":9, // Total count of all events (not particles -> event is one particle or group of coincidence particles)
"CountCoincGroup_Sum_cnt":2, //
"CountCoincGroup_Sum_perc":22.2222, //
"CountParticle_Mean_Event_cnt":1.4444, //
"CountParticle_Sum_NoCoinc_cnt":7, //
"CountParticle_Sum_NoCoinc_perc":53.8462, //
"CountParticle_Sum_CoincGroup_cnt":6, //
"CountParticle_Sum_CoincGroup_perc":46.1538, //
"CountParticle_Mean_CoincGroup_cnt":3, //
"CountParticle_Max_CoincGroup_cnt":4, //
"TimeLive_Int_Sample_s+1":[ 1, 2, 3, 3.0000 ], //
"CountParticle_Mean_Sample_Event_cnt":[ 1, 0, 1, 2 ],//
"CountCoincGroup_Sample_cnt":[ 0, 0, 0, 2 ] //
Graphical Output
Spatial Maps
….
Event Visualization
Another graphical output produced in the engine are visualisations of particle tracks in sensor plane. These plots can be only obtained for raw data containing pixelated information. Example of this output can be seen in the where deposited energy in each pixel is used for this purpose.
The individual events visualisations are completed with so called integrated visualisations. In these plots several particles are shown together to obtain more statistical based information about detected tracks. Two versions of these plots are available: integration of N particles (N is usually around 100) and integration of all detected particles.
Event visualisation in graphical form is completed with text alternative in matrix form.
Error Codes
The engine run can produce an error code into standard error stream (stderr
) and into the standard output stream (stdout
). Successful engine run produces/returns 0.
Any other value signals error in the run. Possible values are listed below (negative numbers with explanation after ///<
):
#define DPE_ERR_NO_ERROR 0 ///< No error occurred.
#define DPE_ERR_NO_PAR_FILE -1000 ///< Missing file with parameters.
#define DPE_ERR_INIT_GEN -1001 ///< Initialization has not been done.
#define DPE_ERR_READ_PAR_GEN -1002 ///< UNUSED
#define DPE_ERR_OPEN_LOG -1003 ///< Can not open log file.
#define DPE_ERR_NO_IN_FILE -1004 ///< Can not find any files for processing.
#define DPE_ERR_NO_IN_FILE_SEL_IGN -1005 ///< No files have passed the selection and ignore criteria.
#define DPE_ERR_NO_INIT_PID -1008 ///< PID init has not been done.
#define DPE_ERR_FIND_CLUSTERER -1009 ///< Can not find clusterer binary.
#define DPE_ERR_FIND_MODELS -1010 ///< Can not find models directory.
#define DPE_ERR_CHECK_PYTHON -1011 ///< Missing python for graphics creation.
#define DPE_ERR_CHECK_PYTHON_MOD -1012 ///< Missing python modules for graphics creation (names are in stderr and stdout).
#define DPE_ERR_IN_FILE_OPEN -1013 ///< Can not open input file.
#define DPE_ERR_IN_SIMPLE_ELIST -1014 ///< Input data recognized as ElistSimple. Can not be processed because the names of variables are missing.
#define DPE_ERR_IN_FILE_CANOP -1015 ///< Can not open input file.
#define DPE_ERR_IN_FILE_ELIST_BAD -1016 ///< Corrupted elist file.
#define DPE_ERR_IN_FILE_UNKNOWN -1017 ///< Unknown in file format.
#define DPE_ERR_IN_FILE_CLOG_BAD -1018 ///< Corrupted clog file.
#define DPE_ERR_IN_FILE_CLOG_EMPTY -1019 ///< Empty clog file.
#define DPE_ERR_IN_FILE_EMPTY -1020 ///< Empty in file format.
#define DPE_ERR_IN_FILE_MATRIX_TOA -1021 ///< Unsupported frame format of ToA.
#define DPE_ERR_FILTER_NO_CONF -1022 ///< Can not open file with filter config.
#define DPE_ERR_INIT_MODULES -1023 ///< Error occurred during module initialization.
#define DPE_ERR_HIST_NO_CONF -1024 ///< Can not open file with histograms config.
#define DPE_ERR_HIST_INIT -1025 ///< Error occurred during histogram module initialization.
#define DPE_ERR_SIGVEC_NO_CONF -1026 ///< Can not open file with sig vec config.
#define DPE_ERR_COMPAR_INIT -1027 ///< Error occurred during comparator module initialization.
#define DPE_ERR_MASK_NO_CONF -1028 ///< Can not open file with mask config.
#define DPE_ERR_MASK_INIT -1029 ///< Error occurred during mask initialization.
#define DPE_ERR_PARAM_VAL -1030 ///< Parameter is missing value in correct format. Check e.g. quotations.
#define DPE_ERR_PARAM_NON_EXIST -1031 ///< Non existing parameter name.
#define DPE_ERR_COMPTON_INIT -1032 ///< Error occurred during Compton camera module initialization.
#define DPE_ERR_DIRECTION_INIT -1033 ///< Error occurred during direction analysis module initialization.
#define DPE_ERR_COMPTON_INDEX -1034 ///< Error occurred during Compton camera indexes mapping.
#define DPE_ERR_FILTER_INIT -1035 ///< Error occurred during filter module initialization.
#define DPE_ERR_INIT -1036 ///< Error occurred during initialization.
#define DPE_ERR_INIT_OUT_DIR -1037 ///< Output directory does not exits.
#define DPE_ERR_INIT_NO_IN_FILES -1038 ///< No in files were found for processing.
#define DPE_ERR_PID_INIT -1039 ///< Error occurred during PID module initialization.
#define DPE_ERR_HISTF_INIT -1040 ///< Error occurred during hist process file initialization.
#define DPE_ERR_SIGVEC_INIT -1041 ///< Error occurred during significant vector initialization.
#define DPE_ERR_COMPAR_NO_DEFAULT -1042 ///< Can not used default .
#define DPE_ERR_DIR_EXPORT -1043 ///< Error occurred during creation of export directories.
#define DPE_ERR_DET_GEO_INIT -1044 ///< Error occurred during detector geometry initialization.
#define DPE_ERR_DIS_GRAPHICS -1045 ///< Creating of graphics disabled, some/all python modules are missing.
#define DPE_ERR_DIS_MULTIPROC -1046 ///< Creating of graphics disabled, some/all python modules are missing.
#define DPE_ERR_PREREQ -1047 ///< Error occurred during checking of prerequisite.
#define DPE_ERR_CAL_MAT -1048 ///< Directory with calibration matrices can not be found. Data will not be calibrated.
#define DPE_ERR_NO_MASK_FILE -2000 ///< File with mask can not be opened.
#define DPE_ERR_MASK_LOAD -2001 ///< Mask load failed.
#define DPE_ERR_LOAD_CLOG_TACQ -2100 ///< Can not load acq time of frames from input clog.
#define DPE_ERR_NO_ELIST -3000 ///< Elist file can not be found/opened.
#define DPE_ERR_NO_EELIST -3001 ///< ExtElist can not be opened.
#define DPE_ERR_RFR -3003 ///< Error occurred during radiation field recognition.
#define DPE_ERR_SIGVEC_HIST -3004 ///< Creating significant vector failed because not all needed histograms were given.
#define DPE_ERR_SIGVEC -3005 ///< Error occurred during significant vector creation.
#define DPE_ERR_PROC -3006 ///< Error occurred during data processing.
#define DPE_ERR_CLPROC_NOCLS -3007 ///< No clusters given for processing.
#define DPE_ERR_CLPROC -3008 ///< Error occurred during clusters processing.
#define DPE_ERR_CLPROC_PID -3009 ///< Error occurred during PID processing.
#define DPE_ERR_CLPROC_FILTER -3010 ///< Error occurred during filtering.
#define DPE_ERR_CLPROC_DIR -3011 ///< Error occurred during dir analysis.
#define DPE_ERR_CLPROC_FRAME -3012 ///< Error occurred during frame analysis.
#define DPE_ERR_ELIST_INIT -3013 ///< Elist init failed.
#define DPE_ERR_PROC_ELIST_init -3014 ///< Error occurred during elist processing init.
#define DPE_ERR_PROC_ELIST -3015 ///< Error occurred during elist processing.
#define DPE_ERR_EXP_NO_DATA -4000 ///< No data was processed. Can not continue with export.
#define DPE_ERR_RELOAD_CLOG -4011 ///< Can not open file with clog for clusters and sensor plots.
#define DPE_ERR_EMPTY_CL -4012 ///< Cluster list is empty for clusters and sensor plots.
#define DPE_ERR_GRAPH_MULTIP_CDIR -4013 ///< Exporting graphics was not successful in the multiprocessing part (open curr dir).
#define DPE_ERR_GRAPH_MULTIP_0FUNC -4014 ///< Exporting graphics was not successful in the multiprocessing part (no sub function for processing).
#define DPE_ERR_GRAPH_MULTIP_RFUNC -4015 ///< Exporting graphics was not successful in the multiprocessing part (can not read sub function file).
#define DPE_ERR_RELOAD_CLOG_OPENF -4016 ///< Can not open file with clog for clusters and sensor plots.
#define DPE_ERR_RELOAD_CLOG_NUPIX -4017 ///< Incorrect number of pixel data in line with cluster.
#define DPE_ERR_CLUSTMATRIX_EXPCL -4018 ///< Can not export coincidence group because the storage is empty.
#define DPE_ERR_EXPORT -4022 ///< Error occurred during export.
#define DPE_ERR_DEL_ELIST -4023 ///< Can not delete elist file.
#define DPE_ERR_FRAME_ANALYSIS_EXP -4024 ///< Error occurred during frame analysis export.
#define DPE_ERR_COINC_ANALYSIS_EXP -4025 ///< Error occurred during coinc event analysis export.
#define DPE_ERR_RFR_EXP -4026 ///< Error occurred during radiation field recognition.
#define DPE_ERR_DIR_ANALYSIS_EXP -4027 ///< Error occurred during directional analysis export.
#define DPE_ERR_COMCAM_ANALYSIS_EXP -4028 ///< Error occurred during Compton analysis export.
Application Programing Interface - C++
…
Application Programing Interface - python
…
Post-processing Library in Python
…
Creators, Contributors and References
List of creators and contributors:
- Lukas Marek
- Carlos Granja
- Jan Jakubek
- Daniel Turecek
- Jan Ingerle
- Cristina Oancea
- Daniela Doubravova
- Jakub Kolarik
- Jacob Miller
Citation of DPE:
…
References:
- ↑ J. Jakubek, ‘Precise energy calibration of pixel detector working in time-over-threshold mode’, Nucl. Instrum. Methods Phys. Res. Sect. Accel. Spectrometers Detect. Assoc. Equip., vol. 633, pp. S262–S266, May 2011, doi: 10.1016/j.nima.2010.06.183.
- ↑ M. Sommer, C. Granja, S. Kodaira, and O. Ploc, ‘High-energy per-pixel calibration of timepix pixel detector with laboratory alpha source’, Nucl. Instrum. Methods Phys. Res. Sect. Accel. Spectrometers Detect. Assoc. Equip., vol. 1022, p. 165957, Jan. 2022, doi: 10.1016/j.nima.2021.165957.
- ↑ B. Hartmann et al., ‘Distortion of the per-pixel signal in the Timepix detector observed in high energy carbon ion beams’, J. Instrum., vol. 9, no. 09, pp. P09006–P09006, Sep. 2014, doi: 10.1088/1748-0221/9/09/P09006.
- ↑ D. Turecek, J. Jakubek, and P. Soukup, ‘USB 3.0 readout and time-walk correction method for Timepix3 detector’, J. Instrum., vol. 11, no. 12, pp. C12065–C12065, Dec. 2016, doi: 10.1088/1748-0221/11/12/C12065.