egammaCellUtils::MaxECell Struct Reference

#include <findMaxECell.h>

Inheritance diagram for egammaCellUtils::MaxECell:

Collaboration diagram for egammaCellUtils::MaxECell:

Public Member Functions
	MaxECell ()=delete
	MaxECell (const xAOD::CaloCluster *clus, const std::string &cellCKey="AllCalo", bool UseWeightForMaxCell=false)
bool	msgLvl (const MSG::Level lvl) const
	Test the output level.
MsgStream &	msg () const
	The standard message stream.
MsgStream &	msg (const MSG::Level lvl) const
	The standard message stream.
void	setLevel (MSG::Level lvl)
	Change the current logging level.

Public Attributes
double	etaCell = 999
double	phiCell = 999
StatusCode	sc = StatusCode::FAILURE

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

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

Detailed Description

Definition at line 14 of file findMaxECell.h.

Constructor & Destructor Documentation

MaxECell() [1/2]

egammaCellUtils::MaxECell::MaxECell ( )

delete

MaxECell() [2/2]

egammaCellUtils::MaxECell::MaxECell	(	const xAOD::CaloCluster *	clus,
		const std::string &	cellCKey = "AllCalo",
		bool	UseWeightForMaxCell = false )

Definition at line 12 of file findMaxECell.cxx.

                                                        :
    AthMessaging("egammaCellUtils::MaxECell") {
    
    const CaloClusterCellLink* cellLinks = clus->getCellLinks();
    if (!cellLinks) {
      ATH_MSG_WARNING("No cell link for cluster. Do nothing");
    } else {
    
      // First check :
      const CaloCellContainer *caloCells = cellLinks->getCellContainer();
      if (!caloCells) {
    ATH_MSG_WARNING("No cells for cluster. Do nothing");
      } else if (cellLinks->getCellContainerLink().dataID() != cellCKey) {
    // Second check :
    // one might have used a custom cell container for the linked cells
    ATH_MSG_ERROR("Wrong key for the calo cells");
      } else {
    
    CaloClusterCellLink::const_iterator it_cell = cellLinks->begin(),
      it_cell_e = cellLinks->end();
    
    // find maximum cell energy in layer 2
    double emax = 0.;
    std::pair<const CaloCell*,double> maxcell{nullptr,0}; //just for debug
    for(; it_cell != it_cell_e; ++it_cell) {
      const CaloCell* cell = (*it_cell);
      if (cell) {
        if (!cell->caloDDE()) {
          ATH_MSG_WARNING("Calo cell without detector element ?? eta = "
                  << cell->eta() << " phi = " << cell->phi());
          continue;
        }
        int layer = cell->caloDDE()->getSampling();
        if (layer == CaloSampling::EMB2 || layer == CaloSampling::EME2) {
          double w     = it_cell.weight();
          double eCell = cell->energy();
          if (UseWeightForMaxCell) eCell *= w;
          if (eCell > emax) {
        emax           = eCell;
        maxcell.first  = cell;
        maxcell.second = w;
          }
        }
      }
    }
    
    if (emax > 0) {
      etaCell = maxcell.first->caloDDE()->eta_raw();
      phiCell = maxcell.first->caloDDE()->phi_raw();
      if (msgLvl(MSG::DEBUG)) {
        CaloSampling::CaloSample sam = maxcell.first->caloDDE()->getSampling();
        double etaAmax = clus->etamax(sam);
        double phiAmax = clus->phimax(sam);
        double vemax   = clus->energy_max(sam);
        ATH_MSG_DEBUG("Cluster energy in sampling 2 = " << clus->energyBE(2)
              << " maximum layer 2 energy cell, E = " << maxcell.first->energy()
              << " check E = " << vemax
              << " w = " << maxcell.second << "\n"
              << " in calo  frame, eta = " << etaCell << " phi = " << phiCell << "\n"
              << " in ATLAS frame, eta = " << etaAmax << " phi = " << phiAmax);
      }
    } else {
      ATH_MSG_WARNING("No layer 2 cell with positive energy ! Should never happen");
      sc = StatusCode::FAILURE;
    }
    sc = StatusCode::SUCCESS;
      } // links OK, good cell container
    }   // no links 
  } // end function

Member Function Documentation

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();
  // If user did not set an explicit level, set a default
  if (m_lvl == MSG::NIL) {
    m_lvl = m_imsg ?
      static_cast<MSG::Level>( m_imsg.load()->outputLevel(m_nm) ) :
      MSG::INFO;
  }
}

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 167 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 182 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 user did not set explicit message level we have to initialize
  // the messaging and retrieve the default via the MessageSvc.
  if (m_lvl==MSG::NIL && !m_initialized.test_and_set()) initMessaging();
 
  if (m_lvl <= lvl) {
    msg() << lvl;
    return true;
  } else {
    return false;
  }
}

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;
}

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.

etaCell

double egammaCellUtils::MaxECell::etaCell = 999

Definition at line 16 of file findMaxECell.h.

m_imsg

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

mutableprivateinherited

MessageSvc pointer.

Definition at line 135 of file AthMessaging.h.

{ nullptr };

m_lvl

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

mutableprivateinherited

Current logging level.

Definition at line 138 of file AthMessaging.h.

{ MSG::NIL };

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.

phiCell

double egammaCellUtils::MaxECell::phiCell = 999

Definition at line 17 of file findMaxECell.h.

sc

StatusCode egammaCellUtils::MaxECell::sc = StatusCode::FAILURE

Definition at line 18 of file findMaxECell.h.

---

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

• findMaxECell.h
• findMaxECell.cxx