|
ATLAS Offline Software
|
Compute Enhanced Bias weights for the events.
More...
#include <EnhancedBiasWeightCompAlg.h>
|
| EnhancedBiasWeightCompAlg (const std::string &name, ISvcLocator *svcLoc) |
|
virtual StatusCode | initialize () override |
|
virtual StatusCode | start () override |
|
virtual StatusCode | execute (const EventContext &context) const override |
|
virtual StatusCode | stop () override |
|
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 |
|
|
SG::ReadHandleKey< TrigConf::HLTMenu > | m_HLTMenuKey {this, "HLTTriggerMenu", "DetectorStore+HLTTriggerMenu", "HLT Menu"} |
|
PublicToolHandle< Trig::TrigDecisionTool > | m_tdt {this, "TrigDecisionTool", "Trig::TrigDecisionTool/TrigDecisionTool", "TrigDecisionTool"} |
|
SG::ReadHandleKey< TrigCompositeUtils::DecisionContainer > | m_finalDecisionKey { this, "FinalDecisionKey", "", "Final stage of all decisions" } |
|
SG::ReadCondHandleKey< TrigConf::HLTPrescalesSet > | m_HLTPrescaleSetInputKey {this, "HLTPrescales", "HLTPrescales", "HLT prescales set"} |
|
SG::ReadCondHandleKey< TrigConf::L1PrescalesSet > | m_L1PrescaleSetInputKey {this, "L1Prescales", "L1Prescales", "L1 prescales set"} |
|
Gaudi::Property< std::map< std::string, std::vector< std::string > > > | m_chainToHLTSeed {this, "ChainToItemMap", {}, "Map for HLT seeded chains" } |
|
std::vector< HLT::Identifier > | m_EBChainIds |
| Available EB chains' IDs from HLT Menu. More...
|
|
std::map< HLT::Identifier, std::vector< std::string > > | m_EBChainIdToItem |
| L1 items for EB chains from HLT Menu. More...
|
|
std::string m_outputFilename | ATLAS_THREAD_SAFE |
| Name of xml output file. More...
|
|
tbb::concurrent_vector< std::pair< double, bool > > m_ebWeights | ATLAS_THREAD_SAFE |
| List of calculated weight + unbiased flag pairs to be saved in xml. More...
|
|
tbb::concurrent_unordered_map< int, double > m_eventToWeight | ATLAS_THREAD_SAFE |
| Map of event number to index in m_ebWeights vector to be saved in xml. More...
|
|
DataObjIDColl | m_extendedExtraObjects |
| Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks. More...
|
|
StoreGateSvc_t | m_evtStore |
| Pointer to StoreGate (event store by default) More...
|
|
StoreGateSvc_t | m_detStore |
| Pointer to StoreGate (detector store by default) More...
|
|
std::vector< SG::VarHandleKeyArray * > | m_vhka |
|
bool | m_varHandleArraysDeclared |
|
Compute Enhanced Bias weights for the events.
Definition at line 28 of file EnhancedBiasWeightCompAlg.h.
◆ StoreGateSvc_t
◆ EnhancedBiasWeightCompAlg()
EnhancedBiasWeightCompAlg::EnhancedBiasWeightCompAlg |
( |
const std::string & |
name, |
|
|
ISvcLocator * |
svcLoc |
|
) |
| |
◆ calculateEBWeight()
◆ 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.
◆ checkIfTriggeredByRandomChain()
bool EnhancedBiasWeightCompAlg::checkIfTriggeredByRandomChain |
( |
const std::vector< EBChainInfo > & |
EBChain | ) |
const |
|
private |
Check if any of random chains passed.
Definition at line 210 of file EnhancedBiasWeightCompAlg.cxx.
212 return std::find_if (EBChains.begin(), EBChains.end(),
213 [](
const EBChainInfo&
chain) ->
bool { return chain.getIsRandom(); }) != EBChains.end();
◆ 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 EnhancedBiasWeightCompAlg::execute |
( |
const EventContext & |
context | ) |
const |
|
overridevirtual |
Definition at line 80 of file EnhancedBiasWeightCompAlg.cxx.
85 ATH_MSG_DEBUG(
"Number of eb chains that passed int this event: " << EBChains.size());
88 if (EBChains.empty()) {
90 return StatusCode::SUCCESS;
101 if ( context.evt() == 0 ){
103 filename <<
"EnhancedBiasWeights_" << context.eventID().run_number() <<
".xml";
105 ATH_MSG_INFO(
"The output file name is " << m_outputFilename);
109 auto resultPair = std::pair<double, bool>(
result.weight,
result.isUnbiased);
110 auto newPair =
std::find(m_ebWeights.begin(), m_ebWeights.end(), resultPair);
111 if (newPair == m_ebWeights.end()){
112 newPair = m_ebWeights.push_back(resultPair);
113 ATH_MSG_DEBUG(
"New weight value: " <<
result.weight <<
" with id " << (m_ebWeights.size()-1));
115 m_eventToWeight[context.eventID().event_number()] =
std::distance(m_ebWeights.begin(), newPair);
119 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();
◆ fillTotalPrescaleForChains()
StatusCode EnhancedBiasWeightCompAlg::fillTotalPrescaleForChains |
( |
const EventContext & |
context, |
|
|
std::vector< EBChainInfo > & |
EBChains |
|
) |
| const |
|
private |
Retrieve total prescales (L1 * HLT) for chains into map.
Definition at line 123 of file EnhancedBiasWeightCompAlg.cxx.
131 for (EBChainInfo&
chain : EBChains) {
132 if ( not HLTPrescalesSet->prescale(
chain.getId()).enabled ){
133 chain.setTotalPrescale(-1);
137 double HLTPrescale = HLTPrescalesSet->prescale(
chain.getId()).prescale;
138 double L1Prescale = 1.0;
139 if (!
chain.getIsNoPS()) {
141 L1Prescale *= L1PrescalesSet->prescale(
item).prescale;
145 chain.setTotalPrescale(HLTPrescale * L1Prescale);
148 return StatusCode::SUCCESS;
◆ filterPassed()
virtual bool AthReentrantAlgorithm::filterPassed |
( |
const EventContext & |
ctx | ) |
const |
|
inlinevirtualinherited |
◆ getPassedEBChains()
Get list of Enhanced Bias chains that cound have passed.
Definition at line 174 of file EnhancedBiasWeightCompAlg.cxx.
176 std::vector<EBChainInfo> passedEBChains;
181 bool ebChainIsSeeded =
false;
185 if (
m_tdt->isPassedBits(l1Item) & TrigDefs::L1_isPassedBeforePrescale) {
186 ebChainIsSeeded =
true;
190 }
else if (
chainName.find(
"L1RD3") != std::string::npos){
195 if (
m_tdt->isPassedBits(l1Item) & TrigDefs::L1_isPassedBeforePrescale) {
196 ebChainIsSeeded =
true;
201 if (ebChainIsSeeded) {
206 return passedEBChains;
◆ initialize()
StatusCode EnhancedBiasWeightCompAlg::initialize |
( |
| ) |
|
|
overridevirtual |
◆ 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.
◆ 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.
◆ parseItems()
std::vector< std::string > EnhancedBiasWeightCompAlg::parseItems |
( |
const std::string & |
itemStr | ) |
|
|
private |
Parse list of l1 items.
Definition at line 216 of file EnhancedBiasWeightCompAlg.cxx.
218 std::vector<std::string>
items;
220 std::stringstream itemStream (itemStr);
223 while (std::getline(itemStream,
item,
',')) {
◆ renounce()
◆ renounceArray()
◆ setFilterPassed()
virtual void AthReentrantAlgorithm::setFilterPassed |
( |
bool |
state, |
|
|
const EventContext & |
ctx |
|
) |
| const |
|
inlinevirtualinherited |
◆ start()
StatusCode EnhancedBiasWeightCompAlg::start |
( |
| ) |
|
|
overridevirtual |
◆ stop()
StatusCode EnhancedBiasWeightCompAlg::stop |
( |
| ) |
|
|
overridevirtual |
Definition at line 51 of file EnhancedBiasWeightCompAlg.cxx.
54 if (!m_outputFilename.empty()){
57 outputStream <<
"<?xml version=\"1.0\" encoding=\"us-ascii\"?>" << std::endl;
61 for (
size_t i = 0;
i < m_ebWeights.size(); ++
i){
62 outputStream <<
"<weight id=\"" <<
i <<
"\" value=\"" << m_ebWeights[
i].first <<
"\" unbiased=\"" << m_ebWeights[
i].second <<
"\"/>" << std::endl;
68 for (
const auto&
event : m_eventToWeight){
69 outputStream <<
"<e n=\"" <<
event.first <<
"\" w=\"" <<
event.second <<
"\"/>" << std::endl;
77 return StatusCode::SUCCESS;
◆ 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()
◆ ATLAS_THREAD_SAFE [1/3]
std::string m_outputFilename EnhancedBiasWeightCompAlg::ATLAS_THREAD_SAFE |
|
mutableprivate |
◆ ATLAS_THREAD_SAFE [2/3]
tbb::concurrent_vector<std::pair<double, bool> > m_ebWeights EnhancedBiasWeightCompAlg::ATLAS_THREAD_SAFE |
|
mutableprivate |
◆ ATLAS_THREAD_SAFE [3/3]
tbb::concurrent_unordered_map<int, double> m_eventToWeight EnhancedBiasWeightCompAlg::ATLAS_THREAD_SAFE |
|
mutableprivate |
◆ m_chainToHLTSeed
Gaudi::Property<std::map<std::string, std::vector<std::string> > > EnhancedBiasWeightCompAlg::m_chainToHLTSeed {this, "ChainToItemMap", {}, "Map for HLT seeded chains" } |
|
private |
◆ m_detStore
◆ m_EBChainIds
◆ m_EBChainIdToItem
std::map<HLT::Identifier, std::vector<std::string> > EnhancedBiasWeightCompAlg::m_EBChainIdToItem |
|
private |
◆ 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_finalDecisionKey
◆ m_HLTMenuKey
◆ m_HLTPrescaleSetInputKey
◆ m_L1PrescaleSetInputKey
◆ m_tdt
◆ m_varHandleArraysDeclared
◆ m_vhka
The documentation for this class was generated from the following files:
std::vector< EBChainInfo > getPassedEBChains() const
Get list of Enhanced Bias chains that cound have passed.
std::vector< std::string > parseItems(const std::string &itemStr)
Parse list of l1 items.
std::vector< HLT::Identifier > m_EBChainIds
Available EB chains' IDs from HLT Menu.
std::string find(const std::string &s)
return a remapped string
std::map< HLT::Identifier, std::vector< std::string > > m_EBChainIdToItem
L1 items for EB chains from HLT Menu.
SG::ReadHandleKey< TrigConf::HLTMenu > m_HLTMenuKey
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 checkIfTriggeredByRandomChain(const std::vector< EBChainInfo > &EBChain) const
Check if any of random chains passed.
StatusCode fillTotalPrescaleForChains(const EventContext &context, std::vector< EBChainInfo > &EBChains) const
Retrieve total prescales (L1 * HLT) for chains into map.
EBResult calculateEBWeight(const std::vector< EBChainInfo > &EBChains) const
Calculate EB result based on total prescales of chains.
virtual void setOwner(IDataHandleHolder *o)=0
SG::ReadCondHandle< T > makeHandle(const SG::ReadCondHandleKey< T > &key, const EventContext &ctx=Gaudi::Hive::currentContext())
AthReentrantAlgorithm()
Default constructor:
POOL::TEvent event(POOL::TEvent::kClassAccess)
virtual std::vector< Gaudi::DataHandle * > outputHandles() const override
Return this algorithm's output handles.
::StatusCode StatusCode
StatusCode definition for legacy code.
Gaudi::Property< std::map< std::string, std::vector< std::string > > > m_chainToHLTSeed
StoreGateSvc_t m_detStore
Pointer to StoreGate (detector store by default)
StatusCode initialize(bool used=true)
If this object is used as a property, then this should be called during the initialize phase.
SG::ReadCondHandleKey< TrigConf::L1PrescalesSet > m_L1PrescaleSetInputKey
virtual bool isValid() override final
Can the handle be successfully dereferenced?
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
virtual StatusCode sysInitialize() override
Override sysInitialize.
StatusCode initialize(bool used=true)
PublicToolHandle< Trig::TrigDecisionTool > m_tdt
#define ATH_MSG_WARNING(x)
std::string name() const
reports human redable name if it is enabled or, empty string
SG::VarHandleKey & vhKey()
Return a non-const reference to the HandleKey.
SG::ReadCondHandleKey< TrigConf::HLTPrescalesSet > m_HLTPrescaleSetInputKey
float distance(const Amg::Vector3D &p1, const Amg::Vector3D &p2)
calculates the distance between two point in 3D space
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>
SG::ReadHandleKey< TrigCompositeUtils::DecisionContainer > m_finalDecisionKey