LArDeltaRespTool Class Reference

#include <LArDeltaRespTool.h>

Inheritance diagram for LArDeltaRespTool:

Collaboration diagram for LArDeltaRespTool:

Public Member Functions
	LArDeltaRespTool (const std::string &type, const std::string &name, const IInterface *parent)
virtual	~LArDeltaRespTool ()
virtual StatusCode	initialize ()
virtual StatusCode	finalize ()
LArCaliWave	makeLArDeltaResp (const LArWFParams &, const LArCaliWave &)
ServiceHandle< StoreGateSvc > &	evtStore ()
	The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.
const ServiceHandle< StoreGateSvc > &	detStore () const
	The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.
virtual StatusCode	sysInitialize () override
	Perform system initialization for an algorithm.
virtual StatusCode	sysStart () override
	Handle START transition.
virtual std::vector< Gaudi::DataHandle * >	inputHandles () const override
	Return this algorithm's input handles.
virtual std::vector< Gaudi::DataHandle * >	outputHandles () const override
	Return this algorithm's output handles.
Gaudi::Details::PropertyBase &	declareProperty (Gaudi::Property< T, V, H > &t)
void	updateVHKA (Gaudi::Details::PropertyBase &)
MsgStream &	msg () const
bool	msgLvl (const MSG::Level lvl) const

Static Public Member Functions
static const InterfaceID &	interfaceID ()

Protected Member Functions
void	renounceArray (SG::VarHandleKeyArray &handlesArray)
	remove all handles from I/O resolution
std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void >	renounce (T &h)
void	extraDeps_update_handler (Gaudi::Details::PropertyBase &ExtraDeps)
	Add StoreName to extra input/output deps as needed.

Private Types
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
void	compute_deltaresp ()
LArWave	deltaResp (const LArWave &) const
LArWave	deltaCorr () const
double	deltaCorr (double t) const
LArWave	injResp (const LArWave &w) const
LArWave	injCorr () const
double	injCorr (double t) const
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey>

Private Attributes
LArCaliWave	m_gCali
LArCaliWave	m_gDelta
bool	m_injPointCorr
bool	m_normalizeCali
bool	m_timeOriginShift
bool	m_subtractBaseline
double	m_Tdrift
double	m_Fstep
double	m_Tcal
double	m_Omega0
double	m_Taur
unsigned	m_Tstart
StoreGateSvc_t	m_evtStore
	Pointer to StoreGate (event store by default)
StoreGateSvc_t	m_detStore
	Pointer to StoreGate (detector store by default)
std::vector< SG::VarHandleKeyArray * >	m_vhka
bool	m_varHandleArraysDeclared

Static Private Attributes
static const int	DEFAULT =-1

Detailed Description

Definition at line 19 of file LArDeltaRespTool.h.

Member Typedef Documentation

StoreGateSvc_t

typedef ServiceHandle<StoreGateSvc> AthCommonDataStore< AthCommonMsg< AlgTool > >::StoreGateSvc_t

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

LArDeltaRespTool()

LArDeltaRespTool::LArDeltaRespTool	(	const std::string &	type,
		const std::string &	name,
		const IInterface *	parent )

Definition at line 11 of file LArDeltaRespTool.cxx.

  : AthAlgTool(type,name,parent),
    m_Tdrift(0),
    m_Fstep(0),
    m_Tcal(0),
    m_Omega0(0),
    m_Taur(0),
    m_Tstart(0)
{
  declareProperty("InjPointCorr",m_injPointCorr=true) ;
  declareProperty("NormalizeCali",m_normalizeCali=true) ;
  declareProperty("TimeOriginShift",m_timeOriginShift=false) ;
  declareProperty("SubtractBaseline",m_subtractBaseline=true) ;
}

~LArDeltaRespTool()

LArDeltaRespTool::~LArDeltaRespTool ( )

virtualdefault

Member Function Documentation

compute_deltaresp()

void LArDeltaRespTool::compute_deltaresp ( )

private

Definition at line 47 of file LArDeltaRespTool.cxx.

{
  LArWaveHelper wHelper;
 
  // calib. signal at Mother Board :
  LArWave gCaliMB(m_gCali), gDelta;  
 
  // shift gCaliMB to start point and remove baseline
  m_Tstart = wHelper.getStart(gCaliMB) ;
  double baseline = wHelper.getBaseline(gCaliMB,m_Tstart) ;
  if ( m_subtractBaseline ) gCaliMB = gCaliMB + (-baseline) ;
  if ( m_timeOriginShift )  gCaliMB = wHelper.translate(gCaliMB,-m_Tstart,baseline) ;
  
  // normalization of calibration pulse
  if ( m_normalizeCali ) {
    ATH_MSG_VERBOSE ( "*** Normalisation \t|-> YES" );
    double peak = gCaliMB.getSample( wHelper.getMax(gCaliMB) ) ;
    ATH_MSG_VERBOSE  ( "peak = " << peak );
    if ( peak <=0 ) {
      ATH_MSG_WARNING ( "Peak value <=0 , cannot normalize!" );
    } else {
      gCaliMB = gCaliMB * (1./peak)  ;
    }
  } else {
    ATH_MSG_VERBOSE ( "*** Normalisation \t|-> NO" );
  }
 
  // delta responde waveform prediction
  ATH_MSG_VERBOSE ( "*** Delta response \t|-> m_Fstep   = " << m_Fstep  << " ns " );
  ATH_MSG_VERBOSE ( "*** Delta response \t|-> m_Tcal    = " << m_Tcal   << " ns " );
  ATH_MSG_VERBOSE ( "*** Delta response \t|-> m_Omega0  = " << m_Omega0 << " GHz" );
  ATH_MSG_VERBOSE ( "*** Delta response \t|-> m_Taur    = " << m_Taur   << " ns " );
  
  if ( ! m_injPointCorr ) {
    // perform only exp->triangle correction
    ATH_MSG_VERBOSE ( "*** Inj.Point Corr. \t|-> NO" );
    gDelta = deltaResp ( gCaliMB ) ;
  } else {
    // perform exp->triangle and then injection point correction
    ATH_MSG_VERBOSE ( "*** Inj.Point Corr. \t|-> YES" );
    //gDelta = deltaResp ( gCaliMB ) ;
    //gDelta = LArPhysWaveTool::injResp( gDelta ); // HOW TO MAKE THIS WORKING???
    gDelta = injResp( deltaResp ( gCaliMB ) ) ;
  }
  
  //int deltaRespDAC = m_gCali.getDAC() ;
  int deltaRespDAC = -3 ;
  
  m_gDelta = LArCaliWave( gDelta.getWave() ,
              m_gCali.getDt() ,
                  deltaRespDAC ,
              m_gCali.getIsPulsedInt(),
                  m_gCali.getFlag() );
             
}

declareGaudiProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty ( Gaudi::Property< T, V, H > & hndl, const SG::VarHandleKeyType &  )

inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleKey>

Definition at line 156 of file AthCommonDataStore.h.

  {
    return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                       hndl.value(), 
                                                       hndl.documentation());
 
  }

declareProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( Gaudi::Property< T, V, H > & t )

inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

                                                                     {
    typedef typename SG::HandleClassifier<T>::type htype;
    return AthCommonDataStore<PBASE>::declareGaudiProperty(t, htype());
  }

deltaCorr() [1/2]

LArWave LArDeltaRespTool::deltaCorr ( ) const

private

Definition at line 112 of file LArDeltaRespTool.cxx.

{
  unsigned N = m_gCali.getSize() ;
  double  dt = m_gCali.getDt() ;
  LArWave w(N,dt) ;
  for ( unsigned i=0 ; i<N ; i++ ) w.setSample(i,deltaCorr(i*dt)) ;
  return w ;
}

deltaCorr() [2/2]

double LArDeltaRespTool::deltaCorr ( double t ) const

private

Definition at line 121 of file LArDeltaRespTool.cxx.

{
  double fstep = m_Fstep ;
  double Tc    = m_Tcal ;
  return ((fstep*fstep-fstep)/(Tc*Tc)) * exp( -fstep*t/Tc );
}

deltaResp()

LArWave LArDeltaRespTool::deltaResp ( const LArWave & w ) const

private

Definition at line 104 of file LArDeltaRespTool.cxx.

{
  double fstep = m_Fstep;
  double Tc    = m_Tcal;
  LArWaveHelper wHelper;
  return ( w % deltaCorr() ) + w*((1.-fstep)/Tc) + wHelper.derive_smooth(w);
}

detStore()

const ServiceHandle< StoreGateSvc > & AthCommonDataStore< AthCommonMsg< AlgTool > >::detStore ( ) const

inlineinherited

The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 95 of file AthCommonDataStore.h.

{ return m_detStore; }

evtStore()

ServiceHandle< StoreGateSvc > & AthCommonDataStore< AthCommonMsg< AlgTool > >::evtStore ( )

inlineinherited

The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 85 of file AthCommonDataStore.h.

{ return m_evtStore; }

extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::extraDeps_update_handler ( Gaudi::Details::PropertyBase & ExtraDeps )

protectedinherited

Add StoreName to extra input/output deps as needed.

use the logic of the VarHandleKey to parse the DataObjID keys supplied via the ExtraInputs and ExtraOuputs Properties to add the StoreName if it's not explicitly given

finalize()

virtual StatusCode LArDeltaRespTool::finalize ( )

inlinevirtual

Definition at line 33 of file LArDeltaRespTool.h.

{return StatusCode::SUCCESS;}

initialize()

virtual StatusCode LArDeltaRespTool::initialize ( )

inlinevirtual

Definition at line 32 of file LArDeltaRespTool.h.

{return StatusCode::SUCCESS;}

injCorr() [1/2]

LArWave LArDeltaRespTool::injCorr ( ) const

private

Definition at line 132 of file LArDeltaRespTool.cxx.

                                        {
  unsigned N = m_gCali.getSize() ;
  double dt = m_gCali.getDt() ;
  LArWave w(N,dt) ;
  for ( unsigned i=0 ; i<N ; i++ ) w.setSample(i,injCorr(i*dt)) ;
  return w ;
}

injCorr() [2/2]

double LArDeltaRespTool::injCorr ( double t ) const

private

Definition at line 140 of file LArDeltaRespTool.cxx.

                                                  {
  double tau0 = 1./m_Omega0;
  double taur = m_Taur;
  double Delta = pow(taur,2.) - pow(2*tau0,2.) ;
  if ( Delta > 0 ) {
    double sqrtDelta = sqrt(Delta) ;
    double taup = 0.5*( taur + sqrtDelta ) ;
    double taum = 0.5*( taur - sqrtDelta ) ;
    return ( exp(-t/taup) - exp(-t/taum) ) / ( taup - taum ) ;
  } else if ( Delta < 0 ) {
    double T = sqrt(-Delta) ;
    double A = 2 * taur / ( pow(taur,2.) - Delta ) ;
    double B = 2 * T / ( pow(taur,2.) - Delta ) ;
    return 2 * exp(-A*t) * sin(B*t) / T ;
  } else {
    double tau = 0.5 * taur ;
    return exp(-t/tau) * t / pow(tau,2.) ;
  }
#if 0
  double taur2 = taur*taur, tau02 = tau0*tau0 ;
  double taua = sqrt( 4.*tau02 - taur2 );
  return (2./taua)*exp(-t*taur/(2.*tau02))*sin(t*taua/(2.*tau02));
#endif
}

injResp()

LArWave LArDeltaRespTool::injResp ( const LArWave & w ) const

private

Definition at line 128 of file LArDeltaRespTool.cxx.

                                                         {
  return  w % injCorr() ;
}

inputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< AlgTool > >::inputHandles ( ) const

overridevirtualinherited

Return this algorithm's input handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

interfaceID()

const InterfaceID & LArDeltaRespTool::interfaceID ( )

inlinestatic

Definition at line 24 of file LArDeltaRespTool.h.

{ return IID_LArDeltaRespTool; }

makeLArDeltaResp()

LArCaliWave LArDeltaRespTool::makeLArDeltaResp	(	const LArWFParams &	larWFParam,
		const LArCaliWave &	larCaliWave )

Definition at line 29 of file LArDeltaRespTool.cxx.

{
 
  m_gCali  = larCaliWave;
 
  m_Tdrift = larWFParam.tdrift();
  m_Tcal   = larWFParam.tcal();
  m_Fstep  = larWFParam.fstep();
  m_Taur   = larWFParam.taur();
  m_Omega0 = larWFParam.omega0();
 
  compute_deltaresp(); 
  
  return m_gDelta; 
}

msg()

MsgStream & AthCommonMsg< AlgTool >::msg ( ) const

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                {
    return this->msgStream();
  }

msgLvl()

bool AthCommonMsg< AlgTool >::msgLvl ( const MSG::Level lvl ) const

inlineinherited

Definition at line 30 of file AthCommonMsg.h.

                                               {
    return this->msgLevel(lvl);
  }

outputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< AlgTool > >::outputHandles ( ) const

overridevirtualinherited

Return this algorithm's output handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

renounce()

std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void > AthCommonDataStore< AthCommonMsg< AlgTool > >::renounce ( T & h )

inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

  {
    h.renounce();
    PBASE::renounce (h);
  }

renounceArray()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::renounceArray ( SG::VarHandleKeyArray & handlesArray )

inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

                                                          {
    handlesArray.renounce();
  }

sysInitialize()

virtual StatusCode AthCommonDataStore< AthCommonMsg< AlgTool > >::sysInitialize ( )

overridevirtualinherited

Perform system initialization for an algorithm.

We override this to declare all the elements of handle key arrays at the end of initialization. See comments on updateVHKA.

Reimplemented in asg::AsgMetadataTool, AthCheckedComponent< AthAlgTool >, AthCheckedComponent<::AthAlgTool >, and DerivationFramework::CfAthAlgTool.

sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< AlgTool > >::sysStart ( )

overridevirtualinherited

Handle START transition.

We override this in order to make sure that conditions handle keys can cache a pointer to the conditions container.

updateVHKA()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::updateVHKA ( Gaudi::Details::PropertyBase & )

inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

                                               {
    // debug() << "updateVHKA for property " << p.name() << " " << p.toString() 
    //         << "  size: " << m_vhka.size() << endmsg;
    for (auto &a : m_vhka) {
      std::vector<SG::VarHandleKey*> keys = a->keys();
      for (auto k : keys) {
        k->setOwner(this);
      }
    }
  }

Member Data Documentation

DEFAULT

const int LArDeltaRespTool::DEFAULT =-1

staticprivate

Definition at line 39 of file LArDeltaRespTool.h.

m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore

privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_evtStore

privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

m_Fstep

double LArDeltaRespTool::m_Fstep

private

Definition at line 45 of file LArDeltaRespTool.h.

m_gCali

LArCaliWave LArDeltaRespTool::m_gCali

private

Definition at line 41 of file LArDeltaRespTool.h.

m_gDelta

LArCaliWave LArDeltaRespTool::m_gDelta

private

Definition at line 41 of file LArDeltaRespTool.h.

m_injPointCorr

bool LArDeltaRespTool::m_injPointCorr

private

Definition at line 43 of file LArDeltaRespTool.h.

m_normalizeCali

bool LArDeltaRespTool::m_normalizeCali

private

Definition at line 43 of file LArDeltaRespTool.h.

m_Omega0

double LArDeltaRespTool::m_Omega0

private

Definition at line 45 of file LArDeltaRespTool.h.

m_subtractBaseline

bool LArDeltaRespTool::m_subtractBaseline

private

Definition at line 43 of file LArDeltaRespTool.h.

m_Taur

double LArDeltaRespTool::m_Taur

private

Definition at line 45 of file LArDeltaRespTool.h.

m_Tcal

double LArDeltaRespTool::m_Tcal

private

Definition at line 45 of file LArDeltaRespTool.h.

m_Tdrift

double LArDeltaRespTool::m_Tdrift

private

Definition at line 45 of file LArDeltaRespTool.h.

m_timeOriginShift

bool LArDeltaRespTool::m_timeOriginShift

private

Definition at line 43 of file LArDeltaRespTool.h.

m_Tstart

unsigned LArDeltaRespTool::m_Tstart

private

Definition at line 46 of file LArDeltaRespTool.h.

m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< AlgTool > >::m_varHandleArraysDeclared

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

m_vhka

std::vector<SG::VarHandleKeyArray*> AthCommonDataStore< AthCommonMsg< AlgTool > >::m_vhka

privateinherited

Definition at line 398 of file AthCommonDataStore.h.

---

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

• LArDeltaRespTool.h
• LArDeltaRespTool.cxx