|
ATLAS Offline Software
|
#include <TCPufitFex.h>
|
| TCPufitFex (const std::string &name, ISvcLocator *pSvcLocator) |
| Constructor. More...
|
|
virtual StatusCode | initialize () override |
| Initialize the fex. More...
|
|
virtual StatusCode | execute (const EventContext &context) const override |
| Run the algorithm. More...
|
|
virtual StatusCode | sysInitialize () override |
| Override sysInitialize. More...
|
|
virtual bool | isClonable () const override |
| Specify if the algorithm is clonable. More...
|
|
virtual unsigned int | cardinality () const override |
| Cardinality (Maximum number of clones that can exist) special value 0 means that algorithm is reentrant. More...
|
|
virtual StatusCode | sysExecute (const EventContext &ctx) override |
| Execute an algorithm. More...
|
|
virtual const DataObjIDColl & | extraOutputDeps () const override |
| Return the list of extra output dependencies. More...
|
|
virtual bool | filterPassed (const EventContext &ctx) const |
|
virtual void | setFilterPassed (bool state, const EventContext &ctx) const |
|
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 | sysStart () override |
| Handle START transition. More...
|
|
virtual std::vector< Gaudi::DataHandle * > | inputHandles () const override |
| Return this algorithm's input handles. More...
|
|
virtual std::vector< Gaudi::DataHandle * > | outputHandles () const override |
| Return this algorithm's output handles. More...
|
|
Gaudi::Details::PropertyBase & | declareProperty (Gaudi::Property< T > &t) |
|
Gaudi::Details::PropertyBase * | declareProperty (const std::string &name, SG::VarHandleKey &hndl, const std::string &doc, const SG::VarHandleKeyType &) |
| Declare a new Gaudi property. More...
|
|
Gaudi::Details::PropertyBase * | declareProperty (const std::string &name, SG::VarHandleBase &hndl, const std::string &doc, const SG::VarHandleType &) |
| Declare a new Gaudi property. More...
|
|
Gaudi::Details::PropertyBase * | declareProperty (const std::string &name, SG::VarHandleKeyArray &hndArr, const std::string &doc, const SG::VarHandleKeyArrayType &) |
|
Gaudi::Details::PropertyBase * | declareProperty (const std::string &name, T &property, const std::string &doc, const SG::NotHandleType &) |
| Declare a new Gaudi property. More...
|
|
Gaudi::Details::PropertyBase * | declareProperty (const std::string &name, T &property, const std::string &doc="none") |
| Declare a new Gaudi property. More...
|
|
void | updateVHKA (Gaudi::Details::PropertyBase &) |
|
MsgStream & | msg () const |
|
MsgStream & | msg (const MSG::Level lvl) const |
|
bool | msgLvl (const MSG::Level lvl) const |
|
Definition at line 40 of file TCPufitFex.h.
◆ StoreGateSvc_t
◆ TCPufitFex()
HLT::MET::TCPufitFex::TCPufitFex |
( |
const std::string & |
name, |
|
|
ISvcLocator * |
pSvcLocator |
|
) |
| |
◆ cardinality()
unsigned int AthReentrantAlgorithm::cardinality |
( |
| ) |
const |
|
overridevirtualinherited |
Cardinality (Maximum number of clones that can exist) special value 0 means that algorithm is reentrant.
Override this to return 0 for reentrant algorithms.
Override this to return 0 for reentrant algorithms.
Definition at line 55 of file AthReentrantAlgorithm.cxx.
◆ declareGaudiProperty() [1/4]
specialization for handling Gaudi::Property<SG::VarHandleKeyArray>
Definition at line 170 of file AthCommonDataStore.h.
175 hndl.documentation());
◆ declareGaudiProperty() [2/4]
specialization for handling Gaudi::Property<SG::VarHandleKey>
Definition at line 156 of file AthCommonDataStore.h.
161 hndl.documentation());
◆ declareGaudiProperty() [3/4]
specialization for handling Gaudi::Property<SG::VarHandleBase>
Definition at line 184 of file AthCommonDataStore.h.
189 hndl.documentation());
◆ declareGaudiProperty() [4/4]
◆ declareProperty() [1/6]
Declare a new Gaudi property.
- Parameters
-
name | Name of the property. |
hndl | Object holding the property value. |
doc | Documentation string for the property. |
This is the version for types that derive from SG::VarHandleBase
. The property value object is put on the input and output lists as appropriate; then we forward to the base class.
Definition at line 245 of file AthCommonDataStore.h.
250 this->declare(hndl.
vhKey());
251 hndl.
vhKey().setOwner(
this);
◆ declareProperty() [2/6]
Declare a new Gaudi property.
- Parameters
-
name | Name of the property. |
hndl | Object holding the property value. |
doc | Documentation string for the property. |
This is the version for types that derive from SG::VarHandleKey
. The property value object is put on the input and output lists as appropriate; then we forward to the base class.
Definition at line 221 of file AthCommonDataStore.h.
◆ declareProperty() [3/6]
◆ declareProperty() [4/6]
Declare a new Gaudi property.
- Parameters
-
name | Name of the property. |
property | Object holding the property value. |
doc | Documentation string for the property. |
This is the generic version, for types that do not derive from SG::VarHandleKey
. It just forwards to the base class version of declareProperty
.
Definition at line 333 of file AthCommonDataStore.h.
◆ declareProperty() [5/6]
Declare a new Gaudi property.
- Parameters
-
name | Name of the property. |
property | Object holding the property value. |
doc | Documentation string for the property. |
This dispatches to either the generic declareProperty
or the one for VarHandle/Key/KeyArray.
Definition at line 352 of file AthCommonDataStore.h.
◆ declareProperty() [6/6]
◆ detStore()
◆ evtStore() [1/2]
◆ evtStore() [2/2]
◆ execute()
StatusCode HLT::MET::FexBase::execute |
( |
const EventContext & |
context | ) |
const |
|
overridevirtualinherited |
Run the algorithm.
Definition at line 60 of file FexBase.cxx.
63 ATH_MSG_ERROR(
"Base class was not initialised! This means that the "
64 <<
"derived class was not correctly written!");
65 return StatusCode::FAILURE;
70 auto metCont = std::make_unique<xAOD::TrigMissingETContainer>();
71 auto metContAux = std::make_unique<xAOD::TrigMissingETAuxContainer>();
72 metCont->setStore(metContAux.get() );
75 metCont->push_back(std::make_unique<xAOD::TrigMissingET>());
91 MonGroupBuilder monitors;
105 for (std::size_t ii = 0; ii <
met->getNumberOfComponents(); ++ii)
107 <<
met->nameOfComponent(ii) <<
", "
108 <<
met->statusComponent(ii) <<
", "
109 << METComponent(ii, *
met) );
114 CHECK( handle.record( std::move(metCont), std::move(metContAux) ) );
115 return StatusCode::SUCCESS;
◆ 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
◆ extraOutputDeps()
const DataObjIDColl & AthReentrantAlgorithm::extraOutputDeps |
( |
| ) |
const |
|
overridevirtualinherited |
Return the list of extra output dependencies.
This list is extended to include symlinks implied by inheritance relations.
Definition at line 79 of file AthReentrantAlgorithm.cxx.
86 return Algorithm::extraOutputDeps();
◆ fillMET()
Calculate and fill the output MET value.
- Parameters
-
met | The object to fill |
context | The event context |
monitors[out] | Any extra monitors to fill |
Implements HLT::MET::FexBase.
Definition at line 44 of file TCPufitFex.cxx.
51 iclus->e(), iclus->eta(), iclus->phi());
62 for (PufitGrid &grid : gridset.grids)
68 const PufitGrid &grid = gridset[
idx];
71 PufitUtils::CovarianceSum pileupSum;
72 std::vector<SignedKinematics> masked;
80 pileupSum.add(tower,
sigma);
83 masked.push_back(tower);
91 sum.fillMETComponent(1,
met);
100 pileupSum.covariance,
107 for (
const SignedKinematics &kin : corrections)
111 return StatusCode::SUCCESS;
◆ filterPassed()
virtual bool AthReentrantAlgorithm::filterPassed |
( |
const EventContext & |
ctx | ) |
const |
|
inlinevirtualinherited |
◆ flagMET()
|
protectedvirtualinherited |
Flag suspicious values in the output MET.
- Parameters
-
Definition at line 118 of file FexBase.cxx.
123 for (std::size_t ii = 0; ii <
met.getNumberOfComponents(); ++ii) {
124 METComponent component(ii,
met);
125 if (component.sumEt > 0 &&
130 met.setStatusComponent(ii, component.status);
132 METComponent total(
met);
133 total.status |= overall;
134 if (total.sumEt > 0 &&
137 met.setFlag(total.status);
138 return StatusCode::SUCCESS;
◆ initialize()
StatusCode HLT::MET::TCPufitFex::initialize |
( |
| ) |
|
|
overridevirtual |
Initialize the fex.
Definition at line 38 of file TCPufitFex.cxx.
42 return StatusCode::FAILURE;
◆ initializeBase()
StatusCode HLT::MET::FexBase::initializeBase |
( |
const std::vector< std::string > & |
componentNames | ) |
|
|
protectedinherited |
Initialize the base class.
- Parameters
-
componentNames | The names of any components in the output MET object. |
The xAOD::TrigMissingET writes out a main 'met' value but can also write out extra MET components - either for debugging, monitoring or further use offline. The object has to be initialized with this names, they cannot be extended after this, so this information must be available for the base class to use.
Definition at line 48 of file FexBase.cxx.
57 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.
◆ isClonable()
bool AthReentrantAlgorithm::isClonable |
( |
| ) |
const |
|
overridevirtualinherited |
Specify if the algorithm is clonable.
Reentrant algorithms are clonable.
Reimplemented in Simulation::BeamEffectsAlg, InDet::SiTrackerSpacePointFinder, InDet::SCT_Clusterization, InDet::SiSPSeededTrackFinder, SCTRawDataProvider, InDet::GNNSeedingTrackMaker, SCT_PrepDataToxAOD, RoIBResultToxAOD, InDet::SiSPGNNTrackMaker, SCT_CablingCondAlgFromCoraCool, SCT_ReadCalibDataTestAlg, SCT_CablingCondAlgFromText, SCT_ReadCalibChipDataTestAlg, SCT_TestCablingAlg, SCT_ConfigurationConditionsTestAlg, ITkPixelCablingAlg, ITkStripCablingAlg, SCTEventFlagWriter, SCT_ConditionsSummaryTestAlg, SCT_ModuleVetoTestAlg, SCT_MonitorConditionsTestAlg, SCT_LinkMaskingTestAlg, SCT_MajorityConditionsTestAlg, SCT_RODVetoTestAlg, SCT_SensorsTestAlg, SCT_TdaqEnabledTestAlg, SCT_SiliconConditionsTestAlg, SCTSiLorentzAngleTestAlg, SCT_ByteStreamErrorsTestAlg, SCT_ConditionsParameterTestAlg, SCT_FlaggedConditionTestAlg, SCT_StripVetoTestAlg, SCT_RawDataToxAOD, and SCTSiPropertiesTestAlg.
Definition at line 44 of file AthReentrantAlgorithm.cxx.
◆ monitor()
|
protectedvirtualinherited |
Add monitor variables from an xAOD::TrigMissingET object.
- Parameters
-
met | The object to fill from |
monitors | Extra variables to be monitored |
Definition at line 141 of file FexBase.cxx.
145 METComponent metComponent(
met);
147 monitors.increaseCapacity(16,
true);
152 monitors.add(
Monitored::Scalar(
"EF_MET_log", toLog10GeV(metComponent.met() ) ) );
153 monitors.add(
Monitored::Scalar(
"EF_ME_log", toLog10GeV(metComponent.magnitude() ) ) );
154 monitors.add(
Monitored::Scalar(
"EF_SumEt_log", toLog10GeV(metComponent.sumEt) ) );
155 monitors.add(
Monitored::Scalar(
"EF_SumE_log", toLog10GeV(metComponent.sumE, -9e9) ) );
160 monitors.add(
Monitored::Scalar(
"EF_ME_lin", toLinGeV(metComponent.magnitude() ) ) );
164 "EF_XS", toLinGeV(metComponent.met() ) / toLinGeV(metComponent.sumEt, 1) ) );
166 return StatusCode::SUCCESS;
◆ msg() [1/2]
◆ msg() [2/2]
◆ 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.
◆ renounce()
◆ renounceArray()
◆ setFilterPassed()
virtual void AthReentrantAlgorithm::setFilterPassed |
( |
bool |
state, |
|
|
const EventContext & |
ctx |
|
) |
| const |
|
inlinevirtualinherited |
◆ sysExecute()
StatusCode AthReentrantAlgorithm::sysExecute |
( |
const EventContext & |
ctx | ) |
|
|
overridevirtualinherited |
Execute an algorithm.
We override this in order to work around an issue with the Algorithm base class storing the event context in a member variable that can cause crashes in MT jobs.
Definition at line 67 of file AthReentrantAlgorithm.cxx.
69 return Gaudi::Algorithm::sysExecute (ctx);
◆ sysInitialize()
StatusCode AthReentrantAlgorithm::sysInitialize |
( |
| ) |
|
|
overridevirtualinherited |
Override sysInitialize.
Override sysInitialize from the base class.
Loop through all output handles, and if they're WriteCondHandles, automatically register them and this Algorithm with the CondSvc
Scan through all outputHandles, and if they're WriteCondHandles, register them with the CondSvc
Reimplemented from AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >.
Reimplemented in InputMakerBase, and HypoBase.
Definition at line 96 of file AthReentrantAlgorithm.cxx.
107 if ( cs.retrieve().isFailure() ) {
109 return StatusCode::SUCCESS;
111 if (cs->regHandle(
this,*
h).isFailure()) {
112 sc = StatusCode::FAILURE;
113 ATH_MSG_ERROR(
"unable to register WriteCondHandle " <<
h->fullKey()
◆ sysStart()
Handle START transition.
We override this in order to make sure that conditions handle keys can cache a pointer to the conditions container.
◆ updateVHKA()
◆ m_baseInitialised
bool HLT::MET::FexBase::m_baseInitialised {false} |
|
privateinherited |
Whether or not this class has been correctly initialised.
Definition at line 143 of file FexBase.h.
◆ m_caloNoise
Gaudi::Property<float> HLT::MET::TCPufitFex::m_caloNoise |
|
private |
Initial value:{
this, "CaloNoise", 50, "The coefficient of the noise term in the "
"calorimeter resolution estimate [MeV]"}
The coefficient of the noise term in the calo resolution estimate.
Definition at line 86 of file TCPufitFex.h.
◆ m_caloStoch
Gaudi::Property<float> HLT::MET::TCPufitFex::m_caloStoch |
|
private |
Initial value:{
this, "CaloStochastic", 15.81, "The coefficient of the stochastic "
"term in the calorimeter resolution estimate [MeV^1/2]"}
The coefficient of the stochastic term in the calo resolution estimate.
Definition at line 90 of file TCPufitFex.h.
◆ m_clusterKey
Initial value:{
this, "ClusterName", "clusters", "Input cluster collection"}
Input clusters.
Definition at line 67 of file TCPufitFex.h.
◆ m_componentNames
std::vector<std::string> HLT::MET::FexBase::m_componentNames |
|
protectedinherited |
The names of the output MET components.
Definition at line 72 of file FexBase.h.
◆ m_constraintImportance
Gaudi::Property<float> HLT::MET::TCPufitFex::m_constraintImportance |
|
private |
Initial value:{
this, "ConstraintImportance", 1,
"The relative importance of the two constraints to the fit"}
The relative constraint importance.
Definition at line 94 of file TCPufitFex.h.
◆ m_detStore
◆ m_evtStore
◆ m_extendedExtraObjects
DataObjIDColl AthReentrantAlgorithm::m_extendedExtraObjects |
|
privateinherited |
Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.
Empty if no symlinks were found.
Definition at line 153 of file AthReentrantAlgorithm.h.
◆ m_maxComponentMETSumEtRatio
Gaudi::Property<float> HLT::MET::FexBase::m_maxComponentMETSumEtRatio |
|
privateinherited |
Initial value:{
"MaxComponentMETSumEtRatio", 1.,
"The maximum MET/SumEt ratio per component"}
Definition at line 132 of file FexBase.h.
◆ m_maxEta
Gaudi::Property<float> HLT::MET::TCPufitFex::m_maxEta |
|
private |
Initial value:{
this, "MaxEta", 5, "The maximum eta range"}
The eta range of the grid.
Definition at line 73 of file TCPufitFex.h.
◆ m_maxGlobalMETSumEtRatio
Gaudi::Property<float> HLT::MET::FexBase::m_maxGlobalMETSumEtRatio |
|
privateinherited |
Initial value:{
"MaxGlobalMETSumEtRatio", 1.,
"The maximum MET/SumEt ratio for the total value"}
Definition at line 135 of file FexBase.h.
◆ m_metContainerKey
Initial value:{
this, "METContainerKey", "HLT_MET", "Name of the output MET object"}
The output MET object.
Definition at line 127 of file FexBase.h.
◆ m_monTool
Initial value:{
this, "MonTool", "", "Monitoring tool"}
The monitoring tool.
Definition at line 130 of file FexBase.h.
◆ m_nEtaBins
Gaudi::Property<std::size_t> HLT::MET::TCPufitFex::m_nEtaBins |
|
private |
Initial value:{
this, "NEtaBins", 14, "The number of eta bins"}
The number of bins in eta.
Definition at line 76 of file TCPufitFex.h.
◆ m_nPhiBins
Gaudi::Property<std::size_t> HLT::MET::TCPufitFex::m_nPhiBins |
|
private |
Initial value:{
this, "NPhiBins", 8, "The number of phi bins"}
The number of bins in phi.
Definition at line 79 of file TCPufitFex.h.
◆ m_nSigma
Gaudi::Property<float> HLT::MET::TCPufitFex::m_nSigma |
|
private |
Initial value:{
this, "NSigma", 5, "Set the threshold at mean + NSigma*variance"}
The sigma threshold.
Definition at line 70 of file TCPufitFex.h.
◆ m_trimFraction
Gaudi::Property<float> HLT::MET::TCPufitFex::m_trimFraction |
|
private |
Initial value:{
this, "TrimFraction", 0.9,
"The fraction of bins to use when calculating the mean and variance"}
The trimming fraction.
Definition at line 82 of file TCPufitFex.h.
◆ m_varHandleArraysDeclared
◆ m_vhka
The documentation for this class was generated from the following files:
Eigen::VectorXd pufit(const Eigen::Vector2d &pileupSum, const Eigen::Matrix2d &pileupCovariance, const Eigen::VectorXd &towerExpectations, const Eigen::VectorXd &towerVariances, const Eigen::VectorXd &correctionDirections, double constraintImportance)
Perform the pile-up fit.
SG::ReadHandleKey< xAOD::CaloClusterContainer > m_clusterKey
Input clusters.
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="")
Gaudi::Property< float > m_caloNoise
The coefficient of the noise term in the calo resolution estimate.
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T > &t)
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)
StoreGateSvc_t m_evtStore
Pointer to StoreGate (event store by default)
std::vector< SG::VarHandleKeyArray * > m_vhka
bool msgLvl(const MSG::Level lvl) const
std::vector< std::string > m_componentNames
The names of the output MET components.
Gaudi::Property< float > m_maxComponentMETSumEtRatio
virtual StatusCode fillMET(xAOD::TrigMissingET &met, const EventContext &context, MonGroupBuilder &monitors) const =0
Calculate and fill the output MET value.
virtual void setOwner(IDataHandleHolder *o)=0
@ NoDisplacement
The grid is not shifted.
SG::ReadCondHandle< T > makeHandle(const SG::ReadCondHandleKey< T > &key, const EventContext &ctx=Gaudi::Hive::currentContext())
GridDisplacement
Enum to describe the positioning of the grid.
bool m_baseInitialised
Whether or not this class has been correctly initialised.
ToolHandle< GenericMonitoringTool > m_monTool
The monitoring tool.
virtual StatusCode monitor(const xAOD::TrigMissingET &met, MonGroupBuilder &monitors) const
Add monitor variables from an xAOD::TrigMissingET object.
virtual std::vector< Gaudi::DataHandle * > outputHandles() const override
Return this algorithm's output handles.
SG::WriteHandleKey< xAOD::TrigMissingETContainer > m_metContainerKey
The output MET object.
::StatusCode StatusCode
StatusCode definition for legacy code.
Gaudi::Property< float > m_maxGlobalMETSumEtRatio
GridDisplacement selectGrid(const PufitGridSet &grids)
Select the grid with the highest masked sumEt.
#define CHECK(...)
Evaluate an expression and check for errors.
virtual StatusCode flagMET(xAOD::TrigMissingET &met) const
Flag suspicious values in the output MET.
StoreGateSvc_t m_detStore
Pointer to StoreGate (detector store by default)
Gaudi::Property< float > m_trimFraction
The trimming fraction.
static SignedKinematics fromEnergyEtaPhi(double energy, double eta, double phi)
Factory function to construct from energy, eta, phi (massless)
DataObjIDColl m_extendedExtraObjects
Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.
virtual void renounce()=0
std::conditional< std::is_base_of< SG::VarHandleKeyArray, T >::value, VarHandleKeyArrayType, type2 >::type type
void trimmedMeanAndVariance(const std::vector< double > &sorted, double trimFraction, double &mean, double &variance)
Calculate the trimmed mean and variance for a vector of tower sumEts.
Gaudi::Property< float > m_constraintImportance
The relative constraint importance.
virtual StatusCode sysInitialize() override
Override sysInitialize.
@ ComponentBigMEtSEtRatio
Gaudi::Property< float > m_caloStoch
The coefficient of the stochastic term in the calo resolution estimate.
#define ATH_MSG_WARNING(x)
Class holding the Missing ET trigger fex results.
Gaudi::Property< float > m_nSigma
The sigma threshold.
SG::VarHandleKey & vhKey()
Return a non-const reference to the HandleKey.
Declare a monitored scalar variable.
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>