TilePulseShape Class Reference

#include <TilePulseShape.h>

Inheritance diagram for TilePulseShape:

Collaboration diagram for TilePulseShape:

Public Member Functions
	TilePulseShape (IMessageSvc *msgSvc, const std::string &name)
	TilePulseShape (IMessageSvc *msgSvc, const std::string &name, const TString &fileName)
	TilePulseShape (IMessageSvc *msgSvc, const std::string &name, const std::vector< double > &shapevec)
virtual	~TilePulseShape ()
void	loadPulseShape (const TString &fileName)
void	setPulseShape (const std::vector< double > &shapevec)
TGraph *	getGraph (double t0=0., double ped=0., double amp=1.)
double	eval (double x, bool useSpline=true, bool useUndershoot=false)
void	resetDeformation ()
int	scalePulse (double leftSF=1., double rightSF=1.)
bool	msgLvl (const MSG::Level lvl) const
	Test the output level. More...
MsgStream &	msg () const
	The standard message stream. More...
MsgStream &	msg (const MSG::Level lvl) const
	The standard message stream. More...
void	setLevel (MSG::Level lvl)
	Change the current logging level. More...

Private Member Functions
void	initMessaging () const
	Initialize our message level and MessageSvc. More...

Private Attributes
TGraph *	m_pulseShape {nullptr}
TGraph *	m_deformedShape {nullptr}
TSpline *	m_deformedSpline {nullptr}
std::string	m_nm
	Message source name. More...
boost::thread_specific_ptr< MsgStream >	m_msg_tls
	MsgStream instance (a std::cout like with print-out levels) More...
std::atomic< IMessageSvc * >	m_imsg { nullptr }
	MessageSvc pointer. More...
std::atomic< MSG::Level >	m_lvl { MSG::NIL }
	Current logging level. More...
std::atomic_flag m_initialized	ATLAS_THREAD_SAFE = ATOMIC_FLAG_INIT
	Messaging initialized (initMessaging) More...

Detailed Description

Definition at line 17 of file TilePulseShape.h.

Constructor & Destructor Documentation

TilePulseShape() [1/3]

TilePulseShape::TilePulseShape	(	IMessageSvc *	msgSvc,
		const std::string &	name
	)

Definition at line 19 of file TilePulseShape.cxx.

  : AthMessaging(msgSvc, name)
{
}

TilePulseShape() [2/3]

TilePulseShape::TilePulseShape	(	IMessageSvc *	msgSvc,
		const std::string &	name,
		const TString &	fileName
	)

Definition at line 25 of file TilePulseShape.cxx.

  : AthMessaging(msgSvc, name)
{
  loadPulseShape(fileName);
} 

TilePulseShape() [3/3]

TilePulseShape::TilePulseShape	(	IMessageSvc *	msgSvc,
		const std::string &	name,
		const std::vector< double > &	shapevec
	)

Definition at line 32 of file TilePulseShape.cxx.

  : AthMessaging(msgSvc, name)
{
  setPulseShape(shapevec);
}

~TilePulseShape()

TilePulseShape::~TilePulseShape ( )

virtual

Definition at line 39 of file TilePulseShape.cxx.

{
  resetDeformation();
  if(m_pulseShape) delete m_pulseShape;
  if(m_deformedSpline) delete m_deformedSpline;
}

Member Function Documentation

eval()

double TilePulseShape::eval	(	double	x,
		bool	useSpline = true,
		bool	useUndershoot = false
	)

Definition at line 75 of file TilePulseShape.cxx.

{
 
  //=== make sure pulseshape is available
  if(!m_deformedShape){
    ATH_MSG_ERROR("No pulseshape loaded!");
    return 0.;
  }
 
  //=== interpolate pulseShape value, TMath::BinarySearch returns 
  //=== index of value smaller or equal to the search value. 
  //=== -> Need to catch boundary values
  double y(0.);
  int n(m_deformedShape->GetN());
  int idx = TMath::BinarySearch(n,m_deformedShape->GetX(),x);
  if(idx<0){
    //=== left out of bounds, return leftmost value
//    y = (_deformedShape->GetY())[0]; 
    y = 0;
    ATH_MSG_DEBUG("Left out of bounds. Replacing y = " << (m_deformedShape->GetY())[0] << " with y = 0. (idx = " << idx << ", x = " << x << ")");
  }
  else if(idx>=n-1){
    //=== right out of bounds, return rightmost value
//    y = (_deformedShape->GetY())[n-1]; 
    if(useUndershoot){
      y = 0.00196 * (1 - exp(- x / 20664.59) ) - 0.00217;
    }
    else{
      y = 0;
      ATH_MSG_DEBUG("Right out of bounds. Replacing y = " << (m_deformedShape->GetY())[0] << " with y = 0. (idx = " << idx << ", x = " << x << ")");
    }
  }
  else{
    //=== linear interpolation
    double xLo = (m_deformedShape->GetX())[idx  ];
    double xHi = (m_deformedShape->GetX())[idx+1];
    double yLo = (m_deformedShape->GetY())[idx  ];
    double yHi = (m_deformedShape->GetY())[idx+1];
    double yLinear = yLo + (yHi-yLo)/(xHi-xLo) * (x-xLo);
    //=== spline interpolation
    double ySpline = m_deformedSpline->Eval(x);
 
    if(useSpline) y = ySpline;
    else          y = yLinear;
  }
 
  return y;
}

getGraph()

TGraph * TilePulseShape::getGraph	(	double	t0 = 0.,
		double	ped = 0.,
		double	amp = 1.
	)

Definition at line 165 of file TilePulseShape.cxx.

{
  TGraph* gr = (TGraph*) m_deformedShape->Clone();
  for(int i=0; i<gr->GetN(); i++){
    double x,y;
    gr->GetPoint(i,x,y);
    y*=amp;
    y+=ped;
    x+=t0;
    gr->SetPoint(i,x,y);
  }
  return gr;
}

initMessaging()

void AthMessaging::initMessaging ( ) const

privateinherited

Initialize our message level and MessageSvc.

This method should only be called once.

Definition at line 39 of file AthMessaging.cxx.

{
  m_imsg = Athena::getMessageSvc();
  m_lvl = m_imsg ?
    static_cast<MSG::Level>( m_imsg.load()->outputLevel(m_nm) ) :
    MSG::INFO;
}

loadPulseShape()

void TilePulseShape::loadPulseShape

(

const TString &

fileName

)

Definition at line 47 of file TilePulseShape.cxx.

{
  resetDeformation();
  if(m_pulseShape) delete m_pulseShape;
  m_pulseShape = new TGraph(fileName.Data());
  if(m_pulseShape->IsZombie()) {
    throw GaudiException(std::string("could not load pulseshape from file: ") + fileName.Data(),
                         "TilePulseShape", StatusCode::FAILURE);
  } else ATH_MSG_INFO("Loaded pulseshape from file: " << fileName);
  resetDeformation();
}

msg() [1/2]

MsgStream & AthMessaging::msg ( ) const

inlineinherited

The standard message stream.

Returns a reference to the default message stream May not be invoked before sysInitialize() has been invoked.

Definition at line 164 of file AthMessaging.h.

{
  MsgStream* ms = m_msg_tls.get();
  if (!ms) {
    if (!m_initialized.test_and_set()) initMessaging();
    ms = new MsgStream(m_imsg,m_nm);
    m_msg_tls.reset( ms );
  }
 
  ms->setLevel (m_lvl);
  return *ms;
}

msg() [2/2]

MsgStream & AthMessaging::msg ( const MSG::Level  lvl ) const

inlineinherited

The standard message stream.

Returns a reference to the default message stream May not be invoked before sysInitialize() has been invoked.

Definition at line 179 of file AthMessaging.h.

{ return msg() << lvl; }

msgLvl()

bool AthMessaging::msgLvl ( const MSG::Level  lvl ) const

inlineinherited

Test the output level.

Parameters

lvl	The message level to test against

Returns

boolean Indicating if messages at given level will be printed

Return values

true	Messages at level "lvl" will be printed

Definition at line 151 of file AthMessaging.h.

{
  if (!m_initialized.test_and_set()) initMessaging();
  if (m_lvl <= lvl) {
    msg() << lvl;
    return true;
  } else {
    return false;
  }
}

resetDeformation()

void TilePulseShape::resetDeformation

(

)

Definition at line 126 of file TilePulseShape.cxx.

{
  if(m_deformedShape==m_pulseShape) return;
  
  delete m_deformedShape;  
  delete m_deformedSpline; 
  m_deformedShape = m_pulseShape;
  m_deformedSpline = new TSpline3("deformedSpline",m_deformedShape);
}

scalePulse()

int TilePulseShape::scalePulse	(	double	leftSF = 1.,
		double	rightSF = 1.
	)

Definition at line 138 of file TilePulseShape.cxx.

{
 
  resetDeformation();
  if(!m_pulseShape) {
    ATH_MSG_WARNING("Attempted pulse shape scaling before loading pulse shape");
    return 1;
   } else {
    m_deformedShape = (TGraph*) m_pulseShape->Clone();
  
    for(int i=0; i<m_deformedShape->GetN(); ++i){
      double x,y;
      m_deformedShape->GetPoint(i,x,y);
      if(x<0.)      x*= leftSF;
      else if(x>0.) x*=rightSF;
      m_deformedShape->SetPoint(i,x,y);
    }
 
    delete m_deformedSpline;
    m_deformedSpline = new TSpline3("deformedSpline",m_deformedShape);
    return 0;
  }   
 
}

setLevel()

void AthMessaging::setLevel ( MSG::Level  lvl )

inherited

Change the current logging level.

Use this rather than msg().setLevel() for proper operation with MT.

Definition at line 28 of file AthMessaging.cxx.

{
  m_lvl = lvl;
}

setPulseShape()

void TilePulseShape::setPulseShape

(

const std::vector< double > &

shapevec

)

Definition at line 60 of file TilePulseShape.cxx.

{
  resetDeformation();
  if(m_pulseShape) delete m_pulseShape;
  m_pulseShape = new TGraph(shapevec.size());
  for(std::vector<double>::size_type i = 0; i != shapevec.size(); i++) {
    m_pulseShape->SetPoint(i, -75.5+i*0.5, shapevec.at(i));
  }
  resetDeformation();
}

Member Data Documentation

ATLAS_THREAD_SAFE

std::atomic_flag m_initialized AthMessaging::ATLAS_THREAD_SAFE = ATOMIC_FLAG_INIT

mutableprivateinherited

Messaging initialized (initMessaging)

Definition at line 141 of file AthMessaging.h.

m_deformedShape

TGraph* TilePulseShape::m_deformedShape {nullptr}

private

Definition at line 36 of file TilePulseShape.h.

m_deformedSpline

TSpline* TilePulseShape::m_deformedSpline {nullptr}

private

Definition at line 37 of file TilePulseShape.h.

m_imsg

std::atomic<IMessageSvc*> AthMessaging::m_imsg { nullptr }

mutableprivateinherited

MessageSvc pointer.

Definition at line 135 of file AthMessaging.h.

m_lvl

std::atomic<MSG::Level> AthMessaging::m_lvl { MSG::NIL }

mutableprivateinherited

Current logging level.

Definition at line 138 of file AthMessaging.h.

m_msg_tls

boost::thread_specific_ptr<MsgStream> AthMessaging::m_msg_tls

mutableprivateinherited

MsgStream instance (a std::cout like with print-out levels)

Definition at line 132 of file AthMessaging.h.

m_nm

std::string AthMessaging::m_nm

privateinherited

Message source name.

Definition at line 129 of file AthMessaging.h.

m_pulseShape

TGraph* TilePulseShape::m_pulseShape {nullptr}

private

Definition at line 35 of file TilePulseShape.h.

---

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

• TilePulseShape.h
• TilePulseShape.cxx