LArAutoCorrDecoderTool Class Reference

#include <LArAutoCorrDecoderTool.h>

Inheritance diagram for LArAutoCorrDecoderTool:

Collaboration diagram for LArAutoCorrDecoderTool:

Public Member Functions
virtual	~LArAutoCorrDecoderTool ()
const Eigen::MatrixXd	AutoCorr (const HWIdentifier &CellID, int gain, unsigned nSamples) const
virtual StatusCode	initialize ()
virtual StatusCode	finalize ()

Static Public Member Functions
static const InterfaceID &	interfaceID ()

Private Member Functions
const Eigen::MatrixXd	ACDiagonal (const HWIdentifier &CellID, int gain, unsigned nSamples) const
const Eigen::MatrixXd	ACPhysics (const HWIdentifier &CellID, int gain, unsigned nSamples) const

Private Attributes
Gaudi::Property< unsigned >	m_decodemode {this,"DecodeMode",0}
Gaudi::Property< bool >	m_alwaysHighGain {this,"UseAlwaysHighGain",false}
Gaudi::Property< bool >	m_isSC {this,"isSC",false}
Gaudi::Property< std::string >	m_keyAutoCorr {this,"KeyAutoCorr","LArAutoCorr"}
const LArOnlineID_Base *	m_onlineID =nullptr

Detailed Description

Definition at line 21 of file LArAutoCorrDecoderTool.h.

Constructor & Destructor Documentation

~LArAutoCorrDecoderTool()

LArAutoCorrDecoderTool::~LArAutoCorrDecoderTool ( )

virtualdefault

Member Function Documentation

ACDiagonal()

const Eigen::MatrixXd LArAutoCorrDecoderTool::ACDiagonal ( const HWIdentifier & CellID, int gain, unsigned nSamples = 5 ) const

private

Definition at line 45 of file LArAutoCorrDecoderTool.cxx.

                                                                                                                         {
 
  if (!m_isSC && m_alwaysHighGain) {
    if (m_onlineID->isFCALchannel(CellID) ||m_onlineID->isHECchannel(CellID))
      gain=1;
    else
      gain=0;
  }
 
  Eigen::MatrixXd AutoCorrMatrix=Eigen::MatrixXd::Zero(nSamples,nSamples);
 
  const ILArAutoCorr* autoCorr=nullptr;
  detStore()->retrieve(autoCorr,m_keyAutoCorr).ignore();
  
 
  if ( autoCorr ) { // LArAutoCorrComplete is loaded in DetStore
 
    ILArAutoCorr::AutoCorrRef_t dbcorr = autoCorr->autoCorr(CellID,gain);
 
    if ( dbcorr.size()== 0 ) { // empty AutoCorr for given channel
      ATH_MSG_WARNING( "Empty AutoCorr vector for channel " <<  m_onlineID->channel_name(CellID) << " in Gain = " << gain);
      nSamples=0;
    }
    else if (dbcorr.size() < nSamples-1 ) {
      ATH_MSG_WARNING( "Not enough samples in AutoCorr vector for channel " <<  m_onlineID->channel_name(CellID) << " in Gain = " << gain);
      nSamples=1+dbcorr.size(); //The remaining values of the eigen matrix are left to 0.0
    } 
  
    // fill diagonal matrix with vector 
    for (unsigned i=0;i<nSamples;i++) {
      AutoCorrMatrix(i,i)= 1 ;
      for (unsigned j=i+1;j<nSamples;j++) {
    AutoCorrMatrix(i,j) = AutoCorrMatrix(j,i) = dbcorr[j-i-1];      
      }
    }
  }//else if m_autoCorr
  else { // no LArAutoCorrComplete loaded in DetStore (e.g. DB problem) :-(
    ATH_MSG_WARNING( "No valid AutoCorr object loaded from DetStore" );
  }
 
  ATH_MSG_DEBUG("AutoCorr Diagonal matrix for channel " <<  m_onlineID->channel_name(CellID) 
        << " in Gain = " << gain
        << ":\n" << AutoCorrMatrix);
 
  return AutoCorrMatrix;
  
}

ACPhysics()

const Eigen::MatrixXd LArAutoCorrDecoderTool::ACPhysics ( const HWIdentifier & CellID, int gain, unsigned nSamples = 5 ) const

private

Definition at line 93 of file LArAutoCorrDecoderTool.cxx.

                                                                                                                        {
 
 
  if (!m_isSC && m_alwaysHighGain) {
    if (m_onlineID->isFCALchannel(CellID) ||m_onlineID->isHECchannel(CellID))
      gain=1;
    else
      gain=0;
  }
 
  Eigen::MatrixXd AutoCorrMatrix=Eigen::MatrixXd::Identity(nSamples,nSamples);
 
  const ILArAutoCorr* autoCorr=nullptr;
  detStore()->retrieve(autoCorr,m_keyAutoCorr).ignore();
  
  if ( autoCorr ) { // LArAutoCorrComplete is loaded in DetStore
 
    ILArAutoCorr::AutoCorrRef_t corrdb = autoCorr->autoCorr(CellID,gain);
   
    if ( corrdb.size()== 0 ) { // empty AutoCorr for given channel
      ATH_MSG_WARNING( "Empty AutoCorr vector for channel " << m_onlineID->channel_name(CellID) << " in Gain = " << gain);
      nSamples=0; //return all-zero matrix
    }
    else  if ( corrdb.size() < nSamples*(nSamples+1)/2 ) {
      ATH_MSG_WARNING( "Not enough samples in AutoCorr vector for channel " <<  m_onlineID->channel_name(CellID) 
            << "in Gain = " << gain << " for AC Physics mode");
      nSamples=0;//return all-zero matrix 
    } 
 
    // Corr size could be bigger, then it's asked now, need remapping:
    const unsigned int nsamples_AC = (-1+((int)(sqrt(1+8*corrdb.size()))))/2;
    unsigned int k=0;
    for (unsigned i=0;i<nSamples;i++) {
      for (unsigned j=i;j<nSamples;j++,k++) {
    if (i<=j) {
      AutoCorrMatrix(i,j) = AutoCorrMatrix(j,i)= corrdb[k];
    }
      }
      k+=nsamples_AC-nSamples;
    }
  } //end if m_autoCorr 
  else { // no LArAutoCorrComplete loaded in DetStore (e.g. DB problem) :-(
    ATH_MSG_WARNING( "No valid AutoCorr object loaded from DetStore" );
  }
   
  ATH_MSG_DEBUG("AutoCorr Physics matrix for channel " <<  m_onlineID->channel_name(CellID) 
        << " in Gain = " << gain
        << ":\n" << AutoCorrMatrix);
  return AutoCorrMatrix;
  
}

AutoCorr()

const Eigen::MatrixXd LArAutoCorrDecoderTool::AutoCorr	(	const HWIdentifier &	CellID,
		int	gain,
		unsigned	nSamples = 5 ) const

Definition at line 36 of file LArAutoCorrDecoderTool.cxx.

{
  if (m_decodemode==1u)
    return ACPhysics(CellID,gain,nSamples);
  else
    return ACDiagonal(CellID,gain,nSamples);
}

finalize()

virtual StatusCode LArAutoCorrDecoderTool::finalize ( )

inlinevirtual

Definition at line 38 of file LArAutoCorrDecoderTool.h.

{return StatusCode::SUCCESS;}

initialize()

StatusCode LArAutoCorrDecoderTool::initialize ( )

virtual

Definition at line 10 of file LArAutoCorrDecoderTool.cxx.

{
  ATH_MSG_DEBUG("LArAutoCorrDecoderTool initialize() begin");
  
  if ( m_isSC ) {
    const LArOnline_SuperCellID* ll;
    ATH_CHECK( detStore()->retrieve(ll, "LArOnline_SuperCellID") );
    m_onlineID = static_cast<const LArOnlineID_Base*>(ll);
    ATH_MSG_DEBUG("Found the LArOnlineID helper");
    
  } else { // m_isSC
    const LArOnlineID* ll;
    ATH_CHECK( detStore()->retrieve(ll, "LArOnlineID") );
    m_onlineID = static_cast<const LArOnlineID_Base*>(ll);
    ATH_MSG_DEBUG(" Found the LArOnlineID helper. ");
  }
 
 
  if (!m_isSC && m_alwaysHighGain)
    ATH_MSG_INFO( "Will always return HIGH gain autocorrelation matrix for EM calo, MEDIUM for HEC and FCAL" );
 
  ATH_MSG_DEBUG("LArAutoCorrDecoderTool initialize() end");
  return StatusCode::SUCCESS;
}

interfaceID()

const InterfaceID & LArAutoCorrDecoderTool::interfaceID ( )

inlinestatic

Definition at line 40 of file LArAutoCorrDecoderTool.h.

                                          { 
    return ILArAutoCorrDecoderTool::interfaceID();
  } 

Member Data Documentation

m_alwaysHighGain

Gaudi::Property<bool> LArAutoCorrDecoderTool::m_alwaysHighGain {this,"UseAlwaysHighGain",false}

private

Definition at line 48 of file LArAutoCorrDecoderTool.h.

{this,"UseAlwaysHighGain",false};

m_decodemode

Gaudi::Property<unsigned> LArAutoCorrDecoderTool::m_decodemode {this,"DecodeMode",0}

private

Definition at line 47 of file LArAutoCorrDecoderTool.h.

{this,"DecodeMode",0};

m_isSC

Gaudi::Property<bool> LArAutoCorrDecoderTool::m_isSC {this,"isSC",false}

private

Definition at line 49 of file LArAutoCorrDecoderTool.h.

{this,"isSC",false};

m_keyAutoCorr

Gaudi::Property<std::string> LArAutoCorrDecoderTool::m_keyAutoCorr {this,"KeyAutoCorr","LArAutoCorr"}

private

Definition at line 50 of file LArAutoCorrDecoderTool.h.

{this,"KeyAutoCorr","LArAutoCorr"};

m_onlineID

const LArOnlineID_Base* LArAutoCorrDecoderTool::m_onlineID =nullptr

private

Definition at line 55 of file LArAutoCorrDecoderTool.h.

---

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

• LArAutoCorrDecoderTool.h
• LArAutoCorrDecoderTool.cxx