DQMHistAnalysisTOPModule Class Reference

Class for TOP histogram analysis. More...

#include <DQMHistAnalysisTOP.h>

Inheritance diagram for DQMHistAnalysisTOPModule:

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
	DQMHistAnalysisTOPModule ()
	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.
void	clearlist (void)
	Clear all static global lists.
TH1 *	findHistInCanvas (const std::string &hname, TCanvas **canvas=nullptr)
	Find histogram in corresponding canvas.
TH1 *	getDelta (const std::string &fullname, int n=0, bool onlyIfUpdated=true)
	Get Delta histogram.
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, 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	setUseEpics (bool flag)
	Setter for EPICS usage.
void	cleanupEpicsPVs (void)
	Unsubscribe from EPICS PVs on terminate.
void	checkPVStatus (void)
	Check the status of all PVs and report if disconnected or not found.
void	CheckEpicsError (int state, const std::string &message, const std::string &name)
	check the return status and check PV in case of error
void	CheckEpicsError (int state, const std::string &message, chid id)
	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 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 &histname, bool onlyIfUpdated=false)
	Get histogram from list (no other search).
static TH1 *	findHist (const std::string &dirname, const std::string &histname, bool onlyIfUpdated=false)
	Find histogram.
static TH1 *	findRefHist (const std::string &histname, ERefScaling scaling=ERefScaling::c_RefScaleNone, const TH1 *hist=nullptr)
	Get referencehistogram from list (no other search).
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	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	printPVStatus (chid pv, bool onlyError=true)
	check the status of a PVs and report if disconnected or not found
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 Types
enum	EAlarmStates {   c_Gray = 0 ,   c_Green = 1 ,   c_Yellow = 2 ,   c_Red = 3 }
	Alarm states enumerator. More...

Private Member Functions
void	updateWindowVsSlotCanvas ()
	Updates canvas of window_vs_slot w/ alarming.
void	updateEventMonitorCanvas ()
	Updates canvas of event desynchronization monitor w/ alarming.
void	updateNGoodHitsCanvas ()
	Updates canvas of number of good hits per event w/ alarming (injection BG)
void	updateEventT0Canvas ()
	Updates canvas of event T0 w/ alarming.
void	updateBunchOffsetCanvas ()
	Updates canvas of bunch offset w/ alarming.
void	updateTimingCanvas ()
	Updates canvas of timing plot w/ alarming.
const TH1F *	makeDeadAndHotFractionsPlot ()
	Makes a plot of dead and hot channel fractions per slot.
void	makePhotonYieldsAndBGRatePlots (const TH1F *activeFraction)
	Make plots of dead-and-hot-channel corrected photon yields and BG rates per slot.
void	makeJunkFractionPlot ()
	Makes a plot of fractions of junk hits per boardstack.
void	makePMTHitRatesPlots ()
	Makes plots of the number of PMT hits per event.
void	makeInjectionBGPlots ()
	Makes projections of injection BG plots.
void	setMiraBelleVariables (const std::string &variableName, const TH1 *histogram)
	Sets MiraBelle variables from the histogram with bins corresponding to slot numbers.
int	getAlarmState (double value, const std::vector< double > &alarmLevels, bool bigRed=true) const
	Returns alarm state.
int	getAlarmColor (unsigned alarmState) const
	Converts alarm state to color.
int	getOffcialAlarmStatus (unsigned alarmState) const
	Converts alarm state to official status (see EStatus of the base class)
void	setAlarmLines ()
	Sets all alarm lines.
void	setEpicsVariables ()
	Calculates and sets epics variables.
void	updateLimits ()
	Updates limits defined by module parameters using EpicsPVs.
void	setIncludedBoardstacks (const std::vector< std::string > &excludedBoardstacks)
	Sets flags for boardstacks to be included in alarming.
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 bool	sameHistDefinition (TH1 *h1, TH1 *h2)
	Checks if histograms are defined in the same way (nbins, xmin, xmax)
static void	setZAxisRange (const std::string &name, double scale)
	Sets z-axis range of 2D histograms.
static void	makeBGSubtractedTimingPlot (const std::string &name, const TH2F *trackHits, int slot)
	Makes background subtracted time distribution plot.
static void	makeFlagFractPlot (const std::string &hname, TH1 *histogram, TCanvas *canvas)
	Makes a plot of fraction of events with the flag is set.
static void	setAlarmLines (const std::vector< double > &alarmLevels, double xmin, double xmax, std::vector< TLine * > &alarmLines, bool bigRed=true)
	Sets alarm lines.
static std::pair< double, double >	getDeadAndHotCuts (const TH1 *h)
	Returns cut levels for dead and hot channels.

Private Attributes
std::vector< int >	m_asicWindowsBand = {215, 235}
	lower and upper bin of a band denoting good windows
std::vector< double >	m_asicWindowsAlarmLevels = {0.002, 0.02}
	alarm levels for fraction of windows outside the band
std::vector< double >	m_eventMonitorAlarmLevels = {1e-4, 2e-3}
	alarm levels for fraction of desynchronized digits
std::vector< double >	m_junkHitsAlarmLevels = {0.05, 0.25}
	alarm levels for the fraction of junk hits
std::vector< double >	m_deadChannelsAlarmLevels = {0.1, 0.35}
	alarm levels for the fraction of dead + hot channels
std::vector< double >	m_backgroundAlarmLevels = {5.0, 10.0}
	alarm levels for background rates [MHz/PMT]
std::vector< double >	m_photonYieldsAlarmLevels = {15.0, 25.0}
	alarm levels for the number of photons per track
std::vector< std::string >	m_excludedBoardstacks
	list of boarstacks to be excluded from alarming
std::string	m_pvPrefix
	Epics PV prefix.
std::vector< double >	m_injectionBGAlarmLevels = {5, 10}
	alarm levels for injection background (in % of events)
std::vector< double >	m_timingAlarmLevels = {0.15, 0.30}
	alarm levels for time distribution (fraction of area difference)
std::vector< double >	m_eventT0MeanAlarmLevels = {8, 20}
	alarm levels for mean of event T0 [ns]
std::vector< double >	m_eventT0RmsAlarmLevels = {10, 20}
	alarm levels for r.m.s.
std::vector< double >	m_offsetMeanAlarmLevels = {0.2, 0.5}
	alarm levels for mean of bunch offset [ns]
std::vector< double >	m_offsetRmsAlarmLevels = {0.25, 0.50}
	alarm levels for r.m.s.
std::vector< int >	m_alarmColors = {c_ColorTooFew, c_ColorGood, c_ColorWarning, c_ColorError}
	alarm colors (see base class)
std::vector< int >	m_officialStates = {c_StatusTooFew, c_StatusGood, c_StatusWarning, c_StatusError}
	official alarm states
std::vector< bool >	m_includedBoardstacks
	boardstacks included in alarming
std::map< std::string, int >	m_bsmap
	a map of boardstack names to ID's
int	m_alarmStateOverall = 0
	overall alarm state of histograms to be sent by EpicsPV
double	m_averageRate = 0
	average BG rate (to pass to EpicsPV)
bool	m_IsNullRun = false
	Run type flag for null runs.
std::string	m_runType
	Run type.
double	m_numEvents = 0
	number of events processed with TOPDQM module
TH1D *	m_photonYields = nullptr
	photon yields per slot
TH1D *	m_backgroundRates = nullptr
	background rates per slot
TCanvas *	m_c_photonYields = nullptr
	Canvas: photon yields per slot.
TCanvas *	m_c_backgroundRates = nullptr
	Canvas: background rates per slot.
TH1F *	m_hotFraction = nullptr
	fraction of hot channels per slot (included boardstacks only)
TH1F *	m_deadFraction = nullptr
	fraction of dead channels per slot (included boardstacks only)
TH1F *	m_excludedFraction = nullptr
	fraction of dead and hot channels per slot in excluded boardstacks only
TH1F *	m_activeFraction = nullptr
	fraction of active channels per slot
THStack *	m_stack = nullptr
	stack for drawing dead, hot and active channel fractions
TLegend *	m_legend = nullptr
	legend for dead and hot channels
TCanvas *	m_c_deadAndHot = nullptr
	Canvas: fractin of dead and hot channels.
TH1F *	m_junkFraction = nullptr
	fraction of junk hits per boardstack
TH1F *	m_excludedBSHisto = nullptr
	histogram to show excluded boardstacks on junk fraction plot
TCanvas *	m_c_junkFraction = nullptr
	Canvas: fraction of junk hits per boardstack.
std::vector< TH1F * >	m_pmtHitRates
	histograms of PMT hits per event (index = slot - 1)
std::vector< TCanvas * >	m_c_pmtHitRates
	Canvases of PMT hits per event (index = slot - 1)
std::map< std::string, TCanvas * >	m_c_injBGs
	Canvases for projections of injection BG histograms.
std::map< std::string, TProfile * >	m_profiles
	profiles of injection BG
std::map< std::string, TH1D * >	m_projections
	projections of injection BG
TH1D *	m_skipProcFlagFract = nullptr
	fraction of events w/ skip processing flag set vs.
TH1D *	m_injVetoFlagFract = nullptr
	fraction of events w/ injection veto flag set vs.
TCanvas *	m_c_skipProcFlagFract = nullptr
	Canvas: fraction of events w/ skip processing flag set vs.
TCanvas *	m_c_injVetoFlagFract = nullptr
	Canvas: fraction of events w/ injection veto flag set vs.
std::vector< TLine * >	m_asicWindowsBandLines
	lines denoting a band of good windows
std::vector< TLine * >	m_verticalLines
	vertical lines splitting slots
std::vector< TLine * >	m_junkHitsAlarmLines
	lines representing alarm levels
std::vector< TLine * >	m_deadChannelsAlarmLines
	lines representing alarm levels
std::vector< TLine * >	m_backgroundAlarmLines
	lines representing alarm levels
std::vector< TLine * >	m_photonYieldsAlarmLines
	lines representing alarm levels
TLine *	m_injBGCutLine = nullptr
	a line denoting the cut on the number of hits for injection BG counting
TPaveText *	m_text1 = nullptr
	text to be written to window_vs_slot
TPaveText *	m_text2 = nullptr
	text to be written to event desynchonization monitor
TPaveText *	m_text3 = nullptr
	text to be written to background rates
TPaveText *	m_text4 = nullptr
	text to be written to number of good hits per event
std::map< std::string, double >	m_mirabelleVariables
	variables for MiraBelle
MonitoringObject *	m_monObj = nullptr
	MiraBelle 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

Class for TOP histogram analysis.

Definition at line 33 of file DQMHistAnalysisTOP.h.

Member Typedef Documentation

CanvasUpdatedList

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

inherited

The type of list of canvas updated status.

Definition at line 91 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 86 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 77 of file DQMHistAnalysis.h.

MonObjList

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

inherited

The type of list of MonitoringObjects.

Definition at line 81 of file DQMHistAnalysis.h.

RefList

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

inherited

The type of list of references.

Definition at line 96 of file DQMHistAnalysis.h.

Member Enumeration Documentation

EAlarmStates

enum EAlarmStates

private

Alarm states enumerator.

Definition at line 73 of file DQMHistAnalysisTOP.h.

                      {
      c_Gray = 0,    // undefined
      c_Green = 1,   // good
      c_Yellow = 2,  // warning
      c_Red = 3      // alarm
    };

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 68 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 46 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 57 of file DQMHistAnalysis.h.

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

Constructor & Destructor Documentation

DQMHistAnalysisTOPModule()

DQMHistAnalysisTOPModule

(

)

Constructor.

Definition at line 34 of file DQMHistAnalysisTOP.cc.

                                                  : DQMHistAnalysisModule()
{
  // Set description
  setDescription("Histogram analysis module for TOP DQM.");
 
  // Add parameters
  addParam("asicWindowsBand", m_asicWindowsBand,
           "lower and upper bin of a band denoting good windows", m_asicWindowsBand);
  addParam("asicWindowsAlarmLevels", m_asicWindowsAlarmLevels,
           "alarm levels for the fraction of windows outside the band (yellow, red)", m_asicWindowsAlarmLevels);
  addParam("eventMonitorAlarmLevels", m_eventMonitorAlarmLevels,
           "alarm levels for the fraction of desynchronized digits (yellow, red)", m_eventMonitorAlarmLevels);
  addParam("junkHitsAlarmLevels", m_junkHitsAlarmLevels,
           "alarm levels for the fraction of junk hits (yellow, red)", m_junkHitsAlarmLevels);
  addParam("deadChannelsAlarmLevels", m_deadChannelsAlarmLevels,
           "alarm levels for the fraction of dead + hot channels (yellow, red)", m_deadChannelsAlarmLevels);
  addParam("backgroundAlarmLevels", m_backgroundAlarmLevels,
           "alarm levels for background rates [MHz/PMT] (yellow, red)", m_backgroundAlarmLevels);
  addParam("photonYieldsAlarmLevels", m_photonYieldsAlarmLevels,
           "alarm levels for the number of photons per track (red, yellow)", m_photonYieldsAlarmLevels);
  addParam("excludedBoardstacks", m_excludedBoardstacks,
           "boarstacks to be excluded from alarming. Names are given like '5c', '13d' etc.", m_excludedBoardstacks);
  addParam("pvPrefix", m_pvPrefix, "Epics PV prefix", std::string("TOP:"));
  addParam("injectionBGAlarmLevels", m_injectionBGAlarmLevels,
           "alarm levels for injection background (in % of events)", m_injectionBGAlarmLevels);
  addParam("timingAlarmLevels", m_timingAlarmLevels,
           "alarm levels for time distribution (residual fraction w.r.t reference plot)", m_timingAlarmLevels);
  addParam("eventT0MeanAlarmLevels", m_eventT0MeanAlarmLevels,
           "alarm levels for mean of event T0 [ns]", m_eventT0MeanAlarmLevels);
  addParam("eventT0RmsAlarmLevels", m_eventT0RmsAlarmLevels,
           "alarm levels for r.m.s. of event T0 [ns]", m_eventT0RmsAlarmLevels);
  addParam("offsetMeanAlarmLevels", m_offsetMeanAlarmLevels,
           "alarm levels for mean of bunch offset [ns]", m_offsetMeanAlarmLevels);
  addParam("offsetRmsAlarmLevels", m_offsetRmsAlarmLevels,
           "alarm levels for r.m.s. of bunch offset [ns]", m_offsetRmsAlarmLevels);
 
  B2DEBUG(20, "DQMHistAnalysisTOP: Constructor done.");
}

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 95 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);
  n.setCanvas(nullptr);
}

beginRun()

void beginRun ( void )

finaloverridevirtual

Called when entering a new run.

Reimplemented from Module.

Definition at line 220 of file DQMHistAnalysisTOP.cc.

{
  m_mirabelleVariables.clear();
 
  B2DEBUG(20, "DQMHistAnalysisTOP: beginRun called.");
}

CheckEpicsError() [1/2]

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

inherited

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 807 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 799 of file DQMHistAnalysis.cc.

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

checkPVStatus()

void checkPVStatus ( void )

inherited

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

Definition at line 763 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 644 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 309 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 337 of file DQMHistAnalysis.cc.

{
  s_histList.clear();
}

clearlist()

void clearlist ( void )

inherited

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 342 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 750 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(); }

endRun()

void endRun ( void )

finaloverridevirtual

This method is called if the current run ends.

Reimplemented from Module.

Definition at line 311 of file DQMHistAnalysisTOP.cc.

{
  // add MiraBelle monitoring
 
  for (const auto& var : m_mirabelleVariables) {
    m_monObj->setVariable(var.first, var.second);
    B2DEBUG(20, var.first << " " << var.second);
  }
 
  B2DEBUG(20, "DQMHistAnalysisTOP : endRun called");
}

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 228 of file DQMHistAnalysisTOP.cc.

{
  // get type of the run (TODO: to be replaced with base class function when fixed)
  auto* rtype = findHist("DQMInfo/rtype");
  m_runType = rtype ? rtype->GetTitle() : "";
  m_IsNullRun = (m_runType == "null");
 
  // get number of events processed with TOPDQM module
  auto* goodHitsPerEvent = findHist("TOP/goodHitsPerEventAll");
  m_numEvents = goodHitsPerEvent ? goodHitsPerEvent->GetEntries() : 0;
 
  bool zeroSupp = gStyle->GetHistMinimumZero();
  gStyle->SetHistMinimumZero(true);
 
  // update alarm levels and other parameters from EpicsPVs
  updateLimits();
 
  // reset overall alarm state
  m_alarmStateOverall = c_Gray;
 
  // Update window_vs_slot canvas w/ alarming
  updateWindowVsSlotCanvas();
 
  // Update event desynchronization monitor w/ alarming
  updateEventMonitorCanvas();
 
  // Update number of good hits per event w/ alarming (injection BG)
  updateNGoodHitsCanvas();
 
  // Update event T0 w/ alarming
  updateEventT0Canvas();
 
  // Update bunch offset w/ alarming
  updateBunchOffsetCanvas();
 
  // Fraction of dead and hot channels
  const auto* activeFraction = makeDeadAndHotFractionsPlot();
 
  // Photon yields and background rates, corrected for dead and hot channels
  makePhotonYieldsAndBGRatePlots(activeFraction);
 
  // Fractions of junk hits
  makeJunkFractionPlot();
 
  // Set z-axis range to 3 times the average for good hits, 30 times the average for junk hits
  setZAxisRange("TOP/good_hits_xy_", 3);
  setZAxisRange("TOP/bad_hits_xy_", 30);
  setZAxisRange("TOP/good_hits_asics_", 3);
  setZAxisRange("TOP/bad_hits_asics_", 30);
 
  // Background subtracted time distributions (only for physics runs)
  if (m_runType == "physics") {
    const auto* trackHits = static_cast<TH2F*>(findHist("TOP/trackHits"));
    makeBGSubtractedTimingPlot("goodHitTimes", trackHits, 0);
    for (int slot = 1; slot <= 16; slot++) {
      makeBGSubtractedTimingPlot("good_timing_" + to_string(slot), trackHits, slot);
    }
  }
 
  // Update timing plot w/ alarming
  updateTimingCanvas();
 
  // PMT hit rates
  makePMTHitRatesPlots();
 
  // Injection BG
  makeInjectionBGPlots();
 
  // normalize histogram for injection veto flags check
  auto* injVetoFlagDiff = static_cast<TH1F*>(findHist("TOP/injVetoFlagDiff"));
  if (injVetoFlagDiff) injVetoFlagDiff->Scale(1 / injVetoFlagDiff->Integral(), "nosw2");
 
  // make flag fraction plots
  makeFlagFractPlot("TOP/skipProcFlag", m_skipProcFlagFract, m_c_skipProcFlagFract);
  makeFlagFractPlot("TOP/injVetoFlag", m_injVetoFlagFract, m_c_injVetoFlagFract);
 
  // Set Epics variables
  setEpicsVariables();
 
  gStyle->SetHistMinimumZero(zeroSupp);
}

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 376 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 364 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");
}

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 146 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() [1/2]

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 174 of file DQMHistAnalysis.cc.

{
  if (dirname.size() > 0) {
    return findHist(dirname + "/" + histname, updated);
  }
  return findHist(histname, updated);
}

findHist() [2/2]

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

staticinherited

Get histogram from list (no other search).

Parameters

histname	The name of the histogram (incl dir).
onlyIfUpdated	req only updated hists, return nullptr otherwise

Returns

The found histogram, or nullptr if not found.

Definition at line 160 of file DQMHistAnalysis.cc.

{
  if (s_histList.find(histname) != s_histList.end()) {
    if (was_updated && !s_histList[histname].isUpdated()) return nullptr;
    if (s_histList[histname].getHist()) {
      return s_histList[histname].getHist();
    } else {
      B2ERROR("Histogram " << histname << " in histogram list but nullptr.");
    }
  }
  B2INFO("Histogram " << histname << " 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 228 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 263 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 284 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() [1/2]

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 219 of file DQMHistAnalysis.cc.

{
  if (dirname.size() > 0) {
    return findRefHist(dirname + "/" + histname, scaling, hist);
  }
  return findRefHist(histname, scaling, hist);
}

findRefHist() [2/2]

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

staticinherited

Get referencehistogram from list (no other search).

Parameters

histname	The name of the histogram (incl dir).
scaling	enum what scaling to use
hist	histogram to scale to

Returns

The found histogram, or nullptr if not found.

Definition at line 208 of file DQMHistAnalysis.cc.

{
  if (s_refList.find(histname) != 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[histname].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;
}

getAlarmColor()

int getAlarmColor ( unsigned alarmState ) const

inlineprivate

Converts alarm state to color.

Parameters

alarmState

alarm state

Returns

alarm color

Definition at line 189 of file DQMHistAnalysisTOP.h.

    {
      if (alarmState < m_alarmColors.size()) return m_alarmColors[alarmState];
      return m_alarmColors[c_Gray];
    }

getAlarmState()

int getAlarmState ( double value, const std::vector< double > & alarmLevels, bool bigRed = true ) const

private

Returns alarm state.

Parameters

value	value to translate to color
alarmLevels	alarm levels
bigRed	on true red color for large values, else red color for small values

Returns

alarm state (to be used as an index of m_alarmColors array)

Definition at line 948 of file DQMHistAnalysisTOP.cc.

{
  if (m_IsNullRun or m_numEvents < 1000) return c_Gray;
 
  if (bigRed) {
    if (value < alarmLevels[0]) return c_Green;
    else if (value < alarmLevels[1]) return c_Yellow;
    else return c_Red;
  } else {
    if (value < alarmLevels[0]) return c_Red;
    else if (value < alarmLevels[1]) return c_Yellow;
    else return c_Green;
  }
}

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 202 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;
}

getDeadAndHotCuts()

std::pair< double, double > getDeadAndHotCuts ( const TH1 * h )

staticprivate

Returns cut levels for dead and hot channels.

Parameters

h	pixel or channel distribution of hits (1D or 2D histogram)

Returns

cut levels (first = dead, second = hot)

Definition at line 1010 of file DQMHistAnalysisTOP.cc.

{
  std::vector<double> binContents;
  for (int k = 1; k <= h->GetNbinsY(); k++) {
    for (int i = 1; i <= h->GetNbinsX(); i++) {
      binContents.push_back(h->GetBinContent(i, k));
    }
  }
 
  double mean = 0;
  double rms = h->GetMaximum();
  for (int iter = 0; iter < 5; iter++) {
    double sumy = 0;
    double sumyy = 0;
    int n = 0;
    for (auto y : binContents) {
      if (y == 0 or fabs(y - mean) > 3 * rms) continue;
      sumy += y;
      sumyy += y * y;
      n++;
    }
    if (n == 0) continue;
    mean = sumy / n;
    rms = sqrt(sumyy / n - mean * mean);
  }
 
  return std::make_pair(mean / 5, std::max(mean * 2, mean + 6 * rms));
}

getDelta() [1/2]

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 118 of file DQMHistAnalysis.cc.

{
  std::string fullname;
  if (dirname.size() > 0) {
    fullname = dirname + "/" + histname;
  } else {
    fullname = histname;
  }
  return getDelta(fullname, n, onlyIfUpdated);
}

getDelta() [2/2]

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

inherited

Get Delta histogram.

Parameters

fullname	directory+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 129 of file DQMHistAnalysis.cc.

{
  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 196 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 507 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], (void*)&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 530 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], (void*)&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 603 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 617 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 553 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 578 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 221 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 184 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 139 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 190 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;}

getOffcialAlarmStatus()

int getOffcialAlarmStatus ( unsigned alarmState ) const

inlineprivate

Converts alarm state to official status (see EStatus of the base class)

Parameters

alarmState

alarm state

Returns

alarm status

Definition at line 200 of file DQMHistAnalysisTOP.h.

    {
      if (alarmState < m_officialStates.size()) return m_officialStates[alarmState];
      return c_StatusDefault;
    }

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

getPVPrefix()

static std::string & getPVPrefix ( void )

inlinestaticinherited

get global prefix for EPICS PVs

Returns

prefix in use

Definition at line 610 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 215 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 293 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 723 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 593 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 599 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 107 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 323 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 74 of file DQMHistAnalysisTOP.cc.

{
 
  // check module parameters
 
  if (m_asicWindowsBand.size() != 2) B2ERROR("Parameter list 'asicWindowsBand' must contain two numbers");
  if (m_asicWindowsAlarmLevels.size() != 2) B2ERROR("Parameter list 'asicWindowsAlarmLevels' must contain two numbers");
  if (m_eventMonitorAlarmLevels.size() != 2) B2ERROR("Parameter list 'eventMonitorAlarmLevels' must contain two numbers");
  if (m_junkHitsAlarmLevels.size() != 2) B2ERROR("Parameter list 'junkHitsAlarmLevels' must contain two numbers");
  if (m_deadChannelsAlarmLevels.size() != 2) B2ERROR("Parameter list 'deadChannelsAlarmLevels' must contain two numbers");
  if (m_backgroundAlarmLevels.size() != 2) B2ERROR("Parameter list 'backgroundAlarmLevels' must contain two numbers");
  if (m_photonYieldsAlarmLevels.size() != 2) B2ERROR("Parameter list 'photonYieldsAlarmLevels' must contain two numbers");
  if (m_injectionBGAlarmLevels.size() != 2) B2ERROR("Parameter list 'injectionBGAlarmLevels' must contain two numbers");
  if (m_timingAlarmLevels.size() != 2) B2ERROR("Parameter list 'timingAlarmLevels' must contain two numbers");
  if (m_eventT0MeanAlarmLevels.size() != 2) B2ERROR("Parameter list 'eventT0MeanAlarmLevels' must contain two numbers");
  if (m_eventT0RmsAlarmLevels.size() != 2) B2ERROR("Parameter list 'eventT0RmsAlarmLevels' must contain two numbers");
  if (m_offsetMeanAlarmLevels.size() != 2) B2ERROR("Parameter list 'offsetMeanAlarmLevels' must contain two numbers");
  if (m_offsetRmsAlarmLevels.size() != 2) B2ERROR("Parameter list 'offsetRmsAlarmLevels' must contain two numbers");
 
  // make a map of boardstack names to ID's
 
  int id = 1;
  for (int slot = 1; slot <= 16; slot++) {
    string slotstr = to_string(slot);
    for (std::string bs : {"a", "b", "c", "d"}) {
      m_bsmap[slotstr + bs] = id;
      id++;
    }
  }
 
  // parse excluded boardstacks
 
  setIncludedBoardstacks(m_excludedBoardstacks);
 
  // MiraBelle monitoring
 
  m_monObj = getMonitoringObject("top");
 
  // Epics used to pass values to shifter's page (output only)
 
  registerEpicsPV(m_pvPrefix + "badBoardstacks", "badBoardstacks");
  registerEpicsPV(m_pvPrefix + "badCarriers", "badCarriers");
  registerEpicsPV(m_pvPrefix + "badAsics", "badAsics");
  registerEpicsPV(m_pvPrefix + "badPMTs", "badPMTs");
  registerEpicsPV(m_pvPrefix + "numExcludedBS", "numExcludedBS");
  registerEpicsPV(m_pvPrefix + "histoAlarmState", "histoAlarmState"); // to pass overall state to central alarm overview panel
 
  // Epics used to get limits from configuration file - override module parameters (input only)
 
  registerEpicsPV(m_pvPrefix + "asicWindowsBand", "asicWindowsBand");
  registerEpicsPV(m_pvPrefix + "asicWindowsAlarmLevels", "asicWindowsAlarmLevels");
  registerEpicsPV(m_pvPrefix + "eventMonitorAlarmLevels", "eventMonitorAlarmLevels");
  registerEpicsPV(m_pvPrefix + "junkHitsAlarmLevels", "junkHitsAlarmLevels");
  registerEpicsPV(m_pvPrefix + "deadChannelsAlarmLevels", "deadChannelsAlarmLevels");
  registerEpicsPV(m_pvPrefix + "backgroundAlarmLevels", "backgroundAlarmLevels"); // also output
  registerEpicsPV(m_pvPrefix + "photonYieldsAlarmLevels", "photonYieldsAlarmLevels");
  registerEpicsPV(m_pvPrefix + "excludedBoardstacks", "excludedBoardstacks");
 
  registerEpicsPV(m_pvPrefix + "injectionBGAlarmLevels", "injectionBGAlarmLevels"); // also output
  registerEpicsPV(m_pvPrefix + "timingAlarmLevels", "timingAlarmLevels");
  registerEpicsPV(m_pvPrefix + "eventT0MeanAlarmLevels", "eventT0MeanAlarmLevels");
  registerEpicsPV(m_pvPrefix + "eventT0RmsAlarmLevels", "eventT0RmsAlarmLevels");
  registerEpicsPV(m_pvPrefix + "offsetMeanAlarmLevels", "offsetMeanAlarmLevels");
  registerEpicsPV(m_pvPrefix + "offsetRmsAlarmLevels", "offsetRmsAlarmLevels");
 
  // new canvases, histograms and graphic primitives
 
  gROOT->cd();
 
  m_c_photonYields = new TCanvas("TOP/c_photonYields", "c_photonYields");
  m_c_backgroundRates = new TCanvas("TOP/c_backgroundRates", "c_backgroundRates");
 
  m_deadFraction = new TH1F("TOP/deadFraction", "Fraction of dead channels in included boardstacks", 16, 0.5, 16.5);
  m_deadFraction->SetXTitle("slot number");
  m_deadFraction->SetYTitle("fraction");
  m_hotFraction = new TH1F("TOP/hotFraction", "Fraction of hot channels in included boardstacks", 16, 0.5, 16.5);
  m_hotFraction->SetXTitle("slot number");
  m_hotFraction->SetYTitle("fraction");
  m_excludedFraction = new TH1F("TOP/excludedFraction", "Fraction of hot and dead channels in excluded bordstacks", 16, 0.5, 16.5);
  m_excludedFraction->SetXTitle("slot number");
  m_excludedFraction->SetYTitle("fraction");
  m_activeFraction = new TH1F("TOP/activeFraction", "Fraction of active channels", 16, 0.5, 16.5);
  m_activeFraction->SetXTitle("slot number");
  m_activeFraction->SetYTitle("fraction");
  m_c_deadAndHot = new TCanvas("TOP/c_deadAndHotChannels", "c_deadAndHotChannels");
 
  m_junkFraction = new TH1F("TOP/junkFraction", "Fraction of junk hits per boardstack", 64, 0.5, 16.5);
  m_junkFraction->SetXTitle("slot number");
  m_junkFraction->SetYTitle("fraction");
  // note: titles are intentionally the same since this one is plotted first
  m_excludedBSHisto = new TH1F("TOP/excludedBSHisto", "Fraction of junk hits per boardstack", 64, 0.5, 16.5);
  m_excludedBSHisto->SetXTitle("slot number");
  m_excludedBSHisto->SetYTitle("fraction");
  m_c_junkFraction = new TCanvas("TOP/c_junkFraction", "c_junkFraction");
 
  for (int slot = 1; slot <= 16; slot++) {
    string hname = "TOP/pmtHitRates_" + to_string(slot);
    string htitle = "PMT hits per event for slot #" + to_string(slot);
    auto* h = new TH1F(hname.c_str(), htitle.c_str(), 32, 0.5, 32.5);
    h->SetXTitle("PMT number");
    h->SetYTitle("Number of good hits per event");
    m_pmtHitRates.push_back(h);
    string cname = "TOP/c_pmtHitRates_" + to_string(slot);
    string ctitle = "c_pmtHitRates_" + to_string(slot);
    m_c_pmtHitRates.push_back(new TCanvas(cname.c_str(), ctitle.c_str()));
  }
 
  for (std::string name : {
         "nhitInjLER", "nhitInjHER", "nhitInjLERcut", "nhitInjHERcut",
         "eventInjLER", "eventInjHER", "eventInjLERcut", "eventInjHERcut"
       }) {
    for (std::string proj : {"_px", "_py"}) {
      std::string cname = "TOP/c_" + name + proj;
      m_c_injBGs[cname] = new TCanvas(cname.c_str(), (name + proj).c_str());
    }
  }
 
  m_c_skipProcFlagFract = new TCanvas("TOP/c_skipProcFlagFract", "c_skipProcFlagFract");
  m_c_injVetoFlagFract = new TCanvas("TOP/c_injVetoFlagFract", "c_injVetoFlagFract");
 
  m_text1 = new TPaveText(0.125, 0.8, 0.675, 0.88, "NDC");
  m_text1->SetFillColorAlpha(kWhite, 0);
  m_text1->SetBorderSize(0);
  m_text2 = new TPaveText(0.55, 0.8, 0.85, 0.89, "NDC");
  m_text2->SetFillColorAlpha(kWhite, 0);
  m_text2->SetBorderSize(0);
  m_text3 = new TPaveText(0.47, 0.8, 0.85, 0.89, "NDC");
  m_text3->SetFillColorAlpha(kWhite, 0);
  m_text3->SetBorderSize(0);
  m_text4 = new TPaveText(0.125, 0.8, 0.675, 0.88, "NDC");
  m_text4->SetFillColorAlpha(kWhite, 0);
  m_text4->SetBorderSize(0);
 
  for (int slot = 1; slot < 16; slot++) {
    auto* line = new TLine(slot + 0.5, 0, slot + 0.5, 1);
    line->SetLineWidth(1);
    line->SetLineStyle(2);
    m_verticalLines.push_back(line);
  }
 
  setAlarmLines();
 
  B2DEBUG(20, "DQMHistAnalysisTOP: initialized.");
}

makeBGSubtractedTimingPlot()

void makeBGSubtractedTimingPlot ( const std::string & name, const TH2F * trackHits, int slot )

staticprivate

Makes background subtracted time distribution plot.

Parameters

name	the name of the histogram
trackHits	histogram used to scale background in case it is available
slot	slot number

Definition at line 796 of file DQMHistAnalysisTOP.cc.

{
  auto* canvas = findCanvas("TOP/c_" + name);
  if (not canvas) return;
 
  auto* h = static_cast<TH1F*>(findHist("TOP/" + name));
  if (not h) return;
 
  auto* hb = static_cast<TH1F*>(findHist("TOP/" + name + "BG"));
  if (not hb) return;
 
  if (trackHits) {
    // use the ratio of events w/ and w/o track in the slot to scale the background
    double s = (slot == 0) ? trackHits->Integral(1, 16, 2, 2) : trackHits->GetBinContent(slot, 2);
    if (s == 0) return;
    double sb = (slot == 0) ? trackHits->Integral(1, 16, 1, 1) : trackHits->GetBinContent(slot, 1);
    if (sb == 0) return;
    h->Add(h, hb, 1, -s / sb);
  } else {
    // use the content of bins at t < 0 to scale the background
    int i0 = h->GetXaxis()->FindBin(0.); // bin at t = 0
    double s = h->Integral(1, i0);
    if (s == 0) return;
    double sb = hb->Integral(1, i0);
    if (sb == 0) return;
    if (s / sb > 1) return; // this can happen due to low statistics and is not reliable
    h->Add(h, hb, 1, -s / sb);
  }
 
  TString title = TString(h->GetTitle()) + " (BG subtracted)";
  h->SetTitle(title);
 
  canvas->Clear();
  canvas->cd();
  h->Draw();
  canvas->Modified();
}

makeDeadAndHotFractionsPlot()

const TH1F * makeDeadAndHotFractionsPlot ( )

private

Makes a plot of dead and hot channel fractions per slot.

Returns

fraction of active channels per slot

Definition at line 543 of file DQMHistAnalysisTOP.cc.

{
  m_deadFraction->Reset();
  m_hotFraction->Reset();
  m_excludedFraction->Reset();
  m_activeFraction->Reset();
  double inactiveFract = 0; // max inactive channel fraction when some boardstacks are excluded from alarming
 
  for (int slot = 1; slot <= 16; slot++) {
    auto* h = static_cast<TH1F*>(findHist("TOP/good_channel_hits_" + std::to_string(slot)));
    if (not h) continue;
 
    auto cuts = getDeadAndHotCuts(h);
    double deadCut = cuts.first;
    double hotCut = cuts.second;
    double deadFract = 0;
    double hotFract = 0;
    double deadFractIncl = 0;
    double hotFractIncl = 0;
    for (int chan = 0; chan < h->GetNbinsX(); chan++) {
      double y = h->GetBinContent(chan + 1);
      int bs = chan / 128 + (slot - 1) * 4;
      bool included = m_includedBoardstacks[bs];
      if (y <= deadCut) {
        deadFract += 1;
        if (included) deadFractIncl += 1;
      } else if (y > hotCut) {
        hotFract += 1;
        if (included) hotFractIncl += 1;
      }
    }
    deadFract /= h->GetNbinsX();
    hotFract /= h->GetNbinsX();
    deadFractIncl /= h->GetNbinsX();
    hotFractIncl /= h->GetNbinsX();
    m_deadFraction->SetBinContent(slot, deadFractIncl);
    m_hotFraction->SetBinContent(slot, hotFractIncl);
    m_excludedFraction->SetBinContent(slot, deadFract - deadFractIncl + hotFract - hotFractIncl);
    m_activeFraction->SetBinContent(slot, 1 - deadFract - hotFract);
    inactiveFract = std::max(inactiveFract, deadFractIncl + hotFractIncl);
  }
 
  setMiraBelleVariables("ActiveChannelFraction_slot", m_activeFraction);
 
  int alarmState = c_Gray;
  if (m_activeFraction->Integral() > 0) {
    alarmState = getAlarmState(inactiveFract, m_deadChannelsAlarmLevels);
  }
 
  m_alarmStateOverall = std::max(m_alarmStateOverall, alarmState);
 
  m_deadFraction->SetFillColor(1);
  m_deadFraction->SetLineColor(1);
  m_deadFraction->GetXaxis()->SetNdivisions(16);
 
  m_hotFraction->SetFillColor(2);
  m_hotFraction->SetLineColor(2);
  m_hotFraction->GetXaxis()->SetNdivisions(16);
 
  m_excludedFraction->SetFillColor(kGray);
  m_excludedFraction->SetLineColor(kGray);
  m_excludedFraction->GetXaxis()->SetNdivisions(16);
 
  m_activeFraction->SetFillColor(0);
  m_activeFraction->GetXaxis()->SetNdivisions(16);
 
  auto* canvas = m_c_deadAndHot;
  canvas->Clear();
  canvas->cd();
  canvas->Pad()->SetFrameFillColor(10);
  if (not m_stack) {
    m_stack = new THStack("TOP/stack", "Fraction of dead and hot channels");
    m_stack->Add(m_deadFraction);
    m_stack->Add(m_hotFraction);
    m_stack->Add(m_excludedFraction);
    m_stack->Add(m_activeFraction);
  }
  m_stack->Draw();
 
  for (auto* line : m_deadChannelsAlarmLines) line->Draw("same");
 
  if (not m_legend) {
    m_legend = new TLegend(0.8, 0.87, 0.99, 0.99);
    m_legend->AddEntry(m_hotFraction, "hot");
    m_legend->AddEntry(m_deadFraction, "dead");
    m_legend->AddEntry(m_excludedFraction, "excluded");
  }
  m_legend->Draw("same");
 
  canvas->Pad()->SetFillColor(getAlarmColor(alarmState));
  canvas->Modified();
 
  return m_activeFraction;
}

makeFlagFractPlot()

void makeFlagFractPlot ( const std::string & hname, TH1 * histogram, TCanvas * canvas )

staticprivate

Makes a plot of fraction of events with the flag is set.

Parameters

hname	name of 2D histogram (y-axis = flag value)
histogram	resulting histogram w/ fractions
canvas	canvas to plot

Definition at line 916 of file DQMHistAnalysisTOP.cc.

{
  if (histogram) delete histogram;
  auto* h = static_cast<TH2F*>(findHist(hname));
  if (not h) return;
 
  histogram = h->ProjectionX((hname + "Fract").c_str(), 2, 2);
  auto* px = h->ProjectionX("tmp");
  histogram->Divide(histogram, px, 1, 1, "B");
  delete px;
  histogram->SetTitle(TString(h->GetTitle()) + " is set");
  histogram->SetYTitle("fraction of events");
  histogram->SetMarkerStyle(24);
  histogram->SetMinimum(0);
 
  if (not canvas) return;
  canvas->Clear();
  canvas->cd();
  histogram->Draw();
  canvas->Modified();
}

makeInjectionBGPlots()

void makeInjectionBGPlots ( )

private

Makes projections of injection BG plots.

Definition at line 868 of file DQMHistAnalysisTOP.cc.

{
  for (std::string name : {"nhitInjLER", "nhitInjHER", "nhitInjLERcut", "nhitInjHERcut"}) {
    std::string hname = "TOP/" + name;
    auto* h = static_cast<TProfile2D*>(findHist(hname));
    if (not h) continue;
    for (std::string proj : {"_px", "_py"}) {
      std::string cname = "TOP/c_" + name + proj;
      auto* canvas = m_c_injBGs[cname];
      if (not canvas) continue;
      canvas->Clear();
      canvas->cd();
      auto& hproj = m_profiles[cname];
      if (hproj) delete hproj;
      hproj = (proj == "_px") ? h->ProfileX((hname + proj).c_str()) : h->ProfileY((hname + proj).c_str());
      std::string xtitle = (proj == "_px") ? h->GetXaxis()->GetTitle() : h->GetYaxis()->GetTitle();
      hproj->SetXTitle(xtitle.c_str());
      hproj->SetYTitle(h->GetZaxis()->GetTitle());
      hproj->SetMinimum(0);
      hproj->Draw("hist");
      canvas->Modified();
    }
  }
 
  for (std::string name : {"eventInjLER", "eventInjHER", "eventInjLERcut", "eventInjHERcut"}) {
    std::string hname = "TOP/" + name;
    auto* h = static_cast<TH2F*>(findHist(hname));
    if (not h) continue;
    for (std::string proj : {"_px", "_py"}) {
      std::string cname = "TOP/c_" + name + proj;
      auto* canvas = m_c_injBGs[cname];
      if (not canvas) continue;
      canvas->Clear();
      canvas->cd();
      auto& hproj = m_projections[cname];
      if (hproj) delete hproj;
      hproj = (proj == "_px") ? h->ProjectionX((hname + proj).c_str()) : h->ProjectionY((hname + proj).c_str());
      std::string xtitle = (proj == "_px") ? h->GetXaxis()->GetTitle() : h->GetYaxis()->GetTitle();
      hproj->SetXTitle(xtitle.c_str());
      hproj->SetYTitle(h->GetZaxis()->GetTitle());
      hproj->SetMinimum(0);
      hproj->Draw("hist");
      canvas->Modified();
    }
  }
 
}

makeJunkFractionPlot()

void makeJunkFractionPlot ( )

private

Makes a plot of fractions of junk hits per boardstack.

Definition at line 728 of file DQMHistAnalysisTOP.cc.

{
  m_junkFraction->Reset();
  m_excludedBSHisto->Reset();
  auto* allHits = static_cast<TH1D*>(m_junkFraction->Clone("tmp"));
  for (int slot = 1; slot <= 16; slot++) {
    auto* good = static_cast<TH1F*>(findHist("TOP/good_channel_hits_" + std::to_string(slot)));
    if (not good) continue;
    auto* bad = static_cast<TH1F*>(findHist("TOP/bad_channel_hits_" + std::to_string(slot)));
    if (not bad) continue;
    for (int i = 0; i < 512; i++) {
      int bs = i / 128;
      allHits->Fill(slot + bs / 4. - 0.5, good->GetBinContent(i + 1) + bad->GetBinContent(i + 1));
      m_junkFraction->Fill(slot + bs / 4. - 0.5, bad->GetBinContent(i + 1));
    }
  }
 
  m_junkFraction->Divide(m_junkFraction, allHits, 1, 1, "B");
 
  int alarmState = c_Gray;
  if (allHits->Integral() > 0) {
    double hmax = 0;
    for (size_t i = 0; i < m_includedBoardstacks.size(); i++) {
      if (m_includedBoardstacks[i]) hmax = std::max(hmax, m_junkFraction->GetBinContent(i + 1));
      else m_excludedBSHisto->SetBinContent(i + 1, 1);
    }
    alarmState = getAlarmState(hmax, m_junkHitsAlarmLevels);
  }
  delete allHits;
  m_alarmStateOverall = std::max(m_alarmStateOverall, alarmState);
 
  auto* canvas = m_c_junkFraction;
  canvas->Clear();
  canvas->cd();
  canvas->Pad()->SetFrameFillColor(10);
  canvas->Pad()->SetFillColor(getAlarmColor(alarmState));
  m_excludedBSHisto->SetFillColor(kGray);
  m_excludedBSHisto->SetLineColor(kGray);
  m_excludedBSHisto->GetXaxis()->SetNdivisions(16);
  m_excludedBSHisto->GetYaxis()->SetRangeUser(0, 1);
  m_excludedBSHisto->Draw();
  m_junkFraction->SetMarkerStyle(24);
  m_junkFraction->GetXaxis()->SetNdivisions(16);
  m_junkFraction->GetYaxis()->SetRangeUser(0, 1); // Note: m_junkFraction->GetMaximum() will now give 1 and not the histogram maximum!
  m_junkFraction->Draw("same");
  for (auto* line : m_verticalLines) line->Draw("same");
  for (auto* line : m_junkHitsAlarmLines) line->Draw("same");
  canvas->Modified();
}

makePhotonYieldsAndBGRatePlots()

void makePhotonYieldsAndBGRatePlots ( const TH1F * activeFraction )

private

Make plots of dead-and-hot-channel corrected photon yields and BG rates per slot.

Parameters

activeFraction

fraction of active channels per slot

Definition at line 639 of file DQMHistAnalysisTOP.cc.

{
  for (auto* canvas : {m_c_photonYields, m_c_backgroundRates}) {
    canvas->Clear();
    canvas->Pad()->SetFrameFillColor(10);
    canvas->Pad()->SetFillColor(getAlarmColor(c_Gray));
    canvas->Modified();
  }
  m_averageRate = 0;
 
  auto* signalHits = static_cast<TProfile*>(findHist("TOP/signalHits"));
  if (not signalHits) return;
 
  auto* backgroundHits = static_cast<TProfile*>(findHist("TOP/backgroundHits"));
  if (not backgroundHits) return;
 
  if (m_photonYields) delete m_photonYields;
  m_photonYields = signalHits->ProjectionX("TOP/photonYields");
  if (m_backgroundRates) delete m_backgroundRates;
  m_backgroundRates = backgroundHits->ProjectionX("TOP/backgroundRates");
  auto* activeFract = static_cast<TH1F*>(activeFraction->Clone("tmp"));
  for (int i = 1; i <= activeFract->GetNbinsX(); i++) activeFract->SetBinError(i, 0);
 
  m_photonYields->Add(m_photonYields, m_backgroundRates, 1, -1);
  m_photonYields->Divide(m_photonYields, activeFract);
  setMiraBelleVariables("PhotonsPerTrack_slot", m_photonYields);
 
  int alarmState = c_Gray;
  if (signalHits->GetEntries() > 0 and activeFraction->Integral() > 0) {
    double hmin = 1000;
    for (int i = 1; i <= m_photonYields->GetNbinsX(); i++) {
      if (signalHits->GetBinEntries(i) < 10) continue;
      hmin = std::min(hmin, m_photonYields->GetBinContent(i) + 3 * m_photonYields->GetBinError(i));
    }
    if (hmin < 1000) alarmState = getAlarmState(hmin, m_photonYieldsAlarmLevels, false);
  }
  m_alarmStateOverall = std::max(m_alarmStateOverall, alarmState);
 
  m_photonYields->SetTitle("Number of photons per track");
  m_photonYields->SetYTitle("photons per track");
  m_photonYields->SetMarkerStyle(24);
  m_photonYields->GetXaxis()->SetNdivisions(16);
 
  auto* canvas = m_c_photonYields;
  canvas->cd();
  m_photonYields->SetMinimum(0);
  m_photonYields->Draw();
  for (auto* line : m_photonYieldsAlarmLines) line->Draw("same");
  canvas->Pad()->SetFillColor(getAlarmColor(alarmState));
  canvas->Modified();
 
  m_backgroundRates->Scale(1.0 / 50.0e-3 / 32);  // measured in 50 ns window, 32 PMT's ==> rate in MHz/PMT
  m_backgroundRates->Divide(m_backgroundRates, activeFract);
  setMiraBelleVariables("BackgroundRate_slot", m_backgroundRates);
 
  alarmState = c_Gray;
  m_text3->Clear();
  if (backgroundHits->GetEntries() > 100 and activeFraction->Integral() > 0) {
    int status = m_backgroundRates->Fit("pol0", "Q0");
    if (status == 0) {
      auto* fun = m_backgroundRates->GetFunction("pol0");
      if (fun) {
        m_averageRate = fun->GetParameter(0);
        double error = fun->GetParError(0);
        alarmState = getAlarmState(m_averageRate - 3 * error, m_backgroundAlarmLevels);
        m_text3->AddText(Form("Average: %.2f MHz/PMT", m_averageRate));
      }
    }
  }
  m_alarmStateOverall = std::max(m_alarmStateOverall, alarmState);
 
  m_backgroundRates->SetTitle("Background rates");
  m_backgroundRates->SetYTitle("background rate [MHz/PMT]");
  m_backgroundRates->SetMarkerStyle(24);
  m_backgroundRates->GetXaxis()->SetNdivisions(16);
 
  canvas = m_c_backgroundRates;
  canvas->cd();
  m_backgroundRates->SetMinimum(0);
  m_backgroundRates->Draw();
  for (auto* line : m_backgroundAlarmLines) line->Draw("same");
  m_text3->Draw();
  canvas->Pad()->SetFillColor(getAlarmColor(alarmState));
  canvas->Modified();
 
  delete activeFract;
}

makePMTHitRatesPlots()

void makePMTHitRatesPlots ( )

private

Makes plots of the number of PMT hits per event.

Definition at line 835 of file DQMHistAnalysisTOP.cc.

{
  auto* h0 = static_cast<TH1F*>(findHist("TOP/goodHitsPerEventAll"));
  if (not h0) return;
  double numEvents = h0->GetEntries();
  if (numEvents == 0) return;
 
  int numSlots = m_pmtHitRates.size();
  for (int slot = 1; slot <= numSlots; slot++) {
    string name = "TOP/good_hits_xy_" + to_string(slot);
    auto* hxy = static_cast<TH2F*>(findHist(name));
    if (not hxy) continue;
    std::vector<double> pmts(32, 0);
    for (int row = 0; row < 8; row++) {
      for (int col = 0; col < 64; col++) {
        int pmt = col / 4 + (row / 4) * 16;
        pmts[pmt] += hxy->GetBinContent(col + 1, row + 1);
      }
    }
    auto* h = m_pmtHitRates[slot - 1];
    for (size_t i = 0; i < pmts.size(); i++) {
      h->SetBinContent(i + 1, pmts[i] / numEvents);
    }
    auto* canvas = m_c_pmtHitRates[slot - 1];
    canvas->Clear();
    canvas->cd();
    h->SetMinimum(0);
    h->Draw();
    canvas->Modified();
  }
}

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 705 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 773 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;
  }
}

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 388 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 398 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 393 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 670 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 658 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 664 of file DQMHistAnalysis.cc.

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

resetDeltaList()

void resetDeltaList ( void )

staticinherited

Reset Delta.

Definition at line 347 of file DQMHistAnalysis.cc.

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

sameHistDefinition()

bool sameHistDefinition ( TH1 * h1, TH1 * h2 )

staticprivate

Checks if histograms are defined in the same way (nbins, xmin, xmax)

Parameters

h1	first histogram
h2	second histogram

Returns

true if (nbins, xmin, xmax) match

Definition at line 535 of file DQMHistAnalysisTOP.cc.

{
  if (h1->GetNbinsX() != h2->GetNbinsX()) return false;
  if (h1->GetXaxis()->GetXmin() != h2->GetXaxis()->GetXmin()) return false;
  if (h1->GetXaxis()->GetXmax() != h2->GetXaxis()->GetXmax()) return false;
  return true;
}

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 182 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));
}

setAlarmLines() [1/2]

void setAlarmLines ( )

private

Sets all alarm lines.

Definition at line 987 of file DQMHistAnalysisTOP.cc.

{
  for (size_t i = 0; i < m_asicWindowsBand.size(); i++) {
    double y = m_asicWindowsBand[i];
    if (i < m_asicWindowsBandLines.size()) {
      auto* line = m_asicWindowsBandLines[i];
      line->SetY1(y);
      line->SetY2(y);
    } else {
      auto* line = new TLine(0.5, y, 16.5, y);
      line->SetLineWidth(2);
      line->SetLineColor(kRed);
      m_asicWindowsBandLines.push_back(line);
    }
  }
 
  setAlarmLines(m_junkHitsAlarmLevels, 0.5, 16.5, m_junkHitsAlarmLines);
  setAlarmLines(m_deadChannelsAlarmLevels, 0.5, 16.5, m_deadChannelsAlarmLines);
  setAlarmLines(m_backgroundAlarmLevels, 0.5, 16.5, m_backgroundAlarmLines);
  setAlarmLines(m_photonYieldsAlarmLevels, 0.5, 16.5, m_photonYieldsAlarmLines, false);
}

setAlarmLines() [2/2]

void setAlarmLines ( const std::vector< double > & alarmLevels, double xmin, double xmax, std::vector< TLine * > & alarmLines, bool bigRed = true )

staticprivate

Sets alarm lines.

Parameters

alarmLevels	alarm levels
xmin	minimal x
xmax	maximal x
alarmLines	alarm lines [out]
bigRed	on true red color for large values, else red color for small values

Definition at line 964 of file DQMHistAnalysisTOP.cc.

{
  std::vector<int> colors = {kOrange, kRed};
  if (not bigRed) std::reverse(colors.begin(), colors.end());
  for (size_t i = 0; i < std::min(colors.size(), alarmLevels.size()); i++) {
    if (i < alarmLines.size()) {
      auto* line = alarmLines[i];
      line->SetX1(xmin);
      line->SetX2(xmax);
      line->SetY1(alarmLevels[i]);
      line->SetY2(alarmLevels[i]);
    } else {
      auto* line = new TLine(xmin, alarmLevels[i], xmax, alarmLevels[i]);
      line->SetLineWidth(2);
      line->SetLineStyle(2);
      line->SetLineColor(colors[i]);
      alarmLines.push_back(line);
    }
  }
}

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 426 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], (void*)&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 438 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], (void*)&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 468 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], (void*)&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 480 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], (void*)&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 450 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, std::string value )

inherited

Write string to a EPICS PV.

Parameters

index	index of PV
value	string to write

Definition at line 492 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
}

setEpicsVariables()

void setEpicsVariables ( )

private

Calculates and sets epics variables.

Definition at line 1040 of file DQMHistAnalysisTOP.cc.

{
  int badBoardstacks = 0;
  int badCarriers = 0;
  int badAsics = 0;
  for (int slot = 1; slot <= 16; slot++) {
    std::string hname = "TOP/good_hits_asics_" + to_string(slot);
    auto* h = static_cast<TH2F*>(findHist(hname));
    if (not h) continue;
 
    auto cuts = getDeadAndHotCuts(h);
    double deadCut = cuts.first;
    double hotCut = cuts.second;
    std::vector<int> asics(64, 0);
    std::vector<int> carriers(16, 0);
    std::vector<int> boardstacks(4, 0);
    for (int asic = 0; asic < 64; asic++) {
      int carrier = asic / 4;
      int boardstack = carrier / 4;
      for (int chan = 0; chan < 8; chan++) {
        double y = h->GetBinContent(asic + 1, chan + 1);
        if (y > deadCut and y <= hotCut) {
          asics[asic]++;
          carriers[carrier]++;
          boardstacks[boardstack]++;
        }
      }
    }
    for (int n : asics) if (n == 0) badAsics++;
    for (int n : carriers) if (n == 0) badCarriers++;
    for (int n : boardstacks) if (n == 0) badBoardstacks++;
  }
  badAsics -= badCarriers * 4;
  badCarriers -= badBoardstacks * 4;
 
  int badPMTs = 0;
  for (int slot = 1; slot <= 16; slot++) {
    std::string hname = "TOP/good_hits_xy_" + to_string(slot);
    auto* h = static_cast<TH2F*>(findHist(hname));
    if (not h) continue;
 
    auto cuts = getDeadAndHotCuts(h);
    double deadCut = cuts.first;
    double hotCut = cuts.second;
    std::vector<int> pmts(32, 0);
    for (int row = 0; row < 8; row++) {
      for (int col = 0; col < 64; col++) {
        double y = h->GetBinContent(col + 1, row + 1);
        if (y > deadCut and y <= hotCut) {
          int pmt = col / 4 + (row / 4) * 16;
          pmts[pmt]++;
        }
      }
    }
    for (int n : pmts) if (n == 0) badPMTs++;
  }
  badPMTs -= badBoardstacks * 8;
 
  setEpicsPV("badBoardstacks", badBoardstacks);
  setEpicsPV("badCarriers", badCarriers);
  setEpicsPV("badAsics", badAsics);
  setEpicsPV("badPMTs", badPMTs);
  int numBS = 0;
  for (auto included : m_includedBoardstacks) if (not included) numBS++;
  setEpicsPV("numExcludedBS", numBS);
  setEpicsPV("histoAlarmState", getOffcialAlarmStatus(m_alarmStateOverall));
  setEpicsPV("backgroundAlarmLevels", m_averageRate);
 
  B2DEBUG(20, "badBoardstacks: " << badBoardstacks);
  B2DEBUG(20, "badCarriers: " << badCarriers);
  B2DEBUG(20, "badAsics: " << badAsics);
  B2DEBUG(20, "badPMTs: " << badPMTs);
  B2DEBUG(20, "excludedBS: " << numBS);
  B2DEBUG(20, "histoAlarmState: " << getOffcialAlarmStatus(m_alarmStateOverall));
  B2DEBUG(20, "backgroundAlarmLevels" << m_averageRate);
}

setEventProcessed()

static void setEventProcessed ( int e )

inlinestaticinherited

Set the number of processed events.

(Attention, asynch histogram updates!)

e Processed events.

Definition at line 233 of file DQMHistAnalysis.h.

{s_eventProcessed = e;};

setIncludedBoardstacks()

void setIncludedBoardstacks ( const std::vector< std::string > & excludedBoardstacks )

private

Sets flags for boardstacks to be included in alarming.

Parameters

excludedBoardstacks

list of boarstacks to be excluded from alarming

Definition at line 1186 of file DQMHistAnalysisTOP.cc.

{
  m_includedBoardstacks.clear();
  m_includedBoardstacks.resize(64, true);
 
  for (const auto& bsname : excludedBoardstacks) {
    int id = m_bsmap[bsname];
    if (id > 0) m_includedBoardstacks[id - 1] = false;
    else B2ERROR("Invalid boardstack name: " << bsname);
  }
}

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));
}

setMiraBelleVariables()

void setMiraBelleVariables ( const std::string & variableName, const TH1 * histogram )

private

Sets MiraBelle variables from the histogram with bins corresponding to slot numbers.

Parameters

variableName	variable name
histogram	histogram with bins corresponding to slot numbers

Definition at line 938 of file DQMHistAnalysisTOP.cc.

{
  for (int slot = 1; slot <= 16; slot++) {
    auto vname = variableName + std::to_string(slot);
    double value = histogram ? histogram->GetBinContent(slot) : 0;
    m_mirabelleVariables[vname] = value;
  }
}

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

t Run type string.

Definition at line 227 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()

void setUseEpics ( bool flag )

inlineinherited

Setter for EPICS usage.

Parameters

flag	set in use

Definition at line 581 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 587 of file DQMHistAnalysis.h.

{m_epicsReadOnly = flag;};

setZAxisRange()

void setZAxisRange ( const std::string & name, double scale )

staticprivate

Sets z-axis range of 2D histograms.

Parameters

name	the name of histograms
scale	scale factor for the histogram z-axis range (maximum = average * scale)

Definition at line 779 of file DQMHistAnalysisTOP.cc.

{
  double totalHits = 0;
  std::vector<TH2F*> histos;
  for (int slot = 1; slot <= 16; slot++) {
    TH2F* h = static_cast<TH2F*>(findHist(name + std::to_string(slot)));
    if (not h) continue;
    histos.push_back(h);
    totalHits += h->Integral();
  }
  if (histos.empty()) return;
  double average = totalHits / 512 / histos.size();  // per pixel or asic channel
 
  for (auto* h : histos) h->GetZaxis()->SetRangeUser(0, std::max(average * scale, 1.0));
}

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 302 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 324 of file DQMHistAnalysisTOP.cc.

{
  B2DEBUG(20, "terminate called");
}

updateBunchOffsetCanvas()

void updateBunchOffsetCanvas ( )

private

Updates canvas of bunch offset w/ alarming.

Definition at line 469 of file DQMHistAnalysisTOP.cc.

{
  int alarmState = c_Gray;
 
  auto* h = static_cast<TH1F*>(findHist("TOP/bunchOffset"));
  if (h) {
    double totalEvts = h->GetEntries();
    if (totalEvts > 100) {
      double mean = h->GetMean();
      double rms = h->GetRMS();
      alarmState = std::max(getAlarmState(fabs(mean), m_offsetMeanAlarmLevels), getAlarmState(rms, m_offsetRmsAlarmLevels));
    }
  }
 
  m_alarmStateOverall = std::max(m_alarmStateOverall, alarmState);
 
  auto* canvas = findCanvas("TOP/c_bunchOffset");
  if (canvas) {
    canvas->cd();
    canvas->Pad()->SetFrameFillColor(10);
    canvas->Pad()->SetFillColor(getAlarmColor(alarmState));
    canvas->Modified();
  }
}

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 354 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 359 of file DQMHistAnalysis.cc.

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

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 631 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;
}

updateEventMonitorCanvas()

void updateEventMonitorCanvas ( )

private

Updates canvas of event desynchronization monitor w/ alarming.

Definition at line 369 of file DQMHistAnalysisTOP.cc.

{
  int alarmState = c_Gray;
  m_text2->Clear();
 
  auto* evtMonitor = static_cast<TH1F*>(findHist("TOP/BoolEvtMonitor"));
  if (evtMonitor) {
    double totalEvts = evtMonitor->Integral();
    double badEvts = evtMonitor->GetBinContent(2);
    if (totalEvts > 0) {
      double badRatio = badEvts / totalEvts;
      alarmState = getAlarmState(badRatio, m_eventMonitorAlarmLevels);
      m_text2->AddText(Form("Fraction: %.4f %%", badRatio * 100.0));
    }
  }
 
  m_alarmStateOverall = std::max(m_alarmStateOverall, alarmState);
 
  auto* canvas = findCanvas("TOP/c_BoolEvtMonitor");
  if (canvas) {
    canvas->cd();
    m_text2->Draw();
    canvas->Pad()->SetFrameFillColor(10);
    canvas->Pad()->SetFillColor(getAlarmColor(alarmState));
    canvas->Modified();
  }
}

updateEventT0Canvas()

void updateEventT0Canvas ( )

private

Updates canvas of event T0 w/ alarming.

Definition at line 443 of file DQMHistAnalysisTOP.cc.

{
  int alarmState = c_Gray;
 
  auto* h = static_cast<TH1F*>(findHist("TOP/eventT0"));
  if (h) {
    double totalEvts = h->GetEntries();
    if (totalEvts > 100) {
      double mean = h->GetMean();
      double rms = h->GetRMS();
      alarmState = std::max(getAlarmState(fabs(mean), m_eventT0MeanAlarmLevels), getAlarmState(rms, m_eventT0RmsAlarmLevels));
    }
  }
 
  m_alarmStateOverall = std::max(m_alarmStateOverall, alarmState);
 
  auto* canvas = findCanvas("TOP/c_eventT0");
  if (canvas) {
    canvas->cd();
    canvas->Pad()->SetFrameFillColor(10);
    canvas->Pad()->SetFillColor(getAlarmColor(alarmState));
    canvas->Modified();
  }
}

updateLimits()

void updateLimits ( )

private

Updates limits defined by module parameters using EpicsPVs.

This includes asicWindowsBand, alarmLevels and excluded boardstacks.

Definition at line 1117 of file DQMHistAnalysisTOP.cc.

{
  double unused = 0;
 
  double yLo = m_asicWindowsBand[0];
  double yHi = m_asicWindowsBand[1];
  requestLimitsFromEpicsPVs("asicWindowsBand", yLo, unused, unused, yHi);
  m_asicWindowsBand[0] = yLo;
  m_asicWindowsBand[1] = yHi;
 
  requestLimitsFromEpicsPVs("asicWindowsAlarmLevels", unused, unused, m_asicWindowsAlarmLevels[0], m_asicWindowsAlarmLevels[1]);
  requestLimitsFromEpicsPVs("eventMonitorAlarmLevels", unused, unused, m_eventMonitorAlarmLevels[0], m_eventMonitorAlarmLevels[1]);
  requestLimitsFromEpicsPVs("junkHitsAlarmLevels", unused, unused, m_junkHitsAlarmLevels[0], m_junkHitsAlarmLevels[1]);
  requestLimitsFromEpicsPVs("deadChannelsAlarmLevels", unused, unused, m_deadChannelsAlarmLevels[0], m_deadChannelsAlarmLevels[1]);
  requestLimitsFromEpicsPVs("backgroundAlarmLevels", unused, unused, m_backgroundAlarmLevels[0], m_backgroundAlarmLevels[1]);
  requestLimitsFromEpicsPVs("photonYieldsAlarmLevels", m_photonYieldsAlarmLevels[0], m_photonYieldsAlarmLevels[1], unused, unused);
 
  requestLimitsFromEpicsPVs("injectionBGAlarmLevels", unused, unused, m_injectionBGAlarmLevels[0], m_injectionBGAlarmLevels[1]);
  requestLimitsFromEpicsPVs("timingAlarmLevels", unused, unused, m_timingAlarmLevels[0], m_timingAlarmLevels[1]);
  requestLimitsFromEpicsPVs("eventT0MeanAlarmLevels", unused, unused, m_eventT0MeanAlarmLevels[0], m_eventT0MeanAlarmLevels[1]);
  requestLimitsFromEpicsPVs("eventT0RmsAlarmLevels", unused, unused, m_eventT0RmsAlarmLevels[0], m_eventT0RmsAlarmLevels[1]);
  requestLimitsFromEpicsPVs("offsetMeanAlarmLevels", unused, unused, m_offsetMeanAlarmLevels[0], m_offsetMeanAlarmLevels[1]);
  requestLimitsFromEpicsPVs("offsetRmsAlarmLevels", unused, unused, m_offsetRmsAlarmLevels[0], m_offsetRmsAlarmLevels[1]);
 
  setAlarmLines();
 
  bool status = false;
  std::string excludedBS = getEpicsStringPV("excludedBoardstacks", status);
 
  if (status) {
    m_excludedBoardstacks.clear();
    std::string name;
    for (auto c : excludedBS) {
      if (isspace(c)) continue;
      else if (ispunct(c)) {
        if (not name.empty()) {
          m_excludedBoardstacks.push_back(name);
          name.clear();
        }
      } else name.push_back(c);
    }
    if (not name.empty()) {
      m_excludedBoardstacks.push_back(name);
    }
    setIncludedBoardstacks(m_excludedBoardstacks);
  }
 
  B2DEBUG(20, "asicWindowsBand:         [" << m_asicWindowsBand[0] << ", " << m_asicWindowsBand[1] << "]");
  B2DEBUG(20, "asicWindowsAlarmLevels:  [" << m_asicWindowsAlarmLevels[0] << ", " << m_asicWindowsAlarmLevels[1] << "]");
  B2DEBUG(20, "eventMonitorAlarmLevels: [" << m_eventMonitorAlarmLevels[0] << ", " << m_eventMonitorAlarmLevels[1] << "]");
  B2DEBUG(20, "junkHitsAlarmLevels:     [" << m_junkHitsAlarmLevels[0] << ", " << m_junkHitsAlarmLevels[1] << "]");
  B2DEBUG(20, "deadChannelsAlarmLevels: [" << m_deadChannelsAlarmLevels[0] << ", " << m_deadChannelsAlarmLevels[1] << "]");
  B2DEBUG(20, "backgroundAlarmLevels:   [" << m_backgroundAlarmLevels[0] << ", " << m_backgroundAlarmLevels[1] << "]");
  B2DEBUG(20, "photonYieldsAlarmLevels: [" << m_photonYieldsAlarmLevels[0] << ", " << m_photonYieldsAlarmLevels[1] << "]");
 
  B2DEBUG(20, "injectionBGAlarmLevels:     [" << m_injectionBGAlarmLevels[0] << ", " << m_injectionBGAlarmLevels[1] << "]");
  B2DEBUG(20, "timingAlarmLevels:     [" << m_timingAlarmLevels[0] << ", " << m_timingAlarmLevels[1] << "]");
  B2DEBUG(20, "eventT0MeanAlarmLevels:     [" << m_eventT0MeanAlarmLevels[0] << ", " << m_eventT0MeanAlarmLevels[1] << "]");
  B2DEBUG(20, "eventT0RmsAlarmLevels:     [" << m_eventT0RmsAlarmLevels[0] << ", " << m_eventT0RmsAlarmLevels[1] << "]");
  B2DEBUG(20, "offsetMeanAlarmLevels:     [" << m_offsetMeanAlarmLevels[0] << ", " << m_offsetMeanAlarmLevels[1] << "]");
  B2DEBUG(20, "offsetRmsAlarmLevels:     [" << m_offsetRmsAlarmLevels[0] << ", " << m_offsetRmsAlarmLevels[1] << "]");
 
  std::string ss;
  for (const auto& s : m_excludedBoardstacks) ss += "'" + s + "', ";
  if (ss.size() > 2)  {ss.pop_back(); ss.pop_back();}
  B2DEBUG(20, "excludedBoardstacks:     [" << ss << "]");
 
}

updateNGoodHitsCanvas()

void updateNGoodHitsCanvas ( )

private

Updates canvas of number of good hits per event w/ alarming (injection BG)

Definition at line 398 of file DQMHistAnalysisTOP.cc.

{
  int alarmState = c_Gray;
  m_text4->Clear();
 
  double fract = 0;
  double xcut = 0;
  double ymax = 0;
  auto* h = static_cast<TH1F*>(findHist("TOP/goodHitsPerEventAll"));
  if (h) {
    double totalEvts = h->GetEntries();
    if (totalEvts > 1000) {
      // fraction of events with more than xcut hits - these are mostly containing injection BG
      xcut = h->GetBinCenter(h->GetMaximumBin()) + 900;
      ymax = h->GetMaximum() / 2;
      fract = h->Integral(h->FindBin(xcut), h->GetNbinsX() + 1) / totalEvts * 100; // in %
      alarmState = getAlarmState(fract, m_injectionBGAlarmLevels);
      m_text4->AddText(Form("Events w/ Injection BG: %.2f %%", fract));
    }
  }
 
  setEpicsPV("injectionBGAlarmLevels", fract);
  m_alarmStateOverall = std::max(m_alarmStateOverall, alarmState);
 
  auto* canvas = findCanvas("TOP/c_goodHitsPerEventAll");
  if (canvas) {
    canvas->cd();
    if (not m_injBGCutLine) {
      m_injBGCutLine = new TLine(xcut, 0, xcut, ymax);
      m_injBGCutLine->SetLineWidth(2);
      m_injBGCutLine->SetLineColor(kRed);
      m_injBGCutLine->Draw("same");
    } else {
      m_injBGCutLine->SetX1(xcut);
      m_injBGCutLine->SetX2(xcut);
      m_injBGCutLine->SetY2(ymax);
    }
    m_text4->Draw();
    canvas->Pad()->SetFrameFillColor(10);
    canvas->Pad()->SetFillColor(getAlarmColor(alarmState));
    canvas->Modified();
  }
}

updateTimingCanvas()

void updateTimingCanvas ( )

private

Updates canvas of timing plot w/ alarming.

Definition at line 495 of file DQMHistAnalysisTOP.cc.

{
  int alarmState = c_Gray;
 
  auto* h = static_cast<TH1F*>(findHist("TOP/goodHitTimes"));
  auto* href = static_cast<TH1F*>(findRefHist("TOP/goodHitTimes"));
  if (h and href) {
    double n = h->Integral();
    double nref = href->Integral();
    if (n > 0 and nref > 0 and sameHistDefinition(h, href)) {
      auto* h_clone = static_cast<TH1F*>(h->Clone("tmp"));
      auto* href_clone = static_cast<TH1F*>(href->Clone("tmpref"));
      h_clone->Scale(1 / n);
      href_clone->Scale(1 / nref);
      h_clone->Add(h_clone, href_clone, 1, -1);
      double sumDiff = 0;
      double errDiff = 0;
      for (int i = 1; i <= h_clone->GetNbinsX(); i++) {
        sumDiff += fabs(h_clone->GetBinContent(i));
        errDiff += pow(h_clone->GetBinError(i), 2);
      }
      errDiff = sqrt(errDiff);
      if (sumDiff < 5 * errDiff) sumDiff = 0;  // difference not significant
      alarmState = getAlarmState(sumDiff, m_timingAlarmLevels);
      delete h_clone;
      delete href_clone;
    }
  }
 
  m_alarmStateOverall = std::max(m_alarmStateOverall, alarmState);
 
  auto* canvas = findCanvas("TOP/c_goodHitTimes");
  if (canvas) {
    canvas->cd();
    canvas->Pad()->SetFrameFillColor(10);
    canvas->Pad()->SetFillColor(getAlarmColor(alarmState));
    canvas->Modified();
  }
}

updateWindowVsSlotCanvas()

void updateWindowVsSlotCanvas ( )

private

Updates canvas of window_vs_slot w/ alarming.

Definition at line 330 of file DQMHistAnalysisTOP.cc.

{
  int alarmState = c_Gray;
  m_text1->Clear();
 
  auto* hraw = static_cast<TH2F*>(findHist("TOP/window_vs_slot"));
  if (hraw) {
    auto* px = hraw->ProjectionX("tmp_px");
    auto* band = hraw->ProjectionX("TOP/windowFractions", m_asicWindowsBand[0], m_asicWindowsBand[1]);
    band->Add(px, band, 1, -1);
    double total = px->Integral();
    double totalWindowFraction = (total != 0) ? band->Integral() / total : 0;
    band->Divide(band, px);
    setMiraBelleVariables("RateBadRaw_slot", band);
    m_mirabelleVariables["RateBadRaw_all"] = totalWindowFraction;
    if (total > 0) {
      alarmState = getAlarmState(totalWindowFraction, m_asicWindowsAlarmLevels);
      m_text1->AddText(Form("Fraction outside red lines: %.2f %%", totalWindowFraction * 100.0));
    }
    delete px;
    delete band;
  }
 
  m_alarmStateOverall = std::max(m_alarmStateOverall, alarmState);
 
  auto* canvas = findCanvas("TOP/c_window_vs_slot");
  if (canvas) {
    canvas->Clear();
    canvas->cd();
    if (hraw) hraw->Draw();
    m_text1->Draw();
    for (auto* line : m_asicWindowsBandLines) line->Draw();
    canvas->Pad()->SetFrameFillColor(10);
    canvas->Pad()->SetFillColor(getAlarmColor(alarmState));
    canvas->Modified();
  }
}

Member Data Documentation

m_activeFraction

TH1F* m_activeFraction = nullptr

private

fraction of active channels per slot

Definition at line 288 of file DQMHistAnalysisTOP.h.

m_alarmColors

std::vector<int> m_alarmColors = {c_ColorTooFew, c_ColorGood, c_ColorWarning, c_ColorError}

private

alarm colors (see base class)

Definition at line 267 of file DQMHistAnalysisTOP.h.

{c_ColorTooFew, c_ColorGood, c_ColorWarning, c_ColorError}; 

m_alarmStateOverall

int m_alarmStateOverall = 0

private

overall alarm state of histograms to be sent by EpicsPV

Definition at line 271 of file DQMHistAnalysisTOP.h.

m_asicWindowsAlarmLevels

std::vector<double> m_asicWindowsAlarmLevels = {0.002, 0.02}

private

alarm levels for fraction of windows outside the band

Definition at line 249 of file DQMHistAnalysisTOP.h.

{0.002, 0.02}; 

m_asicWindowsBand

std::vector<int> m_asicWindowsBand = {215, 235}

private

lower and upper bin of a band denoting good windows

Definition at line 248 of file DQMHistAnalysisTOP.h.

{215, 235}; 

m_asicWindowsBandLines

std::vector<TLine*> m_asicWindowsBandLines

private

lines denoting a band of good windows

Definition at line 311 of file DQMHistAnalysisTOP.h.

m_averageRate

double m_averageRate = 0

private

average BG rate (to pass to EpicsPV)

Definition at line 272 of file DQMHistAnalysisTOP.h.

m_backgroundAlarmLevels

std::vector<double> m_backgroundAlarmLevels = {5.0, 10.0}

private

alarm levels for background rates [MHz/PMT]

Definition at line 253 of file DQMHistAnalysisTOP.h.

{5.0, 10.0}; 

m_backgroundAlarmLines

std::vector<TLine*> m_backgroundAlarmLines

private

lines representing alarm levels

Definition at line 315 of file DQMHistAnalysisTOP.h.

m_backgroundRates

TH1D* m_backgroundRates = nullptr

private

background rates per slot

Definition at line 281 of file DQMHistAnalysisTOP.h.

m_bsmap

std::map<std::string, int> m_bsmap

private

a map of boardstack names to ID's

Definition at line 270 of file DQMHistAnalysisTOP.h.

m_c_backgroundRates

TCanvas* m_c_backgroundRates = nullptr

private

Canvas: background rates per slot.

Definition at line 283 of file DQMHistAnalysisTOP.h.

m_c_deadAndHot

TCanvas* m_c_deadAndHot = nullptr

private

Canvas: fractin of dead and hot channels.

Definition at line 291 of file DQMHistAnalysisTOP.h.

m_c_injBGs

std::map<std::string, TCanvas*> m_c_injBGs

private

Canvases for projections of injection BG histograms.

Definition at line 300 of file DQMHistAnalysisTOP.h.

m_c_injVetoFlagFract

TCanvas* m_c_injVetoFlagFract = nullptr

private

Canvas: fraction of events w/ injection veto flag set vs.

boardstack

Definition at line 307 of file DQMHistAnalysisTOP.h.

m_c_junkFraction

TCanvas* m_c_junkFraction = nullptr

private

Canvas: fraction of junk hits per boardstack.

Definition at line 295 of file DQMHistAnalysisTOP.h.

m_c_photonYields

TCanvas* m_c_photonYields = nullptr

private

Canvas: photon yields per slot.

Definition at line 282 of file DQMHistAnalysisTOP.h.

m_c_pmtHitRates

std::vector<TCanvas*> m_c_pmtHitRates

private

Canvases of PMT hits per event (index = slot - 1)

Definition at line 298 of file DQMHistAnalysisTOP.h.

m_c_skipProcFlagFract

TCanvas* m_c_skipProcFlagFract = nullptr

private

Canvas: fraction of events w/ skip processing flag set vs.

boardstack

Definition at line 306 of file DQMHistAnalysisTOP.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_deadChannelsAlarmLevels

std::vector<double> m_deadChannelsAlarmLevels = {0.1, 0.35}

private

alarm levels for the fraction of dead + hot channels

Definition at line 252 of file DQMHistAnalysisTOP.h.

{0.1, 0.35}; 

m_deadChannelsAlarmLines

std::vector<TLine*> m_deadChannelsAlarmLines

private

lines representing alarm levels

Definition at line 314 of file DQMHistAnalysisTOP.h.

m_deadFraction

TH1F* m_deadFraction = nullptr

private

fraction of dead channels per slot (included boardstacks only)

Definition at line 286 of file DQMHistAnalysisTOP.h.

m_description

std::string m_description

privateinherited

The description of the module.

Definition at line 510 of file Module.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 155 of file DQMHistAnalysis.h.

m_eventMonitorAlarmLevels

std::vector<double> m_eventMonitorAlarmLevels = {1e-4, 2e-3}

private

alarm levels for fraction of desynchronized digits

Definition at line 250 of file DQMHistAnalysisTOP.h.

{1e-4, 2e-3}; 

m_eventT0MeanAlarmLevels

std::vector<double> m_eventT0MeanAlarmLevels = {8, 20}

private

alarm levels for mean of event T0 [ns]

Definition at line 260 of file DQMHistAnalysisTOP.h.

{8, 20}; 

m_eventT0RmsAlarmLevels

std::vector<double> m_eventT0RmsAlarmLevels = {10, 20}

private

alarm levels for r.m.s.

of event T0 [ ns]

Definition at line 261 of file DQMHistAnalysisTOP.h.

{10, 20}; 

m_excludedBoardstacks

std::vector<std::string> m_excludedBoardstacks

private

list of boarstacks to be excluded from alarming

Definition at line 255 of file DQMHistAnalysisTOP.h.

m_excludedBSHisto

TH1F* m_excludedBSHisto = nullptr

private

histogram to show excluded boardstacks on junk fraction plot

Definition at line 294 of file DQMHistAnalysisTOP.h.

m_excludedFraction

TH1F* m_excludedFraction = nullptr

private

fraction of dead and hot channels per slot in excluded boardstacks only

Definition at line 287 of file DQMHistAnalysisTOP.h.

m_hasReturnValue

bool m_hasReturnValue

privateinherited

True, if the return value is set.

Definition at line 517 of file Module.h.

m_hotFraction

TH1F* m_hotFraction = nullptr

private

fraction of hot channels per slot (included boardstacks only)

Definition at line 285 of file DQMHistAnalysisTOP.h.

m_includedBoardstacks

std::vector<bool> m_includedBoardstacks

private

boardstacks included in alarming

Definition at line 269 of file DQMHistAnalysisTOP.h.

m_injBGCutLine

TLine* m_injBGCutLine = nullptr

private

a line denoting the cut on the number of hits for injection BG counting

Definition at line 317 of file DQMHistAnalysisTOP.h.

m_injectionBGAlarmLevels

std::vector<double> m_injectionBGAlarmLevels = {5, 10}

private

alarm levels for injection background (in % of events)

Definition at line 258 of file DQMHistAnalysisTOP.h.

{5, 10}; 

m_injVetoFlagFract

TH1D* m_injVetoFlagFract = nullptr

private

fraction of events w/ injection veto flag set vs.

boardstack

Definition at line 305 of file DQMHistAnalysisTOP.h.

m_IsNullRun

bool m_IsNullRun = false

private

Run type flag for null runs.

Definition at line 274 of file DQMHistAnalysisTOP.h.

m_junkFraction

TH1F* m_junkFraction = nullptr

private

fraction of junk hits per boardstack

Definition at line 293 of file DQMHistAnalysisTOP.h.

m_junkHitsAlarmLevels

std::vector<double> m_junkHitsAlarmLevels = {0.05, 0.25}

private

alarm levels for the fraction of junk hits

Definition at line 251 of file DQMHistAnalysisTOP.h.

{0.05, 0.25}; 

m_junkHitsAlarmLines

std::vector<TLine*> m_junkHitsAlarmLines

private

lines representing alarm levels

Definition at line 313 of file DQMHistAnalysisTOP.h.

m_legend

TLegend* m_legend = nullptr

private

legend for dead and hot channels

Definition at line 290 of file DQMHistAnalysisTOP.h.

m_logConfig

LogConfig m_logConfig

privateinherited

The log system configuration of the module.

Definition at line 513 of file Module.h.

m_mirabelleVariables

std::map<std::string, double> m_mirabelleVariables

private

variables for MiraBelle

Definition at line 326 of file DQMHistAnalysisTOP.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 = nullptr

private

MiraBelle monitoring object.

Definition at line 327 of file DQMHistAnalysisTOP.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_numEvents

double m_numEvents = 0

private

number of events processed with TOPDQM module

Definition at line 276 of file DQMHistAnalysisTOP.h.

m_officialStates

std::vector<int> m_officialStates = {c_StatusTooFew, c_StatusGood, c_StatusWarning, c_StatusError}

private

official alarm states

Definition at line 268 of file DQMHistAnalysisTOP.h.

{c_StatusTooFew, c_StatusGood, c_StatusWarning, c_StatusError}; 

m_offsetMeanAlarmLevels

std::vector<double> m_offsetMeanAlarmLevels = {0.2, 0.5}

private

alarm levels for mean of bunch offset [ns]

Definition at line 262 of file DQMHistAnalysisTOP.h.

{0.2, 0.5}; 

m_offsetRmsAlarmLevels

std::vector<double> m_offsetRmsAlarmLevels = {0.25, 0.50}

private

alarm levels for r.m.s.

of bunch offset [ns]

Definition at line 263 of file DQMHistAnalysisTOP.h.

{0.25, 0.50}; 

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_photonYields

TH1D* m_photonYields = nullptr

private

photon yields per slot

Definition at line 280 of file DQMHistAnalysisTOP.h.

m_photonYieldsAlarmLevels

std::vector<double> m_photonYieldsAlarmLevels = {15.0, 25.0}

private

alarm levels for the number of photons per track

Definition at line 254 of file DQMHistAnalysisTOP.h.

{15.0, 25.0}; 

m_photonYieldsAlarmLines

std::vector<TLine*> m_photonYieldsAlarmLines

private

lines representing alarm levels

Definition at line 316 of file DQMHistAnalysisTOP.h.

m_pmtHitRates

std::vector<TH1F*> m_pmtHitRates

private

histograms of PMT hits per event (index = slot - 1)

Definition at line 297 of file DQMHistAnalysisTOP.h.

m_profiles

std::map<std::string, TProfile*> m_profiles

private

profiles of injection BG

Definition at line 301 of file DQMHistAnalysisTOP.h.

m_projections

std::map<std::string, TH1D*> m_projections

private

projections of injection BG

Definition at line 302 of file DQMHistAnalysisTOP.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 160 of file DQMHistAnalysis.h.

m_pvPrefix

std::string m_pvPrefix

private

Epics PV prefix.

Definition at line 256 of file DQMHistAnalysisTOP.h.

m_returnValue

int m_returnValue

privateinherited

The return value.

Definition at line 518 of file Module.h.

m_runType

std::string m_runType

private

Run type.

Definition at line 275 of file DQMHistAnalysisTOP.h.

m_skipProcFlagFract

TH1D* m_skipProcFlagFract = nullptr

private

fraction of events w/ skip processing flag set vs.

boardstack

Definition at line 304 of file DQMHistAnalysisTOP.h.

m_stack

THStack* m_stack = nullptr

private

stack for drawing dead, hot and active channel fractions

Definition at line 289 of file DQMHistAnalysisTOP.h.

m_text1

TPaveText* m_text1 = nullptr

private

text to be written to window_vs_slot

Definition at line 319 of file DQMHistAnalysisTOP.h.

m_text2

TPaveText* m_text2 = nullptr

private

text to be written to event desynchonization monitor

Definition at line 320 of file DQMHistAnalysisTOP.h.

m_text3

TPaveText* m_text3 = nullptr

private

text to be written to background rates

Definition at line 321 of file DQMHistAnalysisTOP.h.

m_text4

TPaveText* m_text4 = nullptr

private

text to be written to number of good hits per event

Definition at line 322 of file DQMHistAnalysisTOP.h.

m_timingAlarmLevels

std::vector<double> m_timingAlarmLevels = {0.15, 0.30}

private

alarm levels for time distribution (fraction of area difference)

Definition at line 259 of file DQMHistAnalysisTOP.h.

{0.15, 0.30}; 

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 149 of file DQMHistAnalysis.h.

m_verticalLines

std::vector<TLine*> m_verticalLines

private

vertical lines splitting slots

Definition at line 312 of file DQMHistAnalysisTOP.h.

s_canvasUpdatedList

DQMHistAnalysisModule::CanvasUpdatedList s_canvasUpdatedList

staticprivateinherited

The list of canvas updated status.

Definition at line 124 of file DQMHistAnalysis.h.

s_deltaList

DQMHistAnalysisModule::DeltaList s_deltaList

staticinherited

The list of Delta Histograms and settings.

Definition at line 118 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 138 of file DQMHistAnalysis.h.

s_histList

DQMHistAnalysisModule::HistList s_histList

staticprivateinherited

The list of Histograms.

Definition at line 108 of file DQMHistAnalysis.h.

s_monObjList

DQMHistAnalysisModule::MonObjList s_monObjList

staticprivateinherited

The list of MonitoringObjects.

Definition at line 112 of file DQMHistAnalysis.h.

s_refList

DQMHistAnalysisModule::RefList s_refList

staticprivateinherited

The list of references.

Definition at line 129 of file DQMHistAnalysis.h.

s_runType

std::string s_runType = ""

inlinestaticprivateinherited

The Run type.

Definition at line 143 of file DQMHistAnalysis.h.

---

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

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