Make summary of data quality from reconstruction. More...

#include <DQMHistAnalysisMiraBelle.h>

Inheritance diagram for DQMHistAnalysisMiraBelleModule:

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
	DQMHistAnalysisMiraBelleModule ()
	Constructor.

TH1 *	findHistInCanvas (const std::string &hname, TCanvas **canvas=nullptr)
	Find histogram in corresponding canvas.

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

int	registerEpicsPV (const std::string &pvname, const std::string &keyname="")
	EPICS related Functions.

int	registerExternalEpicsPV (const std::string &pvname, const std::string &keyname="")
	Register a PV with its name and a key name.

void	setEpicsPV (const std::string &keyname, double value)
	Write value to a EPICS PV.

void	setEpicsPV (const std::string &keyname, int value)
	Write value to a EPICS PV.

void	setEpicsPV (int index, double value)
	Write value to a EPICS PV.

void	setEpicsPV (int index, int value)
	Write value to a EPICS PV.

void	setEpicsStringPV (const std::string &keyname, const std::string &value)
	Write string to a EPICS PV.

void	setEpicsStringPV (int index, const std::string &value)
	Write string to a EPICS PV.

double	getEpicsPV (const std::string &keyname)
	Read value from a EPICS PV.

double	getEpicsPV (int index)
	Read value from a EPICS PV.

std::string	getEpicsStringPV (const std::string &keyname, bool &status)
	Read value from a EPICS PV.

std::string	getEpicsStringPV (int index, bool &status)
	Read value from a EPICS PV.

int	updateEpicsPVs (float timeout)
	Update all EPICS PV (flush to network)

chid	getEpicsPVChID (const std::string &keyname)
	Get EPICS PV Channel Id.

chid	getEpicsPVChID (int index)
	Get EPICS PV Channel Id.

bool	requestLimitsFromEpicsPVs (chid id, double &lowerAlarm, double &lowerWarn, double &upperWarn, double &upperAlarm)
	Get Alarm Limits from EPICS PV.

bool	requestLimitsFromEpicsPVs (const std::string &keyname, double &lowerAlarm, double &lowerWarn, double &upperWarn, double &upperAlarm)
	Get Alarm Limits from EPICS PV.

bool	requestLimitsFromEpicsPVs (int index, double &lowerAlarm, double &lowerWarn, double &upperWarn, double &upperAlarm)
	Get Alarm Limits from EPICS PV.

void	cleanupEpicsPVs (void)
	Unsubscribe from EPICS PVs on terminate.

void	CheckEpicsError (int state, const std::string &message, const std::string &name)
	check the return status and check PV in case of error

virtual std::vector< std::string >	getFileNames (bool outputFiles)
	Return a list of output filenames for this modules.

const std::string &	getName () const
	Returns the name of the module.

const std::string &	getType () const
	Returns the type of the module (i.e.

const std::string &	getPackage () const
	Returns the package this module is in.

const std::string &	getDescription () const
	Returns the description of the module.

void	setName (const std::string &name)
	Set the name of the module.

void	setPropertyFlags (unsigned int propertyFlags)
	Sets the flags for the module properties.

LogConfig &	getLogConfig ()
	Returns the log system configuration.

void	setLogConfig (const LogConfig &logConfig)
	Set the log system configuration.

void	setLogLevel (int logLevel)
	Configure the log level.

void	setDebugLevel (int debugLevel)
	Configure the debug messaging level.

void	setAbortLevel (int abortLevel)
	Configure the abort log level.

void	setLogInfo (int logLevel, unsigned int logInfo)
	Configure the printed log information for the given level.

void	if_value (const std::string &expression, const std::shared_ptr< Path > &path, EAfterConditionPath afterConditionPath=EAfterConditionPath::c_End)
	Add a condition to the module.

void	if_false (const std::shared_ptr< Path > &path, EAfterConditionPath afterConditionPath=EAfterConditionPath::c_End)
	A simplified version to add a condition to the module.

void	if_true (const std::shared_ptr< Path > &path, EAfterConditionPath afterConditionPath=EAfterConditionPath::c_End)
	A simplified version to set the condition of the module.

bool	hasCondition () const
	Returns true if at least one condition was set for the module.

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

const std::vector< ModuleCondition > &	getAllConditions () const
	Return all set conditions for this module.

bool	evalCondition () const
	If at least one condition was set, it is evaluated and true returned if at least one condition returns true.

std::shared_ptr< Path >	getConditionPath () const
	Returns the path of the last true condition (if there is at least one, else reaturn a null pointer).

Module::EAfterConditionPath	getAfterConditionPath () const
	What to do after the conditional path is finished.

std::vector< std::shared_ptr< Path > >	getAllConditionPaths () const
	Return all condition paths currently set (no matter if the condition is true or not).

bool	hasProperties (unsigned int propertyFlags) const
	Returns true if all specified property flags are available in this module.

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

const ModuleParamList &	getParamList () const
	Return module param list.

template<typename T>
ModuleParam< T > &	getParam (const std::string &name) const
	Returns a reference to a parameter.

bool	hasReturnValue () const
	Return true if this module has a valid return value set.

int	getReturnValue () const
	Return the return value set by this module.

std::shared_ptr< PathElement >	clone () const override
	Create an independent copy of this module.

std::shared_ptr< boost::python::list >	getParamInfoListPython () const
	Returns a python list of all parameters.

Static Public Member Functions
static void	clearlist (void)
	Clear all static global lists.

static HistList &	getHistList ()
	Get the list of the histograms.

static MonObjList &	getMonObjList ()
	Get the list of MonitoringObjects.

static const DeltaList &	getDeltaList ()
	Get the list of the delta histograms.

static const CanvasUpdatedList &	getCanvasUpdatedList ()
	Get the list of the canvas update status.

static const std::string &	getRunType (void)
	Get the list of the reference histograms.

static int	getEventProcessed (void)
	Get the number of processed events.

static void	setRunType (const std::string &t)
	Set the Run Type.

static void	setEventProcessed (int e)
	Set the number of processed events.

static TCanvas *	findCanvas (TString cname)
	Find canvas by name.

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

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

static TH1 *	scaleReference (ERefScaling scaling, const TH1 hist, TH1 ref)
	Using the original and reference, create scaled version.

static TH1 *	findHistInFile (TFile *file, const std::string &histname)
	Find histogram in specific TFile (e.g.

static MonitoringObject *	findMonitoringObject (const std::string &objName)
	Find MonitoringObject.

static double	getSigma68 (TH1 *h)
	Helper function to compute half of the central interval covering 68% of a distribution.

static bool	addHist (const std::string &dirname, const std::string &histname, TH1 *h)
	Add histogram.

static void	addRefHist (const std::string &dirname, TH1 *hist)
	Add reference histogram.

static MonitoringObject *	getMonitoringObject (const std::string &name)
	Get MonitoringObject with given name (new object is created if non-existing)

static void	clearCanvases (void)
	Clear content of all Canvases.

static void	initHistListBeforeEvent (void)
	Reset the list of histograms.

static void	clearHistList (void)
	Clears the list of histograms.

static void	clearRefList (void)
	Clears the list of ref histograms.

static void	resetDeltaList (void)
	Reset Delta.

static void	addDeltaPar (const std::string &dirname, const std::string &histname, HistDelta::EDeltaType t, int p, unsigned int a=1)
	Add Delta histogram parameters.

static bool	hasDeltaPar (const std::string &dirname, const std::string &histname)
	Check if Delta histogram parameters exist for histogram.

static void	UpdateCanvas (const std::string &name, bool updated=true)
	Mark canvas as updated (or not)

static void	UpdateCanvas (TCanvas *canvas, bool updated=true)
	Mark canvas as updated (or not)

static void	ExtractRunType (std::vector< TH1 * > &hs)
	Extract Run Type from histogram title, called from input module.

static void	ExtractNEvent (std::vector< TH1 * > &hs)
	Extract event processed from daq histogram, called from input module.

static void	setUseEpics (bool flag)
	Setter for EPICS usage.

static void	setUseEpicsReadOnly (bool flag)
	Setter EPICS flag in read only mode.

static bool	getUseEpics (void)
	Getter for EPICS usage.

static bool	getUseEpicsReadOnly (void)
	Getter EPICS flag in read only mode.

static std::string &	getPVPrefix (void)
	get global prefix for EPICS PVs

static void	setPVPrefix (const std::string &prefix)
	set global prefix for EPICS PVs

static EStatus	makeStatus (bool enough, bool warn_flag, bool error_flag)
	Helper function to judge the status for coloring and EPICS.

static void	colorizeCanvas (TCanvas *canvas, EStatus status)
	Helper function for Canvas colorization.

static EStatusColor	getStatusColor (EStatus status)
	Return color for canvas state.

static void	checkPVStatus (void)
	Check the status of all PVs and report if disconnected or not found.

static void	printPVStatus (chid pv, bool onlyError=true)
	check the status of a PVs and report if disconnected or not found

static void	CheckEpicsError (int state, const std::string &message, chid id)
	check the return status and check PV in case of error

static std::vector< std::string >	StringSplit (const std::string &s, const char delim)
	Helper function for string token split.

static void	exposePythonAPI ()
	Exposes methods of the Module class to Python.

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

Protected Member Functions
virtual void	def_initialize ()
	Wrappers to make the methods without "def_" prefix callable from Python.

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

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

virtual void	def_endRun ()
	This method can receive that the current run ends as a call from the Python side.

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

void	setDescription (const std::string &description)
	Sets the description of the module.

void	setType (const std::string &type)
	Set the module type.

template<typename T>
void	addParam (const std::string &name, T &paramVariable, const std::string &description, const T &defaultValue)
	Adds a new parameter to the module.

template<typename T>
void	addParam (const std::string &name, T &paramVariable, const std::string &description)
	Adds a new enforced parameter to the module.

void	setReturnValue (int value)
	Sets the return value for this module as integer.

void	setReturnValue (bool value)
	Sets the return value for this module as bool.

void	setParamList (const ModuleParamList &params)
	Replace existing parameter list.

Protected Attributes
double	m_scale_dst
	Scale factor "signal region" / "sideband", assuming uniform events.

double	m_reference_hadb
	Reference value for hadronb2_tight bhabha ratio.

double	m_reference_mumu
	Reference value for mumu inv mass.

Belle2::MonitoringObject *	mon_mumu = nullptr
	MonitoringObject for mumu events.

Belle2::MonitoringObject *	mon_dst = nullptr
	MonitoringObject for D* events.

Belle2::MonitoringObject *	mon_tautau = nullptr
	MonitoringObject for tau pair events.

Belle2::MonitoringObject *	mon_bhabha = nullptr
	MonitoringObject for bhabha events.

Belle2::MonitoringObject *	mon_hadron = nullptr
	MonitoringObject for hadronb2 tight events.

TCanvas *	mumu_main = nullptr
	TCanvas object for main display.

TCanvas *	mumu_resolution = nullptr
	TCanvas object for resolution display.

TCanvas *	mumu_muon_val = nullptr
	TCanvas object for mumu values display.

TCanvas *	dst_mass = nullptr
	TCanvas object for D0 mass, delta mass, and pi0 mass.

TCanvas *	dst_pi_val = nullptr
	TCanvas object for PID distributions of pi.

TCanvas *	dst_k_val = nullptr
	TCanvas object for PID distributions of K.

TCanvas *	bhabha_main = nullptr
	TCanvas object for main display.

TCanvas *	bhabha_resolution = nullptr
	TCanvas object for resolution display.

TCanvas *	hadron_main = nullptr
	TCanvas object for main display.

Private Member Functions
void	initialize () override final
	Initialize the Module.

void	beginRun () override final
	Called when entering a new run.

void	event () override final
	Event processor.

void	endRun () override final
	End-of-run action.

void	terminate () override final
	Termination action.

void	ExtractMumu (void)
	Extract data from histograms related to MuMu.

void	ExtractD0Star (void)
	Extract data from histograms related to Do*.

void	ExtractTaupair (void)
	Extract data from histograms related to TauPairs.

void	ExtractBhabha (void)
	Extract data from histograms related to Bhabha.

void	ExtractHadron (void)
	Extract data from histograms related to Hadron.

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.

Private Attributes
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

Make summary of data quality from reconstruction.

Definition at line 26 of file DQMHistAnalysisMiraBelle.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

◆ 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

◆ DQMHistAnalysisMiraBelleModule()

DQMHistAnalysisMiraBelleModule ( )

Constructor.

Definition at line 24 of file DQMHistAnalysisMiraBelle.cc.

  : DQMHistAnalysisModule()
{
  // set module description (e.g. insert text)
  setDescription("Modify and analyze the data quality histograms of MiraBelle");
  setPropertyFlags(c_ParallelProcessingCertified);
  addParam("scale_dst", m_scale_dst, "Scale factor signal/sideband", 0.09375);
  addParam("m_reference_hadb", m_reference_hadb, "reference for hadronb2_tight bhabha ratio", 0.193);
  addParam("m_reference_mumu", m_reference_mumu, "reference for mumu inv mass", 10.568);
}

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 )

finaloverrideprivatevirtual

Called when entering a new run.

Set run dependent things like run header parameters, alignment, etc.

Reimplemented from Module.

Definition at line 75 of file DQMHistAnalysisMiraBelle.cc.

{
  B2DEBUG(20, "DQMHistAnalysisMiraBelle: beginRun called.");
}

◆ CheckEpicsError() [1/2]

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

staticinherited

check the return status and check PV in case of error

Parameters

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

Definition at line 793 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 785 of file DQMHistAnalysis.cc.

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

◆ checkPVStatus()

void checkPVStatus ( void )

staticinherited

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

Definition at line 749 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 630 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 295 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 323 of file DQMHistAnalysis.cc.

{
  s_histList.clear();
}

◆ clearlist()

void clearlist ( void )

staticinherited

Clear all static global lists.

Definition at line 50 of file DQMHistAnalysis.cc.

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

◆ clearRefList()

void clearRefList ( void )

staticinherited

Clears the list of ref histograms.

Definition at line 328 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 736 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.

425{ 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.

438{ 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.

431{ 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.

419{ 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.

444{ terminate(); }

◆ endRun()

void endRun ( void )

finaloverrideprivatevirtual

End-of-run action.

Save run-related stuff, such as statistics.

Reimplemented from Module.

Definition at line 1159 of file DQMHistAnalysisMiraBelle.cc.

{
  ExtractMumu();
  ExtractD0Star();
  ExtractTaupair();
  ExtractBhabha();
  ExtractHadron();
 
  B2DEBUG(20, "DQMHistAnalysisMiraBelle : 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 )

finaloverrideprivatevirtual

Event processor.

Reimplemented from Module.

Definition at line 80 of file DQMHistAnalysisMiraBelle.cc.

{
  B2DEBUG(20, "DQMHistAnalysisOutputMiraBelle: event called.");
}

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

◆ ExtractBhabha()

void ExtractBhabha ( void )

private

Extract data from histograms related to Bhabha.

Definition at line 870 of file DQMHistAnalysisMiraBelle.cc.

{
  //bhabha,hadrons
  // ========== bhabha_trk_ecl
  // get existing histograms produced by DQM modules
  auto* histbh_npxd = findHist("PhysicsObjectsMiraBelleBhabha/hist_npxd");
  auto* histbh_nsvd = findHist("PhysicsObjectsMiraBelleBhabha/hist_nsvd");
  auto* histbh_ncdc = findHist("PhysicsObjectsMiraBelleBhabha/hist_ncdc");
  auto* histbh_topdig = findHist("PhysicsObjectsMiraBelleBhabha/hist_topdig");
  auto* histbh_DetPhotonARICH = findHist("PhysicsObjectsMiraBelleBhabha/hist_DetPhotonARICH");
  auto* histbh_dD0 = findHist("PhysicsObjectsMiraBelleBhabha/hist_dD0");
  auto* histbh_dZ0 = findHist("PhysicsObjectsMiraBelleBhabha/hist_dZ0");
  auto* histbh_dPtcms = findHist("PhysicsObjectsMiraBelleBhabha/hist_dPtcms");
  auto* histbh_Pval = findHist("PhysicsObjectsMiraBelleBhabha/hist_Pval");
  auto* histbh_nExtraCDCHits = findHist("PhysicsObjectsMiraBelleBhabha/hist_nExtraCDCHits");
  auto* histbh_nECLClusters = findHist("PhysicsObjectsMiraBelleBhabha/hist_nECLClusters");
  auto* histbh_electronid = findHist("PhysicsObjectsMiraBelleBhabha/hist_electronid");
  auto* histbh_inv_p = findHist("PhysicsObjectsMiraBelleBhabha/hist_inv_p");
  auto* histbh_ndf = findHist("PhysicsObjectsMiraBelleBhabha/hist_ndf");
  auto* histbh_D0 = findHist("PhysicsObjectsMiraBelleBhabha/hist_D0");
  auto* histbh_Z0 = findHist("PhysicsObjectsMiraBelleBhabha/hist_Z0");
  auto* histbh_theta = findHist("PhysicsObjectsMiraBelleBhabha/hist_theta");
  auto* histbh_Phi0 = findHist("PhysicsObjectsMiraBelleBhabha/hist_Phi0");
  auto* histbh_Pt = findHist("PhysicsObjectsMiraBelleBhabha/hist_Pt");
  auto* histbh_Mom = findHist("PhysicsObjectsMiraBelleBhabha/hist_Mom");
  auto* histbh_dPhicms = findHist("PhysicsObjectsMiraBelleBhabha/hist_dPhicms");
 
  if (histbh_npxd == nullptr) {
    B2ERROR("Can not find the histbh_npxd histogram!");
    return;
  }
  if (histbh_nsvd == nullptr) {
    B2ERROR("Can not find the histbh_nsvd histogram!");
    return;
  }
  if (histbh_ncdc == nullptr) {
    B2ERROR("Can not find the histbh_ncdc histogram!");
    return;
  }
  if (histbh_topdig == nullptr) {
    B2ERROR("Can not find the histbh_topdig histogram!");
    return;
  }
  if (histbh_DetPhotonARICH == nullptr) {
    B2ERROR("Can not find the histbh_DetPhotonARICH histogram!");
    return;
  }
  if (histbh_dD0 == nullptr) {
    B2ERROR("Can not find the histbh_dD0 histogram!");
    return;
  }
  if (histbh_dZ0 == nullptr) {
    B2ERROR("Can not find the histbh_dZ0 histogram!");
    return;
  }
  if (histbh_dPtcms == nullptr) {
    B2ERROR("Can not find the histbh_dPtcms histogram!");
    return;
  }
  if (histbh_Pval == nullptr) {
    B2ERROR("Can not find the histbh_Pval histogram!");
    return;
  }
  if (histbh_nExtraCDCHits == nullptr) {
    B2ERROR("Can not find the histbh_nExtraCDCHits histogram!");
    return;
  }
  if (histbh_nECLClusters == nullptr) {
    B2ERROR("Can not find the histbh_nECLClusters histogram!");
    return;
  }
  if (histbh_electronid == nullptr) {
    B2ERROR("Can not find the histbh_electronid histogram!");
    return;
  }
  if (histbh_inv_p == nullptr) {
    B2ERROR("Can not find the histbh_inv_p histogram!");
    return;
  }
  if (histbh_ndf == nullptr) {
    B2ERROR("Can not find the histbh_ndf histogram!");
    return;
  }
  if (histbh_D0 == nullptr) {
    B2ERROR("Can not find the histbh_D0 histogram!");
    return;
  }
  if (histbh_Z0 == nullptr) {
    B2ERROR("Can not find the histbh_Z0 histogram!");
    return;
  }
  if (histbh_theta == nullptr) {
    B2ERROR("Can not find the histbh_theta histogram!");
    return;
  }
  if (histbh_Phi0 == nullptr) {
    B2ERROR("Can not find the histbh_Phi0 histogram!");
    return;
  }
  if (histbh_Pt == nullptr) {
    B2ERROR("Can not find the histbh_Pt histogram!");
    return;
  }
  if (histbh_Mom == nullptr) {
    B2ERROR("Can not find the histbh_Mom histogram!");
    return;
  }
  if (histbh_dPhicms == nullptr) {
    B2ERROR("Can not find the histbh_dPhicms histogram!");
    return;
  }
 
  // Make TCanvases
  // --- bhabha_Main
  bhabha_main->Divide(4, 3);
  bhabha_main->cd(1);  histbh_theta->Draw();
  bhabha_main->cd(2);  histbh_Phi0->Draw();
  bhabha_main->cd(3);  histbh_Mom->Draw();
  bhabha_main->cd(4);  histbh_Pt->Draw();
  bhabha_main->cd(5);  histbh_npxd->Draw();
  bhabha_main->cd(6);  histbh_nsvd->Draw();
  bhabha_main->cd(7);  histbh_ncdc->Draw();
  bhabha_main->cd(8);  histbh_topdig->Draw();
  bhabha_main->cd(9);  histbh_DetPhotonARICH->Draw();
  bhabha_main->cd(10);  histbh_nExtraCDCHits->Draw();
  bhabha_main->cd(11);  histbh_nECLClusters->Draw();
  // --- bhabha_Resolution
  bhabha_resolution->Divide(3, 3);
  bhabha_resolution->cd(1);  histbh_inv_p->Draw();
  bhabha_resolution->cd(2);  histbh_dD0->Draw();
  bhabha_resolution->cd(3);  histbh_dZ0->Draw();
  bhabha_resolution->cd(4);  histbh_dPtcms->Draw();
  bhabha_resolution->cd(5);  histbh_dPhicms->Draw();
  bhabha_resolution->cd(6);  histbh_nECLClusters->Draw();
  bhabha_resolution->cd(7);  histbh_nExtraCDCHits->Draw();
  bhabha_resolution->cd(8);  histbh_ndf->Draw();
 
  // calculate the values of monitoring variables
  double bh_mean_npxd = histbh_npxd->GetMean();
  double bh_mean_nsvd = histbh_nsvd->GetMean();
  double bh_mean_ncdc = histbh_ncdc->GetMean();
  double bh_mean_topdig = histbh_topdig->GetMean();
  double bh_mean_parich = histbh_DetPhotonARICH->GetMean();
  double bh_mean_ncdc_ex = histbh_nExtraCDCHits->GetMean();
  double bh_mean_necl = histbh_nECLClusters->GetMean();
  double bh_mean_electronid = histbh_electronid->GetMean();
  double bh_mean_d0 = histbh_D0->GetMean();
  double bh_mean_z0 = histbh_Z0->GetMean();
  double bh_mean_pval = histbh_Pval->GetMean();
  double bh_mean_ndf = histbh_ndf->GetMean();
  double bh_dif_ndf_ncdc = bh_mean_ndf - bh_mean_ncdc;
  double bh_mean_dd0 = histbh_dD0->GetMean();
  double bh_mean_dz0 = histbh_dZ0->GetMean();
  double bh_mean_dpt = histbh_dPtcms->GetMean();
  double bh_rms_dd0 = histbh_dD0->GetRMS();
  double bh_rms_dz0 = histbh_dZ0->GetRMS();
  double bh_rms_dpt = histbh_dPtcms->GetRMS();
  double bh_sigma68_dd0 = getSigma68(histbh_dD0);
  double bh_sigma68_dz0 = getSigma68(histbh_dZ0);
  double bh_sigma68_dpt = getSigma68(histbh_dPtcms);
  double bh_neve_bhabha = histbh_nECLClusters->GetEntries();
  double bh_neve_bhabha_sign = histbh_nsvd->GetEntries();
  double bh_goode_frac = -1.;
  double bh_pval_frac_0 = -1.;
  double bh_pval_frac_1 = -1.;
  double bh_nocdc_frac = -1.;
  double bh_notop_frac = -1.;
  double bh_noarich_frac = -1.;
  double eidcontent = histbh_electronid->GetBinContent(1);
  double bh_goode_o_bade;
  if (eidcontent != 0) {
    bh_goode_o_bade = histbh_electronid->GetBinContent(20) / eidcontent;
  } else {
    bh_goode_o_bade = 0.0;
  }
  double bh_pval_more95 = 0.0;
  double bh_pval_less05 = 0.0;
  for (int i = 95; i < 100; i++) bh_pval_more95 += histbh_Pval->GetBinContent(i + 1);
  for (int i = 0; i < 5; i++) bh_pval_less05 += histbh_Pval->GetBinContent(i + 1);
  if (bh_neve_bhabha_sign != 0) {
    bh_goode_frac = histbh_electronid->GetBinContent(20) / bh_neve_bhabha_sign;
    bh_pval_frac_0 = bh_pval_less05 / bh_neve_bhabha_sign;
    bh_pval_frac_1 = bh_pval_more95 / bh_neve_bhabha_sign;
    bh_nocdc_frac = histbh_ncdc->GetBinContent(1) / bh_neve_bhabha_sign;
    bh_notop_frac = histbh_topdig->GetBinContent(1) / bh_neve_bhabha_sign;
    bh_noarich_frac = histbh_DetPhotonARICH->GetBinContent(1) / bh_neve_bhabha_sign;
  } else {
    bh_goode_frac = 0.0;
    bh_pval_frac_0 = 0.0;
    bh_pval_frac_1 = 0.0;
    bh_nocdc_frac = 0.0;
    bh_notop_frac = 0.0;
    bh_noarich_frac = 0.0;
  }
  // set values
  mon_bhabha->setVariable("bh_mean_npxd", bh_mean_npxd);
  mon_bhabha->setVariable("bh_mean_nsvd", bh_mean_nsvd);
  mon_bhabha->setVariable("bh_mean_ncdc", bh_mean_ncdc);
  mon_bhabha->setVariable("bh_mean_topdig", bh_mean_topdig);
  mon_bhabha->setVariable("bh_mean_parich", bh_mean_parich);
  mon_bhabha->setVariable("bh_mean_ncdc_ex", bh_mean_ncdc_ex);
  mon_bhabha->setVariable("bh_mean_necl", bh_mean_necl);
  mon_bhabha->setVariable("bh_mean_electronid", bh_mean_electronid);
  mon_bhabha->setVariable("bh_mean_d0", bh_mean_d0);
  mon_bhabha->setVariable("bh_mean_z0", bh_mean_z0);
  mon_bhabha->setVariable("bh_mean_pval", bh_mean_pval);
  mon_bhabha->setVariable("bh_mean_ndf", bh_mean_ndf);
  mon_bhabha->setVariable("bh_dif_ndf_ncdc", bh_dif_ndf_ncdc);
  mon_bhabha->setVariable("bh_mean_dd0", bh_mean_dd0);
  mon_bhabha->setVariable("bh_mean_dz0", bh_mean_dz0);
  mon_bhabha->setVariable("bh_mean_dpt", bh_mean_dpt);
  mon_bhabha->setVariable("bh_rms_dd0", bh_rms_dd0);
  mon_bhabha->setVariable("bh_rms_dz0", bh_rms_dz0);
  mon_bhabha->setVariable("bh_rms_dpt", bh_rms_dpt);
  mon_bhabha->setVariable("bh_sigma68_dd0", bh_sigma68_dd0);
  mon_bhabha->setVariable("bh_sigma68_dz0", bh_sigma68_dz0);
  mon_bhabha->setVariable("bh_sigma68_dpt", bh_sigma68_dpt);
  mon_bhabha->setVariable("bh_neve_bhabha", bh_neve_bhabha);
  mon_bhabha->setVariable("bh_neve_bhabha_sign", bh_neve_bhabha_sign);
  mon_bhabha->setVariable("bh_goode_frac", bh_goode_frac);
  mon_bhabha->setVariable("bh_goode_o_bade", bh_goode_o_bade);
  mon_bhabha->setVariable("bh_pval_frac_0", bh_pval_frac_0);
  mon_bhabha->setVariable("bh_pval_frac_1", bh_pval_frac_1);
  mon_bhabha->setVariable("bh_nocdc_frac", bh_nocdc_frac);
  mon_bhabha->setVariable("bh_notop_frac", bh_notop_frac);
  mon_bhabha->setVariable("bh_noarich_frac", bh_noarich_frac);
}

◆ ExtractD0Star()

void ExtractD0Star ( void )

private

Extract data from histograms related to Do*.

Definition at line 376 of file DQMHistAnalysisMiraBelle.cc.

{
  // ========== D*
  // get existing histograms produced by DQM modules
  auto* hist_D0_InvM = findHist("PhysicsObjectsMiraBelleDst/hist_D0_InvM");
  auto* hist_delta_m = findHist("PhysicsObjectsMiraBelleDst/hist_delta_m");
  auto* hist_D0_softpi_PID_ALL_pion = findHist("PhysicsObjectsMiraBelleDst/hist_D0_softpi_PID_ALL_pion");
  auto* hist_D0_softpi_PID_SVD_pion = findHist("PhysicsObjectsMiraBelleDst/hist_D0_softpi_PID_SVD_pion");
  auto* hist_D0_softpi_PID_CDC_pion = findHist("PhysicsObjectsMiraBelleDst/hist_D0_softpi_PID_CDC_pion");
  auto* hist_D0_softpi_PID_TOP_pion = findHist("PhysicsObjectsMiraBelleDst/hist_D0_softpi_PID_TOP_pion");
  auto* hist_D0_softpi_PID_ARICH_pion = findHist("PhysicsObjectsMiraBelleDst/hist_D0_softpi_PID_ARICH_pion");
  auto* hist_D0_softpi_PID_ECL_pion = findHist("PhysicsObjectsMiraBelleDst/hist_D0_softpi_PID_ECL_pion");
  auto* hist_D0_softpi_PID_KLM_pion = findHist("PhysicsObjectsMiraBelleDst/hist_D0_softpi_PID_KLM_pion");
  auto* hist_D0_pi_PID_ALL_pion = findHist("PhysicsObjectsMiraBelleDst/hist_D0_pi_PID_ALL_pion");
  auto* hist_D0_pi_PID_SVD_pion = findHist("PhysicsObjectsMiraBelleDst/hist_D0_pi_PID_SVD_pion");
  auto* hist_D0_pi_PID_CDC_pion = findHist("PhysicsObjectsMiraBelleDst/hist_D0_pi_PID_CDC_pion");
  auto* hist_D0_pi_PID_TOP_pion = findHist("PhysicsObjectsMiraBelleDst/hist_D0_pi_PID_TOP_pion");
  auto* hist_D0_pi_PID_ARICH_pion = findHist("PhysicsObjectsMiraBelleDst/hist_D0_pi_PID_ARICH_pion");
  auto* hist_D0_pi_PID_ECL_pion = findHist("PhysicsObjectsMiraBelleDst/hist_D0_pi_PID_ECL_pion");
  auto* hist_D0_pi_PID_KLM_pion = findHist("PhysicsObjectsMiraBelleDst/hist_D0_pi_PID_KLM_pion");
  auto* hist_D0_K_PID_ALL_kaon = findHist("PhysicsObjectsMiraBelleDst/hist_D0_K_PID_ALL_kaon");
  auto* hist_D0_K_PID_SVD_kaon = findHist("PhysicsObjectsMiraBelleDst/hist_D0_K_PID_SVD_kaon");
  auto* hist_D0_K_PID_CDC_kaon = findHist("PhysicsObjectsMiraBelleDst/hist_D0_K_PID_CDC_kaon");
  auto* hist_D0_K_PID_TOP_kaon = findHist("PhysicsObjectsMiraBelleDst/hist_D0_K_PID_TOP_kaon");
  auto* hist_D0_K_PID_ARICH_kaon = findHist("PhysicsObjectsMiraBelleDst/hist_D0_K_PID_ARICH_kaon");
  auto* hist_D0_K_PID_ECL_kaon = findHist("PhysicsObjectsMiraBelleDst/hist_D0_K_PID_ECL_kaon");
  auto* hist_D0_K_PID_KLM_kaon = findHist("PhysicsObjectsMiraBelleDst/hist_D0_K_PID_KLM_kaon");
  auto* hist_sideband_D0_softpi_PID_ALL_pion = findHist("PhysicsObjectsMiraBelleDst/hist_sideband_D0_softpi_PID_ALL_pion");
  auto* hist_sideband_D0_softpi_PID_SVD_pion = findHist("PhysicsObjectsMiraBelleDst/hist_sideband_D0_softpi_PID_SVD_pion");
  auto* hist_sideband_D0_softpi_PID_CDC_pion = findHist("PhysicsObjectsMiraBelleDst/hist_sideband_D0_softpi_PID_CDC_pion");
  auto* hist_sideband_D0_softpi_PID_TOP_pion = findHist("PhysicsObjectsMiraBelleDst/hist_sideband_D0_softpi_PID_TOP_pion");
  auto* hist_sideband_D0_softpi_PID_ARICH_pion = findHist("PhysicsObjectsMiraBelleDst/hist_sideband_D0_softpi_PID_ARICH_pion");
  auto* hist_sideband_D0_softpi_PID_ECL_pion = findHist("PhysicsObjectsMiraBelleDst/hist_sideband_D0_softpi_PID_ECL_pion");
  auto* hist_sideband_D0_softpi_PID_KLM_pion = findHist("PhysicsObjectsMiraBelleDst/hist_sideband_D0_softpi_PID_KLM_pion");
  auto* hist_sideband_D0_pi_PID_ALL_pion = findHist("PhysicsObjectsMiraBelleDst/hist_sideband_D0_pi_PID_ALL_pion");
  auto* hist_sideband_D0_pi_PID_SVD_pion = findHist("PhysicsObjectsMiraBelleDst/hist_sideband_D0_pi_PID_SVD_pion");
  auto* hist_sideband_D0_pi_PID_CDC_pion = findHist("PhysicsObjectsMiraBelleDst/hist_sideband_D0_pi_PID_CDC_pion");
  auto* hist_sideband_D0_pi_PID_TOP_pion = findHist("PhysicsObjectsMiraBelleDst/hist_sideband_D0_pi_PID_TOP_pion");
  auto* hist_sideband_D0_pi_PID_ARICH_pion = findHist("PhysicsObjectsMiraBelleDst/hist_sideband_D0_pi_PID_ARICH_pion");
  auto* hist_sideband_D0_pi_PID_ECL_pion = findHist("PhysicsObjectsMiraBelleDst/hist_sideband_D0_pi_PID_ECL_pion");
  auto* hist_sideband_D0_pi_PID_KLM_pion = findHist("PhysicsObjectsMiraBelleDst/hist_sideband_D0_pi_PID_KLM_pion");
  auto* hist_sideband_D0_K_PID_ALL_kaon = findHist("PhysicsObjectsMiraBelleDst/hist_sideband_D0_K_PID_ALL_kaon");
  auto* hist_sideband_D0_K_PID_SVD_kaon = findHist("PhysicsObjectsMiraBelleDst/hist_sideband_D0_K_PID_SVD_kaon");
  auto* hist_sideband_D0_K_PID_CDC_kaon = findHist("PhysicsObjectsMiraBelleDst/hist_sideband_D0_K_PID_CDC_kaon");
  auto* hist_sideband_D0_K_PID_TOP_kaon = findHist("PhysicsObjectsMiraBelleDst/hist_sideband_D0_K_PID_TOP_kaon");
  auto* hist_sideband_D0_K_PID_ARICH_kaon = findHist("PhysicsObjectsMiraBelleDst/hist_sideband_D0_K_PID_ARICH_kaon");
  auto* hist_sideband_D0_K_PID_ECL_kaon = findHist("PhysicsObjectsMiraBelleDst/hist_sideband_D0_K_PID_ECL_kaon");
  auto* hist_sideband_D0_K_PID_KLM_kaon = findHist("PhysicsObjectsMiraBelleDst/hist_sideband_D0_K_PID_KLM_kaon");
  auto* hist_D0_pi0_InvM = findHist("PhysicsObjectsMiraBelleDst2/hist_D0_pi0_InvM");
  if (hist_D0_InvM == nullptr) {
    B2ERROR("Can not find the hist_D0_InvM histogram!");
    return;
  }
  if (hist_delta_m == nullptr) {
    B2ERROR("Can not find the hist_delta_m histogram!");
    return;
  }
  if (hist_D0_softpi_PID_ALL_pion == nullptr) {
    B2ERROR("Can not find the hist_D0_softpi_PID_ALL_pion histogram!");
    return;
  }
  if (hist_D0_softpi_PID_SVD_pion == nullptr) {
    B2ERROR("Can not find the hist_D0_softpi_PID_SVD_pion histogram!");
    return;
  }
  if (hist_D0_softpi_PID_CDC_pion == nullptr) {
    B2ERROR("Can not find the hist_D0_softpi_PID_CDC_pion histogram!");
    return;
  }
  if (hist_D0_softpi_PID_TOP_pion == nullptr) {
    B2ERROR("Can not find the hist_D0_softpi_PID_TOP_pion histogram!");
    return;
  }
  if (hist_D0_softpi_PID_ARICH_pion == nullptr) {
    B2ERROR("Can not find the hist_D0_softpi_PID_ARICH_pion histogram!");
    return;
  }
  if (hist_D0_softpi_PID_ECL_pion == nullptr) {
    B2ERROR("Can not find the hist_D0_softpi_PID_ECL_pion histogram!");
    return;
  }
  if (hist_D0_softpi_PID_KLM_pion == nullptr) {
    B2ERROR("Can not find the hist_D0_softpi_PID_KLM_pion histogram!");
    return;
  }
  if (hist_D0_pi_PID_ALL_pion == nullptr) {
    B2ERROR("Can not find the hist_D0_pi_PID_ALL_pion histogram!");
    return;
  }
  if (hist_D0_pi_PID_SVD_pion == nullptr) {
    B2ERROR("Can not find the hist_D0_pi_PID_SVD_pion histogram!");
    return;
  }
  if (hist_D0_pi_PID_CDC_pion == nullptr) {
    B2ERROR("Can not find the hist_D0_pi_PID_CDC_pion histogram!");
    return;
  }
  if (hist_D0_pi_PID_TOP_pion == nullptr) {
    B2ERROR("Can not find the hist_D0_pi_PID_TOP_pion histogram!");
    return;
  }
  if (hist_D0_pi_PID_ARICH_pion == nullptr) {
    B2ERROR("Can not find the hist_D0_pi_PID_ARICH_pion histogram!");
    return;
  }
  if (hist_D0_pi_PID_ECL_pion == nullptr) {
    B2ERROR("Can not find the hist_D0_pi_PID_ECL_pion histogram!");
    return;
  }
  if (hist_D0_pi_PID_KLM_pion == nullptr) {
    B2ERROR("Can not find the hist_D0_pi_PID_KLM_pion histogram!");
    return;
  }
  if (hist_D0_K_PID_ALL_kaon == nullptr) {
    B2ERROR("Can not find the hist_D0_K_PID_ALL_kaon histogram!");
    return;
  }
  if (hist_D0_K_PID_SVD_kaon == nullptr) {
    B2ERROR("Can not find the hist_D0_K_PID_SVD_kaon histogram!");
    return;
  }
  if (hist_D0_K_PID_CDC_kaon == nullptr) {
    B2ERROR("Can not find the hist_D0_K_PID_CDC_kaon histogram!");
    return;
  }
  if (hist_D0_K_PID_TOP_kaon == nullptr) {
    B2ERROR("Can not find the hist_D0_K_PID_TOP_kaon histogram!");
    return;
  }
  if (hist_D0_K_PID_ARICH_kaon == nullptr) {
    B2ERROR("Can not find the hist_D0_K_PID_ARICH_kaon histogram!");
    return;
  }
  if (hist_D0_K_PID_ECL_kaon == nullptr) {
    B2ERROR("Can not find the hist_D0_K_PID_ECL_kaon histogram!");
    return;
  }
  if (hist_D0_K_PID_KLM_kaon == nullptr) {
    B2ERROR("Can not find the hist_D0_K_PID_KLM_kaon histogram!");
    return;
  }
  if (hist_sideband_D0_softpi_PID_ALL_pion == nullptr) {
    B2ERROR("Can not find the hist_sideband_D0_softpi_PID_ALL_pion histogram!");
    return;
  }
  if (hist_sideband_D0_softpi_PID_SVD_pion == nullptr) {
    B2ERROR("Can not find the hist_sideband_D0_softpi_PID_SVD_pion histogram!");
    return;
  }
  if (hist_sideband_D0_softpi_PID_CDC_pion == nullptr) {
    B2ERROR("Can not find the hist_sideband_D0_softpi_PID_CDC_pion histogram!");
    return;
  }
  if (hist_sideband_D0_softpi_PID_TOP_pion == nullptr) {
    B2ERROR("Can not find the hist_sideband_D0_softpi_PID_TOP_pion histogram!");
    return;
  }
  if (hist_sideband_D0_softpi_PID_ARICH_pion == nullptr) {
    B2ERROR("Can not find the hist_sideband_D0_softpi_PID_ARICH_pion histogram!");
    return;
  }
  if (hist_sideband_D0_softpi_PID_ECL_pion == nullptr) {
    B2ERROR("Can not find the hist_sideband_D0_softpi_PID_ECL_pion histogram!");
    return;
  }
  if (hist_sideband_D0_softpi_PID_KLM_pion == nullptr) {
    B2ERROR("Can not find the hist_sideband_D0_softpi_PID_KLM_pion histogram!");
    return;
  }
  if (hist_sideband_D0_pi_PID_ALL_pion == nullptr) {
    B2ERROR("Can not find the hist_sideband_D0_pi_PID_ALL_pion histogram!");
    return;
  }
  if (hist_sideband_D0_pi_PID_SVD_pion == nullptr) {
    B2ERROR("Can not find the hist_sideband_D0_pi_PID_SVD_pion histogram!");
    return;
  }
  if (hist_sideband_D0_pi_PID_CDC_pion == nullptr) {
    B2ERROR("Can not find the hist_sideband_D0_pi_PID_CDC_pion histogram!");
    return;
  }
  if (hist_sideband_D0_pi_PID_TOP_pion == nullptr) {
    B2ERROR("Can not find the hist_sideband_D0_pi_PID_TOP_pion histogram!");
    return;
  }
  if (hist_sideband_D0_pi_PID_ARICH_pion == nullptr) {
    B2ERROR("Can not find the hist_sideband_D0_pi_PID_ARICH_pion histogram!");
    return;
  }
  if (hist_sideband_D0_pi_PID_ECL_pion == nullptr) {
    B2ERROR("Can not find the hist_sideband_D0_pi_PID_ECL_pion histogram!");
    return;
  }
  if (hist_sideband_D0_pi_PID_KLM_pion == nullptr) {
    B2ERROR("Can not find the hist_sideband_D0_pi_PID_KLM_pion histogram!");
    return;
  }
  if (hist_sideband_D0_K_PID_ALL_kaon == nullptr) {
    B2ERROR("Can not find the hist_sideband_D0_K_PID_ALL_kaon histogram!");
    return;
  }
  if (hist_sideband_D0_K_PID_SVD_kaon == nullptr) {
    B2ERROR("Can not find the hist_sideband_D0_K_PID_SVD_kaon histogram!");
    return;
  }
  if (hist_sideband_D0_K_PID_CDC_kaon == nullptr) {
    B2ERROR("Can not find the hist_sideband_D0_K_PID_CDC_kaon histogram!");
    return;
  }
  if (hist_sideband_D0_K_PID_TOP_kaon == nullptr) {
    B2ERROR("Can not find the hist_sideband_D0_K_PID_TOP_kaon histogram!");
    return;
  }
  if (hist_sideband_D0_K_PID_ARICH_kaon == nullptr) {
    B2ERROR("Can not find the hist_sideband_D0_K_PID_ARICH_kaon histogram!");
    return;
  }
  if (hist_sideband_D0_K_PID_ECL_kaon == nullptr) {
    B2ERROR("Can not find the hist_sideband_D0_K_PID_ECL_kaon histogram!");
    return;
  }
  if (hist_sideband_D0_K_PID_KLM_kaon == nullptr) {
    B2ERROR("Can not find the hist_sideband_D0_K_PID_KLM_kaon histogram!");
    return;
  }
  if (hist_D0_pi0_InvM == nullptr) {
    B2ERROR("Can not find the hist_D0_pi0_InvM histogram!");
    return;
  }
  // Fit mass distributions for scale factor
  TF1* f_InvM = new TF1("f_InvM", "[0]*TMath::Gaus(x,[1],[2])+[3]*pow(x-[4],2)+[5]", 1.81, 1.95);
  f_InvM->SetParNames("Height", "#mu", "#sigma", "a", "b", "c");
  f_InvM->SetParameters(hist_D0_InvM->GetMaximum(), 1.86, 5e-3, 0., 1.86, 0.);
  hist_D0_InvM->Fit(f_InvM, "", "", 1.81, 1.95);
  f_InvM->SetLineColor(kRed);
 
  TF1* f_delta_m = new TF1("f_delta_m", "[0]*TMath::Gaus(x,[1],[2])+[3]*pow(x-[4],2)+[5]", 0.14, 0.16);
  f_delta_m->SetParNames("Height", "#mu", "#sigma", "a", "b", "c");
  f_delta_m->SetParameters(hist_delta_m->GetMaximum(), 0.145, 5e-4, 0., 0.145, 0.);
  hist_delta_m->Fit(f_delta_m, "", "", 0.14, 0.16);
  f_delta_m->SetLineColor(kRed);
 
  TF1* f_pi0_InvM = new TF1("f_pi0_InvM", "[0]*TMath::Gaus(x,[1],[2])+[3]*pow(x-[4],2)+[5]", 0.09, 0.17);
  f_pi0_InvM->SetParNames("Height", "#mu", "#sigma", "a", "b", "c");
  f_pi0_InvM->SetParameters(hist_D0_pi0_InvM->GetMaximum(), 0.13, 5e-3, 0., 0.13, 0.);
  hist_D0_pi0_InvM->Fit(f_pi0_InvM, "", "", 0.09, 0.17);
  f_pi0_InvM->SetLineColor(kRed);
 
  // D->KPi and D->KPiPi0 yields
  TF1* f_gaus_InvM = new TF1("f_gaus_InvM", "gaus", 1.81, 1.95);
  f_gaus_InvM->SetParameters(f_InvM->GetParameter(0), f_InvM->GetParameter(1), f_InvM->GetParameter(2));
  double InvM_bin_width = hist_D0_InvM->GetXaxis()->GetBinWidth(1);
  double neve_dst = f_gaus_InvM->Integral(1.81, 1.95) / InvM_bin_width;
 
  TF1* f_gaus_pi0_InvM = new TF1("f_gaus_pi0_InvM", "gaus", 0.09, 0.17);
  f_gaus_pi0_InvM->SetParameters(f_pi0_InvM->GetParameter(0), f_pi0_InvM->GetParameter(1), f_pi0_InvM->GetParameter(2));
  double pi0_InvM_bin_width = hist_D0_pi0_InvM->GetXaxis()->GetBinWidth(1);
  double neve_pi0 = f_gaus_pi0_InvM->Integral(0.09, 0.17) / pi0_InvM_bin_width;
 
  // Sumw2
  hist_D0_softpi_PID_ALL_pion->Sumw2();
  hist_D0_softpi_PID_SVD_pion->Sumw2();
  hist_D0_softpi_PID_CDC_pion->Sumw2();
  hist_D0_softpi_PID_TOP_pion->Sumw2();
  hist_D0_softpi_PID_ARICH_pion->Sumw2();
  hist_D0_softpi_PID_ECL_pion->Sumw2();
  hist_D0_softpi_PID_KLM_pion->Sumw2();
  hist_D0_pi_PID_ALL_pion->Sumw2();
  hist_D0_pi_PID_SVD_pion->Sumw2();
  hist_D0_pi_PID_CDC_pion->Sumw2();
  hist_D0_pi_PID_TOP_pion->Sumw2();
  hist_D0_pi_PID_ARICH_pion->Sumw2();
  hist_D0_pi_PID_ECL_pion->Sumw2();
  hist_D0_pi_PID_KLM_pion->Sumw2();
  hist_D0_K_PID_ALL_kaon->Sumw2();
  hist_D0_K_PID_SVD_kaon->Sumw2();
  hist_D0_K_PID_CDC_kaon->Sumw2();
  hist_D0_K_PID_TOP_kaon->Sumw2();
  hist_D0_K_PID_ARICH_kaon->Sumw2();
  hist_D0_K_PID_ECL_kaon->Sumw2();
  hist_D0_K_PID_KLM_kaon->Sumw2();
  hist_sideband_D0_softpi_PID_ALL_pion->Sumw2();
  hist_sideband_D0_softpi_PID_SVD_pion->Sumw2();
  hist_sideband_D0_softpi_PID_CDC_pion->Sumw2();
  hist_sideband_D0_softpi_PID_TOP_pion->Sumw2();
  hist_sideband_D0_softpi_PID_ARICH_pion->Sumw2();
  hist_sideband_D0_softpi_PID_ECL_pion->Sumw2();
  hist_sideband_D0_softpi_PID_KLM_pion->Sumw2();
  hist_sideband_D0_pi_PID_ALL_pion->Sumw2();
  hist_sideband_D0_pi_PID_SVD_pion->Sumw2();
  hist_sideband_D0_pi_PID_CDC_pion->Sumw2();
  hist_sideband_D0_pi_PID_TOP_pion->Sumw2();
  hist_sideband_D0_pi_PID_ARICH_pion->Sumw2();
  hist_sideband_D0_pi_PID_ECL_pion->Sumw2();
  hist_sideband_D0_pi_PID_KLM_pion->Sumw2();
  hist_sideband_D0_K_PID_ALL_kaon->Sumw2();
  hist_sideband_D0_K_PID_SVD_kaon->Sumw2();
  hist_sideband_D0_K_PID_CDC_kaon->Sumw2();
  hist_sideband_D0_K_PID_TOP_kaon->Sumw2();
  hist_sideband_D0_K_PID_ARICH_kaon->Sumw2();
  hist_sideband_D0_K_PID_ECL_kaon->Sumw2();
  hist_sideband_D0_K_PID_KLM_kaon->Sumw2();
 
  // Scale the distributions in sideband
  hist_sideband_D0_softpi_PID_ALL_pion->Scale(m_scale_dst);
  hist_sideband_D0_softpi_PID_SVD_pion->Scale(m_scale_dst);
  hist_sideband_D0_softpi_PID_CDC_pion->Scale(m_scale_dst);
  hist_sideband_D0_softpi_PID_TOP_pion->Scale(m_scale_dst);
  hist_sideband_D0_softpi_PID_ARICH_pion->Scale(m_scale_dst);
  hist_sideband_D0_softpi_PID_ECL_pion->Scale(m_scale_dst);
  hist_sideband_D0_softpi_PID_KLM_pion->Scale(m_scale_dst);
  hist_sideband_D0_pi_PID_ALL_pion->Scale(m_scale_dst);
  hist_sideband_D0_pi_PID_SVD_pion->Scale(m_scale_dst);
  hist_sideband_D0_pi_PID_CDC_pion->Scale(m_scale_dst);
  hist_sideband_D0_pi_PID_TOP_pion->Scale(m_scale_dst);
  hist_sideband_D0_pi_PID_ARICH_pion->Scale(m_scale_dst);
  hist_sideband_D0_pi_PID_ECL_pion->Scale(m_scale_dst);
  hist_sideband_D0_pi_PID_KLM_pion->Scale(m_scale_dst);
  hist_sideband_D0_K_PID_ALL_kaon->Scale(m_scale_dst);
  hist_sideband_D0_K_PID_SVD_kaon->Scale(m_scale_dst);
  hist_sideband_D0_K_PID_CDC_kaon->Scale(m_scale_dst);
  hist_sideband_D0_K_PID_TOP_kaon->Scale(m_scale_dst);
  hist_sideband_D0_K_PID_ARICH_kaon->Scale(m_scale_dst);
  hist_sideband_D0_K_PID_ECL_kaon->Scale(m_scale_dst);
  hist_sideband_D0_K_PID_KLM_kaon->Scale(m_scale_dst);
 
  // BG subtraction
  hist_D0_softpi_PID_ALL_pion->Add(hist_sideband_D0_softpi_PID_ALL_pion, -1);
  hist_D0_softpi_PID_SVD_pion->Add(hist_sideband_D0_softpi_PID_SVD_pion, -1);
  hist_D0_softpi_PID_CDC_pion->Add(hist_sideband_D0_softpi_PID_CDC_pion, -1);
  hist_D0_softpi_PID_TOP_pion->Add(hist_sideband_D0_softpi_PID_TOP_pion, -1);
  hist_D0_softpi_PID_ARICH_pion->Add(hist_sideband_D0_softpi_PID_ARICH_pion, -1);
  hist_D0_softpi_PID_ECL_pion->Add(hist_sideband_D0_softpi_PID_ECL_pion, -1);
  hist_D0_softpi_PID_KLM_pion->Add(hist_sideband_D0_softpi_PID_KLM_pion, -1);
  hist_D0_pi_PID_ALL_pion->Add(hist_sideband_D0_pi_PID_ALL_pion, -1);
  hist_D0_pi_PID_SVD_pion->Add(hist_sideband_D0_pi_PID_SVD_pion, -1);
  hist_D0_pi_PID_CDC_pion->Add(hist_sideband_D0_pi_PID_CDC_pion, -1);
  hist_D0_pi_PID_TOP_pion->Add(hist_sideband_D0_pi_PID_TOP_pion, -1);
  hist_D0_pi_PID_ARICH_pion->Add(hist_sideband_D0_pi_PID_ARICH_pion, -1);
  hist_D0_pi_PID_ECL_pion->Add(hist_sideband_D0_pi_PID_ECL_pion, -1);
  hist_D0_pi_PID_KLM_pion->Add(hist_sideband_D0_pi_PID_KLM_pion, -1);
  hist_D0_K_PID_ALL_kaon->Add(hist_sideband_D0_K_PID_ALL_kaon, -1);
  hist_D0_K_PID_SVD_kaon->Add(hist_sideband_D0_K_PID_SVD_kaon, -1);
  hist_D0_K_PID_CDC_kaon->Add(hist_sideband_D0_K_PID_CDC_kaon, -1);
  hist_D0_K_PID_TOP_kaon->Add(hist_sideband_D0_K_PID_TOP_kaon, -1);
  hist_D0_K_PID_ARICH_kaon->Add(hist_sideband_D0_K_PID_ARICH_kaon, -1);
  hist_D0_K_PID_ECL_kaon->Add(hist_sideband_D0_K_PID_ECL_kaon, -1);
  hist_D0_K_PID_KLM_kaon->Add(hist_sideband_D0_K_PID_KLM_kaon, -1);
 
  // Make TCanvases
  // --- D* -> D0pi mass resolution
  dst_mass->Divide(3, 1);
  dst_mass->cd(1);
  hist_D0_InvM->Draw();
  f_InvM->Draw("SAME");
  dst_mass->cd(2);
  hist_delta_m->Draw();
  f_delta_m->Draw("SAME");
  dst_mass->cd(3);
  hist_D0_pi0_InvM->Draw();
  f_pi0_InvM->Draw("SAME");
  // --- pi variables for D0 -> K pi
  dst_pi_val->Divide(4, 2);
  dst_pi_val->cd(1);  hist_D0_pi_PID_ALL_pion->Draw();
  dst_pi_val->cd(2);  hist_D0_pi_PID_SVD_pion->Draw();
  dst_pi_val->cd(3);  hist_D0_pi_PID_CDC_pion->Draw();
  dst_pi_val->cd(4);  hist_D0_pi_PID_TOP_pion->Draw();
  dst_pi_val->cd(5);  hist_D0_pi_PID_ARICH_pion->Draw();
  dst_pi_val->cd(6);  hist_D0_pi_PID_ECL_pion->Draw();
  dst_pi_val->cd(7);  hist_D0_pi_PID_KLM_pion->Draw();
  // --- K variables for D0 -> K pi
  dst_k_val->Divide(4, 2);
  dst_k_val->cd(1);  hist_D0_K_PID_ALL_kaon->Draw();
  dst_k_val->cd(2);  hist_D0_K_PID_SVD_kaon->Draw();
  dst_k_val->cd(3);  hist_D0_K_PID_CDC_kaon->Draw();
  dst_k_val->cd(4);  hist_D0_K_PID_TOP_kaon->Draw();
  dst_k_val->cd(5);  hist_D0_K_PID_ARICH_kaon->Draw();
  dst_k_val->cd(6);  hist_D0_K_PID_ECL_kaon->Draw();
  dst_k_val->cd(7);  hist_D0_K_PID_KLM_kaon->Draw();
 
  // calculate the values of monitoring variables
  double mean_D0_InvM = f_InvM->GetParameter(1);
  double width_D0_InvM = f_InvM->GetParameter(2);
  double mean_delta_m = f_delta_m->GetParameter(1);
  double width_delta_m = f_delta_m->GetParameter(2);
  double mean_pi0_InvM = f_pi0_InvM->GetParameter(1);
  double width_pi0_InvM = f_pi0_InvM->GetParameter(2);
  double mean_D0_softpi_PID_ALL_pion = hist_D0_softpi_PID_ALL_pion->GetMean();
  double mean_D0_softpi_PID_SVD_pion = hist_D0_softpi_PID_SVD_pion->GetMean();
  double mean_D0_softpi_PID_CDC_pion = hist_D0_softpi_PID_CDC_pion->GetMean();
  double mean_D0_softpi_PID_TOP_pion = hist_D0_softpi_PID_TOP_pion->GetMean();
  double mean_D0_softpi_PID_ARICH_pion = hist_D0_softpi_PID_ARICH_pion->GetMean();
  double mean_D0_softpi_PID_ECL_pion = hist_D0_softpi_PID_ECL_pion->GetMean();
  double mean_D0_softpi_PID_KLM_pion = hist_D0_softpi_PID_KLM_pion->GetMean();
  double mean_D0_pi_PID_ALL_pion = hist_D0_pi_PID_ALL_pion->GetMean();
  double mean_D0_pi_PID_SVD_pion = hist_D0_pi_PID_SVD_pion->GetMean();
  double mean_D0_pi_PID_CDC_pion = hist_D0_pi_PID_CDC_pion->GetMean();
  double mean_D0_pi_PID_TOP_pion = hist_D0_pi_PID_TOP_pion->GetMean();
  double mean_D0_pi_PID_ARICH_pion = hist_D0_pi_PID_ARICH_pion->GetMean();
  double mean_D0_pi_PID_ECL_pion = hist_D0_pi_PID_ECL_pion->GetMean();
  double mean_D0_pi_PID_KLM_pion = hist_D0_pi_PID_KLM_pion->GetMean();
  double mean_D0_K_PID_ALL_kaon = hist_D0_K_PID_ALL_kaon->GetMean();
  double mean_D0_K_PID_SVD_kaon = hist_D0_K_PID_SVD_kaon->GetMean();
  double mean_D0_K_PID_CDC_kaon = hist_D0_K_PID_CDC_kaon->GetMean();
  double mean_D0_K_PID_TOP_kaon = hist_D0_K_PID_TOP_kaon->GetMean();
  double mean_D0_K_PID_ARICH_kaon = hist_D0_K_PID_ARICH_kaon->GetMean();
  double mean_D0_K_PID_ECL_kaon = hist_D0_K_PID_ECL_kaon->GetMean();
  double mean_D0_K_PID_KLM_kaon = hist_D0_K_PID_KLM_kaon->GetMean();
 
  // set values
  mon_dst->setVariable("neve_dst", neve_dst);
  mon_dst->setVariable("neve_pi0", neve_pi0);
  mon_dst->setVariable("mean_D0_InvM", mean_D0_InvM);
  mon_dst->setVariable("width_D0_InvM", width_D0_InvM);
  mon_dst->setVariable("mean_delta_m", mean_delta_m);
  mon_dst->setVariable("width_delta_m", width_delta_m);
  mon_dst->setVariable("mean_pi0_InvM", mean_pi0_InvM);
  mon_dst->setVariable("width_pi0_InvM", width_pi0_InvM);
  mon_dst->setVariable("mean_D0_softpi_PID_ALL_pion", mean_D0_softpi_PID_ALL_pion);
  mon_dst->setVariable("mean_D0_softpi_PID_SVD_pion", mean_D0_softpi_PID_SVD_pion);
  mon_dst->setVariable("mean_D0_softpi_PID_CDC_pion", mean_D0_softpi_PID_CDC_pion);
  mon_dst->setVariable("mean_D0_softpi_PID_TOP_pion", mean_D0_softpi_PID_TOP_pion);
  mon_dst->setVariable("mean_D0_softpi_PID_ARICH_pion", mean_D0_softpi_PID_ARICH_pion);
  mon_dst->setVariable("mean_D0_softpi_PID_ECL_pion", mean_D0_softpi_PID_ECL_pion);
  mon_dst->setVariable("mean_D0_softpi_PID_KLM_pion", mean_D0_softpi_PID_KLM_pion);
  mon_dst->setVariable("mean_D0_pi_PID_ALL_pion", mean_D0_pi_PID_ALL_pion);
  mon_dst->setVariable("mean_D0_pi_PID_SVD_pion", mean_D0_pi_PID_SVD_pion);
  mon_dst->setVariable("mean_D0_pi_PID_CDC_pion", mean_D0_pi_PID_CDC_pion);
  mon_dst->setVariable("mean_D0_pi_PID_TOP_pion", mean_D0_pi_PID_TOP_pion);
  mon_dst->setVariable("mean_D0_pi_PID_ARICH_pion", mean_D0_pi_PID_ARICH_pion);
  mon_dst->setVariable("mean_D0_pi_PID_ECL_pion", mean_D0_pi_PID_ECL_pion);
  mon_dst->setVariable("mean_D0_pi_PID_KLM_pion", mean_D0_pi_PID_KLM_pion);
  mon_dst->setVariable("mean_D0_K_PID_ALL_kaon", mean_D0_K_PID_ALL_kaon);
  mon_dst->setVariable("mean_D0_K_PID_SVD_kaon", mean_D0_K_PID_SVD_kaon);
  mon_dst->setVariable("mean_D0_K_PID_CDC_kaon", mean_D0_K_PID_CDC_kaon);
  mon_dst->setVariable("mean_D0_K_PID_TOP_kaon", mean_D0_K_PID_TOP_kaon);
  mon_dst->setVariable("mean_D0_K_PID_ARICH_kaon", mean_D0_K_PID_ARICH_kaon);
  mon_dst->setVariable("mean_D0_K_PID_ECL_kaon", mean_D0_K_PID_ECL_kaon);
  mon_dst->setVariable("mean_D0_K_PID_KLM_kaon", mean_D0_K_PID_KLM_kaon);
}

◆ ExtractHadron()

void ExtractHadron ( void )

private

Extract data from histograms related to Hadron.

Definition at line 1098 of file DQMHistAnalysisMiraBelle.cc.

{
  hadron_main->Clear();
  hadron_main->Divide(2, 2);
  // ========== hadronb2 + tight
  // get existing histograms produced by DQM modules
 
  if (auto* histhad_nECLClusters = findHist("PhysicsObjectsMiraBelleHadron/hist_nECLClusters"); histhad_nECLClusters != nullptr) {
    hadron_main->cd(1);
    histhad_nECLClusters->Draw();
  } else {
    B2WARNING("Can not find the histhad_nECLClusters histogram!");
  }
  if (auto* histhad_visibleEnergyCMSnorm = findHist("PhysicsObjectsMiraBelleHadron/hist_visibleEnergyCMSnorm");
      histhad_visibleEnergyCMSnorm != nullptr) {
    hadron_main->cd(2);
    histhad_visibleEnergyCMSnorm->Draw();
  } else {
    B2WARNING("Can not find the histhad_visibleEnergyCMSnorm histogram!");
  }
 
  if (auto* histhad_EsumCMSnorm = findHist("PhysicsObjectsMiraBelleHadron/hist_EsumCMSnorm"); histhad_EsumCMSnorm != nullptr) {
    hadron_main->cd(3);
    histhad_EsumCMSnorm->Draw();
  } else {
    B2WARNING("Can not find the histhad_EsumCMSnorm histogram!");
  }
  if (auto* histhad_R2 = findHist("PhysicsObjectsMiraBelleHadron/hist_R2"); histhad_R2 != nullptr) {
    hadron_main->cd(4);
    histhad_R2->Draw();
  } else {
    B2WARNING("Can not find the histhad_R2 histogram!");
  }
 
  auto* histhad_physicsresultsH = findHist("PhysicsObjectsMiraBelleHadron/hist_physicsresultsH");
  if (histhad_physicsresultsH == nullptr) {
    B2WARNING("Can not find the histhad_physicsresultsH histogram!");
    return;
  }
  auto* histbh_nECLClusters = findHist("PhysicsObjectsMiraBelleBhabha/hist_nECLClusters");
  if (histbh_nECLClusters == nullptr) {
    B2WARNING("Can not find the PhysicsObjectsMiraBelleBhabha/hist_nECLClusters histogram!");
    return;
  }
 
 
  // calculate the values of monitoring variables
  double had_ntot = histhad_physicsresultsH->GetBinContent(3);
  mon_bhabha->setVariable("had_ntot", had_ntot);
  double bh_neve_bhabha = histbh_nECLClusters->GetEntries(); // copy from bhabha
  if (had_ntot != 0 and bh_neve_bhabha != 0) {
    double ratio_hadron_bhabha = had_ntot / bh_neve_bhabha;
    //pull
    double error_ratio = ratio_hadron_bhabha * sqrt((1 / had_ntot) + (1 / bh_neve_bhabha));
    double ratio_pull_hadBhabha = (ratio_hadron_bhabha - m_reference_hadb) / error_ratio;
    mon_hadron->setVariable("ratio_hadron_bhabha", ratio_hadron_bhabha);
    mon_hadron->setVariable("error_ratio", error_ratio);
    mon_hadron->setVariable("ratio_pull_hadBhabha", ratio_pull_hadBhabha);
  }
}

◆ ExtractMumu()

void ExtractMumu ( void )

private

Extract data from histograms related to MuMu.

Definition at line 85 of file DQMHistAnalysisMiraBelle.cc.

{
  // ========== mumutight
  // get existing histograms produced by DQM modules
  auto* hist_npxd = findHist("PhysicsObjectsMiraBelle/hist_npxd");
  auto* hist_nsvd = findHist("PhysicsObjectsMiraBelle/hist_nsvd");
  auto* hist_ncdc = findHist("PhysicsObjectsMiraBelle/hist_ncdc");
  auto* hist_topdig = findHist("PhysicsObjectsMiraBelle/hist_topdig");
  auto* hist_DetPhotonARICH = findHist("PhysicsObjectsMiraBelle/hist_DetPhotonARICH");
  auto* hist_klmTotalHits = findHist("PhysicsObjectsMiraBelle/hist_klmTotalHits");
  auto* hist_klmClusterLayers = findHist("PhysicsObjectsMiraBelle/hist_klmClusterLayers");
  auto* hist_dD0 = findHist("PhysicsObjectsMiraBelle/hist_dD0");
  auto* hist_dZ0 = findHist("PhysicsObjectsMiraBelle/hist_dZ0");
  auto* hist_dPtcms = findHist("PhysicsObjectsMiraBelle/hist_dPtcms");
  auto* hist_Pval = findHist("PhysicsObjectsMiraBelle/hist_Pval");
  auto* hist_nExtraCDCHits = findHist("PhysicsObjectsMiraBelle/hist_nExtraCDCHits");
  auto* hist_nECLClusters = findHist("PhysicsObjectsMiraBelle/hist_nECLClusters");
  auto* hist_muid = findHist("PhysicsObjectsMiraBelle/hist_muid");
  auto* hist_inv_p = findHist("PhysicsObjectsMiraBelle/hist_inv_p");
  auto* hist_ndf = findHist("PhysicsObjectsMiraBelle/hist_ndf");
  auto* hist_D0 = findHist("PhysicsObjectsMiraBelle/hist_D0");
  auto* hist_Z0 = findHist("PhysicsObjectsMiraBelle/hist_Z0");
  auto* hist_theta = findHist("PhysicsObjectsMiraBelle/hist_theta");
  auto* hist_Phi0 = findHist("PhysicsObjectsMiraBelle/hist_Phi0");
  auto* hist_Pt = findHist("PhysicsObjectsMiraBelle/hist_Pt");
  auto* hist_Mom = findHist("PhysicsObjectsMiraBelle/hist_Mom");
  auto* hist_klmTotalBarrelHits = findHist("PhysicsObjectsMiraBelle/hist_klmTotalBarrelHits");
  auto* hist_klmTotalEndcapHits = findHist("PhysicsObjectsMiraBelle/hist_klmTotalEndcapHits");
  auto* hist_dPhicms = findHist("PhysicsObjectsMiraBelle/hist_dPhicms");
  auto* hist_hltEff = findHist("PhysicsObjectsMiraBelle/hist_hltEff");
 
  if (hist_npxd == nullptr) {
    B2ERROR("Can not find the hist_npxd histogram!");
    return;
  }
  if (hist_nsvd == nullptr) {
    B2ERROR("Can not find the hist_nsvd histogram!");
    return;
  }
  if (hist_ncdc == nullptr) {
    B2ERROR("Can not find the hist_ncdc histogram!");
    return;
  }
  if (hist_topdig == nullptr) {
    B2ERROR("Can not find the hist_topdig histogram!");
    return;
  }
  if (hist_DetPhotonARICH == nullptr) {
    B2ERROR("Can not find the hist_DetPhotonARICH histogram!");
    return;
  }
  if (hist_klmTotalHits == nullptr) {
    B2ERROR("Can not find the hist_klmTotalHits histogram!");
    return;
  }
  if (hist_klmClusterLayers == nullptr) {
    B2ERROR("Can not find the hist_klmClusterLayers histogram!");
    return;
  }
  if (hist_dD0 == nullptr) {
    B2ERROR("Can not find the hist_dD0 histogram!");
    return;
  }
  if (hist_dZ0 == nullptr) {
    B2ERROR("Can not find the hist_dZ0 histogram!");
    return;
  }
  if (hist_dPtcms == nullptr) {
    B2ERROR("Can not find the hist_dPtcms histogram!");
    return;
  }
  if (hist_Pval == nullptr) {
    B2ERROR("Can not find the hist_Pval histogram!");
    return;
  }
  if (hist_nExtraCDCHits == nullptr) {
    B2ERROR("Can not find the hist_nExtraCDCHits histogram!");
    return;
  }
  if (hist_nECLClusters == nullptr) {
    B2ERROR("Can not find the hist_nECLClusters histogram!");
    return;
  }
  if (hist_muid == nullptr) {
    B2ERROR("Can not find the hist_muid histogram!");
    return;
  }
  if (hist_inv_p == nullptr) {
    B2ERROR("Can not find the hist_inv_p histogram!");
    return;
  }
  if (hist_ndf == nullptr) {
    B2ERROR("Can not find the hist_ndf histogram!");
    return;
  }
  if (hist_D0 == nullptr) {
    B2ERROR("Can not find the hist_D0 histogram!");
    return;
  }
  if (hist_Z0 == nullptr) {
    B2ERROR("Can not find the hist_Z0 histogram!");
    return;
  }
  if (hist_theta == nullptr) {
    B2ERROR("Can not find the hist_theta histogram!");
    return;
  }
  if (hist_Phi0 == nullptr) {
    B2ERROR("Can not find the hist_Phi0 histogram!");
    return;
  }
  if (hist_Pt == nullptr) {
    B2ERROR("Can not find the hist_Pt histogram!");
    return;
  }
  if (hist_Mom == nullptr) {
    B2ERROR("Can not find the hist_Mom histogram!");
    return;
  }
  if (hist_klmTotalBarrelHits == nullptr) {
    B2ERROR("Can not find the hist_klmTotalBarrelHits histogram!");
    return;
  }
  if (hist_klmTotalEndcapHits == nullptr) {
    B2ERROR("Can not find the hist_klmTotalEndcapHits histogram!");
    return;
  }
  if (hist_dPhicms == nullptr) {
    B2ERROR("Can not find the hist_dPhicms histogram!");
    return;
  }
  if (hist_hltEff == nullptr) {
    B2ERROR("Can not find the hist_hltEff histogram!");
    return;
  }
 
 
  // Make TCanvases
  // --- Mumu_Main
  mumu_main->Divide(4, 3);
  mumu_main->cd(1);  hist_theta->Draw();
  mumu_main->cd(2);  hist_Phi0->Draw();
  mumu_main->cd(3);  hist_Mom->Draw();
  mumu_main->cd(4);  hist_Pt->Draw();
  mumu_main->cd(5);  hist_npxd->Draw();
  mumu_main->cd(6);  hist_nsvd->Draw();
  mumu_main->cd(7);  hist_ncdc->Draw();
  mumu_main->cd(8);  hist_topdig->Draw();
  mumu_main->cd(9);  hist_DetPhotonARICH->Draw();
  mumu_main->cd(10);  hist_klmClusterLayers->Draw();
  mumu_main->cd(11);  hist_nExtraCDCHits->Draw();
  mumu_main->cd(12);  hist_nECLClusters->Draw();
  // --- Mumu_Resolution
  mumu_resolution->Divide(3, 3);
  mumu_resolution->cd(1);  hist_inv_p->Draw();
  mumu_resolution->cd(2);  hist_dD0->Draw();
  mumu_resolution->cd(3);  hist_dZ0->Draw();
  mumu_resolution->cd(4);  hist_dPtcms->Draw();
  mumu_resolution->cd(5);  hist_dPhicms->Draw();
  mumu_resolution->cd(6);  hist_nECLClusters->Draw();
  mumu_resolution->cd(7);  hist_nExtraCDCHits->Draw();
  mumu_resolution->cd(8);  hist_ndf->Draw();
  mumu_muon_val->Divide(2, 2);
  // --- Muon variables
  mumu_muon_val->cd(1);  hist_muid->Draw();
  mumu_muon_val->cd(2);  hist_Pval->Draw();
  mumu_muon_val->cd(3);  hist_theta->Draw();
  mumu_muon_val->cd(4);  hist_Phi0->Draw();
 
  // calculate the values of monitoring variables
  double mean_npxd = hist_npxd->GetMean();
  double mean_nsvd = hist_nsvd->GetMean();
  double mean_ncdc = hist_ncdc->GetMean();
  double mean_nklmlayer = hist_klmClusterLayers->GetMean();
  double mean_nklm = hist_klmTotalHits->GetMean();
  double mean_nbklm = hist_klmTotalBarrelHits->GetMean();
  double mean_neklm = hist_klmTotalEndcapHits->GetMean();
  double mean_topdig = hist_topdig->GetMean();
  double mean_parich = hist_DetPhotonARICH->GetMean();
  double mean_ncdc_ex = hist_nExtraCDCHits->GetMean();
  double mean_necl = hist_nECLClusters->GetMean();
  double mean_muid = hist_muid->GetMean();
  double mean_d0 = hist_D0->GetMean();
  double mean_z0 = hist_Z0->GetMean();
  double mean_pval = hist_Pval->GetMean();
  double mean_ndf = hist_ndf->GetMean();
  double dif_ndf_ncdc = mean_ndf - mean_ncdc;
  double mean_dd0 = hist_dD0->GetMean();
  double mean_dz0 = hist_dZ0->GetMean();
  double mean_dpt = hist_dPtcms->GetMean();
  double rms_dd0 = hist_dD0->GetRMS();
  double rms_dz0 = hist_dZ0->GetRMS();
  double rms_dpt = hist_dPtcms->GetRMS();
  double sigma68_dd0 = getSigma68(hist_dD0);
  double sigma68_dz0 = getSigma68(hist_dZ0);
  double sigma68_dpt = getSigma68(hist_dPtcms);
  int ntot = hist_nsvd->GetEntries();
  double neve_mumu = ntot;
  double goodmu_frac = -1.;
  double pval_frac_0 = -1.;
  double pval_frac_1 = -1.;
  double nocdc_frac = -1.;
  double notop_frac = -1.;
  double noarich_frac = -1.;
  double muidcontent = hist_muid->GetBinContent(1);
  double goodmu_o_badmu = -1.;
  if (muidcontent != 0) {
    goodmu_o_badmu = hist_muid->GetBinContent(20) / muidcontent;
  } else {
    goodmu_o_badmu = 0.0;
  }
  double pval_more95 = 0.0;
  double pval_less05 = 0.0;
  for (int i = 95; i < 100; i++) pval_more95 += hist_Pval->GetBinContent(i + 1);
  for (int i = 0; i < 5; i++) pval_less05 += hist_Pval->GetBinContent(i + 1);
  if (neve_mumu != 0) {
    goodmu_frac = hist_muid->GetBinContent(20) / neve_mumu;
    pval_frac_0 = pval_less05 / neve_mumu;
    pval_frac_1 = pval_more95 / neve_mumu;
    nocdc_frac = hist_ncdc->GetBinContent(1) / neve_mumu;
    notop_frac = hist_topdig->GetBinContent(1) / neve_mumu;
    noarich_frac = hist_DetPhotonARICH->GetBinContent(1) / neve_mumu;
  } else {
    goodmu_frac = 0.0;
    pval_frac_0 = 0.0;
    pval_frac_1 = 0.0;
    nocdc_frac = 0.0;
    notop_frac = 0.0;
    noarich_frac = 0.0;
  }
  //Calculate M(mumu)
  double peak_mumu = hist_inv_p->GetXaxis()->GetBinCenter(hist_inv_p->GetMaximumBin());
  TF1* f_mumuInvM = new TF1("f_mumuInvM", "gaus", peak_mumu - 0.05, peak_mumu + 0.05);
  f_mumuInvM->SetParameters(hist_inv_p->GetMaximum(), peak_mumu, 0.045);
  f_mumuInvM->SetParLimits(1, peak_mumu - 0.05, peak_mumu + 0.05);
  f_mumuInvM->SetParLimits(2, 0.01, 0.08);
  hist_inv_p->Fit(f_mumuInvM, "R");
  double fit_mumumass = f_mumuInvM->GetParameter(1);
  if (fit_mumumass < 9.) fit_mumumass = 9.;
  if (fit_mumumass > 12.) fit_mumumass = 12.;
  double fit_mumumass_error = f_mumuInvM->GetParError(1);
  double pull_mumumass = (fit_mumumass - m_reference_mumu) / fit_mumumass_error;
  double fit_sigma_mumu = f_mumuInvM->GetParameter(2);
 
  // Variables needed for selectmumu efficiency in HLT
  double Nmumu_ECLMuonPair = hist_hltEff->GetBinContent(2);
  double Nmumu_ECLMuonPairSelectmumu = hist_hltEff->GetBinContent(3);
 
  // set values
  mon_mumu->setVariable("mean_npxd", mean_npxd);
  mon_mumu->setVariable("mean_nsvd", mean_nsvd);
  mon_mumu->setVariable("mean_ncdc", mean_ncdc);
  mon_mumu->setVariable("mean_nklmlayer", mean_nklmlayer);
  mon_mumu->setVariable("mean_nklm", mean_nklm);
  mon_mumu->setVariable("mean_nbklm", mean_nbklm);
  mon_mumu->setVariable("mean_neklm", mean_neklm);
  mon_mumu->setVariable("mean_topdig", mean_topdig);
  mon_mumu->setVariable("mean_parich", mean_parich);
  mon_mumu->setVariable("mean_ncdc_ex", mean_ncdc_ex);
  mon_mumu->setVariable("mean_necl", mean_necl);
  mon_mumu->setVariable("mean_muid", mean_muid);
  mon_mumu->setVariable("mean_d0", mean_d0);
  mon_mumu->setVariable("mean_z0", mean_z0);
  mon_mumu->setVariable("mean_pval", mean_pval);
  mon_mumu->setVariable("mean_ndf", mean_ndf);
  mon_mumu->setVariable("dif_ndf_ncdc", dif_ndf_ncdc);
  mon_mumu->setVariable("mean_dd0", mean_dd0);
  mon_mumu->setVariable("mean_dz0", mean_dz0);
  mon_mumu->setVariable("mean_dpt", mean_dpt);
  mon_mumu->setVariable("rms_dd0", rms_dd0);
  mon_mumu->setVariable("rms_dz0", rms_dz0);
  mon_mumu->setVariable("rms_dpt", rms_dpt);
  mon_mumu->setVariable("sigma68_dd0", sigma68_dd0);
  mon_mumu->setVariable("sigma68_dz0", sigma68_dz0);
  mon_mumu->setVariable("sigma68_dpt", sigma68_dpt);
  mon_mumu->setVariable("neve_mumu", neve_mumu);
  mon_mumu->setVariable("goodmu_frac", goodmu_frac);
  mon_mumu->setVariable("goodmu_o_badmu", goodmu_o_badmu);
  mon_mumu->setVariable("pval_frac_0", pval_frac_0);
  mon_mumu->setVariable("pval_frac_1", pval_frac_1);
  mon_mumu->setVariable("nocdc_frac", nocdc_frac);
  mon_mumu->setVariable("notop_frac", notop_frac);
  mon_mumu->setVariable("noarich_frac", noarich_frac);
  mon_mumu->setVariable("fit_mumumass", fit_mumumass);
  mon_mumu->setVariable("fit_mumumass_error", fit_mumumass_error);
  mon_mumu->setVariable("pull_mumumass", pull_mumumass);
  mon_mumu->setVariable("sigma_mumumass", fit_sigma_mumu);
  mon_mumu->setVariable("Nmumu_ECLMuonPair", Nmumu_ECLMuonPair);
  mon_mumu->setVariable("Nmumu_ECLMuonPairSelectmumu", Nmumu_ECLMuonPairSelectmumu);
}

◆ ExtractNEvent()

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

staticinherited

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

Definition at line 362 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 350 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");
}

◆ ExtractTaupair()

void ExtractTaupair ( void )

private

Extract data from histograms related to TauPairs.

Definition at line 817 of file DQMHistAnalysisMiraBelle.cc.

{
  //--- L1 efficiency with taupair
  if (auto* ptr = findHist("PhysicsObjectsMiraBelleTau/hist_L1ECL1x1"); ptr != nullptr) {
    for (int bin = 1; bin <= ptr->GetXaxis()->GetNbins(); bin++) {
      if (!std::string(ptr->GetXaxis()->GetBinLabel(bin)).empty()) {
        std::string label = std::string("ECL_1_1_") + ptr->GetXaxis()->GetBinLabel(bin);
        mon_tautau->setVariable(label, ptr->GetBinContent(bin));
      }
    }
  }
  if (auto* ptr = findHist("PhysicsObjectsMiraBelleTau/hist_L1ECL1x3"); ptr != nullptr) {
    for (int bin = 1; bin <= ptr->GetXaxis()->GetNbins(); bin++) {
      if (!std::string(ptr->GetXaxis()->GetBinLabel(bin)).empty()) {
        std::string label = std::string("ECL_1_3_") + ptr->GetXaxis()->GetBinLabel(bin);
        mon_tautau->setVariable(label, ptr->GetBinContent(bin));
      }
    }
  }
  if (auto* ptr = findHist("PhysicsObjectsMiraBelleTau/hist_L1CDC1x1"); ptr != nullptr) {
    for (int bin = 1; bin <= ptr->GetXaxis()->GetNbins(); bin++) {
      if (!std::string(ptr->GetXaxis()->GetBinLabel(bin)).empty()) {
        std::string label = std::string("CDC_1_1_") + ptr->GetXaxis()->GetBinLabel(bin);
        mon_tautau->setVariable(label, ptr->GetBinContent(bin));
      }
    }
  }
  if (auto* ptr = findHist("PhysicsObjectsMiraBelleTau/hist_L1CDC1x3"); ptr != nullptr) {
    for (int bin = 1; bin <= ptr->GetXaxis()->GetNbins(); bin++) {
      if (!std::string(ptr->GetXaxis()->GetBinLabel(bin)).empty()) {
        std::string label = std::string("CDC_1_3_") + ptr->GetXaxis()->GetBinLabel(bin);
        mon_tautau->setVariable(label, ptr->GetBinContent(bin));
      }
    }
  }
  if (auto* ptr = findHist("PhysicsObjectsMiraBelleTau/hist_L1CDCKLM1x1"); ptr != nullptr) {
    for (int bin = 1; bin <= ptr->GetXaxis()->GetNbins(); bin++) {
      if (!std::string(ptr->GetXaxis()->GetBinLabel(bin)).empty()) {
        std::string label = std::string("CDCKLM_1_1_") + ptr->GetXaxis()->GetBinLabel(bin);
        mon_tautau->setVariable(label, ptr->GetBinContent(bin));
      }
    }
  }
  if (auto* ptr = findHist("PhysicsObjectsMiraBelleTau/hist_L1CDCKLM1x3"); ptr != nullptr) {
    for (int bin = 1; bin <= ptr->GetXaxis()->GetNbins(); bin++) {
      if (!std::string(ptr->GetXaxis()->GetBinLabel(bin)).empty()) {
        std::string label = std::string("CDCKLM_1_3_") + ptr->GetXaxis()->GetBinLabel(bin);
        mon_tautau->setVariable(label, ptr->GetBinContent(bin));
      }
    }
  }
}

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

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

◆ findHist()

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

staticinherited

Find histogram.

Parameters

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

Returns: The found histogram, or nullptr if not found.

Definition at line 154 of file DQMHistAnalysis.cc.

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

◆ findHistInCanvas()

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

inherited

Find histogram in corresponding canvas.

Parameters

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

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

Definition at line 214 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 249 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 270 of file DQMHistAnalysis.cc.

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

◆ findRefHist()

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

staticinherited

Find reference histogram.

Parameters

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

Returns: The found histogram, or nullptr if not found.

Definition at line 198 of file DQMHistAnalysis.cc.

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

◆ getAfterConditionPath()

Module::EAfterConditionPath getAfterConditionPath ( ) const

inherited

What to do after the conditional path is finished.

(defaults to c_End if no condition is set)

Definition at line 133 of file Module.cc.

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

◆ getAllConditionPaths()

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

inherited

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

Definition at line 150 of file Module.cc.

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

◆ getAllConditions()

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

inlineinherited

Return all set conditions for this module.

Definition at line 323 of file Module.h.

    {
      return m_conditions;
    }

◆ getCanvasUpdatedList()

static const CanvasUpdatedList & getCanvasUpdatedList ( )

inlinestaticinherited

Get the list of the canvas update status.

Returns: The list of the canvases.

Definition at line 202 of file DQMHistAnalysis.h.

202{ return s_canvasUpdatedList;};

◆ getCondition()

const ModuleCondition * getCondition ( ) const

inlineinherited

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

Definition at line 313 of file Module.h.

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

◆ getConditionPath()

std::shared_ptr< Path > getConditionPath ( ) const

inherited

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

Definition at line 113 of file Module.cc.

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

◆ getDelta()

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

inherited

Get Delta histogram.

Parameters

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

Returns: delta histogram or nullptr

Definition at line 118 of file DQMHistAnalysis.cc.

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

◆ getDeltaList()

static const DeltaList & getDeltaList ( )

inlinestaticinherited

Get the list of the delta histograms.

Returns: The list of the delta histograms.

Definition at line 196 of file DQMHistAnalysis.h.

196{ return s_deltaList;};

◆ getDescription()

const std::string & getDescription ( ) const

inlineinherited

Returns the description of the module.

Definition at line 201 of file Module.h.

201{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 493 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 516 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 589 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 603 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 539 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 564 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.

221{ 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.

184{ return s_histList;};

◆ getLogConfig()

LogConfig & getLogConfig ( )

inlineinherited

Returns the log system configuration.

Definition at line 224 of file Module.h.

224{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.

505{ 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 132 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.

190{ 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.

186{return m_name;}

◆ getPackage()

const std::string & getPackage ( ) const

inlineinherited

Returns the package this module is in.

Definition at line 196 of file Module.h.

196{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.

362{ 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 584 of file DQMHistAnalysis.h.

584{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.

380{ 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.

215{ 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 279 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 709 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 567 of file DQMHistAnalysis.h.

567{return m_useEpics;};

◆ getUseEpicsReadOnly()

static bool getUseEpicsReadOnly ( void )

inlinestaticinherited

Getter EPICS flag in read only mode.

Returns: flag if read only

Definition at line 573 of file DQMHistAnalysis.h.

573{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.

310{ 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.

377{ 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 309 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 )

finaloverrideprivatevirtual

Initialize the Module.

This method is called at the beginning of data processing.

Reimplemented from Module.

Definition at line 35 of file DQMHistAnalysisMiraBelle.cc.

{
  gROOT->cd();
 
  //mon_mumu = new Belle2::MonitoringObject("mumu");
  mon_mumu = getMonitoringObject("mumu");
  mon_dst = getMonitoringObject("dst");
  mon_tautau = getMonitoringObject("tautau");
  //bhabha,hadrons
  mon_bhabha = getMonitoringObject("bhabha");
  mon_hadron = getMonitoringObject("hadronb2");
 
  // make cavases to be added to MonitoringObject
  mumu_main = new TCanvas("mumu_main", "mumu_main", 0, 0, 800, 600);
  mumu_resolution = new TCanvas("mumu_resolution", "mumu_resolution", 0, 0, 800, 600);
  mumu_muon_val = new TCanvas("mumu_muon_val", "mumu_muon_val", 0, 0, 400, 400);
  dst_mass = new TCanvas("dst_mass", "dst_mass", 0, 0, 1200, 400);
  dst_pi_val = new TCanvas("dst_pi_val", "dst_pi_val", 0, 0, 800, 400);
  dst_k_val = new TCanvas("dst_k_val", "dst_k_val", 0, 0, 800, 400);
  //bhabha,hadrons
  bhabha_main = new TCanvas("bhabha_main", "bhabha_main", 0, 0, 800, 600);
  bhabha_resolution = new TCanvas("bhabha_resolution", "bhabha_resolution", 0, 0, 800, 600);
  hadron_main = new TCanvas("hadron_main", "hadron_main", 0, 0, 800, 600);
 
  // add canvases to MonitoringObject
  mon_mumu->addCanvas(mumu_main);
  mon_mumu->addCanvas(mumu_resolution);
  mon_mumu->addCanvas(mumu_muon_val);
  mon_dst->addCanvas(dst_mass);
  mon_dst->addCanvas(dst_pi_val);
  mon_dst->addCanvas(dst_k_val);
  //bhabha,hadrons
  mon_bhabha->addCanvas(bhabha_main);
  mon_bhabha->addCanvas(bhabha_resolution);
  mon_bhabha->addCanvas(hadron_main);
 
  B2DEBUG(20, "DQMHistAnalysisMiraBelle: initialized.");
 
}

◆ 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 691 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 759 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 374 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 384 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 379 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 656 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 644 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 650 of file DQMHistAnalysis.cc.

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

◆ resetDeltaList()

void resetDeltaList ( void )

staticinherited

Reset Delta.

Definition at line 333 of file DQMHistAnalysis.cc.

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

◆ scaleReference()

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

staticinherited

Using the original and reference, create scaled version.

Parameters

scaling	scaling algorithm
hist	pointer to histogram
ref	pointer to reference

Returns: scaled reference

Definition at line 172 of file DQMHistAnalysis.cc.

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

◆ setAbortLevel()

void setAbortLevel ( int abortLevel )

inherited

Configure the abort log level.

Definition at line 67 of file Module.cc.

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

◆ setDebugLevel()

void setDebugLevel ( int debugLevel )

inherited

Configure the debug messaging level.

Definition at line 61 of file Module.cc.

{
  m_logConfig.setDebugLevel(debugLevel);
}

◆ setDescription()

void setDescription ( const std::string & description )

protectedinherited

Sets the description of the module.

Parameters

description A description of the module.

Definition at line 214 of file Module.cc.

{
  m_description = description;
}

◆ setEpicsPV() [1/4]

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

inherited

Write value to a EPICS PV.

Parameters

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

Definition at line 412 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 424 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 454 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 466 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 436 of file DQMHistAnalysis.cc.

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

◆ setEpicsStringPV() [2/2]

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

inherited

Write string to a EPICS PV.

Parameters

index	index of PV
value	string to write

Definition at line 478 of file DQMHistAnalysis.cc.

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

◆ setEventProcessed()

static void setEventProcessed ( int e )

inlinestaticinherited

Set the number of processed events.

(Attention, asynch histogram updates!)

e Processed events.

Definition at line 233 of file DQMHistAnalysis.h.

233{s_eventProcessed = e;};

◆ setLogConfig()

void setLogConfig ( const LogConfig & logConfig )

inlineinherited

Set the log system configuration.

Definition at line 229 of file Module.h.

229{m_logConfig = logConfig;}

◆ setLogInfo()

void setLogInfo	(	int	logLevel,
		unsigned int	logInfo )

inherited

Configure the printed log information for the given level.

Parameters

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

Definition at line 73 of file Module.cc.

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

◆ setLogLevel()

void setLogLevel ( int logLevel )

inherited

Configure the log level.

Definition at line 55 of file Module.cc.

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

◆ setName()

void setName ( const std::string & name )

inlineinherited

Set the name of the module.

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

Parameters

name	The name of the module

Definition at line 213 of file Module.h.

213{ m_name = name; };

◆ setParamList()

void setParamList ( const ModuleParamList & params )

inlineprotectedinherited

Replace existing parameter list.

Definition at line 500 of file Module.h.

500{ 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 590 of file DQMHistAnalysis.h.

590{ 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.

227{s_runType = t;};

◆ setType()

void setType ( const std::string & type )

protectedinherited

Set the module type.

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

Definition at line 48 of file Module.cc.

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

◆ setUseEpics()

static void setUseEpics ( bool flag )

inlinestaticinherited

Setter for EPICS usage.

Parameters

flag	set in use

Definition at line 555 of file DQMHistAnalysis.h.

555{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 561 of file DQMHistAnalysis.h.

561{m_epicsReadOnly = flag;};

◆ StringSplit()

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

staticinherited

Helper function for string token split.

Parameters

s	String to split
delim	delimiter

Returns: vector of strings

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

finaloverrideprivatevirtual

Termination action.

Clean-up, close files, summarize statistics, etc.

Reimplemented from Module.

Definition at line 1170 of file DQMHistAnalysisMiraBelle.cc.

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

◆ 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 340 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 345 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 617 of file DQMHistAnalysis.cc.

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

Member Data Documentation

◆ bhabha_main

TCanvas* bhabha_main = nullptr

protected

TCanvas object for main display.

Definition at line 110 of file DQMHistAnalysisMiraBelle.h.

◆ bhabha_resolution

TCanvas* bhabha_resolution = nullptr

protected

TCanvas object for resolution display.

Definition at line 111 of file DQMHistAnalysisMiraBelle.h.

◆ dst_k_val

TCanvas* dst_k_val = nullptr

protected

TCanvas object for PID distributions of K.

Definition at line 108 of file DQMHistAnalysisMiraBelle.h.

◆ dst_mass

TCanvas* dst_mass = nullptr

protected

TCanvas object for D0 mass, delta mass, and pi0 mass.

Definition at line 106 of file DQMHistAnalysisMiraBelle.h.

◆ dst_pi_val

TCanvas* dst_pi_val = nullptr

protected

TCanvas object for PID distributions of pi.

Definition at line 107 of file DQMHistAnalysisMiraBelle.h.

◆ hadron_main

TCanvas* hadron_main = nullptr

protected

TCanvas object for main display.

Definition at line 113 of file DQMHistAnalysisMiraBelle.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_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_hasReturnValue

bool m_hasReturnValue

privateinherited

True, if the return value is set.

Definition at line 517 of file Module.h.

◆ m_logConfig

LogConfig m_logConfig

privateinherited

The log system configuration of the module.

Definition at line 513 of file Module.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_name

std::string m_name

privateinherited

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

Definition at line 507 of file Module.h.

◆ m_package

std::string m_package

privateinherited

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

Definition at line 509 of file Module.h.

◆ m_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_reference_hadb

double m_reference_hadb

protected

Reference value for hadronb2_tight bhabha ratio.

Definition at line 93 of file DQMHistAnalysisMiraBelle.h.

◆ m_reference_mumu

double m_reference_mumu

protected

Reference value for mumu inv mass.

Definition at line 94 of file DQMHistAnalysisMiraBelle.h.

◆ m_returnValue

int m_returnValue

privateinherited

The return value.

Definition at line 518 of file Module.h.

◆ m_scale_dst

double m_scale_dst

protected

Scale factor "signal region" / "sideband", assuming uniform events.

Definition at line 92 of file DQMHistAnalysisMiraBelle.h.

◆ m_type

std::string m_type

privateinherited

The type of the module, saved as a string.

Definition at line 508 of file Module.h.

◆ m_useEpics

bool m_useEpics = false

staticprivateinherited

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

Definition at line 149 of file DQMHistAnalysis.h.

◆ mon_bhabha

Belle2::MonitoringObject* mon_bhabha = nullptr

protected

MonitoringObject for bhabha events.

Definition at line 98 of file DQMHistAnalysisMiraBelle.h.

◆ mon_dst

Belle2::MonitoringObject* mon_dst = nullptr

protected

MonitoringObject for D* events.

Definition at line 96 of file DQMHistAnalysisMiraBelle.h.

◆ mon_hadron

Belle2::MonitoringObject* mon_hadron = nullptr

protected

MonitoringObject for hadronb2 tight events.

Definition at line 99 of file DQMHistAnalysisMiraBelle.h.

◆ mon_mumu

Belle2::MonitoringObject* mon_mumu = nullptr

protected

MonitoringObject for mumu events.

Definition at line 95 of file DQMHistAnalysisMiraBelle.h.

◆ mon_tautau

Belle2::MonitoringObject* mon_tautau = nullptr

protected

MonitoringObject for tau pair events.

Definition at line 97 of file DQMHistAnalysisMiraBelle.h.

◆ mumu_main

TCanvas* mumu_main = nullptr

protected

TCanvas object for main display.

Definition at line 102 of file DQMHistAnalysisMiraBelle.h.

◆ mumu_muon_val

TCanvas* mumu_muon_val = nullptr

protected

TCanvas object for mumu values display.

Definition at line 104 of file DQMHistAnalysisMiraBelle.h.

◆ mumu_resolution

TCanvas* mumu_resolution = nullptr

protected

TCanvas object for resolution display.

Definition at line 103 of file DQMHistAnalysisMiraBelle.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/DQMHistAnalysisMiraBelle.h
dqm/analysis/modules/src/DQMHistAnalysisMiraBelle.cc

Public Types

Public Member Functions

Static Public Member Functions

Static Public Attributes

Protected Member Functions

Protected Attributes

Private Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Member Typedef Documentation

◆ CanvasUpdatedList

◆ DeltaList

◆ EAfterConditionPath

◆ HistList

◆ MonObjList

◆ RefList

Member Enumeration Documentation

◆ EModulePropFlags

◆ ERefScaling

◆ EStatus

◆ EStatusColor

Constructor & Destructor Documentation

◆ DQMHistAnalysisMiraBelleModule()

Member Function Documentation

◆ addDeltaPar()

◆ addHist()

◆ addRefHist()

◆ beginRun()

◆ CheckEpicsError() [1/2]

◆ CheckEpicsError() [2/2]

◆ checkPVStatus()

◆ cleanupEpicsPVs()

◆ clearCanvases()

◆ clearHistList()

◆ clearlist()

◆ clearRefList()

◆ clone()

◆ colorizeCanvas()

◆ def_beginRun()

◆ def_endRun()

◆ def_event()

◆ def_initialize()

◆ def_terminate()

◆ endRun()

◆ evalCondition()

◆ event()

◆ exposePythonAPI()

◆ ExtractBhabha()

◆ ExtractD0Star()

◆ ExtractHadron()

◆ ExtractMumu()

◆ ExtractNEvent()

◆ ExtractRunType()

◆ ExtractTaupair()

◆ findCanvas()

◆ findHist()

◆ findHistInCanvas()

◆ findHistInFile()

◆ findMonitoringObject()

◆ findRefHist()

◆ getAfterConditionPath()

◆ getAllConditionPaths()

◆ getAllConditions()

◆ getCanvasUpdatedList()

◆ getCondition()

◆ getConditionPath()

◆ getDelta()

◆ getDeltaList()

◆ getDescription()

◆ getEpicsPV() [1/2]

◆ getEpicsPV() [2/2]

◆ getEpicsPVChID() [1/2]

◆ getEpicsPVChID() [2/2]

◆ getEpicsStringPV() [1/2]

◆ getEpicsStringPV() [2/2]

◆ getEventProcessed()

◆ getFileNames()

◆ getHistList()

◆ getLogConfig()