This object contains ECL DSP coefs – electromagnetic calorimeter digital signal processing coefficients. More...

#include <ECLDspData.h>

Inheritance diagram for ECLDspData:

Collaboration diagram for ECLDspData:

Public Member Functions
	ECLDspData (unsigned char board_number=0)
	Initialize DSP coefficients class from binary file. More...

std::vector< short int >	getExtraData () const

void	getF (std::vector< short int > &dst) const
	Array with tabulated signal waveform.

void	getF1 (std::vector< short int > &dst) const
	Array with tabulated derivative of signal waveform.

void	getF31 (std::vector< short int > &dst) const
	Array FG31, used to estimate signal amplitude. More...

void	getF32 (std::vector< short int > &dst) const
	Array FG32, used to estimate A * delta_t. More...

void	getF33 (std::vector< short int > &dst) const
	Array FG33, used to estimate pedestal height in signal. More...

void	getF41 (std::vector< short int > &dst) const
	Alternative for FG31 for signals with small amplitude.

void	getF43 (std::vector< short int > &dst) const
	Alternative for FG33 for signals with small amplitude.

unsigned char	getverMaj () const

unsigned char	getverMin () const

short int	getaAT () const

short int	getlAT () const

short int	getsT () const

short int	gethT () const

short int	getchiThresh () const
	chi2 threshold for fit quality flag

unsigned char	getk1 () const
	multipliers power of 2 for f, f1

unsigned char	getk2 () const
	multipliers power of 2 for chi2 calculation

unsigned char	getka () const
	Number of bits for FG31, FG41.

unsigned char	getkb () const
	Number of bits for FG32.

unsigned char	getkc () const
	Number of bits for FG33, FG43.

unsigned char	gety0Startr () const
	start point for pedestal calculation

unsigned char	getBoardNumber () const
	Return ShaperDSP board number, 0..11.

void	setExtraData (const std::vector< short int > &newval)
	SETTERS. More...

void	setF (std::vector< short int > &src)
	Array with tabulated signal waveform.

void	setF1 (std::vector< short int > &src)
	Array with tabulated derivative of signal waveform.

void	setF31 (std::vector< short int > &src)
	Array FG31, used to estimate signal amplitude. More...

void	setF32 (std::vector< short int > &src)
	Array FG32, used to estimate A * delta_t. More...

void	setF33 (std::vector< short int > &src)
	Array FG33, used to estimate pedestal height in signal. More...

void	setF41 (std::vector< short int > &src)
	Alternative for FG31 for signals with small amplitude.

void	setF43 (std::vector< short int > &src)
	Alternative for FG33 for signals with small amplitude.

void	setverMaj (unsigned char val)
	Sets major version of DSP coefficients.

void	setverMin (unsigned char val)
	Sets minor version of DSP coefficients.

void	setaAT (short int val)
	Set Low amp threshold (https://confluence.desy.de/display/BI/Electronics+Thresholds)

void	setlAT (short int val)
	Set Low amp threshold (https://confluence.desy.de/display/BI/Electronics+Thresholds)

void	setsT (short int val)
	Set Skip threshold (https://confluence.desy.de/display/BI/Electronics+Thresholds)

void	sethT (short int val)
	Set Hit threshold (https://confluence.desy.de/display/BI/Electronics+Thresholds)

void	setchiThresh (short int val)
	Set chi2 threshold for fit quality flag.

void	setk1 (unsigned char val)
	Set multipliers power of 2 for f, f1.

void	setk2 (unsigned char val)
	Set multipliers power of 2 for chi2 calculation.

void	setka (unsigned char val)
	Set number of bits for FG31, FG41.

void	setkb (unsigned char val)
	Set number of bits for FG32.

void	setkc (unsigned char val)
	Set number of bits for FG33, FG43.

void	sety0Startr (unsigned char val)
	Set start point for pedestal calculation.

Static Public Member Functions
static short int	getPackerVersion ()
	GETTERS. More...

Private Member Functions
void	packCoefVector (const std::vector< short int > &src, std::vector< short int > &dst)
	Convert vector of DSP coefficients (src) to ECLDspData internal format (dst). More...

void	unpackCoefVector (const std::vector< short int > &src, std::vector< short int > &dst) const
	Convert vector of DSP coefficients (src) to ECLDspData internal format (dst). More...

	ClassDef (ECLDspData, 1)
	ClassDef.

Private Attributes
unsigned char	m_boardNumber
	ShaperDSP board number, 0..11.

unsigned char	m_verMaj
	Major version of DSP coefficients.

unsigned char	m_verMin
	Minor version of DSP coefficients.

unsigned char	m_ka
	Number of bits for FG31, FG41.

unsigned char	m_kb
	Number of bits for FG32.

unsigned char	m_kc
	Number of bits for FG33, FG43.

unsigned char	m_y0Startr
	start point for pedestal calculation

short int	m_chiThresh
	chi2 threshold for quality bit More...

unsigned char	m_k1Chi
	multipliers power of 2 for f, f1

unsigned char	m_k2Chi
	multipliers power of 2 for chi2 calculation

short int	m_hitThresh
	See https://confluence.desy.de/display/BI/Electronics+Thresholds.

short int	m_lowAmpThresh
	See https://confluence.desy.de/display/BI/Electronics+Thresholds.

short int	m_skipThresh
	See https://confluence.desy.de/display/BI/Electronics+Thresholds.

short int	m_adcAlwThresh
	See https://confluence.desy.de/display/BI/Electronics+Thresholds.

std::vector< short int >	m_fg31
	See documentation for method 'getF31'.

std::vector< short int >	m_fg32
	See documentation for method 'getF32'.

std::vector< short int >	m_fg33
	See documentation for method 'getF33'.

std::vector< short int >	m_fg41
	See documentation for method 'getF41'.

std::vector< short int >	m_fg43
	See documentation for method 'getF43'.

std::vector< short int >	m_f
	See documentation for method 'getF'.

std::vector< short int >	m_f1
	See documentation for method 'getF1'.

std::vector< short int >	m_extraData
	This vector contains all parameters that didn't exist in the initial version of ECL DSP file format. More...

Static Private Attributes
static const short int	c_packerVersion = 2
	Version of packer used in packCoefVector by default Versioon of unpacker is dynamically determined from the data.

Detailed Description

This object contains ECL DSP coefs – electromagnetic calorimeter digital signal processing coefficients.

Payloads with these objects are used in two use cases.

They are uploaded to ShaperDSP boards and then used in waveform fitting (amplitude and time reconstruction).
ECL DQM uses these payloads with ShaperDSP emulator to quickly detect ShaperDSP boards that started sending corrupted data. (ECL expert then resets these modules)

For detailed overview of ECLDSP file structure see Confluence, https://confluence.desy.de/display/BI/ECL+Technical+Notes Attached file ECL-TN-2013-02 (digi.pdf)

Definition at line 44 of file ECLDspData.h.

Constructor & Destructor Documentation

◆ ECLDspData()

ECLDspData ( unsigned char board_number = 0 )

inlineexplicit

Initialize DSP coefficients class from binary file.

Parameters

board_number Id of specific shaperDSP, 0..11

Definition at line 115 of file ECLDspData.h.

     {}
     ~ECLDspData() {};
  
     /*************/
     /*************/
  

Member Function Documentation

◆ getaAT()

short int getaAT ( ) const

inline

Returns: ADC always threshold (https://confluence.desy.de/display/BI/Electronics+Thresholds)

Definition at line 195 of file ECLDspData.h.

◆ getExtraData()

std::vector<short int> getExtraData ( ) const

inline

Returns: Vector with all parameters that haven't been included in basic ECL DSP file format

Definition at line 152 of file ECLDspData.h.

◆ getF31()

void getF31 ( std::vector< short int > & dst ) const

inline

Array FG31, used to estimate signal amplitude.

Calculated from F and covariance matrix.

Definition at line 166 of file ECLDspData.h.

◆ getF32()

void getF32 ( std::vector< short int > & dst ) const

inline

Array FG32, used to estimate A * delta_t.

(A – amplitude, delta_t – time shift for linearization)

Calculated from F and covariance matrix.

Definition at line 173 of file ECLDspData.h.

◆ getF33()

void getF33 ( std::vector< short int > & dst ) const

inline

Array FG33, used to estimate pedestal height in signal.

Calculated from F and covariance matrix.

Definition at line 179 of file ECLDspData.h.

◆ gethT()

short int gethT ( ) const

inline

Returns: Hit threshold (https://confluence.desy.de/display/BI/Electronics+Thresholds)

Definition at line 207 of file ECLDspData.h.

◆ getlAT()

short int getlAT ( ) const

inline

Returns: Low amp threshold (https://confluence.desy.de/display/BI/Electronics+Thresholds)

Definition at line 199 of file ECLDspData.h.

◆ getPackerVersion()

static short int getPackerVersion ( )

inlinestatic

GETTERS.

Return version of packer algorithm used in [un]packCoefVector

Definition at line 146 of file ECLDspData.h.

◆ getsT()

short int getsT ( ) const

inline

Returns: Skip threshold (https://confluence.desy.de/display/BI/Electronics+Thresholds)

Definition at line 203 of file ECLDspData.h.

◆ getverMaj()

unsigned char getverMaj ( ) const

inline

Returns: Major version of DSP coefficients

Definition at line 189 of file ECLDspData.h.

◆ getverMin()

unsigned char getverMin ( ) const

inline

Returns: Minor version of DSP coefficients

Definition at line 191 of file ECLDspData.h.

◆ packCoefVector()

void packCoefVector	(	const std::vector< short int > &	src,
		std::vector< short int > &	dst
	)

private

Convert vector of DSP coefficients (src) to ECLDspData internal format (dst).

The internal format of coefficient vectors stored in this object is different from their real values: this improves data compression by ~200%. Thus, any accessor methods to m_f* vectors utilize packCoefVector and unpackCoefVector.

Definition at line 5 of file ECLDspData.cc.

 {
   const int N_CHANNELS = 16;
   int size = src.size();
  
   if (getPackerVersion() == 0) {
     dst = src;
   } else if (getPackerVersion() >= 1) {
     // Apply reversible transformation to shrink the range
     // of saved coefficients into much smaller interval.
     // This improves the compression by a factor of ~2.
  
     dst.resize(size);
     for (int i = size - 1; i >= 0; i--) {
       if (i >= 2 * N_CHANNELS)
         dst[i] = src[i] - 2 * src[i - N_CHANNELS] + src[i - 2 * N_CHANNELS];
       else if (i >= N_CHANNELS)
         dst[i] = src[i] - 2 * src[i - N_CHANNELS];
       else
         dst[i] = src[i];
     }
   }
  
   if (getPackerVersion() == 2) {
     // Rearrange the data in such a way that each coeff
     // only takes up 4 bits (special value 0xF indicates
     // coefficients that are too big to be packed)
     //
     // Vector structure after packing:
     // [big_coefs | packed_coefs | packed_coefs.size()]
  
     // Bits allocated for each value
     const int value_bits = 4;
     const long value_max = 0xF;
     const int values_packed = sizeof(short) * 8 / value_bits;
  
     int shift = -6;
  
     std::vector<short> packed(size / values_packed);
     int len = -1;
     for (int i = 0; i < size; i++) {
       if (i % values_packed == 0) packed[i / values_packed] = 0;
  
       short val = dst[i] - shift;
       if (val >= 0 && val < value_max) {
         if (len < 0) len = i;
         packed[i / values_packed] |= val << (value_bits * (i % values_packed));
       } else {
         if (len >= 0 && len < i) {
           dst[len++] = dst[i];
         }
         packed[i / values_packed] |= short(value_max) << (value_bits * (i % values_packed));
       }
     }
  
     int packed_size = 0;
  
     if (len >= 0) {
       dst.resize(len + packed.size());
       packed_size = packed.size();
       for (int i = 0; i < packed_size; i++) {
         dst[len + i] = packed[i];
       }
     }
     dst.push_back(packed_size);
   }
 }

◆ setExtraData()

void setExtraData ( const std::vector< short int > & newval )

inline

SETTERS.

Sets vector with all parameters that haven't been included in basic ECL DSP file format

Definition at line 234 of file ECLDspData.h.

◆ setF31()

void setF31 ( std::vector< short int > & src )

inline

Array FG31, used to estimate signal amplitude.

Calculated from F and covariance matrix.

Definition at line 248 of file ECLDspData.h.

◆ setF32()

void setF32 ( std::vector< short int > & src )

inline

Array FG32, used to estimate A * delta_t.

(A – amplitude, delta_t – time shift for linearization)

Calculated from F and covariance matrix.

Definition at line 255 of file ECLDspData.h.

◆ setF33()

void setF33 ( std::vector< short int > & src )

inline

Array FG33, used to estimate pedestal height in signal.

Calculated from F and covariance matrix.

Definition at line 261 of file ECLDspData.h.

◆ unpackCoefVector()

void unpackCoefVector	(	const std::vector< short int > &	src,
		std::vector< short int > &	dst
	)		const

private

Convert vector of DSP coefficients (src) to ECLDspData internal format (dst).

See documentation for packCoefVector on why this was necessary.

Definition at line 72 of file ECLDspData.cc.

Member Data Documentation

◆ m_chiThresh

short int m_chiThresh

private

chi2 threshold for quality bit

We use condition chi^2 < (Amp^2 / (2^{k2Chi}) + chiThresh)

Definition at line 65 of file ECLDspData.h.

◆ m_extraData

std::vector<short int> m_extraData

private

This vector contains all parameters that didn't exist in the initial version of ECL DSP file format.

This vector allows loose coupling between ECLDspData class and ECL DSP file format.

m_extraData[0] – FORMAT_VERSION m_extraData[1] – PACKER_VERSION

Definition at line 102 of file ECLDspData.h.

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

ecl/dbobjects/include/ECLDspData.h
ecl/dbobjects/src/ECLDspData.cc

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ ECLDspData()

Member Function Documentation

◆ getaAT()

◆ getExtraData()

◆ getF31()

◆ getF32()

◆ getF33()

◆ gethT()

◆ getlAT()

◆ getPackerVersion()

◆ getsT()

◆ getverMaj()

◆ getverMin()

◆ packCoefVector()

◆ setExtraData()

◆ setF31()

◆ setF32()

◆ setF33()

◆ unpackCoefVector()

Member Data Documentation

◆ m_chiThresh

◆ m_extraData