CDCDedxInjectionTime Class Reference

dE/dx injection time calibration constants More...

#include <CDCDedxInjectionTime.h>

Inheritance diagram for CDCDedxInjectionTime:

Collaboration diagram for CDCDedxInjectionTime:

Public Member Functions
	CDCDedxInjectionTime ()
	Default constructor.
	CDCDedxInjectionTime (const std::vector< std::vector< double >> &vinjcalib)
	Constructor.
	~CDCDedxInjectionTime ()
	Destructor.
void	printCorrection (std::string svar, std::string sfx) const
	Return dE/dx mean or norm-reso value for the given time and ring. More...
std::vector< std::vector< double > >	getConstVector () const
	Return vector of all constant vector of payload.
std::vector< double >	getTimeVector (const unsigned int ring) const
	Return time vector. More...
std::vector< double >	getMeanVector (const int ring) const
	Return dedx mean vector. More...
std::vector< double >	getResoVector (const unsigned int ring) const
	Return dedx reso vector. More...
unsigned int	getTimeBin (const std::vector< unsigned int > &array, unsigned int value) const
	Return time bin for the given time array. More...
double	getCorrection (std::string svar, unsigned int ring, unsigned int time) const
	Return dE/dx mean or norm-reso value for the given time and ring. More...

Private Member Functions
	ClassDef (CDCDedxInjectionTime, 1)
	ClassDef.

Private Attributes
std::vector< std::vector< double > >	m_injectionvar
	CDC dE/dx injection time payloads for dEdx mean and reso. More...

Detailed Description

dE/dx injection time calibration constants

Definition at line 29 of file CDCDedxInjectionTime.h.

Member Function Documentation

getCorrection()

double getCorrection	(	std::string	svar,
		unsigned int	ring,
		unsigned int	time
	)		const

Return dE/dx mean or norm-reso value for the given time and ring.

Parameters

svar	is option for mean and reso calibration factor
ring	is injection ring number (0/1 for LER/HER)
time	is injection time (large 0-20sec range)

Definition at line 80 of file CDCDedxInjectionTime.cc.

{
  if (svar != "mean" && svar != "reso") {
    B2ERROR("wrong var input, choose mean or reso");
    return 1.0;
  }
 
  if (std::isnan(ring) || std::isnan(time))return 1.0;
 
  if (ring > 1) {
    B2ERROR("wrong ring input, choose 0 or 1");
    return 1.0;
  }
 
  unsigned int iv = ring * 3 + 1 ;
  if (svar == "reso") iv = ring * 3 + 2 ;
 
  unsigned int sizev = m_injectionvar[iv].size(); //mean or reso
  unsigned int sizet = m_injectionvar[ring * 3].size(); //time
  if (sizet == 0 || sizev == 0) {
    B2ERROR("calibration vectors are empty");
    return 1.0;
  }
 
  std::vector<unsigned int> tedges(sizet); //time edges array
  std::copy(m_injectionvar[ring * 3].begin(), m_injectionvar[ring * 3].end(), tedges.begin());
 
  if (time >= 5e6)time = 5e6 - 10;
 
  unsigned int it = getTimeBin(tedges, time);
  double center = 0.5 * (m_injectionvar[ring * 3].at(it) + m_injectionvar[ring * 3].at(it + 1));
 
  //no corr before veto bin (usually one or two starting bin)
  //intrapolation for entire range except
  //--extrapolation (for first half and last half of intended bin)
  int thisbin = it, nextbin = it;
  if (center != time && it > 0) {
 
    if (time < center) {
      thisbin = it - 1;
    } else {
      if (it < sizet - 2)nextbin = it + 1;
      else thisbin = it - 1;
    }
 
    if (it <= 2) {
      double diff = m_injectionvar[iv].at(2) - m_injectionvar[iv].at(1) ;
      if (diff < -0.015) { //difference above 1.0%
        thisbin = it;
        if (it == 1) nextbin = it;
        else nextbin = it + 1;
      } else {
        if (it == 1) {
          thisbin = it;
          nextbin = it + 1;
        }
      }
    }
  }
 
  double thisdedx = m_injectionvar[iv].at(thisbin);
  double nextdedx = m_injectionvar[iv].at(nextbin);
 
  double thistime = 0.5 * (m_injectionvar[ring * 3].at(thisbin) + m_injectionvar[ring * 3].at(thisbin + 1));
  double nexttime = 0.5 * (m_injectionvar[ring * 3].at(nextbin) + m_injectionvar[ring * 3].at(nextbin + 1));
 
  double newdedx = m_injectionvar[iv].at(it);
  if (thisbin != nextbin)
    newdedx = thisdedx + ((nextdedx - thisdedx) / (nexttime - thistime)) * (time - thistime);
 
  return newdedx;
}

getMeanVector()

std::vector<double> getMeanVector ( const int  ring ) const

inline

Return dedx mean vector.

Parameters

ring	is injection ring number (0/1 for LER/HER)

Definition at line 76 of file CDCDedxInjectionTime.h.

    {
      if (ring < 0 || ring > 1) B2ERROR("wrong index for injection ring ");
      return m_injectionvar[ring * 3 + 1];
    };

getResoVector()

std::vector<double> getResoVector ( const unsigned int  ring ) const

inline

Return dedx reso vector.

Parameters

ring	is injection ring number (0/1 for LER/HER)

Definition at line 85 of file CDCDedxInjectionTime.h.

getTimeBin()

unsigned int getTimeBin ( const std::vector< unsigned int > &  array, unsigned int  value  ) const

inline

Return time bin for the given time array.

Parameters

array	of time
value	of input time

Definition at line 95 of file CDCDedxInjectionTime.h.

getTimeVector()

std::vector<double> getTimeVector ( const unsigned int  ring ) const

inline

Return time vector.

Parameters

ring	is injection ring number (0/1 for LER/HER)

Definition at line 67 of file CDCDedxInjectionTime.h.

printCorrection()

void printCorrection	(	std::string	svar,
		std::string	sfx
	)		const

Return dE/dx mean or norm-reso value for the given time and ring.

Parameters

svar	is option for printing mean and reso calibration
sfx	to add suffix in file save

Definition at line 13 of file CDCDedxInjectionTime.cc.

Member Data Documentation

m_injectionvar

std::vector<std::vector<double> > m_injectionvar

private

CDC dE/dx injection time payloads for dEdx mean and reso.

different for LER and HER vector to store payload values

Definition at line 122 of file CDCDedxInjectionTime.h.

---

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

• reconstruction/dbobjects/include/CDCDedxInjectionTime.h
• reconstruction/dbobjects/src/CDCDedxInjectionTime.cc