Loading [MathJax]/extensions/tex2jax.js
 |
ATLAS Offline Software
|
Implementation of tool to decorate flattened BDT score and working points.
More...
#include <TauWPDecorator.h>
|
| TauWPDecorator (const std::string &name="TauWPDecorator") |
| Constructor. More...
|
|
| ~TauWPDecorator () |
| Destructor. More...
|
|
virtual StatusCode | initialize () override |
| Initialization of this tool. More...
|
|
virtual StatusCode | execute (xAOD::TauJet &tau) const override |
| Executation of this tool. More...
|
|
virtual StatusCode | eventInitialize () override |
| Event initializer - called at the beginning of each event. More...
|
|
virtual StatusCode | executePi0CreateROI (xAOD::TauJet &pTau, CaloConstCellContainer &caloCellContainer, boost::dynamic_bitset<> &map) const override |
|
virtual StatusCode | executeVertexFinder (xAOD::TauJet &pTau, const xAOD::VertexContainer *vertexContainer=nullptr) const override |
|
virtual StatusCode | executeTrackFinder (xAOD::TauJet &pTau, xAOD::TauTrackContainer &tauTrackContainer) const override |
|
virtual StatusCode | executeTrackClassifier (xAOD::TauJet &pTau, xAOD::TauTrackContainer &tauTrackContainer) const override |
|
virtual StatusCode | executeShotFinder (xAOD::TauJet &pTau, xAOD::CaloClusterContainer &shotClusterContainer, xAOD::PFOContainer &PFOContainer) const override |
|
virtual StatusCode | executePi0ClusterCreator (xAOD::TauJet &pTau, xAOD::PFOContainer &neutralPFOContainer, xAOD::PFOContainer &hadronicPFOContainer, const xAOD::CaloClusterContainer &pCaloClusterContainer) const override |
|
virtual StatusCode | executeVertexVariables (xAOD::TauJet &pTau, xAOD::VertexContainer &vertexContainer) const override |
|
virtual StatusCode | executePi0ClusterScaler (xAOD::TauJet &pTau, xAOD::PFOContainer &neutralPFOContainer, xAOD::PFOContainer &chargedPFOContainer) const override |
|
virtual StatusCode | executePi0nPFO (xAOD::TauJet &pTau, xAOD::PFOContainer &neutralPFOContainer) const override |
|
virtual StatusCode | executePanTau (xAOD::TauJet &pTau, xAOD::ParticleContainer &particleContainer, xAOD::PFOContainer &neutralPFOContainer) const override |
|
virtual StatusCode | eventFinalize () override |
| Event finalizer - called at the end of each event. More...
|
|
virtual StatusCode | finalize () override |
| Finalizer. More...
|
|
std::string | find_file (const std::string &fname) const |
|
virtual StatusCode | readConfig () override |
|
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 |
|
|
Gaudi::Property< bool > | m_in_trigger {this, "inTrigger", false, "Indicate if the tool is running on trigger"} |
|
Gaudi::Property< bool > | m_in_AOD {this, "inAOD", false, "Indicate if the tool is running on AOD"} |
|
Gaudi::Property< bool > | m_in_EleRM {this, "inEleRM", false, "Indicate if the tool is running on EleRM routine"} |
|
Gaudi::Property< std::string > | m_tauRecToolsTag {this, "calibFolder", "tauRecTools/R22_preprod", "CVMFS path to the tau calibration folder"} |
|
|
Gaudi::Property< bool > | m_useAbsEta {this, "UseAbsEta", false} |
|
Gaudi::Property< bool > | m_defineWPs {this, "DefineWPs", false} |
|
Gaudi::Property< std::string > | m_scoreName {this, "ScoreName", ""} |
|
Gaudi::Property< std::string > | m_scoreNameTrans {this, "NewScoreName", ""} |
|
Gaudi::Property< std::string > | m_file0p {this, "flatteningFile0Prong", ""} |
|
Gaudi::Property< std::string > | m_file1p {this, "flatteningFile1Prong", ""} |
|
Gaudi::Property< std::string > | m_file2p {this, "flatteningFile2Prong", ""} |
|
Gaudi::Property< std::string > | m_file3p {this, "flatteningFile3Prong", ""} |
|
Gaudi::Property< std::vector< int > > | m_EDMWPs {this, "CutEnumVals", {}} |
|
Gaudi::Property< std::vector< float > > | m_EDMWPEffs0p {this, "SigEff0P", {}} |
|
Gaudi::Property< std::vector< float > > | m_EDMWPEffs1p {this, "SigEff1P", {}} |
|
Gaudi::Property< std::vector< float > > | m_EDMWPEffs2p {this, "SigEff2P", {}} |
|
Gaudi::Property< std::vector< float > > | m_EDMWPEffs3p {this, "SigEff3P", {}} |
|
Gaudi::Property< std::vector< std::string > > | m_decorWPs {this, "DecorWPNames", {}} |
|
Gaudi::Property< std::vector< float > > | m_decorWPEffs0p {this, "DecorWPCutEffs0P", {}} |
|
Gaudi::Property< std::vector< float > > | m_decorWPEffs1p {this, "DecorWPCutEffs1P", {}} |
|
Gaudi::Property< std::vector< float > > | m_decorWPEffs2p {this, "DecorWPCutEffs2P", {}} |
|
Gaudi::Property< std::vector< float > > | m_decorWPEffs3p {this, "DecorWPCutEffs3P", {}} |
|
std::vector< SG::Accessor< char > > | m_charDecors |
|
SG::ReadDecorHandleKey< xAOD::EventInfo > | m_aveIntPerXKey |
|
std::shared_ptr< std::vector< m_pair_t > > | m_hists0p |
| Efficiency and corresponding score distributions of 0-prong taus. More...
|
|
std::shared_ptr< std::vector< m_pair_t > > | m_hists1p |
| Efficiency and corresponding score distributions of 1-prong taus. More...
|
|
std::shared_ptr< std::vector< m_pair_t > > | m_hists2p |
| Efficiency and corresponding score distributions of 2-prong taus. More...
|
|
std::shared_ptr< std::vector< m_pair_t > > | m_hists3p |
| Efficiency and corresponding score distributions of 3-prong taus. More...
|
|
std::map< int, double > | m_xMin |
| Map of n-prong and the minimum value of x variables. More...
|
|
std::map< int, double > | m_yMin |
| Map of n-prong and the minimum value of y variables. More...
|
|
std::map< int, double > | m_xMax |
| Map of n-prong and the maximum value of x variables. More...
|
|
std::map< int, double > | m_yMax |
| Map of n-prong and the maximum value of y variables. More...
|
|
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 |
|
Implementation of tool to decorate flattened BDT score and working points.
Input comes from ROOT files with lists of TH2s containing BDT/RNN score distributions as a function of the dependent variables. For eVeto, the score distributions depend on tau pT and |eta| of the leading track. Otherwise, the score distributions depend on tau pT and pileup.
- Author
- P.O. DeViveiros
-
W. Davey
-
L. Hauswald
Definition at line 35 of file TauWPDecorator.h.
◆ m_pair_t
◆ StoreGateSvc_t
◆ TauWPDecorator()
TauWPDecorator::TauWPDecorator |
( |
const std::string & |
name = "TauWPDecorator" | ) |
|
◆ ~TauWPDecorator()
TauWPDecorator::~TauWPDecorator |
( |
| ) |
|
◆ 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);
253 return PBASE::declareProperty(
name,hndl,
doc);
◆ 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.
229 return PBASE::declareProperty(
name,hndl,
doc);
◆ 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.
338 return PBASE::declareProperty(
name, property,
doc);
◆ 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()
◆ eventFinalize()
StatusCode TauRecToolBase::eventFinalize |
( |
| ) |
|
|
overridevirtualinherited |
◆ eventInitialize()
StatusCode TauRecToolBase::eventInitialize |
( |
| ) |
|
|
overridevirtualinherited |
◆ evtStore() [1/2]
◆ evtStore() [2/2]
◆ execute()
StatusCode TauWPDecorator::execute |
( |
xAOD::TauJet & |
tau | ) |
const |
|
overridevirtual |
Executation of this tool.
Reimplemented from TauRecToolBase.
Definition at line 171 of file TauWPDecorator.cxx.
174 double xVariable = tau.
pt();
177 double yVariable = 0.0;
180 yVariable = std::abs(acc_absEta(tau));
184 if (!eventInfoDecorHandle.isPresent()) {
186 return StatusCode::FAILURE;
188 yVariable = eventInfoDecorHandle(0);
192 int nProng = nTracks;
203 else if (nTracks > 2) {
227 std::shared_ptr<std::vector<m_pair_t>> histArray =
nullptr;
229 else if (nProng == 1) histArray =
m_hists1p;
230 else if (nProng == 2) histArray =
m_hists2p;
233 std::array<double, 2>
cuts = {-1.01, 1.01};
234 std::array<double, 2> effs = {1.0, 0.0};
236 bool gotHigh =
false;
239 double score = acc_score(tau);
242 for (
unsigned int i = 0;
i < histArray->size(); ++
i) {
243 std::shared_ptr<TH2> myHist = histArray->at(
i).second;
244 double myCut = myHist->Interpolate(xVariable, yVariable);
248 effs[0] = histArray->at(
i).first;
251 else if (myCut >
score && ((!gotHigh) || std::abs(myCut-
score) < std::abs(
cuts[1]-
score))) {
253 effs[1] = histArray->at(
i).first;
257 if (gotLow && gotHigh){
263 double scoreTrans = -1111.;
265 scoreTrans = 1 - effs[1];
268 scoreTrans = 1 - effs[0];
275 acc_scoreTrans(tau) = scoreTrans;
279 for (
size_t wpIndex=0; wpIndex <
m_EDMWPs.size(); ++wpIndex) {
283 else if(nProng == 1) {
286 else if(nProng == 2) {
294 for (
size_t wpIndex=0; wpIndex <
m_decorWPs.size(); ++wpIndex) {
300 else if(nProng == 1) {
303 else if(nProng == 2) {
312 return StatusCode::SUCCESS;
◆ executePanTau()
◆ executePi0ClusterCreator()
◆ executePi0ClusterScaler()
◆ executePi0CreateROI()
◆ executePi0nPFO()
◆ executeShotFinder()
◆ executeTrackClassifier()
◆ executeTrackFinder()
◆ executeVertexFinder()
◆ executeVertexVariables()
◆ 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
◆ finalize()
StatusCode TauRecToolBase::finalize |
( |
| ) |
|
|
overridevirtualinherited |
◆ find_file()
std::string TauRecToolBase::find_file |
( |
const std::string & |
fname | ) |
const |
|
inherited |
◆ 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
127 return (
proxy ==
nullptr ? 0 :
proxy->sgkey() );
128 #endif // XAOD_STANDALONE
◆ getName()
const std::string & asg::AsgTool::getName |
( |
const void * |
ptr | ) |
const |
|
inherited |
Get the name of an object that is / should be in the event store.
This is a bit of a special one. StoreGateSvc
and xAOD::TEvent
both provide ways for getting the std::string
name for an object that is in the store, based on a bare pointer. But they provide different interfaces for doing so.
In order to allow tools to efficiently perform this operation, they can use this helper function.
- See also
- asg::AsgTool::getKey
- Parameters
-
ptr | The bare pointer to the object that the event store should know about |
- Returns
- The string name of the object in the store. If not found, an empty string.
Definition at line 106 of file AsgTool.cxx.
108 #ifdef XAOD_STANDALONE
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.
◆ inAOD()
bool TauRecToolBase::inAOD |
( |
| ) |
const |
|
inlineprotectedinherited |
◆ inEleRM()
bool TauRecToolBase::inEleRM |
( |
| ) |
const |
|
inlineprotectedinherited |
◆ initialize()
StatusCode TauWPDecorator::initialize |
( |
| ) |
|
|
overridevirtual |
Initialization of this tool.
Reimplemented from TauRecToolBase.
Definition at line 131 of file TauWPDecorator.cxx.
138 return StatusCode::FAILURE;
143 m_hists0p = std::make_shared<std::vector<m_pair_t>>();
148 m_hists1p = std::make_shared<std::vector<m_pair_t>>();
154 m_hists2p = std::make_shared<std::vector<m_pair_t>>();
159 m_hists3p = std::make_shared<std::vector<m_pair_t>>();
163 for (
size_t wpIndex=0; wpIndex <
m_decorWPs.size(); ++wpIndex) {
167 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.
◆ inTrigger()
bool TauRecToolBase::inTrigger |
( |
| ) |
const |
|
inlineprotectedinherited |
◆ 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()
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, JetDumper, JetBottomUpSoftDrop, JetRecursiveSoftDrop, JetSoftDrop, JetConstituentsRetriever, JetSubStructureMomentToolsBase, JetSplitter, JetToolRunner, JetPruner, JetPseudojetRetriever, JetTrimmer, AsgHelloTool, and KtDeltaRTool.
Definition at line 131 of file AsgTool.cxx.
◆ readConfig()
StatusCode TauRecToolBase::readConfig |
( |
| ) |
|
|
overridevirtualinherited |
Implements ITauToolBase.
Definition at line 27 of file TauRecToolBase.cxx.
32 #ifndef XAOD_STANDALONE
33 bool configPathDeclared =
false;
34 for (Gaudi::Details::PropertyBase* property : getProperties())
36 if (property->name() ==
"ConfigPath")
38 configPathDeclared =
true;
42 if (!configPathDeclared)
43 #elif defined(XAOD_STANDALONE)
44 PropertyMgr::PropMap_t property_map = getPropertyMgr()->getProperties();
45 if (property_map.find(
"ConfigPath") == property_map.end())
47 # error "What environment are we in?!?"
48 #endif // XAOD_STANDALONE
50 ATH_MSG_INFO(
"No config file path property declared yet, this is not recommended");
51 return StatusCode::SUCCESS;
55 const std::string* config_file_path_property;
58 config_file_path_property = getProperty<std::string>(
"ConfigPath");
59 std::string config_file_path =
find_file(*config_file_path_property);
63 THashList* lList =
env.GetTable();
64 for( Int_t
i = 0; lList &&
i < lList->GetEntries(); ++
i )
68 #ifndef XAOD_STANDALONE
70 const std::type_info*
type =
getProperty(lList->At(
i )->GetName()).type_info();
73 if (*
type ==
typeid(
bool))
74 sc = this->setProperty(lList->At(
i )->GetName(),
75 bool(
env.GetValue(lList->At(
i )->GetName(),
bool(
true))));
76 else if (*
type ==
typeid(
int))
77 sc = this->setProperty(lList->At(
i )->GetName(),
78 env.GetValue(lList->At(
i )->GetName(),
int(0)));
79 else if (*
type ==
typeid(
float))
80 sc = this->setProperty(lList->At(
i )->GetName(),
81 env.GetValue(lList->At(
i )->GetName(),
float(0)));
82 else if (*
type ==
typeid(
double))
83 sc = this->setProperty(lList->At(
i )->GetName(),
84 env.GetValue(lList->At(
i )->GetName(),
double(0)));
85 else if (*
type ==
typeid(std::string))
86 sc = this->setProperty(lList->At(
i )->GetName(),
87 env.GetValue(lList->At(
i )->GetName(),
""));
93 sc = this->setProperty(lList->At(
i )->GetName(),
94 bool(
env.GetValue(lList->At(
i )->GetName(),
bool(
true))));
96 sc = this->setProperty(lList->At(
i )->GetName(),
97 env.GetValue(lList->At(
i )->GetName(),
int(0)));
99 sc = this->setProperty(lList->At(
i )->GetName(),
100 env.GetValue(lList->At(
i )->GetName(),
float(0)));
102 sc = this->setProperty(lList->At(
i )->GetName(),
103 env.GetValue(lList->At(
i )->GetName(),
double(0)));
105 sc = this->setProperty(lList->At(
i )->GetName(),
106 env.GetValue(lList->At(
i )->GetName(),
""));
107 #endif // XAOD_STANDALONE
110 #ifndef XAOD_STANDALONE
111 ATH_MSG_FATAL(
"there was a problem to find the correct type enum: "<<
type->name());
114 #endif // XAOD_STANDALONE
115 return StatusCode::FAILURE;
117 if (!
sc.isSuccess()) {
118 ATH_MSG_FATAL(
"failed to set property: " << lList->At(
i )->GetName());
119 return StatusCode::FAILURE;
122 return StatusCode::SUCCESS;
◆ renounce()
◆ renounceArray()
◆ retrieveHistos()
StatusCode TauWPDecorator::retrieveHistos |
( |
int |
nProng | ) |
|
|
private |
Retrieve the histograms containing BDT/RNN score distributions as a function of dependent variables.
- Parameters
-
nProng | Prong of the tau candidate |
Definition at line 24 of file TauWPDecorator.cxx.
27 std::shared_ptr<std::vector<m_pair_t>> histArray =
nullptr;
32 else if (nProng == 1) {
36 else if (nProng == 2) {
40 else if (nProng == 3) {
46 return StatusCode::FAILURE;
50 std::unique_ptr<TFile>
file(TFile::Open(
fullPath.c_str(),
"READ"));
54 return StatusCode::FAILURE;
60 for (
int i = 0;
i < 100; ++
i) {
62 TH2* graph =
dynamic_cast<TH2*
>(
file->Get(Form(
"h2_%02d",
i)));
64 ATH_MSG_ERROR(
"Failed to retrieve Graph " <<
i <<
" named " << Form(
"h2_%02d",
i));
65 return StatusCode::FAILURE;
67 graph->SetDirectory(
nullptr);
68 std::shared_ptr<TH2> sharedGraph(graph);
69 histArray->push_back(
m_pair_t(
float(
i)/100., std::move(sharedGraph)));
74 return StatusCode::SUCCESS;
◆ storeLimits()
StatusCode TauWPDecorator::storeLimits |
( |
int |
nProng | ) |
|
|
private |
Obtain the limit of the dependent variables.
- Parameters
-
nProng | Prong of the tau candidate |
Definition at line 78 of file TauWPDecorator.cxx.
79 std::shared_ptr<std::vector<m_pair_t>> histArray =
nullptr;
83 else if (nProng == 1) {
86 else if (nProng == 2) {
89 else if (nProng == 3) {
94 return StatusCode::FAILURE;
97 std::shared_ptr<TH2> firstHist = histArray->at(0).second;
98 m_xMin[nProng] = firstHist->GetXaxis()->GetXmin();
99 m_xMax[nProng] = firstHist->GetXaxis()->GetBinCenter(firstHist->GetNbinsX());
100 m_yMin[nProng] = firstHist->GetYaxis()->GetXmin();
101 m_yMax[nProng] = firstHist->GetYaxis()->GetBinCenter(firstHist->GetNbinsY());
104 for (
size_t i = 1;
i < histArray->size(); ++
i) {
105 std::shared_ptr<TH2>
hist = histArray->at(
i).second;
107 double xMin =
hist->GetXaxis()->GetXmin();
108 double xMax =
hist->GetXaxis()->GetBinCenter(firstHist->GetNbinsX());
109 double yMin =
hist->GetYaxis()->GetXmin();
110 double yMax =
hist->GetYaxis()->GetBinCenter(firstHist->GetNbinsY());
112 if (std::abs(
m_xMin[nProng] - xMin) > 1
e-5 ||
113 std::abs(
m_xMax[nProng] - xMax) > 1
e-5 ||
114 std::abs(
m_yMin[nProng] - yMin) > 1
e-5 ||
115 std::abs(
m_yMax[nProng] - yMax) > 1
e-5) {
120 return StatusCode::SUCCESS;
◆ 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.
◆ transformScore()
double TauWPDecorator::transformScore |
( |
double |
score, |
|
|
double |
cutLow, |
|
|
double |
effLow, |
|
|
double |
cutHigh, |
|
|
double |
effHigh |
|
) |
| const |
|
private |
Obtain the flattened score.
- Parameters
-
score | Original BDT/RNN score |
cutLow | Lower score cut |
effLow | Efficiency of the lower cut |
cutHigh | Higher score cut |
effHigh | Efficiency of the higher cut |
Definition at line 124 of file TauWPDecorator.cxx.
125 double efficiency = effLow + (
score - cutLow)/(cutHigh - cutLow) * (effHigh - effLow);
◆ updateVHKA()
◆ m_aveIntPerXKey
Initial value:{this,
"averageInteractionsPerCrossingKey",
"EventInfo.averageInteractionsPerCrossing",
"Decoration for Average Interaction Per Crossing"}
Definition at line 98 of file TauWPDecorator.h.
◆ m_charDecors
std::vector<SG::Accessor<char> > TauWPDecorator::m_charDecors |
|
private |
◆ m_decorWPEffs0p
Gaudi::Property<std::vector<float> > TauWPDecorator::m_decorWPEffs0p {this, "DecorWPCutEffs0P", {}} |
|
private |
◆ m_decorWPEffs1p
Gaudi::Property<std::vector<float> > TauWPDecorator::m_decorWPEffs1p {this, "DecorWPCutEffs1P", {}} |
|
private |
◆ m_decorWPEffs2p
Gaudi::Property<std::vector<float> > TauWPDecorator::m_decorWPEffs2p {this, "DecorWPCutEffs2P", {}} |
|
private |
◆ m_decorWPEffs3p
Gaudi::Property<std::vector<float> > TauWPDecorator::m_decorWPEffs3p {this, "DecorWPCutEffs3P", {}} |
|
private |
◆ m_decorWPs
Gaudi::Property<std::vector<std::string> > TauWPDecorator::m_decorWPs {this, "DecorWPNames", {}} |
|
private |
◆ m_defineWPs
Gaudi::Property<bool> TauWPDecorator::m_defineWPs {this, "DefineWPs", false} |
|
private |
◆ m_detStore
◆ m_EDMWPEffs0p
Gaudi::Property<std::vector<float> > TauWPDecorator::m_EDMWPEffs0p {this, "SigEff0P", {}} |
|
private |
◆ m_EDMWPEffs1p
Gaudi::Property<std::vector<float> > TauWPDecorator::m_EDMWPEffs1p {this, "SigEff1P", {}} |
|
private |
◆ m_EDMWPEffs2p
Gaudi::Property<std::vector<float> > TauWPDecorator::m_EDMWPEffs2p {this, "SigEff2P", {}} |
|
private |
◆ m_EDMWPEffs3p
Gaudi::Property<std::vector<float> > TauWPDecorator::m_EDMWPEffs3p {this, "SigEff3P", {}} |
|
private |
◆ m_EDMWPs
Gaudi::Property<std::vector<int> > TauWPDecorator::m_EDMWPs {this, "CutEnumVals", {}} |
|
private |
◆ m_evtStore
◆ m_file0p
Gaudi::Property<std::string> TauWPDecorator::m_file0p {this, "flatteningFile0Prong", ""} |
|
private |
◆ m_file1p
Gaudi::Property<std::string> TauWPDecorator::m_file1p {this, "flatteningFile1Prong", ""} |
|
private |
◆ m_file2p
Gaudi::Property<std::string> TauWPDecorator::m_file2p {this, "flatteningFile2Prong", ""} |
|
private |
◆ m_file3p
Gaudi::Property<std::string> TauWPDecorator::m_file3p {this, "flatteningFile3Prong", ""} |
|
private |
◆ m_hists0p
std::shared_ptr<std::vector<m_pair_t> > TauWPDecorator::m_hists0p |
|
private |
Efficiency and corresponding score distributions of 0-prong taus.
Definition at line 105 of file TauWPDecorator.h.
◆ m_hists1p
std::shared_ptr<std::vector<m_pair_t> > TauWPDecorator::m_hists1p |
|
private |
Efficiency and corresponding score distributions of 1-prong taus.
Definition at line 106 of file TauWPDecorator.h.
◆ m_hists2p
std::shared_ptr<std::vector<m_pair_t> > TauWPDecorator::m_hists2p |
|
private |
Efficiency and corresponding score distributions of 2-prong taus.
Definition at line 107 of file TauWPDecorator.h.
◆ m_hists3p
std::shared_ptr<std::vector<m_pair_t> > TauWPDecorator::m_hists3p |
|
private |
Efficiency and corresponding score distributions of 3-prong taus.
Definition at line 108 of file TauWPDecorator.h.
◆ m_in_AOD
Gaudi::Property<bool> TauRecToolBase::m_in_AOD {this, "inAOD", false, "Indicate if the tool is running on AOD"} |
|
protectedinherited |
◆ m_in_EleRM
Gaudi::Property<bool> TauRecToolBase::m_in_EleRM {this, "inEleRM", false, "Indicate if the tool is running on EleRM routine"} |
|
protectedinherited |
◆ m_in_trigger
Gaudi::Property<bool> TauRecToolBase::m_in_trigger {this, "inTrigger", false, "Indicate if the tool is running on trigger"} |
|
protectedinherited |
◆ m_scoreName
Gaudi::Property<std::string> TauWPDecorator::m_scoreName {this, "ScoreName", ""} |
|
private |
◆ m_scoreNameTrans
Gaudi::Property<std::string> TauWPDecorator::m_scoreNameTrans {this, "NewScoreName", ""} |
|
private |
◆ m_tauRecToolsTag
Gaudi::Property<std::string> TauRecToolBase::m_tauRecToolsTag {this, "calibFolder", "tauRecTools/R22_preprod", "CVMFS path to the tau calibration folder"} |
|
protectedinherited |
◆ m_useAbsEta
Gaudi::Property<bool> TauWPDecorator::m_useAbsEta {this, "UseAbsEta", false} |
|
private |
◆ m_varHandleArraysDeclared
◆ m_vhka
◆ m_xMax
std::map<int, double> TauWPDecorator::m_xMax |
|
private |
Map of n-prong and the maximum value of x variables.
Definition at line 112 of file TauWPDecorator.h.
◆ m_xMin
std::map<int, double> TauWPDecorator::m_xMin |
|
private |
Map of n-prong and the minimum value of x variables.
Definition at line 110 of file TauWPDecorator.h.
◆ m_yMax
std::map<int, double> TauWPDecorator::m_yMax |
|
private |
Map of n-prong and the maximum value of y variables.
Definition at line 113 of file TauWPDecorator.h.
◆ m_yMin
std::map<int, double> TauWPDecorator::m_yMin |
|
private |
Map of n-prong and the minimum value of y variables.
Definition at line 111 of file TauWPDecorator.h.
The documentation for this class was generated from the following files:
Type
Property type enumeration.
std::map< int, double > m_yMax
Map of n-prong and the maximum value of y variables.
Gaudi::Property< std::vector< float > > m_EDMWPEffs2p
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T > &t)
Gaudi::Property< std::vector< std::string > > m_decorWPs
Gaudi::Property< std::vector< float > > m_EDMWPEffs3p
size_t nTracks(TauJetParameters::TauTrackFlag flag=TauJetParameters::TauTrackFlag::classifiedCharged) const
StoreGateSvc_t m_evtStore
Pointer to StoreGate (event store by default)
std::vector< SG::VarHandleKeyArray * > m_vhka
std::map< int, double > m_xMax
Map of n-prong and the maximum value of x variables.
double transformScore(double score, double cutLow, double effLow, double cutHigh, double effHigh) const
Obtain the flattened score.
#define ATH_MSG_VERBOSE(x)
Gaudi::Property< std::vector< float > > m_decorWPEffs0p
Gaudi::Property< std::vector< float > > m_EDMWPEffs0p
void setIsTau(TauJetParameters::IsTauFlag flag, bool value)
Set Flag for tau acceptance based on predefined arbitrary criteria.
std::shared_ptr< std::vector< m_pair_t > > m_hists1p
Efficiency and corresponding score distributions of 1-prong taus.
virtual void setOwner(IDataHandleHolder *o)=0
Gaudi::Property< std::vector< int > > m_EDMWPs
Gaudi::Property< std::vector< float > > m_decorWPEffs1p
virtual double pt() const
The transverse momentum ( ) of the particle.
ServiceHandle< StoreGateSvc > & evtStore()
The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.
SG::ReadDecorHandleKey< xAOD::EventInfo > m_aveIntPerXKey
void efficiency(std::vector< double > &bins, std::vector< double > &values, const std::vector< std::string > &files, const std::string &histname, const std::string &tplotname, const std::string &label="")
Gaudi::Property< std::string > m_scoreName
Handle class for reading a decoration on an object.
::StatusCode StatusCode
StatusCode definition for legacy code.
std::shared_ptr< std::vector< m_pair_t > > m_hists0p
Efficiency and corresponding score distributions of 0-prong taus.
std::vector< SG::Accessor< char > > m_charDecors
std::map< int, double > m_xMin
Map of n-prong and the minimum value of x variables.
StatusCode storeLimits(int nProng)
Obtain the limit of the dependent variables.
Gaudi::Property< std::vector< float > > m_EDMWPEffs1p
std::shared_ptr< std::vector< m_pair_t > > m_hists2p
Efficiency and corresponding score distributions of 2-prong taus.
const std::string & name(Level lvl)
Convenience function for translating message levels to strings.
Gaudi::Property< std::string > m_file3p
StoreGateSvc_t m_detStore
Pointer to StoreGate (detector store by default)
Gaudi::Property< std::string > m_file2p
Gaudi::Property< std::string > m_file0p
Gaudi::Property< std::string > m_file1p
std::map< int, double > m_yMin
Map of n-prong and the minimum value of y variables.
virtual void renounce()=0
std::conditional< std::is_base_of< SG::VarHandleKeyArray, T >::value, VarHandleKeyArrayType, type2 >::type type
Gaudi::Property< bool > m_useAbsEta
std::string PathResolverFindCalibFile(const std::string &logical_file_name)
#define ATH_MSG_WARNING(x)
Gaudi::Property< std::string > m_scoreNameTrans
IsTauFlag
Enum for IsTau flags.
SG::VarHandleKey & vhKey()
Return a non-const reference to the HandleKey.
Gaudi::Property< std::vector< float > > m_decorWPEffs3p
Gaudi::Property< bool > m_defineWPs
std::shared_ptr< std::vector< m_pair_t > > m_hists3p
Efficiency and corresponding score distributions of 3-prong taus.
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>
StatusCode retrieveHistos(int nProng)
Retrieve the histograms containing BDT/RNN score distributions as a function of dependent variables.
Gaudi::Property< std::vector< float > > m_decorWPEffs2p
std::pair< double, std::shared_ptr< TH2 > > m_pair_t