TileCalibDrawerFlt Class Reference

Generic class for storing a number of floats (Flt) for each channel or ADC. More...

#include <TileCalibDrawerFlt.h>

Inheritance diagram for TileCalibDrawerFlt:

Collaboration diagram for TileCalibDrawerFlt:

Public Types
typedef std::vector< std::vector< float > >	DefType
	Object to hold default data used for initialization.

Public Member Functions
virtual	~TileCalibDrawerFlt ()
	Dtor.
virtual uint16_t	getType () const
	Returns TileCalibType::FLT.
float	getCalib (unsigned int channel, unsigned int adc, float energy, bool invert=false) const
	Returns the calibrated energy for a given channel/ADC and input energy.
bool	getYDY (unsigned int channel, unsigned int adc, float x, float &y, float &dy) const
	Returns y and derivative dy for a given x for function blob of type 200.
float	getY (unsigned int channel, unsigned int adc, float x) const
	Returns y for a given x for function blob of type 200.
float	getDY (unsigned int channel, unsigned int adc, float x) const
	Returns y for a given x for function blob of type 200.
float	getData (unsigned int channel, unsigned int adc, unsigned int idx) const
	Returns a single T belonging to a channel/ADC.
void	init (const DefType &def, uint16_t nChans, uint16_t objVers, const std::string &author="", const std::string &comment="", uint64_t timeStamp=0)
	Initializing function.
void	setData (unsigned int channel, unsigned int adc, unsigned int idx, float data)
	Sets a single T belonging to a channel/ADC.
virtual void	dump () const
	Prints out the content of the blob to std::out.
virtual void	dump (std::ostream &stm) const
	Prints the BLOB header summary information.
const void *	getAddress (unsigned int iEle) const
	Returns start address of iEle-th basic unit.
void *	getAddress (unsigned int iEle)
virtual void	clone (const TileCalibDrawerBase &other)
	Initialzes Blob with content of other Blob.
long	getBlobSize () const
	Returns the BLOB size in units of bytes.
long	getBlobSize32 () const
	Returns the BLOB size in units of uint32_t.
uint16_t	getObjType () const
	Returns the BLOB object type.
uint16_t	getObjVersion () const
	Returns the BLOB object version.
uint32_t	getObjSizeUint32 () const
	Returns the size of a data object in units of uint32_t.
uint32_t	getObjSizeByte () const
	Returns the size of a data object in units of bytes.
uint32_t	getNObjs () const
	Returns the number of data objects stored int the BLOB.
uint16_t	getNChans () const
	Returns the number of channels stored in the BLOB.
uint16_t	getNGains () const
	Returns the number of gains stored for each channel.
uint32_t	getCommentSizeUint32 () const
	Returns the space occupied by the comment fields in units of uint32_t.
uint32_t	getCommentSizeChar () const
	Returns the space occupied by the comment fields in units of chars.
std::string	getAuthor () const
	Returns the comment author.
std::string	getComment () const
	Returns the actual comment.
uint64_t	getTimeStamp () const
	Returns the unix timestamp of the comment (seconds since 1.1.1970)
std::string	getDate () const
	Returns the date of the comment as string (derived from timestamp)
std::string	getFullComment () const
	Returns a formated string build from all comment fields.

Static Public Member Functions
static const TileCalibDrawerFlt *	getInstance (const coral::Blob &blob)
	Returns a pointer to a const TileCalibDrawerFlt.
static TileCalibDrawerFlt *	getInstance (coral::Blob &blob, const DefType &def, uint16_t nChans, uint16_t objVers, const std::string &author="", const std::string &comment="", uint64_t timeStamp=0)
	Returns a pointer to a non-const TileCalibDrawerFlt.

Static Public Attributes
static const unsigned int	m_hdrSize32
	The header size in units of uint32_t.

Protected Member Functions
	TileCalibDrawerFlt (const coral::Blob &blob)
	Ctor (const).
	TileCalibDrawerFlt (coral::Blob &blob)
	Ctor (non-const).
const float *	getAddress (unsigned int channel, unsigned int adc) const
	Returns a pointer to the first value for the specified channel & ADC.
uint32_t	createBlob (uint16_t objType, uint16_t objVersion, uint32_t objSizeUint32, uint32_t nObjs, uint16_t nChans, uint16_t nGains, const std::string &author="", const std::string &comment="", uint64_t timeStamp=0)
	(re-)creation of the referenced BLOB object.
void	dumpHeader (std::ostream &stm) const
	Prints the BLOB header summary information.

Private Attributes
coral::Blob *	m_blob_nc
	Non-const reference to the BLOB.
const coral::Blob *	m_blob
	Const reference to the BLOB (always there)
const uint32_t *	m_blobStart32
	Cache blob starting address as uint_32t*.
const uint16_t *	m_blobStart16
	Cache blob starting address as uint_16t*.
uint64_t	m_blobSize32
	Cache blob size in units of uint32_t.
bool	m_isBlobOwner
	Is this TileCalibDrawer owner of the BLOB.

Detailed Description

Generic class for storing a number of floats (Flt) for each channel or ADC.

Author

Nils Gollub nils..nosp@m.goll.nosp@m.ub@ce.nosp@m.rn.c.nosp@m.h

This class implements the infrastructure to store a (variable) number of floats for each channel or ADC. The number of floats to store is encoded in the objSize header variable. The object version header variable is used to encode different functional forms on the stored data, accessed through the getCalib function.

Definition at line 26 of file TileCalibDrawerFlt.h.

Member Typedef Documentation

DefType

typedef std::vector<std::vector<float> > TileCalibDrawerDat< float >::DefType

inherited

Object to hold default data used for initialization.

The length of the first and second vector determines respectively the number of gains and values per gain to be stored. The indexing is thus DefType[iGain][iValue].

Definition at line 37 of file TileCalibDrawerDat.h.

Constructor & Destructor Documentation

~TileCalibDrawerFlt()

virtual TileCalibDrawerFlt::~TileCalibDrawerFlt ( )

inlinevirtual

Dtor.

Definition at line 31 of file TileCalibDrawerFlt.h.

{}

TileCalibDrawerFlt() [1/2]

TileCalibDrawerFlt::TileCalibDrawerFlt ( const coral::Blob & blob )

protected

Ctor (const).

Definition at line 40 of file TileCalibDrawerFlt.cxx.

                                                            : 
  TileCalibDrawerDat<float>(blob)
{
  if(getBlobSize()){
    if(getObjType() != TileCalibDrawerFlt::getType()){
      throw TileCalib::TypeConflict("TileCalibDrawerFlt::Ctor",getObjType(),TileCalibDrawerFlt::getType());
    }
  }
}

TileCalibDrawerFlt() [2/2]

TileCalibDrawerFlt::TileCalibDrawerFlt ( coral::Blob & blob )

protected

Ctor (non-const).

Definition at line 52 of file TileCalibDrawerFlt.cxx.

                                                      : 
  TileCalibDrawerDat<float>(blob)
{
  if(getBlobSize()){
    if(getObjType() != TileCalibDrawerFlt::getType()){
      throw TileCalib::TypeConflict("TileCalibDrawerFlt::Ctor",getObjType(),TileCalibDrawerFlt::getType());
    }
  }
}

Member Function Documentation

clone()

void TileCalibDrawerBase::clone ( const TileCalibDrawerBase & other )

virtualinherited

Initialzes Blob with content of other Blob.

Parameters

other

TileCalibDrawer used for initialization

Definition at line 77 of file TileCalibDrawerBase.cxx.

{
  //=== copy content of other blob
  *m_blob_nc = *other.m_blob;
  //=== and reset cached attributes
  m_blobStart32 = static_cast<const uint32_t*>(m_blob->startingAddress());
  m_blobStart16 = static_cast<const uint16_t*>(m_blob->startingAddress());
  m_blobSize32  = m_blob->size()/sizeof(uint32_t);
}

createBlob()

uint32_t TileCalibDrawerBase::createBlob ( uint16_t objType, uint16_t objVersion, uint32_t objSizeUint32, uint32_t nObjs, uint16_t nChans, uint16_t nGains, const std::string & author = "", const std::string & comment = "", uint64_t timeStamp = 0 )

protectedinherited

(re-)creation of the referenced BLOB object.

Parameters

objType	Object type
objVersion	Object version
objSizeUint32	Size of a data object (in uint32_t)
nObjs	Total number of data objects
nChans	number of stored channels
nGains	number of stored gains
author	comment author
comment	actual comment
timeStamp	unix time stamp (if 0, current time is used)

Definition at line 157 of file TileCalibDrawerBase.cxx.

{
  //=== blob data length including header in bytes
  uint32_t dataSizeByte = (m_hdrSize32+objSizeUint32*nObjs) * sizeof(uint32_t);
  
  //=== calculate comment length, including two ASCII NULLs to end text fields
  uint32_t commentSizeChar(0);
  if(author.size() || comment.size() || timeStamp){
    commentSizeChar += author.size()+comment.size()+sizeof(uint64_t) + 2;
    //=== force comment length to end on 4 byte boundary
    commentSizeChar += (commentSizeChar % sizeof(uint32_t)) ? 
      (sizeof(uint32_t)-(commentSizeChar % sizeof(uint32_t))) : 0;
  }
  
  //=== create blob
  uint32_t blobSizeInBytes = dataSizeByte+commentSizeChar;
  m_blob_nc->resize(blobSizeInBytes);
  uint32_t* blobStart32 = static_cast<uint32_t*>(m_blob_nc->startingAddress());
  uint16_t* blobStart16 = static_cast<uint16_t*>(m_blob_nc->startingAddress());
  m_blobSize32  = m_blob_nc->size()/sizeof(uint32_t);
 
  m_blobStart32 = blobStart32;
  m_blobStart16 = blobStart16;
 
  //=== fill header
  blobStart16[0] = objType;
  blobStart16[1] = objVersion;
  blobStart32[1] = objSizeUint32;
  blobStart32[2] = nObjs;
  blobStart16[6] = nChans;
  blobStart16[7] = nGains;
  blobStart32[4] = commentSizeChar/sizeof(uint32_t);
  
  //==== fill comment fields
  if(commentSizeChar){
    if(!timeStamp) timeStamp = ::time(0);
    uint64_t* pTimeStamp = reinterpret_cast<uint64_t*>(blobStart32+dataSizeByte/sizeof(uint32_t));
    pTimeStamp[0] = timeStamp;
    char* pChar = reinterpret_cast<char*>(++pTimeStamp); 
    std::string::const_iterator iStr = author.begin();
    for(; iStr!=author.end(); ++iStr){ *pChar = *iStr; ++pChar; }
    *pChar = 0; 
    for(iStr=comment.begin(); iStr!=comment.end(); ++iStr){ *(++pChar) = *iStr; }
    *(++pChar) = 0;
  }
 
  return (blobSizeInBytes/sizeof(uint32_t));
}

dump() [1/2]

virtual void TileCalibDrawerBase::dump ( std::ostream & stm ) const

inlinevirtualinherited

Prints the BLOB header summary information.

Parameters

stm	The output stream to use

Reimplemented in TileCalibDrawerCmt, TileCalibDrawerDat< T >, and TileCalibDrawerOfc.

Definition at line 86 of file TileCalibDrawerBase.h.

{dumpHeader(stm);} 

dump() [2/2]

virtual void TileCalibDrawerDat< float >::dump ( ) const

inlinevirtualinherited

Prints out the content of the blob to std::out.

Reimplemented from TileCalibDrawerBase.

Definition at line 81 of file TileCalibDrawerDat.h.

{ dump(std::cout); }

dumpHeader()

void TileCalibDrawerBase::dumpHeader ( std::ostream & stm ) const

protectedinherited

Prints the BLOB header summary information.

Parameters

stm	output stream to use

Definition at line 169 of file TileCalibDrawerBase.cxx.

{
  stm << "This is a " << TileCalibType::getClassName(getObjType()) << std::endl;
  stm << "ObjType        : " << getObjType()     << std::endl;
  stm << "ObjVersion     : " << getObjVersion()  << std::endl;
  stm << "ObjSize [bytes]: " << getObjSizeByte() << std::endl;
  stm << "NObjs          : " << getNObjs()       << std::endl;
  stm << "NChannels      : " << getNChans()      << std::endl;
  stm << "NGains         : " << getNGains()      << std::endl;
  if(!getCommentSizeUint32()){
    stm << "=== No comment available ===" << std::endl;
  }
  else{
    stm << "Author         : " << getAuthor()      << std::endl;
    stm << "Date           : " << getDate() << " ("<< getTimeStamp() << ")" << std::endl;
    stm << "Comment:       : " << getComment()     << std::endl;
  }  
}

getAddress() [1/3]

const float * TileCalibDrawerDat< float >::getAddress ( unsigned int channel, unsigned int adc ) const

protectedinherited

Returns a pointer to the first value for the specified channel & ADC.

Parameters

channel	The channel number; If channel number >= getNChans() it is reset to channel % (maximum number of channels in drawer) if channel number > (maximum number of channels in drawer) otherwise it is reset to 0 without warning (default policy)
adc	The gain index; if >= getNGains() it is reset to 0 without warning (default policy)

Definition at line 100 of file TileCalibDrawerDat.h.

{
  unsigned int idx = channel*getNGains() + adc;
  return static_cast<const T*>(TileCalibDrawerBase::getAddress(idx));
}

getAddress() [2/3]

void * TileCalibDrawerBase::getAddress ( unsigned int iEle )

inherited

getAddress() [3/3]

const void * TileCalibDrawerBase::getAddress ( unsigned int iEle ) const

inherited

Returns start address of iEle-th basic unit.

Parameters

iEle	sequential basic unit number

getAuthor()

std::string TileCalibDrawerBase::getAuthor ( ) const

inherited

Returns the comment author.

Definition at line 118 of file TileCalibDrawerBase.cxx.

{
  if(!getCommentSizeUint32()) return std::string("");
  const char* iBeg =
    reinterpret_cast<const char*>(m_blobStart32 + m_hdrSize32     +
                                  getNObjs()*getObjSizeUint32()   +
                                  sizeof(uint64_t)/sizeof(uint32_t));
  return std::string(iBeg);
}

getBlobSize()

long TileCalibDrawerBase::getBlobSize ( ) const

inlineinherited

Returns the BLOB size in units of bytes.

Definition at line 92 of file TileCalibDrawerBase.h.

{ return m_blobSize32 * sizeof(uint32_t); } 

getBlobSize32()

long TileCalibDrawerBase::getBlobSize32 ( ) const

inlineinherited

Returns the BLOB size in units of uint32_t.

Definition at line 94 of file TileCalibDrawerBase.h.

{ return m_blobSize32; } 

getCalib()

float TileCalibDrawerFlt::getCalib	(	unsigned int	channel,
		unsigned int	adc,
		float	energy,
		bool	invert = false ) const

Returns the calibrated energy for a given channel/ADC and input energy.

Parameters

channel	The channel number; If channel number >= getNChans() it is reset to channel % (maximum number of channels in drawer) if channel number > (maximum number of channels in drawer) otherwise it is reset to 0 without warning (default policy)
adc	The gain index; if >= getNGains() it is reset to 0 without warning (default policy)
energy	The energy that should be calibrated
invert	If true, the calibration is inverted

Definition at line 65 of file TileCalibDrawerFlt.cxx.

{
  //=== interpretation of data depends on the objVersion
  if(getObjVersion()==1){
 
    //=== Application of a single slope
    float slope = TileCalibDrawerFlt::getData(channel,adc,0);
    if(invert) 
      { 
    energy *= slope; 
      }
    else       
      { 
    if(slope==0.) throw TileCalib::InvalidValue("TileCalibDrawerFlt::getCalib(), SLOPE",slope);
    energy /= slope; 
      }
    return energy;
    
  }
  else if(getObjVersion()==100){
    
    //============================================================================
    //=== Look up table (LUT)
    //=== Values are stored in one x and one y column
    //=== Note: when filling the vector for the sqlite input file
    //=== first insert all x-values and then all y-values
    //=== Linear interpolation is done between two x-value entries
    //=== Under- and overflows return the input value (convenient for CIS LUT)
    //=== The return value is the LUT value!
    //============================================================================
 
    //NGO need to do something if "invert" is requested
 
    //NGO change LUT format to slope and residual? Makes it more consistent with rest
    //NGO and easier for visualization??
 
    //=== need at least 2 (x,y) pairs,
    //=== return unit function as default
    if(getObjSizeUint32()<4){
      return energy;
    }
 
    //=== Tile Default Policy has to be enforced here again, 
    //=== since getAddress() below assumes valid channel & adc
 
    if(channel >= getNChans()) {
      if (channel <= TileCalibUtils::MAX_CHAN) {
        channel = 0;
      } else {
        channel = channel % TileCalibUtils::MAX_CHAN;
        if (channel >= getNChans()) channel = 0;
      }
    }
 
    if(adc >= getNGains()) {adc = 0;}
    
    //=== determine position in table
    const float* xBeg   = getAddress(channel,adc);
    const float* xEnd   = xBeg + getObjSizeUint32()/2;
    const float* iX     = std::upper_bound(xBeg,xEnd,energy);
    unsigned int offset = iX-xBeg;
    float y(0.);
    //=== boundary values
    if(iX==xBeg || iX==xEnd){
      y = energy;
    }
    //=== linear interpolation
    else{
      float x1 = *(iX  );
      float x0 = *(iX-1);
      float y1 = *(xEnd+offset  );
      float y0 = *(xEnd+offset-1);
      y = y0+(energy-x0)*(y1-y0)/(x1-x0);  
    }
    return y;
  }
 
  //=== did not recognize the blob version 
  throw TileCalib::VersionConflict("TileCalibDrawerFlt::getCalib",getObjVersion());
}

getComment()

std::string TileCalibDrawerBase::getComment ( ) const

inherited

Returns the actual comment.

Definition at line 120 of file TileCalibDrawerBase.cxx.

{
  if(!getCommentSizeUint32()) return std::string("");
  const char* iBeg =
    reinterpret_cast<const char*>(m_blobStart32 + m_hdrSize32     +
                                  getNObjs()*getObjSizeUint32()   +
                                  sizeof(uint64_t)/sizeof(uint32_t));
  const char* iEnd = iBeg + getCommentSizeChar();
  iBeg = std::find(iBeg,iEnd,0);
  return std::string(++iBeg);
}

getCommentSizeChar()

uint32_t TileCalibDrawerBase::getCommentSizeChar ( ) const

inlineinherited

Returns the space occupied by the comment fields in units of chars.

Definition at line 112 of file TileCalibDrawerBase.h.

{return getCommentSizeUint32()*sizeof(uint32_t);}

getCommentSizeUint32()

uint32_t TileCalibDrawerBase::getCommentSizeUint32 ( ) const

inherited

Returns the space occupied by the comment fields in units of uint32_t.

getData()

float TileCalibDrawerDat< float >::getData ( unsigned int channel, unsigned int adc, unsigned int idx ) const

inherited

Returns a single T belonging to a channel/ADC.

Parameters

adc	The gain index; if >= getNGains() it is reset to 0 without warning (default policy)
idx	The index of the requested value

Definition at line 49 of file TileCalibDrawerDat.h.

{
  if(idx>=getObjSizeUint32()){
    throw TileCalib::IndexOutOfRange("TileCalibDrawerDat::getData",idx,getObjVersion());
  }
 
  //=== Tile Default Policy
 
  if(channel >= getNChans()) {
    if (channel <= TileCalibUtils::MAX_CHAN) {
      channel = 0;
    } else {
      channel = channel % TileCalibUtils::MAX_CHAN;
      if (channel >= getNChans()) channel = 0;
    }
  }
 
  if(adc >= getNGains()) {adc = 0;}
 
  return getAddress(channel,adc)[idx];
}

getDate()

std::string TileCalibDrawerBase::getDate ( ) const

inherited

Returns the date of the comment as string (derived from timestamp)

Definition at line 124 of file TileCalibDrawerBase.cxx.

{
  if(!getCommentSizeUint32()) return std::string("");
  ::time_t timeStamp = getTimeStamp();
  char buf[32];
  char* iBeg = ::ctime_r(&timeStamp, buf);
  char* iEnd = iBeg;
  while(*iEnd!='\n'){++iEnd;}
  return std::string(iBeg,iEnd-iBeg);
}

getDY()

float TileCalibDrawerFlt::getDY	(	unsigned int	channel,
		unsigned int	adc,
		float	x ) const

Returns y for a given x for function blob of type 200.

This function is inefficient and created for python use only. (Python does not like float references)

Definition at line 275 of file TileCalibDrawerFlt.cxx.

{
  float y(-9999.), dy(-9999.);
  (void)getYDY(channel,adc,x,y,dy);
  return dy;
}

getFullComment()

std::string TileCalibDrawerBase::getFullComment ( ) const

inherited

Returns a formated string build from all comment fields.

Definition at line 126 of file TileCalibDrawerBase.cxx.

{
  if(!getCommentSizeUint32()) return std::string("");
  return getAuthor()+" ("+getDate()+"): "+getComment();
}

getInstance() [1/2]

const TileCalibDrawerFlt * TileCalibDrawerFlt::getInstance ( const coral::Blob & blob )

static

Returns a pointer to a const TileCalibDrawerFlt.

Definition at line 13 of file TileCalibDrawerFlt.cxx.

{
  //=== const Blob needs to be correctly formated
  if(blob.size()<static_cast<long>(m_hdrSize32)){
    throw TileCalib::InvalidBlob("TileCalibDrawerFlt::getInstance");
  }
  return (new TileCalibDrawerFlt(blob));
}

getInstance() [2/2]

TileCalibDrawerFlt * TileCalibDrawerFlt::getInstance ( coral::Blob & blob, const DefType & def, uint16_t nChans, uint16_t objVers, const std::string & author = "", const std::string & comment = "", uint64_t timeStamp = 0 )

static

Returns a pointer to a non-const TileCalibDrawerFlt.

Parameters

blob	Reference to a coral::Blob
def	A reference to a DefType object, specifying the inital layout and values
nChans	Number of channels to create
objType	Type of the (derived class) object
objVers	Version of the (derived class) object
author	String containing the author user id
comment	String containing the actual comment
timeStamp	Unix time stamp (if 0, current time is used)

Definition at line 25 of file TileCalibDrawerFlt.cxx.

{
  TileCalibDrawerFlt* calibDrawer = new TileCalibDrawerFlt(blob);
  calibDrawer->init(def,nChans,objVers,author,comment,timeStamp);
  return calibDrawer;
}

getNChans()

uint16_t TileCalibDrawerBase::getNChans ( ) const

inherited

Returns the number of channels stored in the BLOB.

getNGains()

uint16_t TileCalibDrawerBase::getNGains ( ) const

inherited

Returns the number of gains stored for each channel.

getNObjs()

uint32_t TileCalibDrawerBase::getNObjs ( ) const

inherited

Returns the number of data objects stored int the BLOB.

getObjSizeByte()

uint32_t TileCalibDrawerBase::getObjSizeByte ( ) const

inlineinherited

Returns the size of a data object in units of bytes.

Definition at line 102 of file TileCalibDrawerBase.h.

{return getObjSizeUint32()*sizeof(uint32_t);}

getObjSizeUint32()

uint32_t TileCalibDrawerBase::getObjSizeUint32 ( ) const

inherited

Returns the size of a data object in units of uint32_t.

getObjType()

uint16_t TileCalibDrawerBase::getObjType ( ) const

inherited

Returns the BLOB object type.

getObjVersion()

uint16_t TileCalibDrawerBase::getObjVersion ( ) const

inherited

Returns the BLOB object version.

getTimeStamp()

uint64_t TileCalibDrawerBase::getTimeStamp ( ) const

inherited

Returns the unix timestamp of the comment (seconds since 1.1.1970)

getType()

virtual uint16_t TileCalibDrawerFlt::getType ( ) const

inlinevirtual

Returns TileCalibType::FLT.

Reimplemented from TileCalibDrawerBase.

Definition at line 34 of file TileCalibDrawerFlt.h.

{return TileCalibType::FLT;} 

getY()

float TileCalibDrawerFlt::getY	(	unsigned int	channel,
		unsigned int	adc,
		float	x ) const

Returns y for a given x for function blob of type 200.

This function is created as a helper function for TileCalibDrawer::getYDY and for python use. (Python does not like float references)

Definition at line 217 of file TileCalibDrawerFlt.cxx.

{
  float y(-9999.);
  //=== need to have a "function" blob version
  if(getObjVersion()!=200){
    throw TileCalib::VersionConflict("TileCalibDrawerFlt::getY",getObjVersion());
  }
    
  //=== require minimum pulse shape
  if(getObjSizeUint32()<4){
    return false;
  }
  
  //=== Tile Default Policy has to be enforced here again, 
  //=== since getAddress() below assumes valid channel & adc
  
  if(channel >= getNChans()) {
    if (channel <= TileCalibUtils::MAX_CHAN) {
      channel = 0;
    } else {
      channel = channel % TileCalibUtils::MAX_CHAN;
      if (channel >= getNChans()) channel = 0;
    }
  }
  
  if(adc >= getNGains()) {adc = 0;}
  
  //=== determine position in table
  const float* xBeg   = getAddress(channel,adc);
  const float* xEnd   = xBeg + getObjSizeUint32()/2;
  const float* iX     = std::upper_bound(xBeg,xEnd,x);
  unsigned int offset = iX-xBeg;
 
  //=== under-/overflow: return boundary y
  if(iX==xBeg){
    y  = *(xEnd+offset);
    return y;
  }
  else if(iX==xEnd){
    y  = *(xEnd+offset-1);
    return y;
  }
 
  //=== linear interpolation to determine y
  else{
    float x1 = *(iX  );
    float x0 = *(iX-1);
    float y1 = *(xEnd+offset  );
    float y0 = *(xEnd+offset-1);
    float slope = (y1-y0)/(x1-x0);
    y  = y0+(x-x0)*slope;  
  }
  return y;
}

getYDY()

bool TileCalibDrawerFlt::getYDY	(	unsigned int	channel,
		unsigned int	adc,
		float	x,
		float &	y,
		float &	dy ) const

Returns y and derivative dy for a given x for function blob of type 200.

Returns

A bool indicating whether x was within the function range (true) or outside (false)

Parameters

channel	The channel number; If channel number >= getNChans() it is reset to channel % (maximum number of channels in drawer) if channel number > (maximum number of channels in drawer) otherwise it is reset to 0 without warning (default policy)
adc	The gain index; if >= getNGains() it is reset to 0 without warning (default policy)
x	Input x value
y	Y value corresponding to x
dy	Derivative of the function at x

Definition at line 149 of file TileCalibDrawerFlt.cxx.

{
  //=== need to have a "function" blob version
  if(getObjVersion()!=200){
    throw TileCalib::VersionConflict("TileCalibDrawerFlt::getYDY",getObjVersion());
  }
    
  //=== require minimum pulse shape
  if(getObjSizeUint32()<4){
    return false;
  }
  
  //=== Tile Default Policy has to be enforced here again, 
  //=== since getAddress() below assumes valid channel & adc
 
  if(channel >= getNChans()) {
    if (channel <= TileCalibUtils::MAX_CHAN) {
      channel = 0;
    } else {
      channel = channel % TileCalibUtils::MAX_CHAN;
      if (channel >= getNChans()) channel = 0;
    }
  }
 
  if( adc >= getNGains()) {adc = 0;}
  
  //=== determine position in table
  const float* xBeg   = getAddress(channel,adc);
  const float* xEnd   = xBeg + getObjSizeUint32()/2;
  const float* iX     = std::upper_bound(xBeg,xEnd,x);
  unsigned int offset = iX-xBeg;
 
  //=== under-/overflow: return boundary y and dy==0
  //=== offset delta is to avoid float comparison troubles
  const float delta = 1.0e-3;
  if(iX==xBeg) {
    y  = *(xEnd+offset);
    dy = 0.;
    if (x < (*(xBeg) - delta)) return false;
  }
  else if(iX==xEnd) {
    y  = *(xEnd+offset-1);
    dy = 0.;
    if  (x > (*(xEnd-1) + delta)) return false;
  }
 
  //=== linear interpolation to determine y
  else{
    float x1 = *(iX  );
    float x0 = *(iX-1);
    float y1 = *(xEnd+offset  );
    float y0 = *(xEnd+offset-1);
    float slope = (y1-y0)/(x1-x0);
    y  = y0+(x-x0)*slope;  
    
    //=== determine dy from secant near [x,y]
    float dx = x1 - x0;
    float xs1 = x + 0.5 * dx;
    float xs0 = x - 0.5 * dx;
    dy = (getY(channel, adc, xs1) - getY(channel, adc, xs0)) / dx;
  }
  return true;
}

init()

void TileCalibDrawerDat< float >::init ( const DefType & def, uint16_t nChans, uint16_t objVers, const std::string & author = "", const std::string & comment = "", uint64_t timeStamp = 0 )

inherited

Initializing function.

Parameters

def	A reference to a DefType object, specifying the inital layout and values
nChans	Number of channels to create
objType	Type of the (derived class) object
objVers	Version of the (derived class) object
author	String containing the author user id
comment	String containing the actual comment
timeStamp	Unix time stamp (if 0, current time is used)

Definition at line 59 of file TileCalibDrawerDat.h.

{
  uint16_t nGains   = def.size();
  uint32_t objSize  = nGains>0 ? def.begin()->size() * sizeof(T)/sizeof(uint32_t) : 0; 
  uint32_t nObjs    = nChans*nGains;
 
  //=== ensure all gains have the same object size
  for(unsigned int iGain=0; iGain<nGains; ++iGain){
    unsigned int thisSize = def[iGain].size()*sizeof(T)/sizeof(uint32_t);
    if( thisSize != objSize ){
      throw TileCalib::SizeConflict("TileCalibDrawerDat<T>::init",thisSize,objSize);
    }
  }
  
  //=== create blob
  createBlob(getType(), objVers, objSize, nObjs, nChans, nGains, author, comment, timeStamp);
 
  //=== fill the blob with default
  for(unsigned int channel=0; channel<nChans; ++channel){
    for(unsigned int adc=0; adc<nGains; ++adc){
      setData(channel, adc, def[adc]);
    }
  }
}

setData()

void TileCalibDrawerDat< float >::setData ( unsigned int channel, unsigned int adc, unsigned int idx, float data )

inherited

Sets a single T belonging to a channel/ADC.

Parameters

channel	The channel number
adc	The gain index
idx	The index of the value to set
data	The data value to set

Definition at line 71 of file TileCalibDrawerDat.h.

{
  //=== check for out of bounds
  if(channel >= getNChans()){
    throw TileCalib::IndexOutOfRange("TileCalibDrawerDat::getAddress(channel)",channel,getNChans());
  }
  if(adc>=getNGains()){
    throw TileCalib::IndexOutOfRange("CaloCalibDrawerDat::getAddress(gain)",adc,getNGains());
  }
  if(idx>=getObjSizeUint32()){
    throw TileCalib::IndexOutOfRange("TileCalibDrawerDat::getData",idx,getObjVersion());
  }
 
  getAddress(channel,adc)[idx] = data;
}

Member Data Documentation

m_blob

const coral::Blob* TileCalibDrawerBase::m_blob

privateinherited

Const reference to the BLOB (always there)

Definition at line 176 of file TileCalibDrawerBase.h.

m_blob_nc

coral::Blob* TileCalibDrawerBase::m_blob_nc

privateinherited

Non-const reference to the BLOB.

(Only present if we were created with a non-const blob.)

Definition at line 174 of file TileCalibDrawerBase.h.

m_blobSize32

uint64_t TileCalibDrawerBase::m_blobSize32

privateinherited

Cache blob size in units of uint32_t.

Definition at line 182 of file TileCalibDrawerBase.h.

m_blobStart16

const uint16_t* TileCalibDrawerBase::m_blobStart16

privateinherited

Cache blob starting address as uint_16t*.

Definition at line 180 of file TileCalibDrawerBase.h.

m_blobStart32

const uint32_t* TileCalibDrawerBase::m_blobStart32

privateinherited

Cache blob starting address as uint_32t*.

Definition at line 178 of file TileCalibDrawerBase.h.

m_hdrSize32

const unsigned int TileCalibDrawerBase::m_hdrSize32

staticinherited

The header size in units of uint32_t.

Definition at line 138 of file TileCalibDrawerBase.h.

m_isBlobOwner

bool TileCalibDrawerBase::m_isBlobOwner

privateinherited

Is this TileCalibDrawer owner of the BLOB.

Definition at line 184 of file TileCalibDrawerBase.h.

---

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

• TileCalibDrawerFlt.h
• TileCalibDrawerFlt.cxx