Loading [MathJax]/extensions/tex2jax.js
 |
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 | GetQGConstScore (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_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 98 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 395 of file JSSTaggerBase.cxx.
425 float tau21_wta = -999.0;
426 float tau32_wta = -999.0;
427 float tau42_wta = -999.0;
429 float tau1_wta = readTau1WTA(
jet);
430 float tau2_wta = readTau2WTA(
jet);
431 float tau3_wta = readTau3WTA(
jet);
432 float tau4_wta = -999.0;
433 if(readTau4WTA.isAvailable()){
434 tau4_wta = readTau4WTA(
jet);
437 if ( tau1_wta > 1
e-8 ) {
438 tau21_wta = tau2_wta / tau1_wta;
442 if ( tau2_wta > 1
e-8 ) {
443 tau32_wta = tau3_wta / tau2_wta;
444 if(readTau4WTA.isAvailable()){
445 tau42_wta = tau4_wta / tau2_wta;
450 decTau21WTA(
jet) = tau21_wta;
451 decTau32WTA(
jet) = tau32_wta;
452 decTau42WTA(
jet) = tau42_wta;
459 float ECF1 = readECF1(
jet);
460 float ECF2 = readECF2(
jet);
461 float ECF3 = readECF3(
jet);
464 C2 = ECF3 * ECF1 /
std::pow( ECF2, 2.0 );
480 if(!accL2.isAvailable(
jet)){
481 if(readECFG331.isAvailable() && readECFG212.isAvailable()){
482 if(readECFG212(
jet) > 1
e-8){
483 L2 = readECFG331(
jet) /
pow(readECFG212(
jet), (3.0/2.0));
491 if(!accL3.isAvailable(
jet)){
492 if(readECFG331.isAvailable() && readECFG311.isAvailable()){
493 if(readECFG331(
jet) > 1
e-8){
494 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 349 of file JSSTaggerBase.cxx.
383 decValidPtRangeHigh(
jet) = acceptData.
getCutResult(
"ValidPtRangeHigh" );
384 decValidPtRangeLow(
jet) = acceptData.
getCutResult(
"ValidPtRangeLow" );
388 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 508 of file JSSTaggerBase.cxx.
512 bool decor_available = decTau21WTA.isAvailable();
513 if(decor_available)
return;
540 float tau21_wta = -999.0;
541 float tau32_wta = -999.0;
542 float tau42_wta = -999.0;
544 float tau1_wta = readTau1WTA(*
jet);
545 float tau2_wta = readTau2WTA(*
jet);
546 float tau3_wta = readTau3WTA(*
jet);
547 float tau4_wta = -999.0;
548 if(readTau4WTA.isAvailable()){
549 tau4_wta = readTau4WTA(*
jet);
552 if ( tau1_wta > 1
e-8 ) {
553 tau21_wta = tau2_wta / tau1_wta;
556 if ( tau2_wta > 1
e-8 ) {
557 tau32_wta = tau3_wta / tau2_wta;
558 if(readTau4WTA.isAvailable()){
559 tau42_wta = tau4_wta / tau2_wta;
563 decTau21WTA(*
jet) = tau21_wta;
564 decTau32WTA(*
jet) = tau32_wta;
565 decTau42WTA(*
jet) = tau42_wta;
572 float ECF1 = readECF1(*
jet);
573 float ECF2 = readECF2(*
jet);
574 float ECF3 = readECF3(*
jet);
577 C2 = ECF3 * ECF1 /
std::pow( ECF2, 2.0 );
592 if(!accL2.isAvailable(*
jet)){
593 if(readECFG331.isAvailable() && readECFG212.isAvailable()){
594 if(readECFG212(*
jet) > 1
e-8){
595 L2 = readECFG331(*
jet) /
pow(readECFG212(*
jet), (3.0/2.0));
602 if(!accL3.isAvailable(*
jet)){
603 if(readECFG331.isAvailable() && readECFG311.isAvailable()){
604 if(readECFG331(*
jet) > 1
e-8){
605 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 620 of file JSSTaggerBase.cxx.
626 ATH_MSG_WARNING(
"Unable to retrieve primary vertex container PrimaryVertices" );
630 for (
const auto *vx : *vxCont ) {
632 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 278 of file JSSTaggerBase.cxx.
283 std::string configPath;
288 else if (
m_calibArea.find(
"eos") != std::string::npos) {
298 int fSuccess = gSystem->GetPathInfo(configPath.c_str(), fStats);
301 return StatusCode::FAILURE;
304 ATH_MSG_DEBUG(
"Recommendations file was found : " << configPath );
307 if (
m_configReader.ReadFile( configPath.c_str(), EEnvLevel(0) ) ) {
308 ATH_MSG_ERROR(
"Error while reading config file : "<< configPath );
309 return StatusCode::FAILURE;
312 return StatusCode::SUCCESS;
◆ GetConstScore()
Implements IJSSTaggerUtils.
Definition at line 238 of file JSSTaggerUtils.cxx.
250 std::vector<xAOD::JetConstituent> constituents =
jet -> getConstituents().asSTLVector();
253 int MaxConstituents (100);
254 std::vector<xAOD::JetConstituent> constituentsForModel;
256 if( constituents.size() > 100 )
257 constituentsForModel = std::vector<xAOD::JetConstituent> (constituents.begin(), constituents.begin() + MaxConstituents);
259 constituentsForModel = constituents;
262 std::vector<const xAOD::CaloCluster*>
towers;
264 if (towersAcc.isAvailable(*
jet)){
266 std::vector<ElementLink<DataVector<xAOD::IParticle>>> towerLinks = towersAcc(*
jet);
267 for (
auto link_itr : towerLinks){
268 if (!link_itr.isValid())
continue;
276 for(
auto cnst : constituents){
277 m.push_back( cnst ->
m() );
278 pT.push_back( cnst ->
pt() );
279 eta.push_back( cnst ->
eta() );
280 phi.push_back( cnst ->
phi() );
281 E.push_back( cnst ->
e() );
283 std::vector<std::vector<float>> constituents_packed = {
m,
pT,
eta,
phi};
285 m.clear();
pT.clear();
eta.clear();
phi.clear();
287 m.push_back( cnst ->
m() );
288 pT.push_back( cnst ->
pt() );
289 eta.push_back( cnst ->
eta() );
290 phi.push_back( cnst ->
phi() );
291 E.push_back( cnst ->
e() );
293 std::vector<std::vector<float>> towers_packed = {
m,
pT,
eta,
phi};
296 std::vector<std::vector<float>> inputs_packed = {
297 constituents_packed.at(0), constituents_packed.at(1), constituents_packed.at(2), constituents_packed.at(3),
298 towers_packed.at(0), towers_packed.at(1), towers_packed.at(2), towers_packed.at(3),
302 if( (constituents.size() +
towers.size()) > 1 )
309 decNConstituents(*
jet) = constituents.size();
314 return StatusCode::SUCCESS;
◆ GetHLScore()
◆ GetImageScore()
Implements IJSSTaggerUtils.
Definition at line 172 of file JSSTaggerUtils.cxx.
185 std::vector<xAOD::JetConstituent> constituents =
jet -> getConstituents().asSTLVector();
188 int MaxConstituents (100);
189 std::vector<xAOD::JetConstituent> constituentsForModel;
191 if( constituents.size() > 100 )
192 constituentsForModel = std::vector<xAOD::JetConstituent> (constituents.begin(), constituents.begin() + MaxConstituents);
194 constituentsForModel = constituents;
197 std::vector<xAOD::JetConstituent> csts_charged = constituentsForModel;
198 csts_charged.erase( std::remove_if( csts_charged.begin(), csts_charged.end(),
200 const xAOD::FlowElement* ufo = dynamic_cast<const xAOD::FlowElement*>(constituent.rawConstituent());
201 return ufo -> signalType() != xAOD::FlowElement::SignalType::Charged;
202 }), csts_charged.end()) ;
205 std::vector<xAOD::JetConstituent> csts_neutral = constituentsForModel;
206 csts_neutral.erase( std::remove_if( csts_neutral.begin(), csts_neutral.end(),
208 const xAOD::FlowElement* ufo = dynamic_cast<const xAOD::FlowElement*>(constituent.rawConstituent());
209 return ufo -> signalType() != xAOD::FlowElement::SignalType::Neutral;
210 }), csts_neutral.end()) ;
213 std::vector<xAOD::JetConstituent> csts_combined = constituentsForModel;
214 csts_combined.erase( std::remove_if( csts_combined.begin(), csts_combined.end(),
216 const xAOD::FlowElement* ufo = dynamic_cast<const xAOD::FlowElement*>(constituent.rawConstituent());
217 return ufo -> signalType() != xAOD::FlowElement::SignalType::Combined;
218 }), csts_combined.end()) ;
225 std::vector<TH2D> Images = {ImageCharged, ImageNeutral, ImageCombined};
234 return StatusCode::SUCCESS;
◆ GetJSSVars()
| std::map< std::string, double > JSSTaggerUtils::GetJSSVars |
( |
const xAOD::Jet & |
jet | ) |
const |
|
overridevirtual |
Implements IJSSTaggerUtils.
Definition at line 486 of file JSSTaggerUtils.cxx.
488 std::map<std::string, double> JSSVars;
491 int nUngrTracks (-1);
493 if(readNtrk500.isAvailable()){
494 nUngrTracks = readNtrk500(
jet);
497 int pv_location =
findPV();
498 if(pv_location != -1)
503 JSSVars[
"pT"] =
jet.pt();
504 JSSVars[
"nTracks"] = nUngrTracks;
507 JSSVars[
"CNN"] = readConstScore(
jet);
534 JSSVars[
"Split12"] = readSplit12(
jet);
535 JSSVars[
"Split23"] = readSplit23(
jet);
538 JSSVars[
"D2"] = D2Acc(
jet);
541 JSSVars[
"Tau1_wta"] = readTau1_wta(
jet);
542 JSSVars[
"Tau2_wta"] = readTau2_wta(
jet);
543 JSSVars[
"Tau3_wta"] = readTau3_wta(
jet);
546 JSSVars[
"ECF1"] = readECF1(
jet);
547 JSSVars[
"ECF2"] = readECF2(
jet);
548 JSSVars[
"ECF3"] = readECF3(
jet);
551 JSSVars[
"Qw"] = readQw(
jet);
554 JSSVars[
"FoxWolfram0"] = FoxWolfram0Acc.withDefault(
jet, -99.);
555 JSSVars[
"FoxWolfram2"] = FoxWolfram2Acc.withDefault(
jet, -99.);
556 JSSVars[
"PlanarFlow"] = PlanarFlowAcc.withDefault(
jet, -99.);
557 JSSVars[
"Angularity"] = AngularityAcc.withDefault(
jet, -99.);
558 JSSVars[
"Aplanarity"] = AplanarityAcc.withDefault(
jet, -99.);
559 JSSVars[
"ZCut12"] = ZCut12Acc.withDefault(
jet, -99.);
560 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.
◆ GetQGConstScore()
Implements IJSSTaggerUtils.
Definition at line 318 of file JSSTaggerUtils.cxx.
328 std::vector<xAOD::JetConstituent> constituents =
jet -> getConstituents().asSTLVector();
332 std::vector<const xAOD::CaloCluster*>
towers;
334 if (towersAcc.isAvailable(*
jet)){
336 std::vector<ElementLink<DataVector<xAOD::IParticle>>> towerLinks = towersAcc(*
jet);
337 for (
auto link_itr : towerLinks){
338 if (!link_itr.isValid())
continue;
346 for(
auto cnst : constituents){
347 m.push_back( cnst ->
m() );
348 pT.push_back( cnst ->
pt() );
349 eta.push_back( cnst ->
eta() );
350 phi.push_back( cnst ->
phi() );
351 E.push_back( cnst ->
e() );
352 isTower.push_back(0.);
355 pz.push_back( cnst ->
pt() * std::sinh(cnst ->
eta()) );
359 m.push_back( cnst ->
m() );
360 pT.push_back( cnst ->
pt() );
361 eta.push_back( cnst ->
eta() );
362 phi.push_back( cnst ->
phi() );
363 E.push_back( cnst ->
e() );
364 isTower.push_back(1.);
367 pz.push_back( cnst ->
pt() * std::sinh(cnst ->
eta()) );
371 std::vector<std::vector<float>> features_packed;
373 for(
long unsigned int f=0;
f<
pT.size();
f++){
374 std::vector<float> features = {
m.at(
f),
pT.at(
f),
eta.at(
f),
phi.at(
f),
E.at(
f), isTower.at(
f),
px.at(
f),
py.at(
f),
pz.at(
f) };
375 features_packed.push_back(features);
379 std::sort(features_packed.begin(), features_packed.end(), [](
const auto&
i,
const auto& j) { return i.at(1) > j.at(1); });
386 float sum_features_pT = sqrt( sum_features_px*sum_features_px + sum_features_py*sum_features_py);
387 float sum_features_eta =
std::asinh( sum_features_pz /
Clip(sum_features_pT) );
388 float sum_features_phi = std::atan2( sum_features_py, sum_features_px );
391 std::vector<std::vector<float>> const_vars;
392 std::vector<std::vector<std::vector<float>>> inter_vars;
394 for(
auto const &feature_i : features_packed){
397 float log_pT =
log(
Clip(feature_i.at(1) / sum_features_pT) );
398 float log_E =
log(
Clip(feature_i.at(4) / sum_features_pT) );
399 float eta = feature_i.at(2) - sum_features_eta;
400 float phi = feature_i.at(3) - sum_features_phi;
401 float DR = sqrt((feature_i.at(2) - sum_features_eta)*(feature_i.at(2) - sum_features_eta) + (feature_i.at(3) - sum_features_phi)*(feature_i.at(3) - sum_features_phi));
402 float log_m =
log(
Clip(feature_i.at(0)));
403 float type = feature_i.at(5);
405 std::vector<float> vars = {log_pT, log_E,
eta,
phi, DR, log_m,
type};
406 const_vars.push_back(vars);
409 std::vector<std::vector<float>> inter_vars_int;
410 for(
auto const &feature_j : features_packed){
413 float delta = sqrt((feature_i.at(2)-feature_j.at(2))*(feature_i.at(2)-feature_j.at(2)) + (feature_i.at(3)-feature_j.at(3))*(feature_i.at(3)-feature_j.at(3)));
414 float min = feature_i.at(1) != feature_j.at(1) ?
std::min(feature_i.at(1), feature_j.at(1)): 0.;
415 float mass2 = (feature_i.at(4)/sum_features_pT + feature_j.at(4)/sum_features_pT) * (feature_i.at(4)/sum_features_pT + feature_j.at(4)/sum_features_pT);
416 mass2 -= (feature_i.at(6)/sum_features_pT + feature_j.at(6)/sum_features_pT) * (feature_i.at(6)/sum_features_pT + feature_j.at(6)/sum_features_pT);
417 mass2 -= (feature_i.at(7)/sum_features_pT + feature_j.at(7)/sum_features_pT) * (feature_i.at(7)/sum_features_pT + feature_j.at(7)/sum_features_pT);
418 mass2 -= (feature_i.at(8)/sum_features_pT + feature_j.at(8)/sum_features_pT) * (feature_i.at(8)/sum_features_pT + feature_j.at(8)/sum_features_pT);
421 float log_delta =
log(
Clip(delta));
422 float log_mindelta =
log(
Clip(
min * delta / sum_features_pT));
423 float min_over_pT =
min / (feature_i.at(1) + feature_j.at(1));
424 float log_mass =
log(
Clip(mass2));
427 if(feature_i==feature_j){
434 std::vector<float> vars = { log_delta,
439 inter_vars_int.push_back(vars);
443 inter_vars.push_back(inter_vars_int);
447 if( (constituents.size() +
towers.size()) > 1 )
455 return StatusCode::SUCCESS;
◆ 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 756 of file JSSTaggerBase.cxx.
765 const double WtoZmassShift = 10803;
769 if ( logmOverPt > 0 ) logmOverPt = 0;
786 std::string signal_truthLabel=
"";
788 signal_truthLabel=
"t_qqb";
790 signal_truthLabel=
"V_qq";
792 signal_truthLabel=
"t";
799 return std::make_pair( 1.0,
eff );
803 ATH_MSG_DEBUG(
"(pt, m/pt) (" <<
jet.pt()/1.e3 <<
", " <<
jet.m()/
jet.pt() <<
") is out of range for SF calculation. Returning 1.0" );
804 return std::make_pair( 1.0, 1.0 );
807 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 814 of file JSSTaggerBase.cxx.
816 std::string truthLabelStr;
827 truthLabelStr =
"t_qqb";
840 truthLabelStr =
"t_";
844 truthLabelStr =
"V_qq_";
848 truthLabelStr =
"q_";
853 truthLabelStr +=
"passMpassD2";
857 truthLabelStr +=
"failMpassD2";
861 truthLabelStr +=
"passMfailD2";
865 truthLabelStr +=
"failMfailD2";
879 truthLabelStr =
"V_qq";
900 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 644 of file JSSTaggerBase.cxx.
648 const xAOD::Jet * ungroomedJet =
nullptr;
649 if ( readParent.isAvailable() ) {
651 if ( linkToUngroomed.
isValid() ) {
652 ungroomedJet = *linkToUngroomed;
656 if ( acc_Ntrk.isAvailable(*ungroomedJet) ) {
658 const std::vector<int>& NTrkPt500 = acc_Ntrk(*ungroomedJet);
660 int jet_ntrk = NTrkPt500.at(indexPV);
665 ATH_MSG_ERROR(
"WARNING: Unable to retrieve Ntrk of the ungroomed parent jet. Please make sure this variable is in your derivations!!!");
669 ATH_MSG_ERROR(
"WARNING: Unable to retrieve the parent ungroomed jet. Please make sure this variable is in your derivations!!!");
673 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 682 of file JSSTaggerBase.cxx.
683 if ( !
m_calcSF )
return StatusCode::SUCCESS;
687 float sigeffSF = 1.0;
697 sigeffSF =
getSF(
jet,
"t_qqb").first;
699 sigeffSF =
getSF(
jet,
"V_qq").first;
748 decEffSF(
jet) = effSF;
749 decSigeffSF(
jet) = sigeffSF;
751 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 53 of file JSSTaggerUtils.cxx.
71 std::string ConstTaggerFileName =
m_configReader.GetValue(
"ConstTaggerFileName",
"aaa");
72 std::string HLTaggerFileName =
m_configReader.GetValue(
"HLTaggerFileName",
"aaa");
76 std::string ModelPath =
"";
80 else if (
m_calibArea.find(
"eos") != std::string::npos) {
81 ModelPath = (ConstTaggerFileName);
86 ATH_MSG_INFO(
"JSSTaggerUtils::MLBosonTagger()" <<
" + ModelPath " << ModelPath );
92 ATH_MSG_INFO(
"JSSTaggerUtils::MLBosonTagger() read value from config" );
116 std::string ModelPath =
"";
120 else if (
m_calibArea.find(
"eos") != std::string::npos) {
121 ModelPath = (HLTaggerFileName);
126 ATH_MSG_INFO(
"JSSTaggerUtils::MLBosonTagger()" <<
" + ModelPath " << ModelPath );
135 ATH_MSG_INFO(
"JSSTaggerUtils::MLBosonTagger() read value from config" );
143 ATH_MSG_INFO(
"Decorators that will be attached to jet :" );
160 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 610 of file JSSTaggerUtils.cxx.
613 double eta (-99.),
phi (-99.),
pT (-99.),
z (-99.);
614 int BinEta (-99), BinPhi (-99);
617 return sum + cst.pt(); };
618 double SumPT =
std::accumulate( constituents.begin(), constituents.end(), 0., cst_pT) ;
620 auto Image = std::make_unique<TH2D>(
"Image_" + TagImage,
"Image_" + TagImage,
623 for(
auto& cst : constituents ){
641 double current_z = Image -> GetBinContent(
x,
y );
644 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 903 of file JSSTaggerBase.cxx.
904 ATH_MSG_INFO(
"After tagging, you will have access to the following cuts as an asg::AcceptData : (<NCut>) <cut> : <description>)" );
906 for (
int iCut=0; iCut < nCuts; iCut++ ) {
907 std::string cut_string =
"";
◆ ReadScaler()
| StatusCode JSSTaggerUtils::ReadScaler |
( |
| ) |
|
|
overridevirtual |
Implements IJSSTaggerUtils.
Definition at line 566 of file JSSTaggerUtils.cxx.
569 std::vector<std::string> vars_list = {
570 "pT",
"CNN",
"D2",
"nTracks",
"ZCut12",
571 "Tau1_wta",
"Tau2_wta",
"Tau3_wta",
572 "KtDR",
"Split12",
"Split23",
573 "ECF1",
"ECF2",
"ECF3",
574 "Angularity",
"FoxWolfram0",
"FoxWolfram2",
575 "Aplanarity",
"PlanarFlow",
"Qw",
579 for(
const std::string &
var : vars_list){
580 std::string s_mean =
var +
"_mean";
581 std::string s_std =
var +
"_std";
585 if(
mean==-99. && std==-99.){
586 ATH_MSG_ERROR(
"ERROR: one of the parameter for " <<
var <<
" is missing, please, double check the config!!!");
587 return StatusCode::FAILURE;
589 else if(
mean==-99. || std==-99.){
590 ATH_MSG_ERROR(
"ERROR: parameters for " <<
var <<
" are missing, please, double check the config!!!");
591 return StatusCode::FAILURE;
606 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 317 of file JSSTaggerBase.cxx.
330 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 113 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_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 131 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 89 of file JSSTaggerBase.h.
◆ m_readConstScoreKey
◆ 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 137 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 92 of file JSSTaggerBase.h.
◆ m_useScoreCut
| bool JSSTaggerBase::m_useScoreCut = false |
|
protectedinherited |
Flag to indicate if a discriminant score is used.
Definition at line 95 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)
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)
const std::string & key() const
Return the StoreGate ID for the referenced object.
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
def asinh(x)
helper methods ---------------------------------------------------------—
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
accumulate
Update flags based on parser line args.
SG::ReadDecorHandleKey< xAOD::JetContainer > m_readTau1WTAKey
ReadDecorHandle keys for JSS moments.
std::map< std::string, std::unique_ptr< TH2D > > m_weightHistograms
bool empty() const noexcept
void setCutResult(const std::string &cutName, bool cutResult)
Set the result of a cut, based on the cut name (safer)
StatusCode initialize(bool used=true)
If this object is used as a property, then this should be called during the initialize phase.
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_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
1.8.18