CP::IsolationConditionGraph Class Reference

#include <IsolationConditionGraph.h>

Inheritance diagram for CP::IsolationConditionGraph:

Collaboration diagram for CP::IsolationConditionGraph:

Public Member Functions
	IsolationConditionGraph (const std::string &name, const std::vector< std::string > &isoType, std::unique_ptr< TF1 > isoFunction, std::vector< TGraph * > &cutGraphs, std::unique_ptr< TH1F > binning, const std::string &isoDecSuffix="", bool invertCut=false)
virtual	~IsolationConditionGraph () override=default
bool	accept (const xAOD::IParticle &x) const override
bool	accept (const strObj &x) const override
const std::string &	name () const
unsigned int	num_types () const
xAOD::Iso::IsolationType	type (unsigned int n=0) const
const FloatAccessor &	accessor (unsigned int n=0) const
const FloatAccessor &	accessor_noCloseBy (unsigned int n=0) const
void	setLevel (MSG::Level lvl)
	Change the current logging level.
Functions providing the same interface as AthMessaging
bool	msgLvl (const MSG::Level lvl) const
	Test the output level of the object.
MsgStream &	msg () const
	The standard message stream.
MsgStream &	msg (const MSG::Level lvl) const
	The standard message stream.

Protected Attributes
std::string	m_isoDecSuffix {}

Private Member Functions
float	getCutValue (const float pt, const float eta) const
void	initMessaging () const
	Initialize our message level and MessageSvc.

Private Attributes
std::vector< TGraph * >	m_cutGraphs
std::unique_ptr< TF1 >	m_isoFunction
std::unique_ptr< TH1F >	m_binning
bool	m_invertCut {false}
std::string	m_name
std::vector< xAOD::Iso::IsolationType >	m_isolationType
std::vector< FloatAccessor >	m_acc
std::vector< FloatAccessor >	m_acc_noCloseBy
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 15 of file IsolationConditionGraph.h.

Constructor & Destructor Documentation

IsolationConditionGraph()

CP::IsolationConditionGraph::IsolationConditionGraph	(	const std::string &	name,
		const std::vector< std::string > &	isoType,
		std::unique_ptr< TF1 >	isoFunction,
		std::vector< TGraph * > &	cutGraphs,
		std::unique_ptr< TH1F >	binning,
		const std::string &	isoDecSuffix = "",
		bool	invertCut = false )

Definition at line 15 of file IsolationConditionGraph.cxx.

    : IsolationCondition(name, isoTypes, isoDecSuffix),
      m_cutGraphs(cutGraphs),
      m_isoFunction(std::move(isoFunction)),
      m_binning(std::move(binning)),
      m_invertCut(invertCut) {}

~IsolationConditionGraph()

virtual CP::IsolationConditionGraph::~IsolationConditionGraph ( )

overridevirtualdefault

Member Function Documentation

accept() [1/2]

bool CP::IsolationConditionGraph::accept ( const strObj & x ) const

overridevirtual

Implements CP::IsolationCondition.

Definition at line 79 of file IsolationConditionGraph.cxx.

                                                          {
  const float cutValue = getCutValue(x.pt, x.eta);
  std::vector<double> isoVars;
  for (unsigned int itype = 0; itype < num_types(); ++itype)
    isoVars.push_back(x.isolationValues[type(itype)]);
  TF1* f ATLAS_THREAD_SAFE = m_isoFunction.get();
  const float isoValue = f->EvalPar(isoVars.data());
 
  ATH_MSG_DEBUG("PLIV Isolation values of electron %.3f " << isoValue);
 
  if (!m_invertCut)
    return isoValue <= cutValue;
  return isoValue > cutValue;
}

accept() [2/2]

bool CP::IsolationConditionGraph::accept ( const xAOD::IParticle & x ) const

overridevirtual

Implements CP::IsolationCondition.

Definition at line 28 of file IsolationConditionGraph.cxx.

                                                                 {
 
  const float cutValue = getCutValue(x.pt(), x.eta());
  std::vector<double> isoVars(num_types(), 0);
 
  for (unsigned int iacc = 0; iacc < num_types(); ++iacc) {
    const FloatAccessor& acc = accessor(iacc);
 
    if (!acc.isAvailable(x)) {
      // Temporary fix for missing closeByCorr variables if no primary vertex
      // exists for the event If closeByCorr variable does not exist, fallback
      // to the standard isolation variable 2025/02
      const FloatAccessor& acc_noCloseBy = accessor_noCloseBy(iacc);
      if (acc_noCloseBy.isAvailable(x)) {
        isoVars[iacc] = acc_noCloseBy(x);
      } else {
        ATH_MSG_WARNING(__FILE__
                        << ":" << __LINE__ << "Accessor "
                        << SG::AuxTypeRegistry::instance().getName(acc.auxid())
                        << " is not available. Expected when using primary "
                           "AODs, post-p3793 derivations (only for *FixedRad "
                           "or FixedCutPflow* for electrons), "
                        << " pre-p3517 derivations (only for FC*), or "
                           "pre-p3830 derivations (for other electron WPs)");
        if (!m_isoDecSuffix.empty())
          throw std::runtime_error(
              "IsolationConditionCombined: IsolationSelectionTool property "
              "'IsoDecSuffix' is set to " +
              m_isoDecSuffix +
              ". Must run on derivation made with IsolationCloseByCorrection "
              "to create the isolation variables with this suffix, or remove "
              "'IsoDecSuffix'. ");
        isoVars[iacc] = FLT_MAX;
      }
    } else
      isoVars[iacc] = acc(x);
  }
 
  // In general TF1::EvalPar is non-const/not thread-safe. But for the special
  // case of using it as TFormula it can be considered const/safe. It would be
  // safer to change the interface to not allow a generic TF1.
  TF1* f ATLAS_THREAD_SAFE = m_isoFunction.get();
  const float isoValue = f->EvalPar(isoVars.data());
 
  ATH_MSG_DEBUG("PLIV Isolation values of electron %.3f " << isoValue);
 
  if (!m_invertCut)
    return isoValue <= cutValue;
  return isoValue > cutValue;
}

accessor()

const SG::AuxElement::ConstAccessor< float > & CP::IsolationCondition::accessor ( unsigned int n = 0 ) const

inherited

Definition at line 44 of file IsolationCondition.cxx.

{ return m_acc.at(n); }

accessor_noCloseBy()

const SG::AuxElement::ConstAccessor< float > & CP::IsolationCondition::accessor_noCloseBy ( unsigned int n = 0 ) const

inherited

Definition at line 45 of file IsolationCondition.cxx.

{ return m_acc_noCloseBy.at(n); }

getCutValue()

float CP::IsolationConditionGraph::getCutValue ( const float pt, const float eta ) const

private

Definition at line 94 of file IsolationConditionGraph.cxx.

                                                                  {
  //
  int bin = m_binning->GetXaxis()->FindBin(fabs(eta));
  ATH_MSG_DEBUG(Form(" Electron (eta=%.3f) falls into bin %d", eta, bin - 1));
 
  if (bin <= 0 || (unsigned int)bin > m_cutGraphs.size()) {
    ATH_MSG_DEBUG(
        "Stored binning in file does not comply with binning of graphs for "
        "working point "
        << this->name() << "! Please contact the IFF group");
    return -999;
  }
  ATH_MSG_DEBUG(
      Form(" Electron cut value %.3f for pt %.3f as calculated from graph %s",
           (m_cutGraphs[bin - 1])->Eval(pt / 1000.), pt / 1000.,
           (m_cutGraphs[bin - 1])->GetName()));
 
  return (m_cutGraphs[bin - 1])->Eval(pt / 1000.);
}

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 & asg::AsgMessaging::msg ( ) const

inherited

The standard message stream.

Returns

A reference to the default message stream of this object.

Definition at line 49 of file AsgMessaging.cxx.

                                      {
#ifndef XAOD_STANDALONE
      return ::AthMessaging::msg();
#else // not XAOD_STANDALONE
      return m_msg;
#endif // not XAOD_STANDALONE
   }

msg() [2/2]

MsgStream & asg::AsgMessaging::msg ( const MSG::Level lvl ) const

inherited

The standard message stream.

Parameters

lvl	The message level to set the stream to

Returns

A reference to the default message stream, set to level "lvl"

Definition at line 57 of file AsgMessaging.cxx.

                                                          {
#ifndef XAOD_STANDALONE
      return ::AthMessaging::msg( lvl );
#else // not XAOD_STANDALONE
      m_msg << lvl;
      return m_msg;
#endif // not XAOD_STANDALONE
   }

msgLvl()

bool asg::AsgMessaging::msgLvl ( const MSG::Level lvl ) const

inherited

Test the output level of the object.

Parameters

lvl	The message level to test against

Returns

boolean Indicting if messages at given level will be printed

true If messages at level "lvl" will be printed

Definition at line 41 of file AsgMessaging.cxx.

                                                       {
#ifndef XAOD_STANDALONE
      return ::AthMessaging::msgLvl( lvl );
#else // not XAOD_STANDALONE
      return m_msg.msgLevel( lvl );
#endif // not XAOD_STANDALONE
   }

name()

const std::string & CP::IsolationCondition::name ( ) const

inherited

Definition at line 42 of file IsolationCondition.cxx.

{ return m_name; }

num_types()

unsigned int CP::IsolationCondition::num_types ( ) const

inherited

Definition at line 41 of file IsolationCondition.cxx.

{ return m_isolationType.size(); }

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

type()

xAOD::Iso::IsolationType CP::IsolationCondition::type ( unsigned int n = 0 ) const

inherited

Definition at line 43 of file IsolationCondition.cxx.

{ return m_isolationType[n]; }

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_acc

std::vector<FloatAccessor> CP::IsolationCondition::m_acc

privateinherited

Definition at line 54 of file IsolationCondition.h.

m_acc_noCloseBy

std::vector<FloatAccessor> CP::IsolationCondition::m_acc_noCloseBy

privateinherited

Definition at line 55 of file IsolationCondition.h.

m_binning

std::unique_ptr<TH1F> CP::IsolationConditionGraph::m_binning

private

Definition at line 32 of file IsolationConditionGraph.h.

m_cutGraphs

std::vector<TGraph*> CP::IsolationConditionGraph::m_cutGraphs

private

Definition at line 30 of file IsolationConditionGraph.h.

m_imsg

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

mutableprivateinherited

MessageSvc pointer.

Definition at line 135 of file AthMessaging.h.

{ nullptr };

m_invertCut

bool CP::IsolationConditionGraph::m_invertCut {false}

private

Definition at line 33 of file IsolationConditionGraph.h.

{false};

m_isoDecSuffix

std::string CP::IsolationCondition::m_isoDecSuffix {}

protectedinherited

Definition at line 58 of file IsolationCondition.h.

{};

m_isoFunction

std::unique_ptr<TF1> CP::IsolationConditionGraph::m_isoFunction

private

Definition at line 31 of file IsolationConditionGraph.h.

m_isolationType

std::vector<xAOD::Iso::IsolationType> CP::IsolationCondition::m_isolationType

privateinherited

Definition at line 53 of file IsolationCondition.h.

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_name

std::string CP::IsolationCondition::m_name

privateinherited

Definition at line 52 of file IsolationCondition.h.

m_nm

std::string AthMessaging::m_nm

privateinherited

Message source name.

Definition at line 129 of file AthMessaging.h.

---

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

• IsolationConditionGraph.h
• IsolationConditionGraph.cxx