![Logo](../../ATLAS-Logo-Square-Blue-RGB.png) |
ATLAS Offline Software
|
#include <CommonDiTauEfficiencyTool.h>
|
| CommonDiTauEfficiencyTool (const std::string &sName) |
| Create a proper constructor for Athena. More...
|
|
| ~CommonDiTauEfficiencyTool () |
|
virtual StatusCode | initialize () override |
| Dummy implementation of the initialisation function. More...
|
|
virtual CP::CorrectionCode | getEfficiencyScaleFactor (const xAOD::DiTauJet &xDiTau, double &dEfficiencyScaleFactor, unsigned int iRunNumber=0, unsigned int iMu=0) override |
| Get the Efficiency Scale Factor of ditau jet. More...
|
|
virtual CP::CorrectionCode | applyEfficiencyScaleFactor (const xAOD::DiTauJet &xDiTau, unsigned int iRunNumber=0, unsigned int iMu=0) override |
| Get the Efficiency Scale Factor of ditau jet. More...
|
|
void | ReadInputs (std::unique_ptr< TFile > &fFile) |
|
CP::CorrectionCode | getValue (const std::string &sHistName, const xAOD::DiTauJet &xDiTau, double &dEfficiencyScaleFactor) const |
| Get the scale factor from a particular recommendations histogram. More...
|
|
void | generateSystematicSets () |
| generate a set of relevant systematic variations to be applied More...
|
|
virtual CP::CorrectionCode | getEfficiencyScaleFactor (const xAOD::TauJet &tau, double &dEfficiencyScaleFactor, unsigned int iRunNumber=0, unsigned int iMu=0) |
|
virtual CP::CorrectionCode | applyEfficiencyScaleFactor (const xAOD::TauJet &xTau, unsigned int iRunNumber=0, unsigned int iMu=0) |
|
virtual CP::CorrectionCode | getValue (const std::string &sHistName, const xAOD::TauJet &xTau, double &dEfficiencyScaleFactor) const |
|
virtual CP::CorrectionCode | getEfficiencyScaleFactor (const xAOD::TauJet &tau, double &dEfficiencyScaleFactor, unsigned int iRunNumber=0, unsigned int iMu=0) |
| Declare the interface that the class provides. More...
|
|
virtual CP::CorrectionCode | applyEfficiencyScaleFactor (const xAOD::TauJet &xTau, unsigned int iRunNumber=0, unsigned int iMu=0) |
| Decorate the tau with its efficiency. More...
|
|
virtual bool | isAffectedBySystematic (const CP::SystematicVariation &systematic) const |
| returns: whether this tool is affected by the given systematics More...
|
|
virtual CP::SystematicSet | affectingSystematics () const |
| returns: the list of all systematics this tool can be affected by More...
|
|
virtual CP::SystematicSet | recommendedSystematics () const |
| returns: the list of all systematics this tool recommends to use More...
|
|
virtual StatusCode | applySystematicVariation (const CP::SystematicSet &sSystematicSet) |
| configure this tool for the given list of systematic variations. More...
|
|
virtual bool | isSupportedRunNumber (int) const |
| check if run number is supported in recommendations More...
|
|
virtual void | print () const =0 |
| Print the state of the tool. More...
|
|
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 CommonDiTauEfficiencyTool.h.
◆ StoreGateSvc_t
◆ tSFMAP
◆ tTupleObjectFunc
◆ CommonDiTauEfficiencyTool()
CommonDiTauEfficiencyTool::CommonDiTauEfficiencyTool |
( |
const std::string & |
sName | ) |
|
◆ ~CommonDiTauEfficiencyTool()
CommonDiTauEfficiencyTool::~CommonDiTauEfficiencyTool |
( |
| ) |
|
◆ addHistogramToSFMap()
void CommonEfficiencyTool::addHistogramToSFMap |
( |
TKey * |
kKey, |
|
|
const std::string & |
sKeyName |
|
) |
| |
|
protectedinherited |
Definition at line 581 of file CommonEfficiencyTool.cxx.
585 TClass *cClass = gROOT->GetClass(kKey->GetClassName());
586 if (cClass->InheritsFrom(
"TH2"))
588 TH1* oObject = (TH1*)kKey->ReadObj();
589 oObject->SetDirectory(0);
593 else if (cClass->InheritsFrom(
"TH3"))
595 TH1* oObject = (TH1*)kKey->ReadObj();
596 oObject->SetDirectory(0);
599 }
else if (cClass->InheritsFrom(
"TH1"))
601 TH1* oObject = (TH1*)kKey->ReadObj();
602 oObject->SetDirectory(0);
606 else if (cClass->InheritsFrom(
"TF1"))
608 TObject* oObject = kKey->ReadObj();
◆ affectingSystematics()
◆ applyEfficiencyScaleFactor() [1/3]
◆ applyEfficiencyScaleFactor() [2/3]
Definition at line 57 of file CommonEfficiencyTool.cxx.
291 ATH_MSG_DEBUG(
m_sVarName <<
" decoration is available on first tau processed, switched off applyEfficiencyScaleFactor for further taus.");
292 ATH_MSG_DEBUG(
"If an application of efficiency scale factors needs to be redone, please pass a shallow copy of the original tau.");
303 return tmpCorrectionCode;
◆ applyEfficiencyScaleFactor() [3/3]
◆ applySystematicVariation()
configure this tool for the given list of systematic variations.
any requested systematics that are not affecting this tool will be silently ignored
Implements CP::ISystematicsTool.
Reimplemented in TauAnalysisTools::TauEfficiencyTriggerTool.
Definition at line 346 of file CommonEfficiencyTool.cxx.
354 return StatusCode::SUCCESS;
358 double dDirection = 0.;
360 for (
auto sSyst : sSystematicSet)
366 ATH_MSG_VERBOSE(
"unsupported systematic variation: "<< sSyst.basename()<<
"; skipping this one");
371 if (sSyst.parameter() * dDirection < 0)
373 ATH_MSG_ERROR(
"unsupported set of systematic variations, you should either use only \"UP\" or only \"DOWN\" systematics in one set!");
375 return StatusCode::FAILURE;
377 dDirection = sSyst.parameter();
381 ATH_MSG_ERROR(
"unsupported set of systematic variations, you should not combine \"TAUS_{TRUE|FAKE}_EFF_*_TOTAL\" with other systematic variations!");
383 return StatusCode::FAILURE;
387 sSystematicSetAvailable.
insert(sSyst);
393 return StatusCode::SUCCESS;
◆ ConvertDecayModeToString()
std::string CommonEfficiencyTool::ConvertDecayModeToString |
( |
const int |
iDecayMode | ) |
const |
|
protectedinherited |
◆ ConvertMuToString()
std::string CommonEfficiencyTool::ConvertMuToString |
( |
const int |
iMu | ) |
const |
|
protectedinherited |
◆ ConvertProngToString()
std::string CommonEfficiencyTool::ConvertProngToString |
( |
const int |
iProngness | ) |
const |
|
protectedinherited |
◆ 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]
◆ 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
◆ generateSystematicSets()
void CommonDiTauEfficiencyTool::generateSystematicSets |
( |
| ) |
|
generate a set of relevant systematic variations to be applied
Definition at line 244 of file CommonDiTauEfficiencyTool.cxx.
247 std::vector<std::string> vSplitInputFilePath = {};
249 std::string sEfficiencyType = vSplitInputFilePath.at(0);
250 std::string sTruthType = vSplitInputFilePath.at(1);
251 std::transform(sEfficiencyType.begin(), sEfficiencyType.end(), sEfficiencyType.begin(), toupper);
252 std::transform(sTruthType.begin(), sTruthType.end(), sTruthType.begin(), toupper);
253 std::string sSystematicBaseString =
"TAUS_"+sTruthType+
"_EFF_"+sEfficiencyType+
"_";
261 for (
auto mSF : *
m_mSF)
264 std::vector<std::string> vSplitNP = {};
265 split(mSF.first,
'_',vSplitNP);
266 std::string sNP = vSplitNP.at(0);
267 std::string sNPUppercase = vSplitNP.at(0);
269 if (sNP ==
"sf")
continue;
271 bool bIsRecommended =
false;
272 if (isupper(sNP.at(0)))
273 bIsRecommended =
true;
275 std::transform(sNPUppercase.begin(), sNPUppercase.end(), sNPUppercase.begin(), toupper);
276 std::string sSystematicString = sSystematicBaseString+sNPUppercase;
286 ATH_MSG_DEBUG(
"connected base name " << sNP <<
" with systematic " <<sSystematicString);
◆ getEfficiencyScaleFactor() [1/3]
CP::CorrectionCode CommonDiTauEfficiencyTool::getEfficiencyScaleFactor |
( |
const xAOD::DiTauJet & |
xDiTau, |
|
|
double & |
dEfficiencyScaleFactor, |
|
|
unsigned int |
iRunNumber = 0 , |
|
|
unsigned int |
iMu = 0 |
|
) |
| |
|
overridevirtual |
Get the Efficiency Scale Factor of ditau jet.
- Parameters
-
xDiTau | : reco DiTauJet |
dEfficiencyScaleFactor | : reference to output variable where efficiency is returned |
iRunNumber | : run number |
iMu | : number of interactions |
- Returns
- CP::CorrectionCode
Implements TauAnalysisTools::IDiTauEfficiencyCorrectionsTool.
Definition at line 86 of file CommonDiTauEfficiencyTool.cxx.
92 dEfficiencyScaleFactor = 1.;
98 dEfficiencyScaleFactor);
101 return tmpCorrectionCode;
108 double dTotalSystematic2 = 0.;
109 double dDirection = 0.;
117 double dUncertaintySyst = 0.;
120 dDirection = syst.parameter();
123 std::string sHistName =
it->second;
124 if (dDirection>0) sHistName+=
"_up";
125 else sHistName+=
"_down";
129 tmpCorrectionCode =
getValue(sHistName,
135 return tmpCorrectionCode;
138 dUncertaintySyst *= dDirection;
141 dTotalSystematic2 += dUncertaintySyst * dUncertaintySyst;
145 dDirection = (dDirection > 0.) ? 1. : -1.;
148 dEfficiencyScaleFactor *= 1. + dDirection * std::sqrt(dTotalSystematic2);
◆ getEfficiencyScaleFactor() [2/3]
Declare the interface that the class provides.
Get the tau efficiency scale factor
Implements TauAnalysisTools::ITauEfficiencyCorrectionsTool.
Reimplemented in TauAnalysisTools::TauEfficiencyTriggerTool.
Definition at line 158 of file CommonEfficiencyTool.cxx.
164 dEfficiencyScaleFactor = 1.;
171 dEfficiencyScaleFactor = 1.;
184 ATH_MSG_WARNING(
"Found tau with unknown decay mode. Skip efficiency correction.");
192 dEfficiencyScaleFactor = 1.;
199 std::string sMu =
"";
200 std::string sMCCampaign =
"";
204 std::string sHistName =
m_sSFHistName + sMode + sMu + sMCCampaign;
209 dEfficiencyScaleFactor);
212 return tmpCorrectionCode;
219 double dTotalSystematic2 = 0.;
220 double dDirection = 0.;
227 double dUncertaintySyst = 0.;
230 dDirection = syst.parameter();
233 sHistName =
it->second;
234 if (dDirection>0.) sHistName+=
"_up";
235 else sHistName+=
"_down";
237 sHistName += sMode + sMu + sMCCampaign;
240 if( (sHistName.find(
"3P") != std::string::npos && sHistName.find(
"1p") != std::string::npos) ||
241 (sHistName.find(
"1P") != std::string::npos && sHistName.find(
"3p") != std::string::npos))
245 tmpCorrectionCode =
getValue(sHistName,
251 return tmpCorrectionCode;
254 dUncertaintySyst *= dDirection;
257 dTotalSystematic2 += dUncertaintySyst * dUncertaintySyst;
261 dDirection = (dDirection > 0.) ? 1. : -1.;
264 dEfficiencyScaleFactor *= 1. + dDirection * std::sqrt(dTotalSystematic2);
◆ getEfficiencyScaleFactor() [3/3]
Definition at line 54 of file CommonEfficiencyTool.cxx.
164 dEfficiencyScaleFactor = 1.;
171 dEfficiencyScaleFactor = 1.;
184 ATH_MSG_WARNING(
"Found tau with unknown decay mode. Skip efficiency correction.");
191 if( xTau.nTracks() != 1 && xTau.nTracks() != 3) {
192 dEfficiencyScaleFactor = 1.;
199 std::string sMu =
"";
200 std::string sMCCampaign =
"";
204 std::string sHistName =
m_sSFHistName + sMode + sMu + sMCCampaign;
209 dEfficiencyScaleFactor);
212 return tmpCorrectionCode;
219 double dTotalSystematic2 = 0.;
220 double dDirection = 0.;
227 double dUncertaintySyst = 0.;
230 dDirection = syst.parameter();
233 sHistName =
it->second;
234 if (dDirection>0.) sHistName+=
"_up";
235 else sHistName+=
"_down";
237 sHistName += sMode + sMu + sMCCampaign;
240 if( (sHistName.find(
"3P") != std::string::npos && sHistName.find(
"1p") != std::string::npos) ||
241 (sHistName.find(
"1P") != std::string::npos && sHistName.find(
"3p") != std::string::npos))
245 tmpCorrectionCode =
getValue(sHistName,
251 return tmpCorrectionCode;
254 dUncertaintySyst *= dDirection;
257 dTotalSystematic2 += dUncertaintySyst * dUncertaintySyst;
261 dDirection = (dDirection > 0.) ? 1. : -1.;
264 dEfficiencyScaleFactor *= 1. + dDirection * std::sqrt(dTotalSystematic2);
◆ 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
◆ GetMcCampaignString()
std::string CommonEfficiencyTool::GetMcCampaignString |
( |
const int |
iMu | ) |
const |
|
protectedinherited |
◆ 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.
◆ getValue() [1/3]
Get the scale factor from a particular recommendations histogram.
- Parameters
-
sHistName | |
xDiTau | |
dEfficiencyScaleFactor | |
- Returns
- CP::CorrectionCode
Definition at line 296 of file CommonDiTauEfficiencyTool.cxx.
301 auto it = mSF.find (sHistName);
304 ATH_MSG_ERROR(
"Object with name "<<sHistName<<
" was not found in input file.");
306 for (
auto eEntry : mSF)
319 double dVars[3] = {dX, dY, dZ};
322 return (std::get<1>(tTuple))(std::get<0>(tTuple), dEfficiencyScaleFactor, dVars);
◆ getValue() [2/3]
Definition at line 110 of file CommonEfficiencyTool.cxx.
700 auto it = mSF.find (sHistName);
703 ATH_MSG_ERROR(
"Object with name "<<sHistName<<
" was not found in input file.");
705 for (
auto eEntry : mSF)
714 double dPt =
m_fX(xTau);
717 double dVars[2] = {dPt,
dEta};
721 return (std::get<1>(tTuple))(std::get<0>(tTuple), dEfficiencyScaleFactor, dVars);
◆ getValue() [3/3]
|
protectedvirtualinherited |
Definition at line 695 of file CommonEfficiencyTool.cxx.
700 auto it = mSF.find (sHistName);
703 ATH_MSG_ERROR(
"Object with name "<<sHistName<<
" was not found in input file.");
705 for (
auto eEntry : mSF)
714 double dPt =
m_fX(xTau);
717 double dVars[2] = {dPt,
dEta};
721 return (std::get<1>(tTuple))(std::get<0>(tTuple), dEfficiencyScaleFactor, dVars);
◆ getValueTF1()
CP::CorrectionCode CommonEfficiencyTool::getValueTF1 |
( |
const TObject * |
oObject, |
|
|
double & |
dEfficiencyScaleFactor, |
|
|
double |
dVars[] |
|
) |
| |
|
staticprotectedinherited |
Definition at line 827 of file CommonEfficiencyTool.cxx.
830 double dPt = dVars[0];
831 double dEta = dVars[1];
833 const TF1* fFunc =
static_cast<const TF1*
>(oObject);
842 dEfficiencyScaleFactor = fFunc->Eval(dPt,
dEta);
◆ getValueTH1()
CP::CorrectionCode CommonEfficiencyTool::getValueTH1 |
( |
const TObject * |
oObject, |
|
|
double & |
dEfficiencyScaleFactor, |
|
|
double |
dVars[] |
|
) |
| |
|
staticprotectedinherited |
Definition at line 730 of file CommonEfficiencyTool.cxx.
733 double dPt = dVars[0];
735 const TH1* hHist =
dynamic_cast<const TH1*
>(oObject);
744 dPt =
std::max(dPt,hHist->GetXaxis()->GetXmin());
746 dPt =
std::min(dPt,hHist->GetXaxis()->GetXmax() * .999);
749 int iBin = hHist->FindFixBin(dPt);
750 dEfficiencyScaleFactor = hHist->GetBinContent(iBin);
◆ getValueTH2()
CP::CorrectionCode CommonEfficiencyTool::getValueTH2 |
( |
const TObject * |
oObject, |
|
|
double & |
dEfficiencyScaleFactor, |
|
|
double |
dVars[] |
|
) |
| |
|
staticprotectedinherited |
Definition at line 760 of file CommonEfficiencyTool.cxx.
763 double dPt = dVars[0];
764 double dEta = dVars[1];
766 const TH2* hHist =
dynamic_cast<const TH2*
>(oObject);
775 dPt =
std::max(dPt,hHist->GetXaxis()->GetXmin());
778 dPt =
std::min(dPt,hHist->GetXaxis()->GetXmax() * .999);
782 int iBin = hHist->FindFixBin(dPt,
dEta);
783 dEfficiencyScaleFactor = hHist->GetBinContent(iBin);
◆ getValueTH3()
CP::CorrectionCode CommonEfficiencyTool::getValueTH3 |
( |
const TObject * |
oObject, |
|
|
double & |
dEfficiencyScaleFactor, |
|
|
double |
dVars[] |
|
) |
| |
|
staticprotectedinherited |
Definition at line 792 of file CommonEfficiencyTool.cxx.
795 double dX = dVars[0];
796 double dY = dVars[1];
797 double dZ = dVars[2];
799 const TH3* hHist =
dynamic_cast<const TH3*
>(oObject);
808 dX =
std::max(dX,hHist->GetXaxis()->GetXmin());
809 dY =
std::max(dY,hHist->GetYaxis()->GetXmin());
810 dZ =
std::max(dZ,hHist->GetZaxis()->GetXmin());
812 dX =
std::min(dX,hHist->GetXaxis()->GetXmax() * .999);
813 dY =
std::min(dY,hHist->GetYaxis()->GetXmax() * .999);
814 dZ =
std::min(dZ,hHist->GetZaxis()->GetXmax() * .999);
817 int iBin = hHist->FindFixBin(dX,dY,dZ);
818 dEfficiencyScaleFactor = hHist->GetBinContent(iBin);
◆ initialize()
StatusCode CommonDiTauEfficiencyTool::initialize |
( |
| ) |
|
|
overridevirtual |
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 TauAnalysisTools::CommonEfficiencyTool.
Definition at line 43 of file CommonDiTauEfficiencyTool.cxx.
45 ATH_MSG_INFO(
"Initializing CommonDiTauEfficiencyTool" );
51 m_mSF = std::make_unique< tSFMAP >();
52 std::unique_ptr< TFile > fSF( TFile::Open( sInputFilePath.c_str(),
"READ" ) );
55 ATH_MSG_FATAL(
"Could not open file " << sInputFilePath.c_str());
56 return StatusCode::FAILURE;
63 std::vector<std::string> vInputFilePath;
74 return StatusCode::FAILURE;
76 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.
◆ isAffectedBySystematic()
◆ isSupportedRunNumber()
virtual bool TauAnalysisTools::CommonEfficiencyTool::isSupportedRunNumber |
( |
int |
iRunNumber | ) |
const |
|
inlinevirtualinherited |
◆ 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()
◆ 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() [1/2]
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.
◆ print() [2/2]
virtual void asg::IAsgTool::print |
( |
| ) |
const |
|
pure virtualinherited |
Print the state of the tool.
Implemented in JetRecTool, JetFinder, JetModifiedMassDrop, JetFromPseudojet, JetReclusterer, JetReclusteringTool, JetTruthLabelingTool, JetPileupLabelingTool, HI::HIPileupTool, top::TopObjectSelection, asg::AsgTool, JetDumper, JetBottomUpSoftDrop, JetRecursiveSoftDrop, JetSoftDrop, JetConstituentsRetriever, JetSubStructureMomentToolsBase, JetSplitter, JetToolRunner, JetPruner, JetPseudojetRetriever, JetTrimmer, AsgHelloTool, and KtDeltaRTool.
◆ ReadInputs() [1/2]
void CommonEfficiencyTool::ReadInputs |
( |
const TFile & |
fFile | ) |
|
|
protectedinherited |
Definition at line 470 of file CommonEfficiencyTool.cxx.
479 TIter itNext(fFile.GetListOfKeys());
480 while ((kKey = (TKey*)itNext()))
484 std::string sKeyName = kKey->GetName();
485 if (sKeyName ==
"Xaxis")
487 TNamed* tObj = (TNamed*)kKey->ReadObj();
488 std::string sTitle = tObj->GetTitle();
490 if (sTitle ==
"P" || sTitle ==
"PFinalCalib")
495 if (sTitle ==
"TruthDecayMode")
500 if (sTitle ==
"truth pt")
505 if (sTitle ==
"|eta|")
513 else if (sKeyName ==
"Yaxis")
515 TNamed* tObj = (TNamed*)kKey->ReadObj();
516 std::string sTitle = tObj->GetTitle();
518 if (sTitle ==
"track-eta")
523 else if (sTitle ==
"|eta|")
528 else if (sTitle ==
"mu")
532 if (
evtStore()->retrieve(xEventInfo,
"EventInfo").isFailure()) {
547 else if (sTitle ==
"truth |eta|")
555 std::vector<std::string> vSplitName = {};
556 split(sKeyName,
'_',vSplitName);
557 if (vSplitName[0] ==
"sf")
564 if (sKeyName.find(
"_up_") != std::string::npos or sKeyName.find(
"_down_") != std::string::npos)
568 size_t iPos = sKeyName.find(
'_');
◆ ReadInputs() [2/2]
void CommonDiTauEfficiencyTool::ReadInputs |
( |
std::unique_ptr< TFile > & |
fFile | ) |
|
Definition at line 192 of file CommonDiTauEfficiencyTool.cxx.
202 TIter itNext(fFile->GetListOfKeys());
203 while ((kKey = (TKey*)itNext()))
207 std::string sKeyName = kKey->GetName();
209 std::vector<std::string> vSplitName = {};
210 split(sKeyName,
'_',vSplitName);
211 if (vSplitName[0] ==
"sf")
217 if (sKeyName.find(
"_up_") != std::string::npos or sKeyName.find(
"_down_") != std::string::npos)
221 size_t iPos = sKeyName.find(
'_');
◆ recommendedSystematics()
◆ 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_bNoMultiprong
bool TauAnalysisTools::CommonEfficiencyTool::m_bNoMultiprong |
|
protectedinherited |
◆ m_bSFIsAvailable
bool TauAnalysisTools::CommonEfficiencyTool::m_bSFIsAvailable |
|
protectedinherited |
◆ m_bSFIsAvailableChecked
bool TauAnalysisTools::CommonEfficiencyTool::m_bSFIsAvailableChecked |
|
protectedinherited |
◆ m_bSFIsAvailableCheckedDiTau
bool TauAnalysisTools::CommonDiTauEfficiencyTool::m_bSFIsAvailableCheckedDiTau |
◆ m_bSFIsAvailableDiTau
bool TauAnalysisTools::CommonDiTauEfficiencyTool::m_bSFIsAvailableDiTau |
◆ m_bSkipTruthMatchCheck
bool TauAnalysisTools::CommonEfficiencyTool::m_bSkipTruthMatchCheck |
|
protectedinherited |
◆ m_bSplitMCCampaign
bool TauAnalysisTools::CommonEfficiencyTool::m_bSplitMCCampaign |
|
protectedinherited |
◆ m_bSplitMu
bool TauAnalysisTools::CommonEfficiencyTool::m_bSplitMu |
|
protectedinherited |
◆ m_bUseHighPtUncert
bool TauAnalysisTools::CommonEfficiencyTool::m_bUseHighPtUncert |
|
protectedinherited |
◆ m_bUseTauSubstructure
bool TauAnalysisTools::CommonEfficiencyTool::m_bUseTauSubstructure |
|
protectedinherited |
◆ m_detStore
◆ m_eCheckTruth
◆ m_evtStore
◆ m_fX
std::function<double(const xAOD::TauJet& xTau)> TauAnalysisTools::CommonEfficiencyTool::m_fX |
|
protectedinherited |
◆ m_fXDiTau
double(* TauAnalysisTools::CommonDiTauEfficiencyTool::m_fXDiTau) (const xAOD::DiTauJet &xDiTau) |
◆ m_fY
std::function<double(const xAOD::TauJet& xTau)> TauAnalysisTools::CommonEfficiencyTool::m_fY |
|
protectedinherited |
◆ m_fYDiTau
double(* TauAnalysisTools::CommonDiTauEfficiencyTool::m_fYDiTau) (const xAOD::DiTauJet &xDiTau) |
◆ m_fZDiTau
double(* TauAnalysisTools::CommonDiTauEfficiencyTool::m_fZDiTau) (const xAOD::DiTauJet &xDiTau) |
◆ m_iEleIDLevel
int TauAnalysisTools::CommonEfficiencyTool::m_iEleIDLevel |
|
protectedinherited |
◆ m_iJetIDLevel
int TauAnalysisTools::CommonEfficiencyTool::m_iJetIDLevel |
|
protectedinherited |
◆ m_mSF
std::unique_ptr< tSFMAP > TauAnalysisTools::CommonEfficiencyTool::m_mSF |
|
protectedinherited |
◆ m_mSystematics
std::map<std::string, int> TauAnalysisTools::CommonEfficiencyTool::m_mSystematics |
|
protectedinherited |
◆ m_mSystematicSets
std::unordered_map< CP::SystematicSet, std::string > TauAnalysisTools::CommonEfficiencyTool::m_mSystematicSets |
|
protectedinherited |
◆ m_mSystematicsHistNames
std::map<std::string, std::string> TauAnalysisTools::CommonEfficiencyTool::m_mSystematicsHistNames |
|
protectedinherited |
◆ m_sAffectingSystematics
◆ m_sInputFileName
std::string TauAnalysisTools::CommonEfficiencyTool::m_sInputFileName |
|
protectedinherited |
◆ m_sInputFilePath
std::string TauAnalysisTools::CommonEfficiencyTool::m_sInputFilePath |
|
protectedinherited |
◆ m_sMCCampaign
std::string TauAnalysisTools::CommonEfficiencyTool::m_sMCCampaign |
|
protectedinherited |
◆ m_sRecommendedSystematics
◆ m_sSFHistName
std::string TauAnalysisTools::CommonEfficiencyTool::m_sSFHistName |
|
protectedinherited |
◆ m_sSystematicSet
◆ m_sVarName
std::string TauAnalysisTools::CommonEfficiencyTool::m_sVarName |
|
protectedinherited |
◆ m_sWP
std::string TauAnalysisTools::CommonEfficiencyTool::m_sWP |
|
protectedinherited |
◆ m_varHandleArraysDeclared
◆ m_vhka
The documentation for this class was generated from the following files:
size_t size() const
returns: size of the set
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T > &t)
bool empty() const
returns: whether the set is empty
size_t nTracks(TauJetParameters::TauTrackFlag flag=TauJetParameters::TauTrackFlag::classifiedCharged) const
Class to wrap a set of SystematicVariations.
StoreGateSvc_t m_evtStore
Pointer to StoreGate (event store by default)
std::vector< SG::VarHandleKeyArray * > m_vhka
bool panTauDetail(TauJetParameters::PanTauDetails panTauDetail, int &value) const
Get and set values of pantau details variables via enum.
#define ATH_MSG_VERBOSE(x)
uint32_t runNumber() const
The current event's run number.
virtual void setOwner(IDataHandleHolder *o)=0
ServiceHandle< StoreGateSvc > & evtStore()
The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.
@ OutOfValidityRange
Input object is out of validity range.
@ Error
Some error happened during the object correction.
Helper class to provide type-safe access to aux data.
Class describing a tau jet.
Amg::Vector3D transform(Amg::Vector3D &v, Amg::Transform3D &tr)
Transform a point from a Trasformation3D.
const_iterator end() const
description: const iterator to the end of the set
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
float averageInteractionsPerCrossing() const
Average interactions per crossing for all BCIDs - for out-of-time pile-up.
void insert(const SystematicVariation &systematic)
description: insert a systematic into the set
Class describing the basic event information.
std::string PathResolverFindCalibFile(const std::string &logical_file_name)
iterator find(const SystematicVariation &sys) const
description: find an element in the set
@ Ok
The correction was done successfully.
#define ATH_MSG_WARNING(x)
Return value from object correction CP tools.
SG::VarHandleKey & vhKey()
Return a non-const reference to the HandleKey.
bool dEta(const xAOD::TauJet &tau, const xAOD::TauTrack &track, double &out)
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>
float actualInteractionsPerCrossing() const
Average interactions per crossing for the current BCID - for in-time pile-up.