REST-for-Physics  v2.3
Rare Event Searches ToolKit for Physics
Public Member Functions | Private Member Functions | Private Attributes
TRestRawSignalAnalysisProcess Class Reference

Detailed Description

An analysis process to extract valuable information from a TRestRawSignalEvent.

THIS DOCUMENTATION NEEDS REVISION. TODO If possible a figure with the pulse ilustrating the different observables extracted.

The following metadata parameters might be defined at the RML to control the rawsignal analysis is performed:

There are 3 additional parameters that are used in combination to identify the points over threshold at each signal. Those parameters are used, for example, to calculate the number of good signals observable.

Additionaly, there is a metadata parameter,signalsRange that allows to define the signal ids over which this process will have effect. This parameter may allow to define different TRestRawSignalAnalysisProcess instances at the RML operating over different signal ranges. For example:

<TRestRawSignalAnalysisProcess name="rawMM" signalsRange="(0,512)"
baseLineRange="(5,100)" observables="all" />
<TRestRawSignalAnalysisProcess name="rawVETO" signalsRange="(513,1024)" >
<parameter name="pointThreshold" value="5" />
<parameter name="signalThreshold" value="3.5" />
<parameter name="pointsOverThreshold" value="15" />
<observable name="rawAna_max_amplitude_map" type="map<int,double>" />
An analysis process to extract valuable information from a TRestRawSignalEvent.

Observables

Number of signals and base line:

Integrals and energy estimations:

Time observables:

Peak amplitude observables:

Peak time observables:

Observables for individual signal info:

You may add filters to any observable inside the analysis tree. To add a cut, write "cut" sections in your rml file:

<cut name="MeanBaseLineCut" value="(0,4096)" />

Warning
⚠ REST is under continuous development. This documentation is offered to you by the REST community. Your HELP is needed to keep this code up to date. Your feedback will be worth to support this software, please report any problems/suggestions you may find will using it at The REST Framework forum. You are welcome to contribute fixing typos, updating information or adding/proposing new contributions. See also our Contribution Guide.

RESTsoft - Software for Rare Event Searches with TPCs

History of developments:

2017-February: First implementation of raw signal analysis process into REST_v2. Created from TRestDetectorSignalAnalysisProcess

Author
Javier Galan

Definition at line 31 of file TRestRawSignalAnalysisProcess.h.

#include <TRestRawSignalAnalysisProcess.h>

Inheritance diagram for TRestRawSignalAnalysisProcess:
TRestEventProcess TRestMetadata

Public Member Functions

 ClassDefOverride (TRestRawSignalAnalysisProcess, 5)
 
RESTValue GetInputEvent () const override
 Get pointer to input event. Must be implemented in the derived class. More...
 
RESTValue GetOutputEvent () const override
 Get pointer to output event. Must be implemented in the derived class. More...
 
const char * GetProcessName () const override
 
void InitFromConfigFile () override
 To make settings from rml file. This method must be implemented in the derived class. More...
 
void InitProcess () override
 Process initialization. More...
 
void PrintMetadata () override
 Implemented it in the derived metadata class to print out specific metadata information. More...
 
TRestEventProcessEvent (TRestEvent *inputEvent) override
 The main processing event function. More...
 
 TRestRawSignalAnalysisProcess ()
 Default constructor. More...
 
 ~TRestRawSignalAnalysisProcess ()
 Default destructor. More...
 
- Public Member Functions inherited from TRestEventProcess
virtual Bool_t AddInputFile (const std::string &file)
 
bool ApplyCut ()
 
virtual void BeginOfEventProcess (TRestEvent *inputEvent=nullptr)
 Begin of event process, preparation work. Called right before ProcessEvent() More...
 
 ClassDefOverride (TRestEventProcess, 3)
 
virtual void EndOfEventProcess (TRestEvent *inputEvent=nullptr)
 End of event process. Nothing to do. Called directly after ProcessEvent() More...
 
virtual void EndProcess ()
 To be executed at the end of the run (outside event loop) More...
 
TRestAnalysisTreeGetAnalysisTree () const
 Return the local analysis tree (dummy) More...
 
TCanvas * GetCanvas () const
 Get canvas. More...
 
TRestAnalysisTreeGetFullAnalysisTree ()
 
virtual RESTValue GetInputEvent () const =0
 Get pointer to input event. Must be implemented in the derived class. More...
 
std::vector< std::string > GetListOfAddedObservables ()
 
virtual RESTValue GetOutputEvent () const =0
 Get pointer to output event. Must be implemented in the derived class. More...
 
virtual const char * GetProcessName () const =0
 
TRestRunGetRunInfo () const
 Return the pointer of the hosting TRestRun object. More...
 
virtual Long64_t GetTotalBytes () const
 
virtual Long64_t GetTotalBytesRead () const
 Interface to external file reading, get the read bytes. To be implemented in external processes. More...
 
virtual void InitFromConfigFile () override
 To make settings from rml file. This method must be implemented in the derived class. More...
 
virtual void InitProcess ()
 To be executed at the beginning of the run (outside event loop) More...
 
Bool_t isExternal () const
 Return whether this process is external process. More...
 
Int_t LoadSectionMetadata () override
 This method does some preparation of xml section. More...
 
virtual void NotifyAnalysisTreeReset ()
 In case the analysis tree is reset(switched to new file), some process needs to have action. More...
 
virtual Bool_t OpenInputFiles (const std::vector< std::string > &files)
 
virtual TRestEventProcessEvent (TRestEvent *inputEvent)=0
 Process one event. More...
 
std::vector< std::string > ReadObservables ()
 
void RegisterAllObservables (Bool_t value=true)
 
virtual Bool_t ResetEntry ()
 
void SetAnalysisTree (TRestAnalysisTree *tree)
 Set analysis tree of this process, then add observables to it.
 
void SetCanvasSize (Int_t x, Int_t y)
 Set canvas size. More...
 
void SetFriendProcess (TRestEventProcess *p)
 Add friendly process to this process.
 
void SetObservableValidation (bool validate)
 
void SetParallelProcess (TRestEventProcess *p)
 Add parallel process to this process.
 
void SetRunInfo (TRestRun *r)
 Set TRestRun for this process. More...
 
Bool_t singleThreadOnly () const
 Return whether this process is single std::thread only. More...
 
 TRestEventProcess ()
 
 ~TRestEventProcess ()
 

Private Member Functions

void Initialize () override
 Function to initialize input/output event members and define the section name. More...
 

Private Attributes

std::string fBaseLineOption
 Option for calculation of baseline parameters, can be set to "ROBUST". More...
 
TVector2 fBaseLineRange = TVector2(5, 55)
 The range where the baseline range will be calculated. More...
 
std::string fChannelType
 avoid error when loading parameters (not used) More...
 
std::set< std::string > fChannelTypes = {}
 
TRestRawSignalEventfInputEvent
 A pointer to the specific TRestRawSignalEvent input. More...
 
TVector2 fIntegralRange = TVector2(10, 500)
 The range where the observables will be calculated. More...
 
Int_t fPointsOverThreshold = 5
 The minimum number of points over threshold to identify a signal as such. More...
 
Double_t fPointThreshold = 2
 The number of sigmas over baseline fluctuations to identify a point over threshold. More...
 
Bool_t fRangeEnabled = false
 Just a flag to quickly determine if we have to apply the range filter. More...
 
TRestRawReadoutMetadatafReadoutMetadata = nullptr
 
TVector2 fSignalsRange = TVector2(-1, -1)
 It defines the signals id range where analysis is applied. More...
 
Double_t fSignalThreshold = 5
 A parameter to define a minimum signal fluctuation. Measured in sigmas. More...
 

Additional Inherited Members

- Protected Types inherited from TRestEventProcess
enum  REST_Process_Output { No_Output , Observable , Internal_Var , Full_Output }
 
- Protected Member Functions inherited from TRestEventProcess
void BeginPrintProcess ()
 [name, cut range]
 
void CreateCanvas ()
 Create the canvas. More...
 
void EndPrintProcess ()
 
TRestEventProcessGetFriend (const std::string &nameOrType)
 
TRestEventProcessGetFriendLive (const std::string &nameOrType)
 
template<class T >
T * GetMetadata ()
 Get a metadata object from the host TRestRun. More...
 
TRestMetadataGetMetadata (const std::string &nameOrType)
 
size_t GetNumberOfParallelProcesses () const
 
template<class T >
GetObservableValue (const std::string &name)
 
TRestEventProcessGetParallel (int i)
 
template<class T >
std::vector< T > GetParallelDataMembers (T *member_of_process)
 Get a list of data members from parallel processes which is same to this process's certain data member. More...
 
template<class T >
void SetObservableValue (const std::string &name, const T &value)
 Set observable value for AnalysisTree. More...
 
- Protected Attributes inherited from TRestEventProcess
TRestAnalysisTreefAnalysisTree = nullptr
 
TCanvas * fCanvas = nullptr
 < Canvas for some viewer event More...
 
TVector2 fCanvasSize
 Canvas size. More...
 
std::vector< std::pair< std::string, TVector2 > > fCuts
 Stores cut definitions. Any listed observables should be in the range. More...
 
bool fDynamicObs = false
 It defines whether to use added observables only or all the observables appear in the code. More...
 
bool fIsExternal = false
 It defines if the process reads event data from an external source. More...
 
std::map< std::string, int > fObservablesDefined
 Stores the list of all the appeared process observables in the code. More...
 
std::map< std::string, int > fObservablesUpdated
 Stores the list of process observables updated when processing this event. More...
 
bool fReadOnly = false
 not used, keep for compatibility More...
 
TRestRunfRunInfo = nullptr
 < Pointer to TRestRun object where to find metadata. More...
 
bool fSingleThreadOnly = false
 
bool fValidateObservables = false
 It defines if observable names should be added to the validation list. More...
 

Constructor & Destructor Documentation

◆ TRestRawSignalAnalysisProcess()

TRestRawSignalAnalysisProcess::TRestRawSignalAnalysisProcess ( )

Default constructor.

Definition at line 259 of file TRestRawSignalAnalysisProcess.cxx.

◆ ~TRestRawSignalAnalysisProcess()

TRestRawSignalAnalysisProcess::~TRestRawSignalAnalysisProcess ( )

Default destructor.

Definition at line 264 of file TRestRawSignalAnalysisProcess.cxx.

Member Function Documentation

◆ GetInputEvent()

RESTValue TRestRawSignalAnalysisProcess::GetInputEvent ( ) const
inlineoverridevirtual

Get pointer to input event. Must be implemented in the derived class.

Implements TRestEventProcess.

Definition at line 73 of file TRestRawSignalAnalysisProcess.h.

◆ GetOutputEvent()

RESTValue TRestRawSignalAnalysisProcess::GetOutputEvent ( ) const
inlineoverridevirtual

Get pointer to output event. Must be implemented in the derived class.

Implements TRestEventProcess.

Definition at line 75 of file TRestRawSignalAnalysisProcess.h.

◆ GetProcessName()

const char * TRestRawSignalAnalysisProcess::GetProcessName ( ) const
inlineoverridevirtual

Implements TRestEventProcess.

Definition at line 99 of file TRestRawSignalAnalysisProcess.h.

◆ InitFromConfigFile()

void TRestRawSignalAnalysisProcess::InitFromConfigFile ( )
overridevirtual

To make settings from rml file. This method must be implemented in the derived class.

Reimplemented from TRestEventProcess.

Definition at line 286 of file TRestRawSignalAnalysisProcess.cxx.

◆ Initialize()

void TRestRawSignalAnalysisProcess::Initialize ( )
overrideprivatevirtual

Function to initialize input/output event members and define the section name.

Reimplemented from TRestMetadata.

Definition at line 270 of file TRestRawSignalAnalysisProcess.cxx.

◆ InitProcess()

void TRestRawSignalAnalysisProcess::InitProcess ( )
overridevirtual

Process initialization.

Reimplemented from TRestEventProcess.

Definition at line 280 of file TRestRawSignalAnalysisProcess.cxx.

◆ PrintMetadata()

void TRestRawSignalAnalysisProcess::PrintMetadata ( )
inlineoverridevirtual

Implemented it in the derived metadata class to print out specific metadata information.

Prints metadata content on screen. Usually overloaded by the derived metadata class.

Reimplemented from TRestMetadata.

Definition at line 83 of file TRestRawSignalAnalysisProcess.h.

◆ ProcessEvent()

TRestEvent * TRestRawSignalAnalysisProcess::ProcessEvent ( TRestEvent inputEvent)
overridevirtual

The main processing event function.

We define (or re-define) the baseline range and calculation range of our raw-signals.

Important call we need to initialize the points over threshold in a TRestRawSignal

Implements TRestEventProcess.

Definition at line 297 of file TRestRawSignalAnalysisProcess.cxx.

Field Documentation

◆ fBaseLineOption

std::string TRestRawSignalAnalysisProcess::fBaseLineOption
private

Option for calculation of baseline parameters, can be set to "ROBUST".

Definition at line 46 of file TRestRawSignalAnalysisProcess.h.

◆ fBaseLineRange

TVector2 TRestRawSignalAnalysisProcess::fBaseLineRange = TVector2(5, 55)
private

The range where the baseline range will be calculated.

Definition at line 40 of file TRestRawSignalAnalysisProcess.h.

◆ fChannelType

std::string TRestRawSignalAnalysisProcess::fChannelType
private

avoid error when loading parameters (not used)

Definition at line 55 of file TRestRawSignalAnalysisProcess.h.

◆ fChannelTypes

std::set<std::string> TRestRawSignalAnalysisProcess::fChannelTypes = {}
private

Definition at line 67 of file TRestRawSignalAnalysisProcess.h.

◆ fInputEvent

TRestRawSignalEvent* TRestRawSignalAnalysisProcess::fInputEvent
private

A pointer to the specific TRestRawSignalEvent input.

Definition at line 34 of file TRestRawSignalAnalysisProcess.h.

◆ fIntegralRange

TVector2 TRestRawSignalAnalysisProcess::fIntegralRange = TVector2(10, 500)
private

The range where the observables will be calculated.

Definition at line 43 of file TRestRawSignalAnalysisProcess.h.

◆ fPointsOverThreshold

Int_t TRestRawSignalAnalysisProcess::fPointsOverThreshold = 5
private

The minimum number of points over threshold to identify a signal as such.

Definition at line 58 of file TRestRawSignalAnalysisProcess.h.

◆ fPointThreshold

Double_t TRestRawSignalAnalysisProcess::fPointThreshold = 2
private

The number of sigmas over baseline fluctuations to identify a point over threshold.

Definition at line 49 of file TRestRawSignalAnalysisProcess.h.

◆ fRangeEnabled

Bool_t TRestRawSignalAnalysisProcess::fRangeEnabled = false
private

Just a flag to quickly determine if we have to apply the range filter.

Definition at line 37 of file TRestRawSignalAnalysisProcess.h.

◆ fReadoutMetadata

TRestRawReadoutMetadata* TRestRawSignalAnalysisProcess::fReadoutMetadata = nullptr
private

Definition at line 63 of file TRestRawSignalAnalysisProcess.h.

◆ fSignalsRange

TVector2 TRestRawSignalAnalysisProcess::fSignalsRange = TVector2(-1, -1)
private

It defines the signals id range where analysis is applied.

Definition at line 61 of file TRestRawSignalAnalysisProcess.h.

◆ fSignalThreshold

Double_t TRestRawSignalAnalysisProcess::fSignalThreshold = 5
private

A parameter to define a minimum signal fluctuation. Measured in sigmas.

Definition at line 52 of file TRestRawSignalAnalysisProcess.h.


The documentation for this class was generated from the following files: