ATLAS Offline Software
Classes | Public Types | Public Member Functions | Static Public Member Functions | Protected Member Functions | Private Member Functions | Private Attributes | List of all members
TestMuonSF::MuonSFBranches Class Reference

#include <MuonSFTestHelper.h>

Inheritance diagram for TestMuonSF::MuonSFBranches:
Collaboration diagram for TestMuonSF::MuonSFBranches:

Classes

struct  SFSet
 Helper struct to store scale factor, data efficiency + mc efficiency for each systematic variation. More...
 

Public Types

using DataDependency = SG::VarHandleKey *
 

Public Member Functions

 MuonSFBranches (MuonVal::MuonTesterTree &tree, const ToolHandle< CP::IMuonEfficiencyScaleFactors > &handle, const std::string &rel_name="")
 
virtual ~MuonSFBranches ()=default
 
void setMuon (const xAOD::Muon &muon) override final
 
bool fill (const EventContext &ctx) override final
 The fill method checks if enough information is provided such that the branch is cleared from the information of the prevoius event and in cases of that updates did no happen befure a default value is provided. More...
 
bool init () override final
 The init method checks whether the branch name has already registered to the MuonTree and tries then to add the branch to the tree --> return false if the branch has already been defined or the inclusion of the branch to the tree failed or the name of the branch is empty. More...
 
std::string name () const override final
 Returns the name of the branch. More...
 
bool initialized () const
 Returns whether the object has been initialized or not The outcome of the variable is bound to a successful call of the addToTree method. More...
 
TTree * tree () override final
 Returns the underlying TTree object. More...
 
const TTree * tree () const override final
 
std::vector< DataDependencydata_dependencies () override final
 Returns the data dependencies needed by the MuonTesterBranch. More...
 
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...
 

Static Public Member Functions

static std::string systName (const CP::SystematicSet &set)
 Dummy helper function that encodes the systname set. More...
 
static std::string eraseWhiteSpaces (const std::string &In)
 

Protected Member Functions

template<class T >
bool addToTree (T &variable)
 Connects the Memory buffer with the TTree. More...
 
template<class Key >
bool declare_dependency (Key &key)
 Declares the ReadHandle/ ReadCondHandleKey as data dependency of the algorithm. More...
 
MuonTesterTree & parent ()
 Returns the reference to the MuonTesterTree parent. More...
 

Private Member Functions

CP::CorrectionCode fillSystematic (const xAOD::Muon &muon, SFSet &set)
 
void initMessaging () const
 Initialize our message level and MessageSvc. More...
 

Private Attributes

ToolHandle< CP::IMuonEfficiencyScaleFactorsm_handle {}
 
bool m_uncorrelate_sys {false}
 
std::vector< std::unique_ptr< SFSet > > m_SFs {}
 
TTree * m_tree {nullptr}
 
MuonTesterTree * m_parent {nullptr}
 
std::string m_name {}
 
bool m_init {false}
 
std::vector< DataDependencym_dependencies {}
 
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 71 of file MuonSFTestHelper.h.

Member Typedef Documentation

◆ DataDependency

Definition at line 39 of file IMuonTesterBranch.h.

Constructor & Destructor Documentation

◆ MuonSFBranches()

TestMuonSF::MuonSFBranches::MuonSFBranches ( MuonVal::MuonTesterTree tree,
const ToolHandle< CP::IMuonEfficiencyScaleFactors > &  handle,
const std::string &  rel_name = "" 
)

Definition at line 90 of file MuonSFTestHelper.cxx.

92  :
93  MuonVal::MuonTesterBranch{tree, relName + (relName.empty() ? "" : "_") +
94  getProperty<std::string>(handle.get(), "WorkingPoint")},
95  m_handle{handle} {}

◆ ~MuonSFBranches()

virtual TestMuonSF::MuonSFBranches::~MuonSFBranches ( )
virtualdefault

Member Function Documentation

◆ addToTree()

template<class T >
bool MuonVal::MuonTesterBranch::addToTree ( T &  variable)
protectedinherited

Connects the Memory buffer with the TTree.

◆ data_dependencies()

std::vector< MuonTesterBranch::DataDependency > MuonVal::MuonTesterBranch::data_dependencies ( )
finaloverridevirtualinherited

Returns the data dependencies needed by the MuonTesterBranch.

Implements MuonVal::IMuonTesterBranch.

Definition at line 61 of file MuonTesterBranch.cxx.

61 { return m_dependencies;}

◆ declare_dependency()

template<class Key >
bool MuonVal::MuonTesterBranch::declare_dependency ( Key &  key)
protectedinherited

Declares the ReadHandle/ ReadCondHandleKey as data dependency of the algorithm.

◆ eraseWhiteSpaces()

std::string MuonVal::MuonTesterBranch::eraseWhiteSpaces ( const std::string &  In)
staticinherited

Definition at line 56 of file MuonTesterBranch.cxx.

56  {
57  std::string out = In;
58  out.erase(std::remove_if(out.begin(), out.end(), isspace), out.end());
59  return out;
60 }

◆ fill()

bool TestMuonSF::MuonSFBranches::fill ( const EventContext &  ctx)
finaloverridevirtual

The fill method checks if enough information is provided such that the branch is cleared from the information of the prevoius event and in cases of that updates did no happen befure a default value is provided.

--> returns false if no update is called or the dummy value has not been defined

Implements MuonVal::IMuonTesterBranch.

Definition at line 112 of file MuonSFTestHelper.cxx.

112  {
113  for (auto& sys : m_SFs) {
114  if (!sys->fill(ctx)) return false;
115  }
116  return true;
117  }

◆ fillSystematic()

CP::CorrectionCode TestMuonSF::MuonSFBranches::fillSystematic ( const xAOD::Muon muon,
SFSet set 
)
private

No data-mc efficiencies provided for eta beyond 2.5

Definition at line 154 of file MuonSFTestHelper.cxx.

155  {
156  if (m_handle->applySystematicVariation(Syst_SF.cpSet) != StatusCode::SUCCESS) {
157  ATH_MSG_ERROR("Failed to apply variation "<< Syst_SF.cpSet.name()<<" for "<< name());
159  }
160  float sf{1.f};
161  CP::CorrectionCode cc = m_handle->getEfficiencyScaleFactor(muon, sf);
162  if (cc == CP::CorrectionCode::Error) {
163  ATH_MSG_ERROR("Failed to retrieve the scale factor variation "<< Syst_SF.cpSet.name()<<" for "<< name());
165  }
166  Syst_SF.scaleFactor = sf;
167 
169  if (std::abs(muon.eta()) > 2.5) {
170  Syst_SF.dataEff = -1;
171  Syst_SF.mcEff = -1;
172  return CP::CorrectionCode::Ok;
173  }
174  cc = m_handle->getDataEfficiency(muon, sf);
175  if (cc == CP::CorrectionCode::Error) {
176  ATH_MSG_ERROR("Failed to retrieve the data efficiency variation "<< Syst_SF.cpSet.name()<<" for "<< name());
177  return cc;
178  }
179  Syst_SF.dataEff = sf;
180  cc = m_handle->getMCEfficiency(muon, sf);
181  if (cc == CP::CorrectionCode::Error) {
182  ATH_MSG_ERROR("Failed to retrieve the mc efficiency variation "<< Syst_SF.cpSet.name()<<" for "<< name());
184  }
185  Syst_SF.mcEff = sf;
186  return CP::CorrectionCode::Ok;
187  }

◆ init()

bool TestMuonSF::MuonSFBranches::init ( )
finaloverridevirtual

The init method checks whether the branch name has already registered to the MuonTree and tries then to add the branch to the tree --> return false if the branch has already been defined or the inclusion of the branch to the tree failed or the name of the branch is empty.

Implements MuonVal::IMuonTesterBranch.

Definition at line 97 of file MuonSFTestHelper.cxx.

97  {
98  auto mesf = dynamic_cast<const CP::MuonEfficiencyScaleFactors*>(m_handle.get());
99  if (!mesf) {
100  ATH_MSG_ERROR("Given tool not MuonCP::MuonEfficiencyScaleFactors");
101  return false;
102  }
103  m_uncorrelate_sys = mesf->uncorrelate_sys();
104 
105  for (const auto& set : CP::make_systematics_vector(m_handle->recommendedSystematics())) {
106  m_SFs.emplace_back(std::make_unique<SFSet>(set, *this));
107  if (!m_SFs.back()->init()) return false;
108  }
109  return true;
110  }

◆ initialized()

bool MuonVal::MuonTesterBranch::initialized ( ) const
inherited

Returns whether the object has been initialized or not The outcome of the variable is bound to a successful call of the addToTree method.

Definition at line 52 of file MuonTesterBranch.cxx.

52 { return m_init; }

◆ 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.

40 {
42  m_lvl = m_imsg ?
43  static_cast<MSG::Level>( m_imsg.load()->outputLevel(m_nm) ) :
44  MSG::INFO;
45 }

◆ 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.

165 {
166  MsgStream* ms = m_msg_tls.get();
167  if (!ms) {
168  if (!m_initialized.test_and_set()) initMessaging();
169  ms = new MsgStream(m_imsg,m_nm);
170  m_msg_tls.reset( ms );
171  }
172 
173  ms->setLevel (m_lvl);
174  return *ms;
175 }

◆ 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.

180 { return msg() << lvl; }

◆ msgLvl()

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

Test the output level.

Parameters
lvlThe message level to test against
Returns
boolean Indicating if messages at given level will be printed
Return values
trueMessages at level "lvl" will be printed

Definition at line 151 of file AthMessaging.h.

152 {
153  if (!m_initialized.test_and_set()) initMessaging();
154  if (m_lvl <= lvl) {
155  msg() << lvl;
156  return true;
157  } else {
158  return false;
159  }
160 }

◆ name()

std::string MuonVal::MuonTesterBranch::name ( ) const
finaloverridevirtualinherited

Returns the name of the branch.

Implements MuonVal::IMuonTesterBranch.

Definition at line 51 of file MuonTesterBranch.cxx.

51 { return m_name; }

◆ parent()

MuonTesterTree & MuonVal::MuonTesterBranch::parent ( )
protectedinherited

Returns the reference to the MuonTesterTree parent.

Definition at line 38 of file MuonTesterBranch.cxx.

38  {
39  if (!m_parent) {
40  ATH_MSG_WARNING("The parent of " << name() << " is null.");
41  }
42  return *m_parent;
43 }

◆ 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.

29 {
30  m_lvl = lvl;
31 }

◆ setMuon()

void TestMuonSF::MuonSFBranches::setMuon ( const xAOD::Muon muon)
finaloverridevirtual

Only the raw systematic sets have been activated We can loop over each set and fill it properly

Only process the nominal set or the set matching the syst bin

Implements TestMuonSF::MuonEffiBranch.

Definition at line 120 of file MuonSFTestHelper.cxx.

120  {
123  if (!m_uncorrelate_sys) {
124  for (auto& Syst_SF : m_SFs) {
127  return;
128  }
129  }
130  } else {
131  int bin = m_handle->getUnCorrelatedSystBin(muon);
132  if (bin < 0){
133  ATH_MSG_WARNING("Did not find a valid bin for muon with pT: "<<muon.pt() *MeVtoGeV <<" GeV, eta:"<<
134  muon.eta()<<", phi: "<<muon.phi());
135  return;
136  }
137  for (auto& Syst_SF: m_SFs){
139  bool process = Syst_SF->cpSet.name().empty() ||
140  std::find_if(Syst_SF->cpSet.begin(), Syst_SF->cpSet.end(),
141  [bin](const CP::SystematicVariation& t){
142  return t.isToyVariation() && t.getToyVariation().first ==
143  static_cast<unsigned>(bin);
144  }) != Syst_SF->cpSet.end();
145  if (!process) continue;
147  if (cc == CP::CorrectionCode::Error) {
148  return;
149  }
150  }
151  }
152  }

◆ systName()

static std::string TestMuonSF::MuonSFBranches::systName ( const CP::SystematicSet set)
inlinestatic

Dummy helper function that encodes the systname set.

Definition at line 87 of file MuonSFTestHelper.h.

87  {
88  return set.name().empty() ? std::string("") : std::string("__") + set.name();
89  }

◆ tree() [1/2]

const TTree * MuonVal::MuonTesterBranch::tree ( ) const
finaloverridevirtualinherited

Implements MuonVal::IMuonTesterBranch.

Definition at line 54 of file MuonTesterBranch.cxx.

54 { return m_tree; }

◆ tree() [2/2]

TTree * MuonVal::MuonTesterBranch::tree ( )
finaloverridevirtualinherited

Returns the underlying TTree object.

Implements MuonVal::IMuonTesterBranch.

Definition at line 53 of file MuonTesterBranch.cxx.

53 { return m_tree; }

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_dependencies

std::vector<DataDependency> MuonVal::MuonTesterBranch::m_dependencies {}
privateinherited

Definition at line 62 of file MuonTesterBranch.h.

◆ m_handle

ToolHandle<CP::IMuonEfficiencyScaleFactors> TestMuonSF::MuonSFBranches::m_handle {}
private

Definition at line 92 of file MuonSFTestHelper.h.

◆ m_imsg

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

MessageSvc pointer.

Definition at line 135 of file AthMessaging.h.

◆ m_init

bool MuonVal::MuonTesterBranch::m_init {false}
privateinherited

Definition at line 61 of file MuonTesterBranch.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_name

std::string MuonVal::MuonTesterBranch::m_name {}
privateinherited

Definition at line 60 of file MuonTesterBranch.h.

◆ m_nm

std::string AthMessaging::m_nm
privateinherited

Message source name.

Definition at line 129 of file AthMessaging.h.

◆ m_parent

MuonTesterTree* MuonVal::MuonTesterBranch::m_parent {nullptr}
privateinherited

Definition at line 59 of file MuonTesterBranch.h.

◆ m_SFs

std::vector<std::unique_ptr<SFSet> > TestMuonSF::MuonSFBranches::m_SFs {}
private

Definition at line 121 of file MuonSFTestHelper.h.

◆ m_tree

TTree* MuonVal::MuonTesterBranch::m_tree {nullptr}
privateinherited

Definition at line 58 of file MuonTesterBranch.h.

◆ m_uncorrelate_sys

bool TestMuonSF::MuonSFBranches::m_uncorrelate_sys {false}
private

Definition at line 93 of file MuonSFTestHelper.h.


The documentation for this class was generated from the following files:
AthMessaging::m_lvl
std::atomic< MSG::Level > m_lvl
Current logging level.
Definition: AthMessaging.h:138
xAOD::muon
@ muon
Definition: TrackingPrimitives.h:195
CP::make_systematics_vector
std::vector< CP::SystematicSet > make_systematics_vector(const SystematicSet &systematics)
utility functions for working with systematics
Definition: SystematicsUtil.cxx:25
TestMuonSF::MuonSFBranches::m_handle
ToolHandle< CP::IMuonEfficiencyScaleFactors > m_handle
Definition: MuonSFTestHelper.h:92
TestMuonSF::MuonSFBranches::m_uncorrelate_sys
bool m_uncorrelate_sys
Definition: MuonSFTestHelper.h:93
MuonVal::MuonTesterBranch
Definition: MuonTesterBranch.h:21
MuonVal::MuonTesterBranch::m_dependencies
std::vector< DataDependency > m_dependencies
Definition: MuonTesterBranch.h:62
MuonVal::MuonTesterBranch::m_tree
TTree * m_tree
Definition: MuonTesterBranch.h:58
python.AthDsoLogger.out
out
Definition: AthDsoLogger.py:71
CP::MuonEfficiencyScaleFactors
Definition: MuonEfficiencyScaleFactors.h:23
CP::MeVtoGeV
constexpr float MeVtoGeV
Definition: IsolationCloseByCorrectionTool.cxx:33
bin
Definition: BinsDiffFromStripMedian.h:43
read_hist_ntuple.t
t
Definition: read_hist_ntuple.py:5
CP::SystematicVariation
Definition: SystematicVariation.h:47
AthMessaging::m_imsg
std::atomic< IMessageSvc * > m_imsg
MessageSvc pointer.
Definition: AthMessaging.h:135
python.SystemOfUnits.ms
int ms
Definition: SystemOfUnits.py:132
SUSY_SimplifiedModel_PostInclude.process
string process
Definition: SUSY_SimplifiedModel_PostInclude.py:42
Athena::getMessageSvc
IMessageSvc * getMessageSvc(bool quiet=false)
Definition: getMessageSvc.cxx:20
mapkey::sys
@ sys
Definition: TElectronEfficiencyCorrectionTool.cxx:42
MuonVal::MuonTesterBranch::tree
TTree * tree() override final
Returns the underlying TTree object.
Definition: MuonTesterBranch.cxx:53
MuonVal::MuonTesterBranch::m_parent
MuonTesterTree * m_parent
Definition: MuonTesterBranch.h:59
CP::CorrectionCode::Error
@ Error
Some error happened during the object correction.
Definition: CorrectionCode.h:36
TrigConf::MSGTC::Level
Level
Definition: Trigger/TrigConfiguration/TrigConfBase/TrigConfBase/MsgStream.h:21
ATH_MSG_ERROR
#define ATH_MSG_ERROR(x)
Definition: AthMsgStreamMacros.h:33
MuonVal::MuonTesterBranch::m_init
bool m_init
Definition: MuonTesterBranch.h:61
AthMessaging::msg
MsgStream & msg() const
The standard message stream.
Definition: AthMessaging.h:164
MuonVal::MuonTesterBranch::name
std::string name() const override final
Returns the name of the branch.
Definition: MuonTesterBranch.cxx:51
CxxUtils::set
constexpr std::enable_if_t< is_bitmask_v< E >, E & > set(E &lhs, E rhs)
Convenience function to set bits in a class enum bitmask.
Definition: bitmask.h:224
TestMuonSF::MuonSFBranches::fillSystematic
CP::CorrectionCode fillSystematic(const xAOD::Muon &muon, SFSet &set)
Definition: MuonSFTestHelper.cxx:154
CP::CorrectionCode::Ok
@ Ok
The correction was done successfully.
Definition: CorrectionCode.h:38
MuonVal::MuonTesterBranch::m_name
std::string m_name
Definition: MuonTesterBranch.h:60
mapkey::sf
@ sf
Definition: TElectronEfficiencyCorrectionTool.cxx:38
ATH_MSG_WARNING
#define ATH_MSG_WARNING(x)
Definition: AthMsgStreamMacros.h:32
TestMuonSF::MuonSFBranches::m_SFs
std::vector< std::unique_ptr< SFSet > > m_SFs
Definition: MuonSFTestHelper.h:121
AthMessaging::m_nm
std::string m_nm
Message source name.
Definition: AthMessaging.h:129
CP::CorrectionCode
Return value from object correction CP tools.
Definition: CorrectionCode.h:31
AthMessaging::initMessaging
void initMessaging() const
Initialize our message level and MessageSvc.
Definition: AthMessaging.cxx:39
AthMessaging::m_msg_tls
boost::thread_specific_ptr< MsgStream > m_msg_tls
MsgStream instance (a std::cout like with print-out levels)
Definition: AthMessaging.h:132
python.handimod.cc
int cc
Definition: handimod.py:523