Database object for space resolutions. More...

#include <CDCSpaceResols.h>

Inheritance diagram for CDCSpaceResols:

Collaboration diagram for CDCSpaceResols:

[legend]

Public Types
enum	{ c_nSLayers = 56 , c_maxNAlphaBins = 18 , c_maxNThetaBins = 7 , c_maxNSigmaParams = 9 }
	Constants. More...

Public Member Functions
	CDCSpaceResols ()
	Default constructor.

void	setAlphaBin (const array3 &alpha)
	Set alpha-angle bin (rad)

void	setThetaBin (const array3 &theta)
	Set theta-angle bin (rad)

void	setSigmaParamMode (unsigned short mode)
	Set sigma parameterization mode.

void	setMaxSpaceResol (float sigma)
	Set max.

void	setSigmaParams (const SigmaID sigmaID, const std::vector< float > &params)
	Set sigma parameters for the specified id.

void	setSigmaParams (unsigned short iCLayer, unsigned short iLR, unsigned short iAlpha, unsigned short iTheta, const std::vector< float > &params)
	Set sigma parameters for the specified id.

void	addSigmaParams (const SigmaID sigmaID, const std::vector< float > &delta)
	Update parameters for the specified id.

void	addSigmaParams (unsigned short iCLayer, unsigned short iLR, unsigned short iAlpha, unsigned short iTheta, const std::vector< float > &delta)
	Update parameters for the specified id.

unsigned short	getNoOfAlphaBins () const
	Get no.

unsigned short	getNoOfThetaBins () const
	Get no.

const array3 &	getAlphaBin (unsigned short i) const
	Get i-th alpha-angle bin info.

float	getAlphaPoint (unsigned short i) const
	Get i-th alpha-angle point (rad)

const array3 &	getThetaBin (unsigned short i) const
	Get i-th theta-angle bin info.

float	getThetaPoint (unsigned short i) const
	Get i-th theta-angle point (rad)

unsigned short	getSigmaParamMode () const
	Get parameterization mode.

float	getMaxSpaceResol () const
	Get max.

SigmaID	getSigmaID (unsigned short iCLayer, unsigned short iLR, unsigned short iAlpha, unsigned short iTheta) const
	Get id.

SigmaID	getSigmaID (unsigned short iCLayer, unsigned short iLR, float alpha, float theta) const
	Get id.

const std::vector< float > &	getSigmaParams (const SigmaID sigmaID) const
	Get sigma parameters for the specified id.

const std::vector< float > &	getSigmaParams (unsigned short iCLayer, unsigned short iLR, unsigned short iAlpha, unsigned short iTheta) const
	Get sigma parameters for the specified id.

void	dump () const
	Print all contents.

void	outputToFile (std::string fileName) const
	Output the contents in text file format.

Static Public Member Functions
static bool	comp (const array3 &lhs, const array3 &rhs)
	Static function for sorting.

Private Types
typedef std::array< float, 3 >	array3
	angle bin info.

typedef unsigned short	SigmaID
	id.

Private Member Functions
	ClassDef (CDCSpaceResols, 2)
	ClassDef.

Private Attributes
unsigned short	m_sigmaParamMode
	Mode for sigma parameterization.

unsigned short	m_nSigmaParams
	no.

float	m_maxSpaceResol = 0.07
	max.

std::vector< array3 >	m_alphaBins
	alpha bins for sigma (rad)

std::vector< array3 >	m_thetaBins
	theta bins for sigma (rad)

std::map< SigmaID, std::vector< float > >	m_sigmas
	Sigma coefficients for each layer, Left/Right, entrance angle and polar angle.

Detailed Description

Database object for space resolutions.

Definition at line 26 of file CDCSpaceResols.h.

Member Typedef Documentation

◆ array3

typedef std::array<float, 3> array3

private

angle bin info.

Definition at line 28 of file CDCSpaceResols.h.

◆ SigmaID

typedef unsigned short SigmaID

private

id.

for resol.

Definition at line 29 of file CDCSpaceResols.h.

Member Enumeration Documentation

◆ anonymous enum

anonymous enum

Constants.

Enumerator
c_nSLayers	no. of layers
c_maxNAlphaBins	max. no. of alpha angle bins
c_maxNThetaBins	max. no. of theta angle bins
c_maxNSigmaParams	max. no. of params. for sigma

Definition at line 35 of file CDCSpaceResols.h.

         {c_nSLayers      = 56, 
          c_maxNAlphaBins = 18, 
          c_maxNThetaBins =  7, 
          c_maxNSigmaParams =  9  
         };

Constructor & Destructor Documentation

◆ CDCSpaceResols()

CDCSpaceResols ( )

inline

Default constructor.

Definition at line 44 of file CDCSpaceResols.h.

44 : m_sigmaParamMode(0), m_nSigmaParams(0)

45 {}

Belle2::CDCSpaceResols::m_nSigmaParams

unsigned short m_nSigmaParams

no.

Definition: CDCSpaceResols.h:400

Belle2::CDCSpaceResols::m_sigmaParamMode

unsigned short m_sigmaParamMode

Mode for sigma parameterization.

Definition: CDCSpaceResols.h:399

Member Function Documentation

◆ addSigmaParams() [1/2]

void addSigmaParams	(	const SigmaID	sigmaID,
		const std::vector< float > &	delta
	)

inline

Update parameters for the specified id.

(=bin)

Definition at line 127 of file CDCSpaceResols.h.

    {
      std::map<SigmaID, std::vector<float>>::iterator it = m_sigmas.find(sigmaID);
 
      if (it != m_sigmas.end()) {
        for (unsigned short i = 0; i < m_nSigmaParams; ++i) {
          (it->second)[i] += delta[i];
        }
      } else {
        B2FATAL("Specified params not found in addSigmaParams !");
      }
    }

◆ addSigmaParams() [2/2]

void addSigmaParams	(	unsigned short	iCLayer,
		unsigned short	iLR,
		unsigned short	iAlpha,
		unsigned short	iTheta,
		const std::vector< float > &	delta
	)

inline

Update parameters for the specified id.

(=bin)

Definition at line 143 of file CDCSpaceResols.h.

    {
      const SigmaID sigmaID = getSigmaID(iCLayer, iLR, iAlpha, iTheta);
      addSigmaParams(sigmaID, delta);
    }

◆ comp()

static bool comp	(	const array3 &	lhs,
		const array3 &	rhs
	)

inlinestatic

Static function for sorting.

Definition at line 77 of file CDCSpaceResols.h.

    {
      return lhs[0] < rhs[0];
    }

◆ dump()

void dump ( ) const

inline

Print all contents.

Definition at line 307 of file CDCSpaceResols.h.

    {
      std::cout << " " << std::endl;
      std::cout << "Contents of sigma db" << std::endl;
      std::cout << "alpha bins" << std::endl;
 
      const double deg = 180. / M_PI;
 
      unsigned short nAlphaBins = m_alphaBins.size();
      for (unsigned short i = 0; i < nAlphaBins; ++i) {
        std::cout << " " << deg* m_alphaBins[i][0] << " " << deg* m_alphaBins[i][1] << " " << deg* m_alphaBins[i][2] << " " << std::endl;
      }
 
      std::cout << " " << std::endl;
      std::cout << "theta bins" << std::endl;
 
      unsigned short nThetaBins = m_thetaBins.size();
      for (unsigned short i = 0; i < nThetaBins; ++i) {
        std::cout << " " << deg* m_thetaBins[i][0] << " " << deg* m_thetaBins[i][1] << " " << deg* m_thetaBins[i][2] << " " << std::endl;
      }
 
      std::cout << " " << std::endl;
      std::cout << "coefficients for sigma" << std::endl;
 
      for (unsigned short iT = 0; iT < nThetaBins; ++iT) {
        for (unsigned short iA = 0; iA < nAlphaBins; ++iA) {
          for (unsigned short iCL = 0; iCL < c_nSLayers; ++iCL) {
            for (unsigned short iLR = 0; iLR < 2; ++iLR) {
              unsigned short iLRp = abs(iLR - 1);
              std::cout << iCL << "  " << deg* m_thetaBins[iT][2] << "  " << deg* m_alphaBins[iA][2] << "  " << iLRp;
              const std::vector<float> params = getSigmaParams(iCL, iLRp, iA, iT);
              for (unsigned short i = 0; i < m_nSigmaParams; ++i) {
                std::cout << "  " << params[i];
              }
              std::cout << "  " << std::endl;
            }
          }
        }
      }
    }

◆ getAlphaBin()

const array3 & getAlphaBin ( unsigned short i ) const

inline

Get i-th alpha-angle bin info.

(rad)

Definition at line 170 of file CDCSpaceResols.h.

    {
      return m_alphaBins[i];
    }

◆ getAlphaPoint()

float getAlphaPoint ( unsigned short i ) const

inline

Get i-th alpha-angle point (rad)

Definition at line 178 of file CDCSpaceResols.h.

    {
      return m_alphaBins[i][2];
    }

◆ getMaxSpaceResol()

float getMaxSpaceResol ( ) const

inline

Get max.

allowed space resolution

Definition at line 210 of file CDCSpaceResols.h.

    {
      return m_maxSpaceResol;
    }

◆ getNoOfAlphaBins()

unsigned short getNoOfAlphaBins ( ) const

inline

Get no.

of alpha-angle bin

Definition at line 154 of file CDCSpaceResols.h.

    {
      return m_alphaBins.size();
    }

◆ getNoOfThetaBins()

unsigned short getNoOfThetaBins ( ) const

inline

Get no.

of theta-angle bin

Definition at line 162 of file CDCSpaceResols.h.

    {
      return m_thetaBins.size();
    }

◆ getSigmaID() [1/2]

SigmaID getSigmaID	(	unsigned short	iCLayer,
		unsigned short	iLR,
		float	alpha,
		float	theta
	)		const

inline

Get id.

(=bin id.)

Definition at line 231 of file CDCSpaceResols.h.

    {
      /*
      unsigned short iTheta = 999;
      unsigned short ibin = 0;
      for (std::vector<array3>::const_iterator it = m_thetaBins.begin(); it != m_thetaBins.end(); ++it) {
      if ((*it)[0] <= theta && theta <= (*it)[1]) {
      iTheta = ibin;
      break;
      }
      ++ibin;
      }
      */
      unsigned short iTheta = 999;
      unsigned short ibin = 0;
      for (auto const& it : m_thetaBins) {
        if (it[0] <= theta && theta <= it[1]) {
          iTheta = ibin;
          break;
        }
        ++ibin;
      }
      if (iTheta == 999) B2FATAL("Theta bin not found !");
 
      /*
      unsigned short iAlpha = 999;
      ibin = 0;
      for (std::vector<array3>::const_iterator it = m_alphaBins.begin(); it != m_alphaBins.end(); ++it) {
      if ((*it)[0] <= alpha && alpha <= (*it)[1]) {
      iAlpha = ibin;
      break;
      }
      ++ibin;
      }
      */
      unsigned short iAlpha = 999;
      ibin = 0;
      for (auto const& it : m_alphaBins) {
        if (it[0] <= alpha && alpha <= it[1]) {
          iAlpha = ibin;
          break;
        }
        ++ibin;
      }
      if (iAlpha == 999) B2FATAL("Alpha bin not found !");
 
      return getSigmaID(iCLayer, iLR, iAlpha, iTheta);
    }

◆ getSigmaID() [2/2]

SigmaID getSigmaID	(	unsigned short	iCLayer,
		unsigned short	iLR,
		unsigned short	iAlpha,
		unsigned short	iTheta
	)		const

inline

Get id.

(=bin id.) bit#00 - 05: used for layer no. bit#06 - 06: used for left/right bit#07 - 11: used for alpha angle bin bit#12 - 15: used for theta angle bin

Definition at line 222 of file CDCSpaceResols.h.

    {
      SigmaID id = iCLayer + 64 * iLR + 128 * iAlpha + 4096 * iTheta;
      return id;
    }

◆ getSigmaParamMode()

unsigned short getSigmaParamMode ( ) const

inline

Get parameterization mode.

Definition at line 202 of file CDCSpaceResols.h.

    {
      return m_sigmaParamMode;
    }

◆ getSigmaParams() [1/2]

const std::vector< float > & getSigmaParams ( const SigmaID sigmaID ) const

inline

Get sigma parameters for the specified id.

(=bin)

Definition at line 283 of file CDCSpaceResols.h.

    {
      std::map<SigmaID, std::vector<float>>::const_iterator it = m_sigmas.find(sigmaID);
      if (it != m_sigmas.end()) {
        return it->second;
      } else {
        B2FATAL("Specified params. not found in getSigmaParams !");
      }
    }

◆ getSigmaParams() [2/2]

const std::vector< float > & getSigmaParams	(	unsigned short	iCLayer,
		unsigned short	iLR,
		unsigned short	iAlpha,
		unsigned short	iTheta
	)		const

inline

Get sigma parameters for the specified id.

(=bin)

Definition at line 296 of file CDCSpaceResols.h.

    {
      const SigmaID sigmaID = getSigmaID(iCLayer, iLR, iAlpha, iTheta);
      //  std::cout <<"sigmaID in getSigmaParams= " << sigmaID << std::endl;
      return getSigmaParams(sigmaID);
    }

◆ getThetaBin()

const array3 & getThetaBin ( unsigned short i ) const

inline

Get i-th theta-angle bin info.

(rad)

Definition at line 186 of file CDCSpaceResols.h.

    {
      return m_thetaBins[i];
    }

◆ getThetaPoint()

float getThetaPoint ( unsigned short i ) const

inline

Get i-th theta-angle point (rad)

Definition at line 194 of file CDCSpaceResols.h.

    {
      return m_thetaBins[i][2];
    }

◆ outputToFile()

void outputToFile ( std::string fileName ) const

inline

Output the contents in text file format.

Definition at line 351 of file CDCSpaceResols.h.

    {
      std::ofstream fout(fileName);
 
      if (fout.bad()) {
        B2ERROR("Specified output file could not be opened!");
      } else {
        const double deg = 180. / M_PI;
 
        unsigned short nAlphaBins = m_alphaBins.size();
        fout << nAlphaBins << std::endl;
 
        for (unsigned short i = 0; i < nAlphaBins; ++i) {
          fout << deg* m_alphaBins[i][0] << "  " << deg* m_alphaBins[i][1] << "  " << deg* m_alphaBins[i][2] << std::endl;
        }
 
        unsigned short nThetaBins = m_thetaBins.size();
        fout << nThetaBins << std::endl;
 
        for (unsigned short i = 0; i < nThetaBins; ++i) {
          fout << deg* m_thetaBins[i][0] << "  " << deg* m_thetaBins[i][1] << "  " << deg* m_thetaBins[i][2] << std::endl;
        }
 
        fout << m_sigmaParamMode << "  " << m_nSigmaParams << std::endl;
 
        fout << m_maxSpaceResol << std::endl;
 
        for (unsigned short iT = 0; iT < nThetaBins; ++iT) {
          for (unsigned short iA = 0; iA < nAlphaBins; ++iA) {
            for (unsigned short iCL = 0; iCL < c_nSLayers; ++iCL) {
              for (unsigned short iLR = 0; iLR < 2; ++iLR) {
                unsigned short iLRp = abs(iLR - 1);
                fout << std::setw(2) << std::right << std::fixed << iCL << "  " << std::setw(5) << std::setprecision(
                       1) << deg* m_thetaBins[iT][2] << "  " << std::setw(5) << std::right << deg* m_alphaBins[iA][2] << "  " << std::setw(1) << iLRp;
                const std::vector<float> params = getSigmaParams(iCL, iLRp, iA, iT);
                for (unsigned short i = 0; i < m_nSigmaParams; ++i) {
                  fout << "  " << std::setw(15) << std::scientific << std::setprecision(8) << params[i];
                }
                fout << std::endl;
              }
            }
          }
        }
        fout.close();
      }
    }

◆ setAlphaBin()

void setAlphaBin ( const array3 & alpha )

inline

Set alpha-angle bin (rad)

Definition at line 50 of file CDCSpaceResols.h.

    {
      if (m_alphaBins.size() <= c_maxNAlphaBins) {
        m_alphaBins.push_back(alpha);
        sort(m_alphaBins.begin(), m_alphaBins.end(), comp);
      } else {
        //  std::cout<< m_alphaBins.size() <<" "<< c_maxNAlphaBins <<std::endl;
        B2FATAL("The no. of alpha bins > limit !");
      }
    }

◆ setMaxSpaceResol()

void setMaxSpaceResol ( float sigma )

inline

Set max.

allowed space resolution

Definition at line 93 of file CDCSpaceResols.h.

    {
      m_maxSpaceResol = sigma;
    }

◆ setSigmaParamMode()

void setSigmaParamMode ( unsigned short mode )

inline

Set sigma parameterization mode.

Definition at line 85 of file CDCSpaceResols.h.

    {
      m_sigmaParamMode = mode;
    }

◆ setSigmaParams() [1/2]

void setSigmaParams	(	const SigmaID	sigmaID,
		const std::vector< float > &	params
	)

inline

Set sigma parameters for the specified id.

(=bin)

Definition at line 101 of file CDCSpaceResols.h.

    {
      unsigned short nSigmaPars = params.size();
 
      if (nSigmaPars <= c_maxNSigmaParams) {
        m_nSigmaParams = nSigmaPars;
        m_sigmas.insert(std::pair<SigmaID, std::vector<float>>(sigmaID, params));
        //  std::cout <<"sigmaID in setSigmaParams= " << sigmaID << std::endl;
      } else {
        B2FATAL("The no. of sigma params. > limit !");
      }
    }

◆ setSigmaParams() [2/2]

void setSigmaParams	(	unsigned short	iCLayer,
		unsigned short	iLR,
		unsigned short	iAlpha,
		unsigned short	iTheta,
		const std::vector< float > &	params
	)

inline

Set sigma parameters for the specified id.

(=bin)

Definition at line 117 of file CDCSpaceResols.h.

    {
      const SigmaID sigmaID = getSigmaID(iCLayer, iLR, iAlpha, iTheta);
      setSigmaParams(sigmaID, params);
    }

◆ setThetaBin()

void setThetaBin ( const array3 & theta )

inline

Set theta-angle bin (rad)

Definition at line 64 of file CDCSpaceResols.h.

    {
      if (m_thetaBins.size() <= c_maxNThetaBins) {
        m_thetaBins.push_back(theta);
        sort(m_thetaBins.begin(), m_thetaBins.end(), comp);
      } else {
        B2FATAL("The no. of theta bins > limit !");
      }
    }

Member Data Documentation

◆ m_alphaBins

std::vector<array3> m_alphaBins

private

alpha bins for sigma (rad)

Definition at line 402 of file CDCSpaceResols.h.

◆ m_maxSpaceResol

float m_maxSpaceResol = 0.07

private

max.

space resol (cm) allowed

Definition at line 401 of file CDCSpaceResols.h.

◆ m_nSigmaParams

unsigned short m_nSigmaParams

private

no.

of sigma parameters per bin

Definition at line 400 of file CDCSpaceResols.h.

◆ m_sigmaParamMode

unsigned short m_sigmaParamMode

private

Mode for sigma parameterization.

Definition at line 399 of file CDCSpaceResols.h.

◆ m_sigmas

std::map<SigmaID, std::vector<float> > m_sigmas

private

Sigma coefficients for each layer, Left/Right, entrance angle and polar angle.

Definition at line 405 of file CDCSpaceResols.h.

◆ m_thetaBins

std::vector<array3> m_thetaBins

private

theta bins for sigma (rad)

Definition at line 403 of file CDCSpaceResols.h.

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

cdc/dbobjects/include/CDCSpaceResols.h

Public Types

Public Member Functions

Static Public Member Functions

Private Types

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ array3

◆ SigmaID

Member Enumeration Documentation

◆ anonymous enum

Constructor & Destructor Documentation

◆ CDCSpaceResols()

Member Function Documentation

◆ addSigmaParams() [1/2]

◆ addSigmaParams() [2/2]

◆ comp()

◆ dump()

◆ getAlphaBin()

◆ getAlphaPoint()

◆ getMaxSpaceResol()

◆ getNoOfAlphaBins()

◆ getNoOfThetaBins()

◆ getSigmaID() [1/2]

◆ getSigmaID() [2/2]

◆ getSigmaParamMode()

◆ getSigmaParams() [1/2]

◆ getSigmaParams() [2/2]

◆ getThetaBin()

◆ getThetaPoint()

◆ outputToFile()

◆ setAlphaBin()

◆ setMaxSpaceResol()

◆ setSigmaParamMode()

◆ setSigmaParams() [1/2]

◆ setSigmaParams() [2/2]

◆ setThetaBin()

Member Data Documentation

◆ m_alphaBins

◆ m_maxSpaceResol

◆ m_nSigmaParams

◆ m_sigmaParamMode

◆ m_sigmas

◆ m_thetaBins