top

Topics
	top data objects
	top modules

Classes
class	TOPDatabaseImporter
	TOP database importer. More...
class	TOPASICChannel
	Calibration constants of a singe ASIC channel: pedestals, gains and time axis. More...
class	TOPASICGains
	Calibration constants of a single ASIC window: gains. More...
class	TOPASICPedestals
	Calibration constants of a single ASIC window: pedestals. More...
class	TOPCalAsicShift
	Calibration constants for ASIC shifts of all 16 modules. More...
class	TOPCalChannelMask
	Channel status for all 512 channels of 16 modules. More...
class	TOPCalChannelNoise
	r.m.s. More...
class	TOPCalChannelPulseHeight
	Pulse height parameterizations for all 512 channels of 16 modules. More...
class	TOPCalChannelRQE
	Class to store relative quantum efficiency of channels w.r.t initial one measured in PMT QA QE is expected to decrease during the experiment due to aging induced by acuumulated charge. More...
class	TOPCalChannelT0
	Channel T0 calibration constants for all 512 channels of 16 modules. More...
class	TOPCalChannelThreshold
	Value of the threshold (in ADC counts) used for the pulse identification, for all 512 channels of 16 modules. More...
class	TOPCalChannelThresholdEff
	Class to store the threshold efficiency (i.e. More...
class	TOPCalCommonT0
	Common T0 calibration constant. More...
class	TOPCalEventT0Offset
	Class to store the calibrated EventT0 offsets of other detector components The offsets are measured relative to TOP EventT0 (e.g CDC - TOP etc) More...
class	TOPCalFillPatternOffset
	Class to store the offset of reconstructed fill pattern. More...
class	TOPCalIntegratedCharge
	Class to store integrated charge per channel. More...
class	TOPCalModuleAlignment
	Alignment constants for all 16 modules. More...
class	TOPCalModuleT0
	Module T0 calibration constants for all 16 modules. More...
class	TOPCalPhotonYields
	Class to store photon pixel yields for PMT ageing studies, and equalized alpha ratios for finding optically decoupled PMT's. More...
class	TOPCalTimebase
	Sample time calibration constants for all channels. More...
class	TOPCalTimeWalk
	Calibration constants for time-walk correction and for tuning of electronic time resolution in digitization. More...
class	TOPChannelMap
	Map of pixels and channels within the carrier board. More...
class	TOPFrontEndMap
	Mapping of a boardstack number within a module to SCROD and COPPER/Finesse. More...
class	TOPFrontEndSetting
	Front-end settings: storage depths, lookback, readout windows etc. More...
class	TOPGeoBarSegment
	Geometry parameters of a quartz bar segment. More...
class	TOPGeoBase
	Base class for geometry parameters. More...
class	TOPGeoColdPlate
	Geometry parameters of cold plate (simplified) More...
class	TOPGeoEndPlate
	Geometry parameters of forward end plate (simplified) More...
class	TOPGeoFrontEnd
	Geometry parameters of board stack (front-end electronic module) More...
class	TOPGeoHoneycombPanel
	Geometry parameters of honeycomb panel. More...
class	TOPGeometry
	Geometry parameters of TOP. More...
class	TOPGeoMirrorSegment
	Geometry parameters of a mirror segment. More...
class	TOPGeoModule
	Geometry parameters of a module (optical components + positioning) More...
class	TOPGeoModuleDisplacement
	Displacement parameters of a TOP module. More...
class	TOPGeoPMT
	Geometry parameters of MCP-PMT. More...
class	TOPGeoPMTArray
	Geometry parameters of MCP-PMT array. More...
class	TOPGeoPMTArrayDisplacement
	Displacement parameters of MCP-PMT array. More...
class	TOPGeoPrism
	Geometry parameters of prism. More...
class	TOPGeoPrismEnclosure
	Geometry parameters of prism enclosure (simplified) More...
class	TOPGeoQBB
	Geometry parameters of Quartz Bar Box (mother class) More...
class	TOPGeoSideRails
	Geometry parameters of side rails (simplified) More...
class	TOPNominalQE
	Nominal quantum efficiency of PMT. More...
class	TOPNominalTDC
	Nominal time-to-digit conversion parameters (simplified model) More...
class	TOPNominalTTS
	Nominal time transition spread of PMT. More...
class	TOPPmtGainPar
	Parameterized PMT gain (from laser testing at Nagoya) More...
class	TOPPmtInstallation
	PMT position within a module. More...
class	TOPPmtObsoleteData
	PMT specs from Hamamatsu. More...
class	TOPPmtQE
	Nagoya measurements of quantum efficiency. More...
class	TOPPmtTTSHisto
	Raw measurements of TTS (histograms) for each PMT pixel. More...
class	TOPPmtTTSPar
	Parameterized TTS for each PMT pixel. More...
struct	TOPPulseHeightPar
	Struct holding the pulse height parameterizations for the TOP counter. More...
class	TOPSampleTimes
	Calibration constants of a singe ASIC channel: time axis (sample times) More...
class	TOPSignalShape
	Normalized shape of single photon pulse (waveform) Pulse must be positive. More...
class	TOPWavelengthFilter
	Bulk transmittance of wavelength filter. More...

Typedefs
typedef map< unsigned, size_t >::const_iterator	Iterator
	Iteratior for m_map.
typedef std::pair< double, double >	Pair
	Shorthand for std::pair<double, double>

Functions
void	importSampleTimeCalibration (std::string fileNames, int firstExp=0, int firstRun=0, int lastExp=-1, int lastRun=-1)
	Import sample time calibration constants to database Output of TOPTimeBaseCalibrator (root files with histograms)
void	importLocalT0Calibration (std::string fileNames, int firstExp=0, int firstRun=0, int lastExp=-1, int lastRun=-1)
	Import channel-by-channel T0 calibration constants to database The input is the root file with ntuple produced by TOPLaserCalibrator.
void	importChannelT0 (std::string fileName, int expNo, int firstRun, int lastRun)
	Import channel T0 calibration constants The input is a root file with 1D histograms (one per slot, named "channelT0_slot*")
void	importAsicShifts_BS13d (double s0, double s1, double s2, double s3, int expNo, int firstRun, int lastRun)
	Import ASIC shifts of BS13d.
void	importOfflineCommonT0Calibration (std::string fileName, int firstExp=0, int firstRun=0, int lastExp=-1, int lastRun=-1)
	Import common T0 calibration constants derived form the offline data reprocessing to database The input is a root file containing a tree (one per run).
void	importCommonT0 (double value, double error, int expNo, int firstRun, int lastRun, bool roughlyCalibrated=false)
	Import common T0 calibration constants.
void	importModuleT0Calibration (std::string fileName, int firstExp=0, int firstRun=0, int lastExp=-1, int lastRun=-1)
	Import module T0 calibration constants to database The input is the text file.
void	importModuleT0 (std::string fileName, int expNo, int firstRun, int lastRun)
	Import module T0 calibration constants The input is a root file with 1D histogram (name is "moduleT0")
void	printSampleTimeCalibrationInfo ()
	Prints sample time calibration info about constants stored in database.
void	printSampleTimeCalibration ()
	Print sample time calibration constants stored in database.
void	importChannelMask (std::string fileName, int expNo, int firstRun, int lastRun)
	Import channel mask from a root file (given as 1D histograms, one per slot).
void	generateFakeChannelMask (double fractionDead, double fractionHot, int firstExp=0, int firstRun=0, int lastExp=-1, int lastRun=-1)
	Generate and import a (random, fake) channel mask for testing.
void	importPmtQEData (std::string fileName, std::string treeName="qePmtData", int firstExp=0, int firstRun=0, int lastExp=-1, int lastRun=-1)
	Import PMT Quantum Efficiency data to database.
void	importPmtGainData (std::string fileName, std::string treeName="gainPmtData", int firstExp=0, int firstRun=0, int lastExp=-1, int lastRun=-1)
	Import PMT gain parameters data to database.
void	importPmtInstallationData (std::string fileName, std::string treeName="installationPmtData", int firstExp=0, int firstRun=0, int lastExp=-1, int lastRun=-1)
	Import PMT installation data to database.
void	importPmtObsoleteData (std::string fileName, std::string treeName="obsPmtData", int firstExp=0, int firstRun=0, int lastExp=-1, int lastRun=-1)
	Import PMT specifications from Hamamatsu (not to be used!)
void	importPmtTTSPar (std::string fileName, std::string treeName="ttsPmtPar", int firstExp=0, int firstRun=0, int lastExp=-1, int lastRun=-1)
	Import gaussians fitting the TTS distributions.
void	importPmtTTSHisto (std::string fileName, std::string treeName="ttsPmtHisto", int firstExp=0, int firstRun=0, int lastExp=-1, int lastRun=-1)
	Import histograms used for PMT TTS determination.
void	importPmtPulseHeightFitResult (std::string fileName, int firstExp=0, int firstRun=0, int lastExp=-1, int lastRun=-1)
	Import fit results of pulse height disribution for channel gain and threshold efficiency.
void	exportPmtTTSHisto (std::string outFileName="RetrievedHistos.root")
	Example of exporting TTS histograms.
void	importFrontEndSettings (int lookback, int readoutWin, int extraWin, int offset, int expNo, int firstRun, int lastRun)
	Import front-end settings.
void	importDummyCalModuleAlignment (int firstExp=0, int firstRun=0, int lastExp=-1, int lastRun=-1)
	import a dummy payload of TOPCalModuleAlignment DB objects
void	importDummyCalModuleT0 (int firstExp=0, int firstRun=0, int lastExp=-1, int lastRun=-1)
	import a dummy payload of TOPCalCalModuleT0 DB objects
void	importDummyCalChannelT0 (int firstExp=0, int firstRun=0, int lastExp=-1, int lastRun=-1)
	import a dummy payload of TOPCalChannelT0 DB objects
void	importDummyCalTimebase (int firstExp=0, int firstRun=0, int lastExp=-1, int lastRun=-1)
	import a dummy payload of TOPCalTimebase DB objects
void	importDummyCalChannelNoise (int firstExp=0, int firstRun=0, int lastExp=-1, int lastRun=-1)
	import a dummy payload of TOPCalChannelNoise DB objects
void	importDummyCalChannelPulseHeight (int firstExp=0, int firstRun=0, int lastExp=-1, int lastRun=-1)
	import a dummy payload of TOPCalChannelPulseHeight DB objects
void	importDummyCalChannelRQE (int firstExp=0, int firstRun=0, int lastExp=-1, int lastRun=-1)
	import a dummy payload of TOPCalChannelRQE DB objects
void	importDummyCalChannelThresholdEff (int firstExp=0, int firstRun=0, int lastExp=-1, int lastRun=-1)
	import a dummy payload of TOPCalChannelThresholdEff DB objects
void	importDummyCalChannelThreshold (int firstExp=0, int firstRun=0, int lastExp=-1, int lastRun=-1)
	import a dummy payload of TOPCalChannelThreshold DB objects
void	importDummyCalCommonT0 (int firstExp=0, int firstRun=0, int lastExp=-1, int lastRun=-1)
	import a dummy payload of TOPCalCommonT0 DB objects
void	importDummyCalIntegratedCharge (int firstExp=0, int firstRun=0, int lastExp=-1, int lastRun=-1)
	import a dummy payload of TOPCalIntegratedCharge DB objects
void	importDummyCalAsicShift (int firstExp=0, int firstRun=0, int lastExp=-1, int lastRun=-1)
	import a dummy payload of TOPCalAsicShift DB objects
void	correctTOPPmtQE ()
	correct QE values in database for the reflection on window surface be sure that you run the function only once!
void	importTimeWalk (PyObject *list, double a, double b, int firstExp=0, int firstRun=0, int lastExp=-1, int lastRun=-1)
	payload TOPCalTimeWalk import parameters for time-walk correction and electronic time resolution tuning
void	importTest (int runNumber, double syncTimeBase)
	for testing purposes only!
void	importTest ()
	for testing purposes only!
bool	setGains (const std::vector< float > &gains, float error=0)
	Set gains.
int	setPedestals (const TProfile *profile, double average=0)
	Set pedestals from profile histogram with c_WindowSize bins.
unsigned	getOptimizedOffset (const std::vector< unsigned > &values, const std::vector< unsigned > &errors, unsigned maxDif, unsigned maxErr)
	Return the offset that can allow for the maximal number of good pedestal samples.
void	setT0 (int moduleID, unsigned asic, double T0)
	Sets calibration for a single ASIC and switches status to calibrated.
void	setUnusable (int moduleID, unsigned asic)
	Switches calibration status to unusable to flag badly calibrated constant.
double	getT0 (int moduleID, unsigned asic) const
	Returns T0 of a single asic.
bool	isCalibrated (int moduleID, unsigned asic) const
	Returns calibration status.
bool	isDefault (int moduleID, unsigned asic) const
	Returns calibration status.
bool	isUnusable (int moduleID, unsigned asic) const
	Returns calibration status.
void	setStatus (int moduleID, unsigned channel, EStatus status)
	Sets the status for a single channel.
void	setActive (int moduleID, unsigned channel)
	Sets a specific channel as active.
void	setDead (int moduleID, unsigned channel)
	Sets a specific channel as dead.
void	setNoisy (int moduleID, unsigned channel)
	Sets a specific channel as noisy.
EStatus	getStatus (int moduleID, unsigned channel) const
	Returns the status of a single channel.
bool	isActive (int moduleID, unsigned channel) const
	Returns false if the channel is dead or noisy, and true is the channel is active.
bool	check (const int module, const unsigned channel) const
	Check input module and channel arguments are sane.
int	getNumOf (EStatus check) const
	Counts and returns the number of channels having a given status.
int	getNumOf (EStatus check, int moduleID) const
	Counts and returns the number of channels od a module having a given status.
void	setT0 (int moduleID, unsigned channel, double T0, double errT0)
	Sets calibration for a single channel and switches status to calibrated.
void	setUnusable (int moduleID, unsigned channel)
	Switches calibration status to unusable to flag badly calibrated constant.
void	suppressAverage ()
	Subtracts arithmetic average of a module from constants whose status is not c_Default.
double	getT0 (int moduleID, unsigned channel) const
	Returns T0 of a single channel.
double	getT0Error (int moduleID, unsigned channel) const
	Returns error on T0 of a single channel.
bool	isCalibrated (int moduleID, unsigned channel) const
	Returns calibration status.
bool	isDefault (int moduleID, unsigned channel) const
	Returns calibration status.
bool	isUnusable (int moduleID, unsigned channel) const
	Returns calibration status.
void	setAlpha (int moduleID, double alpha, double errAlpha)
	Sets the angle alpha on a single module.
void	setBeta (int moduleID, double beta, double errBeta)
	Sets the angle beta on a single module.
void	setGamma (int moduleID, double gamma, double errGamma)
	Sets the angle gamma on a single module.
void	setX (int moduleID, double x, double errX)
	Sets the displacement x on a single module.
void	setY (int moduleID, double y, double errY)
	Sets the displacement y on a single module.
void	setZ (int moduleID, double z, double errZ)
	Sets the displacement z on a single module.
void	setCalibrated (int moduleID)
	Switches calibration status to calibrated.
void	setUnusable (int moduleID)
	Switches calibration status to unusable to flag badly calibrated constant.
double	getAlpha (int moduleID) const
	Gets the angle alpha on a single module.
double	getBeta (int moduleID) const
	Gets the angle beta on a single module.
double	getGamma (int moduleID) const
	Gets the angle gamma on a single module.
double	getX (int moduleID) const
	Returns the shift x on a single module.
double	getY (int moduleID) const
	Returns the shift y on a single module.
double	getZ (int moduleID) const
	Returns the shift z on a single module.
double	getAlphaErr (int moduleID) const
	Returns the error on alpha on a single module.
double	getBetaErr (int moduleID) const
	Returns the error on beta on a single module.
double	getGammaErr (int moduleID) const
	Returns the error on gamma on a single module.
double	getXErr (int moduleID) const
	Returns the error on x on a single module.
double	getYErr (int moduleID) const
	Returns the error on y on a single module.
double	getZErr (int moduleID) const
	Returns the error on z on a single module.
bool	isCalibrated (int moduleID) const
	Returns calibration status.
bool	isDefault (int moduleID) const
	Returns calibration status.
bool	isUnusable (int moduleID) const
	Returns calibration status.
bool	areAllCalibrated () const
	Returns true if all modules are calibrated.
bool	areAllPrecise (double spatialPrecision, double angularPrecision) const
	Returns true if calibration precision for all modules is within specified values.
void	setTransformations () const
	Sets the transformation cache.
const ROOT::Math::Transform3D &	getTransformation (int moduleID) const
	Returns transformation from local to nominal frame.
void	setT0 (int moduleID, double T0, double errT0)
	Sets calibration for a single module and switches status to calibrated.
void	setUnusable (int moduleID)
	Switches calibration status to unusable to flag badly calibrated constant.
void	suppressAverage ()
	Subtracts arithmetic average from constants whose status is not c_Default.
double	getT0 (int moduleID) const
	Returns T0 of a module.
double	getT0Error (int moduleID) const
	Returns error on T0 of a module.
bool	isCalibrated (int moduleID) const
	Returns calibration status.
bool	isDefault (int moduleID) const
	Returns calibration status.
bool	isUnusable (int moduleID) const
	Returns calibration status.
void	set (int slot, const TH1F *photonYields, const TH1F *backgroundYields, const TH1F *alphaRatio, const TH1F *activePixels, const TH2F *pulseHeights, const TH1F *muonZ)
	Sets the data of a given slot.
void	copyContent (const TH1F *input, TH2F &output)
	Copy content of 1D histogram into 2D histogram.
const TH2F *	getPhotonYields (int slot) const
	Returns a 2D histogram of photon pixel yields.
const TH2F *	getBackgroundYields (int slot) const
	Returns a 2D histogram of background pixel yields.
const TH2F *	getAlphaRatio (int slot) const
	Returns a 2D histogram of equalized pixel alpha ratio.
const TH2F *	getActivePixels (int slot) const
	Returns a 2D histogram of active pixels.
const TH2F *	getPulseHeights (int slot) const
	Returns a 2D histogram of pixel pulse-heights.
const TH1F *	getMuonZ (int slot) const
	Returns z distribution of tracks used to determine pixel yields.
void	append (unsigned scrodID, unsigned channel, const std::vector< double > &sampleTimes, bool replace=true)
	Appends calibration of a single ASIC channel.
const TOPSampleTimes *	getSampleTimes (unsigned scrodID, unsigned channel) const
	Returns calibration constants for given SCROD and channel.
bool	isAvailable (unsigned scrodID, unsigned channel) const
	Checks if calibration is available.
void	createMap () const
	Creates a map in cache.
void	set (const std::vector< double > &timeWalkParams, double noise, double quadratic)
	Sets all the parameters and switches status to calibrated.
double	getTimeWalk (int pulseHeight) const
	Returns time-walk at given pulse height.
double	getSigmaSq (int pulseHeight) const
	Returns an excess of electronic time resolution at given pulse height.
double	getSigma (int pulseHeight) const
	Returns an excess of electronic time resolution at given pulse height.
void	setWindowShifts (std::vector< int > shifts)
	Sets window shifts.
virtual bool	isConsistent () const override
	Check for consistency of data members.
virtual void	print (const std::string &title="Bar segment geometry parameters") const override
	Print the content of the class.
std::vector< std::pair< double, double > >	getBrokenGlueContour () const
	Returns the x-y contour of broken glue.
void	constructContour (double A, double B, double fraction, double angle, std::vector< std::pair< double, double > > &contour) const
	Construct a 2D contour.
virtual void	print (const std::string &title) const
	Print the content of the class.
void	printUnderlined (const std::string &title) const
	Print the content of the class.
virtual void	printSurface (const GeoOpticalSurface &surface) const
	Print the content of optical surface.
double	getReflectivity (const GeoOpticalSurface &surface, double energy) const
	Returns reflectivity of optical surface at given photon energy.
bool	isConsistent () const override
	Check for consistency of data members.
virtual void	print (const std::string &title="QBB cold plate geometry parameters") const override
	Print the content of the class.
bool	isConsistent () const override
	Check for consistency of data members.
virtual void	print (const std::string &title="QBB forward end plate geometry parameters") const override
	Print the content of the class.
bool	isConsistent () const override
	Check for consistency of data members.
virtual void	print (const std::string &title="Front-end geometry parameters") const override
	Print the content of the class.
void	appendContour (std::vector< std::pair< double, double > > &contour, double y0, bool fromLeftToRight) const
	Appends curved surface to contour.
bool	isConsistent () const override
	Check for consistency of data members.
virtual void	print (const std::string &title="QBB honeycomb panel geometry parameters") const override
	Print the content of the class.
void	appendModule (const TOPGeoModule &module)
	Appends module (if its ID differs from already appended modules)
bool	isModuleIDValid (int moduleID) const
	Checks if module exists in m_modules.
const TOPGeoModule &	getModule (int moduleID) const
	Returns module.
const TOPNominalTTS &	getTTS (unsigned type) const
	Returns time transition spread of a given PMT type.
double	getPDETuningFactor (unsigned type) const
	Returns photon detection efficiency tuning factor of a given PMT type.
double	getInnerRadius () const
	Returns inner radius of the volume devoted to TOP counter.
double	getOuterRadius () const
	Returns outer radius of the volume devoted to TOP counter.
double	getRadius () const
	Returns average radius of modules.
double	getBackwardZ () const
	Returns backward z of the volume devoted to TOP counter.
double	getForwardZ () const
	Returns forward z of the volume devoted to TOP counter.
bool	isConsistent () const override
	Check for consistency of data members.
void	print (const std::string &title="TOP geometry parameters") const override
	Print the content of the class.
double	getZc () const
	Returns spherical mirror center of curvature in z (in local frame of this segment)
bool	isConsistent () const override
	Check for consistency of data members.
void	print (const std::string &title="Mirror segment geometry parameters") const override
	Print the content of the class.
	TOPGeoModule (const TOPGeoModule &module)
	Copy constructor.
TOPGeoModule &	operator= (const TOPGeoModule &module)
	Assignment operator.
	~TOPGeoModule ()
	Destructor.
void	setTransformation () const
	Sets transformation cache.
void	setBrokenGlue (int glueID, double fraction, double angle, const std::string &material)
	Sets glue to be broken (delaminated)
void	setPeelOffRegions (double thickness, const std::string &material)
	Sets parameters of the peel-off cookie volumes.
ROOT::Math::XYZPoint	pointToGlobal (const ROOT::Math::XYZPoint &point) const
	Transforms 3D point from module internal (= nominal & displaced) frame to Belle II frame.
ROOT::Math::XYZVector	momentumToGlobal (const ROOT::Math::XYZVector &momentum) const
	Transforms momentum vector from module internal (= nominal & displaced) frame to Belle II frame.
ROOT::Math::XYZPoint	pointToLocal (const ROOT::Math::XYZPoint &point) const
	Transforms 3D point from Belle II to module internal (= nominal & displaced) frame.
ROOT::Math::XYZVector	momentumToLocal (const ROOT::Math::XYZVector &momentum) const
	Transforms momentum vector from Belle II to module internal (= nominal & displaced) frame.
ROOT::Math::XYZPoint	pointNominalToGlobal (const ROOT::Math::XYZPoint &point) const
	Transforms 3D point from module nominal frame to Belle II frame.
ROOT::Math::XYZVector	momentumNominalToGlobal (const ROOT::Math::XYZVector &momentum) const
	Transforms momentum vector from module nominal frame to Belle II frame.
ROOT::Math::XYZPoint	pointGlobalToNominal (const ROOT::Math::XYZPoint &point) const
	Transforms 3D point from Belle II to module nominal frame.
ROOT::Math::XYZVector	momentumGlobalToNominal (const ROOT::Math::XYZVector &momentum) const
	Transforms momentum vector from Belle II to module nominal frame.
bool	isConsistent () const override
	Check for consistency of data members.
void	print (const std::string &title="Module geometry parameters") const override
	Print the content of the class.
void	print (const std::string &title="Module displacement parameters") const override
	Print the content of the class.
ROOT::Math::Transform3D	getTransformation () const
	Returns transformation from local to nominal frame.
unsigned	getPixelColumn (double x) const
	Converts x-coordinate to pixel column (1-based)
unsigned	getPixelRow (double y) const
	Converts y-coordinate to pixel row (1-based)
unsigned	getPixelID (double x, double y) const
	Converts x and y coordinates to PMT pixel ID (1-based)
double	getX (unsigned col) const
	Returns x coordinate of pixel center.
double	getY (unsigned row) const
	Returns y coordinate of pixel center.
bool	isConsistent () const override
	Check for consistency of data members.
void	print (const std::string &title="MCP-PMT geometry parameters") const override
	Print the content of the class.
unsigned	getPmtID (unsigned row, unsigned col) const
	Converts row and column numbers to PMT ID (1-based)
unsigned	getRow (unsigned pmtID) const
	Converts PMT ID to row number (1-based)
unsigned	getColumn (unsigned pmtID) const
	Converts PMT ID to column number (1-based)
int	getPixelID (double x, double y, unsigned pmtID) const
	Converts (digitizes) x, y and PMT ID to pixel ID (1-based)
int	getPixelID (unsigned pmtID, unsigned pmtPixelID) const
	Returns pixel ID (1-based)
void	generateDecoupledPMTs (double fraction)
	Generate randomly a fraction of PMT's to be optically decoupled.
bool	isPMTDecoupled (unsigned pmtID) const
	Checks if PMT is optically decoupled.
bool	isConsistent () const override
	Check for consistency of data members.
void	print (const std::string &title="PMT array geometry parameters") const override
	Print the content of the class.
void	print (const std::string &title="PMT array displacement parameters") const override
	Print the content of the class.
	UnfoldedWindow (const ROOT::Math::XYVector &orig, const ROOT::Math::XYVector &dir, const ROOT::Math::XYVector &norm, const ROOT::Math::XYVector &slanted)
	constructor
bool	isConsistent () const override
	Check for consistency of data members.
void	print (const std::string &title="Prism geometry parameters") const override
	Print the content of the class.
void	appendPeelOffRegion (unsigned ID, double fraction, double angle)
	Appends peel-off cookie region.
std::vector< std::pair< double, double > >	getPeelOffContour (const PeelOffRegion &region) const
	Returns the x-y contour of the peel-off region.
void	unfold () const
	Unfold prism exit window.
void	reflect (const ROOT::Math::XYVector *points, const ROOT::Math::XYVector *normals, const ROOT::Math::XYVector &orig, const ROOT::Math::XYVector &surf, const ROOT::Math::XYVector &norm, const ROOT::Math::XYVector &slanted, int k, std::vector< UnfoldedWindow > &result) const
	Do unfolding.
bool	isConsistent () const override
	Check for consistency of data members.
virtual void	print (const std::string &title="QBB prism enclosure geometry parameters") const override
	Print the content of the class.
std::vector< std::pair< double, double > >	getForwardContour () const
	Returns forward x-y contour.
void	printForwardContour () const
	Print forward contour (polygon) - for debugging.
std::vector< std::pair< double, double > >	getInnerPanelContour () const
	Returns inner honeycomb panel x-y contour.
void	printInnerPanelContour () const
	Print inner honeycomb panel contour (polygon) - for debugging.
std::vector< std::pair< double, double > >	getOuterPanelContour () const
	Returns outer honeycomb panel x-y contour.
void	printOuterPanelContour () const
	Print outer honeycomb panel contour (polygon) - for debugging.
std::vector< std::pair< double, double > >	getBackwardContour () const
	Returns backward x-y contour.
void	printBackwardContour () const
	Print backward contour (polygon) - for debugging.
std::vector< std::pair< double, double > >	getPrismEnclosureContour () const
	Returns prism enclosure wall x-y contour.
void	printPrismEnclosureContour () const
	Print prism enclosure contour (polygon) - for debugging.
std::vector< std::pair< double, double > >	getBackPlateContour () const
	Returns prism enclosure back plate x-y contour.
void	printBackPlateContour () const
	Print prism enclosure back plate contour (polygon) - for debugging.
std::vector< std::pair< double, double > >	getFrontPlateContour () const
	Returns prism enclosure front plate x-y contour.
void	printFrontPlateContour () const
	Print prism enclosure front plate contour (polygon) - for debugging.
bool	isConsistent () const override
	Check for consistency of data members.
virtual void	print (const std::string &title="QBB geometry parameters") const override
	Print the content of the class.
bool	isConsistent () const override
	Check for consistency of data members.
virtual void	print (const std::string &title="QBB side rails geometry parameters") const override
	Print the content of the class.
void	applyFilterTransmission (const TOPWavelengthFilter &filter)
	Multiplies quantum efficiency data points with filter transmission.
double	getQE (double lambda) const
	Returns quantum efficiency at given photon wavelength using linear interpolation.
double	getMinLambda () const
	Returns wavelength of the first nonzero QE data point.
double	getMaxLambda () const
	Returns wavelength of the last nonzero QE data point.
bool	isConsistent () const override
	Check for consistency of data members.
void	print (const std::string &title="Nominal quantum efficiency") const override
	Print the content of the class.
	TOPNominalTDC (int numWindows, int subBits, double syncTimeBase, int numofBunches, double offset, double pileupTime, double doubleHitResolution, double timeJitter, double efficiency, const std::string &name="TOPNominalTDC")
	Useful constructor (new xml file version)
int	getTDCcount (double time) const
	Converts time to TDC count.
int	getSample (double time) const
	Converts time to sample number.
bool	isSampleValid (int sample) const
	Check for the validity of sample number.
bool	isConsistent () const override
	Check for consistency of data members.
void	print (const std::string &title="Nominal time-to-digit conversion parameters") const override
	Print the content of the class.
void	appendGaussian (double norm, double mean, double sigma)
	Append Gaussian.
double	normalize ()
	Normalize the distribution (fractions)
double	generateTTS () const
	Generate time according to TTS distribution.
bool	isConsistent () const override
	Check for consistency of data members.
void	print (const std::string &title="Nominal TTS distribution") const override
	Print the content of the class.
void	print () const
	Print the class content.
const std::vector< float > &	getQE (unsigned pmtPixel) const
	Returns quantum efficiency data points for a given pixel.
double	getQE (unsigned pmtPixel, double lambda) const
	Returns quantum efficiency for a given pixel and wavelength, using linear interpolation.
const std::vector< float > &	getEnvelopeQE () const
	Returns envelope quantum efficiency data points (maximum over pixels)
double	getEnvelopeQE (double lambda) const
	Returns envelope quantum efficiency for a given wavelength (maximum over pixels), using linear interpolation.
double	getEfficiency (unsigned pmtPixel, double lambda, bool BfieldOn) const
	Returns quantum times collection efficiency for a given pixel and wavelength, using linear interpolation.
double	getCE (bool BfieldOn) const
	Returns collection efficiency.
double	getLambdaLast () const
	Returns wavelenght of the last data point (maximal of pixels)
double	getLambdaLast (unsigned pmtPixel) const
	Returns wavelenght of the last data point for a given pixel.
void	setEnvelopeQE () const
	Sets envelope quantum efficiency.
double	interpolate (double lambda, const std::vector< float > &QE) const
	Interpolate between QE datapoints (linear interpolation).
double	getValue (double x) const
	Returns a value of non-normalized pulse height distribution at x Multiply with getNorm() to get normalized one.
double	getNorm () const
	Returns normalization constant.
double	getMean () const
	Returns distribution mean.
double	getRMS () const
	Returns distribution standard deviation.
double	getPeakPosition () const
	Returns the position of distribution maximum.
double	getThresholdEffi (double threshold, double rmsNoise, int n=100) const
	Returns threshold efficiency.
void	setTimeAxis (double syncTimeBase)
	Sets equidistant time axis (uncalibrated).
void	setTimeAxis (const std::vector< double > &sampleTimes, double syncTimeBase)
	Sets time axis from calibration data and switches status to calibrated.
std::vector< double >	getTimeAxis () const
	Returns time axis (sample times)
double	getFullTime (int window, double sample) const
	Returns time with respect to sample 0 of window 0.
double	getSample (int window, double time) const
	Returns sample with respect to sample 0 of the specified ASIC window (inverse of getTime).
double	getTimeBin (int window, int sampleNumber) const
	Returns time bin of a given sample number and window (e.g.
	TOPSignalShape (std::vector< double > shape, double timeBin, double tau, double pole1, double pole2)
	Full constructor.
double	getValue (double t) const
	Returns value at time t of the normalized waveform using interpolator.
double	getPeakingTime () const
	Returns peaking time of the signal.
float	getBulkTransmittance (double lambda) const
	Returns bulk transmittance at given wavelength using linear interpolation.
bool	isConsistent () const override
	Check for consistency of data members.
void	print (const std::string &title="Wavelength filter bulk transmittance") const override
	Print the content of the class.

Variables
static double	s_unit = Unit::cm
	conversion unit for length
static std::string	s_unitName
	conversion unit name

Detailed Description

Typedef Documentation

Iterator

typedef map<unsigned,size_t>::const_iterator Iterator

Iteratior for m_map.

Definition at line 21 of file TOPCalTimebase.cc.

Pair

typedef std::pair< double, double > Pair

Shorthand for std::pair<double, double>

Definition at line 57 of file TOPGeoBarSegment.cc.

Function Documentation

append()

void append	(	unsigned	scrodID,
		unsigned	channel,
		const std::vector< double > &	sampleTimes,
		bool	replace = true )

Appends calibration of a single ASIC channel.

Parameters

scrodID	SCROD ID
channel	channel number within SCROD (0 - 127)
sampleTimes	vector of 256 elements (at least) defining time axis
replace	replace or not if constants for this scrod/channel are already in

Definition at line 24 of file TOPCalTimebase.cc.

  {
    // sanitary checks for NaN's
    for (double t : sampleTimes) {
      if (isnan(t)) {
        B2ERROR("TOPCalTimebase::append: sampleTimes contain NaN's - constants rejected."
                << LogVar("scrodID", scrodID)
                << LogVar("channel", channel));
        return;
      }
    }
 
    // append
    unsigned key = (scrodID << 16) + (channel % 128);
    Iterator it = m_map.find(key);
    if (it == m_map.end()) { // constants not appended yet
      m_sampleTimes.push_back(TOPSampleTimes(scrodID, channel, m_syncTimeBase));
      m_sampleTimes.back().setTimeAxis(sampleTimes, m_syncTimeBase);
      m_map[key] = m_sampleTimes.size() - 1;
    } else { // constants already there
      if (replace) {
        m_sampleTimes[it->second].setTimeAxis(sampleTimes, m_syncTimeBase);
        B2WARNING("TOPCalTimebase::append: constants replaced."
                  << LogVar("scrodID", scrodID)
                  << LogVar("channel", channel));
      } else {
        B2WARNING("TOPCalTimebase::append old constants kept."
                  << LogVar("scrodID", scrodID)
                  << LogVar("channel", channel));
      }
    }
 
  }

appendContour()

void appendContour	(	std::vector< std::pair< double, double > > &	contour,
		double	y0,
		bool	fromLeftToRight ) const

Appends curved surface to contour.

Parameters

contour	2D contour
y0	offset in y0
fromLeftToRight	if true, going from left to right otherwise oposite

Definition at line 21 of file TOPGeoHoneycombPanel.cc.

  {
    if (m_N <= 0) return;
 
    typedef std::pair<double, double> Pair;
 
    double R = getRadius();
    double yc = getMaxThickness() - getMinThickness() - R;
    double x0 = sqrt(R * R - yc * yc);
    if (!fromLeftToRight) x0 = -x0;
    double dx = 2.0 * x0 / m_N;
    double x = -x0;
    contour.push_back(Pair(x, y0));
    for (int i = 1; i < m_N; i++) {
      x += dx;
      double y = y0 + sqrt(R * R - x * x) + yc;
      contour.push_back(Pair(x, y));
    }
    contour.push_back(Pair(x0, y0));
 
  }

appendGaussian()

void appendGaussian	(	double	norm,
		double	mean,
		double	sigma )

Append Gaussian.

Parameters

norm	normalization (area)
mean	mean value [ns]
sigma	sigma (width) [ns]

Definition at line 22 of file TOPNominalTTS.cc.

  {
    if (norm <= 0) {
      B2ERROR("TOPNominalTTS::appendGaussian: normalization is non-positive number");
      return;
    }
    Gauss gauss;
    gauss.fraction = norm;
    gauss.position = mean;
    gauss.sigma = sigma;
    m_tts.push_back(gauss);
    m_normalized = false;
  }

appendModule()

void appendModule

(

const TOPGeoModule &

module

)

Appends module (if its ID differs from already appended modules)

Parameters

module

module geometry parameters

Definition at line 22 of file TOPGeometry.cc.

  {
    if (isModuleIDValid(module.getModuleID())) {
      B2ERROR("TOPGeometry::appendModule: module already appended."
              << LogVar("ID", module.getModuleID()));
      return;
    }
    m_modules.push_back(module);
  }

appendPeelOffRegion()

void appendPeelOffRegion	(	unsigned	ID,
		double	fraction,
		double	angle )

Appends peel-off cookie region.

Parameters

ID	region ID (1-based)
fraction	fraction of the area
angle	angle of the area

Definition at line 94 of file TOPGeoPrism.cc.

  {
    for (const auto& region : m_peelOffRegions) {
      if (region.ID == ID) {
        B2ERROR("TOPGeoPrism::appendPeelOffRegion: region already appended."
                << LogVar("region ID", ID));
        return;
      }
    }
    PeelOffRegion region;
    region.ID = ID;
    region.fraction = fraction;
    region.angle = angle;
    double halfSize = (getWidth() - getPeelOffSize()) / 2;
    if (fabs(getPeelOffCenter(region)) > halfSize) {
      B2ERROR("TOPGeoPrism::appendPeelOffRegion: region doesn't fit into prism."
              << LogVar("region ID", ID));
      return;
    }
    m_peelOffRegions.push_back(region);
  }

applyFilterTransmission()

void applyFilterTransmission

(

const TOPWavelengthFilter &

filter

)

Multiplies quantum efficiency data points with filter transmission.

Parameters

filter

wavelength filter

Definition at line 21 of file TOPNominalQE.cc.

  {
    for (size_t i = 0; i < m_QE.size(); i++) {
      double lambda = m_lambdaFirst + m_lambdaStep * i;
      m_QE[i] *= filter.getBulkTransmittance(lambda);
    }
  }

areAllCalibrated()

bool areAllCalibrated

(

)

const

Returns true if all modules are calibrated.

Definition at line 260 of file TOPCalModuleAlignment.cc.

  {
    for (int i = 0; i <  c_numModules; i++) {
      if (m_status[i] != c_Calibrated) return false;
    }
    return true;
  }

areAllPrecise()

bool areAllPrecise	(	double	spatialPrecision,
		double	angularPrecision ) const

Returns true if calibration precision for all modules is within specified values.

Parameters

spatialPrecision	precision for displacements
angularPrecision	precision for rotations

Definition at line 268 of file TOPCalModuleAlignment.cc.

  {
    for (int i = 0; i <  c_numModules; i++) {
      if (m_errAlpha[i] > angularPrecision) return false;
      if (m_errBeta[i] > angularPrecision) return false;
      if (m_errGamma[i] > angularPrecision) return false;
      if (m_errX[i] > spatialPrecision) return false;
      if (m_errY[i] > spatialPrecision) return false;
      if (m_errZ[i] > spatialPrecision) return false;
    }
    return true;
  }

check()

bool check ( const int module, const unsigned channel ) const

private

Check input module and channel arguments are sane.

Definition at line 82 of file TOPCalChannelMask.cc.

  {
    if (module >= c_numModules) {
      B2ERROR("Invalid module number (" << ClassName() << ")");
      return false;
    }
    if (module < 0) {
      B2ERROR("Invalid module number (" << ClassName() << ")");
      return false;
    }
    if (channel >= c_numChannels) {
      B2ERROR("Invalid channel number (" << ClassName() << ")");
      return false;
    }
    return true;
  }

constructContour()

void constructContour ( double A, double B, double fraction, double angle, std::vector< std::pair< double, double > > & contour ) const

protected

Construct a 2D contour.

Parameters

A	dimension in x
B	dimension in y
fraction	surface fraction
angle	angle
contour	clock-wise polygon to return

Definition at line 70 of file TOPGeoBarSegment.cc.

  {
    if (fraction <= 0.) {
      return;
    }
    if (fraction >= 1.) {
      contour.push_back(Pair(-A, B));
      contour.push_back(Pair(A, B));
      contour.push_back(Pair(A, -B));
      contour.push_back(Pair(-A, -B));
      return;
    }
 
    double alpha0 = atan(B / A);
    double halfPi = alpha0 > 0 ? M_PI / 2 : -M_PI / 2;
    int quadrant = 0;
    while (fabs(angle) > alpha0) {
      angle -= halfPi;
      alpha0 = M_PI / 2 - alpha0;
      double a = A;
      A = B;
      B = a;
      quadrant++;
    }
    quadrant = quadrant % 4;
    if (halfPi < 0) quadrant = (4 - quadrant) % 4;
 
    double tanAngle = tan(fabs(angle));
    double S0 = A / B * tanAngle / 2;
    if (fraction <= S0) {
      double x = sqrt(8 * A * B * fraction / tanAngle);
      double y = x * tanAngle;
      if (angle > 0) {
        contour.push_back(Pair(-A, B - y));
        contour.push_back(Pair(-A, B));
        contour.push_back(Pair(-A + x, B));
      } else {
        contour.push_back(Pair(A, B - y));
        contour.push_back(Pair(A - x, B));
        contour.push_back(Pair(A, B));
      }
    } else if (1 - fraction <= S0) {
      double x = sqrt(8 * A * B * (1 - fraction) / tanAngle);
      double y = x * tanAngle;
      if (angle > 0) {
        contour.push_back(Pair(A, -B + y));
        contour.push_back(Pair(A - x, -B));
        contour.push_back(Pair(-A, -B));
        contour.push_back(Pair(-A, B));
        contour.push_back(Pair(A, B));
      } else {
        contour.push_back(Pair(-A + x, -B));
        contour.push_back(Pair(-A, -B + y));
        contour.push_back(Pair(-A, B));
        contour.push_back(Pair(A, B));
        contour.push_back(Pair(A, -B));
      }
    } else {
      double y = (1 - 2 * fraction) * B;
      double dy = angle > 0 ? A * tanAngle : -A * tanAngle;
      contour.push_back(Pair(-A, y  - dy));
      contour.push_back(Pair(-A, B));
      contour.push_back(Pair(A, B));
      contour.push_back(Pair(A, y  + dy));
    }
 
    switch (quadrant) {
      case 0:
        break;
      case 1:
        for (auto& point : contour) {
          double x = -point.second;
          double y = point.first;
          point.first = x;
          point.second = y;
        }
        break;
      case 2:
        for (auto& point : contour) {
          double x = -point.first;
          double y = -point.second;
          point.first = x;
          point.second = y;
        }
        break;
      case 3:
        for (auto& point : contour) {
          double x = point.second;
          double y = -point.first;
          point.first = x;
          point.second = y;
        }
        break;
      default:
        B2ERROR("TOPGeoBarSegment::constructContour: bug!");
    }
 
  }

copyContent()

void copyContent ( const TH1F * input, TH2F & output )

private

Copy content of 1D histogram into 2D histogram.

Parameters

input	1D histogram
output	2D histogram

Definition at line 51 of file TOPCalPhotonYields.cc.

  {
    for (int bin = 1; bin <= input->GetNbinsX(); bin++) {
      int row = (bin - 1) / 64 + 1;
      int col = (bin - 1) % 64 + 1;
      output.SetBinContent(col, row, input->GetBinContent(bin));
      output.SetBinError(col, row, input->GetBinError(bin));
    }
  }

correctTOPPmtQE()

void correctTOPPmtQE

(

)

correct QE values in database for the reflection on window surface be sure that you run the function only once!

see: BII-4230

Definition at line 1405 of file TOPDatabaseImporter.cc.

  {
    B2ERROR("Function disabled since the corrected payload TOPPmtQEs already imported");
    return;
 
    DBArray<TOPPmtQE> pmtQEData;
    DBImportArray<TOPPmtQE> pmtQECorrected;
 
    for (const auto& pmt : pmtQEData) {
      auto* pmtCorr = pmtQECorrected.appendNew(pmt.getSerialNumber(),
                                               pmt.getLambdaFirst(),
                                               pmt.getLambdaStep(),
                                               pmt.getCE(false),
                                               pmt.getCE(true));
      for (unsigned pmtPixel = 1; pmtPixel <= TOPPmtQE::c_NumPmtPixels; pmtPixel++) {
        auto qeData = pmt.getQE(pmtPixel);
        float lambda = pmt.getLambdaFirst();
        float step = pmt.getLambdaStep();
        for (auto& qe : qeData) {
          double n = TOPGeometryPar::Instance()->getPhaseIndex(TOPGeometryPar::c_hc / lambda); // quartz refractive index
          double reflectance = pow((n - 1) / (n + 1), 2);
          qe /= (1 - reflectance);
          lambda += step;
        }
        pmtCorr->setQE(pmtPixel, qeData);
      }
    }
 
    IntervalOfValidity iov(0, 0, -1, -1);
    pmtQECorrected.import(iov);
 
    B2RESULT("Corrected PMT QE data imported to database for "
             << pmtQECorrected.getEntries() << " PMT's.");
 
  }

createMap()

void createMap ( ) const

private

Creates a map in cache.

Definition at line 87 of file TOPCalTimebase.cc.

  {
    for (size_t i = 0; i < m_sampleTimes.size(); i++) {
      const auto& sampleTime = m_sampleTimes[i];
      unsigned key = (sampleTime.getScrodID() << 16) + sampleTime.getChannel();
      m_map[key] = i;
    }
 
    B2DEBUG(29, "TOPCalTimebase: map created, size = " << m_map.size());
  }

exportPmtTTSHisto()

void exportPmtTTSHisto

(

std::string

outFileName = "RetrievedHistos.root"

)

Example of exporting TTS histograms.

Parameters

outFileName

name of the root file where data will be saved

Definition at line 1219 of file TOPDatabaseImporter.cc.

  {
 
    // this is just an example on how to retrieve TTS histograms
    DBArray<TOPPmtTTSHisto> elements("TOPPmtTTSHistos");
 
    TFile file(outFileName.c_str(), "recreate");
 
    // prints serialNum of PMTs and hv setting used, and saves TTS histograms to root file
    for (const auto& element : elements) {
 
      B2INFO("serialNum = " << element.getSerialNumber() << ", HV = " << element.getHV());
      for (int ic = 0; ic < element.getNumOfPixels(); ic++) {
        const auto* ttsHisto = element.getHistogram(ic + 1);
        if (ttsHisto) ttsHisto->Write();
      }
    }
 
    file.Close();
 
    return;
  }

generateDecoupledPMTs()

void generateDecoupledPMTs

(

double

fraction

)

Generate randomly a fraction of PMT's to be optically decoupled.

Parameters

fraction

decoupled fraction

Definition at line 74 of file TOPGeoPMTArray.cc.

  {
    if (fraction <= 0) return;
 
    for (unsigned i = 0; i < getSize(); i++) {
      if (gRandom->Rndm() < fraction) setDecoupledPMT(i + 1);
    }
 
  }

generateFakeChannelMask()

void generateFakeChannelMask	(	double	fractionDead,
		double	fractionHot,
		int	firstExp = 0,
		int	firstRun = 0,
		int	lastExp = -1,
		int	lastRun = -1 )

Generate and import a (random, fake) channel mask for testing.

Parameters

fractionDead	the fraction of dead PMs to randomly assign
fractionHot	the fraction of noisy PMs to randomly assign
firstExp	first experiment number of IOV
firstRun	first run number of IOV
lastExp	first experiment number of IOV
lastRun	last run number of IOV

Definition at line 673 of file TOPDatabaseImporter.cc.

  {
    // declare db object to be imported -- and construct it
    DBImportObjPtr<TOPCalChannelMask> channelMask;
    channelMask.construct();
 
    // set up for loop channel mapper
    auto& chMapper = TOPGeometryPar::Instance()->getChannelMapper();
    const size_t nModules = TOPGeometryPar::Instance()->getGeometry()->getNumModules();
    unsigned ncall = 0;
    unsigned nall = 0;
 
    // loop over module (1-based)
    for (size_t moduleID = 1; moduleID <= nModules; moduleID++) {
 
      // loop over boardStack*carrierBoard*assic*channel to get channel (0 to 512)
      // TODO: get these loop limits from some sensible enum somewhere
      for (int boardStack = 0; boardStack < 4; boardStack++) {
        for (int carrierBoard = 0; carrierBoard < 4; carrierBoard++) {
          for (int asic = 0; asic < 4; asic++) {
            for (int chan = 0; chan < 8; chan++) {
              auto channel = chMapper.getChannel(boardStack, carrierBoard, asic, chan);
              nall++;
              if (gRandom->Rndm() < fractionDead) {
                channelMask->setDead(moduleID, channel);
                ncall++;
              }
              if (gRandom->Rndm() < fractionHot) {
                channelMask->setNoisy(moduleID, channel);
                ncall++;
              }
            }
          }
        }
      }
    } // module
 
    // declare interval of validity
    IntervalOfValidity iov(firstExp, firstRun, lastExp, lastRun);
    channelMask.import(iov);
 
    B2RESULT("Generated and imported a fake channel mask to database for testing: "
             << ncall << "/" << nall);
    return;
  }

generateTTS()

double generateTTS

(

)

const

Generate time according to TTS distribution.

Returns

time [ns]

Definition at line 46 of file TOPNominalTTS.cc.

  {
    double prob = gRandom->Rndm();
    double s = 0;
    for (const auto& tts : m_tts) {
      s = s + tts.fraction;
      if (prob < s) {
        return gRandom->Gaus(tts.position, tts.sigma);
      }
    }
    return 0;
  }

getActivePixels()

const TH2F * getActivePixels

(

int

slot

)

const

Returns a 2D histogram of active pixels.

Parameters

slot

slot ID

Returns

pixel image of active pixels (nullptr if N/A)

Definition at line 86 of file TOPCalPhotonYields.cc.

  {
    unsigned index = slot - 1;
    if (index < m_activePixels.size()) return &m_activePixels[index];
    return 0;
  }

getAlpha()

double getAlpha

(

int

moduleID

)

const

Gets the angle alpha on a single module.

Parameters

moduleID

module ID (1-based)

Returns

alpha rotation angle around x

Definition at line 119 of file TOPCalModuleAlignment.cc.

  {
    unsigned module = moduleID - 1;
    if (module >= c_numModules) {
      B2WARNING("Invalid module number, returning 0 (" << ClassName() << ")");
      return 0;
    }
    return m_alpha[module];
  }

getAlphaErr()

double getAlphaErr

(

int

moduleID

)

const

Returns the error on alpha on a single module.

Parameters

moduleID

module ID (1-based)

Returns

errAlpha error on the alpha angle

Definition at line 179 of file TOPCalModuleAlignment.cc.

  {
    unsigned module = moduleID - 1;
    if (module >= c_numModules) {
      B2WARNING("Invalid module number, returning 0 (" << ClassName() << ")");
      return 0;
    }
    return m_errAlpha[module];
  }

getAlphaRatio()

const TH2F * getAlphaRatio

(

int

slot

)

const

Returns a 2D histogram of equalized pixel alpha ratio.

For optically well coupled pixels this ratio is around 1 and it is around 0.5 for those which are not.

Parameters

slot

slot ID

Returns

pixel image of alpha ratio (nullptr if N/A)

Definition at line 78 of file TOPCalPhotonYields.cc.

  {
    unsigned index = slot - 1;
    if (index < m_alphaRatio.size()) return &m_alphaRatio[index];
    return 0;
  }

getBackgroundYields()

const TH2F * getBackgroundYields

(

int

slot

)

const

Returns a 2D histogram of background pixel yields.

These are inefficiency corrected number of background hits per track in a pixel.

Parameters

slot

slot ID

Returns

pixel image of photon yields (nullptr if N/A)

Definition at line 70 of file TOPCalPhotonYields.cc.

  {
    unsigned index = slot - 1;
    if (index < m_backgroundYields.size()) return &m_backgroundYields[index];
    return 0;
  }

getBackPlateContour()

std::vector< std::pair< double, double > > getBackPlateContour

(

)

const

Returns prism enclosure back plate x-y contour.

Returns

polygon

Definition at line 174 of file TOPGeoQBB.cc.

  {
    typedef std::pair<double, double> Pair;
    std::vector<Pair> contour;
 
    double x1 = getWidth() / 2;
    double x2 = x1 - m_prismEnclosure.getHeight() * tan(m_prismEnclosure.getAngle());
    double y1 = m_outerPanel.getY() + m_outerPanel.getMinThickness();
    double y2 = y1 - m_sideRails.getHeight();
    double y3 = y2 - m_prismEnclosure.getHeight();
    contour.push_back(Pair(-x1, y1));
    contour.push_back(Pair(x1, y1));
    contour.push_back(Pair(x1, y2));
    contour.push_back(Pair(x2, y3));
    contour.push_back(Pair(-x2, y3));
    contour.push_back(Pair(-x1, y2));
 
    return contour;
  }

getBackwardContour()

std::vector< std::pair< double, double > > getBackwardContour

(

)

const

Returns backward x-y contour.

Returns

polygon

Definition at line 106 of file TOPGeoQBB.cc.

  {
    typedef std::pair<double, double> Pair;
    std::vector<Pair> contour;
 
    double x1 = getWidth() / 2;
    double x2 = x1 - m_prismEnclosure.getHeight() * tan(m_prismEnclosure.getAngle());
    double x3 = getWidth() / 6;
    double y1 = m_outerPanel.getY() + m_outerPanel.getMaxThickness();
    double y2 = m_outerPanel.getY() + m_outerPanel.getMinThickness();
    double y3 = y2 - m_sideRails.getHeight();
    double y4 = y3 - m_prismEnclosure.getHeight();
    contour.push_back(Pair(-x1, y2));
    contour.push_back(Pair(-x3, y1));
    contour.push_back(Pair(x3, y1));
    contour.push_back(Pair(x1, y2));
    contour.push_back(Pair(x1, y3));
    contour.push_back(Pair(x2, y4));
    contour.push_back(Pair(-x2, y4));
    contour.push_back(Pair(-x1, y3));
 
    return contour;
  }

getBackwardZ()

double getBackwardZ

(

)

const

Returns backward z of the volume devoted to TOP counter.

Returns

minimal z

Definition at line 102 of file TOPGeometry.cc.

  {
    if (m_modules.empty()) return 0;
 
    double z = m_modules[0].getBackwardZ();
    for (auto& module : m_modules) {
      double tmp = module.getBackwardZ() - module.getPrism().getFullLength();
      if (tmp < z) z = tmp;
    }
    return z;
  }

getBeta()

double getBeta

(

int

moduleID

)

const

Gets the angle beta on a single module.

Parameters

moduleID

module ID (1-based)

Returns

beta rotation angle around y

Definition at line 129 of file TOPCalModuleAlignment.cc.

  {
    unsigned module = moduleID - 1;
    if (module >= c_numModules) {
      B2WARNING("Invalid module number, returning 0 (" << ClassName() << ")");
      return 0;
    }
    return m_beta[module];
  }

getBetaErr()

double getBetaErr

(

int

moduleID

)

const

Returns the error on beta on a single module.

Parameters

moduleID

module ID (1-based)

Returns

errBeta error on the beta angle

Definition at line 189 of file TOPCalModuleAlignment.cc.

  {
    unsigned module = moduleID - 1;
    if (module >= c_numModules) {
      B2WARNING("Invalid module number, returning 0 (" << ClassName() << ")");
      return 0;
    }
    return m_errBeta[module];
  }

getBrokenGlueContour()

std::vector< std::pair< double, double > > getBrokenGlueContour

(

)

const

Returns the x-y contour of broken glue.

Definition at line 59 of file TOPGeoBarSegment.cc.

  {
 
    std::vector<Pair> contour;
    constructContour(getWidth() / 2, getThickness() / 2, m_brokenFraction, m_brokenAngle,
                     contour);
 
    return contour;
  }

getBulkTransmittance()

float getBulkTransmittance

(

double

lambda

)

const

Returns bulk transmittance at given wavelength using linear interpolation.

Parameters

lambda

photon wavelength in [nm]

Returns

bulk transmittance

Definition at line 21 of file TOPWavelengthFilter.cc.

  {
    if (m_transmittances.empty()) return 1; // old payload from DB (filter included in QE)
 
    double dlam = lambda - m_lambdaFirst;
    if (dlam < 0) return 0;
    if (dlam > (m_transmittances.size() - 1) * m_lambdaStep) return m_transmittances.back();
    unsigned i = int(dlam / m_lambdaStep);
    if (i > m_transmittances.size() - 2) return m_transmittances.back();
    return m_transmittances[i] + (m_transmittances[i + 1] - m_transmittances[i]) / m_lambdaStep * (dlam - i * m_lambdaStep);
  }

getCE()

double getCE

(

bool

BfieldOn

)

const

Returns collection efficiency.

Parameters

BfieldOn

true for magnetic field being ON

Returns

collection efficiency

Definition at line 59 of file TOPPmtQE.cc.

  {
    if (BfieldOn) {
      return m_CE_withB;
    } else {
      return m_CE_noB;
    }
  }

getColumn()

unsigned getColumn

(

unsigned

pmtID

)

const

Converts PMT ID to column number (1-based)

Parameters

pmtID

PMT ID (1-based)

Returns

column number or 0 for invalid PMT ID

Definition at line 34 of file TOPGeoPMTArray.cc.

  {
    if (pmtID == 0 or pmtID > getSize()) return 0;
    return (pmtID - 1) % m_numColumns + 1;
  }

getEfficiency()

double getEfficiency	(	unsigned	pmtPixel,
		double	lambda,
		bool	BfieldOn ) const

Returns quantum times collection efficiency for a given pixel and wavelength, using linear interpolation.

Parameters

pmtPixel	pmtPixel number (1-based)
lambda	wavelength in [nm]
BfieldOn	true for magnetic field being ON

Returns

quantum efficiency

Definition at line 50 of file TOPPmtQE.cc.

  {
    if (BfieldOn) {
      return getQE(pmtPixel, lambda) * m_CE_withB;
    } else {
      return getQE(pmtPixel, lambda) * m_CE_noB;
    }
  }

getEnvelopeQE() [1/2]

const std::vector< float > & getEnvelopeQE

(

)

const

Returns envelope quantum efficiency data points (maximum over pixels)

Returns

envelope quantum efficiency data points

Definition at line 38 of file TOPPmtQE.cc.

  {
    setEnvelopeQE();
    return m_envelopeQE;
  }

getEnvelopeQE() [2/2]

double getEnvelopeQE

(

double

lambda

)

const

Returns envelope quantum efficiency for a given wavelength (maximum over pixels), using linear interpolation.

Parameters

lambda

wavelength in [nm]

Returns

quantum efficiency

Definition at line 44 of file TOPPmtQE.cc.

  {
    setEnvelopeQE();
    return interpolate(lambda, m_envelopeQE);
  }

getForwardContour()

std::vector< std::pair< double, double > > getForwardContour

(

)

const

Returns forward x-y contour.

Returns

polygon

Definition at line 21 of file TOPGeoQBB.cc.

  {
    typedef std::pair<double, double> Pair;
    std::vector<Pair> contour;
 
    double x1 = getWidth() / 2;
    double x2 = x1 - m_sideRails.getThickness();
    double x3 = getWidth() / 6;
    double y1 = m_outerPanel.getY() + m_outerPanel.getMaxThickness();
    double y2 = m_outerPanel.getY() + m_outerPanel.getMinThickness();
    double y3 = y2 - m_sideRails.getHeight();
    double y4 = m_innerPanel.getY() - m_innerPanel.getMaxThickness();
    contour.push_back(Pair(-x1, y2));
    contour.push_back(Pair(-x3, y1));
    contour.push_back(Pair(x3, y1));
    contour.push_back(Pair(x1, y2));
    contour.push_back(Pair(x1, y3));
    contour.push_back(Pair(x2, y4));
    contour.push_back(Pair(-x2, y4));
    contour.push_back(Pair(-x1, y3));
 
    return contour;
  }

getForwardZ()

double getForwardZ

(

)

const

Returns forward z of the volume devoted to TOP counter.

Returns

maximal z

Definition at line 114 of file TOPGeometry.cc.

  {
    if (m_modules.empty()) return 0;
 
    double z = m_modules[0].getForwardZ();
    for (auto& module : m_modules) {
      double tmp = module.getForwardZ();
      if (tmp > z) z = tmp;
    }
    return z;
  }

getFrontPlateContour()

std::vector< std::pair< double, double > > getFrontPlateContour

(

)

const

Returns prism enclosure front plate x-y contour.

Returns

polygon

Definition at line 203 of file TOPGeoQBB.cc.

  {
    typedef std::pair<double, double> Pair;
    std::vector<Pair> contour;
 
    double x1 = getWidth() / 2;
    double x2 = x1 - m_prismEnclosure.getHeight() * tan(m_prismEnclosure.getAngle());
    double y1 = m_outerPanel.getY() + m_outerPanel.getMinThickness();
    double y2 = y1 - m_sideRails.getHeight();
    double y3 = y2 - m_prismEnclosure.getHeight();
    contour.push_back(Pair(x1, y2));
    contour.push_back(Pair(x2, y3));
    contour.push_back(Pair(-x2, y3));
    contour.push_back(Pair(-x1, y2));
 
    return contour;
  }

getFullTime()

double getFullTime	(	int	window,
		double	sample ) const

Returns time with respect to sample 0 of window 0.

Note: sample is float - digits that follow the decimal point are used to interpolate the time btw. two samples

Parameters

window	ASIC window number
sample	sample counted from the first one in the specified ASIC window

Returns

time in [ns]

Definition at line 58 of file TOPSampleTimes.cc.

  {
 
    int sampleNum = int(sample);
    if (sample < 0) sampleNum--;
    double frac = sample - sampleNum;
 
    sampleNum += window * c_WindowSize;  // counted from window 0
    int n = sampleNum / c_TimeAxisSize;
    int k = sampleNum % c_TimeAxisSize;
    if (k < 0) {
      n--;
      k += c_TimeAxisSize;
    }
 
    double time = n * getTimeRange() + m_timeAxis[k]; // from sample 0 window 0
    time += (m_timeAxis[k + 1] - m_timeAxis[k]) * frac; // add fraction
 
    return time;
  }

getGamma()

double getGamma

(

int

moduleID

)

const

Gets the angle gamma on a single module.

Parameters

moduleID

module ID (1-based)

Returns

gamma rotation angle around z

Definition at line 139 of file TOPCalModuleAlignment.cc.

  {
    unsigned module = moduleID - 1;
    if (module >= c_numModules) {
      B2WARNING("Invalid module number, returning 0 (" << ClassName() << ")");
      return 0;
    }
    return m_gamma[module];
  }

getGammaErr()

double getGammaErr

(

int

moduleID

)

const

Returns the error on gamma on a single module.

Parameters

moduleID

module ID (1-based)

Returns

errGamma error on the gamma angle

Definition at line 199 of file TOPCalModuleAlignment.cc.

  {
    unsigned module = moduleID - 1;
    if (module >= c_numModules) {
      B2WARNING("Invalid module number, returning 0 (" << ClassName() << ")");
      return 0;
    }
    return m_errGamma[module];
  }

getInnerPanelContour()

std::vector< std::pair< double, double > > getInnerPanelContour

(

)

const

Returns inner honeycomb panel x-y contour.

Returns

polygon

Definition at line 54 of file TOPGeoQBB.cc.

  {
    typedef std::pair<double, double> Pair;
    std::vector<Pair> contour;
 
    double x1 = getPanelWidth() / 2;
    double y1 = m_innerPanel.getY();
    double y2 = y1 - m_innerPanel.getMaxThickness();
    contour.push_back(Pair(-x1, y1));
    contour.push_back(Pair(x1, y1));
    contour.push_back(Pair(x1, y2));
    m_innerPanel.appendContour(contour, y2, false);
    contour.push_back(Pair(-x1, y2));
 
    return contour;
  }

getInnerRadius()

double getInnerRadius

(

)

const

Returns inner radius of the volume devoted to TOP counter.

Returns

inner radius

Definition at line 64 of file TOPGeometry.cc.

  {
    if (m_modules.empty()) return 0;
 
    double R = m_modules[0].getRadius();
    for (auto& module : m_modules) {
      double tmp = module.getRadius() + module.getBarThickness() / 2
                   - module.getPrism().getExitThickness();
      if (tmp < R) R = tmp;
    }
    return R;
  }

getLambdaLast() [1/2]

double getLambdaLast

(

)

const

Returns wavelenght of the last data point (maximal of pixels)

Returns

wavelength in [nm]

Definition at line 68 of file TOPPmtQE.cc.

  {
    size_t size = 0;
    for (unsigned pixel = 0; pixel < c_NumPmtPixels; pixel++) {
      const auto& QE = m_QE[pixel];
      size = std::max(size, QE.size());
    }
    return m_lambdaFirst + (size - 1) * m_lambdaStep;
  }

getLambdaLast() [2/2]

double getLambdaLast

(

unsigned

pmtPixel

)

const

Returns wavelenght of the last data point for a given pixel.

Parameters

pmtPixel

pmtPixel number (1-based)

Returns

wavelength in [nm]

Definition at line 78 of file TOPPmtQE.cc.

  {
    pmtPixel--;
    if (pmtPixel < c_NumPmtPixels) {
      return m_lambdaFirst + (m_QE[pmtPixel].size() - 1) * m_lambdaStep;
    }
    return 0;
  }

getMaxLambda()

double getMaxLambda

(

)

const

Returns wavelength of the last nonzero QE data point.

Returns

wavelength of the last nonzero QE data point

Definition at line 49 of file TOPNominalQE.cc.

  {
    for (size_t i = m_QE.size() - 1; i < m_QE.size(); i--) {
      if (m_QE[i] > 0) return m_lambdaFirst + m_lambdaStep * i;
    }
    return m_lambdaFirst;
  }

getMean()

double getMean

(

)

const

Returns distribution mean.

Returns

mean [ADC counts]

Definition at line 32 of file TOPPulseHeightPar.cc.

  {
    return std::tgamma((p1 + 2) / p2) / std::tgamma((p1 + 1) / p2) * x0;
  }

getMinLambda()

double getMinLambda

(

)

const

Returns wavelength of the first nonzero QE data point.

Returns

wavelength of the first nonzero QE data point

Definition at line 40 of file TOPNominalQE.cc.

  {
    for (size_t i = 0; i < m_QE.size(); i++) {
      if (m_QE[i] > 0) return m_lambdaFirst + m_lambdaStep * i;
    }
    return getLambdaLast();
  }

getModule()

const TOPGeoModule & getModule

(

int

moduleID

)

const

Returns module.

Parameters

moduleID

valid module ID (1-based)

Returns

module geometry parameters

Definition at line 42 of file TOPGeometry.cc.

  {
    for (const auto& module : m_modules) {
      if (module.getModuleID() == moduleID) return module;
    }
    B2FATAL("TOPGeometry::getModule: invalid module, ID = " << moduleID);
  }

getMuonZ()

const TH1F * getMuonZ

(

int

slot

)

const

Returns z distribution of tracks used to determine pixel yields.

Parameters

slot

slot ID

Returns

z distribution in local frame (nullptr if N/A)

Definition at line 102 of file TOPCalPhotonYields.cc.

  {
    unsigned index = slot - 1;
    if (index < m_muonZ.size()) return &m_muonZ[index];
    return 0;
  }

getNorm()

double getNorm

(

)

const

Returns normalization constant.

Returns

normalization [1/ADC counts]

Definition at line 27 of file TOPPulseHeightPar.cc.

  {
    return p2 / std::tgamma((p1 + 1) / p2) / x0;
  }

getNumOf() [1/2]

int getNumOf ( EStatus check ) const

private

Counts and returns the number of channels having a given status.

Parameters

check

status to be checked

Returns

number of channels of this status

Definition at line 99 of file TOPCalChannelMask.cc.

  {
    int n = 0;
    for (const auto& statuses : m_status) {
      for (const auto status : statuses) {
        if (status == check) n++;
      }
    }
    return n;
  }

getNumOf() [2/2]

int getNumOf ( EStatus check, int moduleID ) const

private

Counts and returns the number of channels od a module having a given status.

Parameters

check	status to be checked
moduleID	module (slot) ID (1-based)

Returns

number of channels of this status

Definition at line 110 of file TOPCalChannelMask.cc.

  {
    unsigned module = moduleID - 1;
    if (module < c_numModules) {
      int n = 0;
      const auto& statuses = m_status[module];
      for (const auto status : statuses) {
        if (status == check) n++;
      }
      return n;
    }
    B2ERROR("Invalid module number (" << ClassName() << ")");
    return 0;
  }

getOptimizedOffset()

unsigned getOptimizedOffset ( const std::vector< unsigned > & values, const std::vector< unsigned > & errors, unsigned maxDif, unsigned maxErr )

private

Return the offset that can allow for the maximal number of good pedestal samples.

Parameters

values	pedestal values
errors	pedestal errors
maxDif	maximal number that can be stored in c_Bits
maxErr	maximal number that can be stored in 16 - c_Bits

Returns

offset value

Definition at line 65 of file TOPASICPedestals.cc.

  {
 
    std::set<unsigned> sortedValues;
    for (unsigned i = 0; i < values.size(); i++) {
      if (errors[i] > maxErr) continue;
      sortedValues.insert(values[i]);
    }
 
    unsigned nprev = sortedValues.size();
    unsigned nl = 0;
    for (auto it = sortedValues.begin(); it != sortedValues.end(); ++it) {
      unsigned nr = 0;
      for (auto back = sortedValues.crbegin(); *back - *it > maxDif; ++back) {
        nr++;
      }
      unsigned nbad = nl + nr;
      if (nbad == 0) return *it;
      if (nbad >= nprev) return *(--it);
      nprev = nbad;
      nl++;
    }
 
    return *sortedValues.begin();
 
  }

getOuterPanelContour()

std::vector< std::pair< double, double > > getOuterPanelContour

(

)

const

Returns outer honeycomb panel x-y contour.

Returns

polygon

Definition at line 80 of file TOPGeoQBB.cc.

  {
    typedef std::pair<double, double> Pair;
    std::vector<Pair> contour;
 
    double x1 = getPanelWidth() / 2;
    double y1 = m_outerPanel.getY();
    double y2 = y1 + m_outerPanel.getMinThickness();
    contour.push_back(Pair(-x1, y2));
    m_outerPanel.appendContour(contour, y2, true);
    contour.push_back(Pair(x1, y2));
    contour.push_back(Pair(x1, y1));
    contour.push_back(Pair(-x1, y1));
 
    return contour;
  }

getOuterRadius()

double getOuterRadius

(

)

const

Returns outer radius of the volume devoted to TOP counter.

Returns

outer radius

Definition at line 77 of file TOPGeometry.cc.

  {
    if (m_modules.empty()) return 0;
 
    double R = m_modules[0].getRadius();
    for (auto& module : m_modules) {
      double x = module.getBarWidth() / 2;
      double y = module.getRadius() + module.getBarThickness() / 2;
      double tmp = sqrt(x * x + y * y);
      if (tmp > R) R = tmp;
    }
    return R;
  }

getPDETuningFactor()

double getPDETuningFactor

(

unsigned

type

)

const

Returns photon detection efficiency tuning factor of a given PMT type.

Parameters

type

PMT type

Returns

tuning factor

Definition at line 57 of file TOPGeometry.cc.

  {
    std::map<unsigned, float>::const_iterator it = m_tuneFactorsPDE.find(type);
    if (it == m_tuneFactorsPDE.end()) return 1.0;
    return it->second;
  }

getPeakingTime()

double getPeakingTime

(

)

const

Returns peaking time of the signal.

Returns

peaking time [ns]

Definition at line 93 of file TOPSignalShape.cc.

  {
    if (m_peakTime == 0) {
      double t1 = 0;
      double t2 = 1;
      while (getValue(t2) > 0.5) t2 += 1;
      for (int i = 0; i < 20; i++) {
        double t = (t1 + t2) / 2;
        if (m_interpolator->Derivative(t) > 0) {
          t1 = t;
        } else {
          t2 = t;
        }
      }
      m_peakTime = (t1 + t2) / 2;
    }
 
    return m_peakTime;
  }

getPeakPosition()

double getPeakPosition

(

)

const

Returns the position of distribution maximum.

Returns

peak position [ADC counts]

Definition at line 44 of file TOPPulseHeightPar.cc.

  {
    return std::pow(p1 / p2, 1 / p2) * x0;
  }

getPeelOffContour()

std::vector< Pair > getPeelOffContour

(

const PeelOffRegion &

region

)

const

Returns the x-y contour of the peel-off region.

Parameters

region

peel-off region

Returns

clock-wise polygon

Definition at line 118 of file TOPGeoPrism.cc.

  {
    std::vector<Pair> contour;
    constructContour(getPeelOffSize() / 2, getExitThickness() / 2,
                     region.fraction, region.angle,
                     contour);
 
    return contour;
  }

getPhotonYields()

const TH2F * getPhotonYields

(

int

slot

)

const

Returns a 2D histogram of photon pixel yields.

These are background subtracted and inefficiency corrected number of photons per track in a pixel.

Parameters

slot

slot ID

Returns

pixel image of photon yields (nullptr if N/A)

Definition at line 62 of file TOPCalPhotonYields.cc.

  {
    unsigned index = slot - 1;
    if (index < m_photonYields.size()) return &m_photonYields[index];
    return 0;
  }

getPixelColumn()

unsigned getPixelColumn

(

double

x

)

const

Converts x-coordinate to pixel column (1-based)

Parameters

x	coordinate

Returns

pixel column or 0 if x outside sensitive area

Definition at line 21 of file TOPGeoPMT.cc.

  {
    x *= s_unit;
    if (fabs(x) >= m_sensSizeX / 2) return 0;
    return m_numColumns - int((x + m_sensSizeX / 2) / m_sensSizeX * m_numColumns);
  }

getPixelID() [1/3]

unsigned getPixelID	(	double	x,
		double	y ) const

Converts x and y coordinates to PMT pixel ID (1-based)

Parameters

x	coordinate
y	coordinate

Returns

pixel ID or 0 if x or y outside sensitive area

Definition at line 36 of file TOPGeoPMT.cc.

  {
    auto col = getPixelColumn(x);
    if (col == 0) return 0;
    auto row = getPixelRow(y);
    if (row == 0) return 0;
 
    return (row - 1) * m_numColumns + col;
  }

getPixelID() [2/3]

int getPixelID	(	double	x,
		double	y,
		unsigned	pmtID ) const

Converts (digitizes) x, y and PMT ID to pixel ID (1-based)

Parameters

x	coordinate in PMT local frame
y	coordinate in PMT local frame
pmtID	PMT ID (1-based)

Returns

pixel ID or 0 if (x,y) outside sensitive area or pmtID invalid

Definition at line 40 of file TOPGeoPMTArray.cc.

  {
    if (pmtID == 0 or pmtID > getSize()) return 0;
    pmtID--;
 
    unsigned col = m_pmt.getPixelColumn(x);
    if (col == 0) return 0;
    col--;
 
    unsigned row = m_pmt.getPixelRow(y);
    if (row == 0) return 0;
    row--;
 
    col += (pmtID % m_numColumns) * m_pmt.getNumColumns();
    row += (pmtID / m_numColumns) * m_pmt.getNumRows();
    return row * m_numColumns * m_pmt.getNumColumns() + col + 1;
  }

getPixelID() [3/3]

int getPixelID	(	unsigned	pmtID,
		unsigned	pmtPixelID ) const

Returns pixel ID (1-based)

Parameters

pmtID	PMT ID (1-based)
pmtPixelID	PMT pixel ID (1-based)

Returns

pixel ID or 0 for invalid input ID's

Definition at line 58 of file TOPGeoPMTArray.cc.

  {
    pmtID--;
    if (pmtID >= getSize()) return 0;
 
    pmtPixelID--;
    if (pmtPixelID >= m_pmt.getNumPixels()) return 0;
 
    unsigned col = pmtPixelID % m_pmt.getNumColumns();
    unsigned row = pmtPixelID / m_pmt.getNumColumns();
 
    col += (pmtID % m_numColumns) * m_pmt.getNumColumns();
    row += (pmtID / m_numColumns) * m_pmt.getNumRows();
    return row * m_numColumns * m_pmt.getNumColumns() + col + 1;
  }

getPixelRow()

unsigned getPixelRow

(

double

y

)

const

Converts y-coordinate to pixel row (1-based)

Parameters

y	coordinate

Returns

pixel row or 0 if y outside sensitive area

Definition at line 29 of file TOPGeoPMT.cc.

  {
    y *= s_unit;
    if (fabs(y) >= m_sensSizeY / 2) return 0;
    return int((y + m_sensSizeY / 2) / m_sensSizeY * m_numRows) + 1;
  }

getPmtID()

unsigned getPmtID	(	unsigned	row,
		unsigned	col ) const

Converts row and column numbers to PMT ID (1-based)

Parameters

row	row number (1-based)
col	column number (1-based)

Returns

PMT ID or 0 for invalid row or column

Definition at line 21 of file TOPGeoPMTArray.cc.

  {
    if (row == 0 or row > m_numRows) return 0;
    if (col == 0 or col > m_numColumns) return 0;
    return (row - 1) * m_numColumns + col;
  }

getPrismEnclosureContour()

std::vector< std::pair< double, double > > getPrismEnclosureContour

(

)

const

Returns prism enclosure wall x-y contour.

Returns

polygon

Definition at line 139 of file TOPGeoQBB.cc.

  {
    typedef std::pair<double, double> Pair;
    std::vector<Pair> contour;
 
    double tanAngle = tan(m_prismEnclosure.getAngle());
    double x1 = getWidth() / 2;
    double x2 = x1 - m_prismEnclosure.getHeight() * tanAngle;
    double x3 = x1 - m_prismEnclosure.getSideThickness();
    double y1 = m_outerPanel.getY() + m_outerPanel.getMinThickness()
                - m_sideRails.getHeight();
    double y2 = y1 - m_prismEnclosure.getHeight();
    double y3 = y2 + m_prismEnclosure.getBottomThickness();
    double x4 = x1 - m_prismEnclosure.getSideThickness() - (y1 - y3) * tanAngle;
    contour.push_back(Pair(x1, y1));
    contour.push_back(Pair(x2, y2));
    contour.push_back(Pair(-x2, y2));
    contour.push_back(Pair(-x1, y1));
    contour.push_back(Pair(-x3, y1));
    contour.push_back(Pair(-x4, y3));
    contour.push_back(Pair(x4, y3));
    contour.push_back(Pair(x3, y1));
 
    return contour;
  }

getPulseHeights()

const TH2F * getPulseHeights

(

int

slot

)

const

Returns a 2D histogram of pixel pulse-heights.

Parameters

slot

slot ID

Returns

pixel pulse-heights (nullptr if N/A)

Definition at line 94 of file TOPCalPhotonYields.cc.

  {
    unsigned index = slot - 1;
    if (index < m_pulseHeights.size()) return &m_pulseHeights[index];
    return 0;
  }

getQE() [1/3]

double getQE

(

double

lambda

)

const

Returns quantum efficiency at given photon wavelength using linear interpolation.

Parameters

lambda

photon wavelength in [nm]

Returns

quantum efficiency

Definition at line 30 of file TOPNominalQE.cc.

  {
    double dlam = lambda - m_lambdaFirst;
    if (dlam < 0 or dlam > (m_QE.size() - 1) * m_lambdaStep) return 0;
    unsigned i = int(dlam / m_lambdaStep);
    if (i > m_QE.size() - 2) return m_QE.back();
    return m_QE[i] + (m_QE[i + 1] - m_QE[i]) / m_lambdaStep * (dlam - i * m_lambdaStep);
  }

getQE() [2/3]

const std::vector< float > & getQE

(

unsigned

pmtPixel

)

const

Returns quantum efficiency data points for a given pixel.

Parameters

pmtPixel

pmtPixel number (1-based)

Returns

quantum efficiency data points (empty vector for invalid pixel)

Definition at line 21 of file TOPPmtQE.cc.

  {
    pmtPixel--;
    if (pmtPixel < c_NumPmtPixels) return m_QE[pmtPixel];
 
    B2ERROR("TOPPmtQE::getQE: invalid pixel ID, returning empty vector");
    m_envelopeQE.clear();
    return m_envelopeQE;
  }

getQE() [3/3]

double getQE	(	unsigned	pmtPixel,
		double	lambda ) const

Returns quantum efficiency for a given pixel and wavelength, using linear interpolation.

Parameters

pmtPixel	pmtPixel number (1-based)
lambda	wavelength in [nm]

Returns

quantum efficiency

Definition at line 31 of file TOPPmtQE.cc.

  {
    pmtPixel--;
    if (pmtPixel < c_NumPmtPixels) return interpolate(lambda, m_QE[pmtPixel]);
    return 0;
  }

getRadius()

double getRadius

(

)

const

Returns average radius of modules.

Returns

radius

Definition at line 91 of file TOPGeometry.cc.

  {
    if (m_modules.empty()) return 0;
 
    double R = 0;
    for (auto& module : m_modules) {
      R += module.getRadius();
    }
    return R / m_modules.size();
  }

getReflectivity()

double getReflectivity	(	const GeoOpticalSurface &	surface,
		double	energy ) const

Returns reflectivity of optical surface at given photon energy.

Parameters

surface	optical surface parameters
energy	photon energy in [eV]

Definition at line 72 of file TOPGeoBase.cc.

  {
    energy *= Unit::eV;
    if (surface.hasProperties()) {
      for (const auto& property : surface.getProperties()) {
        if (property.getName() == "REFLECTIVITY") {
          auto energies = property.getEnergies();
          auto values = property.getValues();
          if (energies.size() < 2) return 0;
          if (energy < energies[0] or energy > energies.back()) {
            B2WARNING("TOPGeoBase::getReflectivity: photon energy out of range - return value not reliable");
          }
          auto spline = TSpline3("tmp", energies.data(), values.data(), energies.size());
          return spline.Eval(energy);
        }
      }
    }
 
    B2ERROR("Optical surface " << surface.getName() << " has no property REFLECTIVITY");
    return 0;
  }

getRMS()

double getRMS

(

)

const

Returns distribution standard deviation.

Returns

standard deviation [ADC counts]

Definition at line 37 of file TOPPulseHeightPar.cc.

  {
    double x2 = std::tgamma((p1 + 3) / p2) / std::tgamma((p1 + 1) / p2) * x0 * x0;
    double x = getMean();
    return std::sqrt(x2 - x * x);
  }

getRow()

unsigned getRow

(

unsigned

pmtID

)

const

Converts PMT ID to row number (1-based)

Parameters

pmtID

PMT ID (1-based)

Returns

row number or 0 for invalid PMT ID

Definition at line 28 of file TOPGeoPMTArray.cc.

  {
    if (pmtID == 0 or pmtID > getSize()) return 0;
    return (pmtID - 1) / m_numColumns + 1;
  }

getSample() [1/2]

int getSample

(

double

time

)

const

Converts time to sample number.

Use isSampleValid(int) to check the validity of sample number

Parameters

time

[ns]

Returns

sample number

Definition at line 67 of file TOPNominalTDC.cc.

  {
    time += m_offset;
    if (time > 0) {
      return int(time / m_sampleWidth);
    } else {
      return int(time / m_sampleWidth) - 1;
    }
  }

getSample() [2/2]

double getSample	(	int	window,
		double	time ) const

Returns sample with respect to sample 0 of the specified ASIC window (inverse of getTime).

Note: sample is float - digits that follow the decimal point are used to interpolate the time btw. two samples

Parameters

window	ASIC window number
time	time with respect to SSTin of specified ASIC window

Returns

sample

Definition at line 80 of file TOPSampleTimes.cc.

  {
    time += window * getTimeRange() / 4.0;
    int n = int(time / getTimeRange());
    if (time < 0) n--;
 
    double t = time - getTimeRange() * n;
    int i1 = 0;
    int i2 = c_TimeAxisSize;
    while (i2 - i1 > 1) {
      int i = (i1 + i2) / 2;
      if (t > m_timeAxis[i]) {
        i1 = i;
      } else {
        i2 = i;
      }
    }
 
    return (n * c_TimeAxisSize - window * c_WindowSize + i1 +
            (t - m_timeAxis[i1]) / (m_timeAxis[i2] - m_timeAxis[i1]));
 
  }

getSampleTimes()

const TOPSampleTimes * getSampleTimes	(	unsigned	scrodID,
		unsigned	channel ) const

Returns calibration constants for given SCROD and channel.

If calibration is not available, return equidistant sample times.

Parameters

scrodID	SCROD ID
channel	channel number (within SCROD or within module)

Returns

calibration constants (always valid pointer)

Definition at line 61 of file TOPCalTimebase.cc.

  {
    if (m_map.empty()) createMap();
 
    unsigned key = (scrodID << 16) + (channel % 128);
    Iterator it = m_map.find(key);
    if (it == m_map.end()) {
      if (!m_sampleTime) m_sampleTime = new TOPSampleTimes(0, 0, m_syncTimeBase);
      return m_sampleTime;
    }
    return &m_sampleTimes[it->second];
  }

getSigma()

double getSigma

(

int

pulseHeight

)

const

Returns an excess of electronic time resolution at given pulse height.

Parameters

pulseHeight

pulse height [ADC counts]

Returns

excess sigma [ns]

Definition at line 48 of file TOPCalTimeWalk.cc.

  {
    return sqrt(getSigmaSq(pulseHeight));
  }

getSigmaSq()

double getSigmaSq

(

int

pulseHeight

)

const

Returns an excess of electronic time resolution at given pulse height.

Parameters

pulseHeight

pulse height [ADC counts]

Returns

excess sigma squared [ns^2]

Definition at line 41 of file TOPCalTimeWalk.cc.

  {
    if (pulseHeight <= 0) return 0;
    double x = pulseHeight;
    return pow(m_noise / x, 2) + pow(m_quadratic * x * x, 2);
  }

getStatus()

TOPCalChannelMask::EStatus getStatus	(	int	moduleID,
		unsigned	channel ) const

Returns the status of a single channel.

Parameters

moduleID	module ID (1-based)
channel	hardware channel number (0-based)

Returns

status channel status (1: active 0: dead -1: noisy )

Definition at line 61 of file TOPCalChannelMask.cc.

  {
    int module = moduleID - 1;
    if (!check(module, channel)) {
      B2WARNING("Returning dead channel value");
      return c_Dead;
    }
    return m_status[module][channel];
  }

getT0() [1/3]

double getT0	(	int	moduleID,
		unsigned	asic ) const

Returns T0 of a single asic.

Parameters

moduleID	module ID (1-based)
asic	ASIC number within a module (0 - 63)

Returns

T0

Definition at line 51 of file TOPCalAsicShift.cc.

  {
    unsigned module = moduleID - 1;
    if (module >= c_numModules) {
      B2WARNING("Invalid module number, returning 0 (" << ClassName() << ")");
      return 0;
    }
    if (asic >= c_numAsics) {
      B2WARNING("Invalid asic number, returning 0 (" << ClassName() << ")");
      return 0;
    }
    return m_T0[module][asic];
  }

getT0() [2/3]

double getT0	(	int	moduleID,
		unsigned	channel ) const

Returns T0 of a single channel.

Parameters

moduleID	module ID (1-based)
channel	hardware channel number (0-based)

Returns

T0

Definition at line 73 of file TOPCalChannelT0.cc.

  {
    unsigned module = moduleID - 1;
    if (module >= c_numModules) {
      B2WARNING("Invalid module number, returning 0 (" << ClassName() << ")");
      return 0;
    }
    if (channel >= c_numChannels) {
      B2WARNING("Invalid channel number, returning 0 (" << ClassName() << ")");
      return 0;
    }
    return m_T0[module][channel];
  }

getT0() [3/3]

double getT0

(

int

moduleID

)

const

Returns T0 of a module.

Parameters

moduleID

module ID (1-based)

Returns

T0

Definition at line 63 of file TOPCalModuleT0.cc.

  {
    unsigned module = moduleID - 1;
    if (module >= c_numModules) {
      B2WARNING("Invalid module number, returning 0 (" << ClassName() << ")");
      return 0;
    }
    return m_T0[module];
  }

getT0Error() [1/2]

double getT0Error	(	int	moduleID,
		unsigned	channel ) const

Returns error on T0 of a single channel.

Parameters

moduleID	module ID (1-based)
channel	hardware channel number (0-based)

Returns

error on T0

Definition at line 88 of file TOPCalChannelT0.cc.

  {
    unsigned module = moduleID - 1;
    if (module >= c_numModules) {
      B2WARNING("Invalid module number, returning 0 (" << ClassName() << ")");
      return 0;
    }
    if (channel >= c_numChannels) {
      B2WARNING("Invalid channel number, returning 0 (" << ClassName() << ")");
      return 0;
    }
    return m_errT0[module][channel];
  }

getT0Error() [2/2]

double getT0Error

(

int

moduleID

)

const

Returns error on T0 of a module.

Parameters

moduleID

module ID (1-based)

Returns

error on T0

Definition at line 74 of file TOPCalModuleT0.cc.

  {
    unsigned module = moduleID - 1;
    if (module >= c_numModules) {
      B2WARNING("Invalid module number, returning 0 (" << ClassName() << ")");
      return 0;
    }
    return m_errT0[module];
  }

getTDCcount()

int getTDCcount

(

double

time

)

const

Converts time to TDC count.

For times being outside TDC range, TDC overflow value is returned.

Parameters

time

[ns]

Returns

TDC count

Definition at line 58 of file TOPNominalTDC.cc.

  {
    time += m_offset;
    int overflow = getOverflowValue();
    if (time < 0) return overflow;
    if (time > overflow * m_binWidth) return overflow;
    return int(time / m_binWidth);
  }

getThresholdEffi()

double getThresholdEffi	(	double	threshold,
		double	rmsNoise,
		int	n = 100 ) const

Returns threshold efficiency.

Parameters

threshold	threshold value [ADC counts]
rmsNoise	r.m.s of electronic noise [ADC counts]
n	number of steps in numerical integration

Returns

threshold efficiency

Definition at line 49 of file TOPPulseHeightPar.cc.

  {
    if (x0 == 0) return 0;
    double s = 0;
    double dx = (threshold + rmsNoise) / n;
    double f0 = 0;
    for (int i = 1; i <= n; i++) {
      double f = getValue(dx * i);
      s += (f0 + f) / 2;
      f0 = f;
    }
    return 1 - s * dx * getNorm();
  }

getTimeAxis()

std::vector< double > getTimeAxis

(

)

const

Returns time axis (sample times)

Returns

vector of sample times

Definition at line 48 of file TOPSampleTimes.cc.

  {
    std::vector<double> timeAxis;
    for (unsigned i = 0; i < c_TimeAxisSize + 1; i++) {
      timeAxis.push_back(m_timeAxis[i]);
    }
    return timeAxis;
  }

getTimeBin()

double getTimeBin	(	int	window,
		int	sampleNumber ) const

Returns time bin of a given sample number and window (e.g.

time interval to the next sample)

Parameters

window	ASIC window number
sampleNumber	sample number counted from begin of the specified ASIC window

Returns

time bin in [ns]

Definition at line 104 of file TOPSampleTimes.cc.

  {
    int i = (window * c_WindowSize + sampleNumber) % c_TimeAxisSize;
    if (i < 0) i += c_TimeAxisSize;
    return m_timeAxis[i + 1] - m_timeAxis[i];
  }

getTimeWalk()

double getTimeWalk

(

int

pulseHeight

)

const

Returns time-walk at given pulse height.

Parameters

pulseHeight

pulse height [ADC counts]

Returns

time-walk [ns]

Definition at line 27 of file TOPCalTimeWalk.cc.

  {
    /* time walk parameters are coefficients in a polynomial
     * that is evaluated at x = pulseHeight here */
 
    double f = 0;
    double x = 1;
    for (auto p : m_timeWalkParams) {
      f += p * x;
      x *= pulseHeight;
    }
    return f;
  }

getTransformation() [1/2]

const Transform3D & getTransformation

(

int

moduleID

)

const

Returns transformation from local to nominal frame.

Transformation is: rotation first then translation.

Parameters

moduleID

module ID (1-based)

Returns

transformation

Definition at line 298 of file TOPCalModuleAlignment.cc.

  {
    if (m_transforms.empty()) setTransformations();
    unsigned module = moduleID - 1;
    if (module >= c_numModules) {
      B2ERROR("Invalid module number, returning identity transformation (" << ClassName() << ")");
      return m_transforms[c_numModules];
    }
    return m_transforms[module];
  }

getTransformation() [2/2]

Transform3D getTransformation

(

)

const

Returns transformation from local to nominal frame.

Translation is always given in basf2 units.

Returns

transformation (rotation and then a translation)

Definition at line 37 of file TOPGeoModuleDisplacement.cc.

  {
    RotationX Rx(m_alpha);
    RotationY Ry(m_beta);
    RotationZ Rz(m_gamma);
    Translation3D t(m_x, m_y, m_z);
    return Transform3D(Rz * Ry * Rx, t);
  }

getTTS()

const TOPNominalTTS & getTTS

(

unsigned

type

)

const

Returns time transition spread of a given PMT type.

Parameters

type

PMT type

Returns

TTS of a given PMT type if found, otherwise nominal TTS

Definition at line 50 of file TOPGeometry.cc.

  {
    std::map<unsigned, TOPNominalTTS>::const_iterator it = m_tts.find(type);
    if (it == m_tts.end()) return m_nominalTTS;
    return it->second;
  }

getValue() [1/2]

double getValue

(

double

x

)

const

Returns a value of non-normalized pulse height distribution at x Multiply with getNorm() to get normalized one.

Parameters

x	argument [ADC counts]

Definition at line 19 of file TOPPulseHeightPar.cc.

  {
    if (x <= 0 or x0 == 0) return 0;
    double xx = x / x0;
    double f = std::pow(xx, p1) * std::exp(-std::pow(xx, p2));
    return f;
  }

getValue() [2/2]

double getValue

(

double

t

)

const

Returns value at time t of the normalized waveform using interpolator.

Parameters

t	time [ns] (t = 0 is at 50% CF leading edge crossing)

Returns

normalized value (normalization: peaking value = 1)

Definition at line 75 of file TOPSignalShape.cc.

  {
    if (m_shape.empty()) {
      B2ERROR("TOPSignalShape::getValue: object not initialized");
      return 0;
    }
    if (isnan(t)) return 0; // to prevent interpolator to crash with seg. fault
    if (t < m_tmin) return 0;
    if (t > m_tmax) return m_shape.back() * exp(-(t - m_tmax) / m_tau);
    if (!m_interpolator) {
      std::vector<double> shape(m_shape); // since argument in TSpline5 is not const!
      m_interpolator = new TSpline5("signalShape", m_tmin, m_tmax,
                                    shape.data(), shape.size());
    }
    return m_interpolator->Eval(t);
  }

getX() [1/2]

double getX

(

int

moduleID

)

const

Returns the shift x on a single module.

Parameters

moduleID

module ID (1-based)

Returns

x shift along the x direction

Definition at line 149 of file TOPCalModuleAlignment.cc.

  {
    unsigned module = moduleID - 1;
    if (module >= c_numModules) {
      B2WARNING("Invalid module number, returning 0 (" << ClassName() << ")");
      return 0;
    }
    return m_x[module];
  }

getX() [2/2]

double getX

(

unsigned

col

)

const

Returns x coordinate of pixel center.

Parameters

col	valid column number (1-based)

Returns

x coordinate of pixel center

Definition at line 47 of file TOPGeoPMT.cc.

  {
    return (int)(m_numColumns + 1 - 2 * col) / 2.0 * getDx();
  }

getXErr()

double getXErr

(

int

moduleID

)

const

Returns the error on x on a single module.

Parameters

moduleID

module ID (1-based)

Returns

errX error on the x shift

Definition at line 209 of file TOPCalModuleAlignment.cc.

  {
    unsigned module = moduleID - 1;
    if (module >= c_numModules) {
      B2WARNING("Invalid module number, returning 0 (" << ClassName() << ")");
      return 0;
    }
    return m_errX[module];
  }

getY() [1/2]

double getY

(

int

moduleID

)

const

Returns the shift y on a single module.

Parameters

moduleID

module ID (1-based)

Returns

y shift along the y direction

Definition at line 159 of file TOPCalModuleAlignment.cc.

  {
    unsigned module = moduleID - 1;
    if (module >= c_numModules) {
      B2WARNING("Invalid module number, returning 0 (" << ClassName() << ")");
      return 0;
    }
    return m_y[module];
  }

getY() [2/2]

double getY

(

unsigned

row

)

const

Returns y coordinate of pixel center.

Parameters

row	valid row number (1-based)

Returns

y coordinate of pixel center

Definition at line 52 of file TOPGeoPMT.cc.

  {
    return (int)(2 * row - m_numRows - 1) / 2.0 * getDy();
  }

getYErr()

double getYErr

(

int

moduleID

)

const

Returns the error on y on a single module.

Parameters

moduleID

module ID (1-based)

Returns

errY error on the y shift

Definition at line 219 of file TOPCalModuleAlignment.cc.

  {
    unsigned module = moduleID - 1;
    if (module >= c_numModules) {
      B2WARNING("Invalid module number, returning 0 (" << ClassName() << ")");
      return 0;
    }
    return m_errY[module];
  }

getZ()

double getZ

(

int

moduleID

)

const

Returns the shift z on a single module.

Parameters

moduleID

module ID (1-based)

Returns

z shift along the z direction

Definition at line 169 of file TOPCalModuleAlignment.cc.

  {
    unsigned module = moduleID - 1;
    if (module >= c_numModules) {
      B2WARNING("Invalid module number, returning 0 (" << ClassName() << ")");
      return 0;
    }
    return m_z[module];
  }

getZc()

double getZc

(

)

const

Returns spherical mirror center of curvature in z (in local frame of this segment)

Returns

center of curvature in z

Definition at line 21 of file TOPGeoMirrorSegment.cc.

  {
    double dx = fabs(getXc()) - getWidth() / 2;
    double dy = fabs(getYc()) - getThickness() / 2;
    double z = getOuterRadius();
    if (dx > 0) z = sqrt(z * z - dx * dx);
    if (dy > 0) z = sqrt(z * z - dy * dy);
 
    return getFullLength() / 2 - z;
  }

getZErr()

double getZErr

(

int

moduleID

)

const

Returns the error on z on a single module.

Parameters

moduleID

module ID (1-based)

Returns

errZ error on the z shift

Definition at line 229 of file TOPCalModuleAlignment.cc.

  {
    unsigned module = moduleID - 1;
    if (module >= c_numModules) {
      B2WARNING("Invalid module number, returning 0 (" << ClassName() << ")");
      return 0;
    }
    return m_errZ[module];
  }

importAsicShifts_BS13d()

void importAsicShifts_BS13d	(	double	s0,
		double	s1,
		double	s2,
		double	s3,
		int	expNo,
		int	firstRun,
		int	lastRun )

Import ASIC shifts of BS13d.

Parameters

s0	shift of carrier 0 [ns]
s1	shift of carrier 1 [ns]
s2	shift of carrier 2 [ns]
s3	shift of carrier 3 [ns]
expNo	experiment number of IOV
firstRun	first run number of IOV
lastRun	last run number of IOV

Definition at line 342 of file TOPDatabaseImporter.cc.

  {
 
    std::vector<double> shifts;
    shifts.push_back(s0);
    shifts.push_back(s1);
    shifts.push_back(s2);
    shifts.push_back(s3);
 
    DBImportObjPtr<TOPCalAsicShift> asicShift;
    asicShift.construct();
 
    int moduleID = 13;
    unsigned bs = 3;
    for (unsigned carrier = 0; carrier < 4; carrier++) {
      for (unsigned a = 0; a < 4; a++) {
        unsigned asic = a + carrier * 4 + bs * 16;
        asicShift->setT0(moduleID, asic, shifts[carrier]);
      }
    }
 
    IntervalOfValidity iov(expNo, firstRun, expNo, lastRun);
    asicShift.import(iov);
 
    B2INFO("ASIC shifts of BS13d imported for exp " << expNo << " run " << firstRun <<
           " to " << lastRun);
  }

importChannelMask()

void importChannelMask	(	std::string	fileName,
		int	expNo,
		int	firstRun,
		int	lastRun )

Import channel mask from a root file (given as 1D histograms, one per slot).

Parameters

fileName	root file name
expNo	experiment number of IOV
firstRun	first run number of IOV
lastRun	last run number of IOV

Definition at line 621 of file TOPDatabaseImporter.cc.

  {
    // declare db object to be imported -- and construct it
    DBImportObjPtr<TOPCalChannelMask> channelMask;
    channelMask.construct();
 
    // open the root file
    TFile* file = TFile::Open(fileName.c_str(), "r");
    if (!file) {
      B2ERROR("openFile: " << fileName << " *** failed to open");
      return;
    }
    B2INFO(fileName << ": open for reading");
 
    // loop over slots and set channel mask
    int nModules = TOPGeometryPar::Instance()->getGeometry()->getNumModules();
    int active = 0, dead = 0, noisy = 0;
    for (int moduleID = 1; moduleID <= nModules; moduleID++) {
      std::string name = "slot_" + std::to_string(moduleID);
      auto* h = (TH1F*) file->Get(name.c_str());
      if (!h) {
        B2ERROR("Histogram with name '" + name + "' not found");
        continue;
      }
      for (int channel = 0; channel < h->GetNbinsX(); channel++) {
        int value = h->GetBinContent(channel + 1);
        if (value == 0) {
          channelMask->setActive(moduleID, channel);
          active++;
        } else if (value == 1) {
          channelMask->setDead(moduleID, channel);
          dead++;
        } else {
          channelMask->setNoisy(moduleID, channel);
          noisy++;
        }
      }
    }
    file->Close();
 
    // import to database
    IntervalOfValidity iov(expNo, firstRun, expNo, lastRun);
    channelMask.import(iov);
 
    B2INFO("Channel mask for exp " << expNo << " run " << firstRun << " to " << lastRun
           << " imported. Active channels: " << active << ", dead: " << dead
           << ", noisy: " << noisy);
 
  }

importChannelT0()

void importChannelT0	(	std::string	fileName,
		int	expNo,
		int	firstRun,
		int	lastRun )

Import channel T0 calibration constants The input is a root file with 1D histograms (one per slot, named "channelT0_slot*")

Parameters

fileName	root file name
expNo	experiment number of IOV
firstRun	first run number of IOV
lastRun	last run number of IOV

Definition at line 290 of file TOPDatabaseImporter.cc.

  {
    // declare db object to be imported -- and construct it
    DBImportObjPtr<TOPCalChannelT0> channelT0;
    channelT0.construct();
 
    // open the root file
    TFile* file = TFile::Open(fileName.c_str(), "r");
    if (!file) {
      B2ERROR("openFile: " << fileName << " *** failed to open");
      return;
    }
    B2INFO(fileName << ": open for reading");
 
    // loop over slots and set channel T0
    int nModules = TOPGeometryPar::Instance()->getGeometry()->getNumModules();
    int count = 0; // counter of calibrated constants
    for (int moduleID = 1; moduleID <= nModules; moduleID++) {
      std::string name = "channelT0_slot";
      if (moduleID < 10) name += "0";
      name += std::to_string(moduleID);
      auto* h = (TH1F*) file->Get(name.c_str());
      if (!h) {
        B2ERROR("Histogram with name '" + name + "' not found");
        continue;
      }
      for (int channel = 0; channel < h->GetNbinsX(); channel++) {
        double value = h->GetBinContent(channel + 1);
        double error = h->GetBinError(channel + 1);
        channelT0->setT0(moduleID, channel, value, error);
        if (error > 0) {
          count++;
        } else {
          channelT0->setUnusable(moduleID, channel);
        }
      }
    }
    file->Close();
 
    channelT0->suppressAverage();
 
    // import to database
    IntervalOfValidity iov(expNo, firstRun, expNo, lastRun);
    channelT0.import(iov);
 
    B2INFO("Channel T0 for exp " << expNo << " run " << firstRun << " to " << lastRun
           << " imported. Calibrated channels: " << count << "/" << nModules * 512);
 
  }

importCommonT0()

void importCommonT0	(	double	value,
		double	error,
		int	expNo,
		int	firstRun,
		int	lastRun,
		bool	roughlyCalibrated = false )

Import common T0 calibration constants.

Parameters

value	central value of T0
error	uncertainty on T0
expNo	experiment number of IOV
firstRun	first run number of IOV
lastRun	last run number of IOV
roughlyCalibrated	if true set payload status to roughly calibrated

Definition at line 446 of file TOPDatabaseImporter.cc.

  {
    DBImportObjPtr<TOPCalCommonT0> commonT0;
    commonT0.construct(value, error);
    if (roughlyCalibrated) commonT0->setRoughlyCalibrated();
 
    IntervalOfValidity iov(expNo, firstRun, expNo, lastRun);
    commonT0.import(iov);
 
    B2INFO("--> constants for exp = " << expNo
           << " run = " << firstRun << " to " << lastRun << " imported");
  }

importDummyCalAsicShift()

void importDummyCalAsicShift	(	int	firstExp = 0,
		int	firstRun = 0,
		int	lastExp = -1,
		int	lastRun = -1 )

import a dummy payload of TOPCalAsicShift DB objects

Parameters

firstExp	first experiment number of IOV
firstRun	first run number of IOV
lastExp	first experiment number of IOV
lastRun	last run number of IOV

Definition at line 1395 of file TOPDatabaseImporter.cc.

  {
    IntervalOfValidity iov(firstExp, firstRun, lastExp, lastRun);
    DBImportObjPtr<TOPCalAsicShift> asicShift;
    asicShift.construct();
    asicShift.import(iov);
    B2INFO("Dummy TOPCalAsicShift imported");
  }

importDummyCalChannelNoise()

void importDummyCalChannelNoise	(	int	firstExp = 0,
		int	firstRun = 0,
		int	lastExp = -1,
		int	lastRun = -1 )

import a dummy payload of TOPCalChannelNoise DB objects

Parameters

firstExp	first experiment number of IOV
firstRun	first run number of IOV
lastExp	first experiment number of IOV
lastRun	last run number of IOV

Definition at line 1318 of file TOPDatabaseImporter.cc.

  {
    IntervalOfValidity iov(firstExp, firstRun, lastExp, lastRun);
    DBImportObjPtr<TOPCalChannelNoise> channelNoise;
    channelNoise.construct();
    channelNoise.import(iov);
    B2INFO("Dummy TOPCalChannelNoise imported");
    return;
  }

importDummyCalChannelPulseHeight()

void importDummyCalChannelPulseHeight	(	int	firstExp = 0,
		int	firstRun = 0,
		int	lastExp = -1,
		int	lastRun = -1 )

import a dummy payload of TOPCalChannelPulseHeight DB objects

Parameters

firstExp	first experiment number of IOV
firstRun	first run number of IOV
lastExp	first experiment number of IOV
lastRun	last run number of IOV

Definition at line 1329 of file TOPDatabaseImporter.cc.

  {
    IntervalOfValidity iov(firstExp, firstRun, lastExp, lastRun);
    DBImportObjPtr<TOPCalChannelPulseHeight> pulseHeight;
    pulseHeight.construct();
    pulseHeight.import(iov);
    B2INFO("Dummy TOPCalChannelPulseHeight imported");
    return;
  }

importDummyCalChannelRQE()

void importDummyCalChannelRQE	(	int	firstExp = 0,
		int	firstRun = 0,
		int	lastExp = -1,
		int	lastRun = -1 )

import a dummy payload of TOPCalChannelRQE DB objects

Parameters

firstExp	first experiment number of IOV
firstRun	first run number of IOV
lastExp	first experiment number of IOV
lastRun	last run number of IOV

Definition at line 1340 of file TOPDatabaseImporter.cc.

  {
    IntervalOfValidity iov(firstExp, firstRun, lastExp, lastRun);
    DBImportObjPtr<TOPCalChannelRQE> channelRQE;
    channelRQE.construct();
    channelRQE.import(iov);
    B2INFO("Dummy TOPCalChannelRQE imported");
    return;
  }

importDummyCalChannelT0()

void importDummyCalChannelT0	(	int	firstExp = 0,
		int	firstRun = 0,
		int	lastExp = -1,
		int	lastRun = -1 )

import a dummy payload of TOPCalChannelT0 DB objects

Parameters

firstExp	first experiment number of IOV
firstRun	first run number of IOV
lastExp	first experiment number of IOV
lastRun	last run number of IOV

Definition at line 1296 of file TOPDatabaseImporter.cc.

  {
    IntervalOfValidity iov(firstExp, firstRun, lastExp, lastRun);
    DBImportObjPtr<TOPCalChannelT0> channelT0;
    channelT0.construct();
    channelT0.import(iov);
    B2INFO("Dummy TOPCalChannelT0 imported");
    return;
  }

importDummyCalChannelThreshold()

void importDummyCalChannelThreshold	(	int	firstExp = 0,
		int	firstRun = 0,
		int	lastExp = -1,
		int	lastRun = -1 )

import a dummy payload of TOPCalChannelThreshold DB objects

Parameters

firstExp	first experiment number of IOV
firstRun	first run number of IOV
lastExp	first experiment number of IOV
lastRun	last run number of IOV

Definition at line 1362 of file TOPDatabaseImporter.cc.

  {
    IntervalOfValidity iov(firstExp, firstRun, lastExp, lastRun);
    DBImportObjPtr<TOPCalChannelThreshold> channelThreshold;
    channelThreshold.construct();
    channelThreshold.import(iov);
    B2INFO("Dummy TOPCalChannelThreshold imported");
    return;
  }

importDummyCalChannelThresholdEff()

void importDummyCalChannelThresholdEff	(	int	firstExp = 0,
		int	firstRun = 0,
		int	lastExp = -1,
		int	lastRun = -1 )

import a dummy payload of TOPCalChannelThresholdEff DB objects

Parameters

firstExp	first experiment number of IOV
firstRun	first run number of IOV
lastExp	first experiment number of IOV
lastRun	last run number of IOV

Definition at line 1351 of file TOPDatabaseImporter.cc.

  {
    IntervalOfValidity iov(firstExp, firstRun, lastExp, lastRun);
    DBImportObjPtr<TOPCalChannelThresholdEff> channelThresholdEff;
    channelThresholdEff.construct();
    channelThresholdEff.import(iov);
    B2INFO("Dummy TOPCalChannelThresholdEff imported");
    return;
  }

importDummyCalCommonT0()

void importDummyCalCommonT0	(	int	firstExp = 0,
		int	firstRun = 0,
		int	lastExp = -1,
		int	lastRun = -1 )

import a dummy payload of TOPCalCommonT0 DB objects

Parameters

firstExp	first experiment number of IOV
firstRun	first run number of IOV
lastExp	first experiment number of IOV
lastRun	last run number of IOV

Definition at line 1373 of file TOPDatabaseImporter.cc.

  {
    IntervalOfValidity iov(firstExp, firstRun, lastExp, lastRun);
    DBImportObjPtr<TOPCalCommonT0> commonT0;
    commonT0.construct();
    commonT0.import(iov);
    B2INFO("Dummy TOPCalCommonT0 imported");
    return;
  }

importDummyCalIntegratedCharge()

void importDummyCalIntegratedCharge	(	int	firstExp = 0,
		int	firstRun = 0,
		int	lastExp = -1,
		int	lastRun = -1 )

import a dummy payload of TOPCalIntegratedCharge DB objects

Parameters

firstExp	first experiment number of IOV
firstRun	first run number of IOV
lastExp	first experiment number of IOV
lastRun	last run number of IOV

Definition at line 1384 of file TOPDatabaseImporter.cc.

  {
    IntervalOfValidity iov(firstExp, firstRun, lastExp, lastRun);
    DBImportObjPtr<TOPCalIntegratedCharge> integratedCharge;
    integratedCharge.construct();
    integratedCharge.import(iov);
    B2INFO("Dummy TOPCalIntegratedCharge imported");
    return;
  }

importDummyCalModuleAlignment()

void importDummyCalModuleAlignment	(	int	firstExp = 0,
		int	firstRun = 0,
		int	lastExp = -1,
		int	lastRun = -1 )

import a dummy payload of TOPCalModuleAlignment DB objects

Parameters

firstExp	first experiment number of IOV
firstRun	first run number of IOV
lastExp	first experiment number of IOV
lastRun	last run number of IOV

Definition at line 1274 of file TOPDatabaseImporter.cc.

  {
    IntervalOfValidity iov(firstExp, firstRun, lastExp, lastRun);
    DBImportObjPtr<TOPCalModuleAlignment> moduleAlignment;
    moduleAlignment.construct();
    moduleAlignment.import(iov);
    B2INFO("Dummy TOPCalModuleAlignment imported");
    return;
  }

importDummyCalModuleT0()

void importDummyCalModuleT0	(	int	firstExp = 0,
		int	firstRun = 0,
		int	lastExp = -1,
		int	lastRun = -1 )

import a dummy payload of TOPCalCalModuleT0 DB objects

Parameters

firstExp	first experiment number of IOV
firstRun	first run number of IOV
lastExp	first experiment number of IOV
lastRun	last run number of IOV

Definition at line 1285 of file TOPDatabaseImporter.cc.

  {
    IntervalOfValidity iov(firstExp, firstRun, lastExp, lastRun);
    DBImportObjPtr<TOPCalModuleT0> moduleT0;
    moduleT0.construct();
    moduleT0.import(iov);
    B2INFO("Dummy TOPCalModuleT0 imported");
    return;
  }

importDummyCalTimebase()

void importDummyCalTimebase	(	int	firstExp = 0,
		int	firstRun = 0,
		int	lastExp = -1,
		int	lastRun = -1 )

import a dummy payload of TOPCalTimebase DB objects

Parameters

firstExp	first experiment number of IOV
firstRun	first run number of IOV
lastExp	first experiment number of IOV
lastRun	last run number of IOV

Definition at line 1307 of file TOPDatabaseImporter.cc.

  {
    IntervalOfValidity iov(firstExp, firstRun, lastExp, lastRun);
    DBImportObjPtr<TOPCalTimebase> timebase;
    timebase.construct();
    timebase.import(iov);
    B2INFO("Dummy TOPCalTimebase imported");
    return;
  }

importFrontEndSettings()

void importFrontEndSettings	(	int	lookback,
		int	readoutWin,
		int	extraWin,
		int	offset,
		int	expNo,
		int	firstRun,
		int	lastRun )

Import front-end settings.

Parameters

lookback	the number of lookback windows
readoutWin	the number of readout windows
extraWin	the number of extra windows btw. lookback and readout window
offset	offset to photon peak [RF clocks]
expNo	experiment number of IOV
firstRun	first run number of IOV
lastRun	last run number of IOV

Definition at line 1242 of file TOPDatabaseImporter.cc.

  {
    DBImportObjPtr<TOPFrontEndSetting> feSetting;
    feSetting.construct();
 
    // write-window depths (write-window is 128 samples)
    std::vector<int> writeDepths;
    for (int i = 0; i < 3; i++) {
      writeDepths.push_back(214);
      writeDepths.push_back(212);
      writeDepths.push_back(214);
    }
    feSetting->setWriteDepths(writeDepths);
    feSetting->setLookbackWindows(lookback);
    feSetting->setReadoutWindows(readoutWin);
    feSetting->setExtraWindows(extraWin);
    feSetting->setOffset(offset);
 
    // window shifts
    std::vector<int> shifts = {0, 0, 1, 1, 1, 2};
    feSetting->setWindowShifts(shifts);
 
    IntervalOfValidity iov(expNo, firstRun, expNo, lastRun);
    feSetting.import(iov);
 
    B2INFO("Front-end settings imported for exp " << expNo << " run " << firstRun <<
           " to " << lastRun);
  }

importLocalT0Calibration()

void importLocalT0Calibration	(	std::string	fileNames,
		int	firstExp = 0,
		int	firstRun = 0,
		int	lastExp = -1,
		int	lastRun = -1 )

Import channel-by-channel T0 calibration constants to database The input is the root file with ntuple produced by TOPLaserCalibrator.

Parameters

fileNames	file names separated by space (since vector doesn't work!)
firstExp	first experiment number of IOV
firstRun	first run number of IOV
lastExp	first experiment number of IOV
lastRun	last run number of IOV

Definition at line 195 of file TOPDatabaseImporter.cc.

  {
    vector<string> fileNames;
    stringstream ss(fNames);
    string fName;
    while (ss >> fName) {
      fileNames.push_back(fName);
    }
 
    const auto* geo = TOPGeometryPar::Instance()->getGeometry();
 
    DBImportObjPtr<TOPCalChannelT0> channelT0;
    channelT0.construct();
 
    int nCal[16] = {0}; // number of calibrated channels per slot
 
    for (const auto& fileName : fileNames) {
      TFile* file = TFile::Open(fileName.c_str(), "r");
      B2INFO("--> Opening constants file " << fileName);
 
      if (!file) {
        B2ERROR("openFile: " << fileName << " *** failed to open");
        continue;
      }
      B2INFO("--> " << fileName << ": open for reading");
 
      TTree* treeCal = (TTree*)file->Get("chT0");
 
      if (!treeCal) {
        B2ERROR("openFile: no tree named chT0 found in " << fileName);
        file->Close();
        continue;
      }
 
      double t0Cal = 0.;
      double t0CalErr = 0.;
      int channelID = 0; // 0-511
      int slotID = 0;    // 1-16
      int fitStatus = 0; // fit status, 0 = OK
 
      treeCal->SetBranchAddress("channel", &channelID);
      treeCal->SetBranchAddress("slot", &slotID);
      treeCal->SetBranchAddress("channelT0", &t0Cal);
      treeCal->SetBranchAddress("channelT0Err", &t0CalErr);
      treeCal->SetBranchAddress("fitStatus", &fitStatus);
 
      B2INFO("--> importing constats");
 
      for (int iCal = 0; iCal < treeCal->GetEntries(); iCal++) {
        treeCal->GetEntry(iCal);
        if (!geo->isModuleIDValid(slotID)) {
          B2ERROR("Slot ID is not valid (fileName = " << fileName
                  << ", SlotID = " << slotID << ", ChannelID = " << channelID <<
                  "). Skipping the entry.");
          continue;
        }
        if (channelID < 0 or channelID > 511) {
          B2ERROR("Channel ID is not valid (fileName = " << fileName
                  << ", SlotID = " << slotID << ", ChannelID = " << channelID <<
                  "). Skipping the entry.");
          continue;
        }
        channelT0->setT0(slotID, channelID, t0Cal, t0CalErr);
        if (fitStatus == 0) {
          nCal[slotID - 1]++;
        } else {
          channelT0->setUnusable(slotID, channelID);
        }
      }
 
      file->Close();
      B2INFO("--> Input file closed");
    }
    channelT0->suppressAverage();
 
    IntervalOfValidity iov(firstExp, firstRun, lastExp, lastRun);
    channelT0.import(iov);
 
    short nCalTot = 0;
    B2INFO("Summary: ");
    for (int iSlot = 1; iSlot < 17; iSlot++) {
      B2INFO("--> Number of calibrated channels on Slot " << iSlot << " : " << nCal[iSlot - 1] << "/512");
      B2INFO("--> Cal on ch 1, 256 and 511:    " << channelT0->getT0(iSlot, 0) << ", " << channelT0->getT0(iSlot,
             257) << ", " << channelT0->getT0(iSlot, 511));
      nCalTot += nCal[iSlot - 1];
    }
 
 
    B2RESULT("Channel T0 calibration constants imported to database, calibrated channels: "
             << nCalTot << "/ 8192");
  }

importModuleT0()

void importModuleT0	(	std::string	fileName,
		int	expNo,
		int	firstRun,
		int	lastRun )

Import module T0 calibration constants The input is a root file with 1D histogram (name is "moduleT0")

Parameters

fileName	root file name
expNo	experiment number of IOV
firstRun	first run number of IOV
lastRun	last run number of IOV

Definition at line 514 of file TOPDatabaseImporter.cc.

  {
 
    // construct DB import object
    DBImportObjPtr<TOPCalModuleT0> moduleT0;
    moduleT0.construct();
 
    // open the root file
    TFile* file = TFile::Open(fileName.c_str(), "r");
    if (!file) {
      B2ERROR("openFile: " << fileName << " *** failed to open");
      return;
    }
    B2INFO(fileName << ": open for reading");
 
    // get histogram and set the DB import object
    auto* h = (TH1F*) file->Get("moduleT0");
    if (not h) {
      B2ERROR("no histogram 'moduleT0' found in the file, nothing imported");
      return;
    }
    int count = 0; // counter of calibrated
    for (int slot = 1; slot <= h->GetNbinsX(); slot++) {
      double value = h->GetBinContent(slot);
      double error = h->GetBinError(slot);
      moduleT0->setT0(slot, value, error);
      if (error > 0) {
        count++;
      } else {
        moduleT0->setUnusable(slot);
      }
    }
    file->Close();
 
    moduleT0->suppressAverage();
 
    // import the object
    IntervalOfValidity iov(expNo, firstRun, expNo, lastRun);
    moduleT0.import(iov);
 
    B2INFO("Module T0 for exp " << expNo << " run " << firstRun << " to " << lastRun
           << " imported. Calibrated modules: " << count << "/" << 16);
 
  }

importModuleT0Calibration()

void importModuleT0Calibration	(	std::string	fileName,
		int	firstExp = 0,
		int	firstRun = 0,
		int	lastExp = -1,
		int	lastRun = -1 )

Import module T0 calibration constants to database The input is the text file.

Parameters

fileName	name of the dat file with constants of all modules
firstExp	first experiment number of IOV
firstRun	first run number of IOV
lastExp	first experiment number of IOV
lastRun	last run number of IOV

Definition at line 461 of file TOPDatabaseImporter.cc.

  {
 
 
    DBImportObjPtr<TOPCalModuleT0> moduleT0;
    moduleT0.construct();
 
 
    ifstream inFile(fileName);
    B2INFO("--> Opening constants file " << fileName);
 
    if (!inFile) {
      B2ERROR("openFile: " << fileName << " *** failed to open");
      return;
    }
    B2INFO("--> " << fileName << ": open for reading");
 
 
    B2INFO("--> importing constants");
 
    while (!inFile.eof()) {
      int slot = 0;
      int dummy = 0;
      double T0 = 0;
      double T0_err = 0;
 
      inFile >> slot >> dummy >> T0 >> T0_err;
      if (slot < 1 or slot > 16) {
        B2ERROR("Module ID is not valid. Skipping the entry.");
        continue;
      }
      moduleT0->setT0(slot, T0, T0_err);
 
    }
    inFile.close();
    B2INFO("--> Input file closed");
 
    moduleT0->suppressAverage();
 
    IntervalOfValidity iov(firstExp, firstRun, lastExp, lastRun);
    moduleT0.import(iov);
 
    B2INFO("Summary: ");
    for (int iSlot = 1; iSlot < 17; iSlot++) {
      B2INFO("--> Time offset of Slot " << iSlot << " = " << moduleT0->getT0(iSlot));
    }
 
 
  }

importOfflineCommonT0Calibration()

void importOfflineCommonT0Calibration	(	std::string	fileName,
		int	firstExp = 0,
		int	firstRun = 0,
		int	lastExp = -1,
		int	lastRun = -1 )

Import common T0 calibration constants derived form the offline data reprocessing to database The input is a root file containing a tree (one per run).

In the future we may add a furter implementation that reads the IOV from the root file itself.

Parameters

fileName	name of the root file with constants of all modules
firstExp	first experiment number of IOV
firstRun	first run number of IOV
lastExp	first experiment number of IOV
lastRun	last run number of IOV

Definition at line 372 of file TOPDatabaseImporter.cc.

  {
    TFile* file = TFile::Open(fileName.c_str(), "r");
    B2INFO("--> Opening constants file " << fileName);
 
    if (!file) {
      B2ERROR("openFile: " << fileName << " *** failed to open");
      return;
    }
    B2INFO("--> " << fileName << ": open for reading");
 
    TTree* treeCal = (TTree*)file->Get("tree");
 
    if (!treeCal) {
      B2ERROR("openFile: no tree named tree found in " << fileName);
      file->Close();
      return;
    }
 
    float t0 = 0.;
    float t0Err = 0;
    float chi2 = 0;
    float integral = 0;
    float sigma = 0;
    int runNum = 0;
    int fitStatus = 0;
 
    treeCal->SetBranchAddress("offset", &t0);
    treeCal->SetBranchAddress("runNum", &runNum);
    treeCal->SetBranchAddress("sigma", &sigma);
    treeCal->SetBranchAddress("offsetErr", &t0Err);
    treeCal->SetBranchAddress("chi2", &chi2);
    treeCal->SetBranchAddress("integral", &integral);
    treeCal->SetBranchAddress("fitStatus", &fitStatus);
 
    treeCal->GetEntry(0);
 
    if (lastRun == -1 and firstRun == -1) {
      lastRun = runNum;
      firstRun = runNum;
      B2INFO("Using the run number from the tree ");
    } else {
      B2INFO("Using the run numbers passed to the importer");
    }
    B2INFO("IOV = (" << firstExp << ", "  << firstRun << ", "
           << lastExp << ", "  << lastRun << ")");
 
    DBImportObjPtr<TOPCalCommonT0> commonT0;
    commonT0.construct(t0, t0Err);
 
    if (fitStatus == 0 and integral > 10 and sigma > 0.05 and sigma < 0.33) {
      B2INFO("Good calibration found ");
      B2INFO("t0 = " << t0 << " +- "  << t0Err);
      B2INFO("sigma = " << sigma);
      B2INFO("chi2 = " << chi2);
    } else {
      B2INFO("BAD calibration found - set calibration to 'unusable'");
      B2INFO("t0 = " << t0 << " +- "  << t0Err);
      B2INFO("sigma = " << sigma);
      B2INFO("chi2 = " << chi2);
      B2INFO("fit status = " << fitStatus);
      commonT0->setUnusable();
    }
 
    file->Close();
 
    IntervalOfValidity iov(firstExp, firstRun, lastExp, lastRun);
    commonT0.import(iov);
    B2INFO("--> constants imported");
    B2INFO("   ");
  }

importPmtGainData()

void importPmtGainData	(	std::string	fileName,
		std::string	treeName = "gainPmtData",
		int	firstExp = 0,
		int	firstRun = 0,
		int	lastExp = -1,
		int	lastRun = -1 )

Import PMT gain parameters data to database.

Parameters

fileName	: name of the root file containing relevant data
treeName	: name of the tree containing relevant data
firstExp	first experiment number of IOV
firstRun	first run number of IOV
lastExp	first experiment number of IOV
lastRun	last run number of IOV

Definition at line 796 of file TOPDatabaseImporter.cc.

  {
 
    // declare db objects to be imported
    DBImportArray<TOPPmtGainPar> pmtGains;
 
    static const int nChann = 16;
    std::string* serialNum = 0;
    float gain_const[nChann], gain_slope[nChann], gain_ratio[nChann];
    float hv_op0, hv_op;
 
    // open root file and get tree
    TFile* file = TFile::Open(fileName.c_str(), "r");
    if (!file) {
      B2ERROR("Cannot open the file " << fileName);
      return;
    }
    TTree* tGainData = (TTree*)file->Get(treeName.c_str());
    if (!tGainData) {
      B2ERROR("No TTree with name " << treeName << " in file " << fileName);
      file->Close();
      return;
    }
 
    tGainData->SetBranchAddress("serialNum", &serialNum);
    tGainData->SetBranchAddress("gain_const", &gain_const);
    tGainData->SetBranchAddress("gain_slope", &gain_slope);
    tGainData->SetBranchAddress("gain_ratio", &gain_ratio);
    tGainData->SetBranchAddress("hv_op0", &hv_op0);
    tGainData->SetBranchAddress("hv_op", &hv_op);
 
 
    // loop on input tree entries and construct the pmtGain objects
    int countPMTs = 0;
 
    for (int ient = 0; ient < tGainData->GetEntries(); ient++) {
      tGainData->GetEntry(ient);
      auto* pmtGain = pmtGains.appendNew(*serialNum);
 
      for (int ic = 0; ic < nChann; ic++) {
        pmtGain->setPmtPixelData(ic + 1, gain_const[ic], gain_slope[ic], gain_ratio[ic]);
        pmtGain->setNominalHV0(-fabs(hv_op0));
        pmtGain->setNominalHV(-fabs(hv_op));
      }
      countPMTs++;
    }
    file->Close();
 
    IntervalOfValidity iov(firstExp, firstRun, lastExp, lastRun);
    pmtGains.import(iov);
 
    B2RESULT("PMT gain data imported to database for " << countPMTs << " PMT's.");
 
    return;
  }

importPmtInstallationData()

void importPmtInstallationData	(	std::string	fileName,
		std::string	treeName = "installationPmtData",
		int	firstExp = 0,
		int	firstRun = 0,
		int	lastExp = -1,
		int	lastRun = -1 )

Import PMT installation data to database.

Parameters

fileName	: name of the root file containing relevant data
treeName	: name of the tree containing relevant data
firstExp	first experiment number of IOV
firstRun	first run number of IOV
lastExp	first experiment number of IOV
lastRun	last run number of IOV

Definition at line 855 of file TOPDatabaseImporter.cc.

  {
 
    // declare db objects to be imported
    DBImportArray<TOPPmtInstallation> pmtInst;
 
    std::string* serialNum = 0;
    int moduleCNum, slotNum, arrayNum, PMTposition;
    TOPPmtObsoleteData::EType type;
 
    // open root file and get tree
    TFile* file = TFile::Open(fileName.c_str(), "r");
    if (!file) {
      B2ERROR("Cannot open the file " << fileName);
      return;
    }
    TTree* tInstData = (TTree*)file->Get(treeName.c_str());
    if (!tInstData) {
      B2ERROR("No TTree with name " << treeName << " in file " << fileName);
      file->Close();
      return;
    }
 
    tInstData->SetBranchAddress("serialNum", &serialNum);
    tInstData->SetBranchAddress("moduleCNum", &moduleCNum);
    tInstData->SetBranchAddress("slotNum", &slotNum);
    tInstData->SetBranchAddress("arrayNum", &arrayNum);
    tInstData->SetBranchAddress("PMTposition", &PMTposition);
    tInstData->SetBranchAddress("type", &type);
 
    // loop on input tree entries and construct the pmtInstallation objects
    int countPMTs = 0;
 
    for (int ient = 0; ient < tInstData->GetEntries(); ient++) {
      tInstData->GetEntry(ient);
      pmtInst.appendNew(*serialNum, moduleCNum, slotNum, arrayNum, PMTposition, type);
      countPMTs++;
    }
    file->Close();
 
    IntervalOfValidity iov(firstExp, firstRun, lastExp, lastRun);
    pmtInst.import(iov);
 
    B2RESULT("PMT installation data imported to database for " << countPMTs << " PMT's.");
 
  }

importPmtObsoleteData()

void importPmtObsoleteData	(	std::string	fileName,
		std::string	treeName = "obsPmtData",
		int	firstExp = 0,
		int	firstRun = 0,
		int	lastExp = -1,
		int	lastRun = -1 )

Import PMT specifications from Hamamatsu (not to be used!)

Parameters

fileName	: name of the root file containing relevant data
treeName	: name of the tree containing relevant data
firstExp	first experiment number of IOV
firstRun	first run number of IOV
lastExp	first experiment number of IOV
lastRun	last run number of IOV

Definition at line 905 of file TOPDatabaseImporter.cc.

  {
 
    // declare db objects to be imported
    DBImportArray<TOPPmtObsoleteData> pmtObsData;
 
    std::string* serialNum = 0;
    std::string* cathode = 0;
    float hv_spec, dark_spec, qe380_spec;
    TOPPmtObsoleteData::EType type;
 
    // open root file and get tree
    TFile* file = TFile::Open(fileName.c_str(), "r");
    if (!file) {
      B2ERROR("Cannot open the file " << fileName);
      return;
    }
    TTree* tObsData = (TTree*)file->Get(treeName.c_str());
    if (!tObsData) {
      B2ERROR("No TTree with name " << treeName << " in file " << fileName);
      file->Close();
      return;
    }
 
    tObsData->SetBranchAddress("serialNum", &serialNum);
    tObsData->SetBranchAddress("cathode", &cathode);
    tObsData->SetBranchAddress("hv_spec", &hv_spec);
    tObsData->SetBranchAddress("dark_spec", &dark_spec);
    tObsData->SetBranchAddress("qe380_spec", &qe380_spec);
    tObsData->SetBranchAddress("type", &type);
 
    // loop on input tree entries and construct the pmt obsolete data objects
    int countPMTs = 0;
 
    for (int ient = 0; ient < tObsData->GetEntries(); ient++) {
      tObsData->GetEntry(ient);
 
      // make sure the HV from specifications is negative
      hv_spec = -fabs(hv_spec);
 
      pmtObsData.appendNew(*serialNum, type, *cathode, hv_spec, dark_spec, qe380_spec);
      countPMTs++;
    }
 
    IntervalOfValidity iov(firstExp, firstRun, lastExp, lastRun);
    pmtObsData.import(iov);
 
    B2RESULT("PMT obsolete data imported to database for " << countPMTs << " PMT's.");
 
    file->Close();
 
    delete serialNum;
    delete cathode;
 
    return;
  }

importPmtPulseHeightFitResult()

void importPmtPulseHeightFitResult	(	std::string	fileName,
		int	firstExp = 0,
		int	firstRun = 0,
		int	lastExp = -1,
		int	lastRun = -1 )

Import fit results of pulse height disribution for channel gain and threshold efficiency.

Parameters

fileName	: name of the root file containing relevant data, which is obtained from TOPGainEfficiencyMonitor
firstExp	first experiment number of IOV
firstRun	first run number of IOV
lastExp	first experiment number of IOV
lastRun	last run number of IOV

Definition at line 1129 of file TOPDatabaseImporter.cc.

  {
    // declare db objects to be imported
    DBImportObjPtr<TOPCalChannelPulseHeight> calChannelPulseHeight;
    DBImportObjPtr<TOPCalChannelThresholdEff> calChannelThresholdEff;
    calChannelPulseHeight.construct();
    calChannelThresholdEff.construct();
 
    TFile* file = TFile::Open(fileName.c_str());
    if (!file) {
      B2ERROR("openFile: " << fileName << " *** failed to open");
      return;
    }
    TTree* tr = (TTree*)file->Get("tree");   // defined in TOPGainEfficiencyCalculatorModule
    if (!tr) {
      B2ERROR("No TTree with name tree found in " << fileName);
      file->Close();
      return;
    }
 
    short slotId = 0;
    short pixelId = 0;
    float p1 = -1;
    float p2 = -1;
    float x0 = -1;
    float threshold = -1;
    float efficiency = -1;
    float chisquare = -1;
    int ndf = 0;
    tr->SetBranchAddress("slotId", &slotId);
    tr->SetBranchAddress("pixelId", &pixelId);
    tr->SetBranchAddress("p1UseIntegral", &p1);
    tr->SetBranchAddress("p2UseIntegral", &p2);
    tr->SetBranchAddress("x0UseIntegral", &x0);
    tr->SetBranchAddress("thresholdForIntegral", &threshold);
    tr->SetBranchAddress("efficiencyUseIntegral", &efficiency);
    tr->SetBranchAddress("chisquareUseIntegral", &chisquare);
    tr->SetBranchAddress("ndfUseIntegral", &ndf);
 
    const auto& channelMapper = TOPGeometryPar::Instance()->getChannelMapper();
    if (!channelMapper.isValid()) {
      B2ERROR("No valid channel mapper found");
      file->Close();
      return;
    }
 
    long nEntries = tr->GetEntries();
    std::map<short, float> reducedChisqMap;
    for (long iEntry = 0 ; iEntry < nEntries ; iEntry++) {
      tr->GetEntry(iEntry);
 
      if (efficiency < 0) continue;
 
      if (!channelMapper.isPixelIDValid(pixelId)) {
        B2ERROR("invalid pixelID" << pixelId);
        continue;
      }
      auto channel = channelMapper.getChannel(pixelId);
      short globalChannelNumber = slotId * 1000 + channel;
      float redChisq = chisquare / ndf;
 
      //in case entries for the same channel appears multiple time, use data with smaller reduced chisquare
      //(This can happen when distribution is fit manually and results are appended for channels with fit failure)
      if (reducedChisqMap.count(globalChannelNumber) == 0
          or reducedChisqMap[globalChannelNumber] > redChisq) {
        reducedChisqMap[globalChannelNumber] = redChisq;
        calChannelPulseHeight->setParameters(slotId, channel, x0, p1, p2);
        calChannelThresholdEff->setThrEff(slotId, channel, efficiency, (short)threshold);
 
        if (redChisq > 10.) {
          calChannelPulseHeight->setUnusable(slotId, channel);
          calChannelThresholdEff->setUnusable(slotId, channel);
        }
      }
    }
    file->Close();
 
    IntervalOfValidity iov(firstExp, firstRun, lastExp, lastRun);
    calChannelPulseHeight.import(iov);
    calChannelThresholdEff.import(iov);
 
    B2RESULT("Imported channel-by-channel gain and efficiency data from fitting of pulse height distribution for "
             << reducedChisqMap.size() << " channels from " << fileName << " file.");
 
    return;
  }

importPmtQEData()

void importPmtQEData	(	std::string	fileName,
		std::string	treeName = "qePmtData",
		int	firstExp = 0,
		int	firstRun = 0,
		int	lastExp = -1,
		int	lastRun = -1 )

Import PMT Quantum Efficiency data to database.

Parameters

fileName	: name of the root file containing relevant data
treeName	: name of the tree containing relevant data
firstExp	first experiment number of IOV
firstRun	first run number of IOV
lastExp	first experiment number of IOV
lastRun	last run number of IOV

Definition at line 723 of file TOPDatabaseImporter.cc.

  {
 
    // declare db objects to be imported
    DBImportArray<TOPPmtQE> pmtQEs;
 
    static const int nChann = 16;
    std::string* serialNum = 0;
    std::vector<float>* QE_data[nChann];
    float lambdaFirst, lambdaStep, collEff0, collEff;
 
    TBranch* bQE_data[nChann];
 
    // open root file and get tree
    TFile* file = TFile::Open(fileName.c_str(), "r");
    if (!file) {
      B2ERROR("Cannot open the file " << fileName);
      return;
    }
    TTree* tQeData = (TTree*)file->Get(treeName.c_str());
    if (!tQeData) {
      B2ERROR("No TTree with name " << treeName << " in file " << fileName);
      file->Close();
      return;
    }
 
    tQeData->SetBranchAddress("serialNum", &serialNum);
    tQeData->SetBranchAddress("lambdaFirst", &lambdaFirst);
    tQeData->SetBranchAddress("lambdaStep", &lambdaStep);
    tQeData->SetBranchAddress("collEff0", &collEff0);
    tQeData->SetBranchAddress("collEff", &collEff);
 
    for (int ic = 0; ic < nChann; ic++) {
      // must initialize vectors and branches
      QE_data[ic] = new std::vector<float>;
      bQE_data[ic] = new TBranch();
 
      TString cString = "QE_ch";
      cString += ic + 1;
      tQeData->SetBranchAddress(cString, &QE_data[ic], &bQE_data[ic]);
    }
 
    // loop on input tree entries and construct the pmtQE objects
    int countPMTs = 0;
 
    for (int ient = 0; ient < tQeData->GetEntries(); ient++) {
 
      tQeData->GetEntry(ient);
 
      auto* pmtQE = pmtQEs.appendNew(*serialNum, lambdaFirst, lambdaStep, collEff0, collEff);
 
      for (int ic = 0; ic < nChann; ic++) {
        int tEntry = tQeData->LoadTree(ient);
        bQE_data[ic]->GetEntry(tEntry);
 
        pmtQE->setQE(ic + 1, *QE_data[ic]);
      }   // end loop on channels
 
      countPMTs++;
    }
    file->Close();
 
    IntervalOfValidity iov(firstExp, firstRun, lastExp, lastRun);
    pmtQEs.import(iov);
 
    B2RESULT("PMT QE data imported to database for " << countPMTs << " PMT's.");
 
    return;
  }

importPmtTTSHisto()

void importPmtTTSHisto	(	std::string	fileName,
		std::string	treeName = "ttsPmtHisto",
		int	firstExp = 0,
		int	firstRun = 0,
		int	lastExp = -1,
		int	lastRun = -1 )

Import histograms used for PMT TTS determination.

Parameters

fileName	: name of the root file containing relevant data
treeName	: name of the tree containing relevant data
firstExp	first experiment number of IOV
firstRun	first run number of IOV
lastExp	first experiment number of IOV
lastRun	last run number of IOV

Definition at line 1066 of file TOPDatabaseImporter.cc.

  {
 
    // declare db objects to be imported
    DBImportArray<TOPPmtTTSHisto> pmtTtsHistos;
 
    static const int nChann = 16;
    std::string* serialNum = 0;
    float hv = 0;
    TH1F* histo[nChann] = {0};
 
    // open root file and get tree
    TFile* file = TFile::Open(fileName.c_str(), "r");
    if (!file) {
      B2ERROR("Cannot open the file " << fileName);
      return;
    }
    TTree* tTtsHisto = (TTree*)file->Get(treeName.c_str());
    if (!tTtsHisto) {
      B2ERROR("No TTree with name " << treeName << " in file " << fileName);
      file->Close();
      return;
    }
 
    tTtsHisto->SetBranchAddress("serialNum", &serialNum);
    tTtsHisto->SetBranchAddress("hv", &hv);
    for (int ic = 0; ic < nChann; ic++) {
      TString hString = "hist_ch";
      hString += ic + 1;
      tTtsHisto->SetBranchAddress(hString, &histo[ic]);
    }
 
    // loop on input tree entries and construct the pmt tts histo objects
    int countHists = 0;
 
    for (int ient = 0; ient < tTtsHisto->GetEntries(); ient++) {
 
      tTtsHisto->GetEntry(ient);
 
      // make sure the HV used in the test is negative
      hv = -fabs(hv);
 
      B2INFO("Saving TTS histograms for PMT " << *serialNum << ", HV = " << hv);
 
      auto* pmtTtsHisto = pmtTtsHistos.appendNew(*serialNum, hv);
      for (int ic = 0; ic < nChann; ic++) {
        pmtTtsHisto->setHistogram(ic + 1, histo[ic]);
      }
      countHists++;
    }
    file->Close();
 
    IntervalOfValidity iov(firstExp, firstRun, lastExp, lastRun);
    pmtTtsHistos.import(iov);
 
    B2RESULT("Imported " << countHists << " sets of TTS histograms from " << fileName << " file.");
 
    return;
  }

importPmtTTSPar()

void importPmtTTSPar	(	std::string	fileName,
		std::string	treeName = "ttsPmtPar",
		int	firstExp = 0,
		int	firstRun = 0,
		int	lastExp = -1,
		int	lastRun = -1 )

Import gaussians fitting the TTS distributions.

Parameters

fileName	: name of the root file containing relevant data
treeName	: name of the tree containing relevant data
firstExp	first experiment number of IOV
firstRun	first run number of IOV
lastExp	first experiment number of IOV
lastRun	last run number of IOV

Definition at line 965 of file TOPDatabaseImporter.cc.

  {
 
    // declare db objects to be imported
    DBImportArray<TOPPmtTTSPar> pmtTtsPars;
 
    static const int nChann = 16;
    std::string* serialNum  = 0;
    std::vector<float>* gausFrac[nChann];
    std::vector<float>* gausMean[nChann];
    std::vector<float>* gausSigma[nChann];
 
    TBranch* bGFrac[nChann];
    TBranch* bGMean[nChann];
    TBranch* bGSigma[nChann];
 
 
    // open root file and get tree
    TFile* file = TFile::Open(fileName.c_str(), "r");
    if (!file) {
      B2ERROR("Cannot open the file " << fileName);
      return;
    }
    TTree* tTtsPar = (TTree*)file->Get(treeName.c_str());
    if (!tTtsPar) {
      B2ERROR("No TTree with name " << treeName << " in file " << fileName);
      file->Close();
      return;
    }
 
    tTtsPar->SetBranchAddress("serialNum", &serialNum);
    for (int ic = 0; ic < nChann; ic++) {
      // must initialize vectors and branches
      gausFrac[ic] = new std::vector<float>;
      gausMean[ic] = new std::vector<float>;
      gausSigma[ic] = new std::vector<float>;
 
      bGFrac[ic] = new TBranch();
      bGMean[ic] = new TBranch();
      bGSigma[ic] = new TBranch();
 
 
      TString cStringF = "gausFrac_ch";
      TString cStringM = "gausMean_ch";
      TString cStringS = "gausSigma_ch";
 
      cStringF += ic + 1;
      cStringM += ic + 1;
      cStringS += ic + 1;
 
      tTtsPar->SetBranchAddress(cStringF, &gausFrac[ic], &bGFrac[ic]);
      tTtsPar->SetBranchAddress(cStringM, &gausMean[ic], &bGMean[ic]);
      tTtsPar->SetBranchAddress(cStringS, &gausSigma[ic], &bGSigma[ic]);
    }
 
    // loop on input tree entries and construct the pmt tts par objects
    int countPMTs = 0;
 
    for (int ient = 0; ient < tTtsPar->GetEntries(); ient++) {
 
      tTtsPar->GetEntry(ient);
 
      auto* pmtTtsPar = pmtTtsPars.appendNew(*serialNum);
 
      for (int ic = 0; ic < nChann; ic++) {
 
        int tEntry = tTtsPar->LoadTree(ient);
        bGFrac[ic]->GetEntry(tEntry);
        bGMean[ic]->GetEntry(tEntry);
        bGSigma[ic]->GetEntry(tEntry);
 
        // check that the vectors have the same size. Otherwise skip the channel
        if ((gausFrac[ic]->size() != gausMean[ic]->size()) ||
            (gausFrac[ic]->size() != gausSigma[ic]->size())) {
 
          B2ERROR("The TTSPar vectors for PMT " << serialNum << ", channel " << ic + 1 << " have different sizes! Skipping channel...");
          continue;
        }
 
        for (uint iv = 0; iv < gausFrac[ic]->size(); iv++) {
          pmtTtsPar->appendGaussian(ic + 1,
                                    gausFrac[ic]->at(iv),
                                    gausMean[ic]->at(iv),
                                    gausSigma[ic]->at(iv));
        }
      }
      countPMTs++;
    }
    file->Close();
 
    IntervalOfValidity iov(firstExp, firstRun, lastExp, lastRun);
    pmtTtsPars.import(iov);
 
    B2RESULT("PMT TTS parameters imported to database for " << countPMTs << " PMT's.");
 
    return;
  }

importSampleTimeCalibration()

void importSampleTimeCalibration	(	std::string	fileNames,
		int	firstExp = 0,
		int	firstRun = 0,
		int	lastExp = -1,
		int	lastRun = -1 )

Import sample time calibration constants to database Output of TOPTimeBaseCalibrator (root files with histograms)

Parameters

fileNames	file names separated by space (since vector doesn't work!)
firstExp	first experiment number of IOV
firstRun	first run number of IOV
lastExp	first experiment number of IOV
lastRun	last run number of IOV

Definition at line 69 of file TOPDatabaseImporter.cc.

  {
 
    // make vector out of files separated with space
 
    vector<string> fileNames;
    stringstream ss(fNames);
    string fName;
    while (ss >> fName) {
      fileNames.push_back(fName);
    }
 
    // prepare what is needed
 
    const auto* geo = TOPGeometryPar::Instance()->getGeometry();
    auto syncTimeBase = geo->getNominalTDC().getSyncTimeBase();
 
    DBImportObjPtr<TOPCalTimebase> timeBase;
    timeBase.construct(syncTimeBase);
 
    set<int> scrodIDs;
 
    // read constants from files and put them to DB object
 
    for (const auto& fileName : fileNames) {
      TFile* file = TFile::Open(fileName.c_str(), "r");
      if (!file) {
        B2ERROR("openFile: " << fileName << " *** failed to open");
        continue;
      }
      B2INFO(fileName << ": open for reading");
 
      TH1F* hsuccess = (TH1F*) file->Get("success");
      if (!hsuccess) {
        B2ERROR("Fit status histogram '" << hsuccess << "' not found");
        file->Close();
        continue;
      }
 
      int goodChannels = 0;
      int numChannels = hsuccess->GetNbinsX();
      for (int channel = 0; channel < numChannels; channel++) {
        if (hsuccess->GetBinContent(channel + 1) == 0) continue;
 
        string hname = "sampleTimes_ch" +  to_string(channel);
 
        TH1F* hsampleTimes = (TH1F*) file->Get(hname.c_str());
        if (!hsampleTimes) {
          B2ERROR("Histogram '" << hname << "' with calibration constants not found");
          continue;
        }
        // parse scrodID from histogram title
        string title = hsampleTimes->GetTitle();
        auto iscrod = title.find("scrod");
        auto ichannel = title.find("channel");
        if (iscrod == string::npos or ichannel == string::npos) {
          B2ERROR("Unsuccessful parsing of scrodID from '" << title << "'");
          continue;
        }
        iscrod += 5;
        int len = ichannel - iscrod;
        if (len < 1) {
          B2ERROR("Unsuccessful parsing of scrodID from '" << title << "'");
          continue;
        }
        int scrodID = stoi(title.substr(iscrod, len));
        scrodIDs.insert(scrodID);
 
        double rescale = 1;
        if (hsampleTimes->GetBinContent(257) > 0)
          rescale = 2 * syncTimeBase / hsampleTimes->GetBinContent(257);
 
        vector<double> sampleTimes;
        for (int isamp = 0; isamp < 256; isamp++) {
          sampleTimes.push_back(hsampleTimes->GetBinContent(isamp + 1) * rescale);
        }
        goodChannels++;
 
        timeBase->append(scrodID, channel, sampleTimes);
      }
 
      file->Close();
      B2INFO("--> number of calibrated channels: " << goodChannels);
      B2INFO("file closed");
    }
 
    // set calibration for missing ones, using previous calibrated channel within asic
 
    B2INFO("set constants for uncalibrated channels using nearest calibrated channel within an ASIC");
    for (auto scrodID : scrodIDs) {
      int nasic = 128 / 8;
      for (int as = 0; as < nasic; as++) {
        const TOPSampleTimes* sampleTimes = 0;
        for (int ch = 0; ch < 15; ch++) {
          int channel = as * 8 + (ch % 8);
          if (timeBase->isAvailable(scrodID, channel)) {
            sampleTimes = timeBase->getSampleTimes(scrodID, channel);
          } else if (sampleTimes) {
            timeBase->append(scrodID, channel, sampleTimes->getTimeAxis());
          }
        }
        if (!sampleTimes) {
          B2INFO("No calibration available for ASIC " << as << " of scrodID " << scrodID);
        }
      }
    }
 
    // import constants
    IntervalOfValidity iov(firstExp, firstRun, lastExp, lastRun);
    timeBase.import(iov);
 
    // final message
 
    int nall = timeBase->getSampleTimes().size();
    int ncal = 0;
    for (const auto& sampleTimes : timeBase->getSampleTimes()) {
      if (sampleTimes.isCalibrated()) ncal++;
    }
 
    B2RESULT("Sample time calibration constants imported to database, calibrated channels: "
             << ncal << "/" << nall);
  }

importTest() [1/2]

void importTest

(

)

for testing purposes only!

• will be removed ...

Definition at line 1506 of file TOPDatabaseImporter.cc.

  {
 
    DBImportArray<TOPPmtGainPar> pmtGains;
 
    auto* pmtGain = pmtGains.appendNew("JT00123");
    pmtGain->setNominalHV(3520);
    for (unsigned channel = 1; channel <= 16; channel++) {
      pmtGain->setPmtPixelData(channel, -13.77, 0.0042, 0.4);
    }
 
    pmtGain = pmtGains.appendNew("JT02135");
    pmtGain->setNominalHV(3450);
    for (unsigned channel = 1; channel <= 16; channel++) {
      pmtGain->setPmtPixelData(channel, -12.77, 0.0045, 0.4);
    }
 
    for (const auto& gain : pmtGains) gain.print();
 
    //  IntervalOfValidity iov(0, 0, -1, -1); // all experiments and runs
    //  pmtGains.import(iov);
 
 
  }

importTest() [2/2]

void importTest	(	int	runNumber,
		double	syncTimeBase )

for testing purposes only!

• will be removed ...

Definition at line 1473 of file TOPDatabaseImporter.cc.

  {
 
    DBImportObjPtr<TOPCalTimebase> timeBase;
    vector<double> timeAxis;
    for (int i = 0; i < 256; i++) {
      timeAxis.push_back(syncTimeBase / 128.0 * i);
    }
 
 
    timeBase.construct(syncTimeBase);
    for (unsigned scrodID = 0; scrodID < 64; scrodID++) {
      for (unsigned channel = 0; channel < 128; channel++) {
        timeBase->append(scrodID, channel, timeAxis);
      }
    }
 
    if (runNumber == 3) {
      timeBase.addEventDependency(10);
      timeBase.construct(syncTimeBase + 100);
      for (unsigned scrodID = 0; scrodID < 64; scrodID++) {
        for (unsigned channel = 0; channel < 128; channel++) {
          timeBase->append(scrodID, channel, timeAxis);
        }
      }
    }
 
    IntervalOfValidity iov(1, runNumber, 1, runNumber);
    timeBase.import(iov);
 
  }

importTimeWalk()

void importTimeWalk	(	PyObject *	list,
		double	a,
		double	b,
		int	firstExp = 0,
		int	firstRun = 0,
		int	lastExp = -1,
		int	lastRun = -1 )

payload TOPCalTimeWalk import parameters for time-walk correction and electronic time resolution tuning

Parameters

list	Python list of parameters of time-walk calibration curve [ns]
a	electronic time resolution: noise term excess parameter [ns]
b	electronic time resolution: quadratic term parameter [ns]
firstExp	first experiment number of IOV
firstRun	first run number of IOV
lastExp	first experiment number of IOV
lastRun	last run number of IOV

Definition at line 1442 of file TOPDatabaseImporter.cc.

  {
 
    std::vector<double> params;
    if (PyList_Check(list)) {
      for (Py_ssize_t i = 0; i < PyList_Size(list); i++) {
        PyObject* value = PyList_GetItem(list, i);
        params.push_back(PyFloat_AsDouble(value));
        B2INFO(i << " " << params.back());
      }
    } else {
      B2ERROR("Input Python object is not a list");
      return;
    }
 
    DBImportObjPtr<TOPCalTimeWalk> timeWalk;
    timeWalk.construct();
    timeWalk->set(params, a, b);
 
    IntervalOfValidity iov(firstExp, firstRun, lastExp, lastRun);
    timeWalk.import(iov);
 
    B2RESULT("Time-walk constants imported");
  }

interpolate()

double interpolate ( double lambda, const std::vector< float > & QE ) const

private

Interpolate between QE datapoints (linear interpolation).

Parameters

lambda	wavelength in [nm]
QE	quantum efficiency data points

Returns

quantum efficiency at lambda

Definition at line 105 of file TOPPmtQE.cc.

  {
    double dlam = lambda - m_lambdaFirst;
    if (dlam < 0 or dlam > (QE.size() - 1) * m_lambdaStep) return 0;
    unsigned i = int(dlam / m_lambdaStep);
    if (i > QE.size() - 2) return QE.back();
    return QE[i] + (QE[i + 1] - QE[i]) / m_lambdaStep * (dlam - i * m_lambdaStep);
  }

isActive()

bool isActive	(	int	moduleID,
		unsigned	channel ) const

Returns false if the channel is dead or noisy, and true is the channel is active.

Parameters

moduleID	module ID (1-based)
channel	hardware channel number (0-based)

Definition at line 72 of file TOPCalChannelMask.cc.

  {
    int module = moduleID - 1;
    if (!check(module, channel)) {
      B2WARNING("Returning false");
      return false;
    }
    return (m_status[module][channel] == c_Active);
  }

isAvailable()

bool isAvailable	(	unsigned	scrodID,
		unsigned	channel ) const

Checks if calibration is available.

Parameters

scrodID	SCROD ID
channel	channel number (within SCROD or within module)

Returns

true if available

Definition at line 76 of file TOPCalTimebase.cc.

  {
    if (m_map.empty()) createMap();
 
    unsigned key = (scrodID << 16) + (channel % 128);
    Iterator it = m_map.find(key);
    return (it != m_map.end());
 
  }

isCalibrated() [1/4]

bool isCalibrated	(	int	moduleID,
		unsigned	asic ) const

Returns calibration status.

Parameters

moduleID	module ID (1-based)
asic	ASIC number within a module (0 - 63)

Returns

true, if good calibrated

Definition at line 66 of file TOPCalAsicShift.cc.

  {
    unsigned module = moduleID - 1;
    if (module >= c_numModules) return false;
    if (asic >= c_numAsics) return false;
    return m_status[module][asic] == c_Calibrated;
  }

isCalibrated() [2/4]

bool isCalibrated	(	int	moduleID,
		unsigned	channel ) const

Returns calibration status.

Parameters

moduleID	module ID (1-based)
channel	hardware channel number (0-based)

Returns

true, if good calibrated

Definition at line 103 of file TOPCalChannelT0.cc.

  {
    unsigned module = moduleID - 1;
    if (module >= c_numModules) return false;
    if (channel >= c_numChannels) return false;
    return m_status[module][channel] == c_Calibrated;
  }

isCalibrated() [3/4]

bool isCalibrated

(

int

moduleID

)

const

Returns calibration status.

Parameters

moduleID

module ID (1-based)

Returns

true, if good calibrated

Definition at line 239 of file TOPCalModuleAlignment.cc.

  {
    unsigned module = moduleID - 1;
    if (module >= c_numModules) return false ;
    return m_status[module] == c_Calibrated;
  }

isCalibrated() [4/4]

bool isCalibrated

(

int

moduleID

)

const

Returns calibration status.

Parameters

moduleID

module ID (1-based)

Returns

true, if good calibrated

Definition at line 85 of file TOPCalModuleT0.cc.

  {
    unsigned module = moduleID - 1;
    if (module >= c_numModules) return false;
    return m_status[module] == c_Calibrated;
  }

isConsistent() [1/18]

bool isConsistent ( ) const

overridevirtual

Check for consistency of data members.

Returns

true if values consistent (valid)

Reimplemented from TOPGeoBase.

Reimplemented in TOPGeoMirrorSegment, and TOPGeoPrism.

Definition at line 23 of file TOPGeoBarSegment.cc.

  {
    if (m_width <= 0) return false;
    if (m_thickness <= 0) return false;
    if (m_length <= 0) return false;
    if (m_material.empty()) return false;
    if (m_glueThickness <= 0) return false;
    if (m_glueMaterial.empty()) return false;
    if (m_surface.getName().empty() and !m_surface.hasProperties()) return false;
    if (m_sigmaAlpha < 0) return false;
    if (m_brokenFraction > 0 and m_brokenGlueMaterial.empty()) return false;
    return true;
  }

isConsistent() [2/18]

bool isConsistent ( ) const

overridevirtual

Check for consistency of data members.

Returns

true if values consistent (valid)

Reimplemented from TOPGeoBase.

Definition at line 20 of file TOPGeoColdPlate.cc.

  {
    if (m_baseThickness <= 0) return false;
    if (m_baseMaterial.empty()) return false;
    if (m_coolThickness <= 0) return false;
    if (m_coolWidth <= 0) return false;
    if (m_coolMaterial.empty()) return false;
    return true;
  }

isConsistent() [3/18]

bool isConsistent ( ) const

overridevirtual

Check for consistency of data members.

Returns

true if values consistent (valid)

Reimplemented from TOPGeoBase.

Definition at line 20 of file TOPGeoEndPlate.cc.

  {
    if (m_thickness <= 0) return false;
    if (m_height <= 0) return false;
    if (m_material.empty()) return false;
    return true;
  }

isConsistent() [4/18]

bool isConsistent ( ) const

overridevirtual

Check for consistency of data members.

Returns

true if values consistent (valid)

Reimplemented from TOPGeoBase.

Definition at line 20 of file TOPGeoFrontEnd.cc.

  {
    if (m_FBWidth <= 0) return false;
    if (m_FBHeight <= 0) return false;
    if (m_FBThickness <= 0) return false;
    if (m_FBGap < 0) return false;
    if (m_FBMaterial.empty()) return false;
 
    if (m_HVWidth <= 0) return false;
    if (m_HVLength <= 0) return false;
    if (m_HVThickness <= 0) return false;
    if (m_HVGap < 0) return false;
    if (m_HVMaterial.empty()) return false;
 
    if (m_BSWidth <= 0) return false;
    if (m_BSHeight <= 0) return false;
    if (m_BSLength <= 0) return false;
    if (m_BSGap < 0) return false;
    if (m_BSMaterial.empty()) return false;
    if (m_spacerWidth <= 0) return false;
    if (m_spacerMaterial.empty()) return false;
 
    return true;
  }

isConsistent() [5/18]

bool isConsistent ( ) const

overridevirtual

Check for consistency of data members.

Returns

true if values consistent (valid)

Reimplemented from TOPGeoBase.

Definition at line 46 of file TOPGeoHoneycombPanel.cc.

  {
    if (m_width <= 0) return false;
    if (m_length <= 0) return false;
    if (m_minThickness <= 0) return false;
    if (m_maxThickness < m_minThickness) return false;
    if (m_radius <= 0) return false;
    if (m_edgeWidth <= 0) return false;
    if (m_material.empty()) return false;
    if (m_edgeMaterial.empty()) return false;
    return true;
  }

isConsistent() [6/18]

bool isConsistent ( ) const

overridevirtual

Check for consistency of data members.

Returns

true if values consistent (valid)

Reimplemented from TOPGeoBase.

Definition at line 126 of file TOPGeometry.cc.

  {
    if (m_modules.empty()) return false;
    for (const auto& module : m_modules) {
      if (!module.isConsistent()) return false;
    }
    if (!m_frontEnd.isConsistent()) return false;
    if (!m_QBB.isConsistent()) return false;
    if (m_numBoardStacks == 0) return false;
    if (!m_nominalQE.isConsistent()) return false;
    if (!m_nominalTTS.isConsistent()) return false;
    if (!m_nominalTDC.isConsistent()) return false;
    if (!m_wavelengthFilter.isConsistent()) return false;
    return true;
  }

isConsistent() [7/18]

bool isConsistent ( ) const

overridevirtual

Check for consistency of data members.

Returns

true if values consistent (valid)

Reimplemented from TOPGeoBarSegment.

Definition at line 33 of file TOPGeoMirrorSegment.cc.

  {
    if (m_radius <= 0) return false;
    if (m_coatingThickness <= 0) return false;
    if (m_coatingMaterial.empty()) return false;
    if (m_coatingSurface.getName().empty() and !m_coatingSurface.hasProperties()) return false;
    if (!TOPGeoBarSegment::isConsistent()) return false;
    return true;
  }

isConsistent() [8/18]

bool isConsistent ( ) const

overridevirtual

Check for consistency of data members.

Returns

true if values consistent (valid)

Reimplemented from TOPGeoBase.

Definition at line 150 of file TOPGeoModule.cc.

  {
    if (m_moduleID <= 0) return false;
    if (!m_bar1.isConsistent()) return false;
    if (!m_bar2.isConsistent()) return false;
    if (!m_mirror.isConsistent()) return false;
    if (!m_prism.isConsistent()) return false;
    if (!m_pmtArray.isConsistent()) return false;
    return true;
  }

isConsistent() [9/18]

bool isConsistent ( ) const

overridevirtual

Check for consistency of data members.

Returns

true if values consistent (valid)

Reimplemented from TOPGeoBase.

Definition at line 58 of file TOPGeoPMT.cc.

  {
    if (m_sizeX <= 0) return false;
    if (m_sizeY <= 0) return false;
    if (m_sizeZ <= 0) return false;
    if (m_wallThickness <= 0) return false;
 
    if (m_sensSizeX <= 0) return false;
    if (m_sensSizeY <= 0) return false;
    if (m_sensThickness <= 0) return false;
    if (m_numColumns == 0) return false;
    if (m_numRows == 0) return false;
    if (m_sensMaterial.empty()) return false;
 
    if (m_winThickness <= 0) return false;
    if (m_winMaterial.empty()) return false;
 
    if (m_botThickness <= 0) return false;
    if (m_botMaterial.empty()) return false;
 
    if (m_reflEdgeWidth <= 0) return false;
    if (m_reflEdgeThickness <= 0) return false;
 
    return true;
  }

isConsistent() [10/18]

bool isConsistent ( ) const

overridevirtual

Check for consistency of data members.

Returns

true if values consistent (valid)

Reimplemented from TOPGeoBase.

Definition at line 93 of file TOPGeoPMTArray.cc.

  {
    if (m_numRows == 0) return false;
    if (m_numColumns == 0) return false;
    if (m_dx <= 0) return false;
    if (m_dy <= 0) return false;
    if (m_material.empty()) return false;
    if (!m_pmt.isConsistent()) return false;
    if (m_cookieThickness <= 0) return false;
    if (m_cookieMaterial.empty()) return false;
    if (m_filterThickness <= 0) return false;
    if (m_filterMaterial.empty()) return false;
    return true;
  }

isConsistent() [11/18]

bool isConsistent ( ) const

overridevirtual

Check for consistency of data members.

Returns

true if values consistent (valid)

Reimplemented from TOPGeoBarSegment.

Definition at line 36 of file TOPGeoPrism.cc.

  {
    if (m_exitThickness <= 0) return false;
    if (m_flatLength < 0) return false;
    if (m_width <= 0) return false;
    if (m_thickness <= 0) return false;
    if (m_length <= 0) return false;
    if (m_material.empty()) return false;
    if (m_surface.getName().empty() and !m_surface.hasProperties()) return false;
    if (m_sigmaAlpha < 0) return false;
    if (m_brokenFraction > 0 and m_brokenGlueMaterial.empty()) return false;
    if (!m_peelOffRegions.empty()) {
      if (m_peelOffSize <= 0) return false;
      if (m_peelOffThickness <= 0) return false;
      if (m_peelOffMaterial.empty()) return false;
      double halfSize = (getWidth() - getPeelOffSize()) / 2;
      for (const auto& region : m_peelOffRegions) {
        if (fabs(getPeelOffCenter(region)) > halfSize) return false;
      }
    }
    return true;
  }

isConsistent() [12/18]

bool isConsistent ( ) const

overridevirtual

Check for consistency of data members.

Returns

true if values consistent (valid)

Reimplemented from TOPGeoBase.

Definition at line 21 of file TOPGeoPrismEnclosure.cc.

  {
    if (m_length <= 0) return false;
    if (m_height <= 0) return false;
    if (m_bottomThickness <= 0) return false;
    if (m_sideThickness <= 0) return false;
    if (m_backThickness <= 0) return false;
    if (m_frontThickness <= 0) return false;
    if (m_extensionThickness <= 0) return false;
    if (m_material.empty()) return false;
    return true;
  }

isConsistent() [13/18]

bool isConsistent ( ) const

overridevirtual

Check for consistency of data members.

Returns

true if values consistent (valid)

Reimplemented from TOPGeoBase.

Definition at line 230 of file TOPGeoQBB.cc.

  {
    if (m_width <= 0) return false;
    if (m_length <= 0) return false;
    if (m_prismPosition <= 0) return false;
    if (m_material.empty()) return false;
    if (!m_innerPanel.isConsistent()) return false;
    if (!m_outerPanel.isConsistent()) return false;
    if (!m_sideRails.isConsistent()) return false;
    if (!m_prismEnclosure.isConsistent()) return false;
    if (!m_endPlate.isConsistent()) return false;
    if (!m_coldPlate.isConsistent()) return false;
    return true;
  }

isConsistent() [14/18]

bool isConsistent ( ) const

overridevirtual

Check for consistency of data members.

Returns

true if values consistent (valid)

Reimplemented from TOPGeoBase.

Definition at line 20 of file TOPGeoSideRails.cc.

  {
    if (m_thickness <= 0) return false;
    if (m_reducedThickness <= 0) return false;
    if (m_height <= 0) return false;
    if (m_material.empty()) return false;
    return true;
  }

isConsistent() [15/18]

bool isConsistent ( ) const

overridevirtual

Check for consistency of data members.

Returns

true if values consistent (valid)

Reimplemented from TOPGeoBase.

Definition at line 58 of file TOPNominalQE.cc.

  {
    if (m_lambdaFirst <= 0) return false;
    if (m_lambdaStep <= 0) return false;
    if (m_CE <= 0) return false;
    if (m_QE.empty()) return false;
    return true;
  }

isConsistent() [16/18]

bool isConsistent ( ) const

overridevirtual

Check for consistency of data members.

Returns

true if values consistent (valid)

Reimplemented from TOPGeoBase.

Definition at line 84 of file TOPNominalTDC.cc.

  {
    if (m_pileupTime < 0) return false;
    if (m_doubleHitResolution < 0) return false;
    if (m_timeJitter < 0) return false;
    if (m_efficiency <= 0 or m_efficiency > 1) return false;
    return true;
  }

isConsistent() [17/18]

bool isConsistent ( ) const

overridevirtual

Check for consistency of data members.

Returns

true if values consistent (valid)

Reimplemented from TOPGeoBase.

Definition at line 60 of file TOPNominalTTS.cc.

  {
    if (m_tts.empty()) return false;
    if (!m_normalized) return false;
    return true;
  }

isConsistent() [18/18]

bool isConsistent ( ) const

overridevirtual

Check for consistency of data members.

Returns

true if values consistent (valid)

Reimplemented from TOPGeoBase.

Definition at line 34 of file TOPWavelengthFilter.cc.

  {
    if (getName().empty()) return true; // old payload from DB (filter included in QE)
 
    if (m_lambdaFirst <= 0) return false;
    if (m_lambdaStep <= 0) return false;
    if (m_transmittances.empty()) return false;
    return true;
  }

isDefault() [1/4]

bool isDefault	(	int	moduleID,
		unsigned	asic ) const

Returns calibration status.

Parameters

moduleID	module ID (1-based)
asic	ASIC number within a module (0 - 63)

Returns

true, if default (not calibrated)

Definition at line 75 of file TOPCalAsicShift.cc.

  {
    unsigned module = moduleID - 1;
    if (module >= c_numModules) return false;
    if (asic >= c_numAsics) return false;
    return m_status[module][asic] == c_Default;
  }

isDefault() [2/4]

bool isDefault	(	int	moduleID,
		unsigned	channel ) const

Returns calibration status.

Parameters

moduleID	module ID (1-based)
channel	hardware channel number (0-based)

Returns

true, if default (not calibrated)

Definition at line 112 of file TOPCalChannelT0.cc.

  {
    unsigned module = moduleID - 1;
    if (module >= c_numModules) return false;
    if (channel >= c_numChannels) return false;
    return m_status[module][channel] == c_Default;
  }

isDefault() [3/4]

bool isDefault

(

int

moduleID

)

const

Returns calibration status.

Parameters

moduleID

module ID (1-based)

Returns

true, if default (not calibrated)

Definition at line 246 of file TOPCalModuleAlignment.cc.

  {
    unsigned module = moduleID - 1;
    if (module >= c_numModules) return false ;
    return m_status[module] == c_Default;
  }

isDefault() [4/4]

bool isDefault

(

int

moduleID

)

const

Returns calibration status.

Parameters

moduleID

module ID (1-based)

Returns

true, if default (not calibrated)

Definition at line 93 of file TOPCalModuleT0.cc.

  {
    unsigned module = moduleID - 1;
    if (module >= c_numModules) return false;
    return m_status[module] == c_Default;
  }

isModuleIDValid()

bool isModuleIDValid

(

int

moduleID

)

const

Checks if module exists in m_modules.

Parameters

moduleID

module ID (1-based)

Returns

true if exists

Definition at line 33 of file TOPGeometry.cc.

  {
    for (const auto& module : m_modules) {
      if (module.getModuleID() == moduleID) return true;
    }
    return false;
  }

isPMTDecoupled()

bool isPMTDecoupled

(

unsigned

pmtID

)

const

Checks if PMT is optically decoupled.

Returns

true if decoupled

Definition at line 84 of file TOPGeoPMTArray.cc.

  {
    for (const auto& pmt : m_decoupledPMTs) {
      if (pmtID == pmt) return true;
    }
    return false;
  }

isSampleValid()

bool isSampleValid

(

int

sample

)

const

Check for the validity of sample number.

Parameters

sample

sample number

Returns

true if sample number is valid

Definition at line 77 of file TOPNominalTDC.cc.

  {
    if (sample < 0) return false;
    if (sample >= (int) m_numWindows * c_WindowSize) return false;
    return true;
  }

isUnusable() [1/4]

bool isUnusable	(	int	moduleID,
		unsigned	asic ) const

Returns calibration status.

Parameters

moduleID	module ID (1-based)
asic	ASIC number within a module (0 - 63)

Returns

true, if bad calibrated

Definition at line 84 of file TOPCalAsicShift.cc.

  {
    unsigned module = moduleID - 1;
    if (module >= c_numModules) return false;
    if (asic >= c_numAsics) return false;
    return m_status[module][asic] == c_Unusable;
  }

isUnusable() [2/4]

bool isUnusable	(	int	moduleID,
		unsigned	channel ) const

Returns calibration status.

Parameters

moduleID	module ID (1-based)
channel	hardware channel number (0-based)

Returns

true, if bad calibrated

Definition at line 121 of file TOPCalChannelT0.cc.

  {
    unsigned module = moduleID - 1;
    if (module >= c_numModules) return false;
    if (channel >= c_numChannels) return false;
    return m_status[module][channel] == c_Unusable;
  }

isUnusable() [3/4]

bool isUnusable

(

int

moduleID

)

const

Returns calibration status.

Parameters

moduleID

module ID (1-based)

Returns

true, if bad calibrated

Definition at line 253 of file TOPCalModuleAlignment.cc.

  {
    unsigned module = moduleID - 1;
    if (module >= c_numModules) return false ;
    return m_status[module] == c_Unusable;
  }

isUnusable() [4/4]

bool isUnusable

(

int

moduleID

)

const

Returns calibration status.

Parameters

moduleID

module ID (1-based)

Returns

true, if bad calibrated

Definition at line 101 of file TOPCalModuleT0.cc.

  {
    unsigned module = moduleID - 1;
    if (module >= c_numModules) return false;
    return m_status[module] == c_Unusable;
  }

momentumGlobalToNominal()

XYZVector momentumGlobalToNominal

(

const ROOT::Math::XYZVector &

momentum

)

const

Transforms momentum vector from Belle II to module nominal frame.

Parameters

momentum

momentum vector in Belle II frame

Returns

momentum vector in module nominal frame

Definition at line 143 of file TOPGeoModule.cc.

  {
    if (not m_transformNominal) setTransformation();
    return m_transformNominal->ApplyInverse(momentum);
  }

momentumNominalToGlobal()

XYZVector momentumNominalToGlobal

(

const ROOT::Math::XYZVector &

momentum

)

const

Transforms momentum vector from module nominal frame to Belle II frame.

Parameters

momentum

momentum vector in module nominal frame

Returns

momentum vector in Belle II frame

Definition at line 131 of file TOPGeoModule.cc.

  {
    if (not m_transformNominal) setTransformation();
    return *m_transformNominal * momentum;
  }

momentumToGlobal()

XYZVector momentumToGlobal

(

const ROOT::Math::XYZVector &

momentum

)

const

Transforms momentum vector from module internal (= nominal & displaced) frame to Belle II frame.

Parameters

momentum

momentum vector in module internal frame

Returns

momentum vector in Belle II frame

Definition at line 107 of file TOPGeoModule.cc.

  {
    if (not m_transform) setTransformation();
    return *m_transform * momentum;
  }

momentumToLocal()

XYZVector momentumToLocal

(

const ROOT::Math::XYZVector &

momentum

)

const

Transforms momentum vector from Belle II to module internal (= nominal & displaced) frame.

Parameters

momentum

momentum vector in Belle II frame

Returns

momentum vector in module internal frame

Definition at line 119 of file TOPGeoModule.cc.

  {
    if (not m_transform) setTransformation();
    return m_transform->ApplyInverse(momentum);
  }

normalize()

double normalize

(

)

Normalize the distribution (fractions)

Returns

value used to normalize fractions

Definition at line 36 of file TOPNominalTTS.cc.

  {
    float sum = 0;
    for (const auto& tts : m_tts) sum += tts.fraction;
    if (sum == 0) return 0;
    for (auto& tts : m_tts) tts.fraction /= sum;
    m_normalized = true;
    return sum;
  }

operator=()

TOPGeoModule & operator=

(

const TOPGeoModule &

module

)

Assignment operator.

Definition at line 33 of file TOPGeoModule.cc.

  {
    if (this != &module) {
      TOPGeoBase::operator=(module);
      m_moduleID = module.getModuleID();
      m_radius = module.getRadius();
      m_phi = module.getPhi();
      m_backwardZ = module.getBackwardZ();
      m_moduleCNumber = module.getModuleCNumber();
      m_bar1 = module.getBarSegment1();
      m_bar2 = module.getBarSegment2();
      m_mirror = module.getMirrorSegment();
      m_prism = module.getPrism();
      m_pmtArray = module.getPMTArray();
      m_arrayDisplacement = module.getPMTArrayDisplacement();
      m_moduleDisplacement = module.getModuleDisplacement();
 
      if (m_transform) delete m_transform;
      if (m_transformNominal) delete m_transformNominal;
      m_transform = 0;
      m_transformNominal = 0;
    }
    return *this;
  }

pointGlobalToNominal()

XYZPoint pointGlobalToNominal

(

const ROOT::Math::XYZPoint &

point

)

const

Transforms 3D point from Belle II to module nominal frame.

Parameters

point

3D point in Belle II frame (basf2 units!)

Returns

3D point in module nominal frame (basf2 units!)

Definition at line 137 of file TOPGeoModule.cc.

  {
    if (not m_transformNominal) setTransformation();
    return m_transformNominal->ApplyInverse(point);
  }

pointNominalToGlobal()

XYZPoint pointNominalToGlobal

(

const ROOT::Math::XYZPoint &

point

)

const

Transforms 3D point from module nominal frame to Belle II frame.

Parameters

point

3D point in module nominal frame (basf2 units!)

Returns

3D point in Belle II frame (basf2 units!)

Definition at line 125 of file TOPGeoModule.cc.

  {
    if (not m_transformNominal) setTransformation();
    return *m_transformNominal * point;
  }

pointToGlobal()

XYZPoint pointToGlobal

(

const ROOT::Math::XYZPoint &

point

)

const

Transforms 3D point from module internal (= nominal & displaced) frame to Belle II frame.

Parameters

point

3D point in module internal frame (basf2 units!)

Returns

3D point in Belle II frame (basf2 units!)

Definition at line 101 of file TOPGeoModule.cc.

  {
    if (not m_transform) setTransformation();
    return *m_transform * point;
  }

pointToLocal()

XYZPoint pointToLocal

(

const ROOT::Math::XYZPoint &

point

)

const

Transforms 3D point from Belle II to module internal (= nominal & displaced) frame.

Parameters

point

3D point in Belle II frame (basf2 units!)

Returns

3D point in module internal frame (basf2 units!)

Definition at line 113 of file TOPGeoModule.cc.

  {
    if (not m_transform) setTransformation();
    return m_transform->ApplyInverse(point);
  }

print() [1/22]

void print ( const std::string & title = "Bar segment geometry parameters" ) const

overridevirtual

Print the content of the class.

Parameters

title

title to be printed

Reimplemented from TOPGeoBase.

Reimplemented in TOPGeoMirrorSegment, and TOPGeoPrism.

Definition at line 38 of file TOPGeoBarSegment.cc.

  {
    TOPGeoBase::print(title);
    cout << " Vendor: " << getVendor() << ", serial number: " << getSerialNumber() << endl;
    cout << " Dimensions: " << getWidth() << " X " << getThickness() << " X " << getLength()
         << " " << s_unitName << endl;
    cout << " Material: " << getMaterial() << endl;
    cout << " Glue: " << getGlueMaterial() << ", thickness = " << getGlueThickness()
         << " " << s_unitName;
    cout << ", broken fraction = " << getBrokenGlueFraction();
    if (getBrokenGlueFraction() > 0) {
      cout << ", angle = " << getBrokenGlueAngle() / Unit::deg << " deg";
      cout << ", material = " << getBrokenGlueMaterial();
    }
    cout << endl;
    printSurface(m_surface);
    cout << "  - sigmaAlpha: " << getSigmaAlpha() << endl;
  }

print() [2/22]

void print ( const std::string & title ) const

virtual

Print the content of the class.

Parameters

title

title to be printed

Reimplemented in TOPGeoBarSegment, TOPGeoColdPlate, TOPGeoEndPlate, TOPGeoFrontEnd, TOPGeoHoneycombPanel, TOPGeometry, TOPGeoMirrorSegment, TOPGeoModule, TOPGeoModuleDisplacement, TOPGeoPMT, TOPGeoPMTArray, TOPGeoPMTArrayDisplacement, TOPGeoPrism, TOPGeoPrismEnclosure, TOPGeoQBB, TOPGeoSideRails, TOPNominalQE, TOPNominalTDC, TOPNominalTTS, and TOPWavelengthFilter.

Definition at line 28 of file TOPGeoBase.cc.

  {
    cout << title << ":" << endl;
    cout << " name: " << m_name << endl;
  }

print() [3/22]

void print ( const std::string & title = "QBB cold plate geometry parameters" ) const

overridevirtual

Print the content of the class.

Parameters

title

title to be printed

Reimplemented from TOPGeoBase.

Definition at line 31 of file TOPGeoColdPlate.cc.

  {
    TOPGeoBase::print(title);
 
    cout << " base plate: thickness = " << getBaseThickness() << " " << s_unitName;
    cout << ", material = " << getBaseMaterial() << endl;
 
    cout << " cooling plate: thickness = " << getCoolThickness() << " " << s_unitName;
    cout << ", width = " << getCoolWidth() << " " << s_unitName;
    cout << ", material = " << getCoolMaterial() << endl;
 
  }

print() [4/22]

void print ( const std::string & title = "QBB forward end plate geometry parameters" ) const

overridevirtual

Print the content of the class.

Parameters

title

title to be printed

Reimplemented from TOPGeoBase.

Definition at line 29 of file TOPGeoEndPlate.cc.

  {
    TOPGeoBase::print(title);
 
    cout << " thickness = " << getThickness() << " " << s_unitName;
    cout << ", height = " << getHeight() << " " << s_unitName;
    cout << ", material = " << getMaterial() << endl;
 
  }

print() [5/22]

void print ( const std::string & title = "Front-end geometry parameters" ) const

overridevirtual

Print the content of the class.

Parameters

title

title to be printed

Reimplemented from TOPGeoBase.

Definition at line 46 of file TOPGeoFrontEnd.cc.

  {
    TOPGeoBase::printUnderlined(title);
 
    cout << " Front board: " << getFrontBoardWidth() << " X " <<
         getFrontBoardHeight() << " X " << getFrontBoardThickness() << " " << s_unitName;
    cout << ", gap = " << getFrontBoardGap() << " " << s_unitName;
    cout << ", y = " << getFrontBoardY() << " " << s_unitName;
    cout << ", material = " << getFrontBoardMaterial() << endl;
 
    cout << " HV board: " << getHVBoardWidth() << " X " <<
         getHVBoardLength() << " X " << getHVBoardThickness() << " " << s_unitName;
    cout << ", gap = " << getHVBoardGap() << " " << s_unitName;
    cout << ", y = " << getHVBoardY() << " " << s_unitName;
    cout << ", material = " << getHVBoardMaterial() << endl;
 
    cout << " Board stack: " << getBoardStackWidth() << " X " <<
         getBoardStackHeight() << " X " << getBoardStackLength() << " " << s_unitName;
    cout << ", gap = " << getBoardStackGap() << " " << s_unitName;
    cout << ", y = " << getBoardStackY() << " " << s_unitName;
    cout << ", material = " << getBoardStackMaterial() << endl;
 
    cout << "              spacers: width = " << getSpacerWidth() << " " << s_unitName;
    cout << ", material = " << getSpacerMaterial() << endl;
 
  }

print() [6/22]

void print ( const std::string & title = "QBB honeycomb panel geometry parameters" ) const

overridevirtual

Print the content of the class.

Parameters

title

title to be printed

Reimplemented from TOPGeoBase.

Definition at line 60 of file TOPGeoHoneycombPanel.cc.

  {
    TOPGeoBase::print(title);
 
    cout << " width = " << getWidth() << " " << s_unitName;
    cout << ", length = " << getLength() << " " << s_unitName;
    cout << ", thickness: min = " << getMinThickness() << " " << s_unitName;
    cout << ", max = " << getMaxThickness() << " " << s_unitName;
    cout << ", radius = " << getRadius() << " " << s_unitName << endl;
 
    cout << " edge width = " << getEdgeWidth() << " " << s_unitName;
    cout << ", y = " << getY() << " " << s_unitName;
    cout << ", N = " << m_N << endl;
 
    cout << " materials: body = " << getMaterial();
    cout << ", edge = " << getEdgeMaterial() << endl;
  }

print() [7/22]

void print ( const std::string & title = "TOP geometry parameters" ) const

overridevirtual

Print the content of the class.

Parameters

title

title to be printed

Reimplemented from TOPGeoBase.

Definition at line 142 of file TOPGeometry.cc.

  {
    cout << endl;
    cout << "Geometry parameters of TOP counter:" << endl;
    cout << "===================================" << endl;
    cout << " name: " << m_name << endl;
    cout << " number of modules: " << m_modules.size() << endl << endl;
 
    for (const auto& module : m_modules) {
      module.print();
      cout << endl;
    }
    m_frontEnd.print();
    cout << " Number of board stacks: " << m_numBoardStacks << endl;
    cout << endl;
    m_QBB.print();
    cout << endl;
    m_nominalQE.print();
    cout << endl;
    cout << "Photo-detection efficiency tuning factors" << endl;
    cout << "-----------------------------------------" << endl;
    if (m_tuneFactorsPDE.empty()) {
      cout << " Not available" << endl;
    } else {
      for (const auto& x : m_tuneFactorsPDE) {
        cout << " PMT type: " << x.first << "  factor = " << x.second << endl;
      }
    }
    cout << endl;
    m_nominalTTS.print();
    cout << endl;
    for (const auto& tts : m_tts) {
      tts.second.print("TTS distribution");
      cout << endl;
    }
    m_nominalTDC.print();
    cout << endl;
    m_wavelengthFilter.print();
    cout << endl;
 
  }

print() [8/22]

void print ( const std::string & title = "Mirror segment geometry parameters" ) const

overridevirtual

Print the content of the class.

Parameters

title

title to be printed

Reimplemented from TOPGeoBarSegment.

Definition at line 44 of file TOPGeoMirrorSegment.cc.

  {
    TOPGeoBarSegment::print(title);
    cout << " Radius of curvature: " << getRadius() << " " << s_unitName << endl;
    cout << " Center of curvature: (" << getXc() << ", " << getYc()
         << ") " << s_unitName << endl;
    cout << " Reflective coating: " << getCoatingMaterial() << ", thickness: " <<
         getCoatingThickness() << " " << s_unitName << endl;
    printSurface(m_coatingSurface);
 
  }

print() [9/22]

void print ( const std::string & title = "Module geometry parameters" ) const

overridevirtual

Print the content of the class.

Parameters

title

title to be printed

Reimplemented from TOPGeoBase.

Definition at line 162 of file TOPGeoModule.cc.

  {
    cout << "Slot " << getModuleID() << " geometry parameters:" << endl;
    cout << "---------------------------" << endl;
    cout << " name: " << m_name << endl;
    cout << " moduleID = " << getModuleID();
    cout << ", radius = " << getRadius() << " " << s_unitName;
    cout << ", phi = " << getPhi() / Unit::deg << " deg";
    cout << ", backward z = " << getBackwardZ() << " " << s_unitName;
    cout << ", construction number = " << getModuleCNumber() << endl;
    cout << endl;
 
    m_prism.print();
    cout << endl;
    m_bar2.print("Bar segment 2 (backward) geometry parameters");
    cout << endl;
    m_bar1.print("Bar segment 1 (forward) geometry parameters");
    cout << endl;
    m_mirror.print();
    cout << endl;
    m_pmtArray.print();
    cout << endl;
    m_arrayDisplacement.print();
    cout << endl;
    m_moduleDisplacement.print();
    cout << endl;
 
  }

print() [10/22]

void print ( const std::string & title = "Module displacement parameters" ) const

overridevirtual

Print the content of the class.

Parameters

title

title to be printed

Reimplemented from TOPGeoBase.

Definition at line 25 of file TOPGeoModuleDisplacement.cc.

  {
    TOPGeoBase::print(title);
    cout << " translation vector: (" << getX() << ", "
         << getY() << ", " << getZ() << ")" << " " << s_unitName << endl;
    cout << " rotation angles: ";
    cout << getAlpha() << " ";
    cout << getBeta() << " ";
    cout << getGamma() << " ";
    cout << endl;
  }

print() [11/22]

void print ( const std::string & title = "MCP-PMT geometry parameters" ) const

overridevirtual

Print the content of the class.

Parameters

title

title to be printed

Reimplemented from TOPGeoBase.

Definition at line 85 of file TOPGeoPMT.cc.

  {
    TOPGeoBase::print(title);
 
    cout << " outer dimensions: " << getSizeX() << " X " << getSizeY()
         << " X " << getSizeZ() << " " << s_unitName << endl;
    cout << " wall thickness:   " << getWallThickness() << " " << s_unitName << endl;
    cout << " window thickness: " << getWinThickness() << " " << s_unitName << endl;
    cout << " bottom thickness: " << getBotThickness() << " " << s_unitName << endl;
    cout << " sensitive area:   " <<  getSensSizeX() << " X " << getSensSizeY()
         << " " << s_unitName;
    cout << ", thickness: " << getSensThickness() << " " << s_unitName << endl;
    cout << " number of pixels: " << getNumColumns() << " X " << getNumRows() << endl;
    cout << " casing material:    " << getWallMaterial() << endl;
    cout << " inner material:     " << getFillMaterial() << endl;
    cout << " window material:    " << getWinMaterial() << endl;
    cout << " bottom material:    " << getBotMaterial() << endl;
    cout << " sensitive material: " << getSensMaterial() << endl;
    cout << " reflective edge: " << endl;
    cout << "  - width: " << getReflEdgeWidth() << " " << s_unitName;
    cout << ", thickness: " << getReflEdgeThickness() << " " << s_unitName << endl;
 
    printSurface(m_reflEdgeSurface);
  }

print() [12/22]

void print ( const std::string & title = "PMT array geometry parameters" ) const

overridevirtual

Print the content of the class.

Parameters

title

title to be printed

Reimplemented from TOPGeoBase.

Definition at line 108 of file TOPGeoPMTArray.cc.

  {
    TOPGeoBase::print(title);
    cout << " size: " << getNumColumns() << " X " << getNumRows();
    cout << ", dimensions: " << getSizeX() << " X " << getSizeY() << " X " << getSizeZ();
    cout << " " << s_unitName << endl;
    cout << " cell: " << getDx() << " X " << getDy() << " " << s_unitName << endl;
    cout << " material: " << getMaterial() << endl;
    if (m_cookieThickness > 0) { // new version of payload
      cout << " silicone cookies:  thickness = " << getCookieThickness() << " "
           << s_unitName << ", material = " << getCookieMaterial() << endl;
      cout << " wavelength filter: thickness = " << getFilterThickness() << " "
           << s_unitName << ", material = " << getFilterMaterial() << endl;
    }
    cout << " air gap (decoupled PMT's): " << getAirGap() << " " << s_unitName << endl;
    cout << " optically decoupled PMT's:";
    if (m_decoupledPMTs.empty()) {
      cout << " None";
    } else {
      for (const auto& pmt : m_decoupledPMTs) cout << " " << pmt;
    }
    cout << endl;
    cout << endl;
    m_pmt.print();
  }

print() [13/22]

void print ( const std::string & title = "PMT array displacement parameters" ) const

overridevirtual

Print the content of the class.

Parameters

title

title to be printed

Reimplemented from TOPGeoBase.

Definition at line 20 of file TOPGeoPMTArrayDisplacement.cc.

  {
    TOPGeoBase::print(title);
    cout << " translation vector: (" << getX() << ", " << getY() << ")"
         << " " << s_unitName << endl;
    cout << " rotation angle: " << m_alpha << endl;
  }

print() [14/22]

void print ( const std::string & title = "Prism geometry parameters" ) const

overridevirtual

Print the content of the class.

Parameters

title

title to be printed

Reimplemented from TOPGeoBarSegment.

Definition at line 60 of file TOPGeoPrism.cc.

  {
    TOPGeoBase::print(title);
    cout << " Vendor: " << getVendor() << ", serial number: " << getSerialNumber() << endl;
    cout << " Dimensions: " << getWidth() << " X " << getThickness() << " X " << getLength()
         << " " << s_unitName << endl;
    cout << " Exit window dimensions: " << getWidth() << " X "
         << getExitThickness() << " " << s_unitName << endl;
    cout << " Prism angle: " << getAngle() / Unit::deg << " deg";
    cout << ", flat surface length: " << getFlatLength() << " " << s_unitName << endl;
    cout << " Material: " << getMaterial() << endl;
    if (getFilterThickness() > 0) { // old payload
      cout << " Wavelenght filter: " << getFilterMaterial()
           << ", thickness: " << getFilterThickness() << " " << s_unitName << endl;
    }
    if (!m_peelOffRegions.empty()) {
      cout << " Peel-off cookie regions: ";
      cout << " size = " << getPeelOffSize() << " " << s_unitName;
      cout << " thickness = " << getPeelOffThickness() << " " << s_unitName;
      cout << " material = " << getPeelOffMaterial() << endl;
      for (const auto& region : m_peelOffRegions) {
        cout << "   ID = " << region.ID << ", fraction = " << region.fraction
             << ", angle = " << region.angle  / Unit::deg << " deg"
             << ", xc = " << getPeelOffCenter(region) << " " << s_unitName << endl;
      }
    } else {
      cout << " Peel-off regions: None" << endl;
    }
    printSurface(m_surface);
    cout << "  - sigmaAlpha: " << getSigmaAlpha() << endl;
 
  }

print() [15/22]

void print ( const std::string & title = "QBB prism enclosure geometry parameters" ) const

overridevirtual

Print the content of the class.

Parameters

title

title to be printed

Reimplemented from TOPGeoBase.

Definition at line 35 of file TOPGeoPrismEnclosure.cc.

  {
    TOPGeoBase::print(title);
 
    cout << " length = " << getLength() << " " << s_unitName;
    cout << ", height = " << getHeight() << " " << s_unitName;
    cout << ", angle = " << getAngle() / Unit::deg << " deg";
    cout << ", material = " << getMaterial() << endl;
 
    cout << " thicknesses: bottom = " << getBottomThickness() << " " << s_unitName;
    cout << ", side = " << getSideThickness() << " " << s_unitName;
    cout << ", back = " << getBackThickness() << " " << s_unitName;
    cout << ", front = " << getFrontThickness() << " " << s_unitName;
    cout << ", extension plate = " << getExtensionThickness() << " " << s_unitName << endl;
  }

print() [16/22]

void print ( const std::string & title = "QBB geometry parameters" ) const

overridevirtual

Print the content of the class.

Parameters

title

title to be printed

Reimplemented from TOPGeoBase.

Definition at line 246 of file TOPGeoQBB.cc.

  {
    TOPGeoBase::printUnderlined(title);
 
    cout << " width = " << getWidth() << " " << s_unitName;
    cout << ", length = " << getLength() << " " << s_unitName;
    cout << ", prism-bar joint = " << getPrismPosition() << " " << s_unitName;
    cout << ", inside material = " << getMaterial() << endl;
 
    m_innerPanel.print("QBB inner honeycomb panel geometry parameters");
    m_outerPanel.print("QBB outer honeycomb panel geometry parameters");
    m_sideRails.print();
    m_prismEnclosure.print();
    m_endPlate.print();
    m_coldPlate.print();
 
  }

print() [17/22]

void print ( const std::string & title = "QBB side rails geometry parameters" ) const

overridevirtual

Print the content of the class.

Parameters

title

title to be printed

Reimplemented from TOPGeoBase.

Definition at line 30 of file TOPGeoSideRails.cc.

  {
    TOPGeoBase::print(title);
 
    cout << " thickness = " << getThickness() << " " << s_unitName;
    cout << ", reduced thickness = " << getReducedThickness() << " " << s_unitName;
    cout << ", height = " << getHeight() << " " << s_unitName;
    cout << ", material = " << getMaterial() << endl;
 
  }

print() [18/22]

void print ( const std::string & title = "Nominal quantum efficiency" ) const

overridevirtual

Print the content of the class.

Parameters

title

title to be printed

Reimplemented from TOPGeoBase.

Definition at line 68 of file TOPNominalQE.cc.

  {
    TOPGeoBase::printUnderlined(title);
    cout << " first point: " << getLambdaFirst() << " nm";
    cout << ", step: " << getLambdaStep() << " nm";
    cout << ", CE: " << getCE() << endl;
    cout << " QE: [";
    for (const auto& qe : m_QE) cout << qe << ", ";
    cout << "]" << endl;
  }

print() [19/22]

void print ( const std::string & title = "Nominal time-to-digit conversion parameters" ) const

overridevirtual

Print the content of the class.

Parameters

title

title to be printed

Reimplemented from TOPGeoBase.

Definition at line 93 of file TOPNominalTDC.cc.

  {
    TOPGeoBase::printUnderlined(title);
    cout << " range: [" << getTimeMin() << ", " << getTimeMax() << "] ns" << endl;
    cout << " offset: " << getOffset() << " ns" << endl;
    cout << " number of bits: " << getNumBits() << endl;
    cout << " bin width: " << getBinWidth() / Unit::ps << " ps" << endl;
    cout << " number of ASIC windows: " << getNumWindows() << endl;
    cout << " number of bits per sample: " << getSubBits() << endl;
    cout << " synchonization time base: " << getSyncTimeBase() << " ns" << endl;
    cout << " pile-up time: " << getPileupTime() << " ns" << endl;
    cout << " double hit resolution: " << getDoubleHitResolution() << " ns" << endl;
    cout << " r.m.s. of time jitter: " << getTimeJitter() / Unit::ps << " ps" << endl;
    cout << " electronic efficiency: " << getEfficiency() << endl;
  }

print() [20/22]

void print ( const std::string & title = "Nominal TTS distribution" ) const

overridevirtual

Print the content of the class.

Parameters

title

title to be printed

Reimplemented from TOPGeoBase.

Definition at line 68 of file TOPNominalTTS.cc.

  {
    TOPGeoBase::printUnderlined(title);
    if (m_tts.empty()) {
      cout << " -> vector is empty!" << endl;
      return;
    }
    cout << " PMT type: " << m_type << endl;
 
    int i = 0;
    for (const auto& tts : m_tts) {
      i++;
      cout << " Gaussian " << i << ": ";
      cout << tts.fraction << ", ";
      cout << tts.position << " ns, ";
      cout << tts.sigma << " ns" << endl;
    }
    if (!m_normalized) cout << " warning: fractions not normalized" << endl;
  }

print() [21/22]

void print

(

)

const

Print the class content.

Definition at line 21 of file TOPPmtGainPar.cc.

  {
    cout << "PMT gain parameters:" << endl;
    cout << " Serial number: " << getSerialNumber() << endl;
    cout << " Nominal HV: " << getNominalHV() << " V" << endl;
    for (unsigned channel = 1; channel <= c_NumPmtPixels; channel++) {
      cout << setw(4) << channel;
      cout << setw(12) << getConstant(channel);
      cout << setw(12) << getSlope(channel);
      cout << setw(12) << getRatio(channel);
      cout << setw(12) << getGain0(channel, getNominalHV());
      cout << endl;
    }
    cout << endl;
  }

print() [22/22]

void print ( const std::string & title = "Wavelength filter bulk transmittance" ) const

overridevirtual

Print the content of the class.

Parameters

title

title to be printed

Reimplemented from TOPGeoBase.

Definition at line 45 of file TOPWavelengthFilter.cc.

  {
    TOPGeoBase::printUnderlined(title);
    cout << " first point: " << getLambdaFirst() << " nm";
    cout << ", step: " << getLambdaStep() << " nm" << endl;
    cout << " bulk transmittances: [";
    for (const auto& trE : m_transmittances) cout << setprecision(3) << trE << ", ";
    cout << "]" << endl;
    if (getName().empty()) {
      cout << "(Old geometry payload, filter transmittance is included in QE)" << endl;
    }
 
  }

printBackPlateContour()

void printBackPlateContour

(

)

const

Print prism enclosure back plate contour (polygon) - for debugging.

Definition at line 194 of file TOPGeoQBB.cc.

  {
    auto contour = getBackPlateContour();
    cout << "Prism enclosure back plate contour [" << s_unitName << "]:" << endl;
    for (auto& xy : contour) cout << xy.first << " " << xy.second << endl;
    cout << endl;
  }

printBackwardContour()

void printBackwardContour

(

)

const

Print backward contour (polygon) - for debugging.

Definition at line 130 of file TOPGeoQBB.cc.

  {
    auto contour = getBackwardContour();
    cout << "Backward contour [" << s_unitName << "]:" << endl;
    for (auto& xy : contour) cout << xy.first << " " << xy.second << endl;
    cout << endl;
  }

printForwardContour()

void printForwardContour

(

)

const

Print forward contour (polygon) - for debugging.

Definition at line 45 of file TOPGeoQBB.cc.

  {
    auto contour = getForwardContour();
    cout << "Forward contour [" << s_unitName << "]:" << endl;
    for (auto& xy : contour) cout << xy.first << " " << xy.second << endl;
    cout << endl;
  }

printFrontPlateContour()

void printFrontPlateContour

(

)

const

Print prism enclosure front plate contour (polygon) - for debugging.

Definition at line 221 of file TOPGeoQBB.cc.

  {
    auto contour = getFrontPlateContour();
    cout << "Prism enclosure front plate contour [" << s_unitName << "]:" << endl;
    for (auto& xy : contour) cout << xy.first << " " << xy.second << endl;
    cout << endl;
  }

printInnerPanelContour()

void printInnerPanelContour

(

)

const

Print inner honeycomb panel contour (polygon) - for debugging.

Definition at line 71 of file TOPGeoQBB.cc.

  {
    auto contour = getInnerPanelContour();
    cout << "Inner panel contour [" << s_unitName << "]:" << endl;
    for (auto& xy : contour) cout << xy.first << " " << xy.second << endl;
    cout << endl;
  }

printOuterPanelContour()

void printOuterPanelContour

(

)

const

Print outer honeycomb panel contour (polygon) - for debugging.

Definition at line 97 of file TOPGeoQBB.cc.

  {
    auto contour = getOuterPanelContour();
    cout << "Outer panel contour [" << s_unitName << "]:" << endl;
    for (auto& xy : contour) cout << xy.first << " " << xy.second << endl;
    cout << endl;
  }

printPrismEnclosureContour()

void printPrismEnclosureContour

(

)

const

Print prism enclosure contour (polygon) - for debugging.

Definition at line 165 of file TOPGeoQBB.cc.

  {
    auto contour = getPrismEnclosureContour();
    cout << "Prism enclosure contour [" << s_unitName << "]:" << endl;
    for (auto& xy : contour) cout << xy.first << " " << xy.second << endl;
    cout << endl;
  }

printSampleTimeCalibration()

void printSampleTimeCalibration

(

)

Print sample time calibration constants stored in database.

Definition at line 601 of file TOPDatabaseImporter.cc.

  {
 
    DBObjPtr<TOPCalTimebase> timeBase;
    if (!timeBase.isValid()) {
      B2ERROR("No time base calibration available");
      return;
    }
 
    for (const auto& sampleTimes : timeBase->getSampleTimes()) {
      cout << sampleTimes.getScrodID() << " " << sampleTimes.getChannel() << endl;
      for (const auto& time : sampleTimes.getTimeAxis()) {
        cout << time << " ";
      }
      cout << endl;
    }
 
  }

printSampleTimeCalibrationInfo()

void printSampleTimeCalibrationInfo

(

)

Prints sample time calibration info about constants stored in database.

Definition at line 561 of file TOPDatabaseImporter.cc.

  {
    DBObjPtr<TOPCalTimebase> timeBase;
    if (!timeBase.isValid()) {
      B2ERROR("No time base calibration available");
      return;
    }
 
    const auto* geo = TOPGeometryPar::Instance()->getGeometry();
    int numModules = geo->getNumModules();
    auto& feMapper = TOPGeometryPar::Instance()->getFrontEndMapper();
 
    cout << "Time base calibration: number of calibrated channels in database" << endl << endl;
    for (int moduleID = 1; moduleID <= numModules; moduleID++) {
      int ncal[4] = {0, 0, 0, 0};
      int scrodID[4] = {0, 0, 0, 0};
      for (int bs = 0; bs < 4; bs++) {
        auto* femap = feMapper.getMap(moduleID, bs);
        if (!femap) {
          B2ERROR("No FrontEnd map available for boardstack " << bs << " of module " << moduleID);
          continue;
        }
        scrodID[bs] = femap->getScrodID();
        for (int channel = 0; channel < 128; channel++) {
          if (timeBase->isAvailable(scrodID[bs], channel)) ncal[bs]++;
        }
      }
 
      cout << "Slot " << moduleID << endl;
      for (int bs = 0; bs < 4; bs++) {
        cout << "  scrodID " << scrodID[bs] << ": " << ncal[bs] << "/128";
        if (ncal[bs] < 128) cout << " (" << 128 - ncal[bs] << " missing)";
        cout << endl;
      }
    }
 
    cout << endl;
  }

printSurface()

void printSurface ( const GeoOpticalSurface & surface ) const

virtual

Print the content of optical surface.

Parameters

surface

optical surface parameters

Definition at line 44 of file TOPGeoBase.cc.

  {
    cout << " Optical surface: ";
    if (surface.getName().empty() and !surface.hasProperties()) {
      cout << "not defined" << endl;
      return;
    }
    cout << surface.getName();
    cout << ", model: " << surface.getModel();
    cout << ", finish: " << surface.getFinish();
    cout << ", type: " << surface.getType();
    cout << ", value: " << surface.getValue();
    cout << endl;
    if (surface.hasProperties()) {
      for (const auto& property : surface.getProperties()) {
        cout << "  - property: ";
        cout << property.getName() << " [";
        for (const auto& value : property.getValues()) cout << value << ", ";
        cout << "], @[";
        for (const auto& value : property.getEnergies()) cout << value / Unit::eV << ", ";
        cout << "] eV" << endl;
      }
    } else {
      cout << "  - properties: None" << endl;
    }
 
  }

printUnderlined()

void printUnderlined ( const std::string & title ) const

protected

Print the content of the class.

Parameters

title

title to be printed

Definition at line 35 of file TOPGeoBase.cc.

  {
    cout << title << ":" << endl;
    for (size_t i = 0; i <= title.length(); i++) cout << "-";
    cout << endl;
    cout << " name: " << m_name << endl;
  }

reflect()

void reflect ( const ROOT::Math::XYVector * points, const ROOT::Math::XYVector * normals, const ROOT::Math::XYVector & orig, const ROOT::Math::XYVector & surf, const ROOT::Math::XYVector & norm, const ROOT::Math::XYVector & slanted, int k, std::vector< UnfoldedWindow > & result ) const

private

Do unfolding.

Parameters

points	points defining position of upper and slanted surfaces
normals	normals of upper and slanted surfaces
orig	origin of exit window
surf	exit window surface direction
norm	exit window surface normal
slanted	slanted surface normal
k	index of the surface to start unfolding (0 or 1)
result	the result of unfolding

Definition at line 152 of file TOPGeoPrism.cc.

  {
    XYVector rp[2] = {points[0], points[1]};
    XYVector n[2] = {normals[0], normals[1]};
    auto r = orig;
    auto s = surf;
    auto q = norm;
    auto sl = slanted;
 
    while (rp[k].Y() < 0) {
      r -= 2 * ((r - rp[k]).Dot(n[k])) * n[k]; // reflect window origin
      s -= 2 * s.Dot(n[k]) * n[k];           // reflect window surface direction
      q -= 2 * q.Dot(n[k]) * n[k];           // reflect window normal
      sl -= 2 * sl.Dot(n[k]) * n[k];         // reflect slanted normal
      result.push_back(UnfoldedWindow(r, s, q, sl));
      if (result.size() > 100) {
        B2ERROR("TOPGeoPrism::reflect: too many reflections -> must be a bug");
        return;
      }
      int i = (k + 1) % 2; // index of opposite surface
      rp[i] -= 2 * (rp[i] - rp[k]).Dot(n[k]) * n[k]; // reflect point on the opposite surface
      n[i] -= 2 * n[i].Dot(n[k]) * n[k];             // reflect normal of opposite surface
      k = i; // toggle the reflection surface
    }
  }

set() [1/2]

void set	(	int	slot,
		const TH1F *	photonYields,
		const TH1F *	backgroundYields,
		const TH1F *	alphaRatio,
		const TH1F *	activePixels,
		const TH2F *	pulseHeights,
		const TH1F *	muonZ )

Sets the data of a given slot.

The data are given as 1D histograms of 512 bins that correspond to pixelID's. Note: the setter must be called in the ascending order of slot ID's, e.g slot = 1, 2, ... , 16

Parameters

slot	slot ID
photonYields	pixel photon yields
backgroundYields	pixel background yields
alphaRatio	equalized alpha ratio
activePixels	active pixels
pulseHeights	pixel pulse-heights
muonZ	z-distribution of tracks

Definition at line 20 of file TOPCalPhotonYields.cc.

  {
    string slotName = (slot < 10) ? "_0" + to_string(slot) : "_" + to_string(slot);
 
    m_photonYields.push_back(TH2F(("photonYields" + slotName).c_str(),
                                  ("Photon yields for slot " + to_string(slot) + "; pixel column; pixel row").c_str(),
                                  c_numCols, 0.5, c_numCols + 0.5, c_numRows, 0.5, c_numRows + 0.5));
    copyContent(photonYields, m_photonYields.back());
 
    m_backgroundYields.push_back(TH2F(("backgroundYields" + slotName).c_str(),
                                      ("Background yields for slot " + to_string(slot) + "; pixel column; pixel row").c_str(),
                                      c_numCols, 0.5, c_numCols + 0.5, c_numRows, 0.5, c_numRows + 0.5));
    copyContent(backgroundYields, m_backgroundYields.back());
 
    m_alphaRatio.push_back(TH2F(("alphaRatio" + slotName).c_str(),
                                ("Equalized alpha ratio for slot " + to_string(slot) + "; pixel column; pixel row").c_str(),
                                c_numCols, 0.5, c_numCols + 0.5, c_numRows, 0.5, c_numRows + 0.5));
    copyContent(alphaRatio, m_alphaRatio.back());
 
    m_activePixels.push_back(TH2F(("activePixels" + slotName).c_str(),
                                  ("Active pixels for slot " + to_string(slot) + "; pixel column; pixel row").c_str(),
                                  c_numCols, 0.5, c_numCols + 0.5, c_numRows, 0.5, c_numRows + 0.5));
    copyContent(activePixels, m_activePixels.back());
    m_activePixels.back().Scale(1 / muonZ->GetEntries());
 
    m_pulseHeights.push_back(*pulseHeights);
    m_muonZ.push_back(*muonZ);
  }

set() [2/2]

void set	(	const std::vector< double > &	timeWalkParams,
		double	noise,
		double	quadratic )

Sets all the parameters and switches status to calibrated.

Parameters

timeWalkParams	polynomial coefficients of time-walk calibration curve [ns]
noise	noise term excess coefficient of electronic time resolution [ns]
quadratic	quadratic term coefficient of electronic time resolution [ns]

Definition at line 18 of file TOPCalTimeWalk.cc.

  {
    m_timeWalkParams = timeWalkParams;
    m_noise = noise;
    m_quadratic = quadratic;
    m_status = c_Calibrated;
  }

setActive()

void setActive	(	int	moduleID,
		unsigned	channel )

Sets a specific channel as active.

Parameters

moduleID	module ID (1-based)
channel	hardware channel number (0-based)

Definition at line 31 of file TOPCalChannelMask.cc.

  {
    int module = moduleID - 1;
    if (!check(module, channel)) {
      B2WARNING("Channel status 'active' not set");
      return;
    }
    m_status[module][channel] = c_Active;
  }

setAlpha()

void setAlpha	(	int	moduleID,
		double	alpha,
		double	errAlpha )

Sets the angle alpha on a single module.

Parameters

moduleID	module ID (1-based)
alpha	rotation angle around x
errAlpha	error on alpha

Definition at line 25 of file TOPCalModuleAlignment.cc.

  {
    unsigned module = moduleID - 1;
    if (module >= c_numModules) {
      B2ERROR("Invalid module number, constant not set (" << ClassName() << ")");
      return;
    }
    m_alpha[module] = alpha;
    m_errAlpha[module] = errAlpha;
    m_transforms.clear();
  }

setBeta()

void setBeta	(	int	moduleID,
		double	beta,
		double	errBeta )

Sets the angle beta on a single module.

Parameters

moduleID	module ID (1-based)
beta	rotation angle around y
errBeta	error on beta

Definition at line 37 of file TOPCalModuleAlignment.cc.

  {
    unsigned module = moduleID - 1;
    if (module >= c_numModules) {
      B2ERROR("Invalid module number, constant not set (" << ClassName() << ")");
      return;
    }
    m_beta[module] = beta;
    m_errBeta[module] = errBeta;
    m_transforms.clear();
  }

setBrokenGlue()

void setBrokenGlue	(	int	glueID,
		double	fraction,
		double	angle,
		const std::string &	material )

Sets glue to be broken (delaminated)

Parameters

glueID	glue ID: bar-mirror (1), bar-bar (2), prism-bar (3)
fraction	fraction of the delaminated surface
angle	angle of the delaminated surface
material	material name to simulate the delaminated glue

Definition at line 72 of file TOPGeoModule.cc.

  {
    switch (glueID) {
      case 1:
        m_mirror.setGlueDelamination(fraction, angle, material);
        break;
      case 2:
        m_bar1.setGlueDelamination(fraction, angle, material);
        break;
      case 3:
        m_bar2.setGlueDelamination(fraction, angle, material);
        break;
      default:
        B2ERROR("TOPGeoModule::setBrokenGlue: invalid glue ID."
                << LogVar("glue ID", glueID));
    }
  }

setCalibrated()

void setCalibrated

(

int

moduleID

)

Switches calibration status to calibrated.

Parameters

moduleID

module ID (1-based)

Definition at line 97 of file TOPCalModuleAlignment.cc.

  {
    unsigned module = moduleID - 1;
    if (module >= c_numModules) {
      B2ERROR("Invalid module number, status not set (" << ClassName() << ")");
      return;
    }
    m_status[module] = c_Calibrated;
    m_transforms.clear();
  }

setDead()

void setDead	(	int	moduleID,
		unsigned	channel )

Sets a specific channel as dead.

Parameters

moduleID	module ID (1-based)
channel	hardware channel number (0-based)

Definition at line 41 of file TOPCalChannelMask.cc.

  {
    int module = moduleID - 1;
    if (!check(module, channel)) {
      B2WARNING("Channel status 'dead' not set");
      return;
    }
    m_status[module][channel] = c_Dead;
  }

setEnvelopeQE()

void setEnvelopeQE ( ) const

private

Sets envelope quantum efficiency.

Definition at line 87 of file TOPPmtQE.cc.

  {
    if (!m_envelopeQE.empty()) return;
 
    size_t size = 0;
    for (unsigned pixel = 0; pixel < c_NumPmtPixels; pixel++) {
      const auto& QE = m_QE[pixel];
      size = std::max(size, QE.size());
    }
    m_envelopeQE.resize(size, 0);
    for (unsigned pixel = 0; pixel < c_NumPmtPixels; pixel++) {
      const auto& QE = m_QE[pixel];
      for (size_t i = 0; i < QE.size(); i++) {
        m_envelopeQE[i] = std::max(m_envelopeQE[i], QE[i]);
      }
    }
  }

setGains()

bool setGains	(	const std::vector< float > &	gains,
		float	error = 0 )

Set gains.

Parameters

gains	vector of relative gains (size must be c_WindowSize)
error	error on gains

Returns

true, on success

Definition at line 20 of file TOPASICGains.cc.

  {
 
    if (gains.size() != c_WindowSize) {
      B2ERROR("TOPASICGains::setGains:  vector with wrong number of elements");
      return false;
    }
 
    for (int i = 0; i < c_WindowSize; i++) {
      float gain = gains[i] * m_unit;
      if (gain > 0 and (gain + 0.5) < 0x10000) {
        m_gains[i] = int(gain + 0.5);
      } else {
        return false;
      }
    }
    error *= m_unit;
    m_gainError = int(error + 0.5);
    return true;
 
  }

setGamma()

void setGamma	(	int	moduleID,
		double	gamma,
		double	errGamma )

Sets the angle gamma on a single module.

Parameters

moduleID	module ID (1-based)
gamma	rotation angle around z
errGamma	error on gamma

Definition at line 49 of file TOPCalModuleAlignment.cc.

  {
    unsigned module = moduleID - 1;
    if (module >= c_numModules) {
      B2ERROR("Invalid module number, constant not set (" << ClassName() << ")");
      return;
    }
    m_gamma[module] = gamma;
    m_errGamma[module] = errGamma;
    m_transforms.clear();
  }

setNoisy()

void setNoisy	(	int	moduleID,
		unsigned	channel )

Sets a specific channel as noisy.

Parameters

moduleID	module ID (1-based)
channel	hardware channel number (0-based)

Definition at line 51 of file TOPCalChannelMask.cc.

  {
    int module = moduleID - 1;
    if (!check(module, channel)) {
      B2WARNING("Channel status 'noisy' not set");
      return;
    }
    m_status[module][channel] = c_Noisy;
  }

setPedestals()

int setPedestals	(	const TProfile *	profile,
		double	average = 0 )

Set pedestals from profile histogram with c_WindowSize bins.

Parameters

profile	profile histogram
average	optional value to add

Returns

number of pedestals that can not be packed into 16-bit word

Definition at line 23 of file TOPASICPedestals.cc.

  {
    int bad = c_WindowSize;
 
    if (!profile) return bad;
 
    if (profile->GetNbinsX() != c_WindowSize) {
      B2ERROR("TOPASICPedestals::setPedestals:  TProfile with wrong number of bins");
      return bad;
    }
 
    if (profile->GetEntries() == 0) {
      B2WARNING("TOPASICPedestals::setPedestals:  TProfile with no entries");
      return bad;
    }
 
    std::vector<unsigned> values;
    std::vector<unsigned> errors;
    for (int i = 0; i < c_WindowSize; i++) {
      values.push_back(int(profile->GetBinContent(i + 1) + average + 0.5));
      errors.push_back(int(profile->GetBinError(i + 1) + 0.5));
    }
 
    unsigned maxDif = 1 << c_Bits; maxDif--;
    unsigned maxErr = 1 << (sizeof(unsigned short) * 8 - c_Bits); maxErr--;
    m_offset = getOptimizedOffset(values, errors, maxDif, maxErr);
 
    for (int i = 0; i < c_WindowSize; i++) {
      unsigned dif = values[i] - m_offset;
      if (dif > maxDif) continue;
      unsigned err = errors[i];
      if (err > maxErr) continue;
      if (err == 0) err = 1;
      m_pedestals[i] = dif + (err << c_Bits);
      bad--;
    }
 
    return bad;
 
  }

setPeelOffRegions()

void setPeelOffRegions	(	double	thickness,
		const std::string &	material )

Sets parameters of the peel-off cookie volumes.

Parameters

thickness	volume thickness
material	material name

Definition at line 92 of file TOPGeoModule.cc.

  {
    double size = 2 * m_pmtArray.getDx();
    double offset = (m_pmtArray.getX(1) + m_pmtArray.getX(2)) / 2 +
                    m_arrayDisplacement.getX();
    m_prism.setPeelOffRegions(size, offset, thickness, material);
  }

setStatus()

void setStatus	(	int	moduleID,
		unsigned	channel,
		EStatus	status )

Sets the status for a single channel.

Parameters

moduleID	module ID (1-based)
channel	hardware channel number (0-based)
status	channel status (0: active 1: dead 2: noisy )

Definition at line 21 of file TOPCalChannelMask.cc.

  {
    int module = moduleID - 1;
    if (!check(module, channel)) {
      B2WARNING("Channel status not set");
      return;
    }
    m_status[module][channel] = status;
  }

setT0() [1/3]

void setT0	(	int	moduleID,
		unsigned	asic,
		double	T0 )

Sets calibration for a single ASIC and switches status to calibrated.

Parameters

moduleID	module ID (1-based)
asic	ASIC number within a module (0 - 63)
T0	shift in time [ns]

Definition at line 20 of file TOPCalAsicShift.cc.

  {
    unsigned module = moduleID - 1;
    if (module >= c_numModules) {
      B2ERROR("Invalid module number, constant not set (" << ClassName() << ")");
      return;
    }
    if (asic >= c_numAsics) {
      B2ERROR("Invalid asic number, constant not set (" << ClassName() << ")");
      return;
    }
    m_T0[module][asic] = T0;
    m_status[module][asic] = c_Calibrated;
  }

setT0() [2/3]

void setT0	(	int	moduleID,
		unsigned	channel,
		double	T0,
		double	errT0 )

Sets calibration for a single channel and switches status to calibrated.

Parameters

moduleID	module ID (1-based)
channel	hardware channel number (0-based)
T0	channel T0
errT0	error on T0

Definition at line 20 of file TOPCalChannelT0.cc.

  {
    unsigned module = moduleID - 1;
    if (module >= c_numModules) {
      B2ERROR("Invalid module number, constant not set (" << ClassName() << ")");
      return;
    }
    if (channel >= c_numChannels) {
      B2ERROR("Invalid channel number, constant not set (" << ClassName() << ")");
      return;
    }
    m_T0[module][channel] = T0;
    m_errT0[module][channel] = errT0;
    m_status[module][channel] = c_Calibrated;
  }

setT0() [3/3]

void setT0	(	int	moduleID,
		double	T0,
		double	errT0 )

Sets calibration for a single module and switches status to calibrated.

Parameters

moduleID	module ID (1-based)
T0	module T0
errT0	error on T0

Definition at line 20 of file TOPCalModuleT0.cc.

  {
    unsigned module = moduleID - 1;
    if (module >= c_numModules) {
      B2ERROR("Invalid module number, constant not set (" << ClassName() << ")");
      return;
    }
    m_T0[module] = T0;
    m_errT0[module] = errT0;
    m_status[module] = c_Calibrated;
  }

setTimeAxis() [1/2]

void setTimeAxis	(	const std::vector< double > &	sampleTimes,
		double	syncTimeBase )

Sets time axis from calibration data and switches status to calibrated.

Parameters

sampleTimes	vector of 256 elements of sample times
syncTimeBase	sinchronization time base (width of 2 ASIC windows)

Definition at line 33 of file TOPSampleTimes.cc.

  {
    if (sampleTimes.size() < c_TimeAxisSize) {
      B2FATAL("TOPSampleTimes::setTimeAxis: vector too short");
      return;
    }
 
    for (unsigned i = 0; i < c_TimeAxisSize; i++) m_timeAxis[i] = sampleTimes[i];
    double DTime = 2 * syncTimeBase;
    m_timeAxis[c_TimeAxisSize] = DTime;
    m_calibrated = c_Calibrated;
  }

setTimeAxis() [2/2]

void setTimeAxis

(

double

syncTimeBase

)

Sets equidistant time axis (uncalibrated).

Parameters

syncTimeBase

sinchronization time base (width of 2 ASIC windows)

Definition at line 23 of file TOPSampleTimes.cc.

  {
    double DTime = 2 * syncTimeBase;
    double timeBin = DTime / static_cast<double>(c_TimeAxisSize);
    for (unsigned i = 0; i < c_TimeAxisSize; i++) m_timeAxis[i] = timeBin * i;
    m_timeAxis[c_TimeAxisSize] = DTime;
    m_calibrated = c_Default;
  }

setTransformation()

void setTransformation ( ) const

private

Sets transformation cache.

Definition at line 64 of file TOPGeoModule.cc.

  {
    RotationZ Rz(m_phi - M_PI / 2);
    m_transformNominal = new Transform3D(Rz, Rz * XYZVector(0, m_radius, getZc() * s_unit));
    m_transform = new Transform3D(*m_transformNominal * m_moduleDisplacement.getTransformation());
  }

setTransformations()

void setTransformations ( ) const

private

Sets the transformation cache.

Definition at line 282 of file TOPCalModuleAlignment.cc.

  {
    for (int i = 0; i < c_numModules; i++) {
      Transform3D T;
      if (m_status[i] == c_Calibrated) {
        RotationX Rx(m_alpha[i]);
        RotationY Ry(m_beta[i]);
        RotationZ Rz(m_gamma[i]);
        Translation3D t(m_x[i], m_y[i], m_z[i]);
        T = Transform3D(Rz * Ry * Rx, t);
      }
      m_transforms.push_back(T);
    }
    m_transforms.push_back(Transform3D()); // for invalid module ID
  }

setUnusable() [1/4]

void setUnusable	(	int	moduleID,
		unsigned	asic )

Switches calibration status to unusable to flag badly calibrated constant.

Parameters

moduleID	module ID (1-based)
asic	ASIC number within a module (0 - 63)

Definition at line 36 of file TOPCalAsicShift.cc.

  {
    unsigned module = moduleID - 1;
    if (module >= c_numModules) {
      B2ERROR("Invalid module number, status not set (" << ClassName() << ")");
      return;
    }
    if (asic >= c_numAsics) {
      B2ERROR("Invalid asic number, status not set (" << ClassName() << ")");
      return;
    }
    m_status[module][asic] = c_Unusable;
  }

setUnusable() [2/4]

void setUnusable	(	int	moduleID,
		unsigned	channel )

Switches calibration status to unusable to flag badly calibrated constant.

Parameters

moduleID	module ID (1-based)
channel	hardware channel number (0-based)

Definition at line 37 of file TOPCalChannelT0.cc.

  {
    unsigned module = moduleID - 1;
    if (module >= c_numModules) {
      B2ERROR("Invalid module number, status not set (" << ClassName() << ")");
      return;
    }
    if (channel >= c_numChannels) {
      B2ERROR("Invalid channel number, status not set (" << ClassName() << ")");
      return;
    }
    m_status[module][channel] = c_Unusable;
  }

setUnusable() [3/4]

void setUnusable

(

int

moduleID

)

Switches calibration status to unusable to flag badly calibrated constant.

Parameters

moduleID

module ID (1-based)

Definition at line 108 of file TOPCalModuleAlignment.cc.

  {
    unsigned module = moduleID - 1;
    if (module >= c_numModules) {
      B2ERROR("Invalid module number, status not set (" << ClassName() << ")");
      return;
    }
    m_status[module] = c_Unusable;
    m_transforms.clear();
  }

setUnusable() [4/4]

void setUnusable

(

int

moduleID

)

Switches calibration status to unusable to flag badly calibrated constant.

Parameters

moduleID

module ID (1-based)

Definition at line 33 of file TOPCalModuleT0.cc.

  {
    unsigned module = moduleID - 1;
    if (module >= c_numModules) {
      B2ERROR("Invalid module number, status not set (" << ClassName() << ")");
      return;
    }
    m_status[module] = c_Unusable;
  }

setWindowShifts()

void setWindowShifts

(

std::vector< int >

shifts

)

Sets window shifts.

Parameters

shifts

vector of size = 6 with window shifts according to revo9count

Definition at line 20 of file TOPFrontEndSetting.cc.

  {
    if (shifts.size() != 6) {
      B2ERROR("TOPFrontEndSetting::setWindowShifts: vector size must be 6");
      return;
    }
    for (int i = 0; i < 6; i++) m_windowShifts[i] = shifts[i];
  }

setX()

void setX	(	int	moduleID,
		double	x,
		double	errX )

Sets the displacement x on a single module.

Parameters

moduleID	module ID (1-based)
x	displacement along the x axis
errX	error on the displacement

Definition at line 61 of file TOPCalModuleAlignment.cc.

  {
    unsigned module = moduleID - 1;
    if (module >= c_numModules) {
      B2ERROR("Invalid module number, constant not set (" << ClassName() << ")");
      return;
    }
    m_x[module] = x;
    m_errX[module] = errX;
    m_transforms.clear();
  }

setY()

void setY	(	int	moduleID,
		double	y,
		double	errY )

Sets the displacement y on a single module.

Parameters

moduleID	module ID (1-based)
y	displacement along the y axis
errY	error on the displacement

Definition at line 73 of file TOPCalModuleAlignment.cc.

  {
    unsigned module = moduleID - 1;
    if (module >= c_numModules) {
      B2ERROR("Invalid module number, constant not set (" << ClassName() << ")");
      return;
    }
    m_y[module] = y;
    m_errY[module] = errY;
    m_transforms.clear();
  }

setZ()

void setZ	(	int	moduleID,
		double	z,
		double	errZ )

Sets the displacement z on a single module.

Parameters

moduleID	module ID (1-based)
z	displacement along the z axis
errZ	error on the displacement

Definition at line 85 of file TOPCalModuleAlignment.cc.

  {
    unsigned module = moduleID - 1;
    if (module >= c_numModules) {
      B2ERROR("Invalid module number, constant not set (" << ClassName() << ")");
      return;
    }
    m_z[module] = z;
    m_errZ[module] = errZ;
    m_transforms.clear();
  }

suppressAverage() [1/2]

void suppressAverage

(

)

Subtracts arithmetic average of a module from constants whose status is not c_Default.

Arithmetic average is calculated from those whose status is c_Calibrated.

Definition at line 52 of file TOPCalChannelT0.cc.

  {
    for (int m = 0; m < c_numModules; m++) {
      float s = 0;
      int n = 0;
      for (int i = 0; i < c_numChannels; i++) {
        if (m_status[m][i] == c_Calibrated) {
          s += m_T0[m][i];
          n++;
        }
      }
      if (n == 0) continue;
      s /= n;
      for (int i = 0; i < c_numChannels; i++) {
        if (m_status[m][i] != c_Default) m_T0[m][i] -= s;
      }
      B2INFO("Slot " << m + 1 << ": average of " << s  << " ns subtracted.");
    }
  }

suppressAverage() [2/2]

void suppressAverage

(

)

Subtracts arithmetic average from constants whose status is not c_Default.

Arithmetic average is calculated from those whose status is c_Calibrated.

Definition at line 44 of file TOPCalModuleT0.cc.

  {
    float s = 0;
    int n = 0;
    for (int i = 0; i < c_numModules; i++) {
      if (m_status[i] == c_Calibrated) {
        s += m_T0[i];
        n++;
      }
    }
    if (n == 0) return;
    s /= n;
    for (int i = 0; i < c_numModules; i++) {
      if (m_status[i] != c_Default) m_T0[i] -= s;
    }
    B2INFO("Average of " << s  << " ns subtracted.");
  }

TOPGeoModule()

TOPGeoModule

(

const TOPGeoModule &

module

)

Copy constructor.

Definition at line 26 of file TOPGeoModule.cc.

                                                      : TOPGeoBase(module.getName())
  {
    *this = module;
    m_transform = 0;
    m_transformNominal = 0;
  }

TOPNominalTDC()

TOPNominalTDC	(	int	numWindows,
		int	subBits,
		double	syncTimeBase,
		int	numofBunches,
		double	offset,
		double	pileupTime,
		double	doubleHitResolution,
		double	timeJitter,
		double	efficiency,
		const std::string &	name = "TOPNominalTDC" )

Useful constructor (new xml file version)

Parameters

numWindows	number of ASIC windows per waveform
subBits	number of bits per sample
syncTimeBase	synchonization time base (time width of c_syncWindows)
numofBunches	number of bunches per syncTimeBase
offset	time offset
pileupTime	pile-up time
doubleHitResolution	double hit resolution time
timeJitter	r.m.s. of time jitter
efficiency	electronic efficiency (fraction of hits above threshold)
name	object name

Definition at line 22 of file TOPNominalTDC.cc.

                                                     :
    TOPGeoBase(name),
    m_numWindows(numWindows), m_subBits(subBits), m_syncTimeBase(syncTimeBase),
    m_numofBunches(numofBunches),
    m_offset(offset), m_pileupTime(pileupTime), m_doubleHitResolution(doubleHitResolution),
    m_timeJitter(timeJitter), m_efficiency(efficiency)
  {
    if (numWindows <= 0)
      B2FATAL("TOPNominalTDC: numWindows must be > 0");
    if (subBits < 0)
      B2FATAL("TOPNominalTDC: subBits must be >= 0");
 
    int numSamples = numWindows * c_WindowSize;
    int numBits = subBits - 1;
    int k = numSamples;
    do {
      numBits++;
      k /= 2;
    } while (k > 0);
    if (numSamples > (1 << (numBits - subBits))) numBits++;
    m_numBits = numBits;
 
    int syncSamples = c_syncWindows * c_WindowSize;
    m_sampleWidth = syncTimeBase / syncSamples;
    m_binWidth = m_sampleWidth / (1 << subBits);
  }

TOPSignalShape()

TOPSignalShape	(	std::vector< double >	shape,
		double	timeBin,
		double	tau,
		double	pole1,
		double	pole2 )

Full constructor.

Parameters

shape	waveform values (unnormalized) of positive pulse
timeBin	time difference between two samples [ns]
tau	time constant of the exponential tail used for extrapolation [ns]
pole1	noise bandwidth: first pole [GHz]
pole2	noise bandwidth: second pole [GHz]

Definition at line 21 of file TOPSignalShape.cc.

                                                            :
    m_shape(shape), m_tau(tau), m_pole1(pole1), m_pole2(pole2)
  {
    m_tmax = m_tmin + (shape.size() - 1) * timeBin;
 
    // find peaking value (positive pulse maximum!)
 
    int samplePeak = 0;
    for (unsigned i = 0; i < shape.size(); i++) {
      if (shape[i] > shape[samplePeak]) samplePeak = i;
    }
 
    TSpline5 spline("spline", m_tmin, m_tmax, shape.data(), shape.size());
    double t1 = m_tmin + (samplePeak - 1) * timeBin;
    double t2 = m_tmin + (samplePeak + 1) * timeBin;
    for (int i = 0; i < 20; i++) {
      double t = (t1 + t2) / 2;
      if (spline.Derivative(t) > 0) {
        t1 = t;
      } else {
        t2 = t;
      }
    }
    double vPeak = spline.Eval((t1 + t2) / 2);
 
    // normalize waveform
 
    for (auto& value : m_shape) value /= vPeak;
 
    // find 50% CF crossing time (positive pulse!)
 
    auto sampleRise = samplePeak;
    while (sampleRise >= 0 and m_shape[sampleRise] > 0.5) sampleRise--;
    t1 = m_tmin + sampleRise * timeBin;
    t2 = m_tmin + (sampleRise + 1) * timeBin;
    for (int i = 0; i < 20; i++) {
      double t = (t1 + t2) / 2;
      if (spline.Eval(t) < vPeak / 2) {
        t1 = t;
      } else {
        t2 = t;
      }
    }
    double crossingTime = (t1 + t2) / 2;
 
    // set 50% CF crossing to happen at t = 0
 
    m_tmin -= crossingTime;
    m_tmax -= crossingTime;
 
  }

unfold()

void unfold ( ) const

private

Unfold prism exit window.

Definition at line 129 of file TOPGeoPrism.cc.

  {
    double z = -(m_length - m_flatLength);
    double yUp = m_thickness / 2;
    double yDown = yUp - m_exitThickness;
    XYVector points[2] = {XYVector(yUp, z), XYVector(yDown, z)}; // points on upper and slanted surfaces
 
    double alpha = getAngle();
    XYVector normals[2] = {XYVector(1, 0), XYVector(-cos(alpha), sin(alpha))}; // normals of upper and slanted surfaces
 
    XYVector orig(0, z); // window origin
    XYVector surf(1, 0); // window surface direction (= upper surface normal)
    XYVector norm(0, -1); // window normal
    auto slanted = normals[1]; // slanted surface normal
 
    reflect(points, normals, orig, surf, norm, slanted, 1, m_unfoldedWindows); // unfolding down
    std::reverse(m_unfoldedWindows.begin(), m_unfoldedWindows.end());
    m_unfoldedWindows.push_back(UnfoldedWindow(orig, surf, norm, slanted)); // true window
    m_k0 = m_unfoldedWindows.size() - 1;
    reflect(points, normals, orig, surf, norm, slanted, 0, m_unfoldedWindows); // unfolding up
  }

UnfoldedWindow()

UnfoldedWindow	(	const ROOT::Math::XYVector &	orig,
		const ROOT::Math::XYVector &	dir,
		const ROOT::Math::XYVector &	norm,
		const ROOT::Math::XYVector &	slanted )

constructor

Parameters

orig	window origin
dir	window surface direction (= upper surface normal)
norm	window surface normal (pointing out of prism)
slanted	slanted surface normal

Definition at line 25 of file TOPGeoPrism.cc.

                                                                                            :
    y0(orig.X()), z0(orig.Y()), sy(dir.X()), sz(dir.Y()), ny(norm.X()), nz(norm.Y())
  {
    nsy[0] = sy;
    nsz[0] = sz;
    nsy[1] = slanted.X();
    nsz[1] = slanted.Y();
  }

~TOPGeoModule()

~TOPGeoModule

(

)

Destructor.

Definition at line 58 of file TOPGeoModule.cc.

  {
    if (m_transform) delete m_transform;
    if (m_transformNominal) delete m_transformNominal;
  }

Variable Documentation

s_unit

double s_unit = Unit::cm

staticprotected

conversion unit for length

Definition at line 86 of file TOPGeoBase.h.

s_unitName

std::string s_unitName

staticprotected

conversion unit name

Definition at line 87 of file TOPGeoBase.h.