ATLAS Offline Software
Public Member Functions | Public Attributes | Protected Types | Protected Member Functions | Static Protected Member Functions | Protected Attributes | Private Types | Private Member Functions | Private Attributes | List of all members
TauAnalysisTools::CommonDiTauEfficiencyTool Class Referenceabstract

#include <CommonDiTauEfficiencyTool.h>

Inheritance diagram for TauAnalysisTools::CommonDiTauEfficiencyTool:
Collaboration diagram for TauAnalysisTools::CommonDiTauEfficiencyTool:

Public Member Functions

 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
 

Public Attributes

double(* m_fXDiTau )(const xAOD::DiTauJet &xDiTau)
 scale factor bin x (e.g. More...
 
double(* m_fYDiTau )(const xAOD::DiTauJet &xDiTau)
 scale factor bin y (e.g. More...
 
double(* m_fZDiTau )(const xAOD::DiTauJet &xDiTau)
 scale factor bin z (e.g. More...
 
bool m_bSFIsAvailableDiTau
 true if scale factor name is already decorated More...
 
bool m_bSFIsAvailableCheckedDiTau
 true if cale factor name is already decorated has already been checked More...
 

Protected Types

typedef std::tuple< TObject *, CP::CorrectionCode(*)(const TObject *oObject, double &dEfficiencyScaleFactor, double dVars[]) > tTupleObjectFunc
 
typedef std::map< std::string, tTupleObjectFunctSFMAP
 

Protected Member Functions

std::string ConvertProngToString (const int iProngness) const
 
std::string ConvertMuToString (const int iMu) const
 
std::string ConvertDecayModeToString (const int iDecayMode) const
 
void ReadInputs (const TFile &fFile)
 
void addHistogramToSFMap (TKey *kKey, const std::string &sKeyName)
 
virtual CP::CorrectionCode getValue (const std::string &sHistName, const xAOD::TauJet &xTau, double &dEfficiencyScaleFactor) const
 
void renounceArray (SG::VarHandleKeyArray &handlesArray)
 remove all handles from I/O resolution More...
 
std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void > renounce (T &h)
 
void extraDeps_update_handler (Gaudi::Details::PropertyBase &ExtraDeps)
 Add StoreName to extra input/output deps as needed. More...
 

Static Protected Member Functions

static CP::CorrectionCode getValueTH1 (const TObject *oObject, double &dEfficiencyScaleFactor, double dVars[])
 
static CP::CorrectionCode getValueTH2 (const TObject *oObject, double &dEfficiencyScaleFactor, double dVars[])
 
static CP::CorrectionCode getValueTH3 (const TObject *oObject, double &dEfficiencyScaleFactor, double dVars[])
 
static CP::CorrectionCode getValueTF1 (const TObject *oObject, double &dEfficiencyScaleFactor, double dVars[])
 

Protected Attributes

std::unique_ptr< tSFMAPm_mSF
 
std::unordered_map< CP::SystematicSet, std::string > m_mSystematicSets
 
const CP::SystematicSetm_sSystematicSet
 
std::map< std::string, int > m_mSystematics
 
std::map< std::string, std::string > m_mSystematicsHistNames
 
std::function< double(const xAOD::TauJet &xTau)> m_fX
 
std::function< double(const xAOD::TauJet &xTau)> m_fY
 
CP::SystematicSet m_sAffectingSystematics
 
CP::SystematicSet m_sRecommendedSystematics
 
std::string m_sInputFilePath
 
std::string m_sInputFileName
 
std::string m_sWP
 
std::string m_sVarName
 
std::string m_sSFHistName
 
bool m_bSkipTruthMatchCheck
 
bool m_bNoMultiprong
 
bool m_bUseTauSubstructure
 
int m_iJetIDLevel
 
int m_iEleIDLevel
 
TruthMatchedParticleType m_eCheckTruth
 
bool m_bSFIsAvailable
 
bool m_bSFIsAvailableChecked
 
bool m_bSplitMu
 

Private Types

typedef ServiceHandle< StoreGateSvcStoreGateSvc_t
 

Private Member Functions

Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
 specialization for handling Gaudi::Property<SG::VarHandleKey> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyArrayType &)
 specialization for handling Gaudi::Property<SG::VarHandleKeyArray> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleType &)
 specialization for handling Gaudi::Property<SG::VarHandleBase> More...
 
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T > &t, const SG::NotHandleType &)
 specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray> More...
 

Private Attributes

StoreGateSvc_t m_evtStore
 Pointer to StoreGate (event store by default) More...
 
StoreGateSvc_t m_detStore
 Pointer to StoreGate (detector store by default) More...
 
std::vector< SG::VarHandleKeyArray * > m_vhka
 
bool m_varHandleArraysDeclared
 

Detailed Description

Definition at line 39 of file CommonDiTauEfficiencyTool.h.

Member Typedef Documentation

◆ StoreGateSvc_t

typedef ServiceHandle<StoreGateSvc> AthCommonDataStore< AthCommonMsg< AlgTool > >::StoreGateSvc_t
privateinherited

Definition at line 388 of file AthCommonDataStore.h.

◆ tSFMAP

typedef std::map<std::string, tTupleObjectFunc > TauAnalysisTools::CommonEfficiencyTool::tSFMAP
protectedinherited

Definition at line 88 of file CommonEfficiencyTool.h.

◆ tTupleObjectFunc

typedef std::tuple<TObject*, CP::CorrectionCode (*)(const TObject* oObject, double& dEfficiencyScaleFactor, double dVars[] ) > TauAnalysisTools::CommonEfficiencyTool::tTupleObjectFunc
protectedinherited

Definition at line 87 of file CommonEfficiencyTool.h.

Constructor & Destructor Documentation

◆ CommonDiTauEfficiencyTool()

CommonDiTauEfficiencyTool::CommonDiTauEfficiencyTool ( const std::string &  sName)

Create a proper constructor for Athena.

Definition at line 21 of file CommonDiTauEfficiencyTool.cxx.

22  : CommonEfficiencyTool( sName )
26  , m_bSFIsAvailableDiTau(false)
28 {
29 }

◆ ~CommonDiTauEfficiencyTool()

CommonDiTauEfficiencyTool::~CommonDiTauEfficiencyTool ( )

Definition at line 31 of file CommonDiTauEfficiencyTool.cxx.

32 {
33 }

Member Function Documentation

◆ addHistogramToSFMap()

void CommonEfficiencyTool::addHistogramToSFMap ( TKey *  kKey,
const std::string &  sKeyName 
)
protectedinherited

Definition at line 555 of file CommonEfficiencyTool.cxx.

556 {
557  // handling for the 3 different input types TH1F/TH1D/TF1, function pointer
558  // handle the access methods for the final scale factor retrieval
559  TClass *cClass = gROOT->GetClass(kKey->GetClassName());
560  if (cClass->InheritsFrom("TH2"))
561  {
562  TH1* oObject = (TH1*)kKey->ReadObj();
563  oObject->SetDirectory(0);
564  (*m_mSF)[sKeyName] = tTupleObjectFunc(oObject,&getValueTH2);
565  ATH_MSG_DEBUG("added histogram with name "<<sKeyName);
566  }
567  else if (cClass->InheritsFrom("TH3"))
568  {
569  TH1* oObject = (TH1*)kKey->ReadObj();
570  oObject->SetDirectory(0);
571  (*m_mSF)[sKeyName] = tTupleObjectFunc(oObject,&getValueTH3);
572  ATH_MSG_DEBUG("added histogram with name "<<sKeyName);
573  }else if (cClass->InheritsFrom("TH1"))
574  {
575  TH1* oObject = (TH1*)kKey->ReadObj();
576  oObject->SetDirectory(0);
577  (*m_mSF)[sKeyName] = tTupleObjectFunc(oObject,&getValueTH1);
578  ATH_MSG_DEBUG("added histogram with name "<<sKeyName);
579  }
580  else if (cClass->InheritsFrom("TF1"))
581  {
582  TObject* oObject = kKey->ReadObj();
583  (*m_mSF)[sKeyName] = tTupleObjectFunc(oObject,&getValueTF1);
584  ATH_MSG_DEBUG("added function with name "<<sKeyName);
585  }
586  else
587  {
588  ATH_MSG_DEBUG("ignored object with name "<<sKeyName);
589  }
590 }

◆ affectingSystematics()

CP::SystematicSet CommonEfficiencyTool::affectingSystematics ( ) const
virtualinherited

returns: the list of all systematics this tool can be affected by

Implements CP::IReentrantSystematicsTool.

Definition at line 317 of file CommonEfficiencyTool.cxx.

318 {
320 }

◆ applyEfficiencyScaleFactor() [1/3]

CP::CorrectionCode CommonDiTauEfficiencyTool::applyEfficiencyScaleFactor ( const xAOD::DiTauJet xDiTau,
unsigned int  iRunNumber = 0,
unsigned int  iMu = 0 
)
overridevirtual

Get the Efficiency Scale Factor of ditau jet.

Parameters
xDiTau
iRunNumber
iMu
Returns
CP::CorrectionCode

Implements TauAnalysisTools::IDiTauEfficiencyCorrectionsTool.

Definition at line 163 of file CommonDiTauEfficiencyTool.cxx.

165 {
166  double dSf = 0.;
167 
170  {
171  m_bSFIsAvailableDiTau = decor.isAvailable(xDiTau);
174  {
175  ATH_MSG_DEBUG(m_sVarName << " decoration is available on first ditau processed, switched of applyEfficiencyScaleFactor for further ditaus.");
176  ATH_MSG_DEBUG("If an application of efficiency scale factors needs to be redone, please pass a shallow copy of the original ditau.");
177  }
178  }
180  return CP::CorrectionCode::Ok;
181 
182  // retrieve scale factor
183  CP::CorrectionCode tmpCorrectionCode = getEfficiencyScaleFactor(xDiTau, dSf, iRunNumber, iMu);
184  // adding scale factor to tau as decoration
185  decor(xDiTau) = dSf;
186 
187  return tmpCorrectionCode;
188 }

◆ applyEfficiencyScaleFactor() [2/3]

CP::CorrectionCode CommonEfficiencyTool::applyEfficiencyScaleFactor

Definition at line 57 of file CommonEfficiencyTool.cxx.

277 {
278  double dSf = 0.;
279 
282  {
283  m_bSFIsAvailable = decor.isAvailable(xTau);
285  if (m_bSFIsAvailable)
286  {
287  ATH_MSG_DEBUG(m_sVarName << " decoration is available on first tau processed, switched off applyEfficiencyScaleFactor for further taus.");
288  ATH_MSG_DEBUG("If an application of efficiency scale factors needs to be redone, please pass a shallow copy of the original tau.");
289  }
290  }
291  if (m_bSFIsAvailable)
292  return CP::CorrectionCode::Ok;
293 
294  // retrieve scale factor
295  CP::CorrectionCode tmpCorrectionCode = getEfficiencyScaleFactor(xTau, dSf, iRunNumber, iMu);
296  // adding scale factor to tau as decoration
297  decor(xTau) = dSf;
298 
299  return tmpCorrectionCode;
300 }

◆ applyEfficiencyScaleFactor() [3/3]

CP::CorrectionCode CommonEfficiencyTool::applyEfficiencyScaleFactor ( const xAOD::TauJet xTau,
unsigned int  iRunNumber = 0,
unsigned int  iMu = 0 
)
virtualinherited

Decorate the tau with its efficiency.

Implements TauAnalysisTools::ITauEfficiencyCorrectionsTool.

Definition at line 275 of file CommonEfficiencyTool.cxx.

277 {
278  double dSf = 0.;
279 
282  {
283  m_bSFIsAvailable = decor.isAvailable(xTau);
285  if (m_bSFIsAvailable)
286  {
287  ATH_MSG_DEBUG(m_sVarName << " decoration is available on first tau processed, switched off applyEfficiencyScaleFactor for further taus.");
288  ATH_MSG_DEBUG("If an application of efficiency scale factors needs to be redone, please pass a shallow copy of the original tau.");
289  }
290  }
291  if (m_bSFIsAvailable)
292  return CP::CorrectionCode::Ok;
293 
294  // retrieve scale factor
295  CP::CorrectionCode tmpCorrectionCode = getEfficiencyScaleFactor(xTau, dSf, iRunNumber, iMu);
296  // adding scale factor to tau as decoration
297  decor(xTau) = dSf;
298 
299  return tmpCorrectionCode;
300 }

◆ applySystematicVariation()

StatusCode CommonEfficiencyTool::applySystematicVariation ( const CP::SystematicSet sSystematicSet)
virtualinherited

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 342 of file CommonEfficiencyTool.cxx.

343 {
344 
345  // first check if we already know this systematic configuration
346  auto itSystematicSet = m_mSystematicSets.find(sSystematicSet);
347  if (itSystematicSet != m_mSystematicSets.end())
348  {
349  m_sSystematicSet = &itSystematicSet->first;
350  return StatusCode::SUCCESS;
351  }
352 
353  // sanity checks if systematic set is supported
354  double dDirection = 0.;
355  CP::SystematicSet sSystematicSetAvailable;
356  for (auto sSyst : sSystematicSet)
357  {
358  // check if systematic is available
359  auto it = m_mSystematicsHistNames.find(sSyst.basename());
360  if (it == m_mSystematicsHistNames.end())
361  {
362  ATH_MSG_VERBOSE("unsupported systematic variation: "<< sSyst.basename()<<"; skipping this one");
363  continue;
364  }
365 
366 
367  if (sSyst.parameter() * dDirection < 0)
368  {
369  ATH_MSG_ERROR("unsupported set of systematic variations, you should either use only \"UP\" or only \"DOWN\" systematics in one set!");
370  ATH_MSG_ERROR("systematic set will not be applied");
371  return StatusCode::FAILURE;
372  }
373  dDirection = sSyst.parameter();
374 
375  if ((m_sRecommendedSystematics.find(sSyst.basename()) != m_sRecommendedSystematics.end()) and sSystematicSet.size() > 1)
376  {
377  ATH_MSG_ERROR("unsupported set of systematic variations, you should not combine \"TAUS_{TRUE|FAKE}_EFF_*_TOTAL\" with other systematic variations!");
378  ATH_MSG_ERROR("systematic set will not be applied");
379  return StatusCode::FAILURE;
380  }
381 
382  // finally add the systematic to the set of systematics to process
383  sSystematicSetAvailable.insert(sSyst);
384  }
385 
386  // store this calibration for future use, and make it current
387  m_sSystematicSet = &m_mSystematicSets.insert(std::pair<CP::SystematicSet,std::string>(sSystematicSetAvailable, sSystematicSet.name())).first->first;
388 
389  return StatusCode::SUCCESS;
390 }

◆ ConvertDecayModeToString()

std::string CommonEfficiencyTool::ConvertDecayModeToString ( const int  iDecayMode) const
protectedinherited

Definition at line 415 of file CommonEfficiencyTool.cxx.

416 {
417  switch(iDecayMode)
418  {
420  return "_r1p0n";
422  return "_r1p1n";
424  return "_r1pXn";
426  return "_r3p0n";
428  return "_r3pXn";
429  default:
430  return "";
431  }
432 }

◆ ConvertMuToString()

std::string CommonEfficiencyTool::ConvertMuToString ( const int  iMu) const
protectedinherited

Definition at line 403 of file CommonEfficiencyTool.cxx.

404 {
405  if (iMu > 35 )
406  return "_highMu";
407 
408  return "_lowMu";
409 }

◆ ConvertProngToString()

std::string CommonEfficiencyTool::ConvertProngToString ( const int  iProngness) const
protectedinherited

Definition at line 393 of file CommonEfficiencyTool.cxx.

394 {
395  return fProngness == 1 ? "_1p" : "_3p";
396 }

◆ declareGaudiProperty() [1/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty ( Gaudi::Property< T > &  hndl,
const SG::VarHandleKeyArrayType  
)
inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleKeyArray>

Definition at line 170 of file AthCommonDataStore.h.

172  {
173  return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
174  hndl.value(),
175  hndl.documentation());
176 
177  }

◆ declareGaudiProperty() [2/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty ( Gaudi::Property< T > &  hndl,
const SG::VarHandleKeyType  
)
inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleKey>

Definition at line 156 of file AthCommonDataStore.h.

158  {
159  return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
160  hndl.value(),
161  hndl.documentation());
162 
163  }

◆ declareGaudiProperty() [3/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty ( Gaudi::Property< T > &  hndl,
const SG::VarHandleType  
)
inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleBase>

Definition at line 184 of file AthCommonDataStore.h.

186  {
187  return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
188  hndl.value(),
189  hndl.documentation());
190  }

◆ declareGaudiProperty() [4/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty ( Gaudi::Property< T > &  t,
const SG::NotHandleType  
)
inlineprivateinherited

specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray>

Definition at line 199 of file AthCommonDataStore.h.

200  {
201  return PBASE::declareProperty(t);
202  }

◆ declareProperty() [1/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( const std::string &  name,
SG::VarHandleBase hndl,
const std::string &  doc,
const SG::VarHandleType  
)
inlineinherited

Declare a new Gaudi property.

Parameters
nameName of the property.
hndlObject holding the property value.
docDocumentation 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.

249  {
250  this->declare(hndl.vhKey());
251  hndl.vhKey().setOwner(this);
252 
253  return PBASE::declareProperty(name,hndl,doc);
254  }

◆ declareProperty() [2/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( const std::string &  name,
SG::VarHandleKey hndl,
const std::string &  doc,
const SG::VarHandleKeyType  
)
inlineinherited

Declare a new Gaudi property.

Parameters
nameName of the property.
hndlObject holding the property value.
docDocumentation 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.

225  {
226  this->declare(hndl);
227  hndl.setOwner(this);
228 
229  return PBASE::declareProperty(name,hndl,doc);
230  }

◆ declareProperty() [3/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( const std::string &  name,
SG::VarHandleKeyArray hndArr,
const std::string &  doc,
const SG::VarHandleKeyArrayType  
)
inlineinherited

Definition at line 259 of file AthCommonDataStore.h.

263  {
264 
265  // std::ostringstream ost;
266  // ost << Algorithm::name() << " VHKA declareProp: " << name
267  // << " size: " << hndArr.keys().size()
268  // << " mode: " << hndArr.mode()
269  // << " vhka size: " << m_vhka.size()
270  // << "\n";
271  // debug() << ost.str() << endmsg;
272 
273  hndArr.setOwner(this);
274  m_vhka.push_back(&hndArr);
275 
276  Gaudi::Details::PropertyBase* p = PBASE::declareProperty(name, hndArr, doc);
277  if (p != 0) {
278  p->declareUpdateHandler(&AthCommonDataStore<PBASE>::updateVHKA, this);
279  } else {
280  ATH_MSG_ERROR("unable to call declareProperty on VarHandleKeyArray "
281  << name);
282  }
283 
284  return p;
285 
286  }

◆ declareProperty() [4/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( const std::string &  name,
T &  property,
const std::string &  doc,
const SG::NotHandleType  
)
inlineinherited

Declare a new Gaudi property.

Parameters
nameName of the property.
propertyObject holding the property value.
docDocumentation 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.

337  {
338  return PBASE::declareProperty(name, property, doc);
339  }

◆ declareProperty() [5/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( const std::string &  name,
T &  property,
const std::string &  doc = "none" 
)
inlineinherited

Declare a new Gaudi property.

Parameters
nameName of the property.
propertyObject holding the property value.
docDocumentation 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.

355  {
356  typedef typename SG::HandleClassifier<T>::type htype;
357  return declareProperty (name, property, doc, htype());
358  }

◆ declareProperty() [6/6]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( Gaudi::Property< T > &  t)
inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

145  {
146  typedef typename SG::HandleClassifier<T>::type htype;
148  }

◆ detStore()

const ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< AlgTool > >::detStore ( ) const
inlineinherited

The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 95 of file AthCommonDataStore.h.

95 { return m_detStore; }

◆ evtStore() [1/2]

ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< AlgTool > >::evtStore ( )
inlineinherited

The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 85 of file AthCommonDataStore.h.

85 { return m_evtStore; }

◆ evtStore() [2/2]

const ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< AlgTool > >::evtStore ( ) const
inlineinherited

The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 90 of file AthCommonDataStore.h.

90 { return m_evtStore; }

◆ extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::extraDeps_update_handler ( Gaudi::Details::PropertyBase &  ExtraDeps)
protectedinherited

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.

245 {
246  // creation of basic string for all NPs, e.g. "TAUS_TRUEHADTAU_EFF_RECO_"
247  std::vector<std::string> vSplitInputFilePath = {};
248  split(m_sInputFileName,'_',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+"_";
254  // set truth type to check for in truth matching
255  if (sTruthType=="TRUEHADTAU") m_eCheckTruth = TauAnalysisTools::TruthHadronicTau;
256  if (sTruthType=="TRUEELECTRON") m_eCheckTruth = TauAnalysisTools::TruthElectron;
257  if (sTruthType=="TRUEMUON") m_eCheckTruth = TauAnalysisTools::TruthMuon;
258  if (sTruthType=="TRUEJET") m_eCheckTruth = TauAnalysisTools::TruthJet;
259  if (sTruthType=="TRUEHADDITAU") m_eCheckTruth = TauAnalysisTools::TruthHadronicDiTau;
260  if (sEfficiencyType=="ELEOLR") m_bNoMultiprong = true;
261  for (auto mSF : *m_mSF)
262  {
263  // parse for nuisance parameter in histogram name
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);
268  // skip nominal scale factors
269  if (sNP == "sf") continue;
270  // test if NP starts with a capital letter indicating that this should be recommended
271  bool bIsRecommended = false;
272  if (isupper(sNP.at(0)))
273  bIsRecommended = true;
274  // make sNP uppercase and build final NP entry name
275  std::transform(sNPUppercase.begin(), sNPUppercase.end(), sNPUppercase.begin(), toupper);
276  std::string sSystematicString = sSystematicBaseString+sNPUppercase;
277  // add all found systematics to the AffectingSystematics
279  m_sAffectingSystematics.insert(CP::SystematicVariation (sSystematicString, -1));
280  // only add found uppercase systematics to the RecommendedSystematics
281  if (bIsRecommended)
282  {
285  }
286  ATH_MSG_DEBUG("connected base name " << sNP << " with systematic " <<sSystematicString);
287  m_mSystematicsHistNames.insert({sSystematicString,sNP});
288  }
289 }

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

88 {
89  // check which true state is requestet
90  if (getTruthParticleType(xDiTau) != m_eCheckTruth)
91  {
92  dEfficiencyScaleFactor = 1.;
94  }
95 
96  CP::CorrectionCode tmpCorrectionCode = getValue(m_sSFHistName,
97  xDiTau,
98  dEfficiencyScaleFactor);
99  // return correction code if histogram is not available
100  if (tmpCorrectionCode != CP::CorrectionCode::Ok)
101  return tmpCorrectionCode;
102 
103  // skip further process if systematic set is empty
104  if (m_sSystematicSet->size() == 0)
105  return CP::CorrectionCode::Ok;
106 
107  // get uncertainties summed in quadrature
108  double dTotalSystematic2 = 0.;
109  double dDirection = 0.;
110  for (auto syst : *m_sSystematicSet)
111  {
112 
113  // check if systematic is available
114  auto it = m_mSystematicsHistNames.find(syst.basename());
115 
116  // get uncertainty value
117  double dUncertaintySyst = 0.;
118 
119  // needed for up/down decision
120  dDirection = syst.parameter();
121 
122  // build up histogram name
123  std::string sHistName = it->second;
124  if (dDirection>0) sHistName+="_up";
125  else sHistName+="_down";
126  if (!m_sWP.empty()) sHistName+="_"+m_sWP;
127 
128  // get the uncertainty from the histogram
129  tmpCorrectionCode = getValue(sHistName,
130  xDiTau,
131  dUncertaintySyst);
132 
133  // return correction code if histogram is not available
134  if (tmpCorrectionCode != CP::CorrectionCode::Ok)
135  return tmpCorrectionCode;
136 
137  // scale uncertainty with direction, i.e. +/- n*sigma
138  dUncertaintySyst *= dDirection;
139 
140  // square uncertainty and add to total uncertainty
141  dTotalSystematic2 += dUncertaintySyst * dUncertaintySyst;
142  }
143 
144  // now use dDirection to use up/down uncertainty
145  dDirection = (dDirection > 0.) ? 1. : -1.;
146 
147  // finally apply uncertainty (eff * ( 1 +/- \sum )
148  dEfficiencyScaleFactor *= 1. + dDirection * std::sqrt(dTotalSystematic2);
149 
150  return CP::CorrectionCode::Ok;
151 }

◆ getEfficiencyScaleFactor() [2/3]

CP::CorrectionCode CommonEfficiencyTool::getEfficiencyScaleFactor ( const xAOD::TauJet xTau,
double &  eff,
unsigned int  iRunNumber = 0,
unsigned int  iMu = 0 
)
virtualinherited

Declare the interface that the class provides.

Get the tau efficiency scale factor

Implements TauAnalysisTools::ITauEfficiencyCorrectionsTool.

Reimplemented in TauAnalysisTools::TauEfficiencyTriggerTool.

Definition at line 155 of file CommonEfficiencyTool.cxx.

157 {
158  // check which true state is requested
160  {
161  dEfficiencyScaleFactor = 1.;
162  return CP::CorrectionCode::Ok;
163  }
164 
165  // check if 1 prong
166  if (m_bNoMultiprong && xTau.nTracks() != 1)
167  {
168  dEfficiencyScaleFactor = 1.;
169  return CP::CorrectionCode::Ok;
170  }
171 
172  // get decay mode or prong extension for histogram name
173  std::string sMode;
175  {
176  int iDecayMode = -1;
178  sMode = ConvertDecayModeToString(iDecayMode);
179  if (sMode.empty())
180  {
181  ATH_MSG_WARNING("Found tau with unknown decay mode. Skip efficiency correction.");
183  }
184  }
185  else
186  {
187  // skip taus which are not 1 or 3 prong
188  if( xTau.nTracks() != 1 && xTau.nTracks() != 3) {
189  dEfficiencyScaleFactor = 1.;
190  return CP::CorrectionCode::Ok;
191  }
192 
193  sMode = ConvertProngToString(xTau.nTracks());
194  }
195 
196  std::string sMu = "";
197  std::string sMCCampaign = "";
198 
199  if (m_bSplitMu) sMu = ConvertMuToString(iMu);
200  std::string sHistName = m_sSFHistName + sMode + sMu;
201 
202  // get standard scale factor
203  CP::CorrectionCode tmpCorrectionCode = getValue(sHistName,
204  xTau,
205  dEfficiencyScaleFactor);
206  // return correction code if histogram is not available
207  if (tmpCorrectionCode != CP::CorrectionCode::Ok)
208  return tmpCorrectionCode;
209 
210  // skip further process if systematic set is empty
211  if (m_sSystematicSet->empty())
212  return CP::CorrectionCode::Ok;
213 
214  // get uncertainties summed in quadrature
215  double dTotalSystematic2 = 0.;
216  double dDirection = 0.;
217  for (auto syst : *m_sSystematicSet)
218  {
219  // check if systematic is available
220  auto it = m_mSystematicsHistNames.find(syst.basename());
221 
222  // get uncertainty value
223  double dUncertaintySyst = 0.;
224 
225  // needed for up/down decision
226  dDirection = syst.parameter();
227 
228  // build up histogram name
229  sHistName = it->second;
230  if (dDirection>0.) sHistName+="_up";
231  else sHistName+="_down";
232  if (!m_sWP.empty()) sHistName+="_"+m_sWP;
233  sHistName += sMode + sMu + sMCCampaign;
234 
235  // filter unwanted combinations
236  if( (sHistName.find("3P") != std::string::npos && sHistName.find("1p") != std::string::npos) ||
237  (sHistName.find("1P") != std::string::npos && sHistName.find("3p") != std::string::npos))
238  continue;
239 
240  // get the uncertainty from the histogram
241  tmpCorrectionCode = getValue(sHistName,
242  xTau,
243  dUncertaintySyst);
244 
245  // return correction code if histogram is not available
246  if (tmpCorrectionCode != CP::CorrectionCode::Ok)
247  return tmpCorrectionCode;
248 
249  // scale uncertainty with direction, i.e. +/- n*sigma
250  dUncertaintySyst *= dDirection;
251 
252  // square uncertainty and add to total uncertainty
253  dTotalSystematic2 += dUncertaintySyst * dUncertaintySyst;
254  }
255 
256  // now use dDirection to use up/down uncertainty
257  dDirection = (dDirection > 0.) ? 1. : -1.;
258 
259  // finally apply uncertainty (eff * ( 1 +/- \sum )
260  dEfficiencyScaleFactor *= 1. + dDirection * std::sqrt(dTotalSystematic2);
261 
262  return CP::CorrectionCode::Ok;
263 }

◆ getEfficiencyScaleFactor() [3/3]

CP::CorrectionCode CommonEfficiencyTool::getEfficiencyScaleFactor

Definition at line 54 of file CommonEfficiencyTool.cxx.

157 {
158  // check which true state is requested
160  {
161  dEfficiencyScaleFactor = 1.;
162  return CP::CorrectionCode::Ok;
163  }
164 
165  // check if 1 prong
166  if (m_bNoMultiprong && xTau.nTracks() != 1)
167  {
168  dEfficiencyScaleFactor = 1.;
169  return CP::CorrectionCode::Ok;
170  }
171 
172  // get decay mode or prong extension for histogram name
173  std::string sMode;
175  {
176  int iDecayMode = -1;
177  xTau.panTauDetail(xAOD::TauJetParameters::PanTau_DecayMode, iDecayMode);
178  sMode = ConvertDecayModeToString(iDecayMode);
179  if (sMode.empty())
180  {
181  ATH_MSG_WARNING("Found tau with unknown decay mode. Skip efficiency correction.");
183  }
184  }
185  else
186  {
187  // skip taus which are not 1 or 3 prong
188  if( xTau.nTracks() != 1 && xTau.nTracks() != 3) {
189  dEfficiencyScaleFactor = 1.;
190  return CP::CorrectionCode::Ok;
191  }
192 
193  sMode = ConvertProngToString(xTau.nTracks());
194  }
195 
196  std::string sMu = "";
197  std::string sMCCampaign = "";
198 
199  if (m_bSplitMu) sMu = ConvertMuToString(iMu);
200  std::string sHistName = m_sSFHistName + sMode + sMu;
201 
202  // get standard scale factor
203  CP::CorrectionCode tmpCorrectionCode = getValue(sHistName,
204  xTau,
205  dEfficiencyScaleFactor);
206  // return correction code if histogram is not available
207  if (tmpCorrectionCode != CP::CorrectionCode::Ok)
208  return tmpCorrectionCode;
209 
210  // skip further process if systematic set is empty
211  if (m_sSystematicSet->empty())
212  return CP::CorrectionCode::Ok;
213 
214  // get uncertainties summed in quadrature
215  double dTotalSystematic2 = 0.;
216  double dDirection = 0.;
217  for (auto syst : *m_sSystematicSet)
218  {
219  // check if systematic is available
220  auto it = m_mSystematicsHistNames.find(syst.basename());
221 
222  // get uncertainty value
223  double dUncertaintySyst = 0.;
224 
225  // needed for up/down decision
226  dDirection = syst.parameter();
227 
228  // build up histogram name
229  sHistName = it->second;
230  if (dDirection>0.) sHistName+="_up";
231  else sHistName+="_down";
232  if (!m_sWP.empty()) sHistName+="_"+m_sWP;
233  sHistName += sMode + sMu + sMCCampaign;
234 
235  // filter unwanted combinations
236  if( (sHistName.find("3P") != std::string::npos && sHistName.find("1p") != std::string::npos) ||
237  (sHistName.find("1P") != std::string::npos && sHistName.find("3p") != std::string::npos))
238  continue;
239 
240  // get the uncertainty from the histogram
241  tmpCorrectionCode = getValue(sHistName,
242  xTau,
243  dUncertaintySyst);
244 
245  // return correction code if histogram is not available
246  if (tmpCorrectionCode != CP::CorrectionCode::Ok)
247  return tmpCorrectionCode;
248 
249  // scale uncertainty with direction, i.e. +/- n*sigma
250  dUncertaintySyst *= dDirection;
251 
252  // square uncertainty and add to total uncertainty
253  dTotalSystematic2 += dUncertaintySyst * dUncertaintySyst;
254  }
255 
256  // now use dDirection to use up/down uncertainty
257  dDirection = (dDirection > 0.) ? 1. : -1.;
258 
259  // finally apply uncertainty (eff * ( 1 +/- \sum )
260  dEfficiencyScaleFactor *= 1. + dDirection * std::sqrt(dTotalSystematic2);
261 
262  return CP::CorrectionCode::Ok;
263 }

◆ getKey()

SG::sgkey_t asg::AsgTool::getKey ( const void *  ptr) const
inherited

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

119  {
120 
121 #ifdef XAOD_STANDALONE
122  // In case we use @c xAOD::TEvent, we have a direct function call
123  // for this.
124  return evtStore()->event()->getKey( ptr );
125 #else
126  const SG::DataProxy* proxy = evtStore()->proxy( ptr );
127  return ( proxy == nullptr ? 0 : proxy->sgkey() );
128 #endif // XAOD_STANDALONE
129  }

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

106  {
107 
108 #ifdef XAOD_STANDALONE
109  // In case we use @c xAOD::TEvent, we have a direct function call
110  // for this.
111  return evtStore()->event()->getName( ptr );
112 #else
113  const SG::DataProxy* proxy = evtStore()->proxy( ptr );
114  static const std::string dummy = "";
115  return ( proxy == nullptr ? dummy : proxy->name() );
116 #endif // XAOD_STANDALONE
117  }

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

CP::CorrectionCode CommonDiTauEfficiencyTool::getValue ( const std::string &  sHistName,
const xAOD::DiTauJet xDiTau,
double &  dEfficiencyScaleFactor 
) const

Get the scale factor from a particular recommendations histogram.

Parameters
sHistName
xDiTau
dEfficiencyScaleFactor
Returns
CP::CorrectionCode

Definition at line 296 of file CommonDiTauEfficiencyTool.cxx.

299 {
300  const tSFMAP& mSF = *m_mSF;
301  auto it = mSF.find (sHistName);
302  if (it == mSF.end())
303  {
304  ATH_MSG_ERROR("Object with name "<<sHistName<<" was not found in input file.");
305  ATH_MSG_DEBUG("Content of input file");
306  for (auto eEntry : mSF)
307  ATH_MSG_DEBUG(" Entry: "<<eEntry.first);
309  }
310 
311  // get a tuple (TObject*,functionPointer) from the scale factor map
312  tTupleObjectFunc tTuple = it->second;
313 
314  // get pt and eta (for x and y axis respectively)
315  double dX = m_fXDiTau(xDiTau);
316  double dY = m_fYDiTau(xDiTau);
317  double dZ = m_fZDiTau(xDiTau);
318 
319  double dVars[3] = {dX, dY, dZ};
320  // finally obtain efficiency scale factor from TH1F/TH1D/TF1, by calling the
321  // function pointer stored in the tuple from the scale factor map
322  return (std::get<1>(tTuple))(std::get<0>(tTuple), dEfficiencyScaleFactor, dVars);
323 }

◆ getValue() [2/3]

CP::CorrectionCode CommonEfficiencyTool::getValue

Definition at line 109 of file CommonEfficiencyTool.cxx.

672 {
673  const tSFMAP& mSF = *m_mSF;
674  auto it = mSF.find (sHistName);
675  if (it == mSF.end())
676  {
677  ATH_MSG_ERROR("Object with name "<<sHistName<<" was not found in input file.");
678  ATH_MSG_DEBUG("Content of input file");
679  for (auto eEntry : mSF)
680  ATH_MSG_DEBUG(" Entry: "<<eEntry.first);
682  }
683 
684  // get a tuple (TObject*,functionPointer) from the scale factor map
685  tTupleObjectFunc tTuple = it->second;
686 
687  // get pt and eta (for x and y axis respectively)
688  double dPt = m_fX(xTau);
689  double dEta = m_fY(xTau);
690 
691  double dVars[2] = {dPt, dEta};
692 
693  // finally obtain efficiency scale factor from TH1F/TH1D/TF1, by calling the
694  // function pointer stored in the tuple from the scale factor map
695  return (std::get<1>(tTuple))(std::get<0>(tTuple), dEfficiencyScaleFactor, dVars);
696 }

◆ getValue() [3/3]

CP::CorrectionCode CommonEfficiencyTool::getValue ( const std::string &  sHistName,
const xAOD::TauJet xTau,
double &  dEfficiencyScaleFactor 
) const
protectedvirtualinherited

Definition at line 669 of file CommonEfficiencyTool.cxx.

672 {
673  const tSFMAP& mSF = *m_mSF;
674  auto it = mSF.find (sHistName);
675  if (it == mSF.end())
676  {
677  ATH_MSG_ERROR("Object with name "<<sHistName<<" was not found in input file.");
678  ATH_MSG_DEBUG("Content of input file");
679  for (auto eEntry : mSF)
680  ATH_MSG_DEBUG(" Entry: "<<eEntry.first);
682  }
683 
684  // get a tuple (TObject*,functionPointer) from the scale factor map
685  tTupleObjectFunc tTuple = it->second;
686 
687  // get pt and eta (for x and y axis respectively)
688  double dPt = m_fX(xTau);
689  double dEta = m_fY(xTau);
690 
691  double dVars[2] = {dPt, dEta};
692 
693  // finally obtain efficiency scale factor from TH1F/TH1D/TF1, by calling the
694  // function pointer stored in the tuple from the scale factor map
695  return (std::get<1>(tTuple))(std::get<0>(tTuple), dEfficiencyScaleFactor, dVars);
696 }

◆ getValueTF1()

CP::CorrectionCode CommonEfficiencyTool::getValueTF1 ( const TObject *  oObject,
double &  dEfficiencyScaleFactor,
double  dVars[] 
)
staticprotectedinherited

Definition at line 801 of file CommonEfficiencyTool.cxx.

803 {
804  double dPt = dVars[0];
805  double dEta = dVars[1];
806 
807  const TF1* fFunc = static_cast<const TF1*>(oObject);
808 
809  if (!fFunc)
810  {
811  // ATH_MSG_ERROR("Problem with casting TObject of type "<<oObject->ClassName()<<" to TF1");
813  }
814 
815  // evaluate TFunction and set scale factor
816  dEfficiencyScaleFactor = fFunc->Eval(dPt, dEta);
817  return CP::CorrectionCode::Ok;
818 }

◆ getValueTH1()

CP::CorrectionCode CommonEfficiencyTool::getValueTH1 ( const TObject *  oObject,
double &  dEfficiencyScaleFactor,
double  dVars[] 
)
staticprotectedinherited

Definition at line 704 of file CommonEfficiencyTool.cxx.

706 {
707  double dPt = dVars[0];
708 
709  const TH1* hHist = dynamic_cast<const TH1*>(oObject);
710 
711  if (!hHist)
712  {
713  // ATH_MSG_ERROR("Problem with casting TObject of type "<<oObject->ClassName()<<" to TH2F");
715  }
716 
717  // protect values from underflow bins
718  dPt = std::max(dPt,hHist->GetXaxis()->GetXmin());
719  // protect values from overflow bins (times .999 to keep it inside last bin)
720  dPt = std::min(dPt,hHist->GetXaxis()->GetXmax() * .999);
721 
722  // get bin from TH2 depending on x and y values; finally set the scale factor
723  int iBin = hHist->FindFixBin(dPt);
724  dEfficiencyScaleFactor = hHist->GetBinContent(iBin);
725  return CP::CorrectionCode::Ok;
726 }

◆ getValueTH2()

CP::CorrectionCode CommonEfficiencyTool::getValueTH2 ( const TObject *  oObject,
double &  dEfficiencyScaleFactor,
double  dVars[] 
)
staticprotectedinherited

Definition at line 734 of file CommonEfficiencyTool.cxx.

736 {
737  double dPt = dVars[0];
738  double dEta = dVars[1];
739 
740  const TH2* hHist = dynamic_cast<const TH2*>(oObject);
741 
742  if (!hHist)
743  {
744  // ATH_MSG_ERROR("Problem with casting TObject of type "<<oObject->ClassName()<<" to TH2F");
746  }
747 
748  // protect values from underflow bins
749  dPt = std::max(dPt,hHist->GetXaxis()->GetXmin());
750  dEta = std::max(dEta,hHist->GetYaxis()->GetXmin());
751  // protect values from overflow bins (times .999 to keep it inside last bin)
752  dPt = std::min(dPt,hHist->GetXaxis()->GetXmax() * .999);
753  dEta = std::min(dEta,hHist->GetYaxis()->GetXmax() * .999);
754 
755  // get bin from TH2 depending on x and y values; finally set the scale factor
756  int iBin = hHist->FindFixBin(dPt,dEta);
757  dEfficiencyScaleFactor = hHist->GetBinContent(iBin);
758  return CP::CorrectionCode::Ok;
759 }

◆ getValueTH3()

CP::CorrectionCode CommonEfficiencyTool::getValueTH3 ( const TObject *  oObject,
double &  dEfficiencyScaleFactor,
double  dVars[] 
)
staticprotectedinherited

Definition at line 766 of file CommonEfficiencyTool.cxx.

768 {
769  double dX = dVars[0];
770  double dY = dVars[1];
771  double dZ = dVars[2];
772 
773  const TH3* hHist = dynamic_cast<const TH3*>(oObject);
774 
775  if (!hHist)
776  {
777  // ATH_MSG_ERROR("Problem with casting TObject of type "<<oObject->ClassName()<<" to TH2D");
779  }
780 
781  // protect values from underflow bins
782  dX = std::max(dX,hHist->GetXaxis()->GetXmin());
783  dY = std::max(dY,hHist->GetYaxis()->GetXmin());
784  dZ = std::max(dZ,hHist->GetZaxis()->GetXmin());
785  // protect values from overflow bins (times .999 to keep it inside last bin)
786  dX = std::min(dX,hHist->GetXaxis()->GetXmax() * .999);
787  dY = std::min(dY,hHist->GetYaxis()->GetXmax() * .999);
788  dZ = std::min(dZ,hHist->GetZaxis()->GetXmax() * .999);
789 
790  // get bin from TH2 depending on x and y values; finally set the scale factor
791  int iBin = hHist->FindFixBin(dX,dY,dZ);
792  dEfficiencyScaleFactor = hHist->GetBinContent(iBin);
793  return CP::CorrectionCode::Ok;
794 }

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

44 {
45  ATH_MSG_INFO( "Initializing CommonDiTauEfficiencyTool" );
46  // only read in histograms once
47  if (m_mSF==nullptr)
48  {
49  std::string sInputFilePath = PathResolverFindCalibFile(m_sInputFilePath);
50 
51  m_mSF = std::make_unique< tSFMAP >();
52  std::unique_ptr< TFile > fSF( TFile::Open( sInputFilePath.c_str(), "READ" ) );
53  if(!fSF)
54  {
55  ATH_MSG_FATAL("Could not open file " << sInputFilePath.c_str());
56  return StatusCode::FAILURE;
57  }
58  ReadInputs(fSF);
59  fSF->Close();
60  }
61 
62  // needed later on in generateSystematicSets(), maybe move it there
63  std::vector<std::string> vInputFilePath;
64  split(m_sInputFilePath,'/',vInputFilePath);
65  m_sInputFileName = vInputFilePath.back();
66 
68 
69  if (m_sWP.length()>0)
70  m_sSFHistName = "sf_"+m_sWP;
71 
72  // load empty systematic variation by default
73  if (applySystematicVariation(CP::SystematicSet()) != StatusCode::SUCCESS )
74  return StatusCode::FAILURE;
75 
76  return StatusCode::SUCCESS;
77 }

◆ inputHandles()

virtual std::vector<Gaudi::DataHandle*> AthCommonDataStore< AthCommonMsg< AlgTool > >::inputHandles ( ) const
overridevirtualinherited

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()

bool CommonEfficiencyTool::isAffectedBySystematic ( const CP::SystematicVariation systematic) const
virtualinherited

returns: whether this tool is affected by the given systematics

Implements CP::ISystematicsTool.

Definition at line 306 of file CommonEfficiencyTool.cxx.

307 {
309  return sys.find(systematic) != sys.end();
310 }

◆ isSupportedRunNumber()

virtual bool TauAnalysisTools::CommonEfficiencyTool::isSupportedRunNumber ( int  iRunNumber) const
inlinevirtualinherited

check if run number is supported in recommendations

Implements TauAnalysisTools::ITauEfficiencyCorrectionsTool.

Reimplemented in TauAnalysisTools::TauEfficiencyTriggerTool.

Definition at line 73 of file CommonEfficiencyTool.h.

74  {
75  return true;
76  };

◆ msg() [1/2]

MsgStream& AthCommonMsg< AlgTool >::msg ( ) const
inlineinherited

Definition at line 24 of file AthCommonMsg.h.

24  {
25  return this->msgStream();
26  }

◆ msg() [2/2]

MsgStream& AthCommonMsg< AlgTool >::msg ( const MSG::Level  lvl) const
inlineinherited

Definition at line 27 of file AthCommonMsg.h.

27  {
28  return this->msgStream(lvl);
29  }

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

101  {
102 
103  return MSG::name( msg().level() );
104  }

◆ msgLvl()

bool AthCommonMsg< AlgTool >::msgLvl ( const MSG::Level  lvl) const
inlineinherited

Definition at line 30 of file AthCommonMsg.h.

30  {
31  return this->msgLevel(lvl);
32  }

◆ outputHandles()

virtual std::vector<Gaudi::DataHandle*> AthCommonDataStore< AthCommonMsg< AlgTool > >::outputHandles ( ) const
overridevirtualinherited

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() [2/2]

virtual void asg::IAsgTool::print ( ) const
pure virtualinherited

◆ ReadInputs() [1/2]

void CommonEfficiencyTool::ReadInputs ( const TFile &  fFile)
protectedinherited

Definition at line 444 of file CommonEfficiencyTool.cxx.

445 {
446  m_mSF->clear();
447 
448  // initialize function pointer
449  m_fX = &finalTauPt;
450  m_fY = &finalTauEta;
451 
452  TKey *kKey;
453  TIter itNext(fFile.GetListOfKeys());
454  while ((kKey = (TKey*)itNext()))
455  {
456  // parse file content for objects of type TNamed, check their title for
457  // known strings and reset funtion pointer
458  std::string sKeyName = kKey->GetName();
459  if (sKeyName == "Xaxis")
460  {
461  TNamed* tObj = (TNamed*)kKey->ReadObj();
462  std::string sTitle = tObj->GetTitle();
463  delete tObj;
464  if (sTitle == "P" || sTitle == "PFinalCalib")
465  {
466  m_fX = &finalTauP;
467  ATH_MSG_DEBUG("using full momentum for x-axis");
468  }
469  if (sTitle == "TruthDecayMode")
470  {
471  m_fX = &truthDecayMode;
472  ATH_MSG_DEBUG("using truth decay mode for x-axis");
473  }
474  if (sTitle == "truth pt")
475  {
476  m_fX = &truthTauPt;
477  ATH_MSG_DEBUG("using truth pT for x-axis");
478  }
479  if (sTitle == "|eta|")
480  {
481  m_fX = &finalTauAbsEta;
482  ATH_MSG_DEBUG("using absolute tau eta for x-axis");
483  }
484 
485  continue;
486  }
487  else if (sKeyName == "Yaxis")
488  {
489  TNamed* tObj = (TNamed*)kKey->ReadObj();
490  std::string sTitle = tObj->GetTitle();
491  delete tObj;
492  if (sTitle == "track-eta")
493  {
495  ATH_MSG_DEBUG("using leading track eta for y-axis");
496  }
497  else if (sTitle == "|eta|")
498  {
499  m_fY = &finalTauAbsEta;
500  ATH_MSG_DEBUG("using absolute tau eta for y-axis");
501  }
502  else if (sTitle == "mu")
503  {
504  m_fY = [this](const xAOD::TauJet&) -> double {
505  const xAOD::EventInfo* xEventInfo = nullptr;
506  if (evtStore()->retrieve(xEventInfo,"EventInfo").isFailure()) {
507  return 0;
508  }
509  if (xEventInfo->runNumber()==284500)
510  {
511  return xEventInfo->averageInteractionsPerCrossing();
512  }
513  else if (xEventInfo->runNumber()==300000 || xEventInfo->runNumber()==310000)
514  {
515  return xEventInfo->actualInteractionsPerCrossing();
516  }
517  return 0;
518  };
519  ATH_MSG_DEBUG("using average mu for y-axis");
520  }
521  else if (sTitle == "truth |eta|")
522  {
523  m_fY = &truthTauAbsEta;
524  ATH_MSG_DEBUG("using absolute truth tau eta for y-axis");
525  }
526  continue;
527  }
528 
529  std::vector<std::string> vSplitName = {};
530  split(sKeyName,'_',vSplitName);
531  if (vSplitName[0] == "sf")
532  {
533  addHistogramToSFMap(kKey, sKeyName);
534  }
535  else
536  {
537  // std::string sDirection = vSplitName[1];
538  if (sKeyName.find("_up_") != std::string::npos or sKeyName.find("_down_") != std::string::npos)
539  addHistogramToSFMap(kKey, sKeyName);
540  else
541  {
542  size_t iPos = sKeyName.find('_');
543  addHistogramToSFMap(kKey, sKeyName.substr(0,iPos)+"_up"+sKeyName.substr(iPos));
544  addHistogramToSFMap(kKey, sKeyName.substr(0,iPos)+"_down"+sKeyName.substr(iPos));
545  }
546  }
547  }
548  ATH_MSG_INFO("data loaded from " << fFile.GetName());
549 }

◆ ReadInputs() [2/2]

void CommonDiTauEfficiencyTool::ReadInputs ( std::unique_ptr< TFile > &  fFile)

Definition at line 192 of file CommonDiTauEfficiencyTool.cxx.

193 {
194  m_mSF->clear();
195 
196  // initialize function pointer
200 
201  TKey *kKey;
202  TIter itNext(fFile->GetListOfKeys());
203  while ((kKey = (TKey*)itNext()))
204  {
205  // parse file content for objects of type TNamed, check their title for
206  // known strings and reset funtion pointer
207  std::string sKeyName = kKey->GetName();
208 
209  std::vector<std::string> vSplitName = {};
210  split(sKeyName,'_',vSplitName);
211  if (vSplitName[0] == "sf")
212  {
213  addHistogramToSFMap(kKey, sKeyName);
214  }
215  else
216  {
217  if (sKeyName.find("_up_") != std::string::npos or sKeyName.find("_down_") != std::string::npos)
218  addHistogramToSFMap(kKey, sKeyName);
219  else
220  {
221  size_t iPos = sKeyName.find('_');
222  addHistogramToSFMap(kKey, sKeyName.substr(0,iPos)+"_up"+sKeyName.substr(iPos));
223  addHistogramToSFMap(kKey, sKeyName.substr(0,iPos)+"_down"+sKeyName.substr(iPos));
224  }
225  }
226  }
227  ATH_MSG_INFO("data loaded from " << fFile->GetName());
228 }

◆ recommendedSystematics()

CP::SystematicSet CommonEfficiencyTool::recommendedSystematics ( ) const
virtualinherited

returns: the list of all systematics this tool recommends to use

Implements CP::IReentrantSystematicsTool.

Definition at line 326 of file CommonEfficiencyTool.cxx.

327 {
329 }

◆ renounce()

std::enable_if_t<std::is_void_v<std::result_of_t<decltype(&T::renounce)(T)> > && !std::is_base_of_v<SG::VarHandleKeyArray, T> && std::is_base_of_v<Gaudi::DataHandle, T>, void> AthCommonDataStore< AthCommonMsg< AlgTool > >::renounce ( T &  h)
inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

381  {
382  h.renounce();
383  PBASE::renounce (h);
384  }

◆ renounceArray()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::renounceArray ( SG::VarHandleKeyArray handlesArray)
inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

364  {
365  handlesArray.renounce();
366  }

◆ sysInitialize()

virtual StatusCode AthCommonDataStore< AthCommonMsg< AlgTool > >::sysInitialize ( )
overridevirtualinherited

Perform system initialization for an algorithm.

We override this to declare all the elements of handle key arrays at the end of initialization. See comments on updateVHKA.

Reimplemented in DerivationFramework::CfAthAlgTool, AthCheckedComponent< AthAlgTool >, AthCheckedComponent<::AthAlgTool >, and asg::AsgMetadataTool.

◆ sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< AlgTool > >::sysStart ( )
overridevirtualinherited

Handle START transition.

We override this in order to make sure that conditions handle keys can cache a pointer to the conditions container.

◆ updateVHKA()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::updateVHKA ( Gaudi::Details::PropertyBase &  )
inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

308  {
309  // debug() << "updateVHKA for property " << p.name() << " " << p.toString()
310  // << " size: " << m_vhka.size() << endmsg;
311  for (auto &a : m_vhka) {
312  std::vector<SG::VarHandleKey*> keys = a->keys();
313  for (auto k : keys) {
314  k->setOwner(this);
315  }
316  }
317  }

Member Data Documentation

◆ m_bNoMultiprong

bool TauAnalysisTools::CommonEfficiencyTool::m_bNoMultiprong
protectedinherited

Definition at line 143 of file CommonEfficiencyTool.h.

◆ m_bSFIsAvailable

bool TauAnalysisTools::CommonEfficiencyTool::m_bSFIsAvailable
protectedinherited

Definition at line 150 of file CommonEfficiencyTool.h.

◆ m_bSFIsAvailableChecked

bool TauAnalysisTools::CommonEfficiencyTool::m_bSFIsAvailableChecked
protectedinherited

Definition at line 151 of file CommonEfficiencyTool.h.

◆ m_bSFIsAvailableCheckedDiTau

bool TauAnalysisTools::CommonDiTauEfficiencyTool::m_bSFIsAvailableCheckedDiTau

true if cale factor name is already decorated has already been checked

Definition at line 108 of file CommonDiTauEfficiencyTool.h.

◆ m_bSFIsAvailableDiTau

bool TauAnalysisTools::CommonDiTauEfficiencyTool::m_bSFIsAvailableDiTau

true if scale factor name is already decorated

Definition at line 106 of file CommonDiTauEfficiencyTool.h.

◆ m_bSkipTruthMatchCheck

bool TauAnalysisTools::CommonEfficiencyTool::m_bSkipTruthMatchCheck
protectedinherited

Definition at line 142 of file CommonEfficiencyTool.h.

◆ m_bSplitMu

bool TauAnalysisTools::CommonEfficiencyTool::m_bSplitMu
protectedinherited

Definition at line 152 of file CommonEfficiencyTool.h.

◆ m_bUseTauSubstructure

bool TauAnalysisTools::CommonEfficiencyTool::m_bUseTauSubstructure
protectedinherited

Definition at line 144 of file CommonEfficiencyTool.h.

◆ m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore
privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

◆ m_eCheckTruth

TruthMatchedParticleType TauAnalysisTools::CommonEfficiencyTool::m_eCheckTruth
protectedinherited

Definition at line 148 of file CommonEfficiencyTool.h.

◆ m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_evtStore
privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

◆ m_fX

std::function<double(const xAOD::TauJet& xTau)> TauAnalysisTools::CommonEfficiencyTool::m_fX
protectedinherited

Definition at line 103 of file CommonEfficiencyTool.h.

◆ m_fXDiTau

double(* TauAnalysisTools::CommonDiTauEfficiencyTool::m_fXDiTau) (const xAOD::DiTauJet &xDiTau)

scale factor bin x (e.g.

lead match pT)

Definition at line 82 of file CommonDiTauEfficiencyTool.h.

◆ m_fY

std::function<double(const xAOD::TauJet& xTau)> TauAnalysisTools::CommonEfficiencyTool::m_fY
protectedinherited

Definition at line 104 of file CommonEfficiencyTool.h.

◆ m_fYDiTau

double(* TauAnalysisTools::CommonDiTauEfficiencyTool::m_fYDiTau) (const xAOD::DiTauJet &xDiTau)

scale factor bin y (e.g.

sublead match pT)

Definition at line 84 of file CommonDiTauEfficiencyTool.h.

◆ m_fZDiTau

double(* TauAnalysisTools::CommonDiTauEfficiencyTool::m_fZDiTau) (const xAOD::DiTauJet &xDiTau)

scale factor bin z (e.g.

dR match particles)

Definition at line 86 of file CommonDiTauEfficiencyTool.h.

◆ m_iEleIDLevel

int TauAnalysisTools::CommonEfficiencyTool::m_iEleIDLevel
protectedinherited

Definition at line 146 of file CommonEfficiencyTool.h.

◆ m_iJetIDLevel

int TauAnalysisTools::CommonEfficiencyTool::m_iJetIDLevel
protectedinherited

Definition at line 145 of file CommonEfficiencyTool.h.

◆ m_mSF

std::unique_ptr< tSFMAP > TauAnalysisTools::CommonEfficiencyTool::m_mSF
protectedinherited

Definition at line 95 of file CommonEfficiencyTool.h.

◆ m_mSystematics

std::map<std::string, int> TauAnalysisTools::CommonEfficiencyTool::m_mSystematics
protectedinherited

Definition at line 100 of file CommonEfficiencyTool.h.

◆ m_mSystematicSets

std::unordered_map< CP::SystematicSet, std::string > TauAnalysisTools::CommonEfficiencyTool::m_mSystematicSets
protectedinherited

Definition at line 98 of file CommonEfficiencyTool.h.

◆ m_mSystematicsHistNames

std::map<std::string, std::string> TauAnalysisTools::CommonEfficiencyTool::m_mSystematicsHistNames
protectedinherited

Definition at line 101 of file CommonEfficiencyTool.h.

◆ m_sAffectingSystematics

CP::SystematicSet TauAnalysisTools::CommonEfficiencyTool::m_sAffectingSystematics
protectedinherited

Definition at line 134 of file CommonEfficiencyTool.h.

◆ m_sInputFileName

std::string TauAnalysisTools::CommonEfficiencyTool::m_sInputFileName
protectedinherited

Definition at line 138 of file CommonEfficiencyTool.h.

◆ m_sInputFilePath

std::string TauAnalysisTools::CommonEfficiencyTool::m_sInputFilePath
protectedinherited

Definition at line 137 of file CommonEfficiencyTool.h.

◆ m_sRecommendedSystematics

CP::SystematicSet TauAnalysisTools::CommonEfficiencyTool::m_sRecommendedSystematics
protectedinherited

Definition at line 135 of file CommonEfficiencyTool.h.

◆ m_sSFHistName

std::string TauAnalysisTools::CommonEfficiencyTool::m_sSFHistName
protectedinherited

Definition at line 141 of file CommonEfficiencyTool.h.

◆ m_sSystematicSet

const CP::SystematicSet* TauAnalysisTools::CommonEfficiencyTool::m_sSystematicSet
protectedinherited

Definition at line 99 of file CommonEfficiencyTool.h.

◆ m_sVarName

std::string TauAnalysisTools::CommonEfficiencyTool::m_sVarName
protectedinherited

Definition at line 140 of file CommonEfficiencyTool.h.

◆ m_sWP

std::string TauAnalysisTools::CommonEfficiencyTool::m_sWP
protectedinherited

Definition at line 139 of file CommonEfficiencyTool.h.

◆ m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< AlgTool > >::m_varHandleArraysDeclared
privateinherited

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

std::vector<SG::VarHandleKeyArray*> AthCommonDataStore< AthCommonMsg< AlgTool > >::m_vhka
privateinherited

Definition at line 398 of file AthCommonDataStore.h.


The documentation for this class was generated from the following files:
xAOD::TauJetParameters::Mode_1p0n
@ Mode_1p0n
Definition: TauDefs.h:386
xAOD::TauJetParameters::PanTau_DecayMode
@ PanTau_DecayMode
Definition: TauDefs.h:360
TauAnalysisTools::CommonDiTauEfficiencyTool::getEfficiencyScaleFactor
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.
Definition: CommonDiTauEfficiencyTool.cxx:86
TauAnalysisTools::TruthElectron
@ TruthElectron
Definition: PhysicsAnalysis/TauID/TauAnalysisTools/TauAnalysisTools/Enums.h:100
TauAnalysisTools::CommonEfficiencyTool::m_eCheckTruth
TruthMatchedParticleType m_eCheckTruth
Definition: CommonEfficiencyTool.h:148
ATH_MSG_FATAL
#define ATH_MSG_FATAL(x)
Definition: AthMsgStreamMacros.h:34
xAOD::TauJetParameters::Mode_1p1n
@ Mode_1p1n
Definition: TauDefs.h:387
StateLessPT_NewConfig.proxy
proxy
Definition: StateLessPT_NewConfig.py:392
max
#define max(a, b)
Definition: cfImp.cxx:41
TauAnalysisTools::truthTauAbsEta
double truthTauAbsEta(const xAOD::TauJet &xTau)
return truth match tau eta (if hadronic truth tau match)
Definition: PhysicsAnalysis/TauID/TauAnalysisTools/Root/HelperFunctions.cxx:172
ATH_MSG_INFO
#define ATH_MSG_INFO(x)
Definition: AthMsgStreamMacros.h:31
TauAnalysisTools::CommonDiTauEfficiencyTool::m_fXDiTau
double(* m_fXDiTau)(const xAOD::DiTauJet &xDiTau)
scale factor bin x (e.g.
Definition: CommonDiTauEfficiencyTool.h:82
TauAnalysisTools::CommonEfficiencyTool::getEfficiencyScaleFactor
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.
Definition: CommonEfficiencyTool.cxx:155
TauAnalysisTools::TruthSubleadPt
double TruthSubleadPt(const xAOD::DiTauJet &xDiTau)
return the truth vis pT of the subleading pT matched particle.
Definition: CommonDiTauEfficiencyTool.cxx:334
TauAnalysisTools::CommonEfficiencyTool::m_bSkipTruthMatchCheck
bool m_bSkipTruthMatchCheck
Definition: CommonEfficiencyTool.h:142
TauAnalysisTools::CommonDiTauEfficiencyTool::generateSystematicSets
void generateSystematicSets()
generate a set of relevant systematic variations to be applied
Definition: CommonDiTauEfficiencyTool.cxx:244
TauAnalysisTools::TruthHadronicTau
@ TruthHadronicTau
Definition: PhysicsAnalysis/TauID/TauAnalysisTools/TauAnalysisTools/Enums.h:97
CP::SystematicSet::size
size_t size() const
returns: size of the set
Definition: SystematicSet.h:71
AthCommonDataStore::declareProperty
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T > &t)
Definition: AthCommonDataStore.h:145
TauAnalysisTools::truthTauPt
double truthTauPt(const xAOD::TauJet &xTau)
return truth match tau pt in GeV (if hadronic truth tau match)
Definition: PhysicsAnalysis/TauID/TauAnalysisTools/Root/HelperFunctions.cxx:158
TauAnalysisTools::CommonEfficiencyTool::m_sSystematicSet
const CP::SystematicSet * m_sSystematicSet
Definition: CommonEfficiencyTool.h:99
CP::SystematicSet::empty
bool empty() const
returns: whether the set is empty
Definition: SystematicSet.h:67
xAOD::TauJet_v3::nTracks
size_t nTracks(TauJetParameters::TauTrackFlag flag=TauJetParameters::TauTrackFlag::classifiedCharged) const
Definition: TauJet_v3.cxx:526
skel.it
it
Definition: skel.GENtoEVGEN.py:396
TauAnalysisTools::CommonEfficiencyTool::ConvertProngToString
std::string ConvertProngToString(const int iProngness) const
Definition: CommonEfficiencyTool.cxx:393
CP::SystematicSet
Class to wrap a set of SystematicVariations.
Definition: SystematicSet.h:31
TauAnalysisTools::CommonEfficiencyTool::tSFMAP
std::map< std::string, tTupleObjectFunc > tSFMAP
Definition: CommonEfficiencyTool.h:88
AthCommonDataStore< AthCommonMsg< AlgTool > >::m_evtStore
StoreGateSvc_t m_evtStore
Pointer to StoreGate (event store by default)
Definition: AthCommonDataStore.h:390
AthCommonDataStore< AthCommonMsg< AlgTool > >::m_vhka
std::vector< SG::VarHandleKeyArray * > m_vhka
Definition: AthCommonDataStore.h:398
TauAnalysisTools::CommonEfficiencyTool::applySystematicVariation
virtual StatusCode applySystematicVariation(const CP::SystematicSet &sSystematicSet)
configure this tool for the given list of systematic variations.
Definition: CommonEfficiencyTool.cxx:342
TauAnalysisTools::CommonDiTauEfficiencyTool::m_fYDiTau
double(* m_fYDiTau)(const xAOD::DiTauJet &xDiTau)
scale factor bin y (e.g.
Definition: CommonDiTauEfficiencyTool.h:84
TauAnalysisTools::CommonEfficiencyTool::m_sVarName
std::string m_sVarName
Definition: CommonEfficiencyTool.h:140
xAOD::TauJet_v3::panTauDetail
bool panTauDetail(TauJetParameters::PanTauDetails panTauDetail, int &value) const
Get and set values of pantau details variables via enum.
Definition: TauJet_v3.cxx:367
read_hist_ntuple.t
t
Definition: read_hist_ntuple.py:5
TauAnalysisTools::finalTauPt
double finalTauPt(const xAOD::TauJet &xTau)
return MVA based tau pt in GeV
Definition: PhysicsAnalysis/TauID/TauAnalysisTools/Root/HelperFunctions.cxx:113
ATH_MSG_VERBOSE
#define ATH_MSG_VERBOSE(x)
Definition: AthMsgStreamMacros.h:28
TauAnalysisTools::TruthJet
@ TruthJet
Definition: PhysicsAnalysis/TauID/TauAnalysisTools/TauAnalysisTools/Enums.h:101
dbg::ptr
void * ptr(T *p)
Definition: SGImplSvc.cxx:74
CP::SystematicVariation
Definition: SystematicVariation.h:47
TauAnalysisTools::CommonEfficiencyTool::getValueTH3
static CP::CorrectionCode getValueTH3(const TObject *oObject, double &dEfficiencyScaleFactor, double dVars[])
Definition: CommonEfficiencyTool.cxx:766
TauAnalysisTools::CommonEfficiencyTool::tTupleObjectFunc
std::tuple< TObject *, CP::CorrectionCode(*)(const TObject *oObject, double &dEfficiencyScaleFactor, double dVars[]) > tTupleObjectFunc
Definition: CommonEfficiencyTool.h:87
TauAnalysisTools::TruthLeadPt
double TruthLeadPt(const xAOD::DiTauJet &xDiTau)
return the truth vis pT of the leading pT matched particle.
Definition: CommonDiTauEfficiencyTool.cxx:326
mapkey::sys
@ sys
Definition: TElectronEfficiencyCorrectionTool.cxx:42
xAOD::EventInfo_v1::runNumber
uint32_t runNumber() const
The current event's run number.
python.iconfTool.models.loaders.level
level
Definition: loaders.py:20
xAOD::TauJetParameters::Mode_1pXn
@ Mode_1pXn
Definition: TauDefs.h:388
SG::VarHandleKeyArray::setOwner
virtual void setOwner(IDataHandleHolder *o)=0
IDTPMcnv.htype
htype
Definition: IDTPMcnv.py:27
AthCommonDataStore< AthCommonMsg< AlgTool > >::evtStore
ServiceHandle< StoreGateSvc > & evtStore()
The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.
Definition: AthCommonDataStore.h:85
TauAnalysisTools::CommonEfficiencyTool::m_mSystematicSets
std::unordered_map< CP::SystematicSet, std::string > m_mSystematicSets
Definition: CommonEfficiencyTool.h:98
TauAnalysisTools::finalTauP
double finalTauP(const xAOD::TauJet &xTau)
return MVA based tau P in GeV
Definition: PhysicsAnalysis/TauID/TauAnalysisTools/Root/HelperFunctions.cxx:134
TauAnalysisTools::CommonEfficiencyTool::getValue
virtual CP::CorrectionCode getValue(const std::string &sHistName, const xAOD::TauJet &xTau, double &dEfficiencyScaleFactor) const
Definition: CommonEfficiencyTool.cxx:669
CP::CorrectionCode::OutOfValidityRange
@ OutOfValidityRange
Input object is out of validity range.
Definition: CorrectionCode.h:37
TauAnalysisTools::truthDecayMode
double truthDecayMode(const xAOD::TauJet &xTau)
return truth decay mode (if hadronic truth tau match)
Definition: PhysicsAnalysis/TauID/TauAnalysisTools/Root/HelperFunctions.cxx:186
CP::CorrectionCode::Error
@ Error
Some error happened during the object correction.
Definition: CorrectionCode.h:36
python.utils.AtlRunQueryDQUtils.p
p
Definition: AtlRunQueryDQUtils.py:210
AthCommonDataStore
Definition: AthCommonDataStore.h:52
ATH_MSG_ERROR
#define ATH_MSG_ERROR(x)
Definition: AthMsgStreamMacros.h:33
SG::Decorator
Helper class to provide type-safe access to aux data.
Definition: Decorator.h:59
TauAnalysisTools::CommonEfficiencyTool::m_sWP
std::string m_sWP
Definition: CommonEfficiencyTool.h:139
ATH_MSG_DEBUG
#define ATH_MSG_DEBUG(x)
Definition: AthMsgStreamMacros.h:29
TauAnalysisTools::CommonEfficiencyTool::m_mSF
std::unique_ptr< tSFMAP > m_mSF
Definition: CommonEfficiencyTool.h:95
xAOD::TauJet_v3
Class describing a tau jet.
Definition: TauJet_v3.h:41
Amg::transform
Amg::Vector3D transform(Amg::Vector3D &v, Amg::Transform3D &tr)
Transform a point from a Trasformation3D.
Definition: GeoPrimitivesHelpers.h:156
TauAnalysisTools::CommonDiTauEfficiencyTool::ReadInputs
void ReadInputs(std::unique_ptr< TFile > &fFile)
Definition: CommonDiTauEfficiencyTool.cxx:192
CP::SystematicSet::end
const_iterator end() const
description: const iterator to the end of the set
Definition: SystematicSet.h:59
python.xAODType.dummy
dummy
Definition: xAODType.py:4
MSG::name
const std::string & name(Level lvl)
Convenience function for translating message levels to strings.
Definition: MsgLevel.cxx:19
TauAnalysisTools::CommonEfficiencyTool::m_fX
std::function< double(const xAOD::TauJet &xTau)> m_fX
Definition: CommonEfficiencyTool.h:103
AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore
StoreGateSvc_t m_detStore
Pointer to StoreGate (detector store by default)
Definition: AthCommonDataStore.h:393
TauAnalysisTools::CommonEfficiencyTool::m_bSFIsAvailable
bool m_bSFIsAvailable
Definition: CommonEfficiencyTool.h:150
TauAnalysisTools::CommonEfficiencyTool::getValueTH2
static CP::CorrectionCode getValueTH2(const TObject *oObject, double &dEfficiencyScaleFactor, double dVars[])
Definition: CommonEfficiencyTool.cxx:734
xAOD::TauJetParameters::Mode_3p0n
@ Mode_3p0n
Definition: TauDefs.h:389
TauAnalysisTools::CommonEfficiencyTool::m_sSFHistName
std::string m_sSFHistName
Definition: CommonEfficiencyTool.h:141
min
#define min(a, b)
Definition: cfImp.cxx:40
TauAnalysisTools::TruthMuon
@ TruthMuon
Definition: PhysicsAnalysis/TauID/TauAnalysisTools/TauAnalysisTools/Enums.h:99
SG::VarHandleKeyArray::renounce
virtual void renounce()=0
SG::HandleClassifier::type
std::conditional< std::is_base_of< SG::VarHandleKeyArray, T >::value, VarHandleKeyArrayType, type2 >::type type
Definition: HandleClassifier.h:54
TauAnalysisTools::CommonEfficiencyTool::m_sInputFilePath
std::string m_sInputFilePath
Definition: CommonEfficiencyTool.h:137
xAOD::EventInfo_v1::averageInteractionsPerCrossing
float averageInteractionsPerCrossing() const
Average interactions per crossing for all BCIDs - for out-of-time pile-up.
Definition: EventInfo_v1.cxx:397
TauAnalysisTools::TruthDeltaR
double TruthDeltaR(const xAOD::DiTauJet &xDiTau)
return the dR of between the leading and subleading pT matched particle.
Definition: CommonDiTauEfficiencyTool.cxx:342
TauAnalysisTools::CommonEfficiencyTool::m_sAffectingSystematics
CP::SystematicSet m_sAffectingSystematics
Definition: CommonEfficiencyTool.h:134
merge_scale_histograms.doc
string doc
Definition: merge_scale_histograms.py:9
TauAnalysisTools::CommonEfficiencyTool::m_bUseTauSubstructure
bool m_bUseTauSubstructure
Definition: CommonEfficiencyTool.h:144
CP::SystematicSet::insert
void insert(const SystematicVariation &systematic)
description: insert a systematic into the set
Definition: SystematicSet.cxx:88
name
std::string name
Definition: Control/AthContainers/Root/debug.cxx:221
TauAnalysisTools::CommonEfficiencyTool::m_fY
std::function< double(const xAOD::TauJet &xTau)> m_fY
Definition: CommonEfficiencyTool.h:104
TauAnalysisTools::CommonEfficiencyTool::ConvertMuToString
std::string ConvertMuToString(const int iMu) const
Definition: CommonEfficiencyTool.cxx:403
xAOD::EventInfo_v1
Class describing the basic event information.
Definition: EventInfo_v1.h:43
TauAnalysisTools::CommonDiTauEfficiencyTool::m_bSFIsAvailableDiTau
bool m_bSFIsAvailableDiTau
true if scale factor name is already decorated
Definition: CommonDiTauEfficiencyTool.h:106
TauAnalysisTools::CommonEfficiencyTool::addHistogramToSFMap
void addHistogramToSFMap(TKey *kKey, const std::string &sKeyName)
Definition: CommonEfficiencyTool.cxx:555
PathResolverFindCalibFile
std::string PathResolverFindCalibFile(const std::string &logical_file_name)
Definition: PathResolver.cxx:431
CP::SystematicSet::find
iterator find(const SystematicVariation &sys) const
description: find an element in the set
Definition: SystematicSet.h:63
TauAnalysisTools::CommonDiTauEfficiencyTool::m_fZDiTau
double(* m_fZDiTau)(const xAOD::DiTauJet &xDiTau)
scale factor bin z (e.g.
Definition: CommonDiTauEfficiencyTool.h:86
CP::CorrectionCode::Ok
@ Ok
The correction was done successfully.
Definition: CorrectionCode.h:38
a
TList * a
Definition: liststreamerinfos.cxx:10
h
TauAnalysisTools::finalTauAbsEta
double finalTauAbsEta(const xAOD::TauJet &xTau)
return MVA based absolute tau eta
Definition: PhysicsAnalysis/TauID/TauAnalysisTools/Root/HelperFunctions.cxx:127
TauAnalysisTools::CommonEfficiencyTool::m_bSFIsAvailableChecked
bool m_bSFIsAvailableChecked
Definition: CommonEfficiencyTool.h:151
ATH_MSG_WARNING
#define ATH_MSG_WARNING(x)
Definition: AthMsgStreamMacros.h:32
TauAnalysisTools::CommonEfficiencyTool::CommonEfficiencyTool
CommonEfficiencyTool(const std::string &sName)
Create a proper constructor for Athena.
Definition: CommonEfficiencyTool.cxx:74
TauAnalysisTools::tauLeadTrackEta
double tauLeadTrackEta(const xAOD::TauJet &xTau)
return leading charge tau track eta
Definition: PhysicsAnalysis/TauID/TauAnalysisTools/Root/HelperFunctions.cxx:141
TauAnalysisTools::CommonEfficiencyTool::ConvertDecayModeToString
std::string ConvertDecayModeToString(const int iDecayMode) const
Definition: CommonEfficiencyTool.cxx:415
TauAnalysisTools::getTruthParticleType
TruthMatchedParticleType getTruthParticleType(const xAOD::TauJet &xTau)
return TauJet match type
Definition: PhysicsAnalysis/TauID/TauAnalysisTools/Root/HelperFunctions.cxx:572
AthCommonMsg< AlgTool >::msg
MsgStream & msg() const
Definition: AthCommonMsg.h:24
CP::CorrectionCode
Return value from object correction CP tools.
Definition: CorrectionCode.h:31
SG::VarHandleBase::vhKey
SG::VarHandleKey & vhKey()
Return a non-const reference to the HandleKey.
Definition: StoreGate/src/VarHandleBase.cxx:623
xAOD::TauJetParameters::Mode_3pXn
@ Mode_3pXn
Definition: TauDefs.h:390
python.Bindings.keys
keys
Definition: Control/AthenaPython/python/Bindings.py:798
TauAnalysisTools::CommonDiTauEfficiencyTool::m_bSFIsAvailableCheckedDiTau
bool m_bSFIsAvailableCheckedDiTau
true if cale factor name is already decorated has already been checked
Definition: CommonDiTauEfficiencyTool.h:108
TauAnalysisTools::CommonEfficiencyTool::m_bNoMultiprong
bool m_bNoMultiprong
Definition: CommonEfficiencyTool.h:143
TauGNNUtils::Variables::Track::dEta
bool dEta(const xAOD::TauJet &tau, const xAOD::TauTrack &track, double &out)
Definition: TauGNNUtils.cxx:527
TauAnalysisTools::CommonEfficiencyTool::m_sInputFileName
std::string m_sInputFileName
Definition: CommonEfficiencyTool.h:138
SG::DataProxy
Definition: DataProxy.h:44
TauAnalysisTools::CommonEfficiencyTool::affectingSystematics
virtual CP::SystematicSet affectingSystematics() const
returns: the list of all systematics this tool can be affected by
Definition: CommonEfficiencyTool.cxx:317
TauAnalysisTools::CommonEfficiencyTool::m_bSplitMu
bool m_bSplitMu
Definition: CommonEfficiencyTool.h:152
AthCommonDataStore::declareGaudiProperty
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>
Definition: AthCommonDataStore.h:156
TauAnalysisTools::CommonEfficiencyTool::m_mSystematicsHistNames
std::map< std::string, std::string > m_mSystematicsHistNames
Definition: CommonEfficiencyTool.h:101
TauAnalysisTools::CommonEfficiencyTool::getValueTH1
static CP::CorrectionCode getValueTH1(const TObject *oObject, double &dEfficiencyScaleFactor, double dVars[])
Definition: CommonEfficiencyTool.cxx:704
TauAnalysisTools::TruthHadronicDiTau
@ TruthHadronicDiTau
Definition: PhysicsAnalysis/TauID/TauAnalysisTools/TauAnalysisTools/Enums.h:102
TauAnalysisTools::CommonEfficiencyTool::getValueTF1
static CP::CorrectionCode getValueTF1(const TObject *oObject, double &dEfficiencyScaleFactor, double dVars[])
Definition: CommonEfficiencyTool.cxx:801
TauAnalysisTools::CommonEfficiencyTool::m_sRecommendedSystematics
CP::SystematicSet m_sRecommendedSystematics
Definition: CommonEfficiencyTool.h:135
TauAnalysisTools::finalTauEta
double finalTauEta(const xAOD::TauJet &xTau)
return MVA based tau eta
Definition: PhysicsAnalysis/TauID/TauAnalysisTools/Root/HelperFunctions.cxx:120
fitman.k
k
Definition: fitman.py:528
xAOD::EventInfo_v1::actualInteractionsPerCrossing
float actualInteractionsPerCrossing() const
Average interactions per crossing for the current BCID - for in-time pile-up.
Definition: EventInfo_v1.cxx:380
TauAnalysisTools::CommonDiTauEfficiencyTool::getValue
CP::CorrectionCode getValue(const std::string &sHistName, const xAOD::DiTauJet &xDiTau, double &dEfficiencyScaleFactor) const
Get the scale factor from a particular recommendations histogram.
Definition: CommonDiTauEfficiencyTool.cxx:296
TauAnalysisTools::split
void split(const std::string &sInput, const char cDelim, std::vector< std::string > &vOut)
Definition: PhysicsAnalysis/TauID/TauAnalysisTools/Root/HelperFunctions.cxx:23