|
ATLAS Offline Software
|
#include <JSSTaggerUtils.h>
|
| JSSTaggerUtils (const std::string &name) |
| Constructor. More...
|
|
virtual StatusCode | initialize () override |
| Run once at the start of the job to setup everything. More...
|
|
virtual StatusCode | tag (const xAOD::Jet &jet) const override |
| IJetSelectorTool interface. More...
|
|
TH2D | MakeJetImage (TString TagImage, const xAOD::Jet *jet, std::vector< xAOD::JetConstituent > constituents) const override |
|
StatusCode | GetImageScore (const xAOD::JetContainer &jets) const override |
|
StatusCode | GetConstScore (const xAOD::JetContainer &jets) const override |
|
StatusCode | GetHLScore (const xAOD::JetContainer &jets) const override |
|
std::map< std::string, double > | GetJSSVars (const xAOD::Jet &jet) const override |
|
StatusCode | ReadScaler () override |
|
virtual void | print () const =0 |
| Print the state of the tool. More...
|
|
virtual StatusCode | decorate (const xAOD::JetContainer &jets) const override |
| Decorate jet collection with tagging info. 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 |
|
virtual StatusCode | modify (xAOD::JetContainer &jets) const override final |
| Concrete implementation of the function inherited from IJetModifier. More...
|
|
|
std::map< std::string, std::vector< double > > | m_scaler |
|
|
StatusCode | getConfigReader () |
| Get configReader StatusCode. More...
|
|
StatusCode | resetCuts (asg::AcceptData &acceptData) const |
| Reset cuts. More...
|
|
bool | passKinRange (const xAOD::Jet &jet) const |
| Check if jet passes kinematic constraints. More...
|
|
StatusCode | checkKinRange (const xAOD::Jet &jet, asg::AcceptData &acceptData) const |
| Check and record if jet passes kinematic constraints. More...
|
|
int | calculateJSSRatios (const xAOD::Jet &jet) const |
| Calculate JSS moment ratios in case they are not already saved TODO: Remove this once JSSMomentTools is modified to take const jets. More...
|
|
void | decorateJSSRatios (const xAOD::JetContainer &jets) const |
|
int | findPV () const |
| Find the PV (to be used for Ntrk) More...
|
|
int | GetUnGroomTracks (const xAOD::Jet &jet, int indexPV) const |
| Retrieve Ntrk variable from the ungroomed parent jet. More...
|
|
StatusCode | getWeight (const xAOD::Jet &jet, bool passSel, asg::AcceptData &acceptData) const |
| Get SF weight. More...
|
|
std::pair< double, double > | getSF (const xAOD::Jet &jet, const std::string &truthLabelStr) const |
| Get scale factor and efficiency. More...
|
|
std::string | getTruthLabelStr (const xAOD::Jet &jet, asg::AcceptData &acceptData) const |
|
void | printCuts () const |
| Print configured cuts. More...
|
|
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...
|
|
|
TEnv | m_configReader |
| TEnv instance to read config files. More...
|
|
asg::AcceptInfo | m_acceptInfo |
| Object that stores the results for a jet. More...
|
|
SG::WriteDecorHandleKey< xAOD::JetContainer > | m_decTaggedKey {this, "TaggedName", "Tagged", "SG key for Tagged"} |
| WriteDecorHandle keys for tagging bools. More...
|
|
SG::WriteDecorHandleKey< xAOD::JetContainer > | m_decValidPtRangeHighKey {this, "ValidPtRangeHighName", "ValidPtRangeHigh", "SG key for ValidPtRangeHigh"} |
|
SG::WriteDecorHandleKey< xAOD::JetContainer > | m_decValidPtRangeLowKey {this, "ValidPtRangeLowName", "ValidPtRangeLow", "SG key for ValidPtRangeLow"} |
|
SG::WriteDecorHandleKey< xAOD::JetContainer > | m_decValidEtaRangeKey {this, "ValidEtaRangeName", "ValidEtaRange", "SG key for ValidEtaRange"} |
|
SG::WriteDecorHandleKey< xAOD::JetContainer > | m_decValidKinRangeKey {this, "ValidKinRangeName", "ValidKinRange", "SG key for ValidKinRange"} |
|
SG::WriteDecorHandleKey< xAOD::JetContainer > | m_decValidJetContentKey {this, "ValidJetContentName", "ValidJetContent", "SG key for ValidJetContent"} |
|
SG::WriteDecorHandleKey< xAOD::JetContainer > | m_decValidEventContentKey {this, "ValidEventContentName", "ValidEventContent", "SG key for ValidEventContent"} |
|
SG::WriteDecorHandleKey< xAOD::JetContainer > | m_decPassMassKey {this, "PassMassName", "PassMass", "SG key for PassMass"} |
|
SG::WriteDecorHandleKey< xAOD::JetContainer > | m_decPassScoreKey {this, "PassScoreName", "PassScore", "SG key for PassScore"} |
|
SG::ReadDecorHandleKey< xAOD::JetContainer > | m_readD2Key {this, "D2Name", "D2", "SG key for D2"} |
|
SG::WriteDecorHandleKey< xAOD::JetContainer > | m_decNtrk500Key {this, "Ntrk500NameOut", "ParentJetNTrkPt500", "SG key for Ntrk500 from ungroomed jet"} |
|
SG::ReadDecorHandleKey< xAOD::JetContainer > | m_readNtrk500Key {this, "Ntrk500Name", "ParentJetNTrkPt500", "SG key for Ntrk500 from ungroomed jet"} |
|
const int | m_nWarnMax = 10 |
| Maximum number of warnings. More...
|
|
std::atomic< int > | m_nWarnVar |
| Warning counters. More...
|
|
bool | m_ptGeV = false |
| Flag to indicate units pT is defined in Set to false by default. More...
|
|
bool | m_useMassCut = false |
| Flag to indicate if mass window cut is used. More...
|
|
bool | m_useScoreCut = false |
| Flag to indicate if a discriminant score is used. More...
|
|
TAGCLASS | m_tagClass {Unknown} |
|
std::string | m_containerName |
| Configurable members. More...
|
|
std::string | m_weightConfigPath |
| Path to the SF configuration root file. More...
|
|
std::string | m_configFile |
| Configuration file name. More...
|
|
std::string | m_calibArea |
| Location where config files live on cvmfs. More...
|
|
std::string | m_kerasConfigFileName |
| Keras configurations for ML taggers. More...
|
|
std::string | m_kerasConfigFilePath |
|
std::string | m_kerasConfigOutputName |
|
std::string | m_kerasCalibArea |
|
std::string | m_tmvaConfigFileName |
| TMVA configurations for BDT taggers. More...
|
|
std::string | m_tmvaConfigFilePath |
|
std::string | m_tmvaCalibArea |
|
std::string | m_wkpt |
| Tagger information. More...
|
|
std::string | m_tagType |
|
float | m_jetPtMin {} |
| Kinematic bounds for the jet - the units are controlled by m_ptGeV. More...
|
|
float | m_jetPtMax {} |
|
float | m_jetEtaMax {} |
|
bool | m_useECFG = false |
| Flags controlling whether generalized ECF moments or L-series ratios are needed TODO: Implement the functionality controlled by these. More...
|
|
bool | m_useLSeries = false |
|
SG::WriteDecorHandleKey< xAOD::JetContainer > | m_decTau21WTAKey {this, "Tau21WTANameOut", "Tau21_wta", "SG key for Tau21_wta"} |
| WriteDecorHandle keys for JSS moments. More...
|
|
SG::WriteDecorHandleKey< xAOD::JetContainer > | m_decTau32WTAKey {this, "Tau32WTANameOut", "Tau32_wta", "SG key for Tau32_wta"} |
|
SG::WriteDecorHandleKey< xAOD::JetContainer > | m_decTau42WTAKey {this, "Tau42WTANameOut", "Tau42_wta", "SG key for Tau42_wta"} |
|
SG::WriteDecorHandleKey< xAOD::JetContainer > | m_decC2Key {this, "C2NameOut", "C2", "SG key for C2"} |
|
SG::WriteDecorHandleKey< xAOD::JetContainer > | m_decD2Key {this, "D2NameOut", "D2", "SG key for D2"} |
|
SG::WriteDecorHandleKey< xAOD::JetContainer > | m_decE3Key {this, "e3NameOut", "e3", "SG key for e3"} |
|
SG::WriteDecorHandleKey< xAOD::JetContainer > | m_decL2Key {this, "L2NameOut", "L2", "SG key for L2"} |
|
SG::WriteDecorHandleKey< xAOD::JetContainer > | m_decL3Key {this, "L3NameOut", "L3", "SG key for L3"} |
|
SG::ReadDecorHandleKey< xAOD::JetContainer > | m_readTau1WTAKey {this, "Tau1WTAName", "Tau1_wta", "SG key for Tau1_wta"} |
| ReadDecorHandle keys for JSS moments. More...
|
|
SG::ReadDecorHandleKey< xAOD::JetContainer > | m_readTau2WTAKey {this, "Tau2WTAName", "Tau2_wta", "SG key for Tau2_wta"} |
|
SG::ReadDecorHandleKey< xAOD::JetContainer > | m_readTau3WTAKey {this, "Tau3WTAName", "Tau3_wta", "SG key for Tau3_wta"} |
|
SG::ReadDecorHandleKey< xAOD::JetContainer > | m_readTau4WTAKey {this, "Tau4WTAName", "Tau4_wta", "SG key for Tau4_wta"} |
|
SG::ReadDecorHandleKey< xAOD::JetContainer > | m_readECF1Key {this, "ECF1Name", "ECF1", "SG key for ECF1"} |
|
SG::ReadDecorHandleKey< xAOD::JetContainer > | m_readECF2Key {this, "ECF2Name", "ECF2", "SG key for ECF2"} |
|
SG::ReadDecorHandleKey< xAOD::JetContainer > | m_readECF3Key {this, "ECF3Name", "ECF3", "SG key for ECF3"} |
|
SG::ReadDecorHandleKey< xAOD::JetContainer > | m_readSplit12Key {this, "Split12Name", "Split12", "SG key for Split12"} |
|
SG::ReadDecorHandleKey< xAOD::JetContainer > | m_readSplit23Key {this, "Split23Name", "Split23", "SG key for Split23"} |
|
SG::ReadDecorHandleKey< xAOD::JetContainer > | m_readQwKey {this, "QwName", "Qw", "SG key for Qw"} |
|
SG::ReadDecorHandleKey< xAOD::JetContainer > | m_readThrustMajKey {this, "ThrustMajName", "ThrustMaj", "SG key for ThrustMaj"} |
|
SG::ReadDecorHandleKey< xAOD::JetContainer > | m_readSphericityKey {this, "SphericityName", "Sphericity", "SG key for Sphericity"} |
|
SG::ReadDecorHandleKey< xAOD::JetContainer > | m_readECFG331Key {this, "ECFG331Name", "ECFG_3_3_1", "SG key for ECFG_3_3_1"} |
|
SG::ReadDecorHandleKey< xAOD::JetContainer > | m_readECFG311Key {this, "ECFG311Name", "ECFG_3_1_1", "SG key for ECFG_3_1_1"} |
|
SG::ReadDecorHandleKey< xAOD::JetContainer > | m_readECFG212Key {this, "ECFG212Name", "ECFG_2_1_2", "SG key for ECFG_2_1_2"} |
|
SG::ReadDecorHandleKey< xAOD::JetContainer > | m_readParentKey {this, "ParentName", "Parent", "SG key for Parent"} |
|
bool | m_suppressOutputDependence {} |
|
std::string | m_strMassCutLow |
| Strings for cut functions. More...
|
|
std::string | m_strMassCutHigh |
|
std::string | m_strScoreCut |
|
std::unique_ptr< TF1 > | m_funcMassCutLow |
| TF1 for cut functions. More...
|
|
std::unique_ptr< TF1 > | m_funcMassCutHigh |
|
std::unique_ptr< TF1 > | m_funcScoreCut |
|
std::string | m_decorationName |
| Decoration name. More...
|
|
bool | m_calcSF {} |
| Flag to calculate scale factor. More...
|
|
bool | m_isMC {} |
|
std::string | m_weightDecorationName |
| String for scale factor decoration names. More...
|
|
std::string | m_weightFileName |
|
std::string | m_weightHistogramName |
|
std::string | m_efficiencyHistogramName |
|
std::string | m_weightFlavors |
|
std::unique_ptr< TFile > | m_weightConfig |
| Histograms for scale factors. More...
|
|
std::map< std::string, std::unique_ptr< TH2D > > | m_weightHistograms |
|
std::map< std::string, std::unique_ptr< TH2D > > | m_efficiencyHistograms |
|
bool | m_truthLabelUseTRUTH3 {} |
| Truth label options. More...
|
|
std::string | m_truthParticleContainerName |
|
std::string | m_truthBosonContainerName |
|
std::string | m_truthTopQuarkContainerName |
|
std::string | m_truthLabelName |
|
SG::ReadDecorHandleKey< xAOD::JetContainer > | m_readTruthLabelKey {this, "truthLabelName", "truthLabel", "SG key for truthLabel"} |
| Truth label ReadDecorHandle key. More...
|
|
SG::WriteDecorHandleKey< xAOD::JetContainer > | m_decCutMLowKey {this, "CutMLowName", "Cut_mlow", "SG key for Cut_mlow"} |
| WriteDecorHandle keys for cut values. More...
|
|
SG::WriteDecorHandleKey< xAOD::JetContainer > | m_decCutMHighKey {this, "CutMHighName", "Cut_mhigh", "SG key for Cut_mhigh"} |
|
SG::WriteDecorHandleKey< xAOD::JetContainer > | m_decScoreCutKey {this, "CutScoreName", "Cut_score", "SG key for Cut_score"} |
|
SG::WriteDecorHandleKey< xAOD::JetContainer > | m_decScoreValueKey {this, "ScoreName", "Score", "SG key for Score"} |
|
SG::WriteDecorHandleKey< xAOD::JetContainer > | m_decWeightKey {this, "weightName", "weight", "SG key for weight"} |
| WriteDecorHandle keys for SF. More...
|
|
SG::WriteDecorHandleKey< xAOD::JetContainer > | m_decEfficiencyKey {this, "efficiencyName", "efficiency", "SG key for efficiency"} |
|
SG::WriteDecorHandleKey< xAOD::JetContainer > | m_decEffSFKey {this, "effSFName", "effSF", "SG key for effSF"} |
|
SG::WriteDecorHandleKey< xAOD::JetContainer > | m_decSigeffSFKey {this, "sigeffSFName", "sigeffSF", "SG key for effSF"} |
|
|
int | m_nbins_eta {} |
|
int | m_nbins_phi {} |
|
double | m_min_eta {} |
|
double | m_max_eta {} |
|
double | m_min_phi {} |
|
double | m_max_phi {} |
|
int | m_ncolors {} |
|
bool | m_dorscaling {} |
|
double | m_rscaling_p0 {} |
|
double | m_rscaling_p1 {} |
|
bool | m_UseConstTagger {} |
|
ToolHandle< AthONNX::IJSSMLTool > | m_MLBosonTagger {"", this} |
|
bool | m_UseHLTagger {} |
|
ToolHandle< AthONNX::IJSSMLTool > | m_MLBosonTagger_HL {"", this} |
|
SG::WriteDecorHandleKey< xAOD::JetContainer > | m_decNConstituentsKey {this, "nConstituentsName", "nConstituents", "SG key for constituents multiplicity"} |
| WriteDecorHandle keys. More...
|
|
SG::WriteDecorHandleKey< xAOD::JetContainer > | m_decNTopoTowersKey {this, "nTopoTowersName", "nTopoTowers", "SG key for towers multiplicity"} |
|
SG::WriteDecorHandleKey< xAOD::JetContainer > | m_decConstScoreKey {this, "ConstScoreNameDec", "ConstScore", "SG key for ConstScore"} |
|
SG::ReadDecorHandleKey< xAOD::JetContainer > | m_readConstScoreKey {this, "ConstScoreNameRead", "ConstScore", "SG key for ConstScore"} |
|
SG::WriteDecorHandleKey< xAOD::JetContainer > | m_decHLScoreKey {this, "HLScoreName", "HLScore", "SG key for HLScore"} |
|
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 |
|
Definition at line 14 of file JSSTaggerUtils.h.
◆ StoreGateSvc_t
◆ TAGCLASS
TAGTYPE enum.
Enumerator |
---|
Unknown | |
WBoson | |
ZBoson | |
TopQuark | |
Definition at line 100 of file JSSTaggerBase.h.
◆ JSSTaggerUtils()
JSSTaggerUtils::JSSTaggerUtils |
( |
const std::string & |
name | ) |
|
◆ calculateJSSRatios()
Calculate JSS moment ratios in case they are not already saved TODO: Remove this once JSSMomentTools is modified to take const jets.
Calculate JSS moment ratios in case they are not already saved These are calculated by hand here because JetSubStructureMomentTools does not operate on const jets.
This should be changed in the future
Create write decor handles
Create read decor handles
WTA N-subjettiness ratios
ECF ratios
Definition at line 399 of file JSSTaggerBase.cxx.
429 float tau21_wta = -999.0;
430 float tau32_wta = -999.0;
431 float tau42_wta = -999.0;
433 float tau1_wta = readTau1WTA(
jet);
434 float tau2_wta = readTau2WTA(
jet);
435 float tau3_wta = readTau3WTA(
jet);
436 float tau4_wta = -999.0;
437 if(readTau4WTA.isAvailable()){
438 tau4_wta = readTau4WTA(
jet);
441 if ( tau1_wta > 1
e-8 ) {
442 tau21_wta = tau2_wta / tau1_wta;
446 if ( tau2_wta > 1
e-8 ) {
447 tau32_wta = tau3_wta / tau2_wta;
448 if(readTau4WTA.isAvailable()){
449 tau42_wta = tau4_wta / tau2_wta;
454 decTau21WTA(
jet) = tau21_wta;
455 decTau32WTA(
jet) = tau32_wta;
456 decTau42WTA(
jet) = tau42_wta;
463 float ECF1 = readECF1(
jet);
464 float ECF2 = readECF2(
jet);
465 float ECF3 = readECF3(
jet);
468 C2 = ECF3 * ECF1 /
std::pow( ECF2, 2.0 );
484 if(!accL2.isAvailable(
jet)){
485 if(readECFG331.isAvailable() && readECFG212.isAvailable()){
486 if(readECFG212(
jet) > 1
e-8){
487 L2 = readECFG331(
jet) /
pow(readECFG212(
jet), (3.0/2.0));
495 if(!accL3.isAvailable(
jet)){
496 if(readECFG331.isAvailable() && readECFG311.isAvailable()){
497 if(readECFG331(
jet) > 1
e-8){
498 L3 = readECFG311(
jet) /
pow(readECFG331(
jet), (1.0/3.0));
◆ checkKinRange()
Check and record if jet passes kinematic constraints.
Check each kinematic constraint Print warnings using counters
Create write decor handles
Decorate kinematic pass information
Definition at line 353 of file JSSTaggerBase.cxx.
387 decValidPtRangeHigh(
jet) = acceptData.
getCutResult(
"ValidPtRangeHigh" );
388 decValidPtRangeLow(
jet) = acceptData.
getCutResult(
"ValidPtRangeLow" );
392 return StatusCode::SUCCESS;
◆ declareGaudiProperty() [1/4]
specialization for handling Gaudi::Property<SG::VarHandleKeyArray>
Definition at line 170 of file AthCommonDataStore.h.
175 hndl.documentation());
◆ declareGaudiProperty() [2/4]
specialization for handling Gaudi::Property<SG::VarHandleKey>
Definition at line 156 of file AthCommonDataStore.h.
161 hndl.documentation());
◆ declareGaudiProperty() [3/4]
specialization for handling Gaudi::Property<SG::VarHandleBase>
Definition at line 184 of file AthCommonDataStore.h.
189 hndl.documentation());
◆ declareGaudiProperty() [4/4]
◆ declareProperty() [1/6]
Declare a new Gaudi property.
- Parameters
-
name | Name of the property. |
hndl | Object holding the property value. |
doc | Documentation string for the property. |
This is the version for types that derive from SG::VarHandleBase
. The property value object is put on the input and output lists as appropriate; then we forward to the base class.
Definition at line 245 of file AthCommonDataStore.h.
250 this->declare(hndl.
vhKey());
251 hndl.
vhKey().setOwner(
this);
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]
◆ decorate()
◆ decorateJSSRatios()
Create write decor handles
Create read decor handles
WTA N-subjettiness ratios
ECF ratios
Definition at line 512 of file JSSTaggerBase.cxx.
516 bool decor_available = decTau21WTA.isAvailable();
517 if(decor_available)
return;
544 float tau21_wta = -999.0;
545 float tau32_wta = -999.0;
546 float tau42_wta = -999.0;
548 float tau1_wta = readTau1WTA(*
jet);
549 float tau2_wta = readTau2WTA(*
jet);
550 float tau3_wta = readTau3WTA(*
jet);
551 float tau4_wta = -999.0;
552 if(readTau4WTA.isAvailable()){
553 tau4_wta = readTau4WTA(*
jet);
556 if ( tau1_wta > 1
e-8 ) {
557 tau21_wta = tau2_wta / tau1_wta;
560 if ( tau2_wta > 1
e-8 ) {
561 tau32_wta = tau3_wta / tau2_wta;
562 if(readTau4WTA.isAvailable()){
563 tau42_wta = tau4_wta / tau2_wta;
567 decTau21WTA(*
jet) = tau21_wta;
568 decTau32WTA(*
jet) = tau32_wta;
569 decTau42WTA(*
jet) = tau42_wta;
576 float ECF1 = readECF1(*
jet);
577 float ECF2 = readECF2(*
jet);
578 float ECF3 = readECF3(*
jet);
581 C2 = ECF3 * ECF1 /
std::pow( ECF2, 2.0 );
596 if(!accL2.isAvailable(*
jet)){
597 if(readECFG331.isAvailable() && readECFG212.isAvailable()){
598 if(readECFG212(*
jet) > 1
e-8){
599 L2 = readECFG331(*
jet) /
pow(readECFG212(*
jet), (3.0/2.0));
606 if(!accL3.isAvailable(*
jet)){
607 if(readECFG331.isAvailable() && readECFG311.isAvailable()){
608 if(readECFG331(*
jet) > 1
e-8){
609 L3 = readECFG311(*
jet) /
pow(readECFG331(*
jet), (1.0/3.0));
◆ detStore()
◆ evtStore() [1/2]
◆ evtStore() [2/2]
◆ extraDeps_update_handler()
Add StoreName to extra input/output deps as needed.
use the logic of the VarHandleKey to parse the DataObjID keys supplied via the ExtraInputs and ExtraOuputs Properties to add the StoreName if it's not explicitly given
◆ findPV()
int JSSTaggerBase::findPV |
( |
| ) |
const |
|
protectedinherited |
Find the PV (to be used for Ntrk)
Definition at line 624 of file JSSTaggerBase.cxx.
630 ATH_MSG_WARNING(
"Unable to retrieve primary vertex container PrimaryVertices" );
634 for (
const auto *vx : *vxCont ) {
636 indexPV = vtx_counter;
◆ getConfigReader()
StatusCode JSSTaggerBase::getConfigReader |
( |
| ) |
|
|
protectedinherited |
Get configReader StatusCode.
Get configReader TEnv.
Check for the existence of the configuration file
https://root.cern.ch/root/roottalk/roottalk02/5332.html
Definition at line 283 of file JSSTaggerBase.cxx.
288 std::string configPath;
293 else if (
m_calibArea.find(
"eos") != std::string::npos) {
302 int fSuccess = gSystem->GetPathInfo(configPath.c_str(), fStats);
305 return StatusCode::FAILURE;
308 ATH_MSG_DEBUG(
"Recommendations file was found : " << configPath );
311 if (
m_configReader.ReadFile( configPath.c_str(), EEnvLevel(0) ) ) {
312 ATH_MSG_ERROR(
"Error while reading config file : "<< configPath );
313 return StatusCode::FAILURE;
316 return StatusCode::SUCCESS;
◆ GetConstScore()
Implements IJSSTaggerUtils.
Definition at line 215 of file JSSTaggerUtils.cxx.
227 std::vector<xAOD::JetConstituent> constituents =
jet -> getConstituents().asSTLVector();
230 int MaxConstituents (100);
231 std::vector<xAOD::JetConstituent> constituentsForModel;
233 if( constituents.size() > 100 )
234 constituentsForModel = std::vector<xAOD::JetConstituent> (constituents.begin(), constituents.begin() + MaxConstituents);
236 constituentsForModel = constituents;
239 std::vector<const xAOD::CaloCluster*>
towers;
241 if (towersAcc.isAvailable(*
jet)){
243 std::vector<ElementLink<DataVector<xAOD::IParticle>>> towerLinks = towersAcc(*
jet);
244 for (
auto link_itr : towerLinks){
245 if (!link_itr.isValid())
continue;
253 for(
auto cnst : constituents){
254 m.push_back( cnst ->
m() );
255 pT.push_back( cnst ->
pt() );
256 eta.push_back( cnst ->
eta() );
257 phi.push_back( cnst ->
phi() );
259 std::vector<std::vector<float>> constituents_packed = {
m,
pT,
eta,
phi};
261 m.clear();
pT.clear();
eta.clear();
phi.clear();
263 m.push_back( cnst ->
m() );
264 pT.push_back( cnst ->
pt() );
265 eta.push_back( cnst ->
eta() );
266 phi.push_back( cnst ->
phi() );
268 std::vector<std::vector<float>> towers_packed = {
m,
pT,
eta,
phi};
271 std::vector<std::vector<float>> inputs_packed = {
272 constituents_packed.at(0), constituents_packed.at(1), constituents_packed.at(2), constituents_packed.at(3),
273 towers_packed.at(0), towers_packed.at(1), towers_packed.at(2), towers_packed.at(3),
278 if( constituents.size()>1 )
285 decNConstituents(*
jet) = constituents.size();
290 return StatusCode::SUCCESS;
◆ GetHLScore()
◆ GetImageScore()
Implements IJSSTaggerUtils.
Definition at line 149 of file JSSTaggerUtils.cxx.
162 std::vector<xAOD::JetConstituent> constituents =
jet -> getConstituents().asSTLVector();
165 int MaxConstituents (100);
166 std::vector<xAOD::JetConstituent> constituentsForModel;
168 if( constituents.size() > 100 )
169 constituentsForModel = std::vector<xAOD::JetConstituent> (constituents.begin(), constituents.begin() + MaxConstituents);
171 constituentsForModel = constituents;
174 std::vector<xAOD::JetConstituent> csts_charged = constituentsForModel;
175 csts_charged.erase( std::remove_if( csts_charged.begin(), csts_charged.end(),
177 const xAOD::FlowElement* ufo = dynamic_cast<const xAOD::FlowElement*>(constituent.rawConstituent());
178 return ufo -> signalType() != xAOD::FlowElement::SignalType::Charged;
179 }), csts_charged.end()) ;
182 std::vector<xAOD::JetConstituent> csts_neutral = constituentsForModel;
183 csts_neutral.erase( std::remove_if( csts_neutral.begin(), csts_neutral.end(),
185 const xAOD::FlowElement* ufo = dynamic_cast<const xAOD::FlowElement*>(constituent.rawConstituent());
186 return ufo -> signalType() != xAOD::FlowElement::SignalType::Neutral;
187 }), csts_neutral.end()) ;
190 std::vector<xAOD::JetConstituent> csts_combined = constituentsForModel;
191 csts_combined.erase( std::remove_if( csts_combined.begin(), csts_combined.end(),
193 const xAOD::FlowElement* ufo = dynamic_cast<const xAOD::FlowElement*>(constituent.rawConstituent());
194 return ufo -> signalType() != xAOD::FlowElement::SignalType::Combined;
195 }), csts_combined.end()) ;
202 std::vector<TH2D> Images = {ImageCharged, ImageNeutral, ImageCombined};
211 return StatusCode::SUCCESS;
◆ GetJSSVars()
std::map< std::string, double > JSSTaggerUtils::GetJSSVars |
( |
const xAOD::Jet & |
jet | ) |
const |
|
overridevirtual |
Implements IJSSTaggerUtils.
Definition at line 321 of file JSSTaggerUtils.cxx.
323 std::map<std::string, double> JSSVars;
326 int nUngrTracks (-1);
328 if(readNtrk500.isAvailable()){
329 nUngrTracks = readNtrk500(
jet);
332 int pv_location =
findPV();
333 if(pv_location != -1)
338 JSSVars[
"pT"] =
jet.pt();
339 JSSVars[
"nTracks"] = nUngrTracks;
342 JSSVars[
"CNN"] = readConstScore(
jet);
370 JSSVars[
"Split12"] = readSplit12(
jet);
371 JSSVars[
"Split23"] = readSplit23(
jet);
374 JSSVars[
"D2"] = readD2(
jet);
377 JSSVars[
"Tau1_wta"] = readTau1_wta(
jet);
378 JSSVars[
"Tau2_wta"] = readTau2_wta(
jet);
379 JSSVars[
"Tau3_wta"] = readTau3_wta(
jet);
382 JSSVars[
"ECF1"] = readECF1(
jet);
383 JSSVars[
"ECF2"] = readECF2(
jet);
384 JSSVars[
"ECF3"] = readECF3(
jet);
387 JSSVars[
"Qw"] = readQw(
jet);
390 JSSVars[
"FoxWolfram0"] = FoxWolfram0Acc.withDefault(
jet, -99.);
391 JSSVars[
"FoxWolfram2"] = FoxWolfram2Acc.withDefault(
jet, -99.);
392 JSSVars[
"PlanarFlow"] = PlanarFlowAcc.withDefault(
jet, -99.);
393 JSSVars[
"Angularity"] = AngularityAcc.withDefault(
jet, -99.);
394 JSSVars[
"Aplanarity"] = AplanarityAcc.withDefault(
jet, -99.);
395 JSSVars[
"ZCut12"] = ZCut12Acc.withDefault(
jet, -99.);
396 JSSVars[
"KtDR"] = KtDRAcc.withDefault(
jet, -99.);
◆ 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.
◆ getSF()
std::pair< double, double > JSSTaggerBase::getSF |
( |
const xAOD::Jet & |
jet, |
|
|
const std::string & |
truthLabelStr |
|
) |
| const |
|
protectedinherited |
Get scale factor and efficiency.
To apply W-tagging efficiency SF to Z-tagger, jet mass is shifted by 10.803 GeV
Definition at line 760 of file JSSTaggerBase.cxx.
769 const double WtoZmassShift = 10803;
773 if ( logmOverPt > 0 ) logmOverPt = 0;
790 std::string signal_truthLabel=
"";
792 signal_truthLabel=
"t_qqb";
794 signal_truthLabel=
"V_qq";
796 signal_truthLabel=
"t";
803 return std::make_pair( 1.0,
eff );
807 ATH_MSG_DEBUG(
"(pt, m/pt) (" <<
jet.pt()/1.e3 <<
", " <<
jet.m()/
jet.pt() <<
") is out of range for SF calculation. Returning 1.0" );
808 return std::make_pair( 1.0, 1.0 );
811 return std::make_pair(
SF,
eff );
◆ getTruthLabelStr()
Truth label string
Truth label value
Contained top tagger
Contained top
QCD
TCC W/Z 2-var tagger
Top
W/Z
QCD
Pass mass and D2
Fail mass, pass D2
Pass mass, fail D2
Fail mass and D2
W/Z tagger
Top
W/Z
QCD
Top
QCD
Definition at line 818 of file JSSTaggerBase.cxx.
820 std::string truthLabelStr;
831 truthLabelStr =
"t_qqb";
844 truthLabelStr =
"t_";
848 truthLabelStr =
"V_qq_";
852 truthLabelStr =
"q_";
857 truthLabelStr +=
"passMpassD2";
861 truthLabelStr +=
"failMpassD2";
865 truthLabelStr +=
"passMfailD2";
869 truthLabelStr +=
"failMfailD2";
883 truthLabelStr =
"V_qq";
904 return truthLabelStr;
◆ GetUnGroomTracks()
int JSSTaggerBase::GetUnGroomTracks |
( |
const xAOD::Jet & |
jet, |
|
|
int |
indexPV |
|
) |
| const |
|
protectedinherited |
Retrieve Ntrk variable from the ungroomed parent jet.
Retrieve the Ntrk variable from the ungroomed parent jet.
Definition at line 648 of file JSSTaggerBase.cxx.
652 const xAOD::Jet * ungroomedJet =
nullptr;
653 if ( readParent.isAvailable() ) {
655 if ( linkToUngroomed.
isValid() ) {
656 ungroomedJet = *linkToUngroomed;
660 if ( acc_Ntrk.isAvailable(*ungroomedJet) ) {
662 const std::vector<int>& NTrkPt500 = acc_Ntrk(*ungroomedJet);
664 int jet_ntrk = NTrkPt500.at(indexPV);
669 ATH_MSG_ERROR(
"WARNING: Unable to retrieve Ntrk of the ungroomed parent jet. Please make sure this variable is in your derivations!!!");
673 ATH_MSG_ERROR(
"WARNING: Unable to retrieve the parent ungroomed jet. Please make sure this variable is in your derivations!!!");
677 ATH_MSG_ERROR(
"WARNING: Unable to retrieve the link to the parent ungroomed jet. Please make sure this variable is in your derivations!!!");
◆ getWeight()
Get SF weight.
Inefficiency SF is directly used
Efficiency SF
Calculate inefficiency SF
If inefficiency SF is not available, SF is always 1.0
Create write decor handles
Decorate values
Definition at line 686 of file JSSTaggerBase.cxx.
687 if ( !
m_calcSF )
return StatusCode::SUCCESS;
691 float sigeffSF = 1.0;
701 sigeffSF =
getSF(
jet,
"t_qqb").first;
703 sigeffSF =
getSF(
jet,
"V_qq").first;
752 decEffSF(
jet) = effSF;
753 decSigeffSF(
jet) = sigeffSF;
755 return StatusCode::SUCCESS;
◆ initialize()
StatusCode JSSTaggerUtils::initialize |
( |
| ) |
|
|
overridevirtual |
Run once at the start of the job to setup everything.
Get configReader
Get the decoration name
Initialize decorators
Reimplemented from JSSTaggerBase.
Definition at line 30 of file JSSTaggerUtils.cxx.
48 std::string ConstTaggerFileName =
m_configReader.GetValue(
"ConstTaggerFileName",
"aaa");
49 std::string HLTaggerFileName =
m_configReader.GetValue(
"HLTaggerFileName",
"aaa");
53 std::string ModelPath =
"";
57 else if (
m_calibArea.find(
"eos") != std::string::npos) {
58 ModelPath = (ConstTaggerFileName);
63 ATH_MSG_INFO(
"JSSTaggerUtils::MLBosonTagger()" <<
" + ModelPath " << ModelPath );
69 ATH_MSG_INFO(
"JSSTaggerUtils::MLBosonTagger() read value from config" );
93 std::string ModelPath =
"";
97 else if (
m_calibArea.find(
"eos") != std::string::npos) {
98 ModelPath = (HLTaggerFileName);
103 ATH_MSG_INFO(
"JSSTaggerUtils::MLBosonTagger()" <<
" + ModelPath " << ModelPath );
112 ATH_MSG_INFO(
"JSSTaggerUtils::MLBosonTagger() read value from config" );
120 ATH_MSG_INFO(
"Decorators that will be attached to jet :" );
137 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.
◆ MakeJetImage()
Implements IJSSTaggerUtils.
Definition at line 446 of file JSSTaggerUtils.cxx.
449 double eta (-99.),
phi (-99.),
pT (-99.),
z (-99.);
450 int BinEta (-99), BinPhi (-99);
453 return sum + cst.pt(); };
454 double SumPT =
std::accumulate( constituents.begin(), constituents.end(), 0., cst_pT) ;
456 auto Image = std::make_unique<TH2D>(
"Image_" + TagImage,
"Image_" + TagImage,
459 for(
auto& cst : constituents ){
477 double current_z = Image -> GetBinContent(
x,
y );
480 Image -> SetBinContent(
x,
y, current_z +=
z );
◆ modify()
|
inlinefinaloverridevirtualinherited |
◆ 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.
◆ passKinRange()
◆ print() [1/2]
void asg::AsgTool::print |
( |
| ) |
const |
|
virtualinherited |
Print the state of the tool.
Implements asg::IAsgTool.
Reimplemented in JetRecTool, JetFinder, JetModifiedMassDrop, JetFromPseudojet, JetReclusterer, JetReclusteringTool, JetTruthLabelingTool, JetPileupLabelingTool, HI::HIPileupTool, JetDumper, JetBottomUpSoftDrop, JetRecursiveSoftDrop, JetSoftDrop, JetConstituentsRetriever, JetSubStructureMomentToolsBase, JetSplitter, JetToolRunner, JetPruner, JetPseudojetRetriever, JetTrimmer, AsgHelloTool, and KtDeltaRTool.
Definition at line 131 of file AsgTool.cxx.
◆ print() [2/2]
virtual void asg::IAsgTool::print |
( |
| ) |
const |
|
pure virtualinherited |
Print the state of the tool.
Implemented in JetRecTool, JetFinder, JetModifiedMassDrop, JetFromPseudojet, JetReclusterer, JetReclusteringTool, JetTruthLabelingTool, JetPileupLabelingTool, HI::HIPileupTool, asg::AsgTool, JetDumper, JetBottomUpSoftDrop, JetRecursiveSoftDrop, JetSoftDrop, JetConstituentsRetriever, JetSubStructureMomentToolsBase, JetSplitter, JetToolRunner, JetPruner, JetPseudojetRetriever, JetTrimmer, AsgHelloTool, and KtDeltaRTool.
◆ printCuts()
void JSSTaggerBase::printCuts |
( |
| ) |
const |
|
protectedinherited |
Print configured cuts.
Definition at line 907 of file JSSTaggerBase.cxx.
908 ATH_MSG_INFO(
"After tagging, you will have access to the following cuts as an asg::AcceptData : (<NCut>) <cut> : <description>)" );
910 for (
int iCut=0; iCut < nCuts; iCut++ ) {
911 std::string cut_string =
"";
◆ ReadScaler()
StatusCode JSSTaggerUtils::ReadScaler |
( |
| ) |
|
|
overridevirtual |
Implements IJSSTaggerUtils.
Definition at line 402 of file JSSTaggerUtils.cxx.
405 std::vector<std::string> vars_list = {
406 "pT",
"CNN",
"D2",
"nTracks",
"ZCut12",
407 "Tau1_wta",
"Tau2_wta",
"Tau3_wta",
408 "KtDR",
"Split12",
"Split23",
409 "ECF1",
"ECF2",
"ECF3",
410 "Angularity",
"FoxWolfram0",
"FoxWolfram2",
411 "Aplanarity",
"PlanarFlow",
"Qw",
415 for(
const std::string &
var : vars_list){
416 std::string s_mean =
var +
"_mean";
417 std::string s_std =
var +
"_std";
421 if(
mean==-99. && std==-99.){
422 ATH_MSG_ERROR(
"ERROR: one of the parameter for " <<
var <<
" is missing, please, double check the config!!!");
423 return StatusCode::FAILURE;
425 else if(
mean==-99. || std==-99.){
426 ATH_MSG_ERROR(
"ERROR: parameters for " <<
var <<
" are missing, please, double check the config!!!");
427 return StatusCode::FAILURE;
442 return StatusCode::SUCCESS;
◆ renounce()
◆ renounceArray()
◆ resetCuts()
Reset cuts.
Reset the AcceptData cut results to false
Initialize common cuts to true by default
Definition at line 321 of file JSSTaggerBase.cxx.
334 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.
◆ tag()
◆ updateVHKA()
◆ m_acceptInfo
Object that stores the results for a jet.
Definition at line 65 of file JSSTaggerBase.h.
◆ m_calcSF
bool JSSTaggerBase::m_calcSF {} |
|
protectedinherited |
◆ m_calibArea
std::string JSSTaggerBase::m_calibArea |
|
protectedinherited |
Location where config files live on cvmfs.
Definition at line 115 of file JSSTaggerBase.h.
◆ m_configFile
std::string JSSTaggerBase::m_configFile |
|
protectedinherited |
◆ m_configReader
TEnv JSSTaggerBase::m_configReader |
|
protectedinherited |
◆ m_containerName
std::string JSSTaggerBase::m_containerName |
|
protectedinherited |
◆ m_decC2Key
◆ m_decConstScoreKey
◆ m_decCutMHighKey
◆ m_decCutMLowKey
◆ m_decD2Key
◆ m_decE3Key
◆ m_decEfficiencyKey
◆ m_decEffSFKey
◆ m_decHLScoreKey
◆ m_decL2Key
◆ m_decL3Key
◆ m_decNConstituentsKey
◆ m_decNTopoTowersKey
◆ m_decNtrk500Key
◆ m_decorationName
std::string JSSTaggerBase::m_decorationName |
|
protectedinherited |
◆ m_decPassMassKey
◆ m_decPassScoreKey
◆ m_decScoreCutKey
◆ m_decScoreValueKey
◆ m_decSigeffSFKey
◆ m_decTaggedKey
WriteDecorHandle keys for tagging bools.
Definition at line 68 of file JSSTaggerBase.h.
◆ m_decTau21WTAKey
◆ m_decTau32WTAKey
◆ m_decTau42WTAKey
◆ m_decValidEtaRangeKey
◆ m_decValidEventContentKey
◆ m_decValidJetContentKey
◆ m_decValidKinRangeKey
◆ m_decValidPtRangeHighKey
◆ m_decValidPtRangeLowKey
◆ m_decWeightKey
◆ m_detStore
◆ m_dorscaling
bool JSSTaggerUtils::m_dorscaling {} |
|
private |
◆ m_efficiencyHistogramName
std::string JSSTaggerBase::m_efficiencyHistogramName |
|
protectedinherited |
◆ m_efficiencyHistograms
std::map<std::string, std::unique_ptr<TH2D> > JSSTaggerBase::m_efficiencyHistograms |
|
protectedinherited |
◆ m_evtStore
◆ m_funcMassCutHigh
std::unique_ptr<TF1> JSSTaggerBase::m_funcMassCutHigh |
|
protectedinherited |
◆ m_funcMassCutLow
std::unique_ptr<TF1> JSSTaggerBase::m_funcMassCutLow |
|
protectedinherited |
◆ m_funcScoreCut
std::unique_ptr<TF1> JSSTaggerBase::m_funcScoreCut |
|
protectedinherited |
◆ m_isMC
bool JSSTaggerBase::m_isMC {} |
|
protectedinherited |
◆ m_jetEtaMax
float JSSTaggerBase::m_jetEtaMax {} |
|
protectedinherited |
◆ m_jetPtMax
float JSSTaggerBase::m_jetPtMax {} |
|
protectedinherited |
◆ m_jetPtMin
float JSSTaggerBase::m_jetPtMin {} |
|
protectedinherited |
Kinematic bounds for the jet - the units are controlled by m_ptGeV.
Definition at line 133 of file JSSTaggerBase.h.
◆ m_kerasCalibArea
std::string JSSTaggerBase::m_kerasCalibArea |
|
protectedinherited |
◆ m_kerasConfigFileName
std::string JSSTaggerBase::m_kerasConfigFileName |
|
protectedinherited |
◆ m_kerasConfigFilePath
std::string JSSTaggerBase::m_kerasConfigFilePath |
|
protectedinherited |
◆ m_kerasConfigOutputName
std::string JSSTaggerBase::m_kerasConfigOutputName |
|
protectedinherited |
◆ m_max_eta
double JSSTaggerUtils::m_max_eta {} |
|
private |
◆ m_max_phi
double JSSTaggerUtils::m_max_phi {} |
|
private |
◆ m_min_eta
double JSSTaggerUtils::m_min_eta {} |
|
private |
◆ m_min_phi
double JSSTaggerUtils::m_min_phi {} |
|
private |
◆ m_MLBosonTagger
◆ m_MLBosonTagger_HL
◆ m_nbins_eta
int JSSTaggerUtils::m_nbins_eta {} |
|
private |
◆ m_nbins_phi
int JSSTaggerUtils::m_nbins_phi {} |
|
private |
◆ m_ncolors
int JSSTaggerUtils::m_ncolors {} |
|
private |
◆ m_nWarnMax
const int JSSTaggerBase::m_nWarnMax = 10 |
|
protectedinherited |
◆ m_nWarnVar
std::atomic<int> JSSTaggerBase::m_nWarnVar |
|
mutableprotectedinherited |
◆ m_ptGeV
bool JSSTaggerBase::m_ptGeV = false |
|
protectedinherited |
Flag to indicate units pT is defined in Set to false by default.
Definition at line 91 of file JSSTaggerBase.h.
◆ m_readConstScoreKey
◆ m_readD2Key
◆ m_readECF1Key
◆ m_readECF2Key
◆ m_readECF3Key
◆ m_readECFG212Key
◆ m_readECFG311Key
◆ m_readECFG331Key
◆ m_readNtrk500Key
◆ m_readParentKey
◆ m_readQwKey
◆ m_readSphericityKey
◆ m_readSplit12Key
◆ m_readSplit23Key
◆ m_readTau1WTAKey
◆ m_readTau2WTAKey
◆ m_readTau3WTAKey
◆ m_readTau4WTAKey
◆ m_readThrustMajKey
◆ m_readTruthLabelKey
◆ m_rscaling_p0
double JSSTaggerUtils::m_rscaling_p0 {} |
|
private |
◆ m_rscaling_p1
double JSSTaggerUtils::m_rscaling_p1 {} |
|
private |
◆ m_scaler
std::map<std::string, std::vector<double> > JSSTaggerUtils::m_scaler |
◆ m_strMassCutHigh
std::string JSSTaggerBase::m_strMassCutHigh |
|
protectedinherited |
◆ m_strMassCutLow
std::string JSSTaggerBase::m_strMassCutLow |
|
protectedinherited |
◆ m_strScoreCut
std::string JSSTaggerBase::m_strScoreCut |
|
protectedinherited |
◆ m_suppressOutputDependence
bool JSSTaggerBase::m_suppressOutputDependence {} |
|
protectedinherited |
◆ m_tagClass
◆ m_tagType
std::string JSSTaggerBase::m_tagType |
|
protectedinherited |
◆ m_tmvaCalibArea
std::string JSSTaggerBase::m_tmvaCalibArea |
|
protectedinherited |
◆ m_tmvaConfigFileName
std::string JSSTaggerBase::m_tmvaConfigFileName |
|
protectedinherited |
◆ m_tmvaConfigFilePath
std::string JSSTaggerBase::m_tmvaConfigFilePath |
|
protectedinherited |
◆ m_truthBosonContainerName
std::string JSSTaggerBase::m_truthBosonContainerName |
|
protectedinherited |
◆ m_truthLabelName
std::string JSSTaggerBase::m_truthLabelName |
|
protectedinherited |
◆ m_truthLabelUseTRUTH3
bool JSSTaggerBase::m_truthLabelUseTRUTH3 {} |
|
protectedinherited |
◆ m_truthParticleContainerName
std::string JSSTaggerBase::m_truthParticleContainerName |
|
protectedinherited |
◆ m_truthTopQuarkContainerName
std::string JSSTaggerBase::m_truthTopQuarkContainerName |
|
protectedinherited |
◆ m_UseConstTagger
bool JSSTaggerUtils::m_UseConstTagger {} |
|
private |
◆ m_useECFG
bool JSSTaggerBase::m_useECFG = false |
|
protectedinherited |
Flags controlling whether generalized ECF moments or L-series ratios are needed TODO: Implement the functionality controlled by these.
Definition at line 139 of file JSSTaggerBase.h.
◆ m_UseHLTagger
bool JSSTaggerUtils::m_UseHLTagger {} |
|
private |
◆ m_useLSeries
bool JSSTaggerBase::m_useLSeries = false |
|
protectedinherited |
◆ m_useMassCut
bool JSSTaggerBase::m_useMassCut = false |
|
protectedinherited |
Flag to indicate if mass window cut is used.
Definition at line 94 of file JSSTaggerBase.h.
◆ m_useScoreCut
bool JSSTaggerBase::m_useScoreCut = false |
|
protectedinherited |
Flag to indicate if a discriminant score is used.
Definition at line 97 of file JSSTaggerBase.h.
◆ m_varHandleArraysDeclared
◆ m_vhka
◆ m_weightConfig
std::unique_ptr<TFile> JSSTaggerBase::m_weightConfig |
|
protectedinherited |
◆ m_weightConfigPath
std::string JSSTaggerBase::m_weightConfigPath |
|
protectedinherited |
◆ m_weightDecorationName
std::string JSSTaggerBase::m_weightDecorationName |
|
protectedinherited |
◆ m_weightFileName
std::string JSSTaggerBase::m_weightFileName |
|
protectedinherited |
◆ m_weightFlavors
std::string JSSTaggerBase::m_weightFlavors |
|
protectedinherited |
◆ m_weightHistogramName
std::string JSSTaggerBase::m_weightHistogramName |
|
protectedinherited |
◆ m_weightHistograms
std::map<std::string, std::unique_ptr<TH2D> > JSSTaggerBase::m_weightHistograms |
|
protectedinherited |
◆ m_wkpt
std::string JSSTaggerBase::m_wkpt |
|
protectedinherited |
The documentation for this class was generated from the following files:
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readParentKey
def retrieve(aClass, aKey=None)
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readTau4WTAKey
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readSplit23Key
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readECF1Key
std::string getTruthLabelStr(const xAOD::Jet &jet, asg::AcceptData &acceptData) const
void mean(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="")
Scalar phi() const
phi method
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readNtrk500Key
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readECFG311Key
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decValidEtaRangeKey
Scalar eta() const
pseudorapidity method
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decHLScoreKey
std::map< std::string, std::unique_ptr< TH2D > > m_efficiencyHistograms
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T > &t)
std::string m_configFile
Configuration file name.
TH2D MakeJetImage(TString TagImage, const xAOD::Jet *jet, std::vector< xAOD::JetConstituent > constituents) const override
TEnv m_configReader
TEnv instance to read config files.
Int_t FindBin(const TAxis *axis, const double x)
bool accumulate(AccumulateMap &map, std::vector< module_t > const &modules, FPGATrackSimMatrixAccumulator const &acc)
Accumulates an accumulator (e.g.
virtual StatusCode tag(const xAOD::Jet &jet) const =0
Decorate single jet with tagging info.
const float SF[NF]
Cross sections for Fluor.
StoreGateSvc_t m_evtStore
Pointer to StoreGate (event store by default)
std::vector< SG::VarHandleKeyArray * > m_vhka
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readTau2WTAKey
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readSplit12Key
float m_jetPtMin
Kinematic bounds for the jet - the units are controlled by m_ptGeV.
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decL3Key
Helper class to provide constant type-safe access to aux data.
ToolHandle< AthONNX::IJSSMLTool > m_MLBosonTagger_HL
bool DescendingPtSorterConstituents(const xAOD::JetConstituent p1, const xAOD::JetConstituent p2)
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decC2Key
std::string m_decorationName
Decoration name.
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readECFG331Key
#define ATH_MSG_VERBOSE(x)
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readConstScoreKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decValidPtRangeHighKey
bool isValid() const
Test to see if the link can be dereferenced.
const std::string & getCutName(unsigned int cutPosition) const
Get the name of a cut, based on the cut position (slow, avoid usage)
std::string m_calibArea
Location where config files live on cvmfs.
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decWeightKey
WriteDecorHandle keys for SF.
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readTruthLabelKey
Truth label ReadDecorHandle key.
virtual void setOwner(IDataHandleHolder *o)=0
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decSigeffSFKey
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decEfficiencyKey
ServiceHandle< StoreGateSvc > & evtStore()
The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decE3Key
Description of a calorimeter cluster.
bool m_ptGeV
Flag to indicate units pT is defined in Set to false by default.
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="")
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decTau32WTAKey
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readQwKey
Handle class for reading a decoration on an object.
std::pair< double, double > getSF(const xAOD::Jet &jet, const std::string &truthLabelStr) const
Get scale factor and efficiency.
std::map< std::string, double > GetJSSVars(const xAOD::Jet &jet) const override
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readECFG212Key
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decNTopoTowersKey
unsigned int getNCuts() const
Get the number of cuts defined.
Handle class for adding a decoration to an object.
std::string m_efficiencyHistogramName
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decTau42WTAKey
const std::string & name(Level lvl)
Convenience function for translating message levels to strings.
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decEffSFKey
StoreGateSvc_t m_detStore
Pointer to StoreGate (detector store by default)
bool m_calcSF
Flag to calculate scale factor.
int findPV() const
Find the PV (to be used for Ntrk)
virtual void renounce()=0
std::conditional< std::is_base_of< SG::VarHandleKeyArray, T >::value, VarHandleKeyArrayType, type2 >::type type
ElementLink implementation for ROOT usage.
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decD2Key
virtual StatusCode initialize() override
Initialize the tool.
asg::AcceptInfo m_acceptInfo
Object that stores the results for a jet.
std::string to_string(const DetectorType &type)
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decTau21WTAKey
WriteDecorHandle keys for JSS moments.
StatusCode ReadScaler() override
void decorateJSSRatios(const xAOD::JetContainer &jets) const
std::string PathResolverFindCalibFile(const std::string &logical_file_name)
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decValidKinRangeKey
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readTau1WTAKey
ReadDecorHandle keys for JSS moments.
std::map< std::string, std::unique_ptr< TH2D > > m_weightHistograms
void setCutResult(const std::string &cutName, bool cutResult)
Set the result of a cut, based on the cut name (safer)
virtual StatusCode decorate(const xAOD::JetContainer &jets) const =0
Decorate a jet collection without otherwise modifying it.
bool passKinRange(const xAOD::Jet &jet) const
Check if jet passes kinematic constraints.
void clear()
Clear all bits.
bool getCutResult(const std::string &cutName) const
Get the result of a cut, based on the cut name (safer)
void sort(typename std::reverse_iterator< DataModel_detail::iterator< DVL > > beg, typename std::reverse_iterator< DataModel_detail::iterator< DVL > > end, const Compare &comp)
Specialization of sort for DataVector/List.
#define ATH_MSG_WARNING(x)
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decValidPtRangeLowKey
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readTau3WTAKey
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readD2Key
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readECF2Key
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readECF3Key
StatusCode getConfigReader()
Get configReader StatusCode.
TypeEnum intToEnum(const int type)
int GetUnGroomTracks(const xAOD::Jet &jet, int indexPV) const
Retrieve Ntrk variable from the ungroomed parent jet.
SG::VarHandleKey & vhKey()
Return a non-const reference to the HandleKey.
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decConstScoreKey
ToolHandle< AthONNX::IJSSMLTool > m_MLBosonTagger
JSSTaggerBase(const std::string &name)
Default constructor - to be used in all derived classes.
4-vector of jet constituent at the scale used during jet finding.
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decL2Key
std::string m_weightFlavors
SG::WriteDecorHandleKey< xAOD::JetContainer > m_decNConstituentsKey
WriteDecorHandle keys.
constexpr int pow(int base, int exp) noexcept
setBGCode setTAP setLVL2ErrorBits bool
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>
const std::string & getCutDescription(const std::string &cutName) const
Get the description of a cut, based on the cut name.
std::string m_containerName
Configurable members.
std::map< std::string, std::vector< double > > m_scaler