Database object for xt-relations. More...

#include <CDCXtRelations.h>

Inheritance diagram for CDCXtRelations:

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

Public Member Functions
	CDCXtRelations ()
	Default constructor.

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

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

void	setXtParamMode (unsigned short mode)
	Set xt parameterization mode.

void	setXtParams (const XtID xtID, const std::vector< float > &params)
	Set xt parameters for the specified id.

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

void	addXTParams (const XtID xtID, const std::vector< float > &delta)
	Update xt parameters for the specified id.

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

void	replaceXTParams (const XtID xtID, const std::vector< float > &param)
	Replace xt parameters for the specified id.

void	replaceXTParams (unsigned short iCLayer, unsigned short iLR, unsigned short iAlpha, unsigned short iTheta, const std::vector< float > &param)
	Replace xt 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	getXtParamMode () const
	Get xt parameterization mode.

XtID	getXtID (unsigned short iCLayer, unsigned short iLR, unsigned short iAlpha, unsigned short iTheta) const
	Get xt id.

XtID	getXtID (unsigned short iCLayer, unsigned short iLR, float alpha, float theta) const
	Get xt id.

const std::vector< float > &	getXtParams (const XtID xtID) const
	Get xt parameters for the specified id.

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

void	dump () const
	Print all contents.

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

double	getGlobalParam (unsigned short xtId, unsigned short i) const
	Get global parameter for i-th component of the specified xtId.

void	setGlobalParam (double value, unsigned short xtId, unsigned short i)
	Set global parameter for i-th component of the specified xtId.

std::vector< std::pair< unsigned short, unsigned short > >	listGlobalParams () const
	list stored global parameters TODO FIXME CDC not ready

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

static unsigned short	getGlobalUniqueID ()
	Get global unique id.

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

typedef unsigned short	XtID
	id.

Private Member Functions
	ClassDef (CDCXtRelations, 2)
	ClassDef.

Private Attributes
unsigned short	m_xtParamMode
	Mode for xt parameterization.

unsigned short	m_nXtParams
	no.

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

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

std::map< XtID, std::vector< float > >	m_xts
	XT-relation coefficients for each layer, Left/Right, entrance angle and polar angle.

Detailed Description

Database object for xt-relations.

Definition at line 28 of file CDCXtRelations.h.

Member Typedef Documentation

◆ array3

typedef std::array<float, 3> array3

private

angle bin info.

Definition at line 30 of file CDCXtRelations.h.

◆ XtID

typedef unsigned short XtID

private

id.

for xt

Definition at line 31 of file CDCXtRelations.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_maxNXtParams	max. no. of params. for xt

Definition at line 37 of file CDCXtRelations.h.

         {c_nSLayers      = 56, 
          c_maxNAlphaBins = 18, 
          c_maxNThetaBins =  7, 
          c_maxNXtParams  =  8  
         };

Constructor & Destructor Documentation

◆ CDCXtRelations()

CDCXtRelations ( )

inline

Default constructor.

Definition at line 46 of file CDCXtRelations.h.

46 : m_xtParamMode(0), m_nXtParams(0)

47 {}

Belle2::CDCXtRelations::m_xtParamMode

unsigned short m_xtParamMode

Mode for xt parameterization.

Definition: CDCXtRelations.h:446

Belle2::CDCXtRelations::m_nXtParams

unsigned short m_nXtParams

no.

Definition: CDCXtRelations.h:447

Member Function Documentation

◆ addXTParams() [1/2]

void addXTParams	(	const XtID	xtID,
		const std::vector< float > &	delta
	)

inline

Update xt parameters for the specified id.

(=bin)

Definition at line 121 of file CDCXtRelations.h.

    {
      std::map<XtID, std::vector<float>>::iterator it = m_xts.find(xtID);
 
      if (it != m_xts.end()) {
        for (unsigned short i = 0; i < m_nXtParams; ++i) {
          (it->second)[i] += delta[i];
        }
      } else {
        B2FATAL("Specified xt params not found in addXTParams !");
      }
    }

◆ addXTParams() [2/2]

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

inline

Update xt parameters for the specified id.

(=bin)

Definition at line 137 of file CDCXtRelations.h.

    {
      const XtID xtID = getXtID(iCLayer, iLR, iAlpha, iTheta);
      addXTParams(xtID, delta);
    }

◆ comp()

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

inlinestatic

Static function for sorting.

Definition at line 79 of file CDCXtRelations.h.

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

◆ dump()

void dump ( ) const

inline

Print all contents.

Definition at line 319 of file CDCXtRelations.h.

    {
      std::cout << " " << std::endl;
      std::cout << "Contents of xt 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 xt" << 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 = getXtParams(iCL, iLRp, iA, iT);
              for (unsigned short i = 0; i < m_nXtParams; ++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 190 of file CDCXtRelations.h.

    {
      return m_alphaBins[i];
    }

◆ getAlphaPoint()

float getAlphaPoint ( unsigned short i ) const

inline

Get i-th alpha-angle point (rad)

Definition at line 198 of file CDCXtRelations.h.

    {
      return m_alphaBins[i][2];
    }

◆ getGlobalParam()

double getGlobalParam	(	unsigned short	xtId,
		unsigned short	i
	)		const

inline

Get global parameter for i-th component of the specified xtId.

Definition at line 414 of file CDCXtRelations.h.

    {
      return getXtParams(xtId).at(i);
    }

◆ getGlobalUniqueID()

static unsigned short getGlobalUniqueID ( )

inlinestatic

Get global unique id.

Definition at line 412 of file CDCXtRelations.h.

412{return 29;}

◆ getNoOfAlphaBins()

unsigned short getNoOfAlphaBins ( ) const

inline

Get no.

of alpha-angle bin

Definition at line 174 of file CDCXtRelations.h.

    {
      return m_alphaBins.size();
    }

◆ getNoOfThetaBins()

unsigned short getNoOfThetaBins ( ) const

inline

Get no.

of theta-angle bin

Definition at line 182 of file CDCXtRelations.h.

    {
      return m_thetaBins.size();
    }

◆ getThetaBin()

const array3 & getThetaBin ( unsigned short i ) const

inline

Get i-th theta-angle bin info.

(rad)

Definition at line 206 of file CDCXtRelations.h.

    {
      return m_thetaBins[i];
    }

◆ getThetaPoint()

float getThetaPoint ( unsigned short i ) const

inline

Get i-th theta-angle point (rad)

Definition at line 214 of file CDCXtRelations.h.

    {
      return m_thetaBins[i][2];
    }

◆ getXtID() [1/2]

XtID getXtID	(	unsigned short	iCLayer,
		unsigned short	iLR,
		float	alpha,
		float	theta
	)		const

inline

Get xt id.

(=bin id.)

Definition at line 243 of file CDCXtRelations.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 ! " << alpha);
 
      return getXtID(iCLayer, iLR, iAlpha, iTheta);
    }

◆ getXtID() [2/2]

XtID getXtID	(	unsigned short	iCLayer,
		unsigned short	iLR,
		unsigned short	iAlpha,
		unsigned short	iTheta
	)		const

inline

Get xt 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 234 of file CDCXtRelations.h.

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

◆ getXtParamMode()

unsigned short getXtParamMode ( ) const

inline

Get xt parameterization mode.

Definition at line 222 of file CDCXtRelations.h.

    {
      return m_xtParamMode;
    }

◆ getXtParams() [1/2]

const std::vector< float > & getXtParams ( const XtID xtID ) const

inline

Get xt parameters for the specified id.

(=bin)

Definition at line 295 of file CDCXtRelations.h.

    {
      std::map<XtID, std::vector<float>>::const_iterator it = m_xts.find(xtID);
      if (it != m_xts.end()) {
        return it->second;
      } else {
        B2FATAL("Specified xt params. not found in getXtParams !");
      }
    }

◆ getXtParams() [2/2]

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

inline

Get xt parameters for the specified id.

(=bin)

Definition at line 308 of file CDCXtRelations.h.

    {
      const XtID xtID = getXtID(iCLayer, iLR, iAlpha, iTheta);
      //  std::cout <<"xtID in getXtParams= " << xtID << std::endl;
      return getXtParams(xtID);
    }

◆ listGlobalParams()

std::vector< std::pair< unsigned short, unsigned short > > listGlobalParams ( ) const

inline

list stored global parameters TODO FIXME CDC not ready

Definition at line 435 of file CDCXtRelations.h.

    {
      std::vector<std::pair<unsigned short, unsigned short>> result;
      for (auto ixt : m_xts) {
        for (int i = 0; i < 8; ++i) {
          result.push_back({ixt.first, i});
        }
      }
      return result;
    }

◆ outputToFile()

void outputToFile ( std::string fileName ) const

inline

Output the contents in test file format.

Definition at line 363 of file CDCXtRelations.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_xtParamMode << "  " << m_nXtParams << std::endl;
 
        signed short phiAngle = 0.;
        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) << std::setprecision(1) << phiAngle << "  " << std::setw(1) << iLRp;
                const std::vector<float> params = getXtParams(iCL, iLRp, iA, iT);
                for (unsigned short i = 0; i < m_nXtParams; ++i) {
                  fout << "  " << std::setw(15) << std::scientific << std::setprecision(8) << params[i];
                }
                fout << std::endl;
              }
            }
          }
        }
        fout.close();
      }
    }

◆ replaceXTParams() [1/2]

void replaceXTParams	(	const XtID	xtID,
		const std::vector< float > &	param
	)

inline

Replace xt parameters for the specified id.

(=bin)

Definition at line 147 of file CDCXtRelations.h.

    {
      std::map<XtID, std::vector<float>>::iterator it = m_xts.find(xtID);
 
      if (it != m_xts.end()) {
        for (unsigned short i = 0; i < m_nXtParams; ++i) {
          (it->second)[i] = param[i];
        }
      } else {
        B2FATAL("Specified xt params not found in replaceXTParams !");
      }
    }

◆ replaceXTParams() [2/2]

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

inline

Replace xt parameters for the specified id.

(=bin)

Definition at line 163 of file CDCXtRelations.h.

    {
      const XtID xtID = getXtID(iCLayer, iLR, iAlpha, iTheta);
      replaceXTParams(xtID, param);
    }

◆ setAlphaBin()

void setAlphaBin ( const array3 & alpha )

inline

Set alpha-angle bin (rad)

Definition at line 52 of file CDCXtRelations.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 !");
      }
    }

◆ setGlobalParam()

void setGlobalParam	(	double	value,
		unsigned short	xtId,
		unsigned short	i
	)

inline

Set global parameter for i-th component of the specified xtId.

Definition at line 419 of file CDCXtRelations.h.

    {
 
      std::map<XtID, std::vector<float>>::const_iterator it = m_xts.find(xtId);
      if (it != m_xts.end()) {
        std::vector<float> allParams =  it->second;
        allParams.at(i) = value;
        setXtParams(xtId, allParams);
      } else {
        B2INFO("Specified xt params. not found in getXtParams.");
        std::vector<float> allParams {0., 0., 0., 0., 0., 0., 0., 0.};
        allParams.at(i) = value;
        setXtParams(xtId, allParams);
      }
    }

◆ setThetaBin()

void setThetaBin ( const array3 & theta )

inline

Set theta-angle bin (rad)

Definition at line 66 of file CDCXtRelations.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 !");
      }
    }

◆ setXtParamMode()

void setXtParamMode ( unsigned short mode )

inline

Set xt parameterization mode.

Definition at line 87 of file CDCXtRelations.h.

    {
      m_xtParamMode = mode;
    }

◆ setXtParams() [1/2]

void setXtParams	(	const XtID	xtID,
		const std::vector< float > &	params
	)

inline

Set xt parameters for the specified id.

(=bin)

Definition at line 95 of file CDCXtRelations.h.

    {
      unsigned short nXtParams = params.size();
 
      if (nXtParams <= c_maxNXtParams) {
        m_nXtParams = nXtParams;
        m_xts.insert(std::pair<XtID, std::vector<float>>(xtID, params));
        //  std::cout <<"xtID in setXtParams= " << xtID << std::endl;
      } else {
        B2FATAL("The no. of xt params. > limit !");
      }
    }

◆ setXtParams() [2/2]

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

inline

Set xt parameters for the specified id.

(=bin)

Definition at line 111 of file CDCXtRelations.h.

    {
      const XtID xtID = getXtID(iCLayer, iLR, iAlpha, iTheta);
      setXtParams(xtID, params);
    }

Member Data Documentation

◆ m_alphaBins

std::vector<array3> m_alphaBins

private

alpha bins for xt (rad)

Definition at line 448 of file CDCXtRelations.h.

◆ m_nXtParams

unsigned short m_nXtParams

private

no.

of xt parameters per bin

Definition at line 447 of file CDCXtRelations.h.

◆ m_thetaBins

std::vector<array3> m_thetaBins

private

theta bins for xt (rad)

Definition at line 449 of file CDCXtRelations.h.

◆ m_xtParamMode

unsigned short m_xtParamMode

private

Mode for xt parameterization.

Definition at line 446 of file CDCXtRelations.h.

◆ m_xts

std::map<XtID, std::vector<float> > m_xts

private

XT-relation coefficients for each layer, Left/Right, entrance angle and polar angle.

Definition at line 451 of file CDCXtRelations.h.

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

cdc/dbobjects/include/CDCXtRelations.h

Public Types

Public Member Functions

Static Public Member Functions

Private Types

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ array3

◆ XtID

Member Enumeration Documentation

◆ anonymous enum

Constructor & Destructor Documentation

◆ CDCXtRelations()

Member Function Documentation

◆ addXTParams() [1/2]

◆ addXTParams() [2/2]

◆ comp()

◆ dump()

◆ getAlphaBin()

◆ getAlphaPoint()

◆ getGlobalParam()

◆ getGlobalUniqueID()

◆ getNoOfAlphaBins()

◆ getNoOfThetaBins()

◆ getThetaBin()

◆ getThetaPoint()

◆ getXtID() [1/2]

◆ getXtID() [2/2]

◆ getXtParamMode()

◆ getXtParams() [1/2]

◆ getXtParams() [2/2]

◆ listGlobalParams()

◆ outputToFile()

◆ replaceXTParams() [1/2]

◆ replaceXTParams() [2/2]

◆ setAlphaBin()

◆ setGlobalParam()

◆ setThetaBin()

◆ setXtParamMode()

◆ setXtParams() [1/2]

◆ setXtParams() [2/2]

Member Data Documentation

◆ m_alphaBins

◆ m_nXtParams

◆ m_thetaBins

◆ m_xtParamMode

◆ m_xts