![Logo](../../ATLAS-Logo-Square-Blue-RGB.png) |
ATLAS Offline Software
|
#include <BTagScaleFactorCalculator.h>
|
| BTagScaleFactorCalculator (const std::string &name) |
|
virtual | ~BTagScaleFactorCalculator () |
|
| BTagScaleFactorCalculator (const BTagScaleFactorCalculator &rhs)=delete |
|
| BTagScaleFactorCalculator (BTagScaleFactorCalculator &&rhs)=delete |
|
BTagScaleFactorCalculator & | operator= (const BTagScaleFactorCalculator &rhs)=delete |
|
StatusCode | initialize () |
| Dummy implementation of the initialisation function. More...
|
|
StatusCode | execute () |
|
StatusCode | apply (const std::shared_ptr< std::unordered_map< std::size_t, std::string > > &jet_syst_collections, bool use_trackjets=false) |
|
StatusCode | debug () |
|
virtual void | print () const |
| Print the state of the tool. More...
|
|
ServiceHandle< StoreGateSvc > & | evtStore () |
| The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc . More...
|
|
const ServiceHandle< StoreGateSvc > & | evtStore () const |
| The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc . More...
|
|
const ServiceHandle< StoreGateSvc > & | detStore () const |
| The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc . More...
|
|
virtual StatusCode | sysInitialize () override |
| Perform system initialization for an algorithm. More...
|
|
virtual StatusCode | sysStart () override |
| Handle START transition. More...
|
|
virtual std::vector< Gaudi::DataHandle * > | inputHandles () const override |
| Return this algorithm's input handles. More...
|
|
virtual std::vector< Gaudi::DataHandle * > | outputHandles () const override |
| Return this algorithm's output handles. More...
|
|
Gaudi::Details::PropertyBase & | declareProperty (Gaudi::Property< T > &t) |
|
Gaudi::Details::PropertyBase * | declareProperty (const std::string &name, SG::VarHandleKey &hndl, const std::string &doc, const SG::VarHandleKeyType &) |
| Declare a new Gaudi property. More...
|
|
Gaudi::Details::PropertyBase * | declareProperty (const std::string &name, SG::VarHandleBase &hndl, const std::string &doc, const SG::VarHandleType &) |
| Declare a new Gaudi property. More...
|
|
Gaudi::Details::PropertyBase * | declareProperty (const std::string &name, SG::VarHandleKeyArray &hndArr, const std::string &doc, const SG::VarHandleKeyArrayType &) |
|
Gaudi::Details::PropertyBase * | declareProperty (const std::string &name, T &property, const std::string &doc, const SG::NotHandleType &) |
| Declare a new Gaudi property. More...
|
|
Gaudi::Details::PropertyBase * | declareProperty (const std::string &name, T &property, const std::string &doc="none") |
| Declare a new Gaudi property. More...
|
|
void | updateVHKA (Gaudi::Details::PropertyBase &) |
|
MsgStream & | msg () const |
|
MsgStream & | msg (const MSG::Level lvl) const |
|
bool | msgLvl (const MSG::Level lvl) const |
|
Definition at line 39 of file BTagScaleFactorCalculator.h.
◆ StoreGateSvc_t
◆ BTagScaleFactorCalculator() [1/3]
top::BTagScaleFactorCalculator::BTagScaleFactorCalculator |
( |
const std::string & |
name | ) |
|
|
explicit |
◆ ~BTagScaleFactorCalculator()
virtual top::BTagScaleFactorCalculator::~BTagScaleFactorCalculator |
( |
| ) |
|
|
inlinevirtual |
◆ BTagScaleFactorCalculator() [2/3]
◆ BTagScaleFactorCalculator() [3/3]
◆ apply()
StatusCode top::BTagScaleFactorCalculator::apply |
( |
const std::shared_ptr< std::unordered_map< std::size_t, std::string > > & |
jet_syst_collections, |
|
|
bool |
use_trackjets = false |
|
) |
| |
– Loop over all jet collections –/// – Lets assume that we're not doing ElectronInJet subtraction –///
– Tell the SF tools to use the nominal systematic –/// – Loop over all jets in each collection –///
– For nominal calibration, vary the SF systematics –///
Definition at line 69 of file BTagScaleFactorCalculator.cxx.
74 for (
auto currentSystematic : *jet_syst_collections) {
80 for (
auto jetPtr : *
jets) {
81 bool passSelection(
false);
82 if (jetPtr->isAvailable<
char>(
"passPreORSelection")) {
83 if (jetPtr->auxdataConst<
char>(
"passPreORSelection") == 1) {
87 if (jetPtr->isAvailable<
char>(
"passPreORSelectionLoose")) {
88 if (jetPtr->auxdataConst<
char>(
"passPreORSelectionLoose") == 1) {
96 for (
const std::string& tagWP : (use_trackjets ?
m_config->bTagWP_calib_trkJet() :
m_config->bTagWP_calib())) {
97 ToolHandle<IBTaggingEfficiencyTool>& btageff =
99 ToolHandle<IBTaggingSelectionTool>& btagsel =
104 unsigned int MapIndex =
m_config->getMapIndex();
106 btageff->setMapIndex(
"Light", MapIndex);
107 btageff->setMapIndex(
"C", MapIndex);
108 btageff->setMapIndex(
"B", MapIndex);
109 btageff->setMapIndex(
"T", MapIndex);
115 if (bTagSystName !=
"") {
118 top::check(btageff->applySystematicVariation(bTagSyst),
119 "Failed to set new b-tagging SF to a shifted systematic set : " + bTagSystName);
122 "Failed to set new b-tagging SF to nominal");
126 bool isTagged =
false;
127 if (std::fabs(jetPtr->eta()) <= 2.5) {
128 if (tagWP.find(
"Continuous") == std::string::npos) {
129 isTagged =
static_cast<bool>(btagsel->accept(*jetPtr));
130 if (isTagged)
top::check(btageff->getScaleFactor(*jetPtr, btag_SF),
131 "Failed to get nominal b-tagging SF");
132 else top::check(btageff->getInefficiencyScaleFactor(*jetPtr, btag_SF),
133 "Failed to get nominal b-tagging SF");
135 top::check(btageff->getScaleFactor(*jetPtr, btag_SF),
136 "Failed to get nominal Continuous b-tagging SF");
139 jetPtr->auxdecor<
float>(
"btag_SF_" + tagWP +
"_nom") = btag_SF;
142 if (currentSystematic.first ==
m_config->nominalHashValue()) {
143 for (
const auto& variation : sysSet) {
146 syst_set.
insert(variation);
147 top::check(btageff->applySystematicVariation(syst_set),
148 "Failed to set new b-tagging systematic variation " + syst_set.
name());
149 if (std::fabs(jetPtr->eta()) <= 2.5) {
150 if (tagWP.find(
"Continuous") == std::string::npos) {
151 if (isTagged)
top::check(btageff->getScaleFactor(*jetPtr, btag_SF),
152 "Failed to get b-tagging SF for variation " + syst_set.
name());
153 else top::check(btageff->getInefficiencyScaleFactor(*jetPtr, btag_SF),
154 "Failed to get b-tagging SF for variation " + syst_set.
name());
156 top::check(btageff->getScaleFactor(*jetPtr, btag_SF),
157 "Failed to get Continuous b-tagging SF for variation " + syst_set.
name());
160 jetPtr->auxdecor<
float>(
"btag_SF_" + tagWP +
"_" + variation.name()) = btag_SF;
169 return StatusCode::SUCCESS;
◆ debug()
StatusCode top::BTagScaleFactorCalculator::debug |
( |
| ) |
|
Definition at line 172 of file BTagScaleFactorCalculator.cxx.
173 ATH_MSG_INFO(
"BTagScaleFactorCalculator::debug function");
175 for (
const std::string& tagWP :
m_config->bTagWP()) {
177 ToolHandle<IBTaggingEfficiencyTool>& btageff =
m_btagEffTools[tagWP];
179 top::check(btageff.retrieve(),
"Failed to retrieve tool");
184 ATH_MSG_INFO(
"-----------------------------------------------------------------------");
187 std::vector<std::string> > allowed_variations = btageff->listSystematics();
189 ATH_MSG_INFO(
"Allowed systematics variations for tool " << btageff->name());
190 for (
auto var : allowed_variations) {
191 std::string flvs =
"";
192 for (
auto flv :
var.second) flvs += flv;
198 ATH_MSG_INFO(
"-----------------------------------------------------------------------");
199 ATH_MSG_INFO(
"List of b-tagging scale factor systematics");
200 std::map<std::string,
201 std::vector<std::string> > listOfScaleFactorSystematics = btageff->listScaleFactorSystematics(
false);
202 for (
auto var : listOfScaleFactorSystematics) {
204 std::vector<std::string> systs =
var.second;
206 for (
const auto&
sys : systs) {
212 ATH_MSG_INFO(
"List of (named) b-tagging scale factor systematics");
213 listOfScaleFactorSystematics = btageff->listScaleFactorSystematics(
true);
214 for (
auto var : listOfScaleFactorSystematics) {
216 std::vector<std::string> systs =
var.second;
218 for (
const auto&
sys : systs) {
223 ATH_MSG_INFO(
"-----------------------------------------------------------------------");
225 return StatusCode::SUCCESS;
◆ declareGaudiProperty() [1/4]
specialization for handling Gaudi::Property<SG::VarHandleKeyArray>
Definition at line 170 of file AthCommonDataStore.h.
175 hndl.documentation());
◆ declareGaudiProperty() [2/4]
specialization for handling Gaudi::Property<SG::VarHandleKey>
Definition at line 156 of file AthCommonDataStore.h.
161 hndl.documentation());
◆ declareGaudiProperty() [3/4]
specialization for handling Gaudi::Property<SG::VarHandleBase>
Definition at line 184 of file AthCommonDataStore.h.
189 hndl.documentation());
◆ declareGaudiProperty() [4/4]
◆ declareProperty() [1/6]
Declare a new Gaudi property.
- Parameters
-
name | Name of the property. |
hndl | Object holding the property value. |
doc | Documentation string for the property. |
This is the version for types that derive from SG::VarHandleBase
. The property value object is put on the input and output lists as appropriate; then we forward to the base class.
Definition at line 245 of file AthCommonDataStore.h.
250 this->declare(hndl.
vhKey());
251 hndl.
vhKey().setOwner(
this);
◆ declareProperty() [2/6]
Declare a new Gaudi property.
- Parameters
-
name | Name of the property. |
hndl | Object holding the property value. |
doc | Documentation string for the property. |
This is the version for types that derive from SG::VarHandleKey
. The property value object is put on the input and output lists as appropriate; then we forward to the base class.
Definition at line 221 of file AthCommonDataStore.h.
◆ declareProperty() [3/6]
◆ declareProperty() [4/6]
Declare a new Gaudi property.
- Parameters
-
name | Name of the property. |
property | Object holding the property value. |
doc | Documentation string for the property. |
This is the generic version, for types that do not derive from SG::VarHandleKey
. It just forwards to the base class version of declareProperty
.
Definition at line 333 of file AthCommonDataStore.h.
◆ declareProperty() [5/6]
Declare a new Gaudi property.
- Parameters
-
name | Name of the property. |
property | Object holding the property value. |
doc | Documentation string for the property. |
This dispatches to either the generic declareProperty
or the one for VarHandle/Key/KeyArray.
Definition at line 352 of file AthCommonDataStore.h.
◆ declareProperty() [6/6]
◆ detStore()
◆ evtStore() [1/2]
◆ evtStore() [2/2]
◆ execute()
StatusCode top::BTagScaleFactorCalculator::execute |
( |
| ) |
|
Definition at line 60 of file BTagScaleFactorCalculator.cxx.
62 "Failed to apply btagging SFs");
64 "Failed to apply track jet btagging SFs");
66 return StatusCode::SUCCESS;
◆ extraDeps_update_handler()
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
◆ getKey()
Get the (hashed) key of an object that is in the event store.
This is a bit of a special one. StoreGateSvc
and xAOD::TEvent
both provide ways for getting the SG::sgkey_t
key for an object that is in the store, based on a bare pointer. But they provide different interfaces for doing so.
In order to allow tools to efficiently perform this operation, they can use this helper function.
- See also
- asg::AsgTool::getName
- Parameters
-
ptr | The bare pointer to the object that the event store should know about |
- Returns
- The hashed key of the object in the store. If not found, an invalid (zero) key.
Definition at line 119 of file AsgTool.cxx.
121 #ifdef XAOD_STANDALONE
124 return evtStore()->event()->getKey( ptr );
127 return (
proxy ==
nullptr ? 0 :
proxy->sgkey() );
128 #endif // XAOD_STANDALONE
◆ getName()
const std::string & asg::AsgTool::getName |
( |
const void * |
ptr | ) |
const |
|
inherited |
Get the name of an object that is / should be in the event store.
This is a bit of a special one. StoreGateSvc
and xAOD::TEvent
both provide ways for getting the std::string
name for an object that is in the store, based on a bare pointer. But they provide different interfaces for doing so.
In order to allow tools to efficiently perform this operation, they can use this helper function.
- See also
- asg::AsgTool::getKey
- Parameters
-
ptr | The bare pointer to the object that the event store should know about |
- Returns
- The string name of the object in the store. If not found, an empty string.
Definition at line 106 of file AsgTool.cxx.
108 #ifdef XAOD_STANDALONE
111 return evtStore()->event()->getName( ptr );
114 static const std::string
dummy =
"";
116 #endif // XAOD_STANDALONE
◆ getProperty()
template<class T >
const T* asg::AsgTool::getProperty |
( |
const std::string & |
name | ) |
const |
|
inherited |
Get one of the tool's properties.
◆ initialize()
StatusCode top::BTagScaleFactorCalculator::initialize |
( |
| ) |
|
|
virtual |
Dummy implementation of the initialisation function.
It's here to allow the dual-use tools to skip defining an initialisation function. Since many are doing so...
Reimplemented from asg::AsgTool.
Definition at line 31 of file BTagScaleFactorCalculator.cxx.
32 ATH_MSG_INFO(
" top::BTagScaleFactorCalculator initialize");
35 for (
const std::string&
WP :
m_config->bTagWP_calib()) {
41 std::set<std::string> base_names =
m_systs[
WP].getBaseNames();
46 for (
const std::string&
WP :
m_config->bTagWP_calib_trkJet()) {
53 m_config->setBTaggingSFSysts(
WP, base_names,
true);
57 return StatusCode::SUCCESS;
◆ inputHandles()
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.
◆ msg() [1/2]
◆ msg() [2/2]
◆ msg_level_name()
const std::string & asg::AsgTool::msg_level_name |
( |
| ) |
const |
|
inherited |
A deprecated function for getting the message level's name.
Instead of using this, weirdly named function, user code should get the string name of the current minimum message level (in case they really need it...), with:
MSG::name( msg().level() )
This function's name doesn't follow the ATLAS coding rules, and as such will be removed in the not too distant future.
- Returns
- The string name of the current minimum message level that's printed
Definition at line 101 of file AsgTool.cxx.
◆ msgLvl()
◆ operator=()
◆ outputHandles()
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.
◆ print()
void asg::AsgTool::print |
( |
| ) |
const |
|
virtualinherited |
Print the state of the tool.
Implements asg::IAsgTool.
Reimplemented in JetRecTool, JetFinder, JetModifiedMassDrop, JetFromPseudojet, JetReclusterer, JetReclusteringTool, JetTruthLabelingTool, JetPileupLabelingTool, HI::HIPileupTool, top::TopObjectSelection, JetDumper, JetBottomUpSoftDrop, JetRecursiveSoftDrop, JetSoftDrop, JetConstituentsRetriever, JetSubStructureMomentToolsBase, JetSplitter, JetToolRunner, JetPruner, JetPseudojetRetriever, JetTrimmer, AsgHelloTool, and KtDeltaRTool.
Definition at line 131 of file AsgTool.cxx.
◆ renounce()
◆ renounceArray()
◆ sysInitialize()
◆ sysStart()
Handle START transition.
We override this in order to make sure that conditions handle keys can cache a pointer to the conditions container.
◆ updateVHKA()
◆ m_btagEffTools
std::unordered_map<std::string, ToolHandle<IBTaggingEfficiencyTool> > top::BTagScaleFactorCalculator::m_btagEffTools |
|
private |
◆ m_btagSelTools
std::unordered_map<std::string, ToolHandle<IBTaggingSelectionTool> > top::BTagScaleFactorCalculator::m_btagSelTools |
|
private |
◆ m_config
std::shared_ptr<top::TopConfig> top::BTagScaleFactorCalculator::m_config |
|
private |
◆ m_detStore
◆ m_evtStore
◆ m_nominal
◆ m_systs
std::unordered_map<std::string, CP::SystematicSet> top::BTagScaleFactorCalculator::m_systs |
|
private |
◆ m_trkjet_btagEffTools
std::unordered_map<std::string, ToolHandle<IBTaggingEfficiencyTool> > top::BTagScaleFactorCalculator::m_trkjet_btagEffTools |
|
private |
◆ m_trkjet_btagSelTools
std::unordered_map<std::string, ToolHandle<IBTaggingSelectionTool> > top::BTagScaleFactorCalculator::m_trkjet_btagSelTools |
|
private |
◆ m_trkjet_systs
std::unordered_map<std::string, CP::SystematicSet> top::BTagScaleFactorCalculator::m_trkjet_systs |
|
private |
◆ m_varHandleArraysDeclared
◆ m_vhka
The documentation for this class was generated from the following files:
def retrieve(aClass, aKey=None)
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T > &t)
Class to wrap a set of SystematicVariations.
StoreGateSvc_t m_evtStore
Pointer to StoreGate (event store by default)
std::vector< SG::VarHandleKeyArray * > m_vhka
std::string name() const
returns: the systematics joined into a single string.
virtual void setOwner(IDataHandleHolder *o)=0
ServiceHandle< StoreGateSvc > & evtStore()
The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.
SystematicVariation getSystematicByBaseName(const std::string &basename) const
description: get the first systematic matching basename
void check(bool thingToCheck, const std::string &usefulFailureMessage)
Print an error message and terminate if thingToCheck is false.
std::string bTagNamedSystCheck(std::shared_ptr< top::TopConfig > config, const std::string &systCollection, const std::string &tagWP, const bool isTrackJet, const bool verbose=false)
const std::string & name(Level lvl)
Convenience function for translating message levels to strings.
StoreGateSvc_t m_detStore
Pointer to StoreGate (detector store by default)
virtual void renounce()=0
std::conditional< std::is_base_of< SG::VarHandleKeyArray, T >::value, VarHandleKeyArrayType, type2 >::type type
void insert(const SystematicVariation &systematic)
description: insert a systematic into the set
void sort(typename std::reverse_iterator< DataModel_detail::iterator< DVL > > beg, typename std::reverse_iterator< DataModel_detail::iterator< DVL > > end, const Compare &comp)
Specialization of sort for DataVector/List.
SG::VarHandleKey & vhKey()
Return a non-const reference to the HandleKey.
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>