DQMHistAnalysisKLMModule Class Reference

Analysis of KLM DQM histograms. More...

#include <DQMHistAnalysisKLM.h>

Inheritance diagram for DQMHistAnalysisKLMModule:

Public Types
enum	EStatus {   c_StatusTooFew = 0 ,   c_StatusDefault = 1 ,   c_StatusGood = 2 ,   c_StatusWarning = 3 ,   c_StatusError = 4 }
	Status flag of histogram/canvas. More...
enum	EStatusColor {   c_ColorTooFew = kGray ,   c_ColorDefault = kWhite ,   c_ColorGood = kGreen ,   c_ColorWarning = kYellow ,   c_ColorError = kRed }
	Status colors of histogram/canvas (corresponding to status) More...
enum	ERefScaling {   c_RefScaleNone = 0 ,   c_RefScaleEntries = 1 ,   c_RefScaleMax = 2 }
	Reference plot scaling type. More...
typedef std::map< std::string, HistObject >	HistList
	The type of list of histograms.
typedef std::map< std::string, MonitoringObject >	MonObjList
	The type of list of MonitoringObjects.
typedef std::map< std::string, HistDelta >	DeltaList
	The type of list of delta settings and histograms.
typedef std::map< std::string, bool >	CanvasUpdatedList
	The type of list of canvas updated status.
typedef std::map< std::string, RefHistObject >	RefList
	The type of list of references.
enum	EModulePropFlags {   c_Input = 1 ,   c_Output = 2 ,   c_ParallelProcessingCertified = 4 ,   c_HistogramManager = 8 ,   c_InternalSerializer = 16 ,   c_TerminateInAllProcesses = 32 ,   c_DontCollectStatistics = 64 }
	Each module can be tagged with property flags, which indicate certain features of the module. More...
typedef ModuleCondition::EAfterConditionPath	EAfterConditionPath
	Forward the EAfterConditionPath definition from the ModuleCondition.

Public Member Functions
	DQMHistAnalysisKLMModule ()
	Constructor.
void	initialize () override final
	Initializer.
void	beginRun () override final
	Called when entering a new run.
void	event () override final
	This method is called for each event.
void	endRun () override final
	This method is called if the current run ends.
void	terminate () override final
	This method is called at the end of the event processing.
TH1 *	findHistInCanvas (const std::string &hname, TCanvas **canvas=nullptr)
	Find histogram in corresponding canvas.
TH1 *	getDelta (const std::string &dirname, const std::string &histname="", int n=0, bool onlyIfUpdated=true)
	Get Delta histogram.
int	registerEpicsPV (const std::string &pvname, const std::string &keyname="")
	EPICS related Functions.
int	registerExternalEpicsPV (const std::string &pvname, const std::string &keyname="")
	Register a PV with its name and a key name.
void	setEpicsPV (const std::string &keyname, double value)
	Write value to a EPICS PV.
void	setEpicsPV (const std::string &keyname, int value)
	Write value to a EPICS PV.
void	setEpicsPV (int index, double value)
	Write value to a EPICS PV.
void	setEpicsPV (int index, int value)
	Write value to a EPICS PV.
void	setEpicsStringPV (const std::string &keyname, const std::string &value)
	Write string to a EPICS PV.
void	setEpicsStringPV (int index, const std::string &value)
	Write string to a EPICS PV.
double	getEpicsPV (const std::string &keyname)
	Read value from a EPICS PV.
double	getEpicsPV (int index)
	Read value from a EPICS PV.
std::string	getEpicsStringPV (const std::string &keyname, bool &status)
	Read value from a EPICS PV.
std::string	getEpicsStringPV (int index, bool &status)
	Read value from a EPICS PV.
int	updateEpicsPVs (float timeout)
	Update all EPICS PV (flush to network)
chid	getEpicsPVChID (const std::string &keyname)
	Get EPICS PV Channel Id.
chid	getEpicsPVChID (int index)
	Get EPICS PV Channel Id.
bool	requestLimitsFromEpicsPVs (chid id, double &lowerAlarm, double &lowerWarn, double &upperWarn, double &upperAlarm)
	Get Alarm Limits from EPICS PV.
bool	requestLimitsFromEpicsPVs (const std::string &keyname, double &lowerAlarm, double &lowerWarn, double &upperWarn, double &upperAlarm)
	Get Alarm Limits from EPICS PV.
bool	requestLimitsFromEpicsPVs (int index, double &lowerAlarm, double &lowerWarn, double &upperWarn, double &upperAlarm)
	Get Alarm Limits from EPICS PV.
void	cleanupEpicsPVs (void)
	Unsubscribe from EPICS PVs on terminate.
void	CheckEpicsError (int state, const std::string &message, const std::string &name)
	check the return status and check PV in case of error
virtual std::vector< std::string >	getFileNames (bool outputFiles)
	Return a list of output filenames for this modules.
const std::string &	getName () const
	Returns the name of the module.
const std::string &	getType () const
	Returns the type of the module (i.e.
const std::string &	getPackage () const
	Returns the package this module is in.
const std::string &	getDescription () const
	Returns the description of the module.
void	setName (const std::string &name)
	Set the name of the module.
void	setPropertyFlags (unsigned int propertyFlags)
	Sets the flags for the module properties.
LogConfig &	getLogConfig ()
	Returns the log system configuration.
void	setLogConfig (const LogConfig &logConfig)
	Set the log system configuration.
void	setLogLevel (int logLevel)
	Configure the log level.
void	setDebugLevel (int debugLevel)
	Configure the debug messaging level.
void	setAbortLevel (int abortLevel)
	Configure the abort log level.
void	setLogInfo (int logLevel, unsigned int logInfo)
	Configure the printed log information for the given level.
void	if_value (const std::string &expression, const std::shared_ptr< Path > &path, EAfterConditionPath afterConditionPath=EAfterConditionPath::c_End)
	Add a condition to the module.
void	if_false (const std::shared_ptr< Path > &path, EAfterConditionPath afterConditionPath=EAfterConditionPath::c_End)
	A simplified version to add a condition to the module.
void	if_true (const std::shared_ptr< Path > &path, EAfterConditionPath afterConditionPath=EAfterConditionPath::c_End)
	A simplified version to set the condition of the module.
bool	hasCondition () const
	Returns true if at least one condition was set for the module.
const ModuleCondition *	getCondition () const
	Return a pointer to the first condition (or nullptr, if none was set)
const std::vector< ModuleCondition > &	getAllConditions () const
	Return all set conditions for this module.
bool	evalCondition () const
	If at least one condition was set, it is evaluated and true returned if at least one condition returns true.
std::shared_ptr< Path >	getConditionPath () const
	Returns the path of the last true condition (if there is at least one, else reaturn a null pointer).
Module::EAfterConditionPath	getAfterConditionPath () const
	What to do after the conditional path is finished.
std::vector< std::shared_ptr< Path > >	getAllConditionPaths () const
	Return all condition paths currently set (no matter if the condition is true or not).
bool	hasProperties (unsigned int propertyFlags) const
	Returns true if all specified property flags are available in this module.
bool	hasUnsetForcedParams () const
	Returns true and prints error message if the module has unset parameters which the user has to set in the steering file.
const ModuleParamList &	getParamList () const
	Return module param list.
template<typename T>
ModuleParam< T > &	getParam (const std::string &name) const
	Returns a reference to a parameter.
bool	hasReturnValue () const
	Return true if this module has a valid return value set.
int	getReturnValue () const
	Return the return value set by this module.
std::shared_ptr< PathElement >	clone () const override
	Create an independent copy of this module.
std::shared_ptr< boost::python::list >	getParamInfoListPython () const
	Returns a python list of all parameters.

Static Public Member Functions
static void	clearlist (void)
	Clear all static global lists.
static HistList &	getHistList ()
	Get the list of the histograms.
static MonObjList &	getMonObjList ()
	Get the list of MonitoringObjects.
static const DeltaList &	getDeltaList ()
	Get the list of the delta histograms.
static const CanvasUpdatedList &	getCanvasUpdatedList ()
	Get the list of the canvas update status.
static const std::string &	getRunType (void)
	Get the list of the reference histograms.
static int	getEventProcessed (void)
	Get the number of processed events.
static void	setRunType (const std::string &t)
	Set the Run Type.
static void	setEventProcessed (int e)
	Set the number of processed events.
static TCanvas *	findCanvas (TString cname)
	Find canvas by name.
static TH1 *	findHist (const std::string &dirname, const std::string &histname="", bool onlyIfUpdated=false)
	Find histogram.
static TH1 *	findRefHist (const std::string &dirname, const std::string &histname="", ERefScaling scaling=ERefScaling::c_RefScaleNone, const TH1 *hist=nullptr)
	Find reference histogram.
static TH1 *	scaleReference (ERefScaling scaling, const TH1 *hist, TH1 *ref)
	Using the original and reference, create scaled version.
static TH1 *	findHistInFile (TFile *file, const std::string &histname)
	Find histogram in specific TFile (e.g.
static MonitoringObject *	findMonitoringObject (const std::string &objName)
	Find MonitoringObject.
static double	getSigma68 (TH1 *h)
	Helper function to compute half of the central interval covering 68% of a distribution.
static bool	addHist (const std::string &dirname, const std::string &histname, TH1 *h)
	Add histogram.
static void	addRefHist (const std::string &dirname, TH1 *hist)
	Add reference histogram.
static MonitoringObject *	getMonitoringObject (const std::string &name)
	Get MonitoringObject with given name (new object is created if non-existing)
static void	clearCanvases (void)
	Clear content of all Canvases.
static void	initHistListBeforeEvent (void)
	Reset the list of histograms.
static void	clearHistList (void)
	Clears the list of histograms.
static void	clearRefList (void)
	Clears the list of ref histograms.
static void	resetDeltaList (void)
	Reset Delta.
static void	addDeltaPar (const std::string &dirname, const std::string &histname, HistDelta::EDeltaType t, int p, unsigned int a=1)
	Add Delta histogram parameters.
static bool	hasDeltaPar (const std::string &dirname, const std::string &histname)
	Check if Delta histogram parameters exist for histogram.
static void	UpdateCanvas (const std::string &name, bool updated=true)
	Mark canvas as updated (or not)
static void	UpdateCanvas (TCanvas *canvas, bool updated=true)
	Mark canvas as updated (or not)
static void	ExtractRunType (std::vector< TH1 * > &hs)
	Extract Run Type from histogram title, called from input module.
static void	ExtractNEvent (std::vector< TH1 * > &hs)
	Extract event processed from daq histogram, called from input module.
static void	setUseEpics (bool flag)
	Setter for EPICS usage.
static void	setUseEpicsReadOnly (bool flag)
	Setter EPICS flag in read only mode.
static bool	getUseEpics (void)
	Getter for EPICS usage.
static bool	getUseEpicsReadOnly (void)
	Getter EPICS flag in read only mode.
static std::string &	getPVPrefix (void)
	get global prefix for EPICS PVs
static void	setPVPrefix (const std::string &prefix)
	set global prefix for EPICS PVs
static EStatus	makeStatus (bool enough, bool warn_flag, bool error_flag)
	Helper function to judge the status for coloring and EPICS.
static void	colorizeCanvas (TCanvas *canvas, EStatus status)
	Helper function for Canvas colorization.
static EStatusColor	getStatusColor (EStatus status)
	Return color for canvas state.
static void	checkPVStatus (void)
	Check the status of all PVs and report if disconnected or not found.
static void	printPVStatus (chid pv, bool onlyError=true)
	check the status of a PVs and report if disconnected or not found
static void	CheckEpicsError (int state, const std::string &message, chid id)
	check the return status and check PV in case of error
static std::vector< std::string >	StringSplit (const std::string &s, const char delim)
	Helper function for string token split.
static void	exposePythonAPI ()
	Exposes methods of the Module class to Python.

Static Public Attributes
static DeltaList	s_deltaList
	The list of Delta Histograms and settings.

Protected Member Functions
virtual void	def_initialize ()
	Wrappers to make the methods without "def_" prefix callable from Python.
virtual void	def_beginRun ()
	Wrapper method for the virtual function beginRun() that has the implementation to be used in a call from Python.
virtual void	def_event ()
	Wrapper method for the virtual function event() that has the implementation to be used in a call from Python.
virtual void	def_endRun ()
	This method can receive that the current run ends as a call from the Python side.
virtual void	def_terminate ()
	Wrapper method for the virtual function terminate() that has the implementation to be used in a call from Python.
void	setDescription (const std::string &description)
	Sets the description of the module.
void	setType (const std::string &type)
	Set the module type.
template<typename T>
void	addParam (const std::string &name, T &paramVariable, const std::string &description, const T &defaultValue)
	Adds a new parameter to the module.
template<typename T>
void	addParam (const std::string &name, T &paramVariable, const std::string &description)
	Adds a new enforced parameter to the module.
void	setReturnValue (int value)
	Sets the return value for this module as integer.
void	setReturnValue (bool value)
	Sets the return value for this module as bool.
void	setParamList (const ModuleParamList &params)
	Replace existing parameter list.

Private Member Functions
double	getProcessedEvents ()
	Get number of processed events.
void	analyseChannelHitHistogram (int subdetector, int section, int sector, int index, TH1 *histogram, TH1 *delta, TCanvas *canvas, TLatex &latex)
	Analyse channel hit histogram.
void	processSpatial2DHitEndcapHistogram (uint16_t section, TH2F *histogram, TCanvas *canvas)
	Process spatial 2D hits histograms for endcap.
void	processTimeHistogram (const std::string &histName)
	Process histogram containing the hit times.
void	processPlaneHistogram (const std::string &histName, TLatex *latex, TH1 *histogram)
	Process histogram containing the number of hits in plane.
void	processPlaneHistogram (const std::string &histName, TLatex &latex)
	Process histogram containing the number of hits in plane.
void	updateCanvasStatus (TCanvas *canvas, const std::vector< KLMModuleNumber > &deadModules, TLatex *latex, int &message_counter, double xAlarm, double yAlarm)
	Helper function to update the canvas status based on dead modules.
void	processFEHistogram (TH1 *feHist, const std::string &histName, TCanvas *canvas)
	Process histogram containing the efficiencies.
void	fillMaskedChannelsHistogram (TH1 *histogram, TCanvas *canvas)
	Fill histogram containing masked channels per sector.
int	registerEpicsPVwithPrefix (const std::string &prefix, const std::string &pvname, const std::string &keyname="")
	Register a PV with its name and a key name.
std::list< ModulePtr >	getModules () const override
	no submodules, return empty list
std::string	getPathString () const override
	return the module name.
void	setParamPython (const std::string &name, const boost::python::object &pyObj)
	Implements a method for setting boost::python objects.
void	setParamPythonDict (const boost::python::dict &dictionary)
	Implements a method for reading the parameter values from a boost::python dictionary.

Static Private Member Functions
static void	deltaDrawer (TH1 *delta, TH1 *histogram, TCanvas *canvas)
	Scales and draws desired delta histogram for current canvas.

Private Attributes
double	m_ProcessedEvents
	Number of processed events.
int	m_ThresholdForMasked
	Threshold for masked channels.
int	m_ThresholdForHot
	Threshold for hot channels.
int	m_ThresholdForLog
	Threshold for log scale.
int	m_MinHitsForFlagging
	Minimal number of hits for flagging.
int	m_MessageThreshold
	Message Threshold for expert pots.
double	m_MinProcessedEventsForMessagesInput
	Input parameter for minimal number of processed events for error messages.
double	m_MinProcessedEventsForMessages
	Minimal number of processed events for error messages.
double	m_minEntries
	Minimal number of entries for delta histogram update.
double	m_minEvents
	Minimal number of processed events for delta histogram update.
std::string	m_histogramDirectoryName
	Name of histogram directory.
std::vector< uint16_t >	m_DeadBarrelModules
	Vector of dead barrel modules.
std::vector< uint16_t >	m_DeadEndcapModules
	Vector of dead endcap modules.
std::vector< uint16_t >	m_MaskedChannels
	Vector of masked channels.
TH1 *	m_MaskedChannelsHist = nullptr
	Histogram for masked channels per sector.
TCanvas *	m_c_masked_channels = nullptr
	Canvas for masked channels per sector summary.
TLine	m_2DHitsLine
	TLine for background region in 2d hits histograms.
TLine	m_PlaneLine
	TLine for boundary in plane histograms.
TText	m_PlaneText
	TText for names in plane histograms.
TH1 *	m_fe_bklm_ratio = nullptr
	Histogram for BKLM plane events fraction w/ FE.
TCanvas *	m_c_fe_bklm_ratio = nullptr
	Canvas for BKLM plane events fraction w/ FE.
TH1 *	m_fe_eklm_ratio = nullptr
	Histogram for EKLM plane events fraction w/ FE.
TCanvas *	m_c_fe_eklm_ratio = nullptr
	Canvas for EKLM plane events fraction w/ FE.
bool	m_IsNullRun
	Run type flag for null runs.
const KLMChannelArrayIndex *	m_ChannelArrayIndex
	KLM channel array index.
const KLMSectorArrayIndex *	m_SectorArrayIndex
	KLM sector array index.
const KLMElementNumbers *	m_ElementNumbers
	KLM element numbers.
const EKLMElementNumbers *	m_EklmElementNumbers
	EKLM element numbers.
DBObjPtr< KLMElectronicsMap >	m_ElectronicsMap
	Electronics map.
MonitoringObject *	m_monObj {}
	Monitoring object.
std::string	m_name
	The name of the module, saved as a string (user-modifiable)
std::string	m_type
	The type of the module, saved as a string.
std::string	m_package
	Package this module is found in (may be empty).
std::string	m_description
	The description of the module.
unsigned int	m_propertyFlags
	The properties of the module as bitwise or (with |) of EModulePropFlags.
LogConfig	m_logConfig
	The log system configuration of the module.
ModuleParamList	m_moduleParamList
	List storing and managing all parameter of the module.
bool	m_hasReturnValue
	True, if the return value is set.
int	m_returnValue
	The return value.
std::vector< ModuleCondition >	m_conditions
	Module condition, only non-null if set.

Static Private Attributes
static HistList	s_histList
	The list of Histograms.
static MonObjList	s_monObjList
	The list of MonitoringObjects.
static CanvasUpdatedList	s_canvasUpdatedList
	The list of canvas updated status.
static RefList	s_refList
	The list of references.
static int	s_eventProcessed = 0
	Number of Events processed to fill histograms.
static std::string	s_runType = ""
	The Run type.
static bool	m_useEpics = false
	Flag if to use EPICS do not set by yourself, use EpicsEnable module to set.
static bool	m_epicsReadOnly
	Flag if to use EPICS in ReadOnly mode (for reading limits) do not set by yourself, use EpicsEnable module to set.
static std::string	m_PVPrefix = "TEST:"
	The Prefix for EPICS PVs.

Detailed Description

Analysis of KLM DQM histograms.

Definition at line 43 of file DQMHistAnalysisKLM.h.

Member Typedef Documentation

CanvasUpdatedList

typedef std::map<std::string, bool> CanvasUpdatedList

inherited

The type of list of canvas updated status.

Definition at line 95 of file DQMHistAnalysis.h.

DeltaList

typedef std::map<std::string, HistDelta> DeltaList

inherited

The type of list of delta settings and histograms.

Definition at line 90 of file DQMHistAnalysis.h.

EAfterConditionPath

typedef ModuleCondition::EAfterConditionPath EAfterConditionPath

inherited

Forward the EAfterConditionPath definition from the ModuleCondition.

Definition at line 88 of file Module.h.

HistList

typedef std::map<std::string, HistObject> HistList

inherited

The type of list of histograms.

Definition at line 81 of file DQMHistAnalysis.h.

MonObjList

typedef std::map<std::string, MonitoringObject> MonObjList

inherited

The type of list of MonitoringObjects.

Definition at line 85 of file DQMHistAnalysis.h.

RefList

typedef std::map<std::string, RefHistObject> RefList

inherited

The type of list of references.

Definition at line 100 of file DQMHistAnalysis.h.

Member Enumeration Documentation

EModulePropFlags

enum EModulePropFlags

inherited

Each module can be tagged with property flags, which indicate certain features of the module.

Enumerator
c_Input	This module is an input module (reads data).
c_Output	This module is an output module (writes data).
c_ParallelProcessingCertified	This module can be run in parallel processing mode safely (All I/O must be done through the data store, in particular, the module must not write any files.)
c_HistogramManager	This module is used to manage histograms accumulated by other modules.
c_InternalSerializer	This module is an internal serializer/deserializer for parallel processing.
c_TerminateInAllProcesses	When using parallel processing, call this module's terminate() function in all processes(). This will also ensure that there is exactly one process (single-core if no parallel modules found) or at least one input, one main and one output process.
c_DontCollectStatistics	No statistics is collected for this module.

Definition at line 77 of file Module.h.

                          {
      c_Input                       = 1,  
      c_Output                      = 2,  
      c_ParallelProcessingCertified = 4,  
      c_HistogramManager            = 8, 
      c_InternalSerializer          = 16,  
      c_TerminateInAllProcesses     = 32,  
      c_DontCollectStatistics       = 64,  
    };

ERefScaling

enum ERefScaling

inherited

Reference plot scaling type.

Enumerator
c_RefScaleNone	no scaling
c_RefScaleEntries	to number of entries (integral)
c_RefScaleMax	to maximum (bin entry)

Definition at line 72 of file DQMHistAnalysis.h.

                     {
      c_RefScaleNone = 0, 
      c_RefScaleEntries = 1, 
      c_RefScaleMax = 2 
    };

EStatus

enum EStatus

inherited

Status flag of histogram/canvas.

Enumerator
c_StatusTooFew	Not enough entries/event to judge.
c_StatusDefault	default for non-coloring
c_StatusGood	Analysis result: Good.
c_StatusWarning	Analysis result: Warning, there may be minor issues.
c_StatusError	Analysis result: Severe issue found.

Definition at line 50 of file DQMHistAnalysis.h.

                 {
      c_StatusTooFew = 0, 
      c_StatusDefault = 1, 
      c_StatusGood = 2, 
      c_StatusWarning = 3, 
      c_StatusError = 4 
    };

EStatusColor

enum EStatusColor

inherited

Status colors of histogram/canvas (corresponding to status)

Enumerator
c_ColorTooFew	Not enough entries/event to judge.
c_ColorDefault	default for non-coloring
c_ColorGood	Analysis result: Good.
c_ColorWarning	Analysis result: Warning, there may be minor issues.
c_ColorError	Analysis result: Severe issue found.

Definition at line 61 of file DQMHistAnalysis.h.

                      {
      c_ColorTooFew = kGray, 
      c_ColorDefault = kWhite, 
      c_ColorGood = kGreen, 
      c_ColorWarning = kYellow, 
      c_ColorError = kRed 
    };

Constructor & Destructor Documentation

DQMHistAnalysisKLMModule()

DQMHistAnalysisKLMModule

(

)

Constructor.

Definition at line 27 of file DQMHistAnalysisKLM.cc.

  : DQMHistAnalysisModule(),
    m_ProcessedEvents{0},
    m_IsNullRun{false},
    m_ChannelArrayIndex{&(KLMChannelArrayIndex::Instance())},
    m_SectorArrayIndex{&(KLMSectorArrayIndex::Instance())},
    m_ElementNumbers{&(KLMElementNumbers::Instance())},
    m_EklmElementNumbers{&(EKLMElementNumbers::Instance())}
{
  setDescription("Module used to analyze KLM DQM histograms.");
  addParam("ThresholdForMasked", m_ThresholdForMasked,
           "Threshold X for masked channels: if a channel has an occupancy X times larger than the average, it will be masked.", 100);
  addParam("ThresholdForHot", m_ThresholdForHot,
           "Threshold Y for hot channels: if a channel has an occupancy Y times larger than the average, it will be marked as hot (but not masked).",
           10);
  addParam("ThresholdForLog", m_ThresholdForLog,
           "Threshold Z for log scale view: if a channel has an occupancy Z times larger than the average, canvas shifts to log scale.",
           20);
  addParam("MinHitsForFlagging", m_MinHitsForFlagging, "Minimal number of hits in a channel required to flag it as 'Masked' or 'Hot'",
           50);
  addParam("MinProcessedEventsForMessages", m_MinProcessedEventsForMessagesInput,
           "Minimal number of processed events required to print error messages", 10000.);
  addParam("MinEntries", m_minEntries,
           "Minimal number of entries for delta histogram updates", 30000.);
  addParam("MinEvents", m_minEvents,
           "Minimal number of processed events for delta histogram updates", 30000.);
  addParam("MessageThreshold", m_MessageThreshold,
           "Max number of messages to show up in channel occupancy plots", 12);
  addParam("HistogramDirectoryName", m_histogramDirectoryName, "Name of histogram directory", std::string("KLM"));
 
  m_MinProcessedEventsForMessages = m_MinProcessedEventsForMessagesInput;
  m_2DHitsLine.SetLineColor(kRed);
  m_2DHitsLine.SetLineWidth(3);
  m_2DHitsLine.SetLineStyle(2); // dashed
  m_PlaneLine.SetLineColor(kMagenta);
  m_PlaneLine.SetLineWidth(1);
  m_PlaneLine.SetLineStyle(2); // dashed
  m_PlaneText.SetTextAlign(22); // centered, middle
  m_PlaneText.SetTextColor(kMagenta);
  m_PlaneText.SetTextFont(42); // Helvetica regular
  m_PlaneText.SetTextSize(0.02); // 2% of TPad's full height
}

Member Function Documentation

addDeltaPar()

void addDeltaPar ( const std::string & dirname, const std::string & histname, HistDelta::EDeltaType t, int p, unsigned int a = 1 )

staticinherited

Add Delta histogram parameters.

Parameters

dirname	directory
histname	name of histogram
t	type of delta histogramming
p	numerical parameter depending on type, e.g. number of entries
a	amount of histograms in the past

Definition at line 94 of file DQMHistAnalysis.cc.

{
  std::string fullname;
  if (dirname.size() > 0) {
    fullname = dirname + "/" + histname;
  } else {
    fullname = histname;
  }
  s_deltaList[fullname].set(t, p, a);
}

addHist()

bool addHist ( const std::string & dirname, const std::string & histname, TH1 * h )

staticinherited

Add histogram.

Parameters

dirname	The name of the directory.
histname	The name of the histogram.
h	The TH1 pointer for the histogram.

Returns

histogram was updated flag

Definition at line 59 of file DQMHistAnalysis.cc.

{
  std::string fullname;
  if (dirname.size() > 0) {
    fullname = dirname + "/" + histname;
  } else {
    fullname = histname;
  }
 
  if (s_histList[fullname].update(h)) {
    // only if histogram changed, check if delta histogram update needed
    auto it = s_deltaList.find(fullname);
    if (it != s_deltaList.end()) {
      B2DEBUG(20, "Found Delta" << fullname);
      it->second.update(h); // update
    }
    return true; // histogram changed
  }
 
  return false; // histogram didn't change
}

addRefHist()

void addRefHist ( const std::string & dirname, TH1 * hist )

staticinherited

Add reference histogram.

Parameters

dirname	The name of the directory.
hist	The TH1 pointer for the reference.

Definition at line 81 of file DQMHistAnalysis.cc.

{
  string histname = hist->GetName();
  std::string name = dirname + "/" + histname;
  auto& n = s_refList[name];
  n.m_orghist_name = name;
  n.m_refhist_name = "ref/" + name;
  hist->SetName((n.m_refhist_name).c_str());
  hist->SetDirectory(0);
  n.setRefHist(hist); // transfer ownership!
  n.setRefCopy(nullptr);
}

analyseChannelHitHistogram()

void analyseChannelHitHistogram ( int subdetector, int section, int sector, int index, TH1 * histogram, TH1 * delta, TCanvas * canvas, TLatex & latex )

private

Analyse channel hit histogram.

Parameters

[in]	subdetector	Subdetector (E or B-KLM).
[in]	section	Section (Forward or Backward).
[in]	sector	Sector.
[in]	index	Histogram Index (The j'th histogram of a given sector).
[in]	histogram	Histogram.
[in]	delta	Delta Histogram.
[in]	canvas	Canvas.
[out]	latex	TLatex to draw messages.

Definition at line 219 of file DQMHistAnalysisKLM.cc.

{
  double x = 0.15;
  double y = 0.85;
  int i, n;
  std::map<KLMModuleNumber, double> moduleHitMap;
  std::map<KLMModuleNumber, double>::iterator it;
  double average = 0;
  int channelSubdetector, channelSection, channelSector;
  int layer, plane, strip;
  std::string str;
  canvas->Clear();
  canvas->cd();
  histogram->SetStats(false);
  canvas->SetLogy(0); //initialize to start without logscale
  histogram->Draw();
  deltaDrawer(delta, histogram, canvas); //draw normalized delta on top
  n = histogram->GetXaxis()->GetNbins();
 
  /* call reference histograms from base class*/
  TH1* ref_histogram = findRefHist(histogram->GetName(), "", ERefScaling::c_RefScaleEntries, histogram);
  if (ref_histogram) {ref_histogram->Draw("hist,same");}
  float ref_average = 0;
 
  if (ref_histogram != nullptr) {
    for (i = 1; i <= n; i++) {
      double nHitsPerModuleRef = ref_histogram->GetBinContent(i);
      ref_average = ref_average + nHitsPerModuleRef;
    }
  }
 
  for (i = 1; i <= n; i++) {
    KLMChannelNumber channelIndex = std::round(histogram->GetBinCenter(i));
    KLMChannelNumber channelNumber =
      m_ChannelArrayIndex->getNumber(channelIndex);
    double nHitsPerModule = histogram->GetBinContent(i);
    average = average + nHitsPerModule;
    m_ElementNumbers->channelNumberToElementNumbers(
      channelNumber, &channelSubdetector, &channelSection, &channelSector,
      &layer, &plane, &strip);
    if ((channelSubdetector != subdetector) ||
        (channelSection != section) ||
        (channelSector != sector))
      B2FATAL("Inconsistent element numbers.");
    KLMModuleNumber module = m_ElementNumbers->moduleNumber(
                               subdetector, section, sector, layer);
    it = moduleHitMap.find(module);
    if (it == moduleHitMap.end()) {
      moduleHitMap.insert(std::pair<KLMModuleNumber, double>(
                            module, nHitsPerModule));
    } else {
      it->second += nHitsPerModule;
    }
  }
  unsigned int activeModuleChannels = 0;
  int message_counter = 0;
  for (it = moduleHitMap.begin(); it != moduleHitMap.end(); ++it) {
    KLMModuleNumber moduleNumber = it->first;
    if (it->second != 0) {
      activeModuleChannels += m_ElementNumbers->getNChannelsModule(moduleNumber);
      continue;
    }
    m_ElementNumbers->moduleNumberToElementNumbers(
      moduleNumber, &channelSubdetector, &channelSection, &channelSector, &layer);
    /* Channel with plane = 1, strip = 1 exists for any BKLM or EKLM module. */
    KLMChannelNumber channel =
      m_ElementNumbers->channelNumber(
        channelSubdetector, channelSection, channelSector, layer, 1, 1);
    const KLMElectronicsChannel* electronicsChannel =
      m_ElectronicsMap->getElectronicsChannel(channel);
    if (electronicsChannel == nullptr)
      B2FATAL("Incomplete KLM electronics map.");
    str = "No data from lane " + std::to_string(electronicsChannel->getLane());
    latex.DrawLatexNDC(x, y, str.c_str());
    y -= 0.05;
    message_counter++;
    /* Store the module number, used later in processPlaneHistogram
     * to color the canvas with red and to raise up an alarm. */
    if (channelSubdetector == KLMElementNumbers::c_BKLM) {
      std::vector<KLMModuleNumber>::iterator ite =
        std::find(m_DeadBarrelModules.begin(),
                  m_DeadBarrelModules.end(),
                  moduleNumber);
      if (ite == m_DeadBarrelModules.end())
        m_DeadBarrelModules.push_back(moduleNumber);
    } else {
      std::vector<KLMModuleNumber>::iterator ite = std::find(m_DeadEndcapModules.begin(),
                                                             m_DeadEndcapModules.end(),
                                                             moduleNumber);
      if (ite == m_DeadEndcapModules.end())
        m_DeadEndcapModules.push_back(moduleNumber);
    }
  }
  if (activeModuleChannels == 0)
    return;
  average /= activeModuleChannels;
  ref_average /= activeModuleChannels;
 
  for (i = 1; i <= n; ++i) {
    KLMChannelNumber channelIndex = std::round(histogram->GetBinCenter(i));
    KLMChannelNumber channelNumber =
      m_ChannelArrayIndex->getNumber(channelIndex);
    double nHits = histogram->GetBinContent(i);
    m_ElementNumbers->channelNumberToElementNumbers(
      channelNumber, &channelSubdetector, &channelSection, &channelSector,
      &layer, &plane, &strip);
    std::string channelStatus = "Normal";
    if ((nHits > average * m_ThresholdForMasked) && (nHits > m_MinHitsForFlagging)) {
      channelStatus = "Masked";
      std::vector<KLMModuleNumber>::iterator ite =
        std::find(m_MaskedChannels.begin(),
                  m_MaskedChannels.end(),
                  channelNumber);
      if (ite == m_MaskedChannels.end())
        m_MaskedChannels.push_back(channelNumber);
      B2DEBUG(20, "KLM@MaskMe " << channelNumber);
    } else if ((nHits > average * m_ThresholdForHot) && (nHits > m_MinHitsForFlagging)) {
      channelStatus = "Hot";
    }
    if (channelStatus != "Normal") {
      const KLMElectronicsChannel* electronicsChannel =
        m_ElectronicsMap->getElectronicsChannel(channelNumber);
      if (electronicsChannel == nullptr)
        B2FATAL("Incomplete BKLM electronics map.");
      if (channelStatus == "Masked") {
        histogram->SetBinContent(i, 0);
        if (delta != nullptr)
          delta->SetBinContent(i, 0);
      }
      str = channelStatus + " channel: ";
      // lane, axis, channel
      str += ("L" + std::to_string(electronicsChannel->getLane()) +
              " A" + std::to_string(electronicsChannel->getAxis()) +
              " Ch" + std::to_string(electronicsChannel->getChannel()));
      message_counter++;
      if (message_counter <= m_MessageThreshold) {
        latex.DrawLatexNDC(x, y, str.c_str());
        y -= 0.05;
      }
    }
  }
  if (message_counter > m_MessageThreshold) {
    std::string verbose_message = " more messages";
    verbose_message = std::to_string(message_counter - m_MessageThreshold) + verbose_message;
    latex.DrawLatexNDC(x, y, verbose_message.c_str());
  }
 
 
  // for hot/masked channels, log scale plots (reference and main)
  if (histogram->GetMaximum()*n > histogram->Integral()*m_ThresholdForLog && average * activeModuleChannels > m_MinHitsForFlagging) {
    histogram->SetMinimum(1);
    canvas->SetLogy();
  } else if (ref_histogram != nullptr) {
    if (ref_histogram->GetMaximum()*n > ref_histogram->Integral()*m_ThresholdForLog
        && ref_average * activeModuleChannels > m_MinHitsForFlagging) {
      histogram->SetMinimum(1);
      canvas->SetLogy();
    } else
      canvas->SetLogy(0);
  } else
    canvas->SetLogy(0);
 
  canvas->Modified();
  canvas->Update();
 
  /* Drawing dividing lines */
  int divisions;
  int bin = 1;
  double xLine;
  // drawing lines for BKLM sectors
  if (subdetector == 1) {
    int shift;
    if (index == 0) {
      divisions = 7;
      shift = 1;
    } else {
      divisions = BKLMElementNumbers::getMaximalSectorNumber();
      shift = BKLMElementNumbers::getMaximalSectorNumber();
    }
    for (int k = 0; k < divisions; k++) {
      xLine = (histogram->GetXaxis()->GetBinLowEdge(bin) - canvas->GetX1()) / (canvas->GetX2() - canvas->GetX1());
      m_PlaneLine.DrawLineNDC(xLine, 0.1, xLine, 0.9);
      bin += BKLMElementNumbers::getNStrips(section, sector, k + shift, 0)
             + BKLMElementNumbers::getNStrips(section, sector, k + shift, 1);
    }
  } else { // drawing lines for EKLM sectors
    if ((section == 2) && (index == 0 || index == 1))
      divisions = 5;
    else
      divisions = 4;
    for (int k = 0; k < divisions; k++) {
      xLine = (histogram->GetXaxis()->GetBinLowEdge(bin) - canvas->GetX1()) / (canvas->GetX2() - canvas->GetX1());
      m_PlaneLine.DrawLineNDC(xLine, 0.1, xLine, 0.9);
      bin += EKLMElementNumbers::getNStripsSector();
    }
  }
  canvas->Modified();
  canvas->Update();
 
}

beginRun()

void beginRun ( void )

finaloverridevirtual

Called when entering a new run.

Reimplemented from Module.

Definition at line 157 of file DQMHistAnalysisKLM.cc.

{
  if (!m_ElectronicsMap.isValid())
    B2FATAL("No KLM electronics map.");
  m_MinProcessedEventsForMessages = m_MinProcessedEventsForMessagesInput;
  m_ProcessedEvents = 0.;
  m_DeadBarrelModules.clear();
  m_DeadEndcapModules.clear();
  m_MaskedChannels.clear();
 
  m_IsNullRun = (getRunType() == "null");
}

CheckEpicsError() [1/2]

void CheckEpicsError ( int state, const std::string & message, chid id = nullptr )

staticinherited

check the return status and check PV in case of error

Parameters

state	return state of epics function
message	message to print out
id	the chid of the affected PV

Definition at line 792 of file DQMHistAnalysis.cc.

{
  if (state != ECA_NORMAL) {
    std::string name;
    if (id) name =  ca_name(id);
    B2WARNING(message << ": " << name);
    printPVStatus(id, false);
  }
}

CheckEpicsError() [2/2]

void CheckEpicsError ( int state, const std::string & message, const std::string & name )

inherited

check the return status and check PV in case of error

Parameters

state	return state of epics function
message	message to print out
name	the (key)name of the affected PV

Definition at line 784 of file DQMHistAnalysis.cc.

{
  if (state != ECA_NORMAL) {
    B2WARNING(message << ": " << name);
    printPVStatus(m_epicsNameToChID[name], false);
  }
}

checkPVStatus()

void checkPVStatus ( void )

staticinherited

Check the status of all PVs and report if disconnected or not found.

Definition at line 748 of file DQMHistAnalysis.cc.

{
  B2INFO("Check PV Connections");
 
  for (const auto& it : m_epicsChID) {
    printPVStatus(it);
  }
  B2INFO("Check PVs done");
}

cleanupEpicsPVs()

void cleanupEpicsPVs ( void )

inherited

Unsubscribe from EPICS PVs on terminate.

Definition at line 629 of file DQMHistAnalysis.cc.

{
  // this should be called in terminate function of analysis modules
#ifdef _BELLE2_EPICS
  if (getUseEpics()) {
    for (auto& it : m_epicsChID) CheckEpicsError(ca_clear_channel(it), "ca_clear_channel failure", it);
    updateEpicsPVs(5.0);
    // Make sure we clean up both afterwards!
    m_epicsChID.clear();
    m_epicsNameToChID.clear();
  }
#endif
}

clearCanvases()

void clearCanvases ( void )

staticinherited

Clear content of all Canvases.

Definition at line 294 of file DQMHistAnalysis.cc.

{
  TIter nextckey(gROOT->GetListOfCanvases());
  TObject* cobj = nullptr;
 
  while ((cobj = dynamic_cast<TObject*>(nextckey()))) {
    if (cobj->IsA()->InheritsFrom("TCanvas")) {
      TCanvas* cnv = dynamic_cast<TCanvas*>(cobj);
      cnv->Clear();
      DQMHistAnalysisModule::colorizeCanvas(cnv, c_StatusDefault);
    }
  }
}

clearHistList()

void clearHistList ( void )

staticinherited

Clears the list of histograms.

Definition at line 322 of file DQMHistAnalysis.cc.

{
  s_histList.clear();
}

clearlist()

void clearlist ( void )

staticinherited

Clear all static global lists.

Definition at line 50 of file DQMHistAnalysis.cc.

{
  s_histList.clear();
  s_refList.clear();
// s_monObjList;
  s_deltaList.clear();
  s_canvasUpdatedList.clear();
}

clearRefList()

void clearRefList ( void )

staticinherited

Clears the list of ref histograms.

Definition at line 327 of file DQMHistAnalysis.cc.

{
  s_refList.clear();
}

clone()

std::shared_ptr< PathElement > clone ( ) const

overridevirtualinherited

Create an independent copy of this module.

Note that parameters are shared, so changing them on a cloned module will also affect the original module.

Implements PathElement.

Definition at line 179 of file Module.cc.

{
  ModulePtr newModule = ModuleManager::Instance().registerModule(getType());
  newModule->m_moduleParamList.setParameters(getParamList());
  newModule->setName(getName());
  newModule->m_package = m_package;
  newModule->m_propertyFlags = m_propertyFlags;
  newModule->m_logConfig = m_logConfig;
  newModule->m_conditions = m_conditions;
 
  return newModule;
}

colorizeCanvas()

void colorizeCanvas ( TCanvas * canvas, EStatus status )

staticinherited

Helper function for Canvas colorization.

Parameters

canvas	Canvas to change
status	status to color

Definition at line 735 of file DQMHistAnalysis.cc.

{
  if (!canvas) return;
  auto color = DQMHistAnalysisModule::getStatusColor(stat);
 
  canvas->Pad()->SetFillColor(color);
 
  canvas->Pad()->SetFrameFillColor(10); // White (kWhite is not used since it results in transparent!)
  canvas->Pad()->SetFrameFillStyle(1001);// White
  canvas->Pad()->Modified();
  canvas->Pad()->Update();
}

def_beginRun()

virtual void def_beginRun ( )

inlineprotectedvirtualinherited

Wrapper method for the virtual function beginRun() that has the implementation to be used in a call from Python.

Reimplemented in PyModule.

Definition at line 425 of file Module.h.

{ beginRun(); }

def_endRun()

virtual void def_endRun ( )

inlineprotectedvirtualinherited

This method can receive that the current run ends as a call from the Python side.

For regular C++-Modules that forwards the call to the regular endRun() method.

Reimplemented in PyModule.

Definition at line 438 of file Module.h.

{ endRun(); }

def_event()

virtual void def_event ( )

inlineprotectedvirtualinherited

Wrapper method for the virtual function event() that has the implementation to be used in a call from Python.

Reimplemented in PyModule.

Definition at line 431 of file Module.h.

{ event(); }

def_initialize()

virtual void def_initialize ( )

inlineprotectedvirtualinherited

Wrappers to make the methods without "def_" prefix callable from Python.

Overridden in PyModule. Wrapper method for the virtual function initialize() that has the implementation to be used in a call from Python.

Reimplemented in PyModule.

Definition at line 419 of file Module.h.

{ initialize(); }

def_terminate()

virtual void def_terminate ( )

inlineprotectedvirtualinherited

Wrapper method for the virtual function terminate() that has the implementation to be used in a call from Python.

Reimplemented in PyModule.

Definition at line 444 of file Module.h.

{ terminate(); }

deltaDrawer()

void deltaDrawer ( TH1 * delta, TH1 * histogram, TCanvas * canvas )

staticprivate

Scales and draws desired delta histogram for current canvas.

Parameters

[in]	delta	Delta histogram.
[in]	histogram	Base histogram (for normalization).
[in]	canvas	Canvas with delta histogram.

Definition at line 205 of file DQMHistAnalysisKLM.cc.

{
  if (delta != nullptr) {
    Double_t scale = (Double_t) histogram->Integral(); //want delta and histo to have same norm
 
    // delta != nullptr should take care of whether update condition is met.
    delta->SetLineColor(kBlackBody); //random choice of not green or blue
    delta->SetLineStyle(4);
    delta->DrawNormalized("SAME", scale); //normalize delta to histo
    canvas->Modified();
    canvas->Update();
  }
}

endRun()

void endRun ( void )

finaloverridevirtual

This method is called if the current run ends.

Reimplemented from Module.

Definition at line 170 of file DQMHistAnalysisKLM.cc.

{
  int hist_max_bin; double max_position;
  TH1* time_rpc = findHist(m_histogramDirectoryName + "/time_rpc");
  if (time_rpc) {
    hist_max_bin = time_rpc->GetMaximumBin();
    max_position = time_rpc->GetXaxis()->GetBinCenter(hist_max_bin);
    m_monObj->setVariable("RPC_Time_Peak", max_position);
  }
 
  TH1* time_scint_bklm = findHist(m_histogramDirectoryName + "/time_scintillator_bklm");
  if (time_scint_bklm) {
    hist_max_bin = time_scint_bklm->GetMaximumBin();
    max_position = time_scint_bklm->GetXaxis()->GetBinCenter(hist_max_bin);
    m_monObj->setVariable("BKLM_Scint_Time_Peak", max_position);
  }
 
  TH1* time_scint_eklm = findHist(m_histogramDirectoryName + "/time_scintillator_eklm");
  if (time_scint_eklm) {
    hist_max_bin = time_scint_eklm->GetMaximumBin();
    max_position = time_scint_eklm->GetXaxis()->GetBinCenter(hist_max_bin);
    m_monObj->setVariable("EKLM_Scint_Time_Peak", max_position);
  }
}

evalCondition()

bool evalCondition ( ) const

inherited

If at least one condition was set, it is evaluated and true returned if at least one condition returns true.

If no condition or result value was defined, the method returns false. Otherwise, the condition is evaluated and true returned, if at least one condition returns true. To speed up the evaluation, the condition strings were already parsed in the method if_value().

Returns

True if at least one condition and return value exists and at least one condition expression was evaluated to true.

Definition at line 96 of file Module.cc.

{
  if (m_conditions.empty()) return false;
 
  //okay, a condition was set for this Module...
  if (!m_hasReturnValue) {
    B2FATAL("A condition was set for '" << getName() << "', but the module did not set a return value!");
  }
 
  for (const auto& condition : m_conditions) {
    if (condition.evaluate(m_returnValue)) {
      return true;
    }
  }
  return false;
}

event()

void event ( void )

finaloverridevirtual

This method is called for each event.

Reimplemented from Module.

Definition at line 714 of file DQMHistAnalysisKLM.cc.

{
  /* If KLM is not included, stop here and return. */
  TH1* daqInclusion = findHist(m_histogramDirectoryName + "/daq_inclusion");
  if (not(daqInclusion == nullptr)) {
    int isKlmIncluded = daqInclusion->GetBinContent(daqInclusion->GetXaxis()->FindBin("Yes"));
    if (isKlmIncluded == 0)
      return;
  }
  /* Make sure that the vectors are cleared at each DQM refresh. */
  m_DeadBarrelModules.clear();
  m_DeadEndcapModules.clear();
  m_MaskedChannels.clear();
  m_ProcessedEvents = getProcessedEvents();
  std::string str, histogramName, canvasName;
  TLatex latex;
  latex.SetTextColor(kRed);
  latex.SetTextAlign(11);
  KLMChannelIndex klmIndex(KLMChannelIndex::c_IndexLevelSector);
  // gathering relevant info for analyseChannelHitHistogram
  for (const KLMChannelIndex& klmSector : klmIndex) {
    int nHistograms;
    if (klmSector.getSubdetector() == KLMElementNumbers::c_BKLM)
      nHistograms = 2;
    else
      nHistograms = 3;
    for (int j = 0; j < nHistograms; j++) {
      str = "strip_hits_subdetector_" +
            std::to_string(klmSector.getSubdetector()) +
            "_section_" + std::to_string(klmSector.getSection()) +
            "_sector_" + std::to_string(klmSector.getSector()) +
            "_" + std::to_string(j);
      histogramName = m_histogramDirectoryName + "/" + str;
      canvasName = m_histogramDirectoryName + "/c_" + str;
      TH1* histogram = findHist(histogramName);
      //get delta histogram (should we work with a clone instead?)
      auto delta = getDelta("", histogramName);
 
      if (histogram == nullptr) {
        B2WARNING("KLM DQM histogram " << histogramName << " is not found.");
        continue;
      }
      TCanvas* canvas = findCanvas(canvasName);
      if (canvas == nullptr) {
        B2WARNING("KLM DQM histogram canvas " << canvasName << " is not found.");
        continue;
      }
      // Add this canvas that it is time to update
      // not sure if this is interfering with the generation of some features
      // after testing, switch condition back to delta != nullptr || histogram != nullptr
      UpdateCanvas(canvas->GetName(), true);
      analyseChannelHitHistogram(
        klmSector.getSubdetector(), klmSector.getSection(),
        klmSector.getSector(), j, histogram, delta, canvas, latex);
 
    }
  }
  /* Temporary change the color palette. */
  gStyle->SetPalette(kLightTemperature);
  klmIndex.setIndexLevel(KLMChannelIndex::c_IndexLevelSection);
  for (const KLMChannelIndex& klmSection : klmIndex) {
    KLMSubdetectorNumber subdetector = klmSection.getSubdetector();
    if (subdetector == KLMElementNumbers::c_EKLM) {
      KLMSubdetectorNumber section = klmSection.getSection();
      int maximalLayerNumber = m_EklmElementNumbers->getMaximalDetectorLayerNumber(section);
      for (int j = 1; j <= maximalLayerNumber; ++j) {
        str = "spatial_2d_hits_subdetector_" + std::to_string(subdetector) +
              "_section_" + std::to_string(section) +
              "_layer_" + std::to_string(j);
        histogramName = m_histogramDirectoryName + "/" + str;
        canvasName = m_histogramDirectoryName + "/c_" + str;
        TH2F* histogram = static_cast<TH2F*>(findHist(histogramName));
        if (histogram == nullptr) {
          B2ERROR("KLM DQM histogram " << histogramName << " is not found.");
          continue;
        }
        TCanvas* canvas = findCanvas(canvasName);
        if (canvas == nullptr) {
          B2ERROR("KLM DQM histogram canvas " << canvasName << " is not found.");
          continue;
        }
        processSpatial2DHitEndcapHistogram(section, histogram, canvas);
      }
    }
  }
 
  /* Reset the color palette to the default one. */
  gStyle->SetPalette(kBird);
  fillMaskedChannelsHistogram(m_MaskedChannelsHist, m_c_masked_channels);
  latex.SetTextColor(kBlue);
  processPlaneHistogram("plane_bklm_phi", latex);
  processPlaneHistogram("plane_bklm_z", latex);
  processPlaneHistogram("plane_eklm", latex);
 
  processFEHistogram(m_fe_bklm_ratio, "feStatus_bklm_scintillator_layers", m_c_fe_bklm_ratio);
  processFEHistogram(m_fe_eklm_ratio, "feStatus_eklm_plane", m_c_fe_eklm_ratio);
  processPlaneHistogram("fe_bklm_ratio", nullptr, m_fe_bklm_ratio);
  processPlaneHistogram("fe_eklm_ratio", nullptr, m_fe_eklm_ratio);
 
 
  processTimeHistogram("time_rpc");
  processTimeHistogram("time_scintillator_bklm");
  processTimeHistogram("time_scintillator_eklm");
 
  B2DEBUG(20, "Updating EPICS PVs for DQMHistAnalysisKLM");
  // only update PVs if there's enough statistics and datasize != 0
  // Check if it's a null run, if so, don't update EPICS PVs
  if (m_IsNullRun) {
    B2INFO("DQMHistAnalysisKLM: Null run detected. No PV Update.");
    return;
  }
  auto* daqDataSize = findHist("DAQ/KLMDataSize");
  double meanDAQDataSize = 0.;
  if (daqDataSize != nullptr) {
    meanDAQDataSize = daqDataSize->GetMean();
  } else
    B2WARNING("DQMHistAnalysisKLM: Cannot find KLMDataSize");
  if ((daqDataSize != nullptr) and (meanDAQDataSize != 0.)) {
    if (m_ProcessedEvents > m_MinProcessedEventsForMessages) {
      setEpicsPV("MaskedChannels", (double)m_MaskedChannels.size());
      setEpicsPV("DeadBarrelModules", (double)m_DeadBarrelModules.size());
      setEpicsPV("DeadEndcapModules", (double)m_DeadEndcapModules.size());
      B2DEBUG(20, "DQMHistAnalysisKLM: MaskedChannels " << m_MaskedChannels.size());
      B2DEBUG(20, "DQMHistAnalysisKLM: DeadBarrelModules " << m_DeadBarrelModules.size());
      B2DEBUG(20, "DQMHistAnalysisKLM: DeadEndcapModules " << m_DeadEndcapModules.size());
    }
  } else
    B2INFO("DQMHistAnalysisKLM: KLM Not included. No PV Update. ");
 
 
}

exposePythonAPI()

void exposePythonAPI ( )

staticinherited

Exposes methods of the Module class to Python.

Definition at line 325 of file Module.cc.

{
  // to avoid confusion between std::arg and boost::python::arg we want a shorthand namespace as well
  namespace bp = boost::python;
 
  docstring_options options(true, true, false); //userdef, py sigs, c++ sigs
 
  void (Module::*setReturnValueInt)(int) = &Module::setReturnValue;
 
  enum_<Module::EAfterConditionPath>("AfterConditionPath",
                                     R"(Determines execution behaviour after a conditional path has been executed:
 
.. attribute:: END
 
  End processing of this path after the conditional path. (this is the default for if_value() etc.)
 
.. attribute:: CONTINUE
 
  After the conditional path, resume execution after this module.)")
  .value("END", Module::EAfterConditionPath::c_End)
  .value("CONTINUE", Module::EAfterConditionPath::c_Continue)
  ;
 
  /* Do not change the names of >, <, ... we use them to serialize conditional paths */
  enum_<Belle2::ModuleCondition::EConditionOperators>("ConditionOperator")
  .value(">", Belle2::ModuleCondition::EConditionOperators::c_GT)
  .value("<", Belle2::ModuleCondition::EConditionOperators::c_ST)
  .value(">=", Belle2::ModuleCondition::EConditionOperators::c_GE)
  .value("<=", Belle2::ModuleCondition::EConditionOperators::c_SE)
  .value("==", Belle2::ModuleCondition::EConditionOperators::c_EQ)
  .value("!=", Belle2::ModuleCondition::EConditionOperators::c_NE)
  ;
 
  enum_<Module::EModulePropFlags>("ModulePropFlags",
                                  R"(Flags to indicate certain low-level features of modules, see :func:`Module.set_property_flags()`, :func:`Module.has_properties()`. Most useful flags are:
 
.. attribute:: PARALLELPROCESSINGCERTIFIED
 
    This module can be run in parallel processing mode safely (All I/O must be done through the data store, in particular, the module must not write any files.)
 
.. attribute:: HISTOGRAMMANAGER
 
  This module is used to manage histograms accumulated by other modules
 
.. attribute:: TERMINATEINALLPROCESSES
 
  When using parallel processing, call this module's terminate() function in all processes. This will also ensure that there is exactly one process (single-core if no parallel modules found) or at least one input, one main and one output process.
)")
  .value("INPUT", Module::EModulePropFlags::c_Input)
  .value("OUTPUT", Module::EModulePropFlags::c_Output)
  .value("PARALLELPROCESSINGCERTIFIED", Module::EModulePropFlags::c_ParallelProcessingCertified)
  .value("HISTOGRAMMANAGER", Module::EModulePropFlags::c_HistogramManager)
  .value("INTERNALSERIALIZER", Module::EModulePropFlags::c_InternalSerializer)
  .value("TERMINATEINALLPROCESSES", Module::EModulePropFlags::c_TerminateInAllProcesses)
  ;
 
  //Python class definition
  class_<Module, PyModule> module("Module", R"(
Base class for Modules.
 
A module is the smallest building block of the framework.
A typical event processing chain consists of a Path containing
modules. By inheriting from this base class, various types of
modules can be created. To use a module, please refer to
:func:`Path.add_module()`.  A list of modules is available by running
``basf2 -m`` or ``basf2 -m package``, detailed information on parameters is
given by e.g. ``basf2 -m RootInput``.
 
The 'Module Development' section in the manual provides detailed information
on how to create modules, setting parameters, or using return values/conditions:
https://xwiki.desy.de/xwiki/rest/p/f4fa4/#HModuleDevelopment
 
)");
  module
  .def("__str__", &Module::getPathString)
  .def("name", &Module::getName, return_value_policy<copy_const_reference>(),
       "Returns the name of the module. Can be changed via :func:`set_name() <Module.set_name()>`, use :func:`type() <Module.type()>` for identifying a particular module class.")
  .def("type", &Module::getType, return_value_policy<copy_const_reference>(),
       "Returns the type of the module (i.e. class name minus 'Module')")
  .def("set_name", &Module::setName, args("name"), R"(
Set custom name, e.g. to distinguish multiple modules of the same type.
 
>>> path.add_module('EventInfoSetter')
>>> ro = path.add_module('RootOutput', branchNames=['EventMetaData'])
>>> ro.set_name('RootOutput_metadata_only')
>>> print(path)
[EventInfoSetter -> RootOutput_metadata_only]
 
)")
  .def("description", &Module::getDescription, return_value_policy<copy_const_reference>(),
       "Returns the description of this module.")
  .def("package", &Module::getPackage, return_value_policy<copy_const_reference>(),
       "Returns the package this module belongs to.")
  .def("available_params", &_getParamInfoListPython,
       "Return list of all module parameters as `ModuleParamInfo` instances")
  .def("has_properties", &Module::hasProperties, (bp::arg("properties")),
       R"DOCSTRING(Allows to check if the module has the given properties out of `ModulePropFlags` set.
 
>>> if module.has_properties(ModulePropFlags.PARALLELPROCESSINGCERTIFIED):
>>>     ...
 
Parameters:
  properties (int): bitmask of `ModulePropFlags` to check for.
)DOCSTRING")
  .def("set_property_flags", &Module::setPropertyFlags, args("property_mask"),
       "Set module properties in the form of an OR combination of `ModulePropFlags`.");
  {
    // python signature is too crowded, make ourselves
    docstring_options subOptions(true, false, false); //userdef, py sigs, c++ sigs
    module
    .def("if_value", &Module::if_value,
         (bp::arg("expression"), bp::arg("condition_path"), bp::arg("after_condition_path")= Module::EAfterConditionPath::c_End),
         R"DOCSTRING(if_value(expression, condition_path, after_condition_path=AfterConditionPath.END)
 
Sets a conditional sub path which will be executed after this
module if the return value set in the module passes the given ``expression``.
 
Modules can define a return value (int or bool) using ``setReturnValue()``,
which can be used in the steering file to split the Path based on this value, for example
 
>>> module_with_condition.if_value("<1", another_path)
 
In case the return value of the ``module_with_condition`` for a given event is
less than 1, the execution will be diverted into ``another_path`` for this event.
 
You could for example set a special return value if an error occurs, and divert
the execution into a path containing :b2:mod:`RootOutput` if it is found;
saving only the data producing/produced by the error.
 
After a conditional path has executed, basf2 will by default stop processing
the path for this event. This behaviour can be changed by setting the
``after_condition_path`` argument.
 
Parameters:
  expression (str): Expression to determine if the conditional path should be executed.
      This should be one of the comparison operators ``<``, ``>``, ``<=``,
      ``>=``, ``==``, or ``!=`` followed by a numerical value for the return value
  condition_path (Path): path to execute in case the expression is fulfilled
  after_condition_path (AfterConditionPath): What to do once the ``condition_path`` has been executed.
)DOCSTRING")
    .def("if_false", &Module::if_false,
         (bp::arg("condition_path"), bp::arg("after_condition_path")= Module::EAfterConditionPath::c_End),
         R"DOC(if_false(condition_path, after_condition_path=AfterConditionPath.END)
 
Sets a conditional sub path which will be executed after this module if
the return value of the module evaluates to False. This is equivalent to
calling `if_value` with ``expression=\"<1\"``)DOC")
    .def("if_true", &Module::if_true,
         (bp::arg("condition_path"), bp::arg("after_condition_path")= Module::EAfterConditionPath::c_End),
         R"DOC(if_true(condition_path, after_condition_path=AfterConditionPath.END)
 
Sets a conditional sub path which will be executed after this module if
the return value of the module evaluates to True. It is equivalent to
calling `if_value` with ``expression=\">=1\"``)DOC");
  }
  module
  .def("has_condition", &Module::hasCondition,
       "Return true if a conditional path has been set for this module "
       "using `if_value`, `if_true` or `if_false`")
  .def("get_all_condition_paths", &_getAllConditionPathsPython,
       "Return a list of all conditional paths set for this module using "
       "`if_value`, `if_true` or `if_false`")
  .def("get_all_conditions", &_getAllConditionsPython,
       "Return a list of all conditional path expressions set for this module using "
       "`if_value`, `if_true` or `if_false`")
  .add_property("logging", make_function(&Module::getLogConfig, return_value_policy<reference_existing_object>()),
                &Module::setLogConfig)

ExtractNEvent()

void ExtractNEvent ( std::vector< TH1 * > & hs )

staticinherited

Extract event processed from daq histogram, called from input module.

Definition at line 361 of file DQMHistAnalysis.cc.

{
  s_eventProcessed = 0;
  for (size_t i = 0; i < hs.size(); i++) {
    if (hs[i]->GetName() == std::string("DAQ/Nevent")) {
      s_eventProcessed = hs[i]->GetEntries();
      return;
    }
  }
  B2ERROR("ExtractEvent: Histogram \"DAQ/Nevent\" missing");
}

ExtractRunType()

void ExtractRunType ( std::vector< TH1 * > & hs )

staticinherited

Extract Run Type from histogram title, called from input module.

Definition at line 349 of file DQMHistAnalysis.cc.

{
  s_runType = "";
  for (size_t i = 0; i < hs.size(); i++) {
    if (hs[i]->GetName() == std::string("DQMInfo/rtype")) {
      s_runType = hs[i]->GetTitle();
      return;
    }
  }
  B2ERROR("ExtractRunType: Histogram \"DQMInfo/rtype\" missing");
}

fillMaskedChannelsHistogram()

void fillMaskedChannelsHistogram ( TH1 * histogram, TCanvas * canvas )

private

Fill histogram containing masked channels per sector.

Parameters

[in]	histogram	Histogram.
[in]	canvas	Canvas.

Definition at line 472 of file DQMHistAnalysisKLM.cc.

{
  histogram->Clear();
  canvas->Clear();
  canvas->cd();
  if (m_MaskedChannels.size() > 0) {
    int channelSubdetector, channelSection, channelSector;
    int layer, plane, strip;
    for (KLMChannelNumber channel : m_MaskedChannels) {
      m_ElementNumbers->channelNumberToElementNumbers(
        channel, &channelSubdetector, &channelSection, &channelSector,
        &layer, &plane, &strip);
      KLMSectorNumber sectorNumber;
      if (channelSubdetector == KLMElementNumbers::c_BKLM)
        sectorNumber = m_ElementNumbers->sectorNumberBKLM(channelSection, channelSector);
      else
        sectorNumber = m_ElementNumbers->sectorNumberEKLM(channelSection, channelSector);
      KLMSectorNumber sectorIndex = m_SectorArrayIndex->getIndex(sectorNumber);
      histogram->Fill(sectorIndex);
    }
  }
  histogram->SetStats(false);
  histogram->SetTitle("Number of masked channels per sector");
  histogram->Draw();
  canvas->Modified();
  canvas->Update();
}

findCanvas()

TCanvas * findCanvas ( TString cname )

staticinherited

Find canvas by name.

Parameters

cname

Name of the canvas

Returns

The pointer to the canvas, or nullptr if not found.

Definition at line 138 of file DQMHistAnalysis.cc.

{
  TIter nextkey(gROOT->GetListOfCanvases());
  TObject* obj{};
 
  while ((obj = dynamic_cast<TObject*>(nextkey()))) {
    if (obj->IsA()->InheritsFrom("TCanvas")) {
      if (obj->GetName() == canvas_name)
        return dynamic_cast<TCanvas*>(obj);
    }
  }
  return nullptr;
}

findHist()

TH1 * findHist ( const std::string & dirname, const std::string & histname = "", bool onlyIfUpdated = false )

staticinherited

Find histogram.

Parameters

dirname	The name of the directory.
histname	The name of the histogram.
onlyIfUpdated	req only updated hists, return nullptr otherwise

Returns

The found histogram, or nullptr if not found.

Definition at line 153 of file DQMHistAnalysis.cc.

{
  std::string fullname = dirname + "/" + histname;
  if (dirname.size() == 0) fullname = histname; // assume contains dirname
  if (histname.size() == 0) fullname = dirname; // assume contains histname
 
  if (s_histList.find(fullname) != s_histList.end()) {
    if (was_updated && !s_histList[fullname].isUpdated()) return nullptr;
    if (s_histList[fullname].getHist()) {
      return s_histList[fullname].getHist();
    } else {
      B2ERROR("Histogram " << fullname << " in histogram list but nullptr.");
    }
  }
  B2INFO("Histogram " << fullname << " not in list.");
  return nullptr;
}

findHistInCanvas()

TH1 * findHistInCanvas ( const std::string & hname, TCanvas ** canvas = nullptr )

inherited

Find histogram in corresponding canvas.

Parameters

hname	Name of the histogram (dir+name)
canvas	ptr to specific canvas ptr or nullptr

Returns

The pointer to the histogram, or nullptr if not found.

Definition at line 213 of file DQMHistAnalysis.cc.

{
  TCanvas* cnv = nullptr;
  // try to get canvas from outside
  if (cobj) cnv = *cobj;
  // if no canvas search for it
  if (cnv == nullptr) {
    // parse the dir+histo name and create the corresponding canvas name
    auto s = StringSplit(histo_name, '/');
    if (s.size() != 2) {
      B2ERROR("findHistInCanvas: histoname not valid (missing dir?), should be 'dirname/histname': " << histo_name);
      return nullptr;
    }
    auto dirname = s.at(0);
    auto hname = s.at(1);
    std::string canvas_name = dirname + "/c_" + hname;
    cnv = findCanvas(canvas_name);
    // set canvas pointer for outside
    if (cnv && cobj) *cobj = cnv;
  }
 
  // get histogram pointer
  if (cnv != nullptr) {
    TIter nextkey(cnv->GetListOfPrimitives());
    TObject* obj{};
    while ((obj = dynamic_cast<TObject*>(nextkey()))) {
      if (obj->IsA()->InheritsFrom("TH1")) {
        if (obj->GetName() == histo_name)
          return  dynamic_cast<TH1*>(obj);
      }
    }
  }
  return nullptr;
}

findHistInFile()

TH1 * findHistInFile ( TFile * file, const std::string & histname )

staticinherited

Find histogram in specific TFile (e.g.

ref file).

Parameters

file	The TFile to search.
histname	The name of the histogram, can incl directory

Returns

The found histogram, or nullptr if not found.

Definition at line 248 of file DQMHistAnalysis.cc.

{
  // find histogram by name in file, histname CAN contain directory!
  // will return nullptr if file is zeroptr, not found or not correct type
  if (file && file->IsOpen()) {
    auto obj = file->Get(histname.data());
    if (obj != nullptr) {
      // check class type
      if (obj->IsA()->InheritsFrom("TH1")) {
        B2DEBUG(20, "Histogram " << histname << " found in file");
        return dynamic_cast<TH1*>(obj);
      } else {
        B2INFO("Found Object " << histname << " in file is not a histogram");
      }
    } else {
      B2INFO("Histogram " << histname << " not found in file");
    }
  }
  return nullptr;
}

findMonitoringObject()

MonitoringObject * findMonitoringObject ( const std::string & objName )

staticinherited

Find MonitoringObject.

Parameters

objName

The name of the MonitoringObject.

Returns

The found MonitoringObject, or nullptr if not found.

Definition at line 269 of file DQMHistAnalysis.cc.

{
  if (s_monObjList.find(objName) != s_monObjList.end()) {
    return &s_monObjList[objName];
  }
  B2INFO("MonitoringObject " << objName << " not in memfile.");
  return nullptr;
}

findRefHist()

TH1 * findRefHist ( const std::string & dirname, const std::string & histname = "", ERefScaling scaling = ERefScaling::c_RefScaleNone, const TH1 * hist = nullptr )

staticinherited

Find reference histogram.

Parameters

dirname	The name of the directory.
histname	The name of the histogram.
scaling	enum what scaling to use
hist	histogram to scale to

Returns

The found histogram, or nullptr if not found.

Definition at line 197 of file DQMHistAnalysis.cc.

{
  std::string fullname = dirname + "/" + histname;
  if (dirname.size() == 0) fullname = histname; // assume contains dirname
  if (histname.size() == 0) fullname = dirname; // assume contains histname
 
  if (s_refList.find(fullname) != s_refList.end()) {
    // get a copy of the reference which we can modify
    // (it is still owned and managed by the framework)
    // then do the scaling
    return scaleReference(scaling, hist, s_refList[fullname].getReference());
  }
  return nullptr;
}

getAfterConditionPath()

Module::EAfterConditionPath getAfterConditionPath ( ) const

inherited

What to do after the conditional path is finished.

(defaults to c_End if no condition is set)

Definition at line 133 of file Module.cc.

{
  if (m_conditions.empty()) return EAfterConditionPath::c_End;
 
  //okay, a condition was set for this Module...
  if (!m_hasReturnValue) {
    B2FATAL("A condition was set for '" << getName() << "', but the module did not set a return value!");
  }
 
  for (const auto& condition : m_conditions) {
    if (condition.evaluate(m_returnValue)) {
      return condition.getAfterConditionPath();
    }
  }
 
  return EAfterConditionPath::c_End;
}

getAllConditionPaths()

std::vector< std::shared_ptr< Path > > getAllConditionPaths ( ) const

inherited

Return all condition paths currently set (no matter if the condition is true or not).

Definition at line 150 of file Module.cc.

{
  std::vector<std::shared_ptr<Path>> allConditionPaths;
  for (const auto& condition : m_conditions) {
    allConditionPaths.push_back(condition.getPath());
  }
 
  return allConditionPaths;
}

getAllConditions()

const std::vector< ModuleCondition > & getAllConditions ( ) const

inlineinherited

Return all set conditions for this module.

Definition at line 323 of file Module.h.

    {
      return m_conditions;
    }

getCanvasUpdatedList()

static const CanvasUpdatedList & getCanvasUpdatedList ( )

inlinestaticinherited

Get the list of the canvas update status.

Returns

The list of the canvases.

Definition at line 206 of file DQMHistAnalysis.h.

{ return s_canvasUpdatedList;};

getCondition()

const ModuleCondition * getCondition ( ) const

inlineinherited

Return a pointer to the first condition (or nullptr, if none was set)

Definition at line 313 of file Module.h.

    {
      if (m_conditions.empty()) {
        return nullptr;
      } else {
        return &m_conditions.front();
      }
    }

getConditionPath()

std::shared_ptr< Path > getConditionPath ( ) const

inherited

Returns the path of the last true condition (if there is at least one, else reaturn a null pointer).

Definition at line 113 of file Module.cc.

{
  PathPtr p;
  if (m_conditions.empty()) return p;
 
  //okay, a condition was set for this Module...
  if (!m_hasReturnValue) {
    B2FATAL("A condition was set for '" << getName() << "', but the module did not set a return value!");
  }
 
  for (const auto& condition : m_conditions) {
    if (condition.evaluate(m_returnValue)) {
      return condition.getPath();
    }
  }
 
  // if none of the conditions were true, return a null pointer.
  return p;
}

getDelta()

TH1 * getDelta ( const std::string & dirname, const std::string & histname = "", int n = 0, bool onlyIfUpdated = true )

inherited

Get Delta histogram.

Parameters

dirname	directory
histname	name of histogram
n	index of delta histogram, 0 is most recent one
onlyIfUpdated	req only updated deltas, return nullptr otherwise

Returns

delta histogram or nullptr

Definition at line 117 of file DQMHistAnalysis.cc.

{
  std::string fullname = dirname + "/" + histname;
  if (dirname.size() == 0) fullname = histname; // assume contains dirname
  if (histname.size() == 0) fullname = dirname; // assume contains histname
 
  auto it = s_deltaList.find(fullname);
  if (it != s_deltaList.end()) {
    return it->second.getDelta(n, onlyIfUpdated);
  }
  B2WARNING("Delta hist " << fullname << " not found");
  return nullptr;
}

getDeltaList()

static const DeltaList & getDeltaList ( )

inlinestaticinherited

Get the list of the delta histograms.

Returns

The list of the delta histograms.

Definition at line 200 of file DQMHistAnalysis.h.

{ return s_deltaList;};

getDescription()

const std::string & getDescription ( ) const

inlineinherited

Returns the description of the module.

Definition at line 201 of file Module.h.

{return m_description;}

getEpicsPV() [1/2]

double getEpicsPV ( const std::string & keyname )

inherited

Read value from a EPICS PV.

Parameters

keyname

key name (or full PV name) of PV

Returns

value or NAN if not existing

Definition at line 492 of file DQMHistAnalysis.cc.

{
  double value{NAN};
  if (!m_useEpics) return value;
#ifdef _BELLE2_EPICS
  if (m_epicsNameToChID[keyname] == nullptr) {
    B2ERROR("Epics PV " << keyname << " not registered!");
    return value;
  }
  // From EPICS doc. When ca_get or ca_array_get are invoked the returned channel value can't be assumed to be stable
  // in the application supplied buffer until after ECA_NORMAL is returned from ca_pend_io. If a connection is lost
  // outstanding get requests are not automatically reissued following reconnect.
  auto r = ca_get(DBR_DOUBLE, m_epicsNameToChID[keyname], &value);
  if (r == ECA_NORMAL) r = ca_pend_io(5.0); // this is needed!
  if (r == ECA_NORMAL) {
    return value;
  } else {
    CheckEpicsError(r, "Read PV failed in ca_get or ca_pend_io failure", keyname);
  }
#endif
  return NAN;
}

getEpicsPV() [2/2]

double getEpicsPV ( int index )

inherited

Read value from a EPICS PV.

Parameters

index

index of PV

Returns

value or NAN if not existing

Definition at line 515 of file DQMHistAnalysis.cc.

{
  double value{NAN};
  if (!m_useEpics) return value;
#ifdef _BELLE2_EPICS
  if (index < 0 || index >= (int)m_epicsChID.size()) {
    B2ERROR("Epics PV with " << index << " not registered!");
    return value;
  }
  // From EPICS doc. When ca_get or ca_array_get are invoked the returned channel value can't be assumed to be stable
  // in the application supplied buffer until after ECA_NORMAL is returned from ca_pend_io. If a connection is lost
  // outstanding get requests are not automatically reissued following reconnect.
  auto r = ca_get(DBR_DOUBLE, m_epicsChID[index], &value);
  if (r == ECA_NORMAL) r = ca_pend_io(5.0); // this is needed!
  if (r == ECA_NORMAL) {
    return value;
  } else {
    CheckEpicsError(r, "Read PV failed in ca_get or ca_pend_io failure", m_epicsChID[index]);
  }
#endif
  return NAN;
}

getEpicsPVChID() [1/2]

chid getEpicsPVChID ( const std::string & keyname )

inherited

Get EPICS PV Channel Id.

Parameters

keyname

key name (or full PV name) of PV

Returns

Channel ID is written on success, otherwise nullptr

Definition at line 588 of file DQMHistAnalysis.cc.

{
#ifdef _BELLE2_EPICS
  if (m_useEpics) {
    if (m_epicsNameToChID[keyname] != nullptr) {
      return m_epicsNameToChID[keyname];
    } else {
      B2ERROR("Epics PV " << keyname << " not registered!");
    }
  }
#endif
  return nullptr;
}

getEpicsPVChID() [2/2]

chid getEpicsPVChID ( int index )

inherited

Get EPICS PV Channel Id.

Parameters

index

index of PV

Returns

Channel ID is written on success, otherwise nullptr

Definition at line 602 of file DQMHistAnalysis.cc.

{
#ifdef _BELLE2_EPICS
  if (m_useEpics) {
    if (index >= 0 && index < (int)m_epicsChID.size()) {
      return m_epicsChID[index];
    } else {
      B2ERROR("Epics PV with " << index << " not registered!");
    }
  }
#endif
  return nullptr;
}

getEpicsStringPV() [1/2]

std::string getEpicsStringPV ( const std::string & keyname, bool & status )

inherited

Read value from a EPICS PV.

Parameters

keyname	key name (or full PV name) of PV
status	return status (true on success)

Returns

string value (empty string if non existing)

Definition at line 538 of file DQMHistAnalysis.cc.

{
  status = false;
  char value[40] = "";
  if (!m_useEpics) return std::string(value);
#ifdef _BELLE2_EPICS
  if (m_epicsNameToChID[keyname] == nullptr) {
    B2ERROR("Epics PV " << keyname << " not registered!");
    return std::string(value);
  }
  // From EPICS doc. When ca_get or ca_array_get are invoked the returned channel value can't be assumed to be stable
  // in the application supplied buffer until after ECA_NORMAL is returned from ca_pend_io. If a connection is lost
  // outstanding get requests are not automatically reissued following reconnect.
  auto r = ca_get(DBR_STRING, m_epicsNameToChID[keyname], value);
  if (r == ECA_NORMAL) r = ca_pend_io(5.0); // this is needed!
  if (r == ECA_NORMAL) {
    status = true;
    return std::string(value);
  } else {
    CheckEpicsError(r, "Read PV (string) failed in ca_get or ca_pend_io failure", keyname);
  }
#endif
  return std::string(value);
}

getEpicsStringPV() [2/2]

std::string getEpicsStringPV ( int index, bool & status )

inherited

Read value from a EPICS PV.

Parameters

index	index of PV
status	return status (true on success)

Returns

string value (empty string if non existing)

Definition at line 563 of file DQMHistAnalysis.cc.

{
  status = false;
  char value[40] = "";
  if (!m_useEpics) return std::string(value);
#ifdef _BELLE2_EPICS
  if (index < 0 || index >= (int)m_epicsChID.size()) {
    B2ERROR("Epics PV with " << index << " not registered!");
    return std::string(value);
  }
  // From EPICS doc. When ca_get or ca_array_get are invoked the returned channel value can't be assumed to be stable
  // in the application supplied buffer until after ECA_NORMAL is returned from ca_pend_io. If a connection is lost
  // outstanding get requests are not automatically reissued following reconnect.
  auto r = ca_get(DBR_STRING, m_epicsChID[index], value);
  if (r == ECA_NORMAL) r = ca_pend_io(5.0); // this is needed!
  if (r == ECA_NORMAL) {
    status = true;
    return std::string(value);
  } else {
    CheckEpicsError(r, "Read PV (string) failed in ca_get or ca_pend_io failure", m_epicsChID[index]);
  }
#endif
  return std::string(value);
}

getEventProcessed()

static int getEventProcessed ( void )

inlinestaticinherited

Get the number of processed events.

(Attention, asynch histogram updates!)

Returns

Processed events.

Definition at line 225 of file DQMHistAnalysis.h.

{ return s_eventProcessed;};

getFileNames()

virtual std::vector< std::string > getFileNames ( bool outputFiles )

inlinevirtualinherited

Return a list of output filenames for this modules.

This will be called when basf2 is run with "--dry-run" if the module has set either the c_Input or c_Output properties.

If the parameter outputFiles is false (for modules with c_Input) the list of input filenames should be returned (if any). If outputFiles is true (for modules with c_Output) the list of output files should be returned (if any).

If a module has sat both properties this member is called twice, once for each property.

The module should return the actual list of requested input or produced output filenames (including handling of input/output overrides) so that the grid system can handle input/output files correctly.

This function should return the same value when called multiple times. This is especially important when taking the input/output overrides from Environment as they get consumed when obtained so the finalized list of output files should be stored for subsequent calls.

Reimplemented in RootInputModule, RootOutputModule, and StorageRootOutputModule.

Definition at line 133 of file Module.h.

    {
      return std::vector<std::string>();
    }

getHistList()

static HistList & getHistList ( )

inlinestaticinherited

Get the list of the histograms.

Returns

The list of the histograms.

Definition at line 188 of file DQMHistAnalysis.h.

{ return s_histList;};

getLogConfig()

LogConfig & getLogConfig ( )

inlineinherited

Returns the log system configuration.

Definition at line 224 of file Module.h.

{return m_logConfig;}

getModules()

std::list< ModulePtr > getModules ( ) const

inlineoverrideprivatevirtualinherited

no submodules, return empty list

Implements PathElement.

Definition at line 505 of file Module.h.

{ return std::list<ModulePtr>(); }

getMonitoringObject()

MonitoringObject * getMonitoringObject ( const std::string & name )

staticinherited

Get MonitoringObject with given name (new object is created if non-existing)

Parameters

name	name of MonitoringObject to get

Returns

The MonitoringObject

Definition at line 131 of file DQMHistAnalysis.cc.

{
  auto obj = &s_monObjList[objName];
  obj->SetName(objName.c_str());
  return obj;
}

getMonObjList()

static MonObjList & getMonObjList ( )

inlinestaticinherited

Get the list of MonitoringObjects.

Returns

The list of the MonitoringObjects.

Definition at line 194 of file DQMHistAnalysis.h.

{ return s_monObjList;};

getName()

const std::string & getName ( ) const

inlineinherited

Returns the name of the module.

This can be changed via e.g. set_name() in the steering file to give more useful names if there is more than one module of the same type.

For identifying the type of a module, using getType() (or type() in Python) is recommended.

Definition at line 186 of file Module.h.

{return m_name;}

getPackage()

const std::string & getPackage ( ) const

inlineinherited

Returns the package this module is in.

Definition at line 196 of file Module.h.

{return m_package;}

getParamInfoListPython()

std::shared_ptr< boost::python::list > getParamInfoListPython ( ) const

inherited

Returns a python list of all parameters.

Each item in the list consists of the name of the parameter, a string describing its type, a python list of all default values and the description of the parameter.

Returns

A python list containing the parameters of this parameter list.

Definition at line 279 of file Module.cc.

{
  return m_moduleParamList.getParamInfoListPython();
}

getParamList()

const ModuleParamList & getParamList ( ) const

inlineinherited

Return module param list.

Definition at line 362 of file Module.h.

{ return m_moduleParamList; }

getPathString()

std::string getPathString ( ) const

overrideprivatevirtualinherited

return the module name.

Implements PathElement.

Definition at line 192 of file Module.cc.

{
 
  std::string output = getName();
 
  for (const auto& condition : m_conditions) {
    output += condition.getString();
  }
 
  return output;
}

getProcessedEvents()

double getProcessedEvents ( )

private

Get number of processed events.

Definition at line 195 of file DQMHistAnalysisKLM.cc.

{
  if (DQMHistAnalysisModule::getEventProcessed() == 0.) {
    B2WARNING("Either DAQ/Nevent is not found or Nevent = 0.");
    /* Set the minimal number of processed events to 0 if we can't determine the processed events. */
    m_MinProcessedEventsForMessages = 0.;
  }
  return (double) DQMHistAnalysisModule::getEventProcessed();
}

getPVPrefix()

static std::string & getPVPrefix ( void )

inlinestaticinherited

get global prefix for EPICS PVs

Returns

prefix in use

Definition at line 588 of file DQMHistAnalysis.h.

{return m_PVPrefix;};

getReturnValue()

int getReturnValue ( ) const

inlineinherited

Return the return value set by this module.

This value is only meaningful if hasReturnValue() is true

Definition at line 380 of file Module.h.

{ return m_returnValue; }

getRunType()

static const std::string & getRunType ( void )

inlinestaticinherited

Get the list of the reference histograms.

Returns

The list of the reference histograms. Get the Run Type.

Run type string.

Definition at line 219 of file DQMHistAnalysis.h.

{ return s_runType;};

getSigma68()

double getSigma68 ( TH1 * h )

staticinherited

Helper function to compute half of the central interval covering 68% of a distribution.

This quantity is an alternative to the standard deviation.

Parameters

h

histogram

Returns

Half of the central interval covering 68% of a distribution.

Definition at line 278 of file DQMHistAnalysis.cc.

{
  double probs[2] = {0.16, 1 - 0.16};
  double quant[2] = {0, 0};
  h->GetQuantiles(2, quant, probs);
  const double sigma68 = (-quant[0] + quant[1]) / 2;
  return sigma68;
}

getStatusColor()

DQMHistAnalysisModule::EStatusColor getStatusColor ( EStatus status )

staticinherited

Return color for canvas state.

Parameters

status

canvas status

Returns

alarm color

Definition at line 708 of file DQMHistAnalysis.cc.

{
  // white color is the default, if no colorize
  EStatusColor color = c_ColorDefault;
  switch (stat) {
    case c_StatusTooFew:
      color = c_ColorTooFew; // Magenta or Gray
      break;
    case c_StatusDefault:
      color = c_ColorDefault; // default no colors
      break;
    case c_StatusGood:
      color = c_ColorGood; // Good
      break;
    case c_StatusWarning:
      color = c_ColorWarning; // Warning
      break;
    case c_StatusError:
      color = c_ColorError; // Severe
      break;
    default:
      color = c_ColorDefault; // default no colors
      break;
  }
  return color;
}

getType()

const std::string & getType ( ) const

inherited

Returns the type of the module (i.e.

class name minus 'Module')

Definition at line 41 of file Module.cc.

{
  if (m_type.empty())
    B2FATAL("Module type not set for " << getName());
  return m_type;
}

getUseEpics()

static bool getUseEpics ( void )

inlinestaticinherited

Getter for EPICS usage.

Returns

flag is in use

Definition at line 571 of file DQMHistAnalysis.h.

{return m_useEpics;};

getUseEpicsReadOnly()

static bool getUseEpicsReadOnly ( void )

inlinestaticinherited

Getter EPICS flag in read only mode.

Returns

flag if read only

Definition at line 577 of file DQMHistAnalysis.h.

{return m_epicsReadOnly;};

hasCondition()

bool hasCondition ( ) const

inlineinherited

Returns true if at least one condition was set for the module.

Definition at line 310 of file Module.h.

{ return not m_conditions.empty(); };

hasDeltaPar()

bool hasDeltaPar ( const std::string & dirname, const std::string & histname )

staticinherited

Check if Delta histogram parameters exist for histogram.

Parameters

dirname	directory
histname	name of histogram

Returns

true if parameters have been set already

Definition at line 106 of file DQMHistAnalysis.cc.

{
  std::string fullname;
  if (dirname.size() > 0) {
    fullname = dirname + "/" + histname;
  } else {
    fullname = histname;
  }
  return s_deltaList.find(fullname) != s_deltaList.end(); // contains() if we switch to C++20
}

hasProperties()

bool hasProperties ( unsigned int propertyFlags ) const

inherited

Returns true if all specified property flags are available in this module.

Parameters

propertyFlags

Ored EModulePropFlags which should be compared with the module flags.

Definition at line 160 of file Module.cc.

{
  return (propertyFlags & m_propertyFlags) == propertyFlags;
}

hasReturnValue()

bool hasReturnValue ( ) const

inlineinherited

Return true if this module has a valid return value set.

Definition at line 377 of file Module.h.

{ return m_hasReturnValue; }

hasUnsetForcedParams()

bool hasUnsetForcedParams ( ) const

inherited

Returns true and prints error message if the module has unset parameters which the user has to set in the steering file.

Definition at line 166 of file Module.cc.

{
  auto missing = m_moduleParamList.getUnsetForcedParams();
  std::string allMissing = "";
  for (const auto& s : missing)
    allMissing += s + " ";
  if (!missing.empty())
    B2ERROR("The following required parameters of Module '" << getName() << "' were not specified: " << allMissing <<
            "\nPlease add them to your steering file.");
  return !missing.empty();
}

if_false()

void if_false ( const std::shared_ptr< Path > & path, EAfterConditionPath afterConditionPath = EAfterConditionPath::c_End )

inherited

A simplified version to add a condition to the module.

Please note that successive calls of this function will add more than one condition to the module. If more than one condition results in true, only the last of them will be used.

Please be careful: Avoid creating cyclic paths, e.g. by linking a condition to a path which is processed before the path where this module is located in.

It is equivalent to the if_value() method, using the expression "<1". This method is meant to be used together with the setReturnValue(bool value) method.

Parameters

path	Shared pointer to the Path which will be executed if the return value is false.
afterConditionPath	What to do after executing 'path'.

Definition at line 85 of file Module.cc.

{
  if_value("<1", path, afterConditionPath);
}

if_true()

void if_true ( const std::shared_ptr< Path > & path, EAfterConditionPath afterConditionPath = EAfterConditionPath::c_End )

inherited

A simplified version to set the condition of the module.

Please note that successive calls of this function will add more than one condition to the module. If more than one condition results in true, only the last of them will be used.

Please be careful: Avoid creating cyclic paths, e.g. by linking a condition to a path which is processed before the path where this module is located in.

It is equivalent to the if_value() method, using the expression ">=1". This method is meant to be used together with the setReturnValue(bool value) method.

Parameters

path	Shared pointer to the Path which will be executed if the return value is true.
afterConditionPath	What to do after executing 'path'.

Definition at line 90 of file Module.cc.

{
  if_value(">=1", path, afterConditionPath);
}

if_value()

void if_value ( const std::string & expression, const std::shared_ptr< Path > & path, EAfterConditionPath afterConditionPath = EAfterConditionPath::c_End )

inherited

Add a condition to the module.

Please note that successive calls of this function will add more than one condition to the module. If more than one condition results in true, only the last of them will be used.

See https://xwiki.desy.de/xwiki/rest/p/a94f2 or ModuleCondition for a description of the syntax.

Please be careful: Avoid creating cyclic paths, e.g. by linking a condition to a path which is processed before the path where this module is located in.

Parameters

expression	The expression of the condition.
path	Shared pointer to the Path which will be executed if the condition is evaluated to true.
afterConditionPath	What to do after executing 'path'.

Definition at line 79 of file Module.cc.

{
  m_conditions.emplace_back(expression, path, afterConditionPath);
}

initHistListBeforeEvent()

void initHistListBeforeEvent ( void )

staticinherited

Reset the list of histograms.

Definition at line 308 of file DQMHistAnalysis.cc.

{
  for (auto& it : s_histList) {
    // attention, we must use reference, otherwise we work on a copy
    it.second.resetBeforeEvent();
  }
  for (auto& it : s_deltaList) {
    // attention, we must use reference, otherwise we work on a copy
    it.second.setNotUpdated();
  }
 
  s_canvasUpdatedList.clear();
}

initialize()

void initialize ( void )

finaloverridevirtual

Initializer.

Reimplemented from Module.

Definition at line 71 of file DQMHistAnalysisKLM.cc.

{
  m_monObj = getMonitoringObject("klm");
 
  if (m_ThresholdForHot > m_ThresholdForMasked)
    B2FATAL("The threshold used for hot channels is larger than the one for masked channels."
            << LogVar("Threshold for hot channels", m_ThresholdForHot)
            << LogVar("Threshold for masked channels", m_ThresholdForMasked));
 
  // register plots for delta histogramming
  addDeltaPar(m_histogramDirectoryName, "time_rpc", HistDelta::c_Entries, m_minEntries, 1);
  addDeltaPar(m_histogramDirectoryName, "time_scintillator_bklm", HistDelta::c_Entries, m_minEntries, 1);
  addDeltaPar(m_histogramDirectoryName, "time_scintillator_eklm", HistDelta::c_Entries, m_minEntries, 1);
 
  // The FE ratio via pair: in c_Entries one fills far more slowly than other,
  // leaving the pair out of sync and now use c_Events to update over the same event window.
  addDeltaPar(m_histogramDirectoryName, "feStatus_bklm_scintillator_layers_0", HistDelta::c_Events, m_minEvents, 1);
  addDeltaPar(m_histogramDirectoryName, "feStatus_bklm_scintillator_layers_1", HistDelta::c_Events, m_minEvents, 1);
  addDeltaPar(m_histogramDirectoryName, "feStatus_eklm_plane_0", HistDelta::c_Events, m_minEvents, 1);
  addDeltaPar(m_histogramDirectoryName, "feStatus_eklm_plane_1", HistDelta::c_Events, m_minEvents, 1);
 
  //register EPICS PVs
  registerEpicsPV("KLM:MaskedChannels", "MaskedChannels");
  registerEpicsPV("KLM:DeadBarrelModules", "DeadBarrelModules");
  registerEpicsPV("KLM:DeadEndcapModules", "DeadEndcapModules");
 
  gROOT->cd();
  std::string c_masked_channels_name = m_histogramDirectoryName + "/c_masked_channels";
  m_c_masked_channels = new TCanvas((c_masked_channels_name).c_str());
  std::string c_fe_bklm_ratio_name = m_histogramDirectoryName + "/c_fe_bklm_ratio";
  m_c_fe_bklm_ratio = new TCanvas((c_fe_bklm_ratio_name).c_str());
  std::string c_fe_eklm_ratio_name = m_histogramDirectoryName + "/c_fe_eklm_ratio";
  m_c_fe_eklm_ratio = new TCanvas((c_fe_eklm_ratio_name).c_str());
  // Feature extraction status histogram for BKLM
  int bklmSectors = BKLMElementNumbers::getMaximalSectorGlobalNumber(); // 16
  int eklmPlanes = EKLMElementNumbers::getMaximalPlaneGlobalNumber(); // 208
  //BKLM
  std::string fe_bklm_ratio_name = m_histogramDirectoryName + "/fe_bklm_ratio";
  m_fe_bklm_ratio = new TH1F((fe_bklm_ratio_name).c_str(),
                             "FE Ratio in BKLM",
                             bklmSectors * 2, 0.5, 0.5 + bklmSectors * 2);
  m_fe_bklm_ratio->GetXaxis()->SetTitle("Scintillator Layer number");
  m_fe_bklm_ratio->SetStats(false);
  m_fe_bklm_ratio->SetOption("HIST");
  //EKLM
  std::string fe_eklm_ratio_name = m_histogramDirectoryName + "/fe_eklm_ratio";
  m_fe_eklm_ratio = new TH1F((fe_eklm_ratio_name).c_str(),
                             "FE Ratio in EKLM",
                             eklmPlanes, 0.5, 0.5 + eklmPlanes);
  m_fe_eklm_ratio->GetXaxis()->SetTitle("Plane number");
  m_fe_eklm_ratio->SetStats(false);
  m_fe_eklm_ratio->SetOption("HIST");
  /* Masked channels per sector. */
  KLMSectorNumber totalSectors = m_SectorArrayIndex->getNElements();
  m_MaskedChannelsHist = new TH1F("masked_channels", "Number of masked channels per sector",
                                  totalSectors, -0.5, totalSectors - 0.5);
 
  std::string str;
  KLMChannelIndex klmIndex(KLMChannelIndex::c_IndexLevelSector);
  for (const KLMChannelIndex& klmSector : klmIndex) {
    std::string label = m_ElementNumbers->getSectorDAQName(klmSector.getSubdetector(), klmSector.getSection(), klmSector.getSector());
    KLMSectorNumber sector = klmSector.getKLMSectorNumber();
    KLMSectorNumber sectorIndex = m_SectorArrayIndex->getIndex(sector);
    m_MaskedChannelsHist->GetXaxis()->SetBinLabel(sectorIndex + 1, label.c_str());
    int nHistograms;
    if (klmSector.getSubdetector() == KLMElementNumbers::c_BKLM)
      nHistograms = 2;
    else
      nHistograms = 3;
    for (int j = 0; j < nHistograms; j++) {
      str = "strip_hits_subdetector_" +
            std::to_string(klmSector.getSubdetector()) +
            "_section_" + std::to_string(klmSector.getSection()) +
            "_sector_" + std::to_string(klmSector.getSector()) +
            "_" + std::to_string(j);
      addDeltaPar(m_histogramDirectoryName, str, HistDelta::c_Entries, m_minEntries, 1);
 
    }
  }
 
}

makeStatus()

DQMHistAnalysisModule::EStatus makeStatus ( bool enough, bool warn_flag, bool error_flag )

staticinherited

Helper function to judge the status for coloring and EPICS.

Parameters

enough	enough events for judging
warn_flag	outside of expected range
error_flag	outside of warning range

Returns

the status

Definition at line 690 of file DQMHistAnalysis.cc.

{
  // white color is the default, if no colorize
  if (!enough) {
    return (c_StatusTooFew);
  } else {
    if (error_flag) {
      return (c_StatusError);
    } else if (warn_flag) {
      return (c_StatusWarning);
    } else {
      return (c_StatusGood);
    }
  }
 
  return (c_StatusDefault); // default, but should not be reached
}

printPVStatus()

void printPVStatus ( chid pv, bool onlyError = true )

staticinherited

check the status of a PVs and report if disconnected or not found

Parameters

pv	the chid of the PV to check
onlyError	print only if in error condition (default)

Definition at line 758 of file DQMHistAnalysis.cc.

{
  if (pv == nullptr) {
    B2WARNING("PV chid was nullptr");
    return;
  }
  auto state = ca_state(pv);
  switch (state) {
    case cs_never_conn: /* valid chid, server not found or unavailable */
      B2WARNING("Channel never connected " << ca_name(pv));
      break;
    case cs_prev_conn:  /* valid chid, previously connected to server */
      B2WARNING("Channel was connected, but now is not " << ca_name(pv));
      break;
    case cs_closed:   /* channel deleted by user */
      B2WARNING("Channel deleted already " << ca_name(pv));
      break;
    case cs_conn:       /* valid chid, connected to server */
      if (!onlyError) B2INFO("Channel connected and OK " << ca_name(pv));
      break;
    default:
      B2WARNING("Undefined status for channel " << ca_name(pv));
      break;
  }
}

processFEHistogram()

void processFEHistogram ( TH1 * feHist, const std::string & histName, TCanvas * canvas )

private

Process histogram containing the efficiencies.

Parameters

[in]	feHist	Histogram itself.
[in]	histName	Name of the histogram.
[in]	canvas	Canvas of interest.

Definition at line 636 of file DQMHistAnalysisKLM.cc.

{
  if (!feHist) {
    B2WARNING("processFEHistogram: feHist is null, exiting function.");
    return;
  }
  if (!canvas) {
    B2WARNING("processFEHistogram: canvas is null, cannot draw histograms.");
    return;
  }
 
  /* Obtain plots necessary for FE Ratio plots */
  auto* numerator = findHist(m_histogramDirectoryName + "/" + histName + "_0");
  auto* denominator = findHist(m_histogramDirectoryName + "/" + histName + "_1");
 
  feHist->Reset();
  std::unique_ptr<TH1> feClone(static_cast<TH1*>(feHist->Clone())); // Clone feHist
 
  canvas->cd();
 
  if (denominator != nullptr && numerator != nullptr) {
    // Clone numerator and denominator into unique_ptrs to ensure they are deleted automatically
    auto tempNumerator = std::unique_ptr<TH1>(static_cast<TH1*>(numerator->Clone()));
    auto tempSum = std::unique_ptr<TH1>(static_cast<TH1*>(denominator->Clone()));
 
    // Add numerator to denominator to get the sum (denominator + numerator)
    tempSum->Add(numerator);
 
    // Divide numerator by (denominator + numerator) and store result in feHist
    feHist->Divide(tempNumerator.get(), tempSum.get(), 1.0, 1.0, "B");
    feHist->Draw();
 
    // Reference check
    TH1* ref = findRefHist(feHist->GetName(), "", ERefScaling::c_RefScaleNone);
    if (ref) {
      ref->Draw("hist,same");
      B2INFO("processFEHistogram: Found and drew reference histogram.");
    } else {
      B2WARNING("processFEHistogram: Reference histogram not found.");
    }
 
    canvas->Modified();
    canvas->Update();
    B2INFO("processFEHistogram: Updated canvas after first draw.");
 
    /* Delta component */
    // Use the latest available deltas, not only "updated in the same event".
    auto deltaDenom = getDelta(m_histogramDirectoryName, histName + "_1", 0, false);
    auto deltaNumer = getDelta(m_histogramDirectoryName, histName + "_0", 0, false);
 
    UpdateCanvas(canvas->GetName(), (feHist != nullptr));
    if (deltaNumer != nullptr && deltaDenom != nullptr) {
      B2INFO("DQMHistAnalysisKLM: FE Ratio Delta Num/Denom Entries is "
             << deltaNumer->GetEntries() << "/" << deltaDenom->GetEntries());
 
      // Clone deltaNumer and deltaDenom into unique_ptrs to create the delta sum
      auto deltaTempNumerator = std::unique_ptr<TH1>(static_cast<TH1*>(deltaNumer->Clone()));
      auto deltaTempSum = std::unique_ptr<TH1>(static_cast<TH1*>(deltaDenom->Clone()));
 
      // Add delta numerator and denominator
      deltaTempSum->Add(deltaNumer);
 
      // Divide deltaNumer by (deltaDenom + deltaNumer) and store in feClone
      feClone->Divide(deltaTempNumerator.get(), deltaTempSum.get(), 1.0, 1.0, "B");
      feClone->SetLineColor(kOrange);
      feClone->DrawCopy("SAME");
 
      canvas->Modified();
      canvas->Update();
    }
    // If the deltas are not available yet (e.g. before the first MinEvents window
    // has closed), silently skip drawing the delta overlay, consistent with the
    // other histogram handlers (e.g. processTimeHistogram).
  } else {
    B2WARNING("processFEHistogram: Skipped histogram processing due to missing numerator/denominator.");
  }
}

processPlaneHistogram() [1/2]

void processPlaneHistogram ( const std::string & histName, TLatex & latex )

private

Process histogram containing the number of hits in plane.

Parameters

[in]	histName	Histogram name.
[out]	latex	TLatex to draw messages.

Definition at line 629 of file DQMHistAnalysisKLM.cc.

{
  // Internally call the main function with a pointer to TLatex
  processPlaneHistogram(histName, &latex, nullptr);
}

processPlaneHistogram() [2/2]

void processPlaneHistogram ( const std::string & histName, TLatex * latex = nullptr, TH1 * histogram = nullptr )

private

Process histogram containing the number of hits in plane.

Parameters

[in]	histName	Histogram name.
[in]	histogram	Histogram.
[out]	latex	TLatex to draw messages.

Definition at line 501 of file DQMHistAnalysisKLM.cc.

{
  // If no histogram is provided, try to retrieve it by name
  if (!histogram) {
    histogram = findHist(m_histogramDirectoryName + "/" + histName);
    if (!histogram) {
      B2WARNING("KLM DQM histogram " + m_histogramDirectoryName + "/" << histName << " is not found.");
      return;
    }
  }
 
  TCanvas* canvas = findCanvas(m_histogramDirectoryName + "/c_" + histName);
  if (!canvas) {
    B2WARNING("KLM DQM histogram canvas " + m_histogramDirectoryName + "/c_" << histName << " is not found.");
    return;
  }
  // Set up the canvas
  canvas->Clear();
  canvas->cd();
  histogram->SetStats(false);
  histogram->Draw();
 
  // Overlay reference histogram if available
  TH1* ref = findRefHist(histogram->GetName(), "", ERefScaling::c_RefScaleEntries, histogram);
  if (ref) {
    ref->Draw("hist,same");
  }
 
  // Define variables for text and alarm positioning
  double xAlarm = 0.15, yAlarm = 0.8;
  int message_counter = 0;
  std::string name;
 
  // Determine if histogram is BKLM or EKLM and process accordingly
  bool isBKLM = histName.find("bklm") != std::string::npos;
 
  // Check if "fe" is in histName to adjust maximalLayer
  bool isFE = histName.find("fe") != std::string::npos;
 
  if (isBKLM) {
    // Use a different maximalLayer if "fe" is found in the name
    const int maximalLayer = isFE ? BKLMElementNumbers::getMaximalSectorGlobalNumber() / 8 :
                             BKLMElementNumbers::getMaximalLayerNumber();
 
    for (int sector = 0; sector < BKLMElementNumbers::getMaximalSectorGlobalNumber(); ++sector) {
      int bin = maximalLayer * sector + 1;
      double xLine = histogram->GetXaxis()->GetBinLowEdge(bin);
      double xText = histogram->GetXaxis()->GetBinLowEdge(bin + maximalLayer / 2);
      double yText = gPad->GetUymin() + 0.98 * (gPad->GetUymax() - gPad->GetUymin());
 
      if (sector > 0) {
        m_PlaneLine.DrawLine(xLine, gPad->GetUymin(), xLine, gPad->GetUymax());
      }
 
      name = (sector < BKLMElementNumbers::getMaximalSectorNumber() ? "BB" : "BF") +
             std::to_string(sector % BKLMElementNumbers::getMaximalSectorNumber());
      m_PlaneText.DrawText(xText, yText, name.c_str());
    }
    // Only update canvas status if latex is provided
    if (latex) {
      updateCanvasStatus(canvas, m_DeadBarrelModules, latex, message_counter, xAlarm, yAlarm);
    }
  } else {
    const double maximalLayer = EKLMElementNumbers::getMaximalLayerGlobalNumber();
    const double maxPlane = EKLMElementNumbers::getMaximalPlaneNumber() * EKLMElementNumbers::getMaximalSectorNumber();
 
    for (int layerGlobal = 1; layerGlobal <= maximalLayer; ++layerGlobal) {
      int bin = maxPlane * layerGlobal + 1;
      double xLine = histogram->GetXaxis()->GetBinLowEdge(bin);
      double xText = histogram->GetXaxis()->GetBinLowEdge(bin - maxPlane / 2);
      double yText = gPad->GetUymin() + 0.98 * (gPad->GetUymax() - gPad->GetUymin());
 
      if (layerGlobal < maximalLayer) {
        m_PlaneLine.DrawLine(xLine, gPad->GetUymin(), xLine, gPad->GetUymax());
      }
 
      int section, layer;
      m_EklmElementNumbers->layerNumberToElementNumbers(layerGlobal, &section, &layer);
      name = (section == EKLMElementNumbers::c_BackwardSection ? "B" : "F") + std::to_string(layer);
      m_PlaneText.DrawText(xText, yText, name.c_str());
    }
 
    // Only update canvas status if latex is provided
    if (latex) {
      updateCanvasStatus(canvas, m_DeadEndcapModules, latex, message_counter, xAlarm, yAlarm);
    }
  }
 
  // Display additional message if message threshold exceeded
  if (latex && message_counter > m_MessageThreshold) {
    std::string verbose_string = std::to_string(message_counter - m_MessageThreshold) + " more messages";
    latex->DrawLatexNDC(xAlarm, yAlarm, verbose_string.c_str());
  }
 
  canvas->Modified();
  canvas->Update();
}

processSpatial2DHitEndcapHistogram()

void processSpatial2DHitEndcapHistogram ( uint16_t section, TH2F * histogram, TCanvas * canvas )

private

Process spatial 2D hits histograms for endcap.

Parameters

[in]	section	Section.
[in]	histogram	Histogram.
[in]	canvas	Canvas.

Definition at line 422 of file DQMHistAnalysisKLM.cc.

{
  canvas->Clear();
  canvas->cd();
  histogram->SetStats(false);
  histogram->Draw("COLZ");
  /* Draw the lines only for the backward layers. */
  if (section == EKLMElementNumbers::c_ForwardSection) {
    m_2DHitsLine.DrawLine(-110, 80, -110, 190);
    m_2DHitsLine.DrawLine(-110, 190, 110, 190);
    m_2DHitsLine.DrawLine(110, 80, 110, 190);
    m_2DHitsLine.DrawLine(-110, 80, 110, 80);
  }
  canvas->Modified();
  canvas->Update();
}

processTimeHistogram()

void processTimeHistogram ( const std::string & histName )

private

Process histogram containing the hit times.

Parameters

[in]

histName

Histogram name.

Definition at line 440 of file DQMHistAnalysisKLM.cc.

{
  TH1* histogram = findHist(m_histogramDirectoryName + "/" + histName);
  if (histogram == nullptr) {
    B2WARNING("KLM DQM histogram " + m_histogramDirectoryName + "/" << histName << " is not found.");
    return;
  }
 
  TCanvas* canvas = findCanvas(m_histogramDirectoryName + "/c_" + histName);
  if (canvas == nullptr) {
    B2WARNING("KLM DQM histogram canvas " + m_histogramDirectoryName + "/c_" << histName << " is not found.");
    return;
  }
 
  else {
    canvas->Clear();
    canvas->cd();
    histogram->Draw();
    /* calling on delta histogram*/
    TH1* delta = getDelta(m_histogramDirectoryName, histName);
    UpdateCanvas(canvas->GetName(), delta != nullptr); //keeping this for testing purposes
    if (delta != nullptr) {
      B2INFO("DQMHistAnalysisKLM: Time Delta Entries is " << delta->GetEntries());
      deltaDrawer(delta, histogram, canvas);
    }
    //reference check
    TH1* ref = findRefHist(histogram->GetName(), "", ERefScaling::c_RefScaleEntries, histogram);
    if (ref) {ref->Draw("hist,same");}
  }
}

registerEpicsPV()

int registerEpicsPV ( const std::string & pvname, const std::string & keyname = "" )

inherited

EPICS related Functions.

Register a PV with its name and a key name

Parameters

pvname	full PV name
keyname	key name for easier access

Returns

an index which can be used to access the PV instead of key name, -1 if failure

Definition at line 373 of file DQMHistAnalysis.cc.

{
  return registerEpicsPVwithPrefix(m_PVPrefix, pvname, keyname);
}

registerEpicsPVwithPrefix()

int registerEpicsPVwithPrefix ( const std::string & prefix, const std::string & pvname, const std::string & keyname = "" )

privateinherited

Register a PV with its name and a key name.

Parameters

prefix	prefix to PV name
pvname	full PV name without prefix
keyname	key name for easier access

Returns

an index which can be used to access the PV instead of key name, -1 if failure

Definition at line 383 of file DQMHistAnalysis.cc.

{
  if (!m_useEpics) return -1;
#ifdef _BELLE2_EPICS
  if (m_epicsNameToChID[pvname] != nullptr) {
    B2ERROR("Epics PV " << pvname << " already registered!");
    return -1;
  }
  if (keyname != "" && m_epicsNameToChID[keyname] != nullptr) {
    B2ERROR("Epics PV with key " << keyname << " already registered!");
    return -1;
  }
 
  m_epicsChID.emplace_back();
  auto ptr = &m_epicsChID.back();
  if (!ca_current_context()) SEVCHK(ca_context_create(ca_disable_preemptive_callback), "ca_context_create");
  // the subscribed name includes the prefix, the map below does *not*
  CheckEpicsError(ca_create_channel((prefix + pvname).data(), NULL, NULL, 10, ptr), "ca_create_channel failure", pvname);
 
  m_epicsNameToChID[pvname] =  *ptr;
  if (keyname != "") m_epicsNameToChID[keyname] =  *ptr;
  return m_epicsChID.size() - 1; // return index to last added item
#else
  return -1;
#endif
}

registerExternalEpicsPV()

int registerExternalEpicsPV ( const std::string & pvname, const std::string & keyname = "" )

inherited

Register a PV with its name and a key name.

Parameters

pvname	full PV name
keyname	key name for easier access

Returns

an index which can be used to access the PV instead of key name, -1 if failure

Definition at line 378 of file DQMHistAnalysis.cc.

{
  return registerEpicsPVwithPrefix(std::string(""), pvname, keyname);
}

requestLimitsFromEpicsPVs() [1/3]

bool requestLimitsFromEpicsPVs ( chid id, double & lowerAlarm, double & lowerWarn, double & upperWarn, double & upperAlarm )

inherited

Get Alarm Limits from EPICS PV.

Parameters

id	Channel ID
&lowerAlarm	return low Alarm limit (lolo) if set, not changed otherwise
&lowerWarn	return low Warning limit (low) if set, not changed otherwise
&upperWarn	return upper Warning limit (high) if set, not changed otherwise
&upperAlarm	return upper Alarm limit (hihi) if set, not changed otherwise

Returns

true if limits could be read (even if there are none set)

Definition at line 655 of file DQMHistAnalysis.cc.

{
  // get warn and error limit only if pv exists
  // overwrite only if limit is defined (not NaN)
  // user should initialize with NaN before calling, unless
  // some "default" values should be set otherwise
  if (pv != nullptr) {
    struct dbr_ctrl_double tPvData;
    // From EPICS doc. When ca_get or ca_array_get are invoked the returned channel value can't be assumed to be stable
    // in the application supplied buffer until after ECA_NORMAL is returned from ca_pend_io. If a connection is lost
    // outstanding get requests are not automatically reissued following reconnect.
    auto r = ca_get(DBR_CTRL_DOUBLE, pv, &tPvData);
    if (r == ECA_NORMAL) r = ca_pend_io(5.0); // this is needed!
    if (r == ECA_NORMAL) {
      if (!std::isnan(tPvData.lower_alarm_limit)) {
        lowerAlarm = tPvData.lower_alarm_limit;
      }
      if (!std::isnan(tPvData.lower_warning_limit)) {
        lowerWarn = tPvData.lower_warning_limit;
      }
      if (!std::isnan(tPvData.upper_warning_limit)) {
        upperWarn = tPvData.upper_warning_limit;
      }
      if (!std::isnan(tPvData.upper_alarm_limit)) {
        upperAlarm = tPvData.upper_alarm_limit;
      }
      return true;
    } else {
      CheckEpicsError(r, "Reading PV Limits failed in ca_get or ca_pend_io failure", pv);
    }
  }
  return false;
}

requestLimitsFromEpicsPVs() [2/3]

bool requestLimitsFromEpicsPVs ( const std::string & keyname, double & lowerAlarm, double & lowerWarn, double & upperWarn, double & upperAlarm )

inherited

Get Alarm Limits from EPICS PV.

Parameters

keyname	key name (or full PV name) of PV
&lowerAlarm	return low Alarm limit (lolo) if set, not changed otherwise
&lowerWarn	return low Warning limit (low) if set, not changed otherwise
&upperWarn	return upper Warning limit (high) if set, not changed otherwise
&upperAlarm	return upper Alarm limit (hihi) if set, not changed otherwise

Returns

true if limits could be read (even if there are none set)

Definition at line 643 of file DQMHistAnalysis.cc.

{
  return requestLimitsFromEpicsPVs(getEpicsPVChID(name), lowerAlarm, lowerWarn, upperWarn, upperAlarm);
}

requestLimitsFromEpicsPVs() [3/3]

bool requestLimitsFromEpicsPVs ( int index, double & lowerAlarm, double & lowerWarn, double & upperWarn, double & upperAlarm )

inherited

Get Alarm Limits from EPICS PV.

Parameters

index	index of PV
&lowerAlarm	return low Alarm limit (lolo) if set, not changed otherwise
&lowerWarn	return low Warning limit (low) if set, not changed otherwise
&upperWarn	return upper Warning limit (high) if set, not changed otherwise
&upperAlarm	return upper Alarm limit (hihi) if set, not changed otherwise

Returns

true if limits could be read (even if there are none set)

Definition at line 649 of file DQMHistAnalysis.cc.

{
  return requestLimitsFromEpicsPVs(getEpicsPVChID(index), lowerAlarm, lowerWarn, upperWarn, upperAlarm);
}

resetDeltaList()

void resetDeltaList ( void )

staticinherited

Reset Delta.

Definition at line 332 of file DQMHistAnalysis.cc.

{
  for (auto& d : s_deltaList) {
    d.second.reset();
  }
}

scaleReference()

TH1 * scaleReference ( ERefScaling scaling, const TH1 * hist, TH1 * ref )

staticinherited

Using the original and reference, create scaled version.

Parameters

scaling	scaling algorithm
hist	pointer to histogram
ref	pointer to reference

Returns

scaled reference

Definition at line 171 of file DQMHistAnalysis.cc.

{
  // if hist/ref is nullptr, nothing to do
  if (!hist || !ref)
    return ref;
 
  switch (scaling) {
    // default: do nothing
    case ERefScaling::c_RefScaleNone: //do nothing
      break;
    case ERefScaling::c_RefScaleEntries: // Integral
      // only if we have entries in reference
      if (hist->Integral() != 0 and ref->Integral() != 0) {
        ref->Scale(hist->Integral() / ref->Integral());
      }
      break;
    case ERefScaling::c_RefScaleMax: // Maximum
      // only if we have entries in reference
      if (hist->GetMaximum() != 0 and ref->GetMaximum() != 0) {
        ref->Scale(hist->GetMaximum() / ref->GetMaximum());
      }
      break;
  }
  return ref;
}

setAbortLevel()

void setAbortLevel ( int abortLevel )

inherited

Configure the abort log level.

Definition at line 67 of file Module.cc.

{
  m_logConfig.setAbortLevel(static_cast<LogConfig::ELogLevel>(abortLevel));
}

setDebugLevel()

void setDebugLevel ( int debugLevel )

inherited

Configure the debug messaging level.

Definition at line 61 of file Module.cc.

{
  m_logConfig.setDebugLevel(debugLevel);
}

setDescription()

void setDescription ( const std::string & description )

protectedinherited

Sets the description of the module.

Parameters

description

A description of the module.

Definition at line 214 of file Module.cc.

{
  m_description = description;
}

setEpicsPV() [1/4]

void setEpicsPV ( const std::string & keyname, double value )

inherited

Write value to a EPICS PV.

Parameters

keyname	key name (or full PV name) of PV
value	value to write

Definition at line 411 of file DQMHistAnalysis.cc.

{
  if (!m_useEpics || m_epicsReadOnly) return;
#ifdef _BELLE2_EPICS
  if (m_epicsNameToChID[keyname] == nullptr) {
    B2ERROR("Epics PV " << keyname << " not registered!");
    return;
  }
  CheckEpicsError(ca_put(DBR_DOUBLE, m_epicsNameToChID[keyname], &value), "ca_set failure", keyname);
#endif
}

setEpicsPV() [2/4]

void setEpicsPV ( const std::string & keyname, int value )

inherited

Write value to a EPICS PV.

Parameters

keyname	key name (or full PV name) of PV
value	value to write

Definition at line 423 of file DQMHistAnalysis.cc.

{
  if (!m_useEpics || m_epicsReadOnly) return;
#ifdef _BELLE2_EPICS
  if (m_epicsNameToChID[keyname] == nullptr) {
    B2ERROR("Epics PV " << keyname << " not registered!");
    return;
  }
  CheckEpicsError(ca_put(DBR_SHORT, m_epicsNameToChID[keyname], &value), "ca_set failure", keyname);
#endif
}

setEpicsPV() [3/4]

void setEpicsPV ( int index, double value )

inherited

Write value to a EPICS PV.

Parameters

index	index of PV
value	value to write

Definition at line 453 of file DQMHistAnalysis.cc.

{
  if (!m_useEpics || m_epicsReadOnly) return;
#ifdef _BELLE2_EPICS
  if (index < 0 || index >= (int)m_epicsChID.size()) {
    B2ERROR("Epics PV with " << index << " not registered!");
    return;
  }
  CheckEpicsError(ca_put(DBR_DOUBLE, m_epicsChID[index], &value), "ca_set failure", m_epicsChID[index]);
#endif
}

setEpicsPV() [4/4]

void setEpicsPV ( int index, int value )

inherited

Write value to a EPICS PV.

Parameters

index	index of PV
value	value to write

Definition at line 465 of file DQMHistAnalysis.cc.

{
  if (!m_useEpics || m_epicsReadOnly) return;
#ifdef _BELLE2_EPICS
  if (index < 0 || index >= (int)m_epicsChID.size()) {
    B2ERROR("Epics PV with " << index << " not registered!");
    return;
  }
  CheckEpicsError(ca_put(DBR_SHORT, m_epicsChID[index], &value), "ca_set failure", m_epicsChID[index]);
#endif
}

setEpicsStringPV() [1/2]

void setEpicsStringPV ( const std::string & keyname, const std::string & value )

inherited

Write string to a EPICS PV.

Parameters

keyname	key name (or full PV name) of PV
value	string to write

Definition at line 435 of file DQMHistAnalysis.cc.

{
  if (!m_useEpics || m_epicsReadOnly) return;
#ifdef _BELLE2_EPICS
  if (m_epicsNameToChID[keyname] == nullptr) {
    B2ERROR("Epics PV " << keyname << " not registered!");
    return;
  }
  if (value.length() > 40) {
    B2ERROR("Epics string PV " << keyname << " too long (>40 characters)!");
    return;
  }
  char text[40];
  strcpy(text, value.c_str());
  CheckEpicsError(ca_put(DBR_STRING, m_epicsNameToChID[keyname], text), "ca_set failure", keyname);
#endif
}

setEpicsStringPV() [2/2]

void setEpicsStringPV ( int index, const std::string & value )

inherited

Write string to a EPICS PV.

Parameters

index	index of PV
value	string to write

Definition at line 477 of file DQMHistAnalysis.cc.

{
  if (!m_useEpics || m_epicsReadOnly) return;
#ifdef _BELLE2_EPICS
  if (index < 0 || index >= (int)m_epicsChID.size()) {
    B2ERROR("Epics PV with " << index << " not registered!");
    return;
  }
  char text[41];
  strncpy(text, value.c_str(), 40);
  text[40] = 0;
  CheckEpicsError(ca_put(DBR_STRING, m_epicsChID[index], text), "ca_set failure", m_epicsChID[index]);
#endif
}

setEventProcessed()

static void setEventProcessed ( int e )

inlinestaticinherited

Set the number of processed events.

(Attention, asynch histogram updates!)

Parameters

e	Processed events.

Definition at line 237 of file DQMHistAnalysis.h.

{s_eventProcessed = e;};

setLogConfig()

void setLogConfig ( const LogConfig & logConfig )

inlineinherited

Set the log system configuration.

Definition at line 229 of file Module.h.

{m_logConfig = logConfig;}

setLogInfo()

void setLogInfo ( int logLevel, unsigned int logInfo )

inherited

Configure the printed log information for the given level.

Parameters

logLevel	The log level (one of LogConfig::ELogLevel)
logInfo	What kind of info should be printed? ORed combination of LogConfig::ELogInfo flags.

Definition at line 73 of file Module.cc.

{
  m_logConfig.setLogInfo(static_cast<LogConfig::ELogLevel>(logLevel), logInfo);
}

setLogLevel()

void setLogLevel ( int logLevel )

inherited

Configure the log level.

Definition at line 55 of file Module.cc.

{
  m_logConfig.setLogLevel(static_cast<LogConfig::ELogLevel>(logLevel));
}

setName()

void setName ( const std::string & name )

inlineinherited

Set the name of the module.

Note

The module name is set when using the REG_MODULE macro, but the module can be renamed before calling process() using the set_name() function in your steering file.

Parameters

name	The name of the module

Definition at line 213 of file Module.h.

{ m_name = name; };

setParamList()

void setParamList ( const ModuleParamList & params )

inlineprotectedinherited

Replace existing parameter list.

Definition at line 500 of file Module.h.

{ m_moduleParamList = params; }

setParamPython()

void setParamPython ( const std::string & name, const boost::python::object & pyObj )

privateinherited

Implements a method for setting boost::python objects.

The method supports the following types: list, dict, int, double, string, bool The conversion of the python object to the C++ type and the final storage of the parameter value is done in the ModuleParam class.

Parameters

name	The unique name of the parameter.
pyObj	The object which should be converted and stored as the parameter value.

Definition at line 234 of file Module.cc.

{
  LogSystem& logSystem = LogSystem::Instance();
  logSystem.updateModule(&(getLogConfig()), getName());
  try {
    m_moduleParamList.setParamPython(name, pyObj);
  } catch (std::runtime_error& e) {
    throw std::runtime_error("Cannot set parameter '" + name + "' for module '"
                             + m_name + "': " + e.what());
  }
 
  logSystem.updateModule(nullptr);
}

setParamPythonDict()

void setParamPythonDict ( const boost::python::dict & dictionary )

privateinherited

Implements a method for reading the parameter values from a boost::python dictionary.

The key of the dictionary has to be the name of the parameter and the value has to be of one of the supported parameter types.

Parameters

dictionary

The python dictionary from which the parameter values are read.

Definition at line 249 of file Module.cc.

{
 
  LogSystem& logSystem = LogSystem::Instance();
  logSystem.updateModule(&(getLogConfig()), getName());
 
  boost::python::list dictKeys = dictionary.keys();
  int nKey = boost::python::len(dictKeys);
 
  //Loop over all keys in the dictionary
  for (int iKey = 0; iKey < nKey; ++iKey) {
    boost::python::object currKey = dictKeys[iKey];
    boost::python::extract<std::string> keyProxy(currKey);
 
    if (keyProxy.check()) {
      const boost::python::object& currValue = dictionary[currKey];
      setParamPython(keyProxy, currValue);
    } else {
      B2ERROR("Setting the module parameters from a python dictionary: invalid key in dictionary!");
    }
  }
 
  logSystem.updateModule(nullptr);
}

setPropertyFlags()

void setPropertyFlags ( unsigned int propertyFlags )

inherited

Sets the flags for the module properties.

Parameters

propertyFlags

bitwise OR of EModulePropFlags

Definition at line 208 of file Module.cc.

{
  m_propertyFlags = propertyFlags;
}

setPVPrefix()

static void setPVPrefix ( const std::string & prefix )

inlinestaticinherited

set global prefix for EPICS PVs

Parameters

prefix

Prefix to set

Definition at line 594 of file DQMHistAnalysis.h.

{ m_PVPrefix = prefix;};

setReturnValue() [1/2]

void setReturnValue ( bool value )

protectedinherited

Sets the return value for this module as bool.

The bool value is saved as an integer with the convention 1 meaning true and 0 meaning false. The value can be used in the steering file to divide the analysis chain into several paths.

Parameters

value

The value of the return value.

Definition at line 227 of file Module.cc.

{
  m_hasReturnValue = true;
  m_returnValue = value;
}

setReturnValue() [2/2]

void setReturnValue ( int value )

protectedinherited

Sets the return value for this module as integer.

The value can be used in the steering file to divide the analysis chain into several paths.

Parameters

value

The value of the return value.

Definition at line 220 of file Module.cc.

{
  m_hasReturnValue = true;
  m_returnValue = value;
}

setRunType()

static void setRunType ( const std::string & t )

inlinestaticinherited

Set the Run Type.

Parameters

t	Run type string.

Definition at line 231 of file DQMHistAnalysis.h.

{s_runType = t;};

setType()

void setType ( const std::string & type )

protectedinherited

Set the module type.

Only for use by internal modules (which don't use the normal REG_MODULE mechanism).

Definition at line 48 of file Module.cc.

{
  if (!m_type.empty())
    B2FATAL("Trying to change module type from " << m_type << " is not allowed, the value is assumed to be fixed.");
  m_type = type;
}

setUseEpics()

static void setUseEpics ( bool flag )

inlinestaticinherited

Setter for EPICS usage.

Parameters

flag	set in use

Definition at line 559 of file DQMHistAnalysis.h.

{m_useEpics = flag;};

setUseEpicsReadOnly()

static void setUseEpicsReadOnly ( bool flag )

inlinestaticinherited

Setter EPICS flag in read only mode.

Parameters

flag	set read only

Definition at line 565 of file DQMHistAnalysis.h.

{m_epicsReadOnly = flag;};

StringSplit()

std::vector< std::string > StringSplit ( const std::string & s, const char delim )

staticinherited

Helper function for string token split.

Parameters

s	String to split
delim	delimiter

Returns

vector of strings

Definition at line 287 of file DQMHistAnalysis.cc.

{
  std::vector <std::string> out;
  boost::split(out, in, [delim](char c) {return c == delim;});
  return out;
}

terminate()

void terminate ( void )

finaloverridevirtual

This method is called at the end of the event processing.

Reimplemented from Module.

Definition at line 153 of file DQMHistAnalysisKLM.cc.

{
}

UpdateCanvas() [1/2]

void UpdateCanvas ( const std::string & name, bool updated = true )

staticinherited

Mark canvas as updated (or not)

Parameters

name	name of Canvas
updated	was updated

Definition at line 339 of file DQMHistAnalysis.cc.

{
  s_canvasUpdatedList[name] = updated;
}

UpdateCanvas() [2/2]

void UpdateCanvas ( TCanvas * canvas, bool updated = true )

staticinherited

Mark canvas as updated (or not)

Parameters

canvas	Canvas from which to take the name for update
updated	was updated

Definition at line 344 of file DQMHistAnalysis.cc.

{
  if (c) UpdateCanvas(c->GetName(), updated);
}

updateCanvasStatus()

void updateCanvasStatus ( TCanvas * canvas, const std::vector< KLMModuleNumber > & deadModules, TLatex * latex, int & message_counter, double xAlarm, double yAlarm )

private

Helper function to update the canvas status based on dead modules.

Parameters

[in]	canvas	Canvas of interest.
[in]	deadModules	List of dead modules.
[in]	latex	LaTeX object to render text annotations.
[out]	message_counter	Counter for how many messages are drawn.
[out]	xAlarm	X-position for the alarm messages.
[out]	yAlarm	Y-position for the alarm messages.

Definition at line 601 of file DQMHistAnalysisKLM.cc.

{
  if (deadModules.empty()) {
    colorizeCanvas(canvas, c_StatusGood);
  } else if (m_ProcessedEvents > m_MinProcessedEventsForMessages) {
    for (const KLMModuleNumber& module : deadModules) {
      int moduleSubdetector, moduleSection, moduleSector, moduleLayer;
      m_ElementNumbers->moduleNumberToElementNumbers(module, &moduleSubdetector, &moduleSection, &moduleSector, &moduleLayer);
 
      std::string alarm = "No data from " + m_ElementNumbers->getSectorDAQName(moduleSubdetector, moduleSection, moduleSector)
                          + ", layer " + std::to_string(moduleLayer);
      message_counter++;
      if (message_counter <= m_MessageThreshold) {
        latex->DrawLatexNDC(xAlarm, yAlarm, alarm.c_str());
        yAlarm -= 0.05;
      }
    }
    if (!m_IsNullRun) {
      colorizeCanvas(canvas, c_StatusError);
    }
  } else {
    colorizeCanvas(canvas, c_StatusTooFew);
  }
}

updateEpicsPVs()

int updateEpicsPVs ( float timeout )

inherited

Update all EPICS PV (flush to network)

Parameters

timeout

maximum time until timeout in s

Returns

status of ca_pend_io

Definition at line 616 of file DQMHistAnalysis.cc.

{
  int state = ECA_NORMAL;
  if (!m_useEpics) return state;
#ifdef _BELLE2_EPICS
  if (wait > 0.) {
    state = ca_pend_io(wait);
    SEVCHK(state, "ca_pend_io failure");
  }
#endif
  return state;
}

Member Data Documentation

m_2DHitsLine

TLine m_2DHitsLine

private

TLine for background region in 2d hits histograms.

Definition at line 217 of file DQMHistAnalysisKLM.h.

m_c_fe_bklm_ratio

TCanvas* m_c_fe_bklm_ratio = nullptr

private

Canvas for BKLM plane events fraction w/ FE.

Definition at line 229 of file DQMHistAnalysisKLM.h.

m_c_fe_eklm_ratio

TCanvas* m_c_fe_eklm_ratio = nullptr

private

Canvas for EKLM plane events fraction w/ FE.

Definition at line 235 of file DQMHistAnalysisKLM.h.

m_c_masked_channels

TCanvas* m_c_masked_channels = nullptr

private

Canvas for masked channels per sector summary.

Definition at line 214 of file DQMHistAnalysisKLM.h.

m_ChannelArrayIndex

const KLMChannelArrayIndex* m_ChannelArrayIndex

private

KLM channel array index.

Definition at line 241 of file DQMHistAnalysisKLM.h.

m_conditions

std::vector<ModuleCondition> m_conditions

privateinherited

Module condition, only non-null if set.

Definition at line 520 of file Module.h.

m_DeadBarrelModules

std::vector<uint16_t> m_DeadBarrelModules

private

Vector of dead barrel modules.

Definition at line 202 of file DQMHistAnalysisKLM.h.

m_DeadEndcapModules

std::vector<uint16_t> m_DeadEndcapModules

private

Vector of dead endcap modules.

Definition at line 205 of file DQMHistAnalysisKLM.h.

m_description

std::string m_description

privateinherited

The description of the module.

Definition at line 510 of file Module.h.

m_EklmElementNumbers

const EKLMElementNumbers* m_EklmElementNumbers

private

EKLM element numbers.

Definition at line 250 of file DQMHistAnalysisKLM.h.

m_ElectronicsMap

DBObjPtr<KLMElectronicsMap> m_ElectronicsMap

private

Electronics map.

Definition at line 253 of file DQMHistAnalysisKLM.h.

m_ElementNumbers

const KLMElementNumbers* m_ElementNumbers

private

KLM element numbers.

Definition at line 247 of file DQMHistAnalysisKLM.h.

m_epicsReadOnly

bool m_epicsReadOnly

staticprivateinherited

Initial value:

=
  false

Flag if to use EPICS in ReadOnly mode (for reading limits) do not set by yourself, use EpicsEnable module to set.

Definition at line 159 of file DQMHistAnalysis.h.

m_fe_bklm_ratio

TH1* m_fe_bklm_ratio = nullptr

private

Histogram for BKLM plane events fraction w/ FE.

Definition at line 226 of file DQMHistAnalysisKLM.h.

m_fe_eklm_ratio

TH1* m_fe_eklm_ratio = nullptr

private

Histogram for EKLM plane events fraction w/ FE.

Definition at line 232 of file DQMHistAnalysisKLM.h.

m_hasReturnValue

bool m_hasReturnValue

privateinherited

True, if the return value is set.

Definition at line 517 of file Module.h.

m_histogramDirectoryName

std::string m_histogramDirectoryName

private

Name of histogram directory.

Definition at line 199 of file DQMHistAnalysisKLM.h.

m_IsNullRun

bool m_IsNullRun

private

Run type flag for null runs.

Definition at line 238 of file DQMHistAnalysisKLM.h.

m_logConfig

LogConfig m_logConfig

privateinherited

The log system configuration of the module.

Definition at line 513 of file Module.h.

m_MaskedChannels

std::vector<uint16_t> m_MaskedChannels

private

Vector of masked channels.

Definition at line 208 of file DQMHistAnalysisKLM.h.

m_MaskedChannelsHist

TH1* m_MaskedChannelsHist = nullptr

private

Histogram for masked channels per sector.

Definition at line 211 of file DQMHistAnalysisKLM.h.

m_MessageThreshold

int m_MessageThreshold

private

Message Threshold for expert pots.

Definition at line 184 of file DQMHistAnalysisKLM.h.

m_minEntries

double m_minEntries

private

Minimal number of entries for delta histogram update.

Definition at line 193 of file DQMHistAnalysisKLM.h.

m_minEvents

double m_minEvents

private

Minimal number of processed events for delta histogram update.

Definition at line 196 of file DQMHistAnalysisKLM.h.

m_MinHitsForFlagging

int m_MinHitsForFlagging

private

Minimal number of hits for flagging.

Definition at line 181 of file DQMHistAnalysisKLM.h.

m_MinProcessedEventsForMessages

double m_MinProcessedEventsForMessages

private

Minimal number of processed events for error messages.

Definition at line 190 of file DQMHistAnalysisKLM.h.

m_MinProcessedEventsForMessagesInput

double m_MinProcessedEventsForMessagesInput

private

Input parameter for minimal number of processed events for error messages.

Definition at line 187 of file DQMHistAnalysisKLM.h.

m_moduleParamList

ModuleParamList m_moduleParamList

privateinherited

List storing and managing all parameter of the module.

Definition at line 515 of file Module.h.

m_monObj

MonitoringObject* m_monObj {}

private

Monitoring object.

Definition at line 256 of file DQMHistAnalysisKLM.h.

{};

m_name

std::string m_name

privateinherited

The name of the module, saved as a string (user-modifiable)

Definition at line 507 of file Module.h.

m_package

std::string m_package

privateinherited

Package this module is found in (may be empty).

Definition at line 509 of file Module.h.

m_PlaneLine

TLine m_PlaneLine

private

TLine for boundary in plane histograms.

Definition at line 220 of file DQMHistAnalysisKLM.h.

m_PlaneText

TText m_PlaneText

private

TText for names in plane histograms.

Definition at line 223 of file DQMHistAnalysisKLM.h.

m_ProcessedEvents

double m_ProcessedEvents

private

Number of processed events.

Definition at line 169 of file DQMHistAnalysisKLM.h.

m_propertyFlags

unsigned int m_propertyFlags

privateinherited

The properties of the module as bitwise or (with |) of EModulePropFlags.

Definition at line 511 of file Module.h.

m_PVPrefix

std::string m_PVPrefix = "TEST:"

staticprivateinherited

The Prefix for EPICS PVs.

Definition at line 164 of file DQMHistAnalysis.h.

m_returnValue

int m_returnValue

privateinherited

The return value.

Definition at line 518 of file Module.h.

m_SectorArrayIndex

const KLMSectorArrayIndex* m_SectorArrayIndex

private

KLM sector array index.

Definition at line 244 of file DQMHistAnalysisKLM.h.

m_ThresholdForHot

int m_ThresholdForHot

private

Threshold for hot channels.

Definition at line 175 of file DQMHistAnalysisKLM.h.

m_ThresholdForLog

int m_ThresholdForLog

private

Threshold for log scale.

Definition at line 178 of file DQMHistAnalysisKLM.h.

m_ThresholdForMasked

int m_ThresholdForMasked

private

Threshold for masked channels.

Definition at line 172 of file DQMHistAnalysisKLM.h.

m_type

std::string m_type

privateinherited

The type of the module, saved as a string.

Definition at line 508 of file Module.h.

m_useEpics

bool m_useEpics = false

staticprivateinherited

Flag if to use EPICS do not set by yourself, use EpicsEnable module to set.

Definition at line 153 of file DQMHistAnalysis.h.

s_canvasUpdatedList

DQMHistAnalysisModule::CanvasUpdatedList s_canvasUpdatedList

staticprivateinherited

The list of canvas updated status.

Definition at line 128 of file DQMHistAnalysis.h.

s_deltaList

DQMHistAnalysisModule::DeltaList s_deltaList

staticinherited

The list of Delta Histograms and settings.

Definition at line 122 of file DQMHistAnalysis.h.

s_eventProcessed

int s_eventProcessed = 0

inlinestaticprivateinherited

Number of Events processed to fill histograms.

Attention: histograms are updates asynchronously Thus the number for a specific histogram may be lower or higher. If you need precise number, you must fill it in the histogram itself (e.g. underflow bin)

Definition at line 142 of file DQMHistAnalysis.h.

s_histList

DQMHistAnalysisModule::HistList s_histList

staticprivateinherited

The list of Histograms.

Definition at line 112 of file DQMHistAnalysis.h.

s_monObjList

DQMHistAnalysisModule::MonObjList s_monObjList

staticprivateinherited

The list of MonitoringObjects.

Definition at line 116 of file DQMHistAnalysis.h.

s_refList

DQMHistAnalysisModule::RefList s_refList

staticprivateinherited

The list of references.

Definition at line 133 of file DQMHistAnalysis.h.

s_runType

std::string s_runType = ""

inlinestaticprivateinherited

The Run type.

Definition at line 147 of file DQMHistAnalysis.h.

---

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

• dqm/analysis/modules/include/DQMHistAnalysisKLM.h
• dqm/analysis/modules/src/DQMHistAnalysisKLM.cc