dE/dx electronic ADC non-linearity correction for highly ionising particles (used in offline hadron saturation calibration) parameters are for X vs Y relation and sep for inner and outer layer vector array 0,1 for inner and 2,3 for outer layers More...

#include <CDCDedxADCNonLinearity.h>

Inheritance diagram for CDCDedxADCNonLinearity:

Public Member Functions
	CDCDedxADCNonLinearity ()
	Default constructor.

	CDCDedxADCNonLinearity (const std::vector< std::vector< double > > &nonlinearity)
	Constructor.

	~CDCDedxADCNonLinearity ()
	Destructor.

unsigned int	getSize (int layer, int axis) const
	Get the number of bins for the non-linearity angle correction.

double	getCorrectedADC (double ADC, int layer) const
	Return corrected ADC with given parameters.

double	getNonLinearityPar (int layer, int axis, unsigned int par) const
	return specific hadron parameter

void	printNonLinearityPars (int layer, int axis) const
	print requested hadron parameter array

void	setNonLinearityPar (unsigned int layer, unsigned int axis, unsigned int par, double value)
	set specific hadron parameter

Private Member Functions
	ClassDef (CDCDedxADCNonLinearity, 1)
	ClassDef.

Private Attributes
std::vector< std::vector< double > >	m_nonlADC
	ADC vs corrected ADC mapping.

Detailed Description

dE/dx electronic ADC non-linearity correction for highly ionising particles (used in offline hadron saturation calibration) parameters are for X vs Y relation and sep for inner and outer layer vector array 0,1 for inner and 2,3 for outer layers

Definition at line 29 of file CDCDedxADCNonLinearity.h.

Constructor & Destructor Documentation

◆ CDCDedxADCNonLinearity() [1/2]

CDCDedxADCNonLinearity ( )

inline

Default constructor.

Definition at line 35 of file CDCDedxADCNonLinearity.h.

35: m_nonlADC() {};

◆ CDCDedxADCNonLinearity() [2/2]

CDCDedxADCNonLinearity ( const std::vector< std::vector< double > > & nonlinearity )

inlineexplicit

Constructor.

Definition at line 40 of file CDCDedxADCNonLinearity.h.

40: m_nonlADC(nonlinearity) {};

◆ ~CDCDedxADCNonLinearity()

~CDCDedxADCNonLinearity ( )

inline

Destructor.

Definition at line 45 of file CDCDedxADCNonLinearity.h.

45{};

Member Function Documentation

◆ getCorrectedADC()

double getCorrectedADC	(	double	ADC,
		int	layer ) const

inline

Return corrected ADC with given parameters.

Parameters

ADC	uncorrected version
layer	outer vs inner 0 to 56 indexing

Definition at line 79 of file CDCDedxADCNonLinearity.h.

    {
 
      unsigned int mylayer = 0;
      if (layer < 0 || layer > 55) {
        B2WARNING("CDCDedxADCNonLinearity:returning uncorrected ADC as input layer is out of range: must be 0-55");
        return ADC;
      } else {
        if (layer >= 8) mylayer = 2;
      }
 
      if (mylayer + 1 >= m_nonlADC.size()) return ADC;
 
      const std::vector<double>& tempX = m_nonlADC[mylayer];//inner or outer X array
      const std::vector<double>& tempY = m_nonlADC[mylayer + 1];//inner or outer Y array
 
      if (tempX.size() < 2 or tempX.size() != tempY.size()) {
        B2WARNING("CDCDedxADCNonLinearity:returning uncorrected ADC as parameters range don't match: X=Y in bins");
        return ADC;
      }
 
      //Find bin for ADC correction
      unsigned int ibin = TMath::BinarySearch(tempY.size(), tempY.data(), double(ADC));
 
      if (ibin >= tempY.size() - 1) ibin = tempY.size() - 2; //overflow to last bin
      if (ibin >= tempY.size() - 1) {
        B2WARNING("CDCDedxADCNonLinearity:returning uncorrected ADC as bins are not in range");
        return ADC;
      }
 
      double slope = (tempY[ibin + 1] - tempY[ibin]) / (tempX[ibin + 1] - tempX[ibin]);
      return  std::round(tempX[ibin] + (ADC - tempY[ibin]) / slope);
    }

◆ getNonLinearityPar()

double getNonLinearityPar	(	int	layer,
		int	axis,
		unsigned int	par ) const

inline

return specific hadron parameter

Parameters

layer	outer(8-55) vs inner(0-7)
axis	0 for X and 1 for Y
par	number starts from 0

Definition at line 118 of file CDCDedxADCNonLinearity.h.

    {
 
      unsigned int mylayer = 0;
      if (layer < 0 || layer > 55) {
        B2ERROR("CDCDedxADCNonLinearity:CDC layer out of range, choose between 0-55");
        return -99.0;
      } else {
        if (layer >= 8) mylayer = 2;
      }
 
      if (axis < 0 || axis > 1) {
        B2ERROR("CDCDedxADCNonLinearity:Problem with X/Y axis: choose 0/1 for X/Y");
        return -99.0;
      } else {
        mylayer = mylayer + axis;
      }
 
      if (mylayer >= m_nonlADC.size()) {
        B2ERROR("CDCDedxADCNonLinearity: vector-of-vectors too short");
        return -99.0;
      }
 
      if (par >= m_nonlADC[mylayer].size()) {
        B2ERROR("CDCDedxADCNonLinearity:Problem with par index: choose 0 and " << m_nonlADC[mylayer].size()); //
        return -99.0;
      }
 
      return m_nonlADC[mylayer][par];
    };

◆ getSize()

unsigned int getSize	(	int	layer,
		int	axis ) const

inline

Get the number of bins for the non-linearity angle correction.

Parameters

layer	outer vs inner 0 to 56 indexing
axis	0 for X and 1 for Y

Definition at line 52 of file CDCDedxADCNonLinearity.h.

    {
 
      unsigned int mylayer = 0;
      if (layer < 0 || layer > 55) {
        B2ERROR("CDCDedxADCNonLinearity:CDC layer out of range, choose between 0-55");
        return 0;
      } else {
        if (layer >= 8) mylayer = 2;
      }
 
      if (axis < 0 || axis > 1) {
        B2ERROR("CDCDedxADCNonLinearity:Problem with X/Y axis: choose 0/1 for X/Y");
        return 0;
      } else {
        mylayer = mylayer + axis;
      }
 
      if (mylayer < m_nonlADC.size()) return m_nonlADC[mylayer].size();
      return 0;
    };

◆ printNonLinearityPars()

void printNonLinearityPars	(	int	layer,
		int	axis ) const

inline

print requested hadron parameter array

Parameters

layer	outer vs inner
axis	0 for X and 1 for Y

Definition at line 153 of file CDCDedxADCNonLinearity.h.

    {
 
      unsigned int mylayer = 0;
      if (layer < 0 || layer > 55) {
        B2ERROR("CDCDedxADCNonLinearity:CDC layer out of range, choose between 0-55");
        return;
      } else {
        if (layer >= 8) mylayer = 2;
      }
 
      if (axis < 0 || axis > 1) {
        B2ERROR("CDCDedxADCNonLinearity:Problem with X/Y axis: choose 0/1 for X/Y");
        return;
      } else {
        mylayer = mylayer + axis;
      }
 
      if (mylayer >= m_nonlADC.size()) {
        B2ERROR("CDCDedxADCNonLinearity: vector-of-vectors too short");
        return;
      }
 
      B2INFO("Printing parameters for layer: " << layer << ", axis: " << axis << ", nPars:" << m_nonlADC[mylayer].size());
      for (unsigned int iPar = 0; iPar < m_nonlADC[mylayer].size(); iPar++)
        std::cout << "Par # " << iPar << ": " << m_nonlADC[mylayer][iPar] << std::endl;
 
    };

◆ setNonLinearityPar()

void setNonLinearityPar	(	unsigned int	layer,
		unsigned int	axis,
		unsigned int	par,
		double	value )

inline

set specific hadron parameter

Parameters

layer	outer(8-55) vs inner(0-7)
axis	0 for X and 1 for Y
par	number starts from 0
value	of parameter to set

Definition at line 188 of file CDCDedxADCNonLinearity.h.

    {
      unsigned int mylayer = 0;
      if (layer > 55) {
        B2ERROR("CDCDedxADCNonLinearity:CDC layer out of range, choose between 0-55");
        return;
      } else {
        if (layer >= 8) mylayer = 2;
      }
 
      if (axis > 1) {
        B2ERROR("CDCDedxADCNonLinearity:Problem with X/Y axis: choose 0/1 for X/Y");
        return;
      } else {
        mylayer = mylayer + axis;
      }
 
      if (mylayer >= m_nonlADC.size()) {
        B2ERROR("CDCDedxADCNonLinearity: vector-of-vectors too short");
        return;
      }
      if (par >= m_nonlADC[mylayer].size()) {
        B2ERROR("CDCDedxADCNonLinearity:Problem with par index: choose 0 and " << m_nonlADC[mylayer].size()); //
        return;
      }
 
      //here set parameter to requested value
      m_nonlADC[mylayer][par] = value;
 
    };

Member Data Documentation

◆ m_nonlADC

std::vector<std::vector<double> > m_nonlADC

private

ADC vs corrected ADC mapping.

Definition at line 221 of file CDCDedxADCNonLinearity.h.

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

cdc/dbobjects/include/CDCDedxADCNonLinearity.h

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ CDCDedxADCNonLinearity() [1/2]

◆ CDCDedxADCNonLinearity() [2/2]

◆ ~CDCDedxADCNonLinearity()

Member Function Documentation

◆ getCorrectedADC()

◆ getNonLinearityPar()

◆ getSize()

◆ printNonLinearityPars()

◆ setNonLinearityPar()

Member Data Documentation

◆ m_nonlADC