ATLAS Offline Software
Loading...
Searching...
No Matches
TrigBphysMonitorAlgorithm Class Reference

#include <TrigBphysMonitorAlgorithm.h>

Inheritance diagram for TrigBphysMonitorAlgorithm:
Collaboration diagram for TrigBphysMonitorAlgorithm:

Public Types

enum class  Environment_t {
  user = 0 , online , tier0 , tier0Raw ,
  tier0ESD , AOD , altprod
}
 Specifies the processing environment. More...
enum class  DataType_t {
  userDefined = 0 , monteCarlo , collisions , cosmics ,
  heavyIonCollisions
}
 Specifies what type of input data is being monitored. More...

Public Member Functions

 TrigBphysMonitorAlgorithm (const std::string &name, ISvcLocator *pSvcLocator)
virtual ~TrigBphysMonitorAlgorithm ()
virtual StatusCode initialize () override
 initialize
virtual StatusCode fillHistograms (const EventContext &ctx) const override
 adds event to the monitoring histograms
virtual StatusCode execute (const EventContext &ctx) const override
 Applies filters and trigger requirements.
void fill (const ToolHandle< GenericMonitoringTool > &groupHandle, std::vector< std::reference_wrapper< Monitored::IMonitoredVariable > > &&variables) const
 Fills a vector of variables to a group by reference.
void fill (const ToolHandle< GenericMonitoringTool > &groupHandle, const std::vector< std::reference_wrapper< Monitored::IMonitoredVariable > > &variables) const
 Fills a vector of variables to a group by reference.
template<typename... T>
void fill (const ToolHandle< GenericMonitoringTool > &groupHandle, T &&... variables) const
 Fills a variadic list of variables to a group by reference.
void fill (const std::string &groupName, std::vector< std::reference_wrapper< Monitored::IMonitoredVariable > > &&variables) const
 Fills a vector of variables to a group by name.
void fill (const std::string &groupName, const std::vector< std::reference_wrapper< Monitored::IMonitoredVariable > > &variables) const
 Fills a vector of variables to a group by name.
template<typename... T>
void fill (const std::string &groupName, T &&... variables) const
 Fills a variadic list of variables to a group by name.
Environment_t environment () const
 Accessor functions for the environment.
Environment_t envStringToEnum (const std::string &str) const
 Convert the environment string from the python configuration to an enum object.
DataType_t dataType () const
 Accessor functions for the data type.
DataType_t dataTypeStringToEnum (const std::string &str) const
 Convert the data type string from the python configuration to an enum object.
const ToolHandle< GenericMonitoringTool > & getGroup (const std::string &name) const
 Get a specific monitoring tool from the tool handle array.
const ToolHandle< Trig::TrigDecisionTool > & getTrigDecisionTool () const
 Get the trigger decision tool member.
bool trigChainsArePassed (const std::vector< std::string > &vTrigNames) const
 Check whether triggers are passed.
SG::ReadHandle< xAOD::EventInfoGetEventInfo (const EventContext &) const
 Return a ReadHandle for an EventInfo object (get run/event numbers, etc.)
virtual float lbAverageInteractionsPerCrossing (const EventContext &ctx=Gaudi::Hive::currentContext()) const
 Calculate the average mu, i.e.
virtual float lbInteractionsPerCrossing (const EventContext &ctx=Gaudi::Hive::currentContext()) const
 Calculate instantaneous number of interactions, i.e.
virtual float lbAverageLuminosity (const EventContext &ctx=Gaudi::Hive::currentContext()) const
 Calculate average luminosity (in ub-1 s-1 => 10^30 cm-2 s-1).
virtual float lbLuminosityPerBCID (const EventContext &ctx=Gaudi::Hive::currentContext()) const
 Calculate the instantaneous luminosity per bunch crossing.
virtual double lbDuration (const EventContext &ctx=Gaudi::Hive::currentContext()) const
 Calculate the duration of the luminosity block (in seconds)
virtual float lbAverageLivefraction (const EventContext &ctx=Gaudi::Hive::currentContext()) const
 Calculate the average luminosity livefraction.
virtual float livefractionPerBCID (const EventContext &ctx=Gaudi::Hive::currentContext()) const
 Calculate the live fraction per bunch crossing ID.
virtual double lbLumiWeight (const EventContext &ctx=Gaudi::Hive::currentContext()) const
 Calculate the average integrated luminosity multiplied by the live fraction.
virtual StatusCode parseList (const std::string &line, std::vector< std::string > &result) const
 Parse a string into a vector.
virtual StatusCode sysInitialize () override
 Override sysInitialize.
virtual bool isClonable () const override
 Specify if the algorithm is clonable.
virtual unsigned int cardinality () const override
 Cardinality (Maximum number of clones that can exist) special value 0 means that algorithm is reentrant.
virtual StatusCode sysExecute (const EventContext &ctx) override
 Execute an algorithm.
virtual const DataObjIDColl & extraOutputDeps () const override
 Return the list of extra output dependencies.
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.
const ServiceHandle< StoreGateSvc > & detStore () const
 The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.
virtual StatusCode sysStart () override
 Handle START transition.
virtual std::vector< Gaudi::DataHandle * > inputHandles () const override
 Return this algorithm's input handles.
virtual std::vector< Gaudi::DataHandle * > outputHandles () const override
 Return this algorithm's output handles.
Gaudi::Details::PropertyBase & declareProperty (Gaudi::Property< T, V, H > &t)
void updateVHKA (Gaudi::Details::PropertyBase &)
MsgStream & msg () const
bool msgLvl (const MSG::Level lvl) const

Public Attributes

str nightly = '/cvmfs/atlas-nightlies.cern.ch/repo/data/data-art/CommonInputs/'
str file = 'data16_13TeV.00311321.physics_Main.recon.AOD.r9264/AOD.11038520._000001.pool.root.1'
 flags = initConfigFlags()
 Files
 isMC
 HISTFileName
 cfg = MainServicesCfg(flags)
 trigBphysMonitorAcc = TrigBphysMonConfig(flags)
 withDetails

Protected Member Functions

void renounceArray (SG::VarHandleKeyArray &handlesArray)
 remove all handles from I/O resolution
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.

Protected Attributes

ToolHandleArray< GenericMonitoringToolm_tools {this,"GMTools",{}}
 Array of Generic Monitoring Tools.
PublicToolHandle< Trig::TrigDecisionToolm_trigDecTool
 Tool to tell whether a specific trigger is passed.
ToolHandleArray< IDQFilterToolm_DQFilterTools {this,"FilterTools",{}}
 Array of Data Quality filter tools.
SG::ReadCondHandleKey< LuminosityCondDatam_lumiDataKey {this,"LuminosityCondDataKey","LuminosityCondData","SG Key of LuminosityCondData object"}
SG::ReadCondHandleKey< LBDurationCondDatam_lbDurationDataKey {this,"LBDurationCondDataKey","LBDurationCondData","SG Key of LBDurationCondData object"}
SG::ReadCondHandleKey< TrigLiveFractionCondDatam_trigLiveFractionDataKey {this,"TrigLiveFractionCondDataKey","TrigLiveFractionCondData", "SG Key of TrigLiveFractionCondData object"}
AthMonitorAlgorithm::Environment_t m_environment
 Instance of the Environment_t enum.
AthMonitorAlgorithm::DataType_t m_dataType
 Instance of the DataType_t enum.
Gaudi::Property< std::string > m_environmentStr {this,"Environment","user"}
 Environment string pulled from the job option and converted to enum.
Gaudi::Property< std::string > m_dataTypeStr {this,"DataType","userDefined"}
 DataType string pulled from the job option and converted to enum.
Gaudi::Property< std::string > m_triggerChainString {this,"TriggerChain",""}
 Trigger chain string pulled from the job option and parsed into a vector.
std::vector< std::string > m_vTrigChainNames
 Vector of trigger chain names parsed from trigger chain string.
Gaudi::Property< std::string > m_fileKey {this,"FileKey",""}
 Internal Athena name for file.
Gaudi::Property< bool > m_useLumi {this,"EnableLumi",false}
 Allows use of various luminosity functions.
Gaudi::Property< float > m_defaultLBDuration {this,"DefaultLBDuration",60.}
 Default duration of one lumi block.
Gaudi::Property< int > m_detailLevel {this,"DetailLevel",0}
 Sets the level of detail used in the monitoring.
SG::ReadHandleKey< xAOD::EventInfom_EventInfoKey {this,"EventInfoKey","EventInfo"}
 Key for retrieving EventInfo from StoreGate.

Private Types

typedef std::vector< std::reference_wrapper< Monitored::IMonitoredVariable > > MonVarVec_t
typedef ServiceHandle< StoreGateSvcStoreGateSvc_t

Private Member Functions

StatusCode fillContainers (const EventContext &ctx) const
 For lxy etc.
StatusCode fillContainerHists (const EventContext &ctx, const SG::ReadHandleKey< xAOD::TrigBphysContainer > &trigBphysContainerKey) const
StatusCode fillChains (const EventContext &ctx) const
StatusCode fillDimuonChainHists (const EventContext &ctx, const std::string &chainName) const
StatusCode fillBmumuxChainHists (const EventContext &ctx, const std::string &chainName) const
StatusCode fillDielectronChainHists (const EventContext &ctx, const std::string &chainName) const
StatusCode fillBhhChainHists (const EventContext &ctx, const std::string &chainName) const
StatusCode fillChainGenericHists (const EventContext &, const ToolHandle< GenericMonitoringTool > &currentMonGroup, const std::string &chainName) const
StatusCode fillBphysObjectHists (const ToolHandle< GenericMonitoringTool > &currentMonGroup, const ElementLink< xAOD::TrigBphysContainer > &bphysLink, const std::string &objStr) const
StatusCode fillTrigLeptonHists (const ToolHandle< GenericMonitoringTool > &currentMonGroup, const ElementLink< xAOD::TrigBphysContainer > &bphysLink, const std::string &name_prefix) const
StatusCode fillTrigBmumuxTrkHists (const ToolHandle< GenericMonitoringTool > &currentMonGroup, const ElementLink< xAOD::TrigBphysContainer > &bphysLink, UInt_t tracksStartFrom=2) const
StatusCode fillDiTracksHists (const ToolHandle< GenericMonitoringTool > &currentMonGroup, const std::vector< ElementLink< xAOD::TrackParticleContainer > > &tpLinkVector, const std::string &name_prefix) const
StatusCode fillTracksHists (const ToolHandle< GenericMonitoringTool > &currentMonGroup, const std::vector< ElementLink< xAOD::TrackParticleContainer > > &tpLinkVector, const std::string &prefix, bool separateHists=false, UInt_t offsetIndex=0) const
StatusCode fillTrkHists (const ToolHandle< GenericMonitoringTool > &currentMonGroup, const xAOD::TrackParticle *trk, const std::string &name_prefix) const
StatusCode fillOfflineDimuons (const EventContext &ctx, const std::vector< std::unique_ptr< xAOD::Vertex > > &dimuonContainer) const
StatusCode fillOfflineDimuonHists (const EventContext &, const std::string &dimuonMonGroupName, const std::vector< std::unique_ptr< xAOD::Vertex > > &dimuonContainer) const
StatusCode fillVertexHists (const ToolHandle< GenericMonitoringTool > &currentMonGroup, const xAOD::Vertex *vertex, const std::string &objStr) const
bool isChainPassed (const std::string &chain) const
StatusCode buildDimuons (const EventContext &ctx, std::vector< std::unique_ptr< xAOD::Vertex > > &vxContainer) const
std::unique_ptr< xAOD::VertexdimuonFit (const xAOD::TrackParticle *mu1, const xAOD::TrackParticle *mu2) const
double dimuonMass (const xAOD::TrackParticle *mu1, const xAOD::TrackParticle *mu2) const
std::vector< const xAOD::Vertex * > GetGoodPVs (const xAOD::VertexContainer *pvContainer) const
const xAOD::VertexgetPvForDimuon_lowestA0 (const xAOD::Vertex *vtx, const std::vector< const xAOD::Vertex * > &PVs) const
bool matchDimuon (const xAOD::Vertex *dimuonVertex, const std::string &chainName) const
bool matchDimuon (const xAOD::Vertex *dimuonVertex, const ElementLink< xAOD::TrigBphysContainer > &bphysLink) const
bool matchTrackParticles (const xAOD::TrackParticle *trk1, const xAOD::TrackParticle *trk2) const
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
 specialization for handling Gaudi::Property<SG::VarHandleKey>

Private Attributes

SG::ReadHandleKeyArray< xAOD::TrigBphysContainerm_TrigBphysContainerKeys {this, "ContainerNames", {}}
SG::ReadHandleKey< xAOD::MuonContainerm_offlineMuonCollectionKey {this, "offlineMuonCollectionKey", "Muons"}
SG::ReadHandleKey< xAOD::TrackParticleContainerm_offlineIDTrackCollectionKey {this, "offlineIDTrackCollectionKey", "InDetTrackParticles"}
SG::ReadHandleKey< xAOD::VertexContainerm_offlinePvCollectionKey {this,"offlinePvCollectionKey","PrimaryVertices"}
Gaudi::Property< std::vector< std::string > > m_ChainNames_MuMu {this, "ChainNames_MuMu", {}}
Gaudi::Property< std::vector< std::string > > m_ChainNames_MuMuX {this, "ChainNames_MuMuX", {}}
Gaudi::Property< std::vector< std::string > > m_ChainNames_ElEl {this, "ChainNames_ElEl", {}}
Gaudi::Property< std::vector< std::string > > m_ChainNames_Bhh {this, "ChainNames_Bhh", {}}
Gaudi::Property< bool > m_requireExplicitESDecision {this, "requireExplicitESDecision", false}
Gaudi::Property< double > m_dimuMassLower_prefit {this, "dimuMassLower_prefit", 0.}
Gaudi::Property< double > m_dimuMassUpper_prefit {this, "dimuMassUpper_prefit", 16000.}
Gaudi::Property< double > m_dimuMassLower_postfit {this, "dimuMassLower_postfit", 0.}
Gaudi::Property< double > m_dimuMassUpper_postfit {this, "dimuMassUpper_postfit", 15000.}
Gaudi::Property< double > m_dimuChi2Cut {this, "dimuChi2Cut", 20.}
Gaudi::Property< double > m_deltaRMatchingThreshold {this, "deltaRMatchingThreshold", 0.05}
ToolHandle< InDet::VertexPointEstimatorm_vertexPointEstimator {this, "VertexPointEstimator", "", "tool to find starting point for the vertex fitter"}
ToolHandle< Trk::TrkVKalVrtFitterm_vertexFitter {this, "VertexFitter", "", "VKalVrtFitter tool to fit tracks into the common vertex"}
ToolHandle< Trk::V0Toolsm_v0Tools {this, "V0Tools", "", "V0 tools to calculate things like Lxy"}
std::string m_name
std::unordered_map< std::string, size_t > m_toolLookupMap
const ToolHandle< GenericMonitoringToolm_dummy
Gaudi::Property< bool > m_enforceExpressTriggers
DataObjIDColl m_extendedExtraObjects
 Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.
StoreGateSvc_t m_evtStore
 Pointer to StoreGate (event store by default)
StoreGateSvc_t m_detStore
 Pointer to StoreGate (detector store by default)
std::vector< SG::VarHandleKeyArray * > m_vhka
bool m_varHandleArraysDeclared

Detailed Description

Definition at line 20 of file TrigBphysMonitorAlgorithm.h.

Member Typedef Documentation

◆ MonVarVec_t

typedef std::vector<std::reference_wrapper<Monitored::IMonitoredVariable> > AthMonitorAlgorithm::MonVarVec_t
privateinherited

Definition at line 370 of file AthMonitorAlgorithm.h.

◆ StoreGateSvc_t

typedef ServiceHandle<StoreGateSvc> AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::StoreGateSvc_t
privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Member Enumeration Documentation

◆ DataType_t

enum class AthMonitorAlgorithm::DataType_t
stronginherited

Specifies what type of input data is being monitored.

An enumeration of the different types of data the monitoring application may be running over. This can be used to select which histograms to produce, e.g. to prevent the production of colliding-beam histograms when running on cosmic-ray data. Strings of the same names may be given as jobOptions.

Enumerator
userDefined 
monteCarlo 
collisions 
cosmics 
heavyIonCollisions 

Definition at line 194 of file AthMonitorAlgorithm.h.

194 {
195 userDefined = 0,
196 monteCarlo,
197 collisions,
198 cosmics,
199 heavyIonCollisions,
200 };
AthConfigFlags cosmics(AthConfigFlags flags, str instanceName, str recoMode)

◆ Environment_t

enum class AthMonitorAlgorithm::Environment_t
stronginherited

Specifies the processing environment.

The running environment may be used to select which histograms are produced, and where they are located in the output. For example, the output paths of the histograms are different for the "user", "online" and the various offline flags. Strings of the same names may be given as jobOptions.

Enumerator
user 
online 
tier0 
tier0Raw 
tier0ESD 
AOD 
altprod 

Definition at line 175 of file AthMonitorAlgorithm.h.

175 {
176 user = 0,
177 online,
178 tier0,
179 tier0Raw,
180 tier0ESD,
181 AOD,
182 altprod,
183 };

Constructor & Destructor Documentation

◆ TrigBphysMonitorAlgorithm()

TrigBphysMonitorAlgorithm::TrigBphysMonitorAlgorithm ( const std::string & name,
ISvcLocator * pSvcLocator )

Definition at line 16 of file TrigBphysMonitorAlgorithm.cxx.

17 : AthMonitorAlgorithm(name,pSvcLocator)
18{}
Base class for Athena Monitoring Algorithms.

◆ ~TrigBphysMonitorAlgorithm()

TrigBphysMonitorAlgorithm::~TrigBphysMonitorAlgorithm ( )
virtual

Definition at line 21 of file TrigBphysMonitorAlgorithm.cxx.

21{}

Member Function Documentation

◆ buildDimuons()

StatusCode TrigBphysMonitorAlgorithm::buildDimuons ( const EventContext & ctx,
std::vector< std::unique_ptr< xAOD::Vertex > > & vxContainer ) const
private

Definition at line 481 of file TrigBphysMonitorAlgorithm.cxx.

481 {
482 ATH_MSG_DEBUG( "TrigBphysMonitorAlgorithm::buildDimuons" );
483
484 // Get the muons
485 SG::ReadHandle<xAOD::MuonContainer> muonHandle(m_offlineMuonCollectionKey,ctx);
486 ATH_CHECK(muonHandle.isValid());
487 const xAOD::MuonContainer* muonContainer = muonHandle.cptr();
488 ATH_MSG_DEBUG("Muon container size " << muonContainer->size());
489
490 // Muon preselection
491 std::vector<const xAOD::Muon*> selectedMuons;
492 for (const auto mu : *muonContainer) {
493 if ( mu == nullptr ) continue;
494 if ( mu->muonType() != xAOD::Muon::Combined ) continue; // require combined muons
495 if (!mu->inDetTrackParticleLink()) continue; // No muons without ID tracks
496 if (!mu->inDetTrackParticleLink().isValid()) continue; // No muons without ID tracks
497 selectedMuons.push_back(mu);
498 }
499 if(selectedMuons.size() < 2) {
500 ATH_MSG_DEBUG("Only " << selectedMuons.size() << "muons pass preselection");
501 return StatusCode::SUCCESS;
502 }
503 std::sort(selectedMuons.begin(), selectedMuons.end(), [](const auto mu1, const auto mu2){ return mu1->pt() > mu2->pt(); });
504
505 // Build dimuons
506 for(auto outerItr=selectedMuons.begin(); outerItr<selectedMuons.end(); ++outerItr){
507 for(auto innerItr=(outerItr+1); innerItr!=selectedMuons.end(); ++innerItr){
508 const auto muon1 = *outerItr;
509 const auto muon2 = *innerItr;
510 const auto trackParticle1 = muon1->trackParticle( xAOD::Muon::InnerDetectorTrackParticle );
511 const auto trackParticle2 = muon2->trackParticle( xAOD::Muon::InnerDetectorTrackParticle );
512 // Charge selection
513 if(trackParticle1->qOverP() * trackParticle2->qOverP() > 0.)
514 continue;
515 // Mass selection
516 double dimu_momentum_prefit = dimuonMass(trackParticle1, trackParticle2);
517 if( !(dimu_momentum_prefit > m_dimuMassLower_prefit && dimu_momentum_prefit < m_dimuMassUpper_prefit) )
518 continue;
519 // Fit
520 std::unique_ptr<xAOD::Vertex> dimuon = dimuonFit(trackParticle1, trackParticle2);
521 if(!dimuon) continue;
522 if(dimuon->chiSquared() > m_dimuChi2Cut) continue;
523 vxContainer.push_back(std::move(dimuon));
524 }
525 }
526 ATH_MSG_DEBUG("Found " << vxContainer.size() << " good dimuons");
527
528 // Augment the dimuon vertices
529 if( !vxContainer.empty() ) {
530 SG::ReadHandle<xAOD::VertexContainer> pvHandle(m_offlinePvCollectionKey,ctx);
531 ATH_CHECK(pvHandle.isValid());
532 const xAOD::VertexContainer* pvContainer = pvHandle.cptr();
533 std::vector<const xAOD::Vertex*> goodPVs = GetGoodPVs(pvContainer);
534 ATH_MSG_DEBUG("Found " << goodPVs.size() << " good PVs");
535 for(auto& dimuon : vxContainer) {
536 xAOD::BPhysHelper jpsiHelper(dimuon.get());
537 jpsiHelper.setRefTrks();
538 const xAOD::Vertex* lowestA0Pv = getPvForDimuon_lowestA0(dimuon.get(), goodPVs);
539 if(lowestA0Pv) {
540 jpsiHelper.setLxy ( m_v0Tools->lxy ( dimuon.get(),lowestA0Pv ), xAOD::BPhysHelper::PV_MIN_A0 );
541 jpsiHelper.setLxyErr( m_v0Tools->lxyError( dimuon.get(),lowestA0Pv ), xAOD::BPhysHelper::PV_MIN_A0 );
542 ATH_MSG_VERBOSE("Lxy = " << m_v0Tools->lxy ( dimuon.get(),lowestA0Pv ));
543 ATH_MSG_VERBOSE("LxyErr = " << m_v0Tools->lxyError( dimuon.get(),lowestA0Pv ));
544 }
545 }
546 }
547
548 return StatusCode::SUCCESS;
549}
#define ATH_CHECK
Evaluate an expression and check for errors.
#define ATH_MSG_VERBOSE(x)
#define ATH_MSG_DEBUG(x)
xAOD::MuonContainer * muonContainer
ToolHandle< Trk::V0Tools > m_v0Tools
std::vector< const xAOD::Vertex * > GetGoodPVs(const xAOD::VertexContainer *pvContainer) const
const xAOD::Vertex * getPvForDimuon_lowestA0(const xAOD::Vertex *vtx, const std::vector< const xAOD::Vertex * > &PVs) const
std::unique_ptr< xAOD::Vertex > dimuonFit(const xAOD::TrackParticle *mu1, const xAOD::TrackParticle *mu2) const
double dimuonMass(const xAOD::TrackParticle *mu1, const xAOD::TrackParticle *mu2) const
void sort(typename DataModel_detail::iterator< DVL > beg, typename DataModel_detail::iterator< DVL > end)
Specialization of sort for DataVector/List.
VertexContainer_v1 VertexContainer
Definition of the current "Vertex container version".
Vertex_v1 Vertex
Define the latest version of the vertex class.
MuonContainer_v1 MuonContainer
Definition of the current "Muon container version".

◆ cardinality()

unsigned int AthCommonReentrantAlgorithm< Gaudi::Algorithm >::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.

Definition at line 75 of file AthCommonReentrantAlgorithm.cxx.

64{
65 return 0;
66}

◆ dataType()

DataType_t AthMonitorAlgorithm::dataType ( ) const
inlineinherited

Accessor functions for the data type.

Returns
the current value of the class's DataType_t instance.

Definition at line 224 of file AthMonitorAlgorithm.h.

224{ return m_dataType; }
AthMonitorAlgorithm::DataType_t m_dataType
Instance of the DataType_t enum.

◆ dataTypeStringToEnum()

AthMonitorAlgorithm::DataType_t AthMonitorAlgorithm::dataTypeStringToEnum ( const std::string & str) const
inherited

Convert the data type string from the python configuration to an enum object.

Returns
a value in the DataType_t enumeration which matches the input string.

Definition at line 144 of file AthMonitorAlgorithm.cxx.

144 {
145 // convert the string to all lowercase
146 std::string lowerCaseStr = str;
147 std::transform(lowerCaseStr.begin(), lowerCaseStr.end(), lowerCaseStr.begin(), ::tolower);
148
149 // check if it matches one of the enum choices
150 if( lowerCaseStr == "userdefined" ) {
152 } else if( lowerCaseStr == "montecarlo" ) {
154 } else if( lowerCaseStr == "collisions" ) {
156 } else if( lowerCaseStr == "cosmics" ) {
157 return DataType_t::cosmics;
158 } else if( lowerCaseStr == "heavyioncollisions" ) {
160 } else { // otherwise, warn the user and return "userDefined"
161 ATH_MSG_WARNING("AthMonitorAlgorithm::dataTypeStringToEnum(): Unknown data type "
162 <<str<<", returning userDefined.");
164 }
165}
#define ATH_MSG_WARNING(x)
void tolower(std::string &s)

◆ declareGaudiProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareGaudiProperty ( Gaudi::Property< T, V, H > & hndl,
const SG::VarHandleKeyType &  )
inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleKey>

Definition at line 156 of file AthCommonDataStore.h.

158 {
160 hndl.value(),
161 hndl.documentation());
162
163 }
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T, V, H > &t)

◆ declareProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareProperty ( Gaudi::Property< T, V, H > & t)
inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

145 {
146 typedef typename SG::HandleClassifier<T>::type htype;
148 }
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>

◆ detStore()

const ServiceHandle< StoreGateSvc > & AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::detStore ( ) const
inlineinherited

The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 95 of file AthCommonDataStore.h.

◆ dimuonFit()

std::unique_ptr< xAOD::Vertex > TrigBphysMonitorAlgorithm::dimuonFit ( const xAOD::TrackParticle * mu1,
const xAOD::TrackParticle * mu2 ) const
private

Definition at line 551 of file TrigBphysMonitorAlgorithm.cxx.

551 {
552
553 const Trk::Perigee& mu1Perigee = mu1->perigeeParameters();
554 const Trk::Perigee& mu2Perigee = mu2->perigeeParameters();
555 int sflag = 0; int errorcode = 0;
556 Amg::Vector3D startingPoint = m_vertexPointEstimator->getCirclesIntersectionPoint(&mu1Perigee,&mu2Perigee,sflag,errorcode);
557 if (errorcode != 0) {startingPoint(0) = 0.0; startingPoint(1) = 0.0; startingPoint(2) = 0.0;}
558 const std::vector<const xAOD::TrackParticle*> trackPair = {mu1, mu2};
559 std::unique_ptr<xAOD::Vertex> myVxCandidate(m_vertexFitter->fit(trackPair, startingPoint));
560
561 return myVxCandidate;
562}
ToolHandle< Trk::TrkVKalVrtFitter > m_vertexFitter
const Trk::Perigee & perigeeParameters() const
Returns the Trk::MeasuredPerigee track parameters.
Eigen::Matrix< double, 3, 1 > Vector3D
ParametersT< TrackParametersDim, Charged, PerigeeSurface > Perigee

◆ dimuonMass()

double TrigBphysMonitorAlgorithm::dimuonMass ( const xAOD::TrackParticle * mu1,
const xAOD::TrackParticle * mu2 ) const
private

Definition at line 615 of file TrigBphysMonitorAlgorithm.cxx.

615 {
616 using namespace TrigVtx;
621 xAOD::TrackParticle::GenVecFourMom_t dimu_mom = mom1 + mom2;
622 return dimu_mom.M();
623}
ROOT::Math::LorentzVector< ROOT::Math::PxPyPzM4D< double > > GenVecFourMom_t
Base 4 Momentum type for TrackParticle.
GenVecFourMom_t genvecP4() const
The full 4-momentum of the particle : GenVector form.

◆ environment()

Environment_t AthMonitorAlgorithm::environment ( ) const
inlineinherited

Accessor functions for the environment.

Returns
the current value of the class's Environment_t instance.

Definition at line 208 of file AthMonitorAlgorithm.h.

208{ return m_environment; }
AthMonitorAlgorithm::Environment_t m_environment
Instance of the Environment_t enum.

◆ envStringToEnum()

AthMonitorAlgorithm::Environment_t AthMonitorAlgorithm::envStringToEnum ( const std::string & str) const
inherited

Convert the environment string from the python configuration to an enum object.

Returns
a value in the Environment_t enumeration which matches the input string.

Definition at line 116 of file AthMonitorAlgorithm.cxx.

116 {
117 // convert the string to all lowercase
118 std::string lowerCaseStr = str;
119 std::transform(lowerCaseStr.begin(), lowerCaseStr.end(), lowerCaseStr.begin(), ::tolower);
120
121 // check if it matches one of the enum choices
122 if( lowerCaseStr == "user" ) {
123 return Environment_t::user;
124 } else if( lowerCaseStr == "online" ) {
126 } else if( lowerCaseStr == "tier0" ) {
128 } else if( lowerCaseStr == "tier0raw" ) {
130 } else if( lowerCaseStr == "tier0esd" ) {
132 } else if( lowerCaseStr == "aod" ) {
133 return Environment_t::AOD;
134 } else if( lowerCaseStr == "altprod" ) {
136 } else { // otherwise, warn the user and return "user"
137 ATH_MSG_WARNING("AthMonitorAlgorithm::envStringToEnum(): Unknown environment "
138 <<str<<", returning user.");
139 return Environment_t::user;
140 }
141}

◆ evtStore()

ServiceHandle< StoreGateSvc > & AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::evtStore ( )
inlineinherited

The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 85 of file AthCommonDataStore.h.

◆ execute()

StatusCode AthMonitorAlgorithm::execute ( const EventContext & ctx) const
overridevirtualinherited

Applies filters and trigger requirements.

Then, calls fillHistograms().

Parameters
ctxevent context for reentrant Athena call
Returns
StatusCode

Definition at line 77 of file AthMonitorAlgorithm.cxx.

77 {
78
79 // Checks that all of the DQ filters are passed. If any one of the filters
80 // fails, return SUCCESS code and do not fill the histograms with the event.
81 for ( const auto& filterItr : m_DQFilterTools ) {
82 if (!filterItr->accept()) {
83 ATH_MSG_DEBUG("Event rejected due to filter tool.");
84 return StatusCode::SUCCESS;
85 }
86 }
87
88 // Trigger: If there is a decision tool and the chains fail, skip the event.
90 ATH_MSG_DEBUG("Event rejected due to trigger filter.");
91 return StatusCode::SUCCESS;
92 }
93
94 ATH_MSG_DEBUG("Event accepted!");
95 return fillHistograms(ctx);
96}
virtual StatusCode fillHistograms(const EventContext &ctx) const =0
adds event to the monitoring histograms
bool trigChainsArePassed(const std::vector< std::string > &vTrigNames) const
Check whether triggers are passed.
std::vector< std::string > m_vTrigChainNames
Vector of trigger chain names parsed from trigger chain string.
PublicToolHandle< Trig::TrigDecisionTool > m_trigDecTool
Tool to tell whether a specific trigger is passed.
ToolHandleArray< IDQFilterTool > m_DQFilterTools
Array of Data Quality filter tools.

◆ extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::extraDeps_update_handler ( Gaudi::Details::PropertyBase & ExtraDeps)
protectedinherited

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 & AthCommonReentrantAlgorithm< Gaudi::Algorithm >::extraOutputDeps ( ) const
overridevirtualinherited

Return the list of extra output dependencies.

This list is extended to include symlinks implied by inheritance relations.

Definition at line 94 of file AthCommonReentrantAlgorithm.cxx.

90{
91 // If we didn't find any symlinks to add, just return the collection
92 // from the base class. Otherwise, return the extended collection.
93 if (!m_extendedExtraObjects.empty()) {
95 }
97}
An algorithm that can be simultaneously executed in multiple threads.

◆ fillBhhChainHists()

StatusCode TrigBphysMonitorAlgorithm::fillBhhChainHists ( const EventContext & ctx,
const std::string & chainName ) const
private

Definition at line 226 of file TrigBphysMonitorAlgorithm.cxx.

226 {
227 ATH_MSG_DEBUG("Filling " << chainName << " chain histograms");
228
229 std::string monGroupName = std::string("Chain_")+chainName;
230 auto monGroup = getGroup(monGroupName);
231
232 if (fillChainGenericHists(ctx, monGroup, chainName).isFailure()) {
233 ATH_MSG_ERROR(Form("Problems filling generic histograms for %s chain",chainName.c_str()));
234 }
235
236 std::vector< TrigCompositeUtils::LinkInfo<xAOD::TrigBphysContainer> > chainFeatureContainer = getTrigDecisionTool()->features<xAOD::TrigBphysContainer>(chainName, TrigDefs::Physics);
237 for (const auto& featLinkInfo: chainFeatureContainer){
238 ATH_CHECK(featLinkInfo.isValid());
239 const auto featLink = featLinkInfo.link;
240 if (fillBphysObjectHists(monGroup, featLink, "B").isFailure()) {
241 ATH_MSG_ERROR(Form("Problems filling bphys object histograms for %s chain",chainName.c_str()));
242 }
243 if (fillTrigBmumuxTrkHists(monGroup, featLink, 0).isFailure()) { // same logic as for Bmumux, but start from 1st track - there are no muons
244 ATH_MSG_ERROR(Form("Problems filling muon histograms for %s chain",chainName.c_str()));
245 }
246 }
247
248 return StatusCode::SUCCESS;
249}
#define ATH_MSG_ERROR(x)
StatusCode fillTrigBmumuxTrkHists(const ToolHandle< GenericMonitoringTool > &currentMonGroup, const ElementLink< xAOD::TrigBphysContainer > &bphysLink, UInt_t tracksStartFrom=2) const
TrigBphysContainer_v1 TrigBphysContainer

◆ fillBmumuxChainHists()

StatusCode TrigBphysMonitorAlgorithm::fillBmumuxChainHists ( const EventContext & ctx,
const std::string & chainName ) const
private

Definition at line 168 of file TrigBphysMonitorAlgorithm.cxx.

168 {
169 ATH_MSG_DEBUG("Filling " << chainName << " chain histograms");
170
171 std::string monGroupName = std::string("Chain_")+chainName;
172 auto monGroup = getGroup(monGroupName);
173
174 if (fillChainGenericHists(ctx, monGroup, chainName).isFailure()) {
175 ATH_MSG_ERROR(Form("Problems filling generic histograms for %s chain",chainName.c_str()));
176 }
177
178 std::vector< TrigCompositeUtils::LinkInfo<xAOD::TrigBphysContainer> > chainFeatureContainer = getTrigDecisionTool()->features<xAOD::TrigBphysContainer>(chainName, TrigDefs::Physics);
179 for (const auto& featLinkInfo: chainFeatureContainer){
180 ATH_CHECK(featLinkInfo.isValid());
181 const auto featLink = featLinkInfo.link;
182 if (fillBphysObjectHists(monGroup, featLink, "B").isFailure()) {
183 ATH_MSG_ERROR(Form("Problems filling bphys object histograms for %s chain",chainName.c_str()));
184 }
185 if (fillTrigBmumuxTrkHists(monGroup, featLink).isFailure()) {
186 ATH_MSG_ERROR(Form("Problems filling muon histograms for %s chain",chainName.c_str()));
187 }
188
189 auto dimuonLink = (*featLink)->lowerChainLink();
190 ATH_CHECK(dimuonLink.isValid());
191 if (fillBphysObjectHists(monGroup, dimuonLink, "dimu").isFailure()) {
192 ATH_MSG_ERROR(Form("Problems filling dimuon object histograms for %s chain",chainName.c_str()));
193 }
194 if (fillTrigLeptonHists(monGroup, dimuonLink, "mu").isFailure()) {
195 ATH_MSG_ERROR(Form("Problems filling muon histograms for %s chain",chainName.c_str()));
196 }
197 }
198
199 return StatusCode::SUCCESS;
200}
StatusCode fillBphysObjectHists(const ToolHandle< GenericMonitoringTool > &currentMonGroup, const ElementLink< xAOD::TrigBphysContainer > &bphysLink, const std::string &objStr) const
StatusCode fillTrigLeptonHists(const ToolHandle< GenericMonitoringTool > &currentMonGroup, const ElementLink< xAOD::TrigBphysContainer > &bphysLink, const std::string &name_prefix) const

◆ fillBphysObjectHists()

StatusCode TrigBphysMonitorAlgorithm::fillBphysObjectHists ( const ToolHandle< GenericMonitoringTool > & currentMonGroup,
const ElementLink< xAOD::TrigBphysContainer > & bphysLink,
const std::string & objStr ) const
private

Definition at line 266 of file TrigBphysMonitorAlgorithm.cxx.

266 {
267
268 auto dimu_mass = Monitored::Scalar<float>(objStr+"_mass",-999.);
269 auto dimu_fitmass = Monitored::Scalar<float>(objStr+"_fitmass",-999.);
270 auto dimu_pt = Monitored::Scalar<float>(objStr+"_pt",-999.);
271 auto dimu_y = Monitored::Scalar<float>(objStr+"_y",-999.);
272 auto dimu_chi2 = Monitored::Scalar<float>(objStr+"_chi2",-999.);
273
274 dimu_mass = (*bphysLink)->mass() / GeV;
275 dimu_fitmass = (*bphysLink)->fitmass() / GeV;
276 dimu_pt = (*bphysLink)->pt() / GeV;
277 dimu_chi2 = (*bphysLink)->fitchi2();
278 dimu_y = (*bphysLink)->eta();
279
280 fill(currentMonGroup, dimu_mass, dimu_fitmass, dimu_pt, dimu_y, dimu_chi2);
281
282 return StatusCode::SUCCESS;
283}
Declare a monitored scalar variable.
void fill(H5::Group &out_file, size_t iterations)

◆ fillChainGenericHists()

StatusCode TrigBphysMonitorAlgorithm::fillChainGenericHists ( const EventContext & ,
const ToolHandle< GenericMonitoringTool > & currentMonGroup,
const std::string & chainName ) const
private

Definition at line 253 of file TrigBphysMonitorAlgorithm.cxx.

253 {
254
255 std::vector< TrigCompositeUtils::LinkInfo<xAOD::TrigBphysContainer> > chainFeatureContainer = getTrigDecisionTool()->features<xAOD::TrigBphysContainer>(chainName, TrigDefs::Physics);
256 auto ncandidates = Monitored::Scalar<int>("ncandidates",-999);
257 ncandidates = chainFeatureContainer.size();
258
259 fill(currentMonGroup, ncandidates);
260
261 return StatusCode::SUCCESS;
262}

◆ fillChains()

StatusCode TrigBphysMonitorAlgorithm::fillChains ( const EventContext & ctx) const
private

Definition at line 94 of file TrigBphysMonitorAlgorithm.cxx.

94 {
95
96 for(const auto& monitoredChain : m_ChainNames_MuMu) {
97 ATH_MSG_DEBUG("Process chain " << monitoredChain);
98 if( !isChainPassed(monitoredChain) ) {
99 ATH_MSG_DEBUG("Chain " << monitoredChain << " is not passed");
100 continue;
101 }
102 if(fillDimuonChainHists(ctx, monitoredChain).isFailure()) {
103 ATH_MSG_ERROR(Form("Problems filling %s chain histograms", monitoredChain.c_str()));
104 }
105 }
106
107 for(const auto& monitoredChain : m_ChainNames_MuMuX) {
108 ATH_MSG_DEBUG("Process chain " << monitoredChain);
109 if( !isChainPassed(monitoredChain) ) {
110 ATH_MSG_DEBUG("Chain " << monitoredChain << " is not passed");
111 continue;
112 }
113 if(fillBmumuxChainHists(ctx, monitoredChain).isFailure()) {
114 ATH_MSG_ERROR(Form("Problems filling %s chain histograms", monitoredChain.c_str()));
115 }
116 }
117
118 for(const auto& monitoredChain : m_ChainNames_ElEl) {
119 ATH_MSG_DEBUG("Process chain " << monitoredChain);
120 if( !isChainPassed(monitoredChain) ) {
121 ATH_MSG_DEBUG("Chain " << monitoredChain << " is not passed");
122 continue;
123 }
124 if(fillDielectronChainHists(ctx, monitoredChain).isFailure()) {
125 ATH_MSG_ERROR(Form("Problems filling %s chain histograms", monitoredChain.c_str()));
126 }
127 }
128
129 for(const auto& monitoredChain : m_ChainNames_Bhh) {
130 ATH_MSG_DEBUG("Process chain " << monitoredChain);
131 if( !isChainPassed(monitoredChain) ) {
132 ATH_MSG_DEBUG("Chain " << monitoredChain << " is not passed");
133 continue;
134 }
135 if(fillBhhChainHists(ctx, monitoredChain).isFailure()) {
136 ATH_MSG_ERROR(Form("Problems filling %s chain histograms", monitoredChain.c_str()));
137 }
138 }
139
140 return StatusCode::SUCCESS;
141}
StatusCode fillBhhChainHists(const EventContext &ctx, const std::string &chainName) const
StatusCode fillBmumuxChainHists(const EventContext &ctx, const std::string &chainName) const
bool isChainPassed(const std::string &chain) const
StatusCode fillDielectronChainHists(const EventContext &ctx, const std::string &chainName) const

◆ fillContainerHists()

StatusCode TrigBphysMonitorAlgorithm::fillContainerHists ( const EventContext & ctx,
const SG::ReadHandleKey< xAOD::TrigBphysContainer > & trigBphysContainerKey ) const
private

Definition at line 72 of file TrigBphysMonitorAlgorithm.cxx.

72 {
73 SG::ReadHandle<xAOD::TrigBphysContainer> trigBphysContainer(trigBphysContainerKey, ctx);
74 if( !trigBphysContainer.isValid() ) {
75 ATH_MSG_WARNING("No valid TrigBphysContainer with tag: " << trigBphysContainerKey);
76 }
77 else {
78 // Fill container-wise histograms
79 std::string monGroupName = std::string("Container_")+trigBphysContainerKey.key();
80 auto monGroup = getGroup(monGroupName);
81
82 auto ncandidates = Monitored::Scalar<int>("ncandidates",-999);
83 ncandidates = trigBphysContainer->size();
84
85 auto bphys_mass = Monitored::Collection("bphys_mass", (*trigBphysContainer), [](const xAOD::TrigBphys* bphys){ return bphys->fitmass() / GeV;});
86
87 fill(monGroup, ncandidates, bphys_mass);
88
89 }
90 return StatusCode::SUCCESS;
91}
const std::string & key() const
Return the StoreGate ID for the referenced object.
float fitmass() const
accessor method: mass from vertex fit
ValuesCollection< T > Collection(std::string name, const T &collection)
Declare a monitored (double-convertible) collection.
TrigBphys_v1 TrigBphys
Definition TrigBphys.h:18

◆ fillContainers()

StatusCode TrigBphysMonitorAlgorithm::fillContainers ( const EventContext & ctx) const
private

For lxy etc.

Definition at line 60 of file TrigBphysMonitorAlgorithm.cxx.

60 {
61
62 for(const SG::ReadHandleKey<xAOD::TrigBphysContainer>& monitoredContainerKey : m_TrigBphysContainerKeys) {
63 ATH_MSG_DEBUG("Process container " << monitoredContainerKey.key().c_str());
64 if(fillContainerHists(ctx, monitoredContainerKey).isFailure()) {
65 ATH_MSG_ERROR(Form("Problems filling %s container histograms", monitoredContainerKey.key().c_str()));
66 }
67 }
68
69 return StatusCode::SUCCESS;
70}
Property holding a SG store/key/clid from which a ReadHandle is made.

◆ fillDielectronChainHists()

StatusCode TrigBphysMonitorAlgorithm::fillDielectronChainHists ( const EventContext & ctx,
const std::string & chainName ) const
private

Definition at line 202 of file TrigBphysMonitorAlgorithm.cxx.

202 {
203 ATH_MSG_DEBUG("Filling " << chainName << " chain histograms");
204
205 std::string monGroupName = std::string("Chain_")+chainName;
206 auto monGroup = getGroup(monGroupName);
207
208 if (fillChainGenericHists(ctx, monGroup, chainName).isFailure()) {
209 ATH_MSG_ERROR(Form("Problems filling generic histograms for %s chain",chainName.c_str()));
210 }
211
212 std::vector< TrigCompositeUtils::LinkInfo<xAOD::TrigBphysContainer> > chainFeatureContainer = getTrigDecisionTool()->features<xAOD::TrigBphysContainer>(chainName, TrigDefs::Physics);
213 for (const auto& featLinkInfo: chainFeatureContainer){
214 ATH_CHECK(featLinkInfo.isValid());
215 const auto featLink = featLinkInfo.link;
216 if (fillBphysObjectHists(monGroup, featLink, "diel").isFailure()) {
217 ATH_MSG_ERROR(Form("Problems filling bphys object histograms for %s chain",chainName.c_str()));
218 }
219 if (fillTrigLeptonHists(monGroup, featLink, "el").isFailure()) {
220 ATH_MSG_ERROR(Form("Problems filling muon histograms for %s chain",chainName.c_str()));
221 }
222 }
223 return StatusCode::SUCCESS;
224}

◆ fillDimuonChainHists()

StatusCode TrigBphysMonitorAlgorithm::fillDimuonChainHists ( const EventContext & ctx,
const std::string & chainName ) const
private

Definition at line 143 of file TrigBphysMonitorAlgorithm.cxx.

143 {
144 ATH_MSG_DEBUG("Filling " << chainName << " chain histograms");
145
146 std::string monGroupName = std::string("Chain_")+chainName;
147 auto monGroup = getGroup(monGroupName);
148
149 if (fillChainGenericHists(ctx, monGroup, chainName).isFailure()) {
150 ATH_MSG_ERROR(Form("Problems filling generic histograms for %s chain",chainName.c_str()));
151 }
152
153 std::vector< TrigCompositeUtils::LinkInfo<xAOD::TrigBphysContainer> > chainFeatureContainer = getTrigDecisionTool()->features<xAOD::TrigBphysContainer>(chainName, TrigDefs::Physics);
154 for (const auto& featLinkInfo: chainFeatureContainer){
155 ATH_CHECK(featLinkInfo.isValid());
156 const auto featLink = featLinkInfo.link;
157 if (fillBphysObjectHists(monGroup, featLink, "dimu").isFailure()) {
158 ATH_MSG_ERROR(Form("Problems filling bphys object histograms for %s chain",chainName.c_str()));
159 }
160 if (fillTrigLeptonHists(monGroup, featLink, "mu").isFailure()) {
161 ATH_MSG_ERROR(Form("Problems filling muon histograms for %s chain",chainName.c_str()));
162 }
163 }
164 return StatusCode::SUCCESS;
165}

◆ fillDiTracksHists()

StatusCode TrigBphysMonitorAlgorithm::fillDiTracksHists ( const ToolHandle< GenericMonitoringTool > & currentMonGroup,
const std::vector< ElementLink< xAOD::TrackParticleContainer > > & tpLinkVector,
const std::string & name_prefix ) const
private

Definition at line 338 of file TrigBphysMonitorAlgorithm.cxx.

340 {
341 // Use first two tracks
342 if (tpLinkVector.size() <2) {
343 ATH_MSG_ERROR("Unexpected number of tracks in a dimuon: " << tpLinkVector.size());
344 return StatusCode::FAILURE;
345 }
346 ATH_CHECK(tpLinkVector.at(0).isValid());
347 ATH_CHECK(tpLinkVector.at(1).isValid());
348 const xAOD::TrackParticle* trk1 = *tpLinkVector.at(0);
349 const xAOD::TrackParticle* trk2 = *tpLinkVector.at(1);
350
351 auto ditrk_dR = Monitored::Scalar<float>(name_prefix+"_dR",-999.);
352 auto ditrk_deta = Monitored::Scalar<float>(name_prefix+"_deta",-999.);
353 auto ditrk_dphi = Monitored::Scalar<float>(name_prefix+"_dphi",-999.);
354
355 ditrk_dR = xAOD::P4Helpers::deltaR(*trk1, *trk2, false); // false for pseudo, not true rapidity
356 ditrk_deta = std::abs(trk1->eta()-trk2->eta());
357 ditrk_dphi = std::abs(xAOD::P4Helpers::deltaPhi(*trk1, *trk2));
358
359 fill(currentMonGroup, ditrk_dR, ditrk_deta, ditrk_dphi);
360
361 return StatusCode::SUCCESS;
362}
virtual double eta() const override final
The pseudorapidity ( ) of the particle.
double deltaPhi(double phiA, double phiB)
delta Phi in range [-pi,pi[
double deltaR(double rapidity1, double phi1, double rapidity2, double phi2)
from bare bare rapidity,phi
TrackParticle_v1 TrackParticle
Reference the current persistent version:

◆ fillHistograms()

StatusCode TrigBphysMonitorAlgorithm::fillHistograms ( const EventContext & ctx) const
overridevirtual

adds event to the monitoring histograms

User will overwrite this function. Histogram booking is no longer done in C++. This function is called in execute once the filters are all passed.

Parameters
ctxforwarded from execute
Returns
StatusCode

Implements AthMonitorAlgorithm.

Definition at line 38 of file TrigBphysMonitorAlgorithm.cxx.

38 {
39
40 // fill the container hists
41 ATH_MSG_DEBUG ("fill... Containers");
42 if (fillContainers(ctx).isFailure()) {ATH_MSG_ERROR("Problems filling Container histograms");}
43
44 // fill the chain hists
45 ATH_MSG_DEBUG ("fill... Chains");
46 if (fillChains(ctx).isFailure()) {ATH_MSG_ERROR("Problems filling Chain histograms");}
47
48 // build offline dimuons and fill offline histograms
49 std::vector<std::unique_ptr<xAOD::Vertex>> offlineDimuons;
50 if (buildDimuons(ctx, offlineDimuons).isFailure()) {
51 ATH_MSG_ERROR("Problems building offline dimuons, won't fill corresponding histograms");
52 }
53 else {
54 if (fillOfflineDimuons(ctx, offlineDimuons).isFailure()) {ATH_MSG_ERROR("Problems filling Offline dimuon histograms");}
55 }
56
57 return StatusCode::SUCCESS;
58}
StatusCode fillOfflineDimuons(const EventContext &ctx, const std::vector< std::unique_ptr< xAOD::Vertex > > &dimuonContainer) const
StatusCode buildDimuons(const EventContext &ctx, std::vector< std::unique_ptr< xAOD::Vertex > > &vxContainer) const

◆ fillOfflineDimuonHists()

StatusCode TrigBphysMonitorAlgorithm::fillOfflineDimuonHists ( const EventContext & ,
const std::string & dimuonMonGroupName,
const std::vector< std::unique_ptr< xAOD::Vertex > > & dimuonContainer ) const
private

Definition at line 406 of file TrigBphysMonitorAlgorithm.cxx.

406 {
407 ATH_MSG_DEBUG("Filling " << dimuonMonGroupName << " offline dimuons histograms");
408
409 // Do matching
410 std::vector<const xAOD::Vertex*> matchedDimuons;
411 for(auto& offlineDimuon : dimuonContainer) {
412 if( dimuonMonGroupName == "Any" || matchDimuon(offlineDimuon.get(), dimuonMonGroupName) )
413 matchedDimuons.push_back(offlineDimuon.get());
414 }
415
416 std::string monGroupName = std::string("OfflineDimu_")+dimuonMonGroupName;
417 auto monGroup = getGroup(monGroupName);
418
419 auto ncandidates = Monitored::Scalar<int>("ncandidates",-999);
420 ncandidates = matchedDimuons.size();
421
422 fill(monGroup, ncandidates);
423
424 for(auto matchedDimuon : matchedDimuons) {
425 if (fillVertexHists(monGroup, matchedDimuon, "dimu").isFailure()) {
426 ATH_MSG_ERROR(Form("Problems filling histogram for an offline dimuon vertex in %s", dimuonMonGroupName.c_str()));
427 }
428 if (fillTracksHists(monGroup, matchedDimuon->trackParticleLinks(), "mu", true).isFailure()) {
429 ATH_MSG_ERROR(Form("Problems filling histogram for offline dimuon muons in %s", dimuonMonGroupName.c_str()));
430 }
431 if (fillDiTracksHists(monGroup, matchedDimuon->trackParticleLinks(), "dimu").isFailure()) {
432 ATH_MSG_ERROR(Form("Problems filling histogram for offline dimuon muon pairs in %s", dimuonMonGroupName.c_str()));
433 }
434 }
435
436 return StatusCode::SUCCESS;
437
438}
StatusCode fillDiTracksHists(const ToolHandle< GenericMonitoringTool > &currentMonGroup, const std::vector< ElementLink< xAOD::TrackParticleContainer > > &tpLinkVector, const std::string &name_prefix) const
StatusCode fillTracksHists(const ToolHandle< GenericMonitoringTool > &currentMonGroup, const std::vector< ElementLink< xAOD::TrackParticleContainer > > &tpLinkVector, const std::string &prefix, bool separateHists=false, UInt_t offsetIndex=0) const

◆ fillOfflineDimuons()

StatusCode TrigBphysMonitorAlgorithm::fillOfflineDimuons ( const EventContext & ctx,
const std::vector< std::unique_ptr< xAOD::Vertex > > & dimuonContainer ) const
private

Definition at line 383 of file TrigBphysMonitorAlgorithm.cxx.

383 {
384
385 std::vector<std::string> dimuonMonGroupNames = {"Any"};
386 dimuonMonGroupNames.insert( dimuonMonGroupNames.end(), m_ChainNames_MuMu.begin(), m_ChainNames_MuMu.end() );
387 dimuonMonGroupNames.insert( dimuonMonGroupNames.end(), m_ChainNames_MuMuX.begin(), m_ChainNames_MuMuX.end() );
388
389 for(const auto& dimuonMonGroupName : dimuonMonGroupNames) {
390 ATH_MSG_DEBUG("Process dimuons for " << dimuonMonGroupName);
391 if(dimuonMonGroupName != "Any") {
392 auto& monitoredChain = dimuonMonGroupName;
393 if( !(dimuonMonGroupName == "Any") && !isChainPassed(monitoredChain) ) {
394 ATH_MSG_DEBUG("Chain " << monitoredChain << " is not passed");
395 continue;
396 }
397 }
398 if(fillOfflineDimuonHists(ctx, dimuonMonGroupName, dimuonContainer).isFailure()) {
399 ATH_MSG_ERROR(Form("Problems filling offline dimuon histograms for %s", dimuonMonGroupName.c_str()));
400 }
401 }
402
403 return StatusCode::SUCCESS;
404}
StatusCode fillOfflineDimuonHists(const EventContext &, const std::string &dimuonMonGroupName, const std::vector< std::unique_ptr< xAOD::Vertex > > &dimuonContainer) const

◆ fillTracksHists()

StatusCode TrigBphysMonitorAlgorithm::fillTracksHists ( const ToolHandle< GenericMonitoringTool > & currentMonGroup,
const std::vector< ElementLink< xAOD::TrackParticleContainer > > & tpLinkVector,
const std::string & prefix,
bool separateHists = false,
UInt_t offsetIndex = 0 ) const
private

Definition at line 320 of file TrigBphysMonitorAlgorithm.cxx.

324 {
325 for(UInt_t i = offsetIndex; i < tpLinkVector.size(); ++i) {
326 ATH_CHECK(tpLinkVector.at(i).isValid());
327 std::string curPrefix = prefix;
328 if(separateHists)
329 curPrefix += std::to_string(i+1);
330 if (fillTrkHists(currentMonGroup, *tpLinkVector.at(i), curPrefix).isFailure()) {
331 ATH_MSG_ERROR(Form("Problems filling track histograms"));
332 return StatusCode::FAILURE;
333 }
334 }
335 return StatusCode::SUCCESS;
336}

◆ fillTrigBmumuxTrkHists()

StatusCode TrigBphysMonitorAlgorithm::fillTrigBmumuxTrkHists ( const ToolHandle< GenericMonitoringTool > & currentMonGroup,
const ElementLink< xAOD::TrigBphysContainer > & bphysLink,
UInt_t tracksStartFrom = 2 ) const
private

Definition at line 304 of file TrigBphysMonitorAlgorithm.cxx.

304 {
305
306 const std::vector<ElementLink<xAOD::TrackParticleContainer> > trackVector = (*bphysLink)->trackParticleLinks();
307 if (trackVector.size() < tracksStartFrom) {
308 ATH_MSG_WARNING("Unexpected number of tracks: "<< trackVector.size() << " found, while expected at least " << tracksStartFrom);
309 return StatusCode::SUCCESS;
310 }
311
312 if (fillTracksHists(currentMonGroup, trackVector, "trk", false, tracksStartFrom).isFailure()) {
313 ATH_MSG_ERROR(Form("Problems filling track histograms for a BMuMuX chain"));
314 return StatusCode::FAILURE;
315 }
316
317 return StatusCode::SUCCESS;
318}

◆ fillTrigLeptonHists()

StatusCode TrigBphysMonitorAlgorithm::fillTrigLeptonHists ( const ToolHandle< GenericMonitoringTool > & currentMonGroup,
const ElementLink< xAOD::TrigBphysContainer > & bphysLink,
const std::string & name_prefix ) const
private

Definition at line 287 of file TrigBphysMonitorAlgorithm.cxx.

287 {
288
289 const std::vector<ElementLink<xAOD::TrackParticleContainer> > trackVector = (*bphysLink)->trackParticleLinks();
290 ATH_MSG_DEBUG("fillTrigLeptonHists: number of lepton tracks: " << trackVector.size());
291 if( fillTracksHists(currentMonGroup, trackVector, name_prefix, true).isFailure() ) {
292 ATH_MSG_ERROR(Form("Problems filling muon histograms for a chain"));
293 return StatusCode::FAILURE;
294 }
295 if( fillDiTracksHists(currentMonGroup, trackVector, std::string("di")+name_prefix).isFailure() ) {
296 ATH_MSG_ERROR(Form("Problems filling two-muon histograms for a chain"));
297 return StatusCode::FAILURE;
298 }
299
300 return StatusCode::SUCCESS;
301}

◆ fillTrkHists()

StatusCode TrigBphysMonitorAlgorithm::fillTrkHists ( const ToolHandle< GenericMonitoringTool > & currentMonGroup,
const xAOD::TrackParticle * trk,
const std::string & name_prefix ) const
private

Definition at line 365 of file TrigBphysMonitorAlgorithm.cxx.

365 {
366 if (!trk) {
367 ATH_MSG_ERROR("Null pointer for a track");
368 return StatusCode::FAILURE;
369 }
370 auto trk_pt = Monitored::Scalar<float>(name_prefix+"_pt",-999.);
371 auto trk_eta = Monitored::Scalar<float>(name_prefix+"_eta",-999.);
372 auto trk_d0 = Monitored::Scalar<float>(name_prefix+"_d0",-999.);
373
374 trk_pt = trk->pt() / GeV;
375 trk_eta = trk->eta();
376 trk_d0 = trk->d0();
377
378 fill(currentMonGroup, trk_pt, trk_eta, trk_d0);
379
380 return StatusCode::SUCCESS;
381}
float d0() const
Returns the parameter.
virtual double pt() const override final
The transverse momentum ( ) of the particle.

◆ fillVertexHists()

StatusCode TrigBphysMonitorAlgorithm::fillVertexHists ( const ToolHandle< GenericMonitoringTool > & currentMonGroup,
const xAOD::Vertex * vertex,
const std::string & objStr ) const
private

Definition at line 440 of file TrigBphysMonitorAlgorithm.cxx.

440 {
441
442 xAOD::BPhysHelper dimuonVertexHelper(vertex);
443 if(dimuonVertexHelper.nRefTrks() != 2) {
444 ATH_MSG_WARNING("Unexpected number of refitted tracks at vertex is " << dimuonVertexHelper.nRefTrks());
445 return StatusCode::SUCCESS;
446 }
447
448 using namespace TrigVtx;
449 TLorentzVector dimuonMom = dimuonVertexHelper.totalP(std::array<double,2>{TrigParticleMasses().mass[TrigParticleName::muon], TrigParticleMasses().mass[TrigParticleName::muon]});
450
451 auto dimu_mass = Monitored::Scalar<float>(objStr+"_mass",-999.);
452 auto dimu_pt = Monitored::Scalar<float>(objStr+"_pt",-999.);
453 auto dimu_y = Monitored::Scalar<float>(objStr+"_y",-999.);
454 auto dimu_chi2 = Monitored::Scalar<float>(objStr+"_chi2",-999.);
455 auto dimu_Lxy = Monitored::Scalar<float>(objStr+"_Lxy",-999.);
456 auto dimu_LxySig = Monitored::Scalar<float>(objStr+"_LxySig",-999.);
457
458 dimu_mass = dimuonMom.M() / GeV;
459 dimu_pt = dimuonMom.Pt() / GeV;
460 dimu_chi2 = vertex->chiSquared();
461 dimu_y = dimuonMom.Rapidity();
462 dimu_Lxy = dimuonVertexHelper.lxy(xAOD::BPhysHelper::PV_MIN_A0);
463 dimu_LxySig = dimuonVertexHelper.lxy(xAOD::BPhysHelper::PV_MIN_A0)/dimuonVertexHelper.lxyErr(xAOD::BPhysHelper::PV_MIN_A0);
464
465 fill(currentMonGroup, dimu_mass, dimu_pt, dimu_y, dimu_chi2, dimu_Lxy, dimu_LxySig);
466
467 return StatusCode::SUCCESS;
468}
std::array< double, 6 > mass

◆ filterPassed()

virtual bool AthCommonReentrantAlgorithm< Gaudi::Algorithm >::filterPassed ( const EventContext & ctx) const
inlinevirtualinherited

Definition at line 96 of file AthCommonReentrantAlgorithm.h.

96 {
97 return execState( ctx ).filterPassed();
98 }
virtual bool filterPassed(const EventContext &ctx) const

◆ GetEventInfo()

SG::ReadHandle< xAOD::EventInfo > AthMonitorAlgorithm::GetEventInfo ( const EventContext & ctx) const
inherited

Return a ReadHandle for an EventInfo object (get run/event numbers, etc.)

Parameters
ctxEventContext for the event
Returns
a SG::ReadHandle<xAOD::EventInfo>

Definition at line 111 of file AthMonitorAlgorithm.cxx.

111 {
112 return SG::ReadHandle<xAOD::EventInfo>(m_EventInfoKey, ctx);
113}
SG::ReadHandleKey< xAOD::EventInfo > m_EventInfoKey
Key for retrieving EventInfo from StoreGate.

◆ GetGoodPVs()

std::vector< const xAOD::Vertex * > TrigBphysMonitorAlgorithm::GetGoodPVs ( const xAOD::VertexContainer * pvContainer) const
private

Definition at line 625 of file TrigBphysMonitorAlgorithm.cxx.

625 {
626 std::vector<const xAOD::Vertex*> goodPrimaryVertices;
627 goodPrimaryVertices.reserve(pvContainer->size());
628 for (auto pv : *pvContainer) {
629 xAOD::VxType::VertexType pvType = pv->vertexType();
630 if ( pvType == xAOD::VxType::PriVtx || pvType == xAOD::VxType::PileUp ) {
631 goodPrimaryVertices.push_back(pv);
632 }
633 }
634 return goodPrimaryVertices;
635}
size_type size() const noexcept
Returns the number of elements in the collection.
VertexType
Vertex types.
@ PileUp
Pile-up vertex.
@ PriVtx
Primary vertex.

◆ getGroup()

const ToolHandle< GenericMonitoringTool > & AthMonitorAlgorithm::getGroup ( const std::string & name) const
inherited

Get a specific monitoring tool from the tool handle array.

Finds a specific GenericMonitoringTool instance from the list of monitoring tools (a ToolHandleArray). Throws a FATAL warning if the object found is empty.

Parameters
namestring name of the desired tool
Returns
reference to the desired monitoring tool

Definition at line 168 of file AthMonitorAlgorithm.cxx.

168 {
169 // get the pointer to the tool, and check that it exists
170 auto idx = m_toolLookupMap.find(name);
171 if (ATH_LIKELY(idx != m_toolLookupMap.end())) {
172 return m_tools[idx->second];
173 }
174 else {
175 // treat empty tool handle case as in Monitored::Group
176 if (m_toolLookupMap.empty()) {
177 return m_dummy;
178 }
179
180 if (!isInitialized()) {
182 "It seems that the AthMonitorAlgorithm::initialize was not called "
183 "in derived class initialize method");
184 } else {
185 std::string available = std::accumulate(
186 m_toolLookupMap.begin(), m_toolLookupMap.end(), std::string(""),
187 [](const std::string& s, auto h) { return s + "," + h.first; });
188 ATH_MSG_FATAL("The tool " << name << " could not be found in the tool array of the "
189 << "monitoring algorithm " << m_name << ". This probably reflects a discrepancy between "
190 << "your python configuration and c++ filling code. Note: your available groups are {"
191 << available << "}.");
192 }
193 }
194 return m_dummy;
195}
#define ATH_MSG_FATAL(x)
#define ATH_LIKELY(x)
std::unordered_map< std::string, size_t > m_toolLookupMap
const ToolHandle< GenericMonitoringTool > m_dummy
ToolHandleArray< GenericMonitoringTool > m_tools
Array of Generic Monitoring Tools.

◆ getPvForDimuon_lowestA0()

const xAOD::Vertex * TrigBphysMonitorAlgorithm::getPvForDimuon_lowestA0 ( const xAOD::Vertex * vtx,
const std::vector< const xAOD::Vertex * > & PVs ) const
private

Definition at line 637 of file TrigBphysMonitorAlgorithm.cxx.

637 {
638 std::vector<const xAOD::Vertex*>::const_iterator pv = std::min_element(PVs.begin(), PVs.end(),
639 [&, vtx](const xAOD::Vertex* pv1, const xAOD::Vertex* pv2)
640 { return (std::abs(m_v0Tools->a0(vtx, pv1)) < std::abs(m_v0Tools->a0(vtx, pv2)));}
641 );
642 if(pv == PVs.end()) {
643 return nullptr;
644 }
645 ATH_MSG_VERBOSE("Min-a0 PV has index " << std::distance(PVs.begin(), pv) << ", a0 = " << m_v0Tools->a0(vtx, *pv));
646 return *pv;
647}

◆ getTrigDecisionTool()

const ToolHandle< Trig::TrigDecisionTool > & AthMonitorAlgorithm::getTrigDecisionTool ( ) const
inherited

Get the trigger decision tool member.

The trigger decision tool is used to check whether a specific trigger is passed by an event.

Returns
m_trigDecTool

Definition at line 198 of file AthMonitorAlgorithm.cxx.

198 {
199 return m_trigDecTool;
200}

◆ initialize()

StatusCode TrigBphysMonitorAlgorithm::initialize ( )
overridevirtual

initialize

Returns
StatusCode

Reimplemented from AthMonitorAlgorithm.

Definition at line 24 of file TrigBphysMonitorAlgorithm.cxx.

24 {
25
26 ATH_CHECK( m_TrigBphysContainerKeys.initialize() );
27 ATH_CHECK( m_offlineMuonCollectionKey.initialize() );
28 ATH_CHECK( m_offlineIDTrackCollectionKey.initialize() );
29 ATH_CHECK( m_offlinePvCollectionKey.initialize() );
30
32
33 return StatusCode::SUCCESS;
34
35}
virtual StatusCode initialize() override
initialize
StatusCode initialize(bool used=true)
If this object is used as a property, then this should be called during the initialize phase.

◆ inputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::inputHandles ( ) const
overridevirtualinherited

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.

◆ isChainPassed()

bool TrigBphysMonitorAlgorithm::isChainPassed ( const std::string & chain) const
private

Definition at line 470 of file TrigBphysMonitorAlgorithm.cxx.

470 {
471 // Check if a chain is passed after prescale, accounting for Express Stream prescales if necessary
472 if(!m_requireExplicitESDecision) {
473 return getTrigDecisionTool()->isPassed(chain);
474 }
475 else {
476 const unsigned int passBits = getTrigDecisionTool()->isPassedBits(chain);
477 return passBits & TrigDefs::Express_passed;
478 }
479}

◆ isClonable()

◆ matchDimuon() [1/2]

bool TrigBphysMonitorAlgorithm::matchDimuon ( const xAOD::Vertex * dimuonVertex,
const ElementLink< xAOD::TrigBphysContainer > & bphysLink ) const
private

Definition at line 576 of file TrigBphysMonitorAlgorithm.cxx.

576 {
577 auto offlineTrk1 = dimuonVertex->trackParticle(0);
578 auto offlineTrk2 = dimuonVertex->trackParticle(1);
579 if(!offlineTrk1 || !offlineTrk2) {
580 ATH_MSG_DEBUG("TrackParticle from dimuon is null, won't match");
581 return false;
582 }
583
584 const std::vector<ElementLink<xAOD::TrackParticleContainer> > trackVector = (*bphysLink)->trackParticleLinks();
585 if (trackVector.size() < 2) {
586 ATH_MSG_ERROR("Unexpected number of tracks in a dimuon: " << trackVector.size() << ", won't match");
587 return false;
588 }
589
590 const xAOD::TrackParticle *hltTrk1(nullptr),*hltTrk2(nullptr);
591 // Assume the first two tracks are always muons'
592 if( !trackVector.at(0).isValid() || !trackVector.at(1).isValid() )
593 return false;
594 hltTrk1 = *trackVector.at(0);
595 hltTrk2 = *trackVector.at(1);
596 if (!hltTrk1 || !hltTrk2) {
597 ATH_MSG_ERROR("Null pointer for track in a dimuon!");
598 return false;
599 }
600 if( ( matchTrackParticles(offlineTrk1, hltTrk1) && matchTrackParticles(offlineTrk2, hltTrk2) ) ||
601 ( matchTrackParticles(offlineTrk1, hltTrk2) && matchTrackParticles(offlineTrk2, hltTrk1) ) )
602 return true;
603 else
604 return false;
605}
bool matchTrackParticles(const xAOD::TrackParticle *trk1, const xAOD::TrackParticle *trk2) const
const TrackParticle * trackParticle(size_t i) const
Get the pointer to a given track that was used in vertex reco.

◆ matchDimuon() [2/2]

bool TrigBphysMonitorAlgorithm::matchDimuon ( const xAOD::Vertex * dimuonVertex,
const std::string & chainName ) const
private

Definition at line 564 of file TrigBphysMonitorAlgorithm.cxx.

564 {
565 std::vector< TrigCompositeUtils::LinkInfo<xAOD::TrigBphysContainer> > chainFeatureContainer = getTrigDecisionTool()->features<xAOD::TrigBphysContainer>(chainName, TrigDefs::Physics);
566 for (const auto& featLinkInfo: chainFeatureContainer){
567 if(!featLinkInfo.isValid())
568 return false;
569 const auto featLink = featLinkInfo.link;
570 if( matchDimuon(dimuonVertex, featLink) )
571 return true;
572 }
573 return false;
574}

◆ matchTrackParticles()

bool TrigBphysMonitorAlgorithm::matchTrackParticles ( const xAOD::TrackParticle * trk1,
const xAOD::TrackParticle * trk2 ) const
private

Definition at line 607 of file TrigBphysMonitorAlgorithm.cxx.

607 {
608 double deltaR = xAOD::P4Helpers::deltaR(*trk1, *trk2);
609 if( deltaR < m_deltaRMatchingThreshold)
610 return true;
611 else
612 return false;
613}
Scalar deltaR(const MatrixBase< Derived > &vec) const

◆ msg()

MsgStream & AthCommonMsg< Gaudi::Algorithm >::msg ( ) const
inlineinherited

Definition at line 24 of file AthCommonMsg.h.

24 {
25 return this->msgStream();
26 }

◆ msgLvl()

bool AthCommonMsg< Gaudi::Algorithm >::msgLvl ( const MSG::Level lvl) const
inlineinherited

Definition at line 30 of file AthCommonMsg.h.

30 {
31 return this->msgLevel(lvl);
32 }

◆ outputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::outputHandles ( ) const
overridevirtualinherited

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.

◆ parseList()

StatusCode AthMonitorAlgorithm::parseList ( const std::string & line,
std::vector< std::string > & result ) const
virtualinherited

Parse a string into a vector.

The input string is a single long string of all of the trigger names. It parses this string and turns it into a vector, where each element is one trigger or trigger category.

Parameters
lineThe input string.
resultThe parsed output vector of strings.
Returns
StatusCode

Definition at line 345 of file AthMonitorAlgorithm.cxx.

345 {
346 std::string item;
347 std::stringstream ss(line);
348
349 ATH_MSG_DEBUG( "AthMonitorAlgorithm::parseList()" );
350
351 while ( std::getline(ss, item, ',') ) {
352 std::stringstream iss(item); // remove whitespace
353 iss >> item;
354 result.push_back(item);
355 }
356
357 return StatusCode::SUCCESS;
358}
static Double_t ss

◆ renounce()

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 > AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::renounce ( T & h)
inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

381 {
382 h.renounce();
384 }
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)

◆ renounceArray()

void AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::renounceArray ( SG::VarHandleKeyArray & handlesArray)
inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

364 {
366 }

◆ setFilterPassed()

virtual void AthCommonReentrantAlgorithm< Gaudi::Algorithm >::setFilterPassed ( bool state,
const EventContext & ctx ) const
inlinevirtualinherited

Definition at line 100 of file AthCommonReentrantAlgorithm.h.

100 {
102 }
virtual void setFilterPassed(bool state, const EventContext &ctx) const

◆ sysExecute()

StatusCode AthCommonReentrantAlgorithm< Gaudi::Algorithm >::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 85 of file AthCommonReentrantAlgorithm.cxx.

77{
78 return BaseAlg::sysExecute (ctx);
79}

◆ sysInitialize()

StatusCode AthCommonReentrantAlgorithm< Gaudi::Algorithm >::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 HypoBase, and InputMakerBase.

Definition at line 61 of file AthCommonReentrantAlgorithm.cxx.

107 {
109
110 if (sc.isFailure()) {
111 return sc;
112 }
113
114 ServiceHandle<ICondSvc> cs("CondSvc",name());
115 for (auto h : outputHandles()) {
116 if (h->isCondition() && h->mode() == Gaudi::DataHandle::Writer) {
117 // do this inside the loop so we don't create the CondSvc until needed
118 if ( cs.retrieve().isFailure() ) {
119 ATH_MSG_WARNING("no CondSvc found: won't autoreg WriteCondHandles");
120 return StatusCode::SUCCESS;
121 }
122 if (cs->regHandle(this,*h).isFailure()) {
124 ATH_MSG_ERROR("unable to register WriteCondHandle " << h->fullKey()
125 << " with CondSvc");
126 }
127 }
128 }
129 return sc;
130}
virtual std::vector< Gaudi::DataHandle * > outputHandles() const override

◆ sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::sysStart ( )
overridevirtualinherited

Handle START transition.

We override this in order to make sure that conditions handle keys can cache a pointer to the conditions container.

◆ trigChainsArePassed()

bool AthMonitorAlgorithm::trigChainsArePassed ( const std::vector< std::string > & vTrigNames) const
inherited

Check whether triggers are passed.

For the event, use the trigger decision tool to check that at least one of the triggers listed in the supplied vector is passed.

Parameters
vTrigNamesList of trigger names.
Returns
If empty input, default to true. If at least one trigger is specified, returns whether at least one trigger was passed.

Definition at line 203 of file AthMonitorAlgorithm.cxx.

203 {
204
205
206 // If no triggers were given, return true.
207 if (vTrigNames.empty()) return true;
208
209
210 // Trigger: Check if this Algorithm is being run as an Express Stream job.
211 // Events are entering the express stream are chosen randomly, and by chain,
212 // Hence an additional check should be aplied to see if the chain(s)
213 // monitored here are responsible for the event being selected for
214 // the express stream.
215
216 const auto group = m_trigDecTool->getChainGroup(vTrigNames);
218 const auto passedBits = m_trigDecTool->isPassedBits(group);
219 bool expressPass = passedBits & TrigDefs::Express_passed; //bitwise AND
220 if(!expressPass) {
221 return false;
222 }
223 }
224
225 // monitor the event if any of the chains in the chain group passes the event.
226 return group->isPassed();
227
228}
Gaudi::Property< bool > m_enforceExpressTriggers

◆ updateVHKA()

void AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::updateVHKA ( Gaudi::Details::PropertyBase & )
inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

308 {
309 // debug() << "updateVHKA for property " << p.name() << " " << p.toString()
310 // << " size: " << m_vhka.size() << endmsg;
311 for (auto &a : m_vhka) {
313 for (auto k : keys) {
314 k->setOwner(this);
315 }
316 }
317 }

Member Data Documentation

◆ cfg

TrigBphysMonitorAlgorithm.cfg = MainServicesCfg(flags)

Definition at line 38 of file TrigBphysMonitorAlgorithm.py.

◆ file

str TrigBphysMonitorAlgorithm.file = 'data16_13TeV.00311321.physics_Main.recon.AOD.r9264/AOD.11038520._000001.pool.root.1'

Definition at line 26 of file TrigBphysMonitorAlgorithm.py.

◆ Files

TrigBphysMonitorAlgorithm.Files

Definition at line 29 of file TrigBphysMonitorAlgorithm.py.

◆ flags

TrigBphysMonitorAlgorithm.flags = initConfigFlags()

Definition at line 28 of file TrigBphysMonitorAlgorithm.py.

◆ HISTFileName

TrigBphysMonitorAlgorithm.HISTFileName

Definition at line 31 of file TrigBphysMonitorAlgorithm.py.

◆ isMC

TrigBphysMonitorAlgorithm.isMC

Definition at line 30 of file TrigBphysMonitorAlgorithm.py.

◆ m_ChainNames_Bhh

Gaudi::Property<std::vector<std::string> > TrigBphysMonitorAlgorithm::m_ChainNames_Bhh {this, "ChainNames_Bhh", {}}
private

Definition at line 38 of file TrigBphysMonitorAlgorithm.h.

38{this, "ChainNames_Bhh", {}};

◆ m_ChainNames_ElEl

Gaudi::Property<std::vector<std::string> > TrigBphysMonitorAlgorithm::m_ChainNames_ElEl {this, "ChainNames_ElEl", {}}
private

Definition at line 37 of file TrigBphysMonitorAlgorithm.h.

37{this, "ChainNames_ElEl", {}};

◆ m_ChainNames_MuMu

Gaudi::Property<std::vector<std::string> > TrigBphysMonitorAlgorithm::m_ChainNames_MuMu {this, "ChainNames_MuMu", {}}
private

Definition at line 35 of file TrigBphysMonitorAlgorithm.h.

35{this, "ChainNames_MuMu", {}};

◆ m_ChainNames_MuMuX

Gaudi::Property<std::vector<std::string> > TrigBphysMonitorAlgorithm::m_ChainNames_MuMuX {this, "ChainNames_MuMuX", {}}
private

Definition at line 36 of file TrigBphysMonitorAlgorithm.h.

36{this, "ChainNames_MuMuX", {}};

◆ m_dataType

AthMonitorAlgorithm::DataType_t AthMonitorAlgorithm::m_dataType
protectedinherited

Instance of the DataType_t enum.

Definition at line 356 of file AthMonitorAlgorithm.h.

◆ m_dataTypeStr

Gaudi::Property<std::string> AthMonitorAlgorithm::m_dataTypeStr {this,"DataType","userDefined"}
protectedinherited

DataType string pulled from the job option and converted to enum.

Definition at line 358 of file AthMonitorAlgorithm.h.

358{this,"DataType","userDefined"};

◆ m_defaultLBDuration

Gaudi::Property<float> AthMonitorAlgorithm::m_defaultLBDuration {this,"DefaultLBDuration",60.}
protectedinherited

Default duration of one lumi block.

Definition at line 365 of file AthMonitorAlgorithm.h.

365{this,"DefaultLBDuration",60.};

◆ m_deltaRMatchingThreshold

Gaudi::Property<double> TrigBphysMonitorAlgorithm::m_deltaRMatchingThreshold {this, "deltaRMatchingThreshold", 0.05}
private

Definition at line 48 of file TrigBphysMonitorAlgorithm.h.

48{this, "deltaRMatchingThreshold", 0.05};

◆ m_detailLevel

Gaudi::Property<int> AthMonitorAlgorithm::m_detailLevel {this,"DetailLevel",0}
protectedinherited

Sets the level of detail used in the monitoring.

Definition at line 366 of file AthMonitorAlgorithm.h.

366{this,"DetailLevel",0};

◆ m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_detStore
privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

◆ m_dimuChi2Cut

Gaudi::Property<double> TrigBphysMonitorAlgorithm::m_dimuChi2Cut {this, "dimuChi2Cut", 20.}
private

Definition at line 46 of file TrigBphysMonitorAlgorithm.h.

46{this, "dimuChi2Cut", 20.};

◆ m_dimuMassLower_postfit

Gaudi::Property<double> TrigBphysMonitorAlgorithm::m_dimuMassLower_postfit {this, "dimuMassLower_postfit", 0.}
private

Definition at line 44 of file TrigBphysMonitorAlgorithm.h.

44{this, "dimuMassLower_postfit", 0.};

◆ m_dimuMassLower_prefit

Gaudi::Property<double> TrigBphysMonitorAlgorithm::m_dimuMassLower_prefit {this, "dimuMassLower_prefit", 0.}
private

Definition at line 42 of file TrigBphysMonitorAlgorithm.h.

42{this, "dimuMassLower_prefit", 0.};

◆ m_dimuMassUpper_postfit

Gaudi::Property<double> TrigBphysMonitorAlgorithm::m_dimuMassUpper_postfit {this, "dimuMassUpper_postfit", 15000.}
private

Definition at line 45 of file TrigBphysMonitorAlgorithm.h.

45{this, "dimuMassUpper_postfit", 15000.};

◆ m_dimuMassUpper_prefit

Gaudi::Property<double> TrigBphysMonitorAlgorithm::m_dimuMassUpper_prefit {this, "dimuMassUpper_prefit", 16000.}
private

Definition at line 43 of file TrigBphysMonitorAlgorithm.h.

43{this, "dimuMassUpper_prefit", 16000.};

◆ m_DQFilterTools

ToolHandleArray<IDQFilterTool> AthMonitorAlgorithm::m_DQFilterTools {this,"FilterTools",{}}
protectedinherited

Array of Data Quality filter tools.

Definition at line 346 of file AthMonitorAlgorithm.h.

346{this,"FilterTools",{}};

◆ m_dummy

const ToolHandle<GenericMonitoringTool> AthMonitorAlgorithm::m_dummy
privateinherited

Definition at line 374 of file AthMonitorAlgorithm.h.

◆ m_enforceExpressTriggers

Gaudi::Property<bool> AthMonitorAlgorithm::m_enforceExpressTriggers
privateinherited
Initial value:
{this,
"EnforceExpressTriggers", false,
"Requires that matched triggers made the event enter the express stream"}

Definition at line 377 of file AthMonitorAlgorithm.h.

377 {this,
378 "EnforceExpressTriggers", false,
379 "Requires that matched triggers made the event enter the express stream"};

◆ m_environment

AthMonitorAlgorithm::Environment_t AthMonitorAlgorithm::m_environment
protectedinherited

Instance of the Environment_t enum.

Definition at line 355 of file AthMonitorAlgorithm.h.

◆ m_environmentStr

Gaudi::Property<std::string> AthMonitorAlgorithm::m_environmentStr {this,"Environment","user"}
protectedinherited

Environment string pulled from the job option and converted to enum.

Definition at line 357 of file AthMonitorAlgorithm.h.

357{this,"Environment","user"};

◆ m_EventInfoKey

SG::ReadHandleKey<xAOD::EventInfo> AthMonitorAlgorithm::m_EventInfoKey {this,"EventInfoKey","EventInfo"}
protectedinherited

Key for retrieving EventInfo from StoreGate.

Definition at line 367 of file AthMonitorAlgorithm.h.

367{this,"EventInfoKey","EventInfo"};

◆ m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_evtStore
privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

◆ m_extendedExtraObjects

DataObjIDColl AthCommonReentrantAlgorithm< Gaudi::Algorithm >::m_extendedExtraObjects
privateinherited

Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.

Empty if no symlinks were found.

Definition at line 114 of file AthCommonReentrantAlgorithm.h.

◆ m_fileKey

Gaudi::Property<std::string> AthMonitorAlgorithm::m_fileKey {this,"FileKey",""}
protectedinherited

Internal Athena name for file.

Definition at line 363 of file AthMonitorAlgorithm.h.

363{this,"FileKey",""};

◆ m_lbDurationDataKey

SG::ReadCondHandleKey<LBDurationCondData> AthMonitorAlgorithm::m_lbDurationDataKey {this,"LBDurationCondDataKey","LBDurationCondData","SG Key of LBDurationCondData object"}
protectedinherited

Definition at line 350 of file AthMonitorAlgorithm.h.

351{this,"LBDurationCondDataKey","LBDurationCondData","SG Key of LBDurationCondData object"};

◆ m_lumiDataKey

SG::ReadCondHandleKey<LuminosityCondData> AthMonitorAlgorithm::m_lumiDataKey {this,"LuminosityCondDataKey","LuminosityCondData","SG Key of LuminosityCondData object"}
protectedinherited

Definition at line 348 of file AthMonitorAlgorithm.h.

349{this,"LuminosityCondDataKey","LuminosityCondData","SG Key of LuminosityCondData object"};

◆ m_name

std::string AthMonitorAlgorithm::m_name
privateinherited

Definition at line 371 of file AthMonitorAlgorithm.h.

◆ m_offlineIDTrackCollectionKey

SG::ReadHandleKey<xAOD::TrackParticleContainer> TrigBphysMonitorAlgorithm::m_offlineIDTrackCollectionKey {this, "offlineIDTrackCollectionKey", "InDetTrackParticles"}
private

Definition at line 32 of file TrigBphysMonitorAlgorithm.h.

32{this, "offlineIDTrackCollectionKey", "InDetTrackParticles"};

◆ m_offlineMuonCollectionKey

SG::ReadHandleKey<xAOD::MuonContainer> TrigBphysMonitorAlgorithm::m_offlineMuonCollectionKey {this, "offlineMuonCollectionKey", "Muons"}
private

Definition at line 31 of file TrigBphysMonitorAlgorithm.h.

31{this, "offlineMuonCollectionKey", "Muons"};

◆ m_offlinePvCollectionKey

SG::ReadHandleKey<xAOD::VertexContainer> TrigBphysMonitorAlgorithm::m_offlinePvCollectionKey {this,"offlinePvCollectionKey","PrimaryVertices"}
private

Definition at line 33 of file TrigBphysMonitorAlgorithm.h.

33{this,"offlinePvCollectionKey","PrimaryVertices"};

◆ m_requireExplicitESDecision

Gaudi::Property<bool> TrigBphysMonitorAlgorithm::m_requireExplicitESDecision {this, "requireExplicitESDecision", false}
private

Definition at line 40 of file TrigBphysMonitorAlgorithm.h.

40{this, "requireExplicitESDecision", false};

◆ m_toolLookupMap

std::unordered_map<std::string, size_t> AthMonitorAlgorithm::m_toolLookupMap
privateinherited

Definition at line 372 of file AthMonitorAlgorithm.h.

◆ m_tools

ToolHandleArray<GenericMonitoringTool> AthMonitorAlgorithm::m_tools {this,"GMTools",{}}
protectedinherited

Array of Generic Monitoring Tools.

Definition at line 341 of file AthMonitorAlgorithm.h.

341{this,"GMTools",{}};

◆ m_TrigBphysContainerKeys

SG::ReadHandleKeyArray<xAOD::TrigBphysContainer> TrigBphysMonitorAlgorithm::m_TrigBphysContainerKeys {this, "ContainerNames", {}}
private

Definition at line 29 of file TrigBphysMonitorAlgorithm.h.

29{this, "ContainerNames", {}};

◆ m_trigDecTool

PublicToolHandle<Trig::TrigDecisionTool> AthMonitorAlgorithm::m_trigDecTool
protectedinherited

Tool to tell whether a specific trigger is passed.

Definition at line 345 of file AthMonitorAlgorithm.h.

◆ m_triggerChainString

Gaudi::Property<std::string> AthMonitorAlgorithm::m_triggerChainString {this,"TriggerChain",""}
protectedinherited

Trigger chain string pulled from the job option and parsed into a vector.

Definition at line 360 of file AthMonitorAlgorithm.h.

360{this,"TriggerChain",""};

◆ m_trigLiveFractionDataKey

SG::ReadCondHandleKey<TrigLiveFractionCondData> AthMonitorAlgorithm::m_trigLiveFractionDataKey {this,"TrigLiveFractionCondDataKey","TrigLiveFractionCondData", "SG Key of TrigLiveFractionCondData object"}
protectedinherited

Definition at line 352 of file AthMonitorAlgorithm.h.

353{this,"TrigLiveFractionCondDataKey","TrigLiveFractionCondData", "SG Key of TrigLiveFractionCondData object"};

◆ m_useLumi

Gaudi::Property<bool> AthMonitorAlgorithm::m_useLumi {this,"EnableLumi",false}
protectedinherited

Allows use of various luminosity functions.

Definition at line 364 of file AthMonitorAlgorithm.h.

364{this,"EnableLumi",false};

◆ m_v0Tools

ToolHandle<Trk::V0Tools> TrigBphysMonitorAlgorithm::m_v0Tools {this, "V0Tools", "", "V0 tools to calculate things like Lxy"}
private

Definition at line 52 of file TrigBphysMonitorAlgorithm.h.

52{this, "V0Tools", "", "V0 tools to calculate things like Lxy"};

◆ m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_varHandleArraysDeclared
privateinherited

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vertexFitter

ToolHandle<Trk::TrkVKalVrtFitter> TrigBphysMonitorAlgorithm::m_vertexFitter {this, "VertexFitter", "", "VKalVrtFitter tool to fit tracks into the common vertex"}
private

Definition at line 51 of file TrigBphysMonitorAlgorithm.h.

51{this, "VertexFitter", "", "VKalVrtFitter tool to fit tracks into the common vertex"};

◆ m_vertexPointEstimator

ToolHandle<InDet::VertexPointEstimator> TrigBphysMonitorAlgorithm::m_vertexPointEstimator {this, "VertexPointEstimator", "", "tool to find starting point for the vertex fitter"}
private

Definition at line 50 of file TrigBphysMonitorAlgorithm.h.

50{this, "VertexPointEstimator", "", "tool to find starting point for the vertex fitter"};

◆ m_vhka

std::vector<SG::VarHandleKeyArray*> AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_vhka
privateinherited

Definition at line 398 of file AthCommonDataStore.h.

◆ m_vTrigChainNames

std::vector<std::string> AthMonitorAlgorithm::m_vTrigChainNames
protectedinherited

Vector of trigger chain names parsed from trigger chain string.

Definition at line 361 of file AthMonitorAlgorithm.h.

◆ nightly

str TrigBphysMonitorAlgorithm.nightly = '/cvmfs/atlas-nightlies.cern.ch/repo/data/data-art/CommonInputs/'

Definition at line 25 of file TrigBphysMonitorAlgorithm.py.

◆ trigBphysMonitorAcc

TrigBphysMonitorAlgorithm.trigBphysMonitorAcc = TrigBphysMonConfig(flags)

Definition at line 41 of file TrigBphysMonitorAlgorithm.py.

◆ withDetails

TrigBphysMonitorAlgorithm.withDetails

Definition at line 46 of file TrigBphysMonitorAlgorithm.py.


The documentation for this class was generated from the following files: