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

#include <CscClusterValMonAlg.h>

Inheritance diagram for CscClusterValMonAlg:
Collaboration diagram for CscClusterValMonAlg:

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

 CscClusterValMonAlg (const std::string &name, ISvcLocator *pSvcLocator)
virtual ~CscClusterValMonAlg ()=default
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

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

bool evtSelTriggersPassed () const
Gaudi::Details::PropertyBase & declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
 specialization for handling Gaudi::Property<SG::VarHandleKey>

Private Attributes

ServiceHandle< Muon::IMuonIdHelperSvcm_idHelperSvc {this, "MuonIdHelperSvc", "Muon::MuonIdHelperSvc/MuonIdHelperSvc"}
ToolHandle< ICscStripFitterm_stripFitter {this, "CSCStripFitter", "CSCStripFitter", "Strip fitter Tool"}
ToolHandle< ICscCalibToolm_cscCalibTool {this, "CscCalibTool", "CscCalib", "Calib Tool"}
SG::ReadHandleKey< Muon::CscPrepDataContainerm_cscClusterKey {this,"CSCClusterKey","CSC_Clusters","CSC clusters"}
SG::ReadHandleKey< Muon::CscStripPrepDataContainerm_cscPRDKey {this,"CSCPrepRawDataKey", "CSC_Measurements","CSC PRDs"}
Gaudi::Property< std::vector< std::string > > m_sampSelTriggers {this,"CSCEventSelTriggers",{}}
Gaudi::Property< bool > m_doEvtSel {this,"CSCDoEventSelection",false}
Gaudi::Property< double > m_qmaxADCCut {this,"CSCQmaxCutADC",100}
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 22 of file CscClusterValMonAlg.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

◆ CscClusterValMonAlg()

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

Definition at line 32 of file CscClusterValMonAlg.cxx.

32 :
33 AthMonitorAlgorithm(name,pSvcLocator)
34 { }
AthMonitorAlgorithm(const std::string &name, ISvcLocator *pSvcLocator)
Constructor.

◆ ~CscClusterValMonAlg()

virtual CscClusterValMonAlg::~CscClusterValMonAlg ( )
virtualdefault

Member Function Documentation

◆ 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.

◆ 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}

◆ evtSelTriggersPassed()

bool CscClusterValMonAlg::evtSelTriggersPassed ( ) const
private

Definition at line 452 of file CscClusterValMonAlg.cxx.

452 {
453
454 if(!m_doEvtSel) return true;
455
456 for(const auto& trig : m_sampSelTriggers) {
457 if(m_trigDecTool->isPassed(trig,TrigDefs::eventAccepted)){
458 return true;
459 }
460 }
461 return false;
462
463} // end evtSelTriggersPassed
PublicToolHandle< Trig::TrigDecisionTool > m_trigDecTool
Tool to tell whether a specific trigger is passed.
Gaudi::Property< bool > m_doEvtSel
Gaudi::Property< std::vector< std::string > > m_sampSelTriggers

◆ 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}
#define ATH_MSG_DEBUG(x)
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.
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.

◆ fillHistograms()

StatusCode CscClusterValMonAlg::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 58 of file CscClusterValMonAlg.cxx.

58 {
59
60
61 StatusCode sc = StatusCode::SUCCESS;
62
63 // check if event passed sample-selection triggers
64 if(m_doEvtSel) { if(!evtSelTriggersPassed()) return sc; }
65
66 // retrieve cluster / strip collection
67 SG::ReadHandle<CscPrepDataContainer> cscCluster(m_cscClusterKey, ctx);
68 SG::ReadHandle<CscStripPrepDataContainer> cscStrip(m_cscPRDKey, ctx);
69
70 // we can do (some) monitoring plots with just the cluster
71 // ideally we need both the cluster and the strips that make up that cluster
72 //FillCSCClusters(*cscCluster.cptr(), *cscStrip.cptr());
73
74 if(!(cscCluster.isValid())) {
75 ATH_MSG_ERROR("evtStore() does not contain csc prd Collection with name "<< m_cscClusterKey);
76 return StatusCode::FAILURE;
77 }
78
79 ATH_MSG_DEBUG ( " Size of Cluster Collection : " << cscCluster->size() );
80 ATH_MSG_DEBUG ( " Size of Strip Collection : " << cscStrip->size() );
81
82 MonStruct monstruct;
83
84 for ( CscPrepDataContainer::const_iterator Icol = cscCluster->begin(); Icol != cscCluster->end(); ++Icol )
85 {
86 const CscPrepDataCollection& clus = **Icol;
87
88 // arrays to hold cluster-count
89 // 32 chambers and 8 layers (each has one extra - index '0' is not counted)
90 int clusCount[33][9], sigclusCount[33][9];
91 unsigned int nEtaClusWidthCnt = 0, nPhiClusWidthCnt = 0;
92 for(unsigned int kl = 0; kl < 33; kl++ ) {
93 for(unsigned int km = 0; km < 9; km++ ) {
94 clusCount[kl][km] = 0;
95 sigclusCount[kl][km] = 0;
96 }
97 }
98
99 float stripsSum_EA = 0.;
100 float stripsSum_EAtest = -50.;
101 float stripsSum_EC = 0.;
102 float stripsSum_ECtest = -50.;
103
104 ATH_MSG_DEBUG ( " Begin loop over clusters ============================");
105 for ( CscPrepDataCollection::const_iterator Itclu = clus.begin(); Itclu != clus.end(); ++Itclu )
106 {
107 const CscPrepData& iClus = **Itclu;
108 const std::vector<Identifier>& stripIds = iClus.rdoList();
109 float clu_charge = iClus.charge();
110 auto clu_time = Monitored::Scalar<float>("clu_time", (iClus.time()));
111
112 ATH_MSG_DEBUG(" cluster charge = " << clu_charge << "\t cluster time = " << clu_time );
113
114 unsigned int noStrips = stripIds.size(); // no. of strips in this cluster = stripIds.size()
115 auto noStrips_mon = Monitored::Scalar<int> ("noStrips_mon",noStrips);
116 Identifier clusId = iClus.identify();
117
118 // get the cluster coordinates
119 int stationName = m_idHelperSvc->cscIdHelper().stationName(clusId);
120 std::string stationString = m_idHelperSvc->cscIdHelper().stationNameString(stationName);
121 int chamberType = stationString == "CSS" ? 0 : 1;
122 int stationEta = m_idHelperSvc->cscIdHelper().stationEta(clusId);
123 int stationPhi = m_idHelperSvc->cscIdHelper().stationPhi(clusId);
124 int wireLayer = m_idHelperSvc->cscIdHelper().wireLayer(clusId);
125 int measuresPhi = m_idHelperSvc->cscIdHelper().measuresPhi(clusId);
126
127 auto x = Monitored::Scalar<float> ("x",iClus.globalPosition().x());
128 auto y = Monitored::Scalar<float> ("y",iClus.globalPosition().y());
129 auto z = Monitored::Scalar<float> ("z",iClus.globalPosition().z());
130 auto r = Monitored::Scalar<float> ("r",std::hypot(x,y));
131
132 fill("CscClusMonitor",z,r);
133 fill("CscClusMonitor",y,x);
134
135 // convert to my coordinates
136 int sectorNo = stationEta * (2 * stationPhi - chamberType); // [-16 -> -1] and [+1 -> +16]
137 auto secLayer = Monitored::Scalar<float> ("secLayer", (sectorNo + 0.2 * (wireLayer - 1) + 0.1) );
138 int xfac = measuresPhi ? -1 : 1; // [-1 -> -48] / [+1 -> +192]
139
140 //total cluster width (EA and EC) calculation
141 if(secLayer > 0.) {
142 stripsSum_EA = stripsSum_EA + noStrips;
143 }
144 if(stripsSum_EA > stripsSum_EAtest) {
145 stripsSum_EAtest = stripsSum_EA;
146 }
147
148 if(secLayer < 0. || secLayer == 0.) {
149 stripsSum_EC = stripsSum_EC + noStrips;
150 }
151 if(stripsSum_EC > stripsSum_ECtest) {
152 stripsSum_ECtest = stripsSum_EC;
153 }
154
155 // compute the indices to store cluster count
156 int ns = sectorNo < 0 ? sectorNo*(-1) : sectorNo+16; // [-16 -> -1] shifted to [1 -> 16] and [+1 -> +16] shifted to [+17 -> +32]
157 int nl = (measuresPhi ? wireLayer : wireLayer+4); // [ 1 -> 4] (phi-layers) and [5 -> 8] (eta-layers)
158
159 // increment the cluster-count for this layer
160 clusCount[ns][nl]++;
161
162 // indices for ns = [+1 -> +32]; 32 places (index '0' is not counted); allocated 33 places
163 // indices for nl = [+1 -> +8]; 8 places (index '0' is not counted); allocated 9 places
164 ATH_MSG_DEBUG(" ns = " << ns << "\tm_nl = " << nl << "\tm_sec = " << sectorNo << "\t m_lay= " << wireLayer << "\tmPhi = " << measuresPhi);
165
166
167 // check the cluster status; probably need to read status info from jobOptions - not done for the moment
168 // status = Muon::CscStatusUnspoiled (i.e 0) or Muon::CscStatusSplitUnspoiled (i.e 10) are considered good for precision clusters
169 // status = Muon::CscStatusSimple (i.e 1) could be good for non-precision clusters (i.e for phi-layers)
170 std::string stat = Muon::toString(iClus.status());
171 bool cluster_status = ( (stat == "unspoiled") ||
172 (stat == "unspoiled with split") ||
173 (stat == "simple")
174 ) ? true : false;
175
176 // Also need at least three strips in an eta-cluster to compute Q_max, Q_left and Q_right
177 bool eta_cluster_status = cluster_status && ( noStrips > 2 ) && (measuresPhi == 0);
178
179 // Also need at least one strip in a phi-cluster to compute Q_max = Q_sum
180 bool phi_cluster_status = cluster_status && ( noStrips > 0 ) && (measuresPhi == 1);
181
182 ATH_MSG_DEBUG ( " ClusterStatus eta = " << eta_cluster_status << " ,phi = " << phi_cluster_status);
183 ATH_MSG_DEBUG ( " ClusterID (eta:" << stationEta << ",phi:" << stationPhi << ",type:" << chamberType << ", measPhi: "
184 << measuresPhi << ",wire:" << wireLayer << ") = " << secLayer << " status = "
185 << stat << " #of strips = " << noStrips );
186
187 // if cluster is okay get Qmax, Qleft, Qright and Qsum = (Qmax + Qleft + Qright)
188 if(eta_cluster_status || phi_cluster_status ) {
189 const CscStripPrepDataCollection* pcol(nullptr);
190 bool found_id = true;
191 std::vector <const CscStripPrepData*> stripVec;
192 std::vector <float> fStripIDs;
193
194 float maxStripCharge = 0., maxStipId = 0.;
195 int sIdx = 0, mxIdx = 0; // index-counter and index of max strip in the vector of Id's
196
197 auto clus_phi_mon = Monitored::Scalar<int>("clus_phi_mon", (int)measuresPhi);
198 auto clus_eta_mon = Monitored::Scalar<int>("clus_eta_mon", (int)(!measuresPhi));
199 fill("CscClusMonitor",noStrips_mon, secLayer,clus_phi_mon, clus_eta_mon);
200
201 // fill cluster width (no. of strips per cluster)
202 if(measuresPhi) {
203 nPhiClusWidthCnt++;
204 } else {
205 nEtaClusWidthCnt++;
206 }
207
208 // Loop over strip id's vector / strip collection and match the id's from vector with strips in collection
209 for ( std::vector<Identifier>::const_iterator sId = stripIds.begin(); sId != stripIds.end(); ++sId, sIdx++ ) {
210 Identifier id = *sId; // for strip Id's
211 int thisStrip = m_idHelperSvc->cscIdHelper().strip(id);
212 auto stripid = Monitored::Scalar<int> ("stripid", thisStrip*xfac);// x-axis fill value
213 fStripIDs.push_back(stripid);
214
215 fill("CscClusMonitor",stripid, secLayer);
216
217 if(!pcol) {
218 const CscStripPrepDataCollection* icol = cscStrip->indexFindPtr(clus.identifyHash());
219 if ( icol == nullptr ) {
220 found_id = false;
221 break; // could not identify the strips
222 } else {
223 pcol = icol;
224 }
225 } // end if !pcol
226
227 bool found_strip = false;
228 float maxsampChVal = 0.;
229 if(found_id) {
230 for ( CscStripPrepDataCollection::const_iterator istrip= pcol->begin(); istrip != pcol->end(); ++ istrip ) {
231 found_strip = ( *istrip )->identify() == id ;
232 if(found_strip) {
233 stripVec.push_back(*istrip);
234 std::vector<float> samp_charges = ( *istrip )->sampleCharges();
235 for(unsigned int i = 0; i < samp_charges.size(); i++ ) {
236 if(samp_charges[i] > maxsampChVal) maxsampChVal = samp_charges[i];
237 }
238 if(maxsampChVal > maxStripCharge ) {
239 maxStripCharge = maxsampChVal;
240 maxStipId = stripid;
241 mxIdx = sIdx;
242 }
243 break; // break from inner loop
244 }
245 } // end for loop on strip collection
246 ATH_MSG_DEBUG ( " " << (found_strip? "FoundStrip " : "NoStripFound ") << " with max sampling = " << maxsampChVal);
247 } // end if found_id
248 } // end for loop over strips
249
250 auto fStripIDs_col = Monitored::Collection("fStripIDs_col",fStripIDs);
251 ATH_MSG_DEBUG ( " Max Strip charge = " << maxStripCharge << " and strip Id = " << maxStipId << " and index = " << mxIdx);
252 float qmax = 0., qleft = 0., qright = 0., qsum = 0.;
253 // if we are here and loop over strips is successful we should have found_id = true
254 // and the size of strip-ID-vector == size of strips-vector
255 bool size_ids_coll = (noStrips == stripVec.size() ? true : false) ;
256
257 if(found_id && size_ids_coll ) {
258 // store results of three strips (Qmax, Qleft, Qright)
259 std::vector<ICscStripFitter::Result> res;
260 res.resize(3);
261 bool range_check = (mxIdx > -1) && (mxIdx < int(noStrips));
262
263 ATH_MSG_DEBUG ( " Range check = (" << mxIdx << " > -1 ) && (" << mxIdx << " < " << noStrips << " ) = " << range_check
264 << "\t size of vec check " << noStrips << " == " << stripVec.size());
265
266 if( range_check ) {
267 // fit Q_left fit
268 if(mxIdx-1 >= 0 ) {
269 res[0] = m_stripFitter->fit(*stripVec[mxIdx-1]);
270 qleft = res[0].charge;
271 qsum += qleft;
272 ATH_MSG_DEBUG ( " Left Strip q +- dq = " << res[0].charge << " +- " << res[0].dcharge << "\t t +- dt = "
273 << res[0].time << " +- " << res[0].dtime << "\t w +- dw = " << res[0].width << " +- "
274 << res[0].dwidth << "\t status= " << res[0].status << "\t chisq= " << res[0].chsq);
275 }// end if q_left
276 // fit Q_max strip
277 res[1] = m_stripFitter->fit(*stripVec[mxIdx]);
278 qmax = res[1].charge;
279 qsum += qmax;
280 ATH_MSG_DEBUG ( " Peak Strip q +- dq = " << res[1].charge << " +- " << res[1].dcharge << "\t t +- dt = "
281 << res[1].time << " +- " << res[1].dtime << "\t w +- dw = " << res[1].width << " +- "
282 << res[1].dwidth << "\t status= " << res[1].status << "\t chisq= " << res[1].chsq);
283 // fit Q_right strip
284 if(mxIdx+1 < int(noStrips)) {
285 res[2] = m_stripFitter->fit(*stripVec[mxIdx+1]);
286 qright = res[2].charge;
287 qsum += qright;
288 ATH_MSG_DEBUG ( " Right Strip q +- dq = " << res[2].charge << " +- " << res[2].dcharge << "\t t +- dt = "
289 << res[2].time << " +- " << res[2].dtime << "\t w +- dw = " << res[2].width << " +- "
290 << res[2].dwidth << "\t status= " << res[2].status << "\t chisq= " << res[2].chsq);
291 } // end if q_right
292 } // end if range_check
293
294 // not used at the moment
295 // 1 e = 1.602176487 10^{-19} C = 1.6022 x 10^{-4} fC
296 // float m_fCperElectron = 1.6022e-4; // multiply # of electrons by this number to get fC
297
298 float kiloele = 1.0e-3; // multiply # of electrons by this number to get kiloElectrons (1 ke = 1 ADC)
299
300
301 // Assume 1000 e = 1 ADC for now = 1000 x 1.6022 x 10^{-4} fC = 0.16022 fC
302 // convert qmax, qleft, qright into ADC
303
304 auto QmaxADC = Monitored::Scalar<float>("QmaxADC", (qmax * kiloele));
305 auto QsumADC = Monitored::Scalar<float>("QsumADC", (qsum * kiloele));
306
307
308 // check if signal or noise
309 // QmaxADC > m_qmaxADCCut is signal
310 bool signal = QmaxADC > m_qmaxADCCut;
311
312 // fill signal/noise histograms
313 auto signal_mon = Monitored::Scalar<int>("signal_mon",(int)signal);
314 auto noise_mon = Monitored::Scalar<int>("noise_mon",(int)!(signal));
315 auto clus_phi = Monitored::Scalar<int>("clus_phi", (int)measuresPhi );
316 auto clus_eta = Monitored::Scalar<int>("clus_eta", (int)(!measuresPhi) );
317 auto clus_phiSig = Monitored::Scalar<int>("clus_phiSig", (int)measuresPhi && (signal));
318 auto clus_etaSig = Monitored::Scalar<int>("clus_etaSig", (int)(!measuresPhi) && (signal));
319 auto clus_phiNoise = Monitored::Scalar<int>("clus_phiNoise", (int)measuresPhi && !(signal));
320 auto clus_etaNoise = Monitored::Scalar<int>("clus_etaNoise", (int)(!measuresPhi) && !(signal));
321 auto sideA = Monitored::Scalar<int>("sideA",(int)((stationEta==1) && (signal)));
322 auto sideC = Monitored::Scalar<int>("sideC",(int)((stationEta==-1) && (signal)));
323 auto sideA_phiSig = Monitored::Scalar<int>("sideA_phiSig", (int)(signal && stationEta==1 && !measuresPhi) );
324 auto sideC_phiSig = Monitored::Scalar<int>("sideC_phiSig", (int)(signal && stationEta==-1 && !measuresPhi) );
325
326 if(signal) sigclusCount[ns][nl]++;
327
328 auto clu_charge_kiloele = Monitored::Scalar<float>("clu_charge_kiloele", (iClus.charge()*kiloele));
329
330 fill("CscClusMonitor",fStripIDs_col, QmaxADC, secLayer, noStrips_mon, QsumADC, clu_time, clu_charge_kiloele, clus_phi, clus_eta, clus_phiSig, clus_etaSig, clus_phiNoise, clus_etaNoise, signal_mon, noise_mon, sideA, sideC, sideA_phiSig, sideC_phiSig);
331
332 ATH_MSG_DEBUG ( " End of strip fits " );
333
334 } // if found_id
335 } // end if cluster_status
336
337 auto stripsSum_EA_mon = Monitored::Scalar<float> ("stripsSum_EA_mon",stripsSum_EA);
338 auto stripsSum_EC_mon = Monitored::Scalar<float> ("stripsSum_EC_mon",stripsSum_EC);
339 fill("CscClusMonitor",stripsSum_EA_mon, stripsSum_EC_mon);
340
341 } // end for loop over prep-data collection
342 ATH_MSG_DEBUG ( " End loop over clusters ============================");
343 auto nPhiClusWidthCnt_mon = Monitored::Scalar<int> ("nPhiClusWidthCnt_mon",nPhiClusWidthCnt);
344 auto nEtaClusWidthCnt_mon = Monitored::Scalar<int> ("nEtaClusWidthCnt_mon",nEtaClusWidthCnt);
345 fill("CscClusMonitor",nPhiClusWidthCnt_mon,nEtaClusWidthCnt_mon);
346
347 // Fill cluster counts
348 int numeta = 0, numphi = 0;
349 int numetasignal = 0, numphisignal = 0;
350
351 //loop over chambers
352 for(int kl = 1; kl < 33; kl++ ) {
353
354 // loop over layers
355 int eta_hits[4] = {0,0,0,0};
356 bool chamber_empty = true;
357 int sec = kl < 17 ? kl*(-1) : kl; // [1->16](-side) [17-32] (+side)
358 for(int km = 1; km < 9; km++ ) {
359 int lay = (km > 4 && km < 9) ? km-4 : km; // 1,2,3,4 (phi-layers) 5-4, 6-4, 7-4, 8-4 (eta-layers)
360 bool mphi = (km > 0 && km < 5) ? true : false; // 1,2,3,4 (phi-layers) 5,6,7,8 (eta-layers)
361 std::string wlay = mphi ? "Phi-Layer " : "Eta-Layer: ";
362 int count = clusCount[kl][km];
363 int scount = sigclusCount[kl][km];
364
365 if(count){
366 ATH_MSG_DEBUG ("sec[" << sec << "]\t" << wlay << "[" << lay << "] = " << monstruct.secLayer.back() << "= " << "\tNsig = " << scount << ", Ntot = " << count);
367 if(mphi){
368 numphi += count;
369 if(scount){
370 chamber_empty = false;
371 numphisignal += scount;
372 }
373 }
374 else{
375 numeta += count;
376 if(scount){
377 eta_hits[lay-1]++;
378 chamber_empty = false;
379 numetasignal +=scount;
380 }
381 }
382 ATH_MSG_DEBUG ( wlay << "Counts sec: [" << kl-16 << "]\tlayer: [" << km << "] = " << monstruct.secLayer.back() << "\t = " << count << "\t" << scount);
383 }
384
385 monstruct.count_mon.push_back(count);
386 monstruct.scount_mon.push_back(scount);
387 monstruct.count_diff.push_back(count-scount);
388 monstruct.mphi_true.push_back((int)mphi && count);
389 monstruct.mphi_false.push_back((int)!(mphi) && count);
390 monstruct.scount_phi_true.push_back((int)mphi && count && scount);
391 monstruct.scount_phi_false.push_back((int)mphi && count && !scount);
392 monstruct.scount_eta_true.push_back((int)!(mphi) && count && scount);
393 monstruct.scount_eta_false.push_back((int)!(mphi) && count && !scount);
394 monstruct.secLayer.push_back((sec + 0.2 * (lay - 1) + 0.1));
395
396 }// end loop over layers
397
398 int segNum_new = -999.;
399 if(!chamber_empty){
400 std::ostringstream nseglist;
401 std::bitset<4> segNum;
402 for(unsigned int mm = 0; mm < 4; mm++) {
403 bool set = (eta_hits[mm] > 0 ? true : false);
404 if(set) segNum.set(mm);
405 nseglist << (set ? "1" : "0");
406 }
407 segNum_new = segNum.to_ulong();
408 ATH_MSG_DEBUG("segments= " << nseglist.str() << "\t = " << segNum.to_ulong());
409 }
410 auto segNum_mon = Monitored::Scalar<int>("segNum_mon", segNum_new);
411 auto sec_mon = Monitored::Scalar<float>("sec_mon",sec+0.3);
412 fill("CscClusMonitor", segNum_mon, sec_mon);
413 } // end loop over chambers
414
415 ATH_MSG_DEBUG(" numphi = " << numphi << "\t numeta = " << numeta << "\tm_sphi = " << numphisignal << "\t m_seta = " << numetasignal);
416 auto numphi_mon = Monitored::Scalar<int>("numphi_mon", numphi);
417 auto numeta_mon = Monitored::Scalar<int>("numeta_mon", numeta);
418 auto numphi_sig_mon = Monitored::Scalar<int>("numphi_sig_mon", numphisignal);
419 auto numeta_sig_mon = Monitored::Scalar<int>("numeta_sig_mon", numetasignal);
420 auto numphi_numeta_mon = Monitored::Scalar<int>("numphi_numeta_mon", numphi+numeta);
421 auto numphi_numeta_sig_mon = Monitored::Scalar<int>("numphi_numeta_sig_mon", numphisignal+numetasignal);
422 auto num_num_noise_mon = Monitored::Scalar<int>("num_num_noise_mon", (numphi-numphisignal)+(numeta-numetasignal));
423 auto numphi_diff_mon = Monitored::Scalar<int>("numphi_diff_mon", numphi-numphisignal);
424 auto numeta_diff_mon = Monitored::Scalar<int>("numeta_diff_mon", numeta-numetasignal);
425
426 fill("CscClusMonitor",numphi_mon,numeta_mon,numphi_sig_mon,numeta_sig_mon,numphi_numeta_mon,numphi_numeta_sig_mon,num_num_noise_mon,numphi_diff_mon,numeta_diff_mon);
427
428 } // end for loop over prep-data container
429
430 auto count_mon = Monitored::Collection("count_mon", monstruct.count_mon);
431 auto scount_mon = Monitored::Collection("scount_mon", monstruct.scount_mon);
432 auto count_diff = Monitored::Collection("count_diff", monstruct.count_diff);
433 auto mphi_true = Monitored::Collection("mphi_true", monstruct.mphi_true);
434 auto mphi_false = Monitored::Collection("mphi_false", monstruct.mphi_false);
435 auto scount_phi_true = Monitored::Collection("scount_phi_true", monstruct.scount_phi_true);
436 auto scount_phi_false = Monitored::Collection("scount_phi_false", monstruct.scount_phi_false);
437 auto scount_eta_true = Monitored::Collection("scount_eta_true", monstruct.scount_eta_true);
438 auto scount_eta_false = Monitored::Collection("scount_eta_false", monstruct.scount_eta_false);
439 auto secLayer = Monitored::Collection("secLayer", monstruct.secLayer);
440 auto tmp_val_mon = Monitored::Collection("tmp_val_mon", monstruct.tmp_val_mon);
441 fill("CscClusMonitor", count_mon, scount_mon, count_diff, secLayer, mphi_true, mphi_false, scount_phi_true, scount_phi_false, scount_eta_true, scount_eta_false);
442
443 ATH_MSG_DEBUG ( " END EVENT ============================");
444
445 return sc;
446
447}
#define ATH_MSG_ERROR(x)
double charge(const T &p)
Definition AtlasPID.h:997
std::pair< std::vector< unsigned int >, bool > res
static Double_t sc
const double width
#define y
#define x
#define z
SG::ReadHandleKey< Muon::CscPrepDataContainer > m_cscClusterKey
Gaudi::Property< double > m_qmaxADCCut
ServiceHandle< Muon::IMuonIdHelperSvc > m_idHelperSvc
ToolHandle< ICscStripFitter > m_stripFitter
SG::ReadHandleKey< Muon::CscStripPrepDataContainer > m_cscPRDKey
DataModel_detail::const_iterator< DataVector > const_iterator
Definition DataVector.h:838
const_iterator end() const noexcept
const_iterator begin() const noexcept
virtual const Amg::Vector3D & globalPosition() const override final
Returns the global position.
double time() const
Returns the time.
int charge() const
Returns the charge.
CscClusterStatus status() const
Returns the Csc status (position measurement) flag.
virtual IdentifierHash identifyHash() const override final
Identifier identify() const
return the identifier
const std::vector< Identifier > & rdoList() const
return the List of rdo identifiers (pointers)
int r
Definition globals.cxx:22
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.
int count(std::string s, const std::string &regx)
count how many occurances of a regx are in a string
Definition hcg.cxx:146
constexpr std::enable_if_t< is_bitmask_v< E >, E & > set(E &lhs, E rhs)
Convenience function to set bits in a class enum bitmask.
Definition bitmask.h:232
::StatusCode StatusCode
StatusCode definition for legacy code.
ValuesCollection< T > Collection(std::string name, const T &collection)
Declare a monitored (double-convertible) collection.
MuonPrepDataCollection< CscPrepData > CscPrepDataCollection
MuonPrepDataCollection< CscStripPrepData > CscStripPrepDataCollection
std::string toString(CscStripStatus cstat)
Return a string description of a CSC cluster status flag.
unsigned int constexpr sideC
Definition RPDUtils.h:15
unsigned int constexpr sideA
Definition RPDUtils.h:16
count_diff(count_in, count_ref, total_in, total_ref, thr_frac, thr_num)
Definition chainDump.py:217

◆ 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.

◆ 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.

◆ 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 CscClusterValMonAlg::initialize ( )
overridevirtual

initialize

Returns
StatusCode

Reimplemented from AthMonitorAlgorithm.

Definition at line 36 of file CscClusterValMonAlg.cxx.

36 {
37
38
39 ATH_MSG_INFO ( "Initializing : " << name() );
40 ATH_MSG_INFO ( "CSCClusterKey : " << m_cscClusterKey );
41 ATH_MSG_INFO ( "CSCPrepRawDataKey : " << m_cscPRDKey );
42
43 ATH_CHECK(m_idHelperSvc.retrieve());
44
45 ATH_CHECK(m_stripFitter.retrieve());
46 ATH_MSG_INFO ( "CSCStripFitter : " << "Using Fitter with name \"" << m_stripFitter->name() << "\"" );
47
48 if( m_doEvtSel ) ATH_CHECK(m_trigDecTool.retrieve());
49
50 ATH_CHECK(m_cscCalibTool.retrieve());
51 ATH_CHECK(m_cscClusterKey.initialize((m_idHelperSvc->hasCSC())));
52 ATH_CHECK(m_cscPRDKey.initialize((m_idHelperSvc->hasCSC())));
53
55}
#define ATH_CHECK
Evaluate an expression and check for errors.
#define ATH_MSG_INFO(x)
virtual StatusCode initialize() override
initialize
ToolHandle< ICscCalibTool > m_cscCalibTool

◆ 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.

◆ isClonable()

◆ 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

◆ m_cscCalibTool

ToolHandle<ICscCalibTool> CscClusterValMonAlg::m_cscCalibTool {this, "CscCalibTool", "CscCalib", "Calib Tool"}
private

Definition at line 40 of file CscClusterValMonAlg.h.

40{this, "CscCalibTool", "CscCalib", "Calib Tool"};

◆ m_cscClusterKey

SG::ReadHandleKey<Muon::CscPrepDataContainer> CscClusterValMonAlg::m_cscClusterKey {this,"CSCClusterKey","CSC_Clusters","CSC clusters"}
private

Definition at line 43 of file CscClusterValMonAlg.h.

43{this,"CSCClusterKey","CSC_Clusters","CSC clusters"};

◆ m_cscPRDKey

SG::ReadHandleKey<Muon::CscStripPrepDataContainer> CscClusterValMonAlg::m_cscPRDKey {this,"CSCPrepRawDataKey", "CSC_Measurements","CSC PRDs"}
private

Definition at line 44 of file CscClusterValMonAlg.h.

44{this,"CSCPrepRawDataKey", "CSC_Measurements","CSC PRDs"};

◆ 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_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_doEvtSel

Gaudi::Property<bool> CscClusterValMonAlg::m_doEvtSel {this,"CSCDoEventSelection",false}
private

Definition at line 47 of file CscClusterValMonAlg.h.

47{this,"CSCDoEventSelection",false};

◆ 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_idHelperSvc

ServiceHandle<Muon::IMuonIdHelperSvc> CscClusterValMonAlg::m_idHelperSvc {this, "MuonIdHelperSvc", "Muon::MuonIdHelperSvc/MuonIdHelperSvc"}
private

Definition at line 35 of file CscClusterValMonAlg.h.

35{this, "MuonIdHelperSvc", "Muon::MuonIdHelperSvc/MuonIdHelperSvc"};

◆ 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_qmaxADCCut

Gaudi::Property<double> CscClusterValMonAlg::m_qmaxADCCut {this,"CSCQmaxCutADC",100}
private

Definition at line 48 of file CscClusterValMonAlg.h.

48{this,"CSCQmaxCutADC",100};

◆ m_sampSelTriggers

Gaudi::Property<std::vector<std::string> > CscClusterValMonAlg::m_sampSelTriggers {this,"CSCEventSelTriggers",{}}
private

Definition at line 46 of file CscClusterValMonAlg.h.

46{this,"CSCEventSelTriggers",{}};

◆ m_stripFitter

ToolHandle<ICscStripFitter> CscClusterValMonAlg::m_stripFitter {this, "CSCStripFitter", "CSCStripFitter", "Strip fitter Tool"}
private

Definition at line 37 of file CscClusterValMonAlg.h.

37{this, "CSCStripFitter", "CSCStripFitter", "Strip fitter Tool"};

◆ 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_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_varHandleArraysDeclared

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

Definition at line 399 of file AthCommonDataStore.h.

◆ 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.


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