top modules

Collaboration diagram for top modules:

Classes
class	TOPAlignmentCollectorModule
	Collector for geometrical alignment of a TOP module with dimuons or Bhabhas. More...
class	TOPAsicShiftsBS13dCollectorModule
	Collector for carrier shifts of BS13d. More...
class	TOPChannelMaskCollectorModule
	Collector for preparing masks of hot and dead channels. More...
class	TOPCommonT0BFCollectorModule
	Collector for common T0 calibration with a fit of bunch finder residuals (method BF) More...
class	TOPCommonT0LLCollectorModule
	Collector for common T0 calibration with neg. More...
class	TOPModuleT0DeltaTCollectorModule
	Collector for module T0 calibration with chi2 minimization of time differences between slots (method DeltaT). More...
class	TOPModuleT0LLCollectorModule
	Collector for module T0 calibration with neg. More...
class	TOPOffsetCollectorModule
	Collector for eventT0 and fill pattern offset calibrations. More...
class	TOPPhotonYieldsCollectorModule
	Collector for photon pixel yields aimed for PMT ageing studies and for finding optically decoupled PMT's. More...
class	TOPPulseHeightCollectorModule
	Collector for channel pulse-height distributions. More...
class	TOPValidationCollectorModule
	Collector for automatic validation of calibration with dimuon events. More...
class	OpticalGunModule
	Source of optical photons for the simulation of the TOP laser system. More...
class	TOPAlignerModule
	Alignment of TOP. More...
class	TOPBackgroundModule
	TOP backgound module. More...
class	TOPBunchFinderModule
	Bunch finder: searches for the bunch crossing where the interaction happened using track-based TOP likelihood. More...
class	TOPChannelMaskerModule
	Masks dead PMs from the reconstruction. More...
class	TOPChannelT0CalibratorModule
	A module for alternative channel T0 calibration with collision data Note: after this kind of calibration one cannot do the geometrical alignment This module can also be used to check the calibration. More...
class	TOPChannelT0MCModule
	TOP Channel T0 MC Extraction module (under development) More...
class	TOPCommonT0CalibratorModule
	A module for common T0 calibration with collision data (dimuons or bhabhas) More...
class	TOPCosmicT0FinderModule
	Event T0 finder for global cosmic runs. More...
class	TOPDigitizerModule
	TOP digitizer. More...
class	TOPTriggerDigitizerModule
	Digitizer that provides time stamps for TOP trigger. More...
class	TOPCalPulseGeneratorModule
	Generator of calibration pulses Output to TOPSimCalPulses. More...
class	TOPDoublePulseGeneratorModule
	Generator of double calibration pulses Output to TOPDigits. More...
class	TOPDQMModule
	TOP DQM histogrammer. More...
class	TOPGainEfficiencyCalculatorModule
	Module for channel-by-channel gain/efficiency analysis. More...
class	TOPLaserHitSelectorModule
	Module for pixel-by-pixel gain/efficiency analysis. More...
class	TOPGeometryParInitializerModule
	Class for initializing TOPGeometryPar. More...
class	TOPInterimFENtupleModule
	Module to produce ntuple from TOPDigits and TOPRawDigits. More...
class	TOPLaserCalibratorModule
	T0 Laser calibration module (under development) More...
class	TOPLaserCalibratorCollectorModule
	Collector module for the TOP channelT0 calibration and, more in general, for the time resolution studies using the laser and pulser data. More...
class	TOPLLScannerModule
	A module to perform the TOP PID likelihood scan and find the actual minimum as function of the mass. More...
class	TOPMCTrackMakerModule
	Constructs Tracks and ExtHits from MCParticles and TOPBarHits Utility needed for testing and debugging of TOP reconstruction. More...
class	TOPModuleT0CalibratorModule
	A module for module T0 calibration with collision data (dimuons or bhabhas) Useful when the geometrical alignment need not to be repeated. More...
class	TOPNtupleModule
	Module to write out a ntuple summarizing TOP reconstruction output. More...
class	TOPPackerModule
	Raw data packer. More...
class	TOPPDFCheckerModule
	Module for checking analytic PDF used in likelihood calculation. More...
class	TOPPDFDebuggerModule
	TOP reconstruction module. More...
class	TOPRawDigitConverterModule
	TOPRawDigits to TOPDigits converter. More...
class	TOPReconstructorModule
	TOP reconstruction module. More...
class	TOPRingPlotterModule
	A module to plot the x-t images from TOP, and in general study the distribution of the digits associated to the particles in a particleList. More...
class	TOPTBCComparatorModule
	Module for the comparison of different sets of time base correction (TBC) constants and to produce monitoring plots out of a given set. More...
struct	Hit
	Structure to hold some of the calpulse data. More...
struct	TwoTimes
	Structure to hold calpulse raw times expressed in samples since sample 0 of window 0. More...
class	TOPTimeBaseCalibratorModule
	Time base calibrator. More...
class	TOPTimeRecalibratorModule
	Utility module for re-calibrating time of TOPDigits pulseWidth and timeError are not changed although they may depend no calibration! More...
class	TOPUnpackerModule
	Raw data unpacker. More...
class	TOPWaveformFeatureExtractorModule
	Waveform feature extractor: module adds rawDigits that are found in waveforms by feature extraction but are not already present in RawDigits. More...
class	TOPWaveformQualityPlotterModule
	Plots and histograms of waveforms and feature extracted parameters. More...
class	TOPXTalkChargeShareSetterModule
	Crosstalk & chargeshare flag setter. More...

Enumerations
enum	{   c_NChannelPerAsic = 8 ,   c_NModule = 16 ,   c_NChannelPerPMT = 16 ,   c_NChannelPerPMTRow = 4 ,   c_NPMTPerRow = 16 ,   c_NPMTPerModule = 32 ,   c_NPixelPerRow = 64 ,   c_NPixelPerModule = 512 }
	enum for maximum number of array elements (# of hits per event)

Functions
	REG_MODULE (TOPAlignmentCollector)
	Register module.
	REG_MODULE (TOPAsicShiftsBS13dCollector)
	Register module.
	REG_MODULE (TOPChannelMaskCollector)
	Register module.
	REG_MODULE (TOPCommonT0BFCollector)
	Register module.
	REG_MODULE (TOPCommonT0LLCollector)
	Register module.
	REG_MODULE (TOPModuleT0DeltaTCollector)
	Register module.
	REG_MODULE (TOPModuleT0LLCollector)
	Register module.
	REG_MODULE (TOPOffsetCollector)
	Register module.
	REG_MODULE (TOPPhotonYieldsCollector)
	Register module.
	REG_MODULE (TOPPulseHeightCollector)
	Register module.
	REG_MODULE (TOPValidationCollector)
	Register module.
	REG_MODULE (OpticalGun)
	Register module.
	REG_MODULE (TOPAligner)
	Register module.
	REG_MODULE (TOPBackground)
	Register the Module.
	REG_MODULE (TOPBunchFinder)
	Register module.
	REG_MODULE (TOPChannelMasker)
	Register the Module.
	REG_MODULE (TOPChannelT0Calibrator)
	Register module.
	REG_MODULE (TOPChannelT0MC)
	Register module.
	REG_MODULE (TOPCommonT0Calibrator)
	Register module.
	REG_MODULE (TOPCosmicT0Finder)
	Register module.
	REG_MODULE (TOPDigitizer)
	Register the Module.
	REG_MODULE (TOPTriggerDigitizer)
	Register the Module.
	REG_MODULE (TOPCalPulseGenerator)
	Register module.
	REG_MODULE (TOPDoublePulseGenerator)
	Register module.
	REG_MODULE (TOPDQM)
	Register module.
	REG_MODULE (TOPGainEfficiencyCalculator)
	Register the Module.
	REG_MODULE (TOPLaserHitSelector)
	Register the Module.
	REG_MODULE (TOPGeometryParInitializer)
	Register module.
	REG_MODULE (TOPInterimFENtuple)
	Register module.
	REG_MODULE (TOPLaserCalibrator)
	Register module.
	REG_MODULE (TOPMCTrackMaker)
	Register module.
	REG_MODULE (TOPModuleT0Calibrator)
	Register module.
	REG_MODULE (TOPNtuple)
	Register module.
	REG_MODULE (TOPPacker)
	Register module.
	REG_MODULE (TOPPDFChecker)
	Register module.
	REG_MODULE (TOPPDFDebugger)
	Register the Module.
	REG_MODULE (TOPRawDigitConverter)
	Register module.
	REG_MODULE (TOPReconstructor)
	Register the Module.
	REG_MODULE (TOPTBCComparator)
	Register module.
	REG_MODULE (TOPTimeBaseCalibrator)
	Register module.
	REG_MODULE (TOPTimeRecalibrator)
	Register module.
	REG_MODULE (TOPUnpacker)
	Register module.
	REG_MODULE (TOPWaveformFeatureExtractor)
	Register module.
	REG_MODULE (TOPWaveformQualityPlotter)
	Register module.
	TOPAlignmentCollectorModule ()
	Constructor.
virtual void	prepare () final
	Replacement for initialize(). More...
virtual void	collect () final
	Replacement for event(). More...
	TOPAsicShiftsBS13dCollectorModule ()
	Constructor.
virtual void	prepare () final
	Replacement for initialize(). More...
virtual void	collect () final
	Replacement for event(). More...
	TOPChannelMaskCollectorModule ()
	Constructor.
virtual void	prepare () final
	Replacement for initialize(). More...
virtual void	collect () final
	Replacement for event(). More...
	TOPCommonT0BFCollectorModule ()
	Constructor.
virtual void	prepare () final
	Replacement for initialize(). More...
virtual void	collect () final
	Replacement for event(). More...
	TOPCommonT0LLCollectorModule ()
	Constructor.
virtual void	prepare () final
	Replacement for initialize(). More...
virtual void	collect () final
	Replacement for event(). More...
	TOPModuleT0DeltaTCollectorModule ()
	Constructor.
virtual void	prepare () final
	Replacement for initialize(). More...
virtual void	collect () final
	Replacement for event(). More...
	TOPModuleT0LLCollectorModule ()
	Constructor.
virtual void	prepare () final
	Replacement for initialize(). More...
virtual void	collect () final
	Replacement for event(). More...
	TOPOffsetCollectorModule ()
	Constructor.
virtual void	prepare () final
	Replacement for initialize(). More...
virtual void	collect () final
	Replacement for event(). More...
	TOPPhotonYieldsCollectorModule ()
	Constructor.
virtual void	prepare () final
	Replacement for initialize(). More...
virtual void	collect () final
	Replacement for event(). More...
	TOPPulseHeightCollectorModule ()
	Constructor.
virtual void	prepare () final
	Replacement for initialize(). More...
virtual void	collect () final
	Replacement for event(). More...
	TOPValidationCollectorModule ()
	Constructor.
virtual void	prepare () final
	Replacement for initialize(). More...
virtual void	startRun () final
	Replacement for beginRun(). More...
virtual void	collect () final
	Replacement for event(). More...
virtual void	closeRun () final
	Replacement for endRun(). More...
	OpticalGunModule ()
	Constructor.
virtual	~OpticalGunModule ()
	Destructor.
virtual void	initialize () override
	Initialize the Module. More...
virtual void	event () override
	Event processor.
bool	isInsideSlit (const ROOT::Math::XYZPoint &point, const ROOT::Math::XYZVector &direction) const
	Checks if photon passes the slit. More...
ROOT::Math::XYZVector	getDirectionGaussian () const
	Return photon direction according to a projected 2D gaussian distribution based on numerical aperture NA. More...
ROOT::Math::XYZVector	getDirectionUniform () const
	Return photon direction according to a projected uniform distribution with opening angle alpha. More...
ROOT::Math::XYZVector	getDirectionLambertian () const
	Return photon direction according to a lambertian distribution with opening angle alpha. More...
ROOT::Math::XYZVector	getDirectionCustom () const
	Return photon direction according to a custom angular distribution given by TFormula. More...
	TOPAlignerModule ()
	Constructor.
virtual void	initialize () override
	Initialize the Module. More...
virtual void	event () override
	Event processor.
virtual void	terminate () override
	Termination action. More...
	TOPBackgroundModule ()
	Constructor.
virtual	~TOPBackgroundModule ()
	Destructor.
virtual void	initialize () override
	Initialize the Module. More...
virtual void	beginRun () override
	Called when entering a new run. More...
virtual void	event () override
	Event processor. More...
void	myprint (TH1F *histo, const char *path, const char *xtit, const char *ytit, double tresh)
	Print histogram 1D, helper function.
virtual void	endRun () override
	End-of-run action. More...
virtual void	terminate () override
	Termination action. More...
void	printModuleParams () const
	Prints module parameters.
	TOPBunchFinderModule ()
	Constructor.
virtual void	initialize () override
	Initialize the Module. More...
virtual void	beginRun () override
	Called when entering a new run.
virtual void	event () override
	Event processor.
virtual void	terminate () override
	Termination action. More...
Const::ChargedStable	getMostProbable (const Track &track)
	Returns most probable charged stable particle according to dEdx and predefined prior probabilities. More...
int	setFinder (TOP::Chi2MinimumFinder1D &finder, const TOP::PDFConstructor &reco, double timeMin, double timeMax)
	Sets finder object with chi2 values. More...
TimeSeed	getTimeSeed ()
	Returns a time seed. More...
bool	isBucketFilled (int bunchNo)
	Does reconstructed bunch number correspond to filled bucket. More...
	TOPChannelMaskerModule ()
	Constructor: Sets the description of the module.
virtual void	initialize () override
	initialize method: registers datastore objects (the TOP hits)
virtual void	beginRun () override
	Called when entering a new run. More...
virtual void	event () override
	event method: removes channels from the reconstruction pdf, flags hits from noisy channels as junk
	TOPChannelT0CalibratorModule ()
	Constructor.
virtual void	initialize () override
	Initialize the Module. More...
virtual void	event () override
	Event processor.
virtual void	terminate () override
	Termination action. More...
	TOPChannelT0MCModule ()
	Constructor.
virtual	~TOPChannelT0MCModule ()
	Destructor.
virtual void	initialize () override
	Initialize the Module. More...
virtual void	beginRun () override
	Called when entering a new run.
virtual void	event () override
	Event processor.
virtual void	endRun () override
	End-of-run action.
virtual void	terminate () override
	Termination action. More...
	TOPCommonT0CalibratorModule ()
	Constructor.
virtual void	initialize () override
	Initialize the Module. More...
virtual void	event () override
	Event processor.
virtual void	terminate () override
	Termination action. More...
bool	isRunningOffsetSubtracted ()
	Checks if running offset is subtracted in TOPDigits. More...
	TOPCosmicT0FinderModule ()
	Constructor.
virtual void	initialize () override
	Initialize the Module. More...
virtual void	event () override
	Event processor.
virtual void	terminate () override
	Termination action. More...
	TOPDigitizerModule ()
	Constructor.
virtual void	initialize () override
	Initialize the Module. More...
virtual void	beginRun () override
	Called when entering a new run.
virtual void	event () override
	Event processor.
TimeOffset	getTimeOffset (double trgOffset, int moduleID, int pixelID)
	Returns a complete time offset by adding time mis-calibration to trgOffset. More...
double	generatePulseHeight (int moduleID, int pixelID) const
	Generates and returns pulse height. More...
	TOPTriggerDigitizerModule ()
	Constructor.
virtual	~TOPTriggerDigitizerModule ()
	Destructor.
virtual void	initialize () override
	Initialize the Module. More...
virtual void	beginRun () override
	Called when entering a new run.
virtual void	event () override
	Event processor.
virtual void	endRun () override
	End-of-run action.
virtual void	terminate () override
	Termination action. More...
	TOPCalPulseGeneratorModule ()
	Constructor.
virtual void	initialize () override
	Initialize the Module. More...
virtual void	event () override
	Event processor.
	TOPDoublePulseGeneratorModule ()
	Constructor.
virtual	~TOPDoublePulseGeneratorModule ()
	Destructor.
virtual void	initialize () override
	Initialize the Module. More...
virtual void	beginRun () override
	Called when entering a new run. More...
virtual void	event () override
	Event processor.
virtual void	endRun () override
	End-of-run action. More...
virtual void	terminate () override
	Termination action. More...
void	storeSampleTimes (std::string fileName)
	Optionally store sample times used by the generator as root histograms fileName root output file name.
void	saveAsHistogram (const std::vector< double > &vec, const std::string &name, const std::string &title, const std::string &xTitle="", const std::string &yTitle="") const
	Save vector to histogram and write it out. More...
	TOPDQMModule ()
	Constructor.
virtual	~TOPDQMModule ()
	Destructor.
virtual void	defineHisto () override
	Histogram definitions such as TH1(), TH2(), TNtuple(), TTree().... More...
virtual void	initialize () override
	Initialize the Module. More...
virtual void	beginRun () override
	Called when entering a new run. More...
virtual void	event () override
	Event processor.
int	getModuleID (const Track &track) const
	Returns slot ID of the module that is hit by the track. More...
	TOPGainEfficiencyCalculatorModule ()
	Constructor.
virtual	~TOPGainEfficiencyCalculatorModule ()
	Destructor.
virtual void	initialize () override
	Load time vs charge 2D histogram from a given input file (paramter "inputFile") and prepare hit timing and pulse charge distribution for each channel.
virtual void	defineHisto () override
	Define TTree branches to store fit results for each channel This TTree is saved in an output file given by "histoFileName" parameter of "HistoManager" module.
virtual void	beginRun () override
	The main processes, fitting charge distribution and calculating gain/efficiency, are done in this function.
virtual void	event () override
	This will be empty as the all the processes are done in beginRun() function thus input file can be a dummy file.
virtual void	endRun () override
	Draw plots to show fitting results for each channel and save them into a given PDF file (outputPDFFile).
virtual void	terminate () override
	Termination action. More...
void	LoadHistograms (const std::string &histotype)
	Load 2D histograms from a given input file (output of TOPLaserHitSelector) and create timing and charge distribution as projection histograms for the x- and y-axis, respectively. More...
void	FitHistograms (EHistogramType LoadHisto)
	Fit charge (or integrated charged) distribution to calculate gain and efficiency for each channel.
void	DummyFillBranch (EHistogramType LoadHisto)
	Fill Dummy for Branch. More...
void	DrawResult (const std::string &histotype, EHistogramType LoadHisto)
	Draw results of gain/efficiency calculation for each channel to a given output file.
static double	TOPGainFunc (double *var, double *par)
	Fit function of pulse charge (or charnge) distribution for channel(pixel)-by-channel gain extraction, given by "[0]*pow(x-[4],[1])*exp(-pow(x-[4],[2])/[3])" smeared by Gaussian with a constant sigma to consider baseline fluctuation.
static double	FindPeakForSmallerXThan (TH1 *histo, double xmax=0)
	Find peak and return its position for a limited range of x (x smaller than the given value (xmax))
	TOPLaserHitSelectorModule ()
	Constructor.
virtual	~TOPLaserHitSelectorModule ()
	Destructor.
virtual void	initialize () override
	Initialize the Module. More...
virtual void	defineHisto () override
	create timing-height 2D histograms for all 8192 pixels
virtual void	beginRun () override
	Called when entering a new run. More...
virtual void	event () override
	Event processor.
virtual void	endRun () override
	End-of-run action. More...
virtual void	terminate () override
	Termination action. More...
	TOPGeometryParInitializerModule ()
	Constructor.
virtual void	initialize () override
	Initialize the Module. More...
	TOPInterimFENtupleModule ()
	Constructor.
virtual	~TOPInterimFENtupleModule ()
	Destructor.
virtual void	initialize () override
	Initialize the Module. More...
virtual void	defineHisto () override
	Module funcions to define histograms.
virtual void	beginRun () override
	Called when entering a new run. More...
virtual void	event () override
	Event processor.
virtual void	endRun () override
	End-of-run action. More...
virtual void	terminate () override
	Termination action. More...
void	getReferenceTiming ()
	Find reference timing. More...
	TOPLaserCalibratorModule ()
	Constructor.
virtual	~TOPLaserCalibratorModule ()
	Destructor.
virtual void	initialize () override
	Initialize the Module. More...
virtual void	beginRun () override
	Called when entering a new run. More...
virtual void	event () override
	Event processor.
virtual void	endRun () override
	End-of-run action. More...
virtual void	terminate () override
	Termination action. More...
	TOPMCTrackMakerModule ()
	Constructor.
virtual	~TOPMCTrackMakerModule ()
	Destructor.
virtual void	initialize () override
	Initialize the Module. More...
virtual void	beginRun () override
	Called when entering a new run. More...
virtual void	event () override
	Event processor.
virtual void	endRun () override
	End-of-run action. More...
virtual void	terminate () override
	Termination action. More...
	TOPModuleT0CalibratorModule ()
	Constructor.
virtual void	initialize () override
	Initialize the Module. More...
virtual void	beginRun () override
	Called when entering a new run. More...
virtual void	event () override
	Event processor.
virtual void	terminate () override
	Termination action. More...
bool	isRunningOffsetSubtracted ()
	Checks if running offset is subtracted in TOPDigits. More...
	TOPNtupleModule ()
	Constructor.
virtual	~TOPNtupleModule ()
	Destructor.
virtual void	initialize () override
	Initialize the Module. More...
virtual void	beginRun () override
	Called when entering a new run. More...
virtual void	event () override
	Event processor.
virtual void	endRun () override
	End-of-run action. More...
virtual void	terminate () override
	Termination action. More...
	TOPPackerModule ()
	Constructor.
virtual	~TOPPackerModule ()
	Destructor.
virtual void	initialize () override
	Initialize the Module. More...
virtual void	beginRun () override
	Called when entering a new run. More...
virtual void	event () override
	Event processor.
virtual void	endRun () override
	End-of-run action. More...
virtual void	terminate () override
	Termination action. More...
void	packProductionDraft ()
	Pack in format: c_Draft (tentative production format) this format was never implemented in firmware!
void	packType0Ver16 ()
	Pack in format: c_Type0Ver16 (Feature-extracted data) this format was never implemented in firmware!
void	packProductionDebug ()
	Pack in format: Production Debugging Data Format 01.
	TOPPDFCheckerModule ()
	Constructor.
virtual void	defineHisto () override
	Histogram definitions such as TH1(), TH2(), TNtuple(), TTree().... More...
virtual void	initialize () override
	Initialize the Module. More...
virtual void	event () override
	Event processor.
virtual void	terminate () override
	Termination action. More...
bool	isFromThisParticle (const TOPDigit &digit, const MCParticle *particle)
	Checks if digit comes from given MC particle. More...
	TOPPDFDebuggerModule ()
	Constructor.
virtual void	initialize () override
	Initialize the Module. More...
virtual void	event () override
	Event processor.
void	associatePDFPeaks (const TOP::PDFConstructor &pdfConstructor)
	Associate PDF peaks with photons using S-plot technique. More...
	TOPRawDigitConverterModule ()
	Constructor.
virtual	~TOPRawDigitConverterModule ()
	Destructor.
virtual void	initialize () override
	Initialize the Module. More...
virtual void	beginRun () override
	Called when entering a new run. More...
virtual void	event () override
	Event processor.
virtual void	endRun () override
	End-of-run action. More...
virtual void	terminate () override
	Termination action. More...
	TOPReconstructorModule ()
	Constructor.
virtual void	initialize () override
	Initialize the Module. More...
virtual void	event () override
	Event processor.
	TOPTBCComparatorModule ()
	Constructor.
void	defineHisto () override
	Defining the histograms. More...
int	analyzeCalFile ()
	Analyzes the calibrations stored in the file m_calSetFile. More...
int	makeComparisons ()
	Last function to be called, compared the histograms of different datasets filled by analyzeCalFile() Every new comparison histogram added to the module has to be filled here.
void	initialize () override
	Initialize the module.
void	beginRun () override
	Called when entering a new run.
void	event () override
	Event processor.
void	endRun () override
	End-of-run action. More...
void	terminate () override
	Termination action. More...
int	parseInputDirectoryLine (std::string)
	Utility function to get the directory name and the label from a line of the m_inputDirectoryList file Sets the values of m_calSetDirectory and m_calSetLabel.
int	parseSlotAndScrodIDfromFileName (std::string)
	Utility function to parse the slot and BS id from the calibration file names.
TH1F *	calculateHistoRatio (TH1F *, TH1F *, TH1F *)
	Utility function to take the ratio of two histograms using TH1::Divide(), without overwriting the output name and title initialized in defineHisto().
TH2F *	calculateHistoRatio (TH2F *, TH2F *, TH2F *)
	Utility function to take the ratio of two histograms using TH2::Divide(), without overwriting the output name and title initialized in defineHisto().
	TOPTimeBaseCalibratorModule ()
	Constructor.
virtual	~TOPTimeBaseCalibratorModule ()
	Destructor.
virtual void	initialize () override
	Initialize the Module. More...
virtual void	beginRun () override
	Called when entering a new run. More...
virtual void	event () override
	Event processor.
virtual void	endRun () override
	End-of-run action. More...
virtual void	terminate () override
	Termination action. More...
bool	calibrateChannel (std::vector< TwoTimes > &ntuple, unsigned scrodID, unsigned scrodChannel, TH1F &Hchi2, TH1F &Hndf, TH1F &HDeltaT)
	calibrate single channel More...
bool	matrixInversion (const std::vector< TwoTimes > &ntuple, unsigned scrodID, unsigned scrodChannel, double meanTimeDifference, TH1F &Hchi2, TH1F &Hndf, TH1F &HDeltaT)
	Method by matrix inversion. More...
bool	iterativeTBC (const std::vector< TwoTimes > &ntuple, unsigned scrodID, unsigned scrodChannel, double meanTimeDifference, TH1F &Hchi2, TH1F &Hndf, TH1F &HDeltaT)
	Method by iteration of chi2 minimization. More...
void	Iteration (const std::vector< TwoTimes > &ntuple, std::vector< double > &xval)
	Iteration function called by iterativeTBC() More...
double	Chisq (const std::vector< TwoTimes > &ntuple, const std::vector< double > &xxval)
	Return the chisqure of one set of TBC constants (xval) in iTBC calculaton. More...
void	saveAsHistogram (const std::vector< double > &vec, const std::string &name, const std::string &title, const std::string &xTitle="", const std::string &yTitle="") const
	Save vector to histogram and write it out. More...
void	saveAsHistogram (const std::vector< double > &vec, const std::vector< double > &err, const std::string &name, const std::string &title, const std::string &xTitle="", const std::string &yTitle="") const
	Save vector and errors to histogram and write it out. More...
void	saveAsHistogram (const TMatrixDSym &M, const std::string &name, const std::string &title) const
	Save matrix to histogram and write it out. More...
	TOPTimeRecalibratorModule ()
	Constructor.
virtual void	initialize () override
	Initialize the Module. More...
virtual void	beginRun () override
	Called when entering a new run. More...
virtual void	event () override
	Event processor.
	TOPUnpackerModule ()
	Constructor.
virtual	~TOPUnpackerModule ()
	Destructor.
virtual void	initialize () override
	Initialize the Module. More...
virtual void	beginRun () override
	Called when entering a new run. More...
virtual void	event () override
	Event processor.
std::string	getFrontEndName (RawTOP &raw, int finesse) const
	Returns the name of the front-end. More...
bool	printTheError ()
	Error messages suppression logic. More...
void	unpackProductionDraft (const int *buffer, int bufferSize)
	Unpack raw data given in a tentative production format (will vanish in future) More...
void	unpackType0Ver16 (const int *buffer, int bufferSize)
	Unpack raw data given in feature-extraction production format. More...
bool	unpackHeadersInterimFEVer01 (const int *buffer, int bufferSize, bool swapBytes)
	Tries to unpack raw data assuming it is in feature-extraction interim format. More...
int	unpackInterimFEVer01 (const int *buffer, int bufferSize, bool pedestalSubtracted)
	Unpack raw data given in feature-extraction interim format. More...
int	unpackProdDebug (const int *buffer, int bufferSize, TOP::RawDataType dataFormat, bool pedestalSubtracted)
	Unpack raw data given in production debugging format. More...
virtual void	endRun () override
	End-of-run action. More...
virtual void	terminate () override
	Termination action. More...
	TOPWaveformFeatureExtractorModule ()
	Constructor.
virtual	~TOPWaveformFeatureExtractorModule ()
	Destructor.
virtual void	initialize () override
	Initialize the Module. More...
virtual void	beginRun () override
	Called when entering a new run. More...
virtual void	event () override
	Event processor.
virtual void	endRun () override
	End-of-run action. More...
virtual void	terminate () override
	Termination action. More...
	TOPWaveformQualityPlotterModule ()
	Constructor.
void	defineHisto () override
	Books the empty histograms.
void	initialize () override
	Module initialization, calls defineHisto and gets waveform.
void	basicDebuggingPlots (const TOPRawWaveform &rawwave)
	Fills the debugging 1D histograms and hitmaps. More...
void	drawWaveforms (const TOPRawWaveform &rawwave)
	Draws the full waveforms onto the TProfiles. More...
void	event () override
	Event processor.
void	endRun () override
	End-of-run action. More...

Detailed Description

Function Documentation

analyzeCalFile()

int analyzeCalFile

(

)

Analyzes the calibrations stored in the file m_calSetFile.

This is the main function in which the analysis and the hisogram filling takes place. Every new monitoring histogram added to the module has to be filled here.

Definition at line 207 of file TOPTBCComparatorModule.cc.

  {
    // Here the histograms are filled
    // WARNING! What root calls "RMS" is actually a standard deviation
 
    // Sanity check on the parsed parameters
    if (m_slotID < 0 || m_boardstackID < 0 || m_scrodID < 0) {
      B2WARNING("Negative slot, BS or scrod ID found while calling analyzeCalFile(). Looks like they have not been initialized, or that a function re-initialized them");
      return 0;
    }
 
    // Loop over all the histograms that should be found in the file
    for (short iChannel = 0; iChannel < 128; iChannel++) {
 
      // Pics up one histogram just to check that the channel iChannel was actually used to calculate the calibrations
      if (!m_calSetFile->Get(str(format("timeDiff_ch%1%") % (iChannel)).c_str())) {
        continue;
      }
 
 
      // ---------
      // 1) Define some channel numbering quantities
      // Watchout: m_slotID is in [1..16] so slot m_slotID is on the element m_slotID-1 of the array
      // ---------
 
      short hardwareChannel = iChannel + 128 * m_boardstackID; // 0-511 across the whole module, BS by BS
      auto& chMapper = TOP::TOPGeometryPar::Instance()->getChannelMapper();
      short pixelID = chMapper.getPixelID(hardwareChannel); // 1-512 across the whole module, in rows
      short colNum = (pixelID - 1) % 64 + 1; // 1- 64  column ID (right to left looking from the quartz to the PMTs)
      short rowNum = (pixelID - 1) / 64 + 1 ; // 1- 8 row ID (bottom to top looking from the quartz to the PMTs)
      short globalChannel = hardwareChannel + 512 * (m_slotID -
                                                     1); // channel number across the whole detector, 0-8191. Used for cal monitorin only
 
 
      // ---------
      // 2) Channel-by-channel DeltaT summaries (average on the 256 samples of each set)
      // ---------
 
      // Checks that the histogram needed here is not corrupted before using it
      if (!m_calSetFile->Get(str(format("timeDiffcal_ch%1%") % (iChannel)).c_str())) {
        B2WARNING("Error opening " << str(format("timeDiffcal_ch%1%") % (iChannel)));
      } else {
        TH2F* h_timeDiffcal = (TH2F*)m_calSetFile->Get(str(format("timeDiffcal_ch%1%") % (iChannel)).c_str());
        TH1D* h_projection = h_timeDiffcal->ProjectionY("h_projection", 1, m_numSamples); // full projection
 
        m_slotAverageDeltaT[m_slotID - 1][m_calSetID]->SetBinContent(hardwareChannel + 1, h_projection->GetMean());
        m_slotAverageDeltaT[m_slotID - 1][m_calSetID]->SetBinError(hardwareChannel + 1,
                                                                   h_projection->GetMeanError()); // Do we trust root on this?
        m_slotSigmaDeltaT[m_slotID - 1][m_calSetID]->SetBinContent(hardwareChannel + 1,
                                                                   h_projection->GetRMS()); // WARNING! What root calls "RMS" is actually a standard deviation
        m_slotSigmaDeltaT[m_slotID - 1][m_calSetID]->SetBinError(hardwareChannel + 1,
                                                                 h_projection->GetRMSError()); // WARNING! What root calls "RMS" is actually a standard deviation
 
        m_slotAverageDeltaTMap[m_slotID - 1][m_calSetID]->SetBinContent(colNum, rowNum, h_projection->GetMean());
        m_slotSigmaDeltaTMap[m_slotID - 1][m_calSetID]->SetBinContent(colNum, rowNum,
            h_projection->GetRMS());
 
        m_topAverageDeltaT[m_calSetID]->SetBinContent(globalChannel + 1, h_projection->GetMean());
        m_topSigmaDeltaT[m_calSetID]->SetBinContent(globalChannel + 1, h_projection->GetRMS());
      }
 
      // ---------
      // 3) Channel-by-channel average  occupancy
      // ---------
 
      // Checks that the histogram needed here is not corrupted before using it
      if (!m_calSetFile->Get(str(format("sampleOccup_ch%1%") % (iChannel)).c_str())) {
        B2WARNING("Error opening " << str(format("sampleOccup_ch%1%") % (iChannel)));
      } else {
        TH1F* h_sampleOccup = (TH1F*)m_calSetFile->Get(str(format("sampleOccup_ch%1%") % (iChannel)).c_str());
 
        // reads the occupancy histogram bin-by-by to look for (almost) empty samples
        int nEmpty = 0;
        for (int iSample = 1; iSample < m_numSamples + 1 ; iSample++) {
          if (h_sampleOccup->GetBinContent(iSample) < m_minCalPulses) nEmpty++;
        }
 
        m_slotSampleOccupancy[m_slotID - 1][m_calSetID]->SetBinContent(hardwareChannel + 1,  h_sampleOccup->Integral() / m_numSamples);
        m_slotEmptySamples[m_slotID - 1][m_calSetID]->SetBinContent(hardwareChannel + 1,  nEmpty);
 
        m_slotSampleOccupancyMap[m_slotID - 1][m_calSetID]->SetBinContent(colNum, rowNum,  h_sampleOccup->Integral() / m_numSamples);
        m_slotEmptySamplesMap[m_slotID - 1][m_calSetID]->SetBinContent(colNum, rowNum,  nEmpty);
 
        m_topSampleOccupancy[m_calSetID]->SetBinContent(globalChannel + 1, h_sampleOccup->Integral() / m_numSamples);
 
      }
 
 
 
    }
    return 1;
  }

associatePDFPeaks()

void associatePDFPeaks ( const TOP::PDFConstructor &  pdfConstructor )

private

Associate PDF peaks with photons using S-plot technique.

Parameters

pdfConstructor

reconstruction object for given track and hypothesis

Definition at line 192 of file TOPPDFDebuggerModule.cc.

  {
 
    const auto& signalPDFs = pdfConstructor.getSignalPDF();
    auto signalPhotons = pdfConstructor.getExpectedSignalPhotons();
    auto bkgPhotons = pdfConstructor.getExpectedBkgPhotons();
    auto deltaPhotons = pdfConstructor.getExpectedDeltaPhotons();
    int PDGCode = pdfConstructor.getHypothesis().getPDGCode();
 
    for (const auto& digit : m_digits) {
      if (digit.getModuleID() != pdfConstructor.getModuleID()) continue;
      if (digit.getHitQuality() != TOPDigit::c_Good) continue;
      int pixelID = digit.getPixelID();
      unsigned index = pixelID - 1;
      if (index >= signalPDFs.size()) continue;
      const auto& signalPDF = signalPDFs[index];
      const auto* tts = signalPDF.getTTS();
      const auto& pdfPeaks = signalPDF.getPDFPeaks();
      const auto& extraInfos = signalPDF.getPDFExtraInfo();
      if (extraInfos.size() != pdfPeaks.size()) {
        B2ERROR("associatePDFPeaks: sizes of PDFPeaks and ExtraInfo's differ");
        continue;
      }
 
      auto* associatedPDF = m_associatedPDFs.appendNew(PDGCode);
      digit.addRelationTo(associatedPDF);
 
      double bkgLevel = pdfConstructor.getBackgroundPDF()->getPDFValue(pixelID) * bkgPhotons;
      associatedPDF->setBackgroundWeight(bkgLevel);
      double deltaLevel = pdfConstructor.getDeltaRayPDF().getPDFValue(pixelID, digit.getTime()) * deltaPhotons;
      associatedPDF->setDeltaRayWeight(deltaLevel);
 
      float time = digit.getTime();
      float timeErr = digit.getTimeError();
 
      for (size_t k = 0; k < pdfPeaks.size(); k++) {
        const auto& pdfPeak = pdfPeaks[k];
        const auto& pdfExtra = extraInfos[k];
        TOPAssociatedPDF::PDFPeak peak;
        peak.position = pdfPeak.t0;
        peak.width = sqrt(pdfPeak.wid);
        peak.numPhotons = pdfPeak.nph * signalPhotons;
        float wt = 0;
        for (const auto& gaus : tts->getTTS()) {
          float sig2 = peak.width * peak.width + gaus.sigma * gaus.sigma + timeErr * timeErr;
          float x = pow(time - peak.position - gaus.position, 2) / sig2;
          if (x > 20) continue;
          wt += peak.numPhotons * gaus.fraction / sqrt(2 * M_PI * sig2) * exp(-x / 2);
        }
        if (wt > 0) {
          peak.fic = pdfPeak.fic;
          peak.e = pdfExtra.e;
          peak.sige = pdfExtra.sige;
          peak.nx = abs(pdfExtra.Nxm) + abs(pdfExtra.Nxb) + abs(pdfExtra.Nxe);
          peak.ny = abs(pdfExtra.Nym) + abs(pdfExtra.Nyb) + abs(pdfExtra.Nye);
          peak.nxm = abs(pdfExtra.Nxm);
          peak.nym = abs(pdfExtra.Nym);
          peak.nxe = abs(pdfExtra.Nxe);
          peak.nye = abs(pdfExtra.Nye);
          peak.xd = pdfExtra.xD;
          peak.yd = pdfExtra.yD;
          peak.type = pdfExtra.type;
          peak.kxe = pdfExtra.kxE;
          peak.kye = pdfExtra.kyE;
          peak.kze = pdfExtra.kzE;
          peak.kxd = pdfExtra.kxD;
          peak.kyd = pdfExtra.kyD;
          peak.kzd = pdfExtra.kzD;
          associatedPDF->appendPeak(peak, wt);
        }
      }
    }
 
  }

basicDebuggingPlots()

void basicDebuggingPlots

(

const TOPRawWaveform &

rawwave

)

Fills the debugging 1D histograms and hitmaps.

Parameters

rawwave

the raw waveform

Definition at line 81 of file TOPWaveformQualityPlotterModule.cc.

  {
    int scrodid = v.getScrodID();
    int asicid = v.getASICNumber();
    int channelid = v.getASICChannel();
    int carrierid = v.getCarrierNumber();
    m_scrod_id->Fill(scrodid);
    m_asic_ch->Fill(channelid);
    m_asic->Fill(asicid);
    m_carrier->Fill(carrierid);
    m_asic_win->Fill(v.getStorageWindow());
    m_entries->Fill(v.getWaveform().size());
    m_moduleID->Fill(v.getModuleID());
    m_pixelID->Fill(v.getPixelID());
 
    if (m_hitmap.find(scrodid) == m_hitmap.end()) {
      m_hitmap[scrodid] = new TH2F((string("scrod ") + to_string(scrodid) + string("Hitmap")).c_str(),
                                   (string("scrod ") + to_string(scrodid) + string("carrier vs. asic;asic;carrier")).c_str(), 4, 0, 4, 4, 0, 4);
    }
    m_hitmap[scrodid]->Fill(asicid, carrierid);
    const vector<short>& waveform = v.getWaveform();
    if (waveform.empty()) {
      return;
    }
    for (short adc : waveform) {
      m_samples->Fill(adc);
    }
  }

beginRun() [1/16]

void beginRun ( void  )

overridevirtual

Called when entering a new run.

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

Reimplemented from Module.

Definition at line 195 of file TOPBackgroundModule.cc.

  {
    // Print run number
    B2INFO("TOPBackground: Processing:");
 
  }

beginRun() [2/16]

void beginRun ( void  )

overridevirtual

Called when entering a new run.

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

Reimplemented from Module.

Definition at line 54 of file TOPChannelMaskerModule.cc.

  {
 
    if (not m_channelMask.isValid()) {
      B2FATAL("channel mask not available");
    }
    if (not m_channelT0.isValid()) {
      B2FATAL("channel T0 calibration not available");
    }
    if (not m_timebase.isValid()) {
      B2FATAL("timebase calibration not available");
    }
 
  }

beginRun() [3/16]

void beginRun ( void  )

overridevirtual

Called when entering a new run.

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

Reimplemented from Module.

Definition at line 148 of file TOPDoublePulseGeneratorModule.cc.

  {
    StoreObjPtr<EventMetaData> evtMetaData;
    if (m_useDatabase) {
      if (!(*m_timebase).isValid()) {
        B2FATAL("Sample time calibration requested but not available for run "
                << evtMetaData->getRun()
                << " of experiment " << evtMetaData->getExperiment());
      }
    }
 
  }

beginRun() [4/16]

void beginRun ( void  )

overridevirtual

Called when entering a new run.

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

Reimplemented from HistoModule.

Definition at line 307 of file TOPDQMModule.cc.

  {
    m_BoolEvtMonitor->Reset();
    m_window_vs_slot->Reset();
    m_eventT0->Reset();
    m_bunchOffset->Reset();
    m_time->Reset();
    m_timeBG->Reset();
    m_signalHits->Reset();
    m_backgroundHits->Reset();
    m_trackHits->Reset();
    m_goodTDCAll->Reset();
    m_badTDCAll->Reset();
    m_goodHitsPerEventAll->Reset();
    m_badHitsPerEventAll->Reset();
    m_TOPOccAfterInjLER->Reset();
    m_TOPOccAfterInjHER->Reset();
    m_TOPEOccAfterInjLER->Reset();
    m_TOPEOccAfterInjHER->Reset();
    m_nhitInjLER->Reset();
    m_nhitInjHER->Reset();
    m_nhitInjLERcut->Reset();
    m_nhitInjHERcut->Reset();
    m_eventInjLER->Reset();
    m_eventInjHER->Reset();
    m_eventInjLERcut->Reset();
    m_eventInjHERcut->Reset();
 
    for (int i = 0; i < m_numModules; i++) {
      m_window_vs_asic[i]->Reset();
      m_goodHitsXY[i]->Reset();
      m_badHitsXY[i]->Reset();
      m_goodHitsAsics[i]->Reset();
      m_badHitsAsics[i]->Reset();
      m_goodTDC[i]->Reset();
      m_badTDC[i]->Reset();
      m_goodTiming[i]->Reset();
      m_goodTimingBG[i]->Reset();
      m_goodChannelHits[i]->Reset();
      m_badChannelHits[i]->Reset();
      m_pulseHeights[i]->Reset();
    }
  }

beginRun() [5/16]

void beginRun ( void  )

overridevirtual

Called when entering a new run.

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

Reimplemented from HistoModule.

Definition at line 150 of file TOPLaserHitSelectorModule.cc.

  {
  }

beginRun() [6/16]

void beginRun ( void  )

overridevirtual

Called when entering a new run.

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

Reimplemented from HistoModule.

Definition at line 159 of file TOPInterimFENtupleModule.cc.

  {
  }

beginRun() [7/16]

void beginRun ( void  )

overridevirtual

Called when entering a new run.

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

Reimplemented from Module.

Definition at line 81 of file TOPLaserCalibratorModule.cc.

  {
  }

beginRun() [8/16]

void beginRun ( void  )

overridevirtual

Called when entering a new run.

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

Reimplemented from Module.

Definition at line 88 of file TOPMCTrackMakerModule.cc.

  {
  }

beginRun() [9/16]

void beginRun ( void  )

overridevirtual

Called when entering a new run.

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

Reimplemented from Module.

Definition at line 160 of file TOPModuleT0CalibratorModule.cc.

  {
    if (m_moduleT0.hasChanged()) {
      StoreObjPtr<EventMetaData> evtMetaData;
      B2WARNING("ModuleT0 payload has changed. Calibration results may not be correct."
                << LogVar("run", evtMetaData->getRun()));
    }
  }

beginRun() [10/16]

void beginRun ( void  )

overridevirtual

Called when entering a new run.

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

Reimplemented from Module.

Definition at line 122 of file TOPNtupleModule.cc.

  {
  }

beginRun() [11/16]

void beginRun ( void  )

overridevirtual

Called when entering a new run.

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

Reimplemented from Module.

Definition at line 99 of file TOPPackerModule.cc.

  {
  }

beginRun() [12/16]

void beginRun ( void  )

overridevirtual

Called when entering a new run.

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

Reimplemented from Module.

Definition at line 132 of file TOPRawDigitConverterModule.cc.

  {
 
    StoreObjPtr<EventMetaData> evtMetaData;
 
    // check if calibrations are available when needed - if not, terminate
 
    if (m_useSampleTimeCalibration) {
      if (not m_timebase.isValid()) {
        B2FATAL("Sample time calibration requested but not available for run "
                << evtMetaData->getRun()
                << " of experiment " << evtMetaData->getExperiment());
      }
    }
    if (m_useChannelT0Calibration) {
      if (not m_channelT0.isValid()) {
        B2FATAL("Channel T0 calibration requested but not available for run "
                << evtMetaData->getRun()
                << " of experiment " << evtMetaData->getExperiment());
      }
    }
    if (m_useAsicShiftCalibration) {
      if (not m_asicShift.isValid()) {
        B2FATAL("ASIC shifts calibration requested but not available for run "
                << evtMetaData->getRun()
                << " of experiment " << evtMetaData->getExperiment());
      }
    }
    if (m_useModuleT0Calibration) {
      if (not m_moduleT0.isValid()) {
        B2FATAL("Module T0 calibration requested but not available for run "
                << evtMetaData->getRun()
                << " of experiment " << evtMetaData->getExperiment());
      }
    }
    if (m_useCommonT0Calibration) {
      if (not m_commonT0.isValid()) {
        B2FATAL("Common T0 calibration requested but not available for run "
                << evtMetaData->getRun()
                << " of experiment " << evtMetaData->getExperiment());
      }
    }
    if (not m_timeWalk.isValid()) {
      // B2FATAL("Time-walk calibration is not available for run "
      B2WARNING("Time-walk calibration is not available for run "
                << evtMetaData->getRun()
                << " of experiment " << evtMetaData->getExperiment());
    }
    if (m_pedestalRMS > 0 and not m_noises.isValid()) {
      B2FATAL("Channel noise levels are not available for run "
              << evtMetaData->getRun()
              << " of experiment " << evtMetaData->getExperiment());
    }
 
    if (not m_feSetting.isValid()) {
      B2FATAL("Front-end settings are not available for run "
              << evtMetaData->getRun()
              << " of experiment " << evtMetaData->getExperiment());
    }
 
  }

beginRun() [13/16]

void beginRun ( void  )

overridevirtual

Called when entering a new run.

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

Reimplemented from Module.

Definition at line 143 of file TOPTimeBaseCalibratorModule.cc.

  {
  }

beginRun() [14/16]

void beginRun ( void  )

overridevirtual

Called when entering a new run.

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

Reimplemented from Module.

Definition at line 85 of file TOPTimeRecalibratorModule.cc.

  {
    StoreObjPtr<EventMetaData> evtMetaData;
 
    // check if calibrations are available when needed - if not, terminate
 
    if (m_useSampleTimeCalibration) {
      if (not m_timebase.isValid()) {
        B2FATAL("Sample time calibration requested but not available for run "
                << evtMetaData->getRun()
                << " of experiment " << evtMetaData->getExperiment());
      }
    }
    if (m_useChannelT0Calibration) {
      if (not m_channelT0.isValid()) {
        B2FATAL("Channel T0 calibration requested but not available for run "
                << evtMetaData->getRun()
                << " of experiment " << evtMetaData->getExperiment());
      }
    }
    if (m_useAsicShiftCalibration) {
      if (not m_asicShift.isValid()) {
        B2FATAL("ASIC shifts calibration requested but not available for run "
                << evtMetaData->getRun()
                << " of experiment " << evtMetaData->getExperiment());
      }
    }
    if (m_useModuleT0Calibration) {
      if (not m_moduleT0.isValid()) {
        B2FATAL("Module T0 calibration requested but not available for run "
                << evtMetaData->getRun()
                << " of experiment " << evtMetaData->getExperiment());
      }
    }
    if (m_useCommonT0Calibration) {
      if (not m_commonT0.isValid()) {
        B2FATAL("Common T0 calibration requested but not available for run "
                << evtMetaData->getRun()
                << " of experiment " << evtMetaData->getExperiment());
      }
    }
    if (m_useTimeWalkCalibration) {
      if (not m_timeWalk.isValid()) {
        // B2FATAL("Time-walk calibration requested but not available for run "
        B2WARNING("Time-walk calibration is not available for run "
                  << evtMetaData->getRun()
                  << " of experiment " << evtMetaData->getExperiment());
      }
    }
 
    if (not m_feSetting.isValid()) {
      B2FATAL("Front-end settings are not available for run "
              << evtMetaData->getRun()
              << " of experiment " << evtMetaData->getExperiment());
    }
 
  }

beginRun() [15/16]

void beginRun ( void  )

overridevirtual

Called when entering a new run.

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

Reimplemented from Module.

Definition at line 113 of file TOPUnpackerModule.cc.

  {
    m_channelStatistics.clear();
  }

beginRun() [16/16]

void beginRun ( void  )

overridevirtual

Called when entering a new run.

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

Reimplemented from Module.

Definition at line 78 of file TOPWaveformFeatureExtractorModule.cc.

  {
  }

calibrateChannel()

bool calibrateChannel ( std::vector< TwoTimes > &  ntuple, unsigned  scrodID, unsigned  scrodChannel, TH1F &  Hchi2, TH1F &  Hndf, TH1F &  HDeltaT  )

private

calibrate single channel

Parameters

ntuple	ntuple data
scrodID	SCROD ID
scrodChannel	channel number within SCROD (0 - 127)
Hchi2	histogram to store normalized chi^2
Hndf	histogram to store degrees of freedom
HDeltaT	histogram to store fittet double pulse delay

Returns

true on success

Definition at line 311 of file TOPTimeBaseCalibratorModule.cc.

Chisq()

double Chisq ( const std::vector< TwoTimes > &  ntuple, const std::vector< double > &  xxval  )

private

Return the chisqure of one set of TBC constants (xval) in iTBC calculaton.

Parameters

ntuple	ntuple data
xxval	TBC constants of 256 samples, and xxval[0]=0, xxval[256]=2*FTSW

Definition at line 669 of file TOPTimeBaseCalibratorModule.cc.

closeRun()

void closeRun ( )

finalprivatevirtual

Replacement for endRun().

Do anything you would normally do in endRun here.

Reimplemented from CalibrationCollectorModule.

Definition at line 288 of file TOPValidationCollectorModule.cc.

  {
 
    // module T0 pulls
 
    for (int module = 0; module < c_numModules; module++) {
      std::vector<double> pos, err;
      for (int set = 0; set < c_numSets; set++) {
        const auto& minimum = m_finders[set][module].getMinimum();
        if (minimum.valid) {
          pos.push_back(minimum.position);
          err.push_back(minimum.error);
        }
      }
      auto h_pulls = getObjectPtr<TH1F>("moduleT0_pulls");
      for (unsigned i = 0; i < pos.size(); i++) {
        for (unsigned j = i + 1; j < pos.size(); j++) {
          double pull = (pos[i] - pos[j]) / sqrt(err[i] * err[i] + err[j] * err[j]);
          h_pulls->Fill(pull);
        }
      }
    }
 
    // module T0 residuals
 
    for (int module = 0; module < c_numModules; module++) {
      auto& finder = m_finders[0][module];
      for (int set = 1; set < c_numSets; set++) {
        finder.add(m_finders[set][module]);
      }
      const auto& minimum = finder.getMinimum();
      if (minimum.valid) {
        m_treeEntry.moduleT0[module] = minimum.position;
        m_treeEntry.moduleT0Err[module] = minimum.error;
      }
    }
 
    // common T0 residual
 
    auto& finder = m_finders[0][0];
    for (int module = 1; module < c_numModules; module++) {
      finder.add(m_finders[0][module]);
    }
    const auto& minimum = finder.getMinimum();
    if (minimum.valid) {
      m_treeEntry.commonT0 = minimum.position;
      m_treeEntry.commonT0Err = minimum.error;
    }
 
    // fill the tree
 
    auto tree = getObjectPtr<TTree>("tree");
    tree->Fill();
  }

collect() [1/11]

void collect ( )

finalprivatevirtual

Replacement for event().

Fill your calibration data objects here

Reimplemented from CalibrationCollectorModule.

Definition at line 216 of file TOPAlignmentCollectorModule.cc.

  {
    // bunch must be reconstructed
 
    if (not m_recBunch.isValid()) return;
    if (not m_recBunch->isReconstructed()) return;
 
    // track-by-track iterations
 
    for (const auto& track : m_tracks) {
 
      // construct TOPTrack from mdst track
      TOPTrack trk(track);
      if (not trk.isValid()) continue;
 
      // skip if track not hitting target module
      if (trk.getModuleID() != m_targetMid) continue;
 
      // track selection
      if (not m_selector.isSelected(trk)) continue;
 
      // generate sub-sample number
      int sub = gRandom->Integer(c_numSets);
      auto& align = m_align[sub];
      auto& countFails = m_countFails[sub];
      const auto& name = m_treeNames[sub];
      auto h1 = getObjectPtr<TH2F>("tracks_per_slot");
      h1->Fill(trk.getModuleID(), sub);
 
      // do an iteration
      int err = align.iterate(trk, m_selector.getChargedStable());
      m_iter++;
 
      // check number of consecutive failures, and in case reset
      if (err == 0) {
        countFails = 0;
      } else if (countFails <= m_maxFails) {
        countFails++;
      } else {
        B2INFO("Reached maximum allowed number of failed iterations. "
               "Resetting TOPalign object of set = " << sub << " at iter = " << m_iter);
        align.reset();
        countFails = 0;
      }
 
      // get new parameter values and estimated errors
      m_vAlignPars = align.getParameters();
      m_vAlignParsErr = align.getErrors();
      m_ntrk = align.getNumUsedTracks();
      m_errorCode = err;
      m_valid = align.isValid();
      m_numPhot = align.getNumOfPhotons();
 
      // set other ntuple variables
      const auto& localPosition = m_selector.getLocalPosition();
      m_x = localPosition.X();
      m_y = localPosition.Y();
      m_z = localPosition.Z();
      const auto& localMomentum = m_selector.getLocalMomentum();
      m_p = localMomentum.R();
      m_theta = localMomentum.Theta();
      m_phi = localMomentum.Phi();
      const auto& pocaPosition = m_selector.getPOCAPosition();
      m_pocaR = pocaPosition.Rho();
      m_pocaZ = pocaPosition.Z();
      m_pocaX = pocaPosition.X();
      m_pocaY = pocaPosition.Y();
      m_cmsE = m_selector.getCMSEnergy();
      m_charge = trk.getCharge();
      m_PDG = trk.getPDGCode();
 
      // fill output tree
      auto alignTree = getObjectPtr<TTree>(name);
      alignTree->Fill();
 
      // fill control histograms
      std::string slotName = "_s" + to_string(m_targetMid);
      auto h2 = getObjectPtr<TH1F>("local_z" + slotName);
      h2->Fill(m_z);
      auto h3 = getObjectPtr<TH2F>("cth_vs_p" + slotName);
      h3->Fill(m_p, cos(m_theta));
      auto h4 = getObjectPtr<TH2F>("poca_xy" + slotName);
      h4->Fill(m_pocaX, m_pocaY);
      auto h5 = getObjectPtr<TH1F>("poca_z" + slotName);
      h5->Fill(m_pocaZ);
      auto h6 = getObjectPtr<TH1F>("Ecms" + slotName);
      h6->Fill(m_cmsE);
      auto h7 = getObjectPtr<TH1F>("charge" + slotName);
      h7->Fill(m_charge);
      auto h8 = getObjectPtr<TH2F>("timeHits" + slotName);
      for (const auto& digit : m_digits) {
        if (digit.getHitQuality() != TOPDigit::c_Good) continue;
        if (digit.getModuleID() != m_targetMid) continue;
        h8->Fill(digit.getChannel(), digit.getTime());
      }
      auto h9 = getObjectPtr<TH1F>("numPhot" + slotName);
      h9->Fill(m_numPhot);
 
    } // tracks
 
  }

collect() [2/11]

void collect ( )

finalprivatevirtual

Replacement for event().

Fill your calibration data objects here.

Reimplemented from CalibrationCollectorModule.

Definition at line 91 of file TOPAsicShiftsBS13dCollectorModule.cc.

  {
 
    if (m_requireRecBunch) {
      if (not m_recBunch.isValid()) return;
      if (not m_recBunch->isReconstructed()) return;
    }
 
    auto time_vs_BS = getObjectPtr<TH2F>("time_vs_BS");
    auto timeReference = getObjectPtr<TH1F>("time_reference");
    std::vector<TH1F*> timeCarriers;
    for (unsigned i = 0; i < 4; i++) {
      string name = "time_carr_" + to_string(i);
      auto h = getObjectPtr<TH1F>(name);
      timeCarriers.push_back(h);
    }
 
    for (const auto& digit : m_topDigits) {
      if (digit.getModuleID() != 13) continue;
      if (digit.getHitQuality() != TOPDigit::c_Good) continue;
      double time = digit.getTime() - m_timeOffset;
      time_vs_BS->Fill(digit.getChannel(), time);
      if (digit.getBoardstackNumber() < 3) {
        timeReference->Fill(time);
      } else {
        auto c = digit.getCarrierNumber();
        timeCarriers[c]->Fill(time);
      }
    }
 
  }

collect() [3/11]

void collect ( )

finalprivatevirtual

Replacement for event().

Fill your calibration data objects here

Reimplemented from CalibrationCollectorModule.

Definition at line 70 of file TOPChannelMaskCollectorModule.cc.

  {
 
    std::vector<TH1F*> slots;
    for (const auto& name : m_names) {
      auto h = getObjectPtr<TH1F>(name);
      slots.push_back(h);
    }
 
    std::vector<TH2F*> asics;
    for (const auto& name : m_asicNames) {
      auto h = getObjectPtr<TH2F>(name);
      asics.push_back(h);
    }
 
    int NGood = 0;
    for (const auto& digit : m_digits) {
      if (digit.getHitQuality() == TOPDigit::c_Junk) continue;
      NGood++;
      unsigned m = digit.getModuleID() - 1;
      if (m < slots.size()) slots[m]->Fill(digit.getChannel());
      if (m < asics.size()) asics[m]->Fill(digit.getChannel() / 8, digit.getRawTime() / 64 + 220);
    }
 
    auto h = getObjectPtr<TH1F>("nhits");
    h->Fill(NGood);
 
  }

collect() [4/11]

void collect ( )

finalprivatevirtual

Replacement for event().

Fill your calibration data objects here

Reimplemented from CalibrationCollectorModule.

Definition at line 68 of file TOPCommonT0BFCollectorModule.cc.

  {
 
    if (not m_recBunch.isValid()) return;
    if (not m_recBunch->isReconstructed()) return;
    if (m_recBunch->getNumTracks() != 2) return;
 
    auto h1a = getObjectPtr<TH1F>("offset_a");
    auto h1b = getObjectPtr<TH1F>("offset_b");
    auto offset = m_recBunch->getCurrentOffset();
    if (m_commonT0->isCalibrated()) offset += m_commonT0->getT0();
 
    // wrap-around into [-1/2, 1/2] of bunch cycle
    auto a = offset - round(offset / m_bunchTimeSep) * m_bunchTimeSep;
    h1a->Fill(a);
 
    // wrap-around into [0, 1] of bunch cycle
    auto b = offset - round(offset / m_bunchTimeSep - 0.5) * m_bunchTimeSep;
    h1b->Fill(b);
 
  }

collect() [5/11]

void collect ( )

finalprivatevirtual

Replacement for event().

Fill your calibration data objects here

Reimplemented from CalibrationCollectorModule.

Definition at line 148 of file TOPCommonT0LLCollectorModule.cc.

  {
    // bunch must be reconstructed
 
    if (not m_recBunch.isValid()) return;
    if (not m_recBunch->isReconstructed()) return;
 
    TOPRecoManager::setDefaultTimeWindow();
    double timeMin = TOPRecoManager::getMinTime();
    double timeMax = TOPRecoManager::getMaxTime();
 
    // correct-back digits for common T0
 
    double T0 = 0;
    if (m_commonT0->isCalibrated()) T0 = m_commonT0->getT0();
    for (auto& digit : m_digits) {
      if (digit.isCommonT0Calibrated()) digit.subtractT0(-T0);
      if (digit.hasStatus(TOPDigit::c_BunchOffsetSubtracted)) digit.subtractT0(-m_recBunch->getAverageOffset());
    }
 
    // loop over reconstructed tracks, make a selection and accumulate log likelihoods
 
    int ntra = 0;
    for (const auto& track : m_tracks) {
 
      // track selection
      TOPTrack trk(track);
      if (not trk.isValid()) continue;
 
      if (not m_selector.isSelected(trk)) continue;
 
      // construct PDF
      PDFConstructor pdfConstructor(trk, m_selector.getChargedStable(), m_PDFOption);
      if (not pdfConstructor.isValid()) continue;
 
      // minimization procedure: accumulate
      int sub = gRandom->Integer(c_numSets); // generate sub-sample number
      auto h = getObjectPtr<TH1D>(m_names[sub]);
      for (int ibin = 0; ibin < h->GetNbinsX(); ibin++) {
        double t0 = h->GetBinCenter(ibin + 1);
        double chi = h->GetBinContent(ibin + 1);
        chi += -2 * pdfConstructor.getLogL(t0, m_sigmaSmear).logL;
        h->SetBinContent(ibin + 1, chi);
      }
      auto h1 = getObjectPtr<TH1F>("tracks_per_set");
      h1->Fill(sub);
 
      // fill histograms of hits
      auto h2 = getObjectPtr<TH1F>("numHits");
      auto h3 = getObjectPtr<TH2F>("timeHits");
      for (const auto& digit : m_digits) {
        if (digit.getHitQuality() != TOPDigit::c_Good) continue;
        if (digit.getModuleID() != trk.getModuleID()) continue;
        if (digit.getTime() < timeMin) continue;
        if (digit.getTime() > timeMax) continue;
        h2->Fill(digit.getModuleID());
        int bs = digit.getBoardstackNumber();
        h3->Fill((digit.getModuleID() * 4 + bs - 1.5) / 4.0, digit.getTime());
      }
      ntra++;
    }
 
    // fill just another control histogram
 
    if (ntra > 0) {
      auto h4 = getObjectPtr<TH1F>("offset");
      h4->Fill(m_recBunch->getCurrentOffset());
    }
 
  }

collect() [6/11]

void collect ( )

finalprivatevirtual

Replacement for event().

Fill your calibration data objects here

Reimplemented from CalibrationCollectorModule.

Definition at line 78 of file TOPModuleT0DeltaTCollectorModule.cc.

  {
 
    if (m_timeZeros.getEntries() != 2) return;
 
    const auto* timeZero1 = m_timeZeros[0];
    const auto* timeZero2 = m_timeZeros[1];
    if (timeZero1->getModuleID() > timeZero2->getModuleID()) {
      timeZero1 = m_timeZeros[1];
      timeZero2 = m_timeZeros[0];
    }
 
    int slot1 = timeZero1->getModuleID();
    int slot2 = timeZero2->getModuleID();
    if (slot1 > 9 or slot2 - slot1 < 7 or slot2 - slot1 > 9) return;
    string name = "deltaT0_" + to_string(slot1) + "-" + to_string(slot2);
    auto h = getObjectPtr<TH1F>(name);
    if (not h) return;
 
    double t1 = timeZero1->getTime();
    if (m_moduleT0->isCalibrated(slot1)) t1 += m_moduleT0->getT0(slot1);
    double t2 = timeZero2->getTime();
    if (m_moduleT0->isCalibrated(slot2)) t2 += m_moduleT0->getT0(slot2);
    h->Fill(t1 - t2);
 
    auto slotPairs = getObjectPtr<TH2F>("slots");
    slotPairs->Fill(slot1, slot2);
  }

collect() [7/11]

void collect ( )

finalprivatevirtual

Replacement for event().

Fill your calibration data objects here

Reimplemented from CalibrationCollectorModule.

Definition at line 150 of file TOPModuleT0LLCollectorModule.cc.

  {
    // bunch must be reconstructed
 
    if (not m_recBunch.isValid()) return;
    if (not m_recBunch->isReconstructed()) return;
 
    TOPRecoManager::setDefaultTimeWindow();
    double timeMin = TOPRecoManager::getMinTime();
    double timeMax = TOPRecoManager::getMaxTime();
 
    // correct-back digits for module T0
 
    for (auto& digit : m_digits) {
      int slot = digit.getModuleID();
      if (digit.isModuleT0Calibrated()) digit.subtractT0(-m_moduleT0->getT0(slot));
      if (digit.hasStatus(TOPDigit::c_BunchOffsetSubtracted)) digit.subtractT0(-m_recBunch->getAverageOffset());
    }
 
    // loop over reconstructed tracks, make a selection and accumulate log likelihoods
 
    int ntra = 0;
    for (const auto& track : m_tracks) {
 
      // track selection
      TOPTrack trk(track);
      if (not trk.isValid()) continue;
 
      if (not m_selector.isSelected(trk)) continue;
 
      // construct PDF
      PDFConstructor pdfConstructor(trk, m_selector.getChargedStable(), m_PDFOption);
      if (not pdfConstructor.isValid()) continue;
 
      // minimization procedure: accumulate
      int sub = gRandom->Integer(c_numSets); // generate sub-sample number
      unsigned module = trk.getModuleID() - 1;
      if (module >= m_names[sub].size()) continue;
      auto h = getObjectPtr<TH1D>(m_names[sub][module]);
      for (int ibin = 0; ibin < h->GetNbinsX(); ibin++) {
        double t0 = h->GetBinCenter(ibin + 1);
        double chi = h->GetBinContent(ibin + 1);
        chi += -2 * pdfConstructor.getLogL(t0, m_sigmaSmear).logL;
        h->SetBinContent(ibin + 1, chi);
      }
      auto h1 = getObjectPtr<TH2F>("tracks_per_slot");
      h1->Fill(trk.getModuleID(), sub);
 
      // fill histograms of hits
      auto h2 = getObjectPtr<TH1F>("numHits");
      auto h3 = getObjectPtr<TH2F>("timeHits");
      for (const auto& digit : m_digits) {
        if (digit.getHitQuality() != TOPDigit::c_Good) continue;
        if (digit.getModuleID() != trk.getModuleID()) continue;
        if (digit.getTime() < timeMin) continue;
        if (digit.getTime() > timeMax) continue;
        h2->Fill(digit.getModuleID());
        int bs = digit.getBoardstackNumber();
        h3->Fill((digit.getModuleID() * 4 + bs - 1.5) / 4.0, digit.getTime());
      }
      ntra++;
    }
 
    // fill just another control histogram
 
    if (ntra > 0) {
      auto h4 = getObjectPtr<TH1F>("offset");
      h4->Fill(m_recBunch->getCurrentOffset());
    }
 
  }

collect() [8/11]

void collect ( )

finalprivatevirtual

Replacement for event().

Fill your calibration data objects here

Reimplemented from CalibrationCollectorModule.

Definition at line 61 of file TOPOffsetCollectorModule.cc.

  {
    if (m_firstEvent) {
      m_firstEvent = false;
      if (m_bunchStructure->isSet()) {
        auto h = getObjectPtr<TH1F>("fillPattern");
        int RFBuckets = m_bunchStructure->getRFBucketsPerRevolution();
        for (int i = 0; i < RFBuckets; i++) {
          if (m_bunchStructure->getBucket(i)) h->SetBinContent(i + 1, 1);
        }
      }
    }
 
    if (not m_recBunch->isReconstructed()) return;
 
    for (const auto& x : m_names) {
      const auto& detector = x.first;
      if (m_eventT0->hasTemporaryEventT0(detector)) {
        auto eventT0s = m_eventT0->getTemporaryEventT0s(detector);
        if (eventT0s.empty()) continue;
        if (detector == Const::CDC and eventT0s.back().algorithm != "chi2") continue;
        double t0 = eventT0s.back().eventT0;
        auto h = getObjectPtr<TH1F>(x.second);
        h->Fill(t0 - m_recBunch->getTime());
      }
    }
 
    auto h = getObjectPtr<TH1F>("recBuckets");
    h->Fill(m_recBunch->getBucketNumber(0, m_bunchStructure->getRFBucketsPerRevolution()));
  }

collect() [9/11]

void collect ( )

finalprivatevirtual

Replacement for event().

Fill your calibration data objects here

Reimplemented from CalibrationCollectorModule.

Definition at line 171 of file TOPPhotonYieldsCollectorModule.cc.

  {
    // bunch must be reconstructed
 
    if (not m_recBunch->isReconstructed()) return;
 
    // loop over reconstructed tracks, make a selection and fill histograms
 
    for (const auto& track : m_tracks) {
      // track selection
      TOPTrack trk(track);
      if (not trk.isValid()) continue;
      if (not m_selector.isSelected(trk)) continue;
 
      // fill histograms
      auto timeStamp = getObjectPtr<TProfile>("timeStamp");
      timeStamp->Fill(0.5, m_eventMetaData->getTime() / 1000000000);
 
      int slot = trk.getModuleID();
      auto numTracks = getObjectPtr<TH1F>("numTracks");
      numTracks->Fill(slot);
 
      auto muonZ = getObjectPtr<TH1F>(m_muonZNames[slot - 1]);
      muonZ->Fill(m_selector.getLocalPosition().Z());
 
      auto signalHits = getObjectPtr<TH1F>(m_signalNames[slot - 1]);
      auto bkgHits = getObjectPtr<TH1F>(m_bkgNames[slot - 1]);
      auto pulseHeight = getObjectPtr<TH2F>(m_pulseHeightNames[slot - 1]);
      for (const auto& digit : m_digits) {
        if (digit.getModuleID() != slot) continue;
        if (digit.getHitQuality() != TOPDigit::c_Good) continue;  // junk hit or pixel masked-out
        if (not m_thresholdEff->isCalibrated(slot, digit.getChannel())) continue;  // threshold effi. not calibrated
        double effi = m_thresholdEff->getThrEff(slot, digit.getChannel());
        if (effi < m_minThresholdEffi) continue; // to suppress possibly unreliable calibration
        if (std::abs(digit.getTime()) > m_timeWindow) continue;
        // fill signal and background hits with weight=1/effi to correct for threshold efficiency
        if (digit.getTime() > 0) {
          signalHits->Fill(digit.getPixelID(), 1 / effi);
          pulseHeight->Fill(digit.getPixelID(), digit.getPulseHeight());
        } else {
          bkgHits->Fill(digit.getPixelID(), 1 / effi);
        }
      }
 
      auto activePixels = getObjectPtr<TH1F>(m_activeNames[slot - 1]);
      const auto& chMapper = TOPGeometryPar::Instance()->getChannelMapper();
      for (int pixel = 1; pixel <= activePixels->GetNbinsX(); pixel++) {
        unsigned channel = chMapper.getChannel(pixel);
        if (not m_thresholdEff->isCalibrated(slot, channel)) continue; // pixel excluded in counting hits
        if (m_thresholdEff->getThrEff(slot, channel) < m_minThresholdEffi) continue; // pixel excluded in counting hits
        if (m_channelMask->isActive(slot, channel) and m_asicMask->isActive(slot, channel)) activePixels->Fill(pixel);
      }
 
      if (std::abs(m_selector.getLocalPosition().Z()) > m_excludedZ) {
        auto alphaLow = getObjectPtr<TH1F>(m_alphaLowNames[slot - 1]);
        auto alphaHigh = getObjectPtr<TH1F>(m_alphaHighNames[slot - 1]);
        for (const auto& digit : m_digits) {
          if (digit.getModuleID() != slot) continue;
          if (digit.getHitQuality() != TOPDigit::c_Good) continue;  // junk hit or pixel masked-out
          const auto* pdf = digit.getRelated<TOPAssociatedPDF>();
          if (not pdf) continue;
          const auto* peak = pdf->getSinglePeak();
          if (not peak) continue;  // hit associated with background
          double alpha = acos(std::abs(peak->kzd)) / Unit::deg;
          if (alpha > 60) continue;
          if (alpha < 30) alphaLow->Fill(digit.getPixelID());
          else alphaHigh->Fill(digit.getPixelID());
        }
      }
    }
 
  }

collect() [10/11]

void collect ( )

finalprivatevirtual

Replacement for event().

Fill your calibration data objects here

Reimplemented from CalibrationCollectorModule.

Definition at line 90 of file TOPPulseHeightCollectorModule.cc.

  {
 
    std::vector<TH2F*> slots;
    for (const auto& name : m_names) {
      auto h = getObjectPtr<TH2F>(name);
      slots.push_back(h);
    }
 
    auto h1a = getObjectPtr<TH1F>("time");
    auto h1b = getObjectPtr<TH1F>("time_sel");
    auto h2a = getObjectPtr<TH2F>("ph_vs_width");
    auto h2b = getObjectPtr<TH2F>("ph_vs_width_sel");
 
    for (const auto& digit : m_digits) {
      if (digit.getHitQuality() != TOPDigit::c_Good) continue;
      auto t = digit.getTime();
      auto w = digit.getPulseWidth();
      auto ph = digit.getPulseHeight();
      h1a->Fill(t);
      h2a->Fill(w, ph);
      if (m_widthWindow.size() == 2) {
        if (w < m_widthWindow[0] or w > m_widthWindow[1]) continue;
      }
      if (m_timeWindow.size() == 2) {
        if (t < m_timeWindow[0] or t > m_timeWindow[1]) continue;
      }
      h1b->Fill(t);
      h2b->Fill(w, ph);
      unsigned m = digit.getModuleID() - 1;
      if (m < slots.size()) slots[m]->Fill(digit.getChannel(), ph);
    }
 
  }

collect() [11/11]

void collect ( )

finalprivatevirtual

Replacement for event().

Fill your calibration data objects here

Reimplemented from CalibrationCollectorModule.

Definition at line 233 of file TOPValidationCollectorModule.cc.

defineHisto() [1/3]

void defineHisto ( )

overridevirtual

Histogram definitions such as TH1(), TH2(), TNtuple(), TTree()....

are supposed to be placed in this function.

Reimplemented from HistoModule.

Definition at line 59 of file TOPDQMModule.cc.

defineHisto() [2/3]

void defineHisto ( )

overridevirtual

Histogram definitions such as TH1(), TH2(), TNtuple(), TTree()....

are supposed to be placed in this function.

Reimplemented from HistoModule.

Definition at line 69 of file TOPPDFCheckerModule.cc.

  {
 
    // time vs. pixel ID
    m_hits = new TH2F("hits", "photon hits", 512, 0.5, 512.5,
                      m_numBins, m_minTime, m_maxTime);
    m_hits->SetXTitle("pixel ID");
    m_hits->SetYTitle("time [ns]");
 
    m_pdf = new TH2F("pdf", "PDF", 512, 0.5, 512.5,
                     m_numBins, m_minTime, m_maxTime);
    m_pdf->SetXTitle("pixel ID");
    m_pdf->SetYTitle("time [ns]");
 
    // time vs pixel column
    m_hitsCol = new TH2F("hitsCol", "photon hits", 64, 0.5, 64.5,
                         m_numBins, m_minTime, m_maxTime);
    m_hitsCol->SetXTitle("pixel column");
    m_hitsCol->SetYTitle("time [ns]");
 
    m_pdfCol = new TH2F("pdfCol", "PDF", 64, 0.5, 64.5,
                        m_numBins, m_minTime, m_maxTime);
    m_pdfCol->SetXTitle("pixel column");
    m_pdfCol->SetYTitle("time [ns]");
 
  }

defineHisto() [3/3]

void defineHisto ( )

overridevirtual

Defining the histograms.

Reads once the m_inputDirectoryList to initialize the proper amount of histograms. Every new histogram added to the module has to be initialized here.

Reimplemented from HistoModule.

Definition at line 61 of file TOPTBCComparatorModule.cc.

drawWaveforms()

void drawWaveforms

(

const TOPRawWaveform &

rawwave

)

Draws the full waveforms onto the TProfiles.

Parameters

rawwave

the raw waveform

Definition at line 111 of file TOPWaveformQualityPlotterModule.cc.

DummyFillBranch()

void DummyFillBranch

(

EHistogramType

LoadHisto

)

Fill Dummy for Branch.

Use it when there aren't 2D-Histogram.

Definition at line 432 of file TOPGainEfficiencyCalculatorModule.cc.

  {
    m_fitMax = -1;
    m_hitTiming = -1;
    m_hitTimingSigma = -1;
    m_nEntries = -1;
    m_nCalPulse = -1;
    m_nOverflowEvents = -1;
    m_meanPulseHeight = -1;
    m_meanPulseHeightError = -1;
    m_gain = -1;
    m_efficiency = -1;
    m_p0 = -1;
    m_p1 = -1;
    m_p2 = -1;
    m_x0 = -1;
    m_p0Error = -1;
    m_p1Error = -1;
    m_p2Error = -1;
    m_x0Error = -1;
    m_chisquare = -1;
    m_ndf = -1;
    m_funcFullRangeIntegral = -1;
    m_funcFitRangeIntegral = -1;
    m_histoFitRangeIntegral = -1;
    m_histoMeanAboveThre = -1;
 
    for (auto itr = m_branch[LoadHisto].begin(); itr != m_branch[LoadHisto].end(); ++itr) {
      (*itr)->Fill();
    }
  }

endRun() [1/14]

void endRun ( void  )

overridevirtual

End-of-run action.

Save run-related stuff, such as statistics.

Reimplemented from Module.

Definition at line 400 of file TOPBackgroundModule.cc.

endRun() [2/14]

void endRun ( void  )

overridevirtual

End-of-run action.

Save run-related stuff, such as statistics.

Reimplemented from Module.

Definition at line 219 of file TOPDoublePulseGeneratorModule.cc.

endRun() [3/14]

void endRun ( void  )

overridevirtual

End-of-run action.

Save run-related stuff, such as statistics.

Reimplemented from HistoModule.

Definition at line 256 of file TOPLaserHitSelectorModule.cc.

endRun() [4/14]

void endRun ( void  )

overridevirtual

End-of-run action.

Save run-related stuff, such as statistics.

Reimplemented from HistoModule.

Definition at line 331 of file TOPInterimFENtupleModule.cc.

endRun() [5/14]

void endRun ( void  )

overridevirtual

End-of-run action.

Save run-related stuff, such as statistics.

Reimplemented from Module.

Definition at line 109 of file TOPLaserCalibratorModule.cc.

endRun() [6/14]

void endRun ( void  )

overridevirtual

End-of-run action.

Save run-related stuff, such as statistics.

Reimplemented from Module.

Definition at line 142 of file TOPMCTrackMakerModule.cc.

endRun() [7/14]

void endRun ( void  )

overridevirtual

End-of-run action.

Save run-related stuff, such as statistics.

Reimplemented from Module.

Definition at line 246 of file TOPNtupleModule.cc.

endRun() [8/14]

void endRun ( void  )

overridevirtual

End-of-run action.

Save run-related stuff, such as statistics.

Reimplemented from Module.

Definition at line 123 of file TOPPackerModule.cc.

endRun() [9/14]

void endRun ( void  )

overridevirtual

End-of-run action.

Save run-related stuff, such as statistics.

Reimplemented from Module.

Definition at line 477 of file TOPRawDigitConverterModule.cc.

endRun() [10/14]

void endRun ( void  )

overridevirtual

End-of-run action.

The main analysis loop over the calibration sets happens here. This function has to be modified only if the directory structure of the TOPTimeBaseCalibrator module output is changed Both the core functions analyzeCalFile() and makeComparisons() are called here.

Reimplemented from HistoModule.

Definition at line 356 of file TOPTBCComparatorModule.cc.

endRun() [11/14]

void endRun ( void  )

overridevirtual

End-of-run action.

Save run-related stuff, such as statistics.

Reimplemented from Module.

Definition at line 216 of file TOPTimeBaseCalibratorModule.cc.

endRun() [12/14]

void endRun ( void  )

overridevirtual

End-of-run action.

Save run-related stuff, such as statistics.

Reimplemented from Module.

Definition at line 1245 of file TOPUnpackerModule.cc.

endRun() [13/14]

void endRun ( void  )

overridevirtual

End-of-run action.

Save run-related stuff, such as statistics.

Reimplemented from Module.

Definition at line 140 of file TOPWaveformFeatureExtractorModule.cc.

endRun() [14/14]

void endRun ( void  )

overridevirtual

End-of-run action.

Save run-related stuff, such as statistics.

Reimplemented from HistoModule.

Definition at line 181 of file TOPWaveformQualityPlotterModule.cc.

event()

void event ( void  )

overridevirtual

Event processor.

Convert TOPSimHits of the event to TOPHits.

Reimplemented from Module.

Definition at line 202 of file TOPBackgroundModule.cc.

generatePulseHeight()

double generatePulseHeight ( int  moduleID, int  pixelID  ) const

private

Generates and returns pulse height.

Parameters

moduleID	module ID (1-based)
pixelID	pixel ID (1-based)

Returns

pulse height [ADC counts]

Definition at line 533 of file TOPDigitizerModule.cc.

  {
    if (m_useDatabase) {
      const auto& channelMapper = TOPGeometryPar::Instance()->getChannelMapper();
      auto channel = channelMapper.getChannel(pixelID);
      if (m_pulseHeights->isCalibrated(moduleID, channel)) {
        const auto& par = m_pulseHeights->getParameters(moduleID, channel);
        PulseHeightGenerator generator(par.x0, par.p1, par.p2, m_ADCmax);
        return generator.generate();
      }
    }
 
    return m_pulseHeightGenerator.generate();
  }

getDirectionCustom()

XYZVector getDirectionCustom ( ) const

private

Return photon direction according to a custom angular distribution given by TFormula.

Returns

photon direction (unit vector)

Definition at line 269 of file OpticalGunModule.cc.

  {
    double alpha = m_customDistribution->GetRandom();
    double phi = 2.0 * M_PI * gRandom->Rndm();
    return XYZVector(cos(phi) * sin(alpha), sin(phi) * sin(alpha), cos(alpha));
  }

getDirectionGaussian()

XYZVector getDirectionGaussian ( ) const

private

Return photon direction according to a projected 2D gaussian distribution based on numerical aperture NA.

Returns

photon direction (unit vector)

Definition at line 237 of file OpticalGunModule.cc.

getDirectionLambertian()

XYZVector getDirectionLambertian ( ) const

private

Return photon direction according to a lambertian distribution with opening angle alpha.

Returns

photon direction (unit vector)

Definition at line 260 of file OpticalGunModule.cc.

getDirectionUniform()

XYZVector getDirectionUniform ( ) const

private

Return photon direction according to a projected uniform distribution with opening angle alpha.

Be careful.

Returns

photon direction (unit vector)

Definition at line 252 of file OpticalGunModule.cc.

getFrontEndName()

std::string getFrontEndName ( RawTOP &  raw, int  finesse  ) const

private

Returns the name of the front-end.

Parameters

raw	raw data
finesse	finesse number

Returns

front-end name

Definition at line 219 of file TOPUnpackerModule.cc.

getModuleID()

int getModuleID ( const Track &  track ) const

private

Returns slot ID of the module that is hit by the track.

Parameters

track

charged track

Returns

slotID or 0 if track does not hit any module

Definition at line 515 of file TOPDQMModule.cc.

getMostProbable()

Const::ChargedStable getMostProbable ( const Track &  track )

private

Returns most probable charged stable particle according to dEdx and predefined prior probabilities.

Parameters

track

reconstructed track

Returns

charged stable

Definition at line 564 of file TOPBunchFinderModule.cc.

  {
 
    std::vector<double> logL;
    std::vector<double> priors;
 
    if (m_usePIDLikelihoods) {
      const auto* pid = track.getRelated<PIDLikelihood>();
      if (not pid) {
        m_nodEdxCount++;
        return Const::pion;
      }
      auto subset = Const::PIDDetectorSet(Const::SVD);
      subset += Const::PIDDetectorSet(Const::CDC);
      for (const auto& type : Const::chargedStableSet) {
        logL.push_back(pid->getLogL(type, subset));
        priors.push_back(m_priors[abs(type.getPDGCode())]);
      }
    } else {
      const auto* cdcdedx = track.getRelated<CDCDedxLikelihood>();
      const auto* vxddedx = track.getRelated<VXDDedxLikelihood>();
      if (not cdcdedx and not vxddedx) {
        m_nodEdxCount++;
        return Const::pion;
      }
      for (const auto& type : Const::chargedStableSet) {
        if (cdcdedx and vxddedx) {
          logL.push_back(cdcdedx->getLogL(type) + vxddedx->getLogL(type));
        } else if (cdcdedx) {
          logL.push_back(cdcdedx->getLogL(type));
        } else {
          logL.push_back(vxddedx->getLogL(type));
        }
        priors.push_back(m_priors[abs(type.getPDGCode())]);
      }
    }
 
    // get maximal logL
    auto logL_max = logL[0];
    for (auto x : logL) {
      if (x > logL_max) logL_max = x;
    }
 
    // calculate probabilities, normalizaton is not needed
    std::vector<double> probability(logL.size());
    for (unsigned i = 0; i < logL.size(); ++i) {
      probability[i] = exp(logL[i] - logL_max) * priors[i];
    }
 
    // find most probable
    unsigned i0 = 0;
    for (unsigned i = 0; i < probability.size(); ++i) {
      if (probability[i] > probability[i0]) i0 = i;
    }
    return Const::chargedStableSet.at(i0);
 
  }

getReferenceTiming()

void getReferenceTiming ( )

private

Find reference timing.

In case that the waveform analysis is enabled, try to find a pair of calibration pulses and timing of the first one is used as reference timing for the correcponding asic. If not enabled, feature extracted timing for a calibration channel is used.

Definition at line 339 of file TOPInterimFENtupleModule.cc.

getTimeOffset()

TOPDigitizerModule::TimeOffset getTimeOffset ( double  trgOffset, int  moduleID, int  pixelID  )

private

Returns a complete time offset by adding time mis-calibration to trgOffset.

Parameters

trgOffset	trigger related time offset
moduleID	slot ID
pixelID	pixel ID

Returns

time offset and its error squared

Definition at line 496 of file TOPDigitizerModule.cc.

getTimeSeed()

TOPBunchFinderModule::TimeSeed getTimeSeed ( )

private

Returns a time seed.

Returns

time seed

Definition at line 643 of file TOPBunchFinderModule.cc.

initialize() [1/29]

void initialize ( void  )

overridevirtual

Initialize the Module.

This method is called at the beginning of data processing.

Reimplemented from Module.

Definition at line 90 of file OpticalGunModule.cc.

initialize() [2/29]

void initialize ( void  )

overridevirtual

Initialize the Module.

This method is called at the beginning of data processing.

Reimplemented from Module.

Definition at line 87 of file TOPAlignerModule.cc.

  {
    // check if target module ID is valid
 
    const auto* geo = TOPGeometryPar::Instance()->getGeometry();
    if (not geo->isModuleIDValid(m_targetMid)) {
      B2ERROR("Target module ID = " << m_targetMid << " is invalid.");
    }
 
    // check if sample type is valid
 
    if (not(m_sample == "dimuon" or m_sample == "bhabha" or m_sample == "cosmics")) {
      B2ERROR("Invalid sample type '" << m_sample << "'");
    }
    if (m_sample == "bhabha") m_chargedStable = Const::electron;
 
    // set track selector
 
    m_selector = TrackSelector(m_sample);
    m_selector.setMinMomentum(m_minMomentum);
    m_selector.setDeltaEcms(m_deltaEcms);
    m_selector.setCutOnPOCA(m_dr, m_dz);
    m_selector.setCutOnLocalZ(m_minZ, m_maxZ);
 
    // set alignment object
 
    m_align.setModuleID(m_targetMid);
    m_align.setSteps(m_stepPosition, m_stepAngle, m_stepTime);
    m_align.setParameters(m_parInit);
    for (const auto& parName : m_parFixed) {
      m_align.fixParameter(parName);
    }
 
    // input
 
    m_digits.isRequired();
    m_tracks.isRequired();
    m_extHits.isRequired();
    m_recBunch.isOptional();
 
    // open output file
 
    m_file = TFile::Open(m_outFileName.c_str(), "RECREATE");
    if (m_file->IsZombie()) {
      B2FATAL("Couldn't open file '" << m_outFileName << "' for writing!");
      return;
    }
 
    // create output tree
 
    m_alignTree = new TTree("alignTree", "TOP alignment results");
    m_alignTree->Branch("ModuleId", &m_targetMid);
    m_alignTree->Branch("iter", &m_iter);
    m_alignTree->Branch("ntrk", &m_ntrk);
    m_alignTree->Branch("errorCode", &m_errorCode);
    m_alignTree->Branch("iterPars", &m_vAlignPars);
    m_alignTree->Branch("iterParsErr", &m_vAlignParsErr);
    m_alignTree->Branch("valid", &m_valid);
    m_alignTree->Branch("numPhot", &m_numPhot);
    m_alignTree->Branch("x", &m_x);
    m_alignTree->Branch("y", &m_y);
    m_alignTree->Branch("z", &m_z);
    m_alignTree->Branch("p", &m_p);
    m_alignTree->Branch("theta", &m_theta);
    m_alignTree->Branch("phi", &m_phi);
    m_alignTree->Branch("r_poca", &m_pocaR);
    m_alignTree->Branch("z_poca", &m_pocaZ);
    m_alignTree->Branch("x_poca", &m_pocaX);
    m_alignTree->Branch("y_poca", &m_pocaY);
    m_alignTree->Branch("Ecms", &m_cmsE);
    m_alignTree->Branch("charge", &m_charge);
    m_alignTree->Branch("PDG", &m_PDG);
 
  }

initialize() [3/29]

void initialize ( void  )

overridevirtual

Initialize the Module.

This method is called at the beginning of data processing.

Reimplemented from Module.

Definition at line 116 of file TOPBackgroundModule.cc.

initialize() [4/29]

void initialize ( void  )

overridevirtual

Initialize the Module.

This method is called at the beginning of data processing.

Reimplemented from Module.

Definition at line 111 of file TOPBunchFinderModule.cc.

initialize() [5/29]

void initialize ( void  )

overridevirtual

Initialize the Module.

This method is called at the beginning of data processing.

Reimplemented from Module.

Definition at line 76 of file TOPChannelT0CalibratorModule.cc.

  {
    // input collections
    m_digits.isRequired();
    m_tracks.isRequired();
    m_extHits.isRequired();
    m_recBunch.isOptional();
 
    // Parse PDF option
    if (m_pdfOption == "rough") {
      m_PDFOption = PDFConstructor::c_Rough;
    } else if (m_pdfOption == "fine") {
      m_PDFOption = PDFConstructor::c_Fine;
    } else if (m_pdfOption == "optimal") {
      m_PDFOption = PDFConstructor::c_Optimal;
    } else {
      B2ERROR("TOPPDFDebuggerModule: unknown PDF option '" << m_pdfOption << "'");
    }
 
    // set track selector
    m_selector = TrackSelector(m_sample);
    m_selector.setMinMomentum(m_minMomentum);
    m_selector.setDeltaEcms(m_deltaEcms);
    m_selector.setCutOnPOCA(m_dr, m_dz);
    m_selector.setCutOnLocalZ(m_minZ, m_maxZ);
 
    // minimum finders
    double tmin = -m_timeRange / 2;
    double tmax =  m_timeRange / 2;
    for (unsigned i = 0; i < 2; i++) {
      for (unsigned m = 0; m < c_numModules; m++) {
        for (unsigned c = 0; c < c_numChannels; c++) {
          m_finders[i][m][c] = Chi2MinimumFinder1D(m_numBins, tmin, tmax);
        }
      }
    }
 
    // if file name includes *'s replace them with a run number
    auto pos = m_outFileName.find("*");
    if (pos != std::string::npos) {
      StoreObjPtr<EventMetaData> evtMetaData;
      auto run = std::to_string(evtMetaData->getRun());
      while (run.size() < 5) run = "0" + run;
      while (pos != std::string::npos) {
        m_outFileName.replace(pos, 1, run);
        pos = m_outFileName.find("*");
      }
    }
 
    // open root file for ntuple and histogram output
    m_file = TFile::Open(m_outFileName.c_str(), "RECREATE");
    if (not m_file) {
      B2ERROR("Cannot open output file '" << m_outFileName << "'");
      return;
    }
 
    // histograms
    for (unsigned module = 0; module < c_numModules; module++) {
      int moduleID = module + 1;
 
      std::string slotNum = std::to_string(moduleID);
      if (moduleID < 10) slotNum = "0" + slotNum;
 
      std::string name = "numHits_slot" + slotNum;
      std::string title = "Number of hits per channel for slot " + slotNum;
      TH1F h1(name.c_str(), title.c_str(), c_numChannels, 0, c_numChannels);
      h1.SetXTitle("channel number");
      h1.SetYTitle("hits per channel");
      m_hits1D.push_back(h1);
 
      name = "timeHits_slot" + slotNum;
      title = "hit time vs. channel for slot " + slotNum;
      TH2F h2(name.c_str(), title.c_str(), c_numChannels, 0, c_numChannels,
              200, 0.0, 20.0);
      h2.SetXTitle("channel number");
      h2.SetYTitle("time [ns]");
      m_hits2D.push_back(h2);
    }
 
    // create output tree
 
    m_tree = new TTree("tree", "Channel T0 calibration results");
    m_tree->Branch("slot", &m_moduleID);
    m_tree->Branch("numPhotons", &m_numPhotons);
    m_tree->Branch("x", &m_x);
    m_tree->Branch("y", &m_y);
    m_tree->Branch("z", &m_z);
    m_tree->Branch("p", &m_p);
    m_tree->Branch("theta", &m_theta);
    m_tree->Branch("phi", &m_phi);
    m_tree->Branch("r_poca", &m_pocaR);
    m_tree->Branch("z_poca", &m_pocaZ);
    m_tree->Branch("x_poca", &m_pocaX);
    m_tree->Branch("y_poca", &m_pocaY);
    m_tree->Branch("Ecms", &m_cmsE);
    m_tree->Branch("charge", &m_charge);
    m_tree->Branch("PDG", &m_PDG);
 
  }

initialize() [6/29]

void initialize ( void  )

overridevirtual

Initialize the Module.

This method is called at the beginning of data processing.

Reimplemented from Module.

Definition at line 56 of file TOPChannelT0MCModule.cc.

  {
 
    m_digits.isRequired();
 
  }

initialize() [7/29]

void initialize ( void  )

overridevirtual

Initialize the Module.

This method is called at the beginning of data processing.

Reimplemented from Module.

Definition at line 73 of file TOPCommonT0CalibratorModule.cc.

  {
    // input collections
    m_digits.isRequired();
    m_tracks.isRequired();
    m_extHits.isRequired();
    m_recBunch.isOptional();
 
    // bunch separation in time
 
    const auto* geo = TOPGeometryPar::Instance()->getGeometry();
    m_bunchTimeSep = geo->getNominalTDC().getSyncTimeBase() / 24;
 
    // Parse PDF option
    if (m_pdfOption == "rough") {
      m_PDFOption = PDFConstructor::c_Rough;
    } else if (m_pdfOption == "fine") {
      m_PDFOption = PDFConstructor::c_Fine;
    } else if (m_pdfOption == "optimal") {
      m_PDFOption = PDFConstructor::c_Optimal;
    } else {
      B2ERROR("TOPPDFDebuggerModule: unknown PDF option '" << m_pdfOption << "'");
    }
 
    // set track selector
    m_selector = TrackSelector(m_sample);
    m_selector.setMinMomentum(m_minMomentum);
    m_selector.setDeltaEcms(m_deltaEcms);
    m_selector.setCutOnPOCA(m_dr, m_dz);
    m_selector.setCutOnLocalZ(m_minZ, m_maxZ);
 
    // Chi2 minimum finders
    double tmin = -m_timeRange / 2;
    double tmax =  m_timeRange / 2;
    for (unsigned i = 0; i < c_numSets; i++) {
      m_finders[i] = Chi2MinimumFinder1D(m_numBins, tmin, tmax);
    }
 
    // if file name includes *'s replace them with a run number
    auto pos = m_outFileName.find("*");
    if (pos != std::string::npos) {
      StoreObjPtr<EventMetaData> evtMetaData;
      auto run = std::to_string(evtMetaData->getRun());
      while (run.size() < 5) run = "0" + run;
      while (pos != std::string::npos) {
        m_outFileName.replace(pos, 1, run);
        pos = m_outFileName.find("*");
      }
    }
 
    // open root file for ntuple and histogram output
    m_file = TFile::Open(m_outFileName.c_str(), "RECREATE");
    if (not m_file) {
      B2ERROR("Cannot open output file '" << m_outFileName << "'");
      return;
    }
 
    // control histograms
    m_hits1D = TH1F("numHits", "Number of photons per slot",
                    c_numModules, 0.5, static_cast<float>(c_numModules) + 0.5);
    m_hits1D.SetXTitle("slot number");
    m_hits1D.SetYTitle("hits per slot");
 
    m_hits2D = TH2F("timeHits", "Photon times vs. boardstacks",
                    c_numModules * 4, 0.5, static_cast<float>(c_numModules) + 0.5, 200, 0.0, 20.0);
    m_hits2D.SetXTitle("slot number");
    m_hits2D.SetYTitle("time [ns]");
 
    // create output tree
    m_tree = new TTree("tree", "Channel T0 calibration results");
    m_tree->Branch("slot", &m_moduleID);
    m_tree->Branch("numPhotons", &m_numPhotons);
    m_tree->Branch("x", &m_x);
    m_tree->Branch("y", &m_y);
    m_tree->Branch("z", &m_z);
    m_tree->Branch("p", &m_p);
    m_tree->Branch("theta", &m_theta);
    m_tree->Branch("phi", &m_phi);
    m_tree->Branch("r_poca", &m_pocaR);
    m_tree->Branch("z_poca", &m_pocaZ);
    m_tree->Branch("x_poca", &m_pocaX);
    m_tree->Branch("y_poca", &m_pocaY);
    m_tree->Branch("Ecms", &m_cmsE);
    m_tree->Branch("charge", &m_charge);
    m_tree->Branch("PDG", &m_PDG);
 
  }

initialize() [8/29]

void initialize ( void  )

overridevirtual

Initialize the Module.

This method is called at the beginning of data processing.

Reimplemented from Module.

Definition at line 78 of file TOPCosmicT0FinderModule.cc.

  {
    // input
 
    StoreArray<TOPDigit> topDigits;
    topDigits.isRequired();
 
    StoreArray<Track> tracks;
    tracks.isRequired();
 
    StoreArray<ExtHit> extHits;
    extHits.isRequired();
 
    StoreArray<TOPBarHit> barHits;
    barHits.isOptional();
 
    // output
 
    StoreArray<TOPTimeZero> timeZeros;
    timeZeros.registerInDataStore();
    timeZeros.registerRelationTo(tracks);
    timeZeros.registerRelationTo(extHits);
    timeZeros.registerRelationTo(barHits);
  }

initialize() [9/29]

void initialize ( void  )

overridevirtual

Initialize the Module.

This method is called at the beginning of data processing.

Reimplemented from Module.

Definition at line 103 of file TOPDigitizerModule.cc.

initialize() [10/29]

void initialize ( void  )

overridevirtual

Initialize the Module.

This method is called at the beginning of data processing.

Reimplemented from Module.

Definition at line 68 of file TOPTriggerDigitizerModule.cc.

  {
    // input
    StoreArray<TOPRawWaveform> waveforms;
    waveforms.isRequired();
 
    // output
    StoreArray<TOPTriggerDigit> digits;
    digits.registerInDataStore();
    digits.registerRelationTo(waveforms);
    StoreObjPtr<TOPTriggerMCInfo> mcInfo;
    mcInfo.registerInDataStore();
 
    if (m_samplingPhase < 0 or m_samplingPhase >= c_SamplingCycle)
      B2ERROR("samplingPhase must be positive and less than " << c_SamplingCycle);
 
  }

initialize() [11/29]

void initialize ( void  )

overridevirtual

Initialize the Module.

This method is called at the beginning of data processing.

Reimplemented from Module.

Definition at line 68 of file TOPCalPulseGeneratorModule.cc.

  {
    // Output
 
    m_calPulses.registerInDataStore();
 
    // prepare vectors to loop on
 
    const auto* geo = TOPGeometryPar::Instance()->getGeometry();
 
    if (m_moduleIDs.empty()) {
      for (const auto& module : geo->getModules()) {
        m_moduleIDs.push_back(module.getModuleID());
      }
    } else {
      for (auto moduleID : m_moduleIDs) {
        if (!geo->isModuleIDValid(moduleID))
          B2ERROR("Invalid module ID found in input list: " << moduleID);
      }
    }
 
    if (m_asicChannels.empty()) {
      for (unsigned ch = 0; ch < 8; ch++) m_asicChannels.push_back(ch);
    } else {
      for (unsigned ch : m_asicChannels) {
        if (ch > 7)
          B2ERROR("Invalid ASIC channel found in input list: " << ch);
      }
    }
 
  }

initialize() [12/29]

void initialize ( void  )

overridevirtual

Initialize the Module.

This method is called at the beginning of data processing.

Reimplemented from Module.

Definition at line 90 of file TOPDoublePulseGeneratorModule.cc.

initialize() [13/29]

void initialize ( void  )

overridevirtual

Initialize the Module.

This method is called at the beginning of data processing.

better use isRequired(), but RawFTSW is not in sim

Reimplemented from HistoModule.

Definition at line 291 of file TOPDQMModule.cc.

initialize() [14/29]

void initialize ( void  )

overridevirtual

Initialize the Module.

This method is called at the beginning of data processing.

Reimplemented from HistoModule.

Definition at line 94 of file TOPLaserHitSelectorModule.cc.

initialize() [15/29]

void initialize ( void  )

overridevirtual

Initialize the Module.

This method is called at the beginning of data processing.

Reimplemented from Module.

Definition at line 53 of file TOPGeometryParInitializerModule.cc.

  {
    auto* gp = TOPGeometryPar::Instance();
    if (gp->isValid()) {
      B2WARNING("TOPGeometryPar has already been initialized");
      return;
    }
 
    if (m_useDB) {
      gp->Initialize();
    } else {
      bool found = false;
      GearDir detectorDir("/Detector");
      for (const GearDir& component : detectorDir.getNodes("DetectorComponent")) {
        if (component.getString("@name") == "TOP") {
          gp->Initialize(GearDir(component, "Content"));
          found = true;
          break;
        }
      }
      if (not found) B2ERROR("Component TOP does not exist or is empty");
    }
 
  }

initialize() [16/29]

void initialize ( void  )

overridevirtual

Initialize the Module.

This method is called at the beginning of data processing.

Reimplemented from HistoModule.

Definition at line 80 of file TOPInterimFENtupleModule.cc.

initialize() [17/29]

void initialize ( void  )

overridevirtual

Initialize the Module.

This method is called at the beginning of data processing.

Reimplemented from Module.

Definition at line 75 of file TOPLaserCalibratorModule.cc.

initialize() [18/29]

void initialize ( void  )

overridevirtual

Initialize the Module.

This method is called at the beginning of data processing.

Reimplemented from Module.

Definition at line 61 of file TOPMCTrackMakerModule.cc.

initialize() [19/29]

void initialize ( void  )

overridevirtual

Initialize the Module.

This method is called at the beginning of data processing.

Reimplemented from Module.

Definition at line 72 of file TOPModuleT0CalibratorModule.cc.

initialize() [20/29]

void initialize ( void  )

overridevirtual

Initialize the Module.

This method is called at the beginning of data processing.

Reimplemented from Module.

Definition at line 70 of file TOPNtupleModule.cc.

initialize() [21/29]

void initialize ( void  )

overridevirtual

Initialize the Module.

This method is called at the beginning of data processing.

Reimplemented from Module.

Definition at line 70 of file TOPPackerModule.cc.

initialize() [22/29]

void initialize ( void  )

overridevirtual

Initialize the Module.

This method is called at the beginning of data processing.

Reimplemented from HistoModule.

Definition at line 96 of file TOPPDFCheckerModule.cc.

initialize() [23/29]

void initialize ( void  )

overridevirtual

Initialize the Module.

This method is called at the beginning of data processing.

Reimplemented from Module.

Definition at line 72 of file TOPPDFDebuggerModule.cc.

initialize() [24/29]

void initialize ( void  )

overridevirtual

Initialize the Module.

This method is called at the beginning of data processing.

Reimplemented from Module.

Definition at line 114 of file TOPRawDigitConverterModule.cc.

initialize() [25/29]

void initialize ( void  )

overridevirtual

Initialize the Module.

This method is called at the beginning of data processing.

Reimplemented from Module.

Definition at line 69 of file TOPReconstructorModule.cc.

  {
    // input
 
    m_digits.isRequired(m_topDigitCollectionName);
    m_tracks.isRequired();
    m_extHits.isRequired();
    m_barHits.isOptional();
    m_recBunch.isOptional();
 
    // output
 
    m_likelihoods.registerInDataStore(m_topLikelihoodCollectionName);
    m_likelihoods.registerRelationTo(m_extHits);
    m_likelihoods.registerRelationTo(m_barHits);
    m_tracks.registerRelationTo(m_likelihoods);
 
    m_topPulls.registerInDataStore(m_topPullCollectionName, DataStore::c_DontWriteOut);
    m_tracks.registerRelationTo(m_topPulls, DataStore::c_Event, DataStore::c_DontWriteOut);
  }

initialize() [26/29]

void initialize ( void  )

overridevirtual

Initialize the Module.

This method is called at the beginning of data processing.

Reimplemented from Module.

Definition at line 106 of file TOPTimeBaseCalibratorModule.cc.

initialize() [27/29]

void initialize ( void  )

overridevirtual

Initialize the Module.

This method is called at the beginning of data processing.

Reimplemented from Module.

Definition at line 71 of file TOPTimeRecalibratorModule.cc.

initialize() [28/29]

void initialize ( void  )

overridevirtual

Initialize the Module.

This method is called at the beginning of data processing.

Reimplemented from Module.

Definition at line 82 of file TOPUnpackerModule.cc.

initialize() [29/29]

void initialize ( void  )

overridevirtual

Initialize the Module.

This method is called at the beginning of data processing.

Reimplemented from Module.

Definition at line 70 of file TOPWaveformFeatureExtractorModule.cc.

isBucketFilled()

bool isBucketFilled ( int  bunchNo )

private

Does reconstructed bunch number correspond to filled bucket.

Parameters

bunchNo

reconstructed relative bunch number

Returns

true if filled

Definition at line 670 of file TOPBunchFinderModule.cc.

isFromThisParticle()

bool isFromThisParticle ( const TOPDigit &  digit, const MCParticle *  particle  )

private

Checks if digit comes from given MC particle.

Parameters

digit	TOP digit
particle	MC particle

Definition at line 207 of file TOPPDFCheckerModule.cc.

isInsideSlit()

bool isInsideSlit ( const ROOT::Math::XYZPoint &  point, const ROOT::Math::XYZVector &  direction  ) const

private

Checks if photon passes the slit.

Parameters

point	photon emission point
direction	photon emission direction

Returns

true on success

Definition at line 218 of file OpticalGunModule.cc.

isRunningOffsetSubtracted() [1/2]

bool isRunningOffsetSubtracted ( )

private

Checks if running offset is subtracted in TOPDigits.

Returns

true if subtracted at least in one digit

Definition at line 305 of file TOPCommonT0CalibratorModule.cc.

isRunningOffsetSubtracted() [2/2]

bool isRunningOffsetSubtracted ( )

private

Checks if running offset is subtracted in TOPDigits.

Returns

true if subtracted at least in one digit

Definition at line 333 of file TOPModuleT0CalibratorModule.cc.

Iteration()

void Iteration ( const std::vector< TwoTimes > &  ntuple, std::vector< double > &  xval  )

private

Iteration function called by iterativeTBC()

Parameters

ntuple	ntuple data
xval	TBC constants of 256 samples, time interval is the difference of nearby xvals, xval[0]=0 and xval[256]=2*m_syncTimeBase

Definition at line 623 of file TOPTimeBaseCalibratorModule.cc.

iterativeTBC()

bool iterativeTBC ( const std::vector< TwoTimes > &  ntuple, unsigned  scrodID, unsigned  scrodChannel, double  meanTimeDifference, TH1F &  Hchi2, TH1F &  Hndf, TH1F &  HDeltaT  )

private

Method by iteration of chi2 minimization.

Parameters

ntuple	ntuple data
scrodID	SCROD ID
scrodChannel	channel number within SCROD (0 - 127)
meanTimeDifference	average time difference [samples]
Hchi2	histogram to store normalized chi^2
Hndf	histogram to store degrees of freedom
HDeltaT	histogram to store fittet double pulse delay

Returns

true on success

Definition at line 538 of file TOPTimeBaseCalibratorModule.cc.

LoadHistograms()

void LoadHistograms

(

const std::string &

histotype

)

Load 2D histograms from a given input file (output of TOPLaserHitSelector) and create timing and charge distribution as projection histograms for the x- and y-axis, respectively.

Timing cut is also applied for charge distributiion

Definition at line 218 of file TOPGainEfficiencyCalculatorModule.cc.

matrixInversion()

bool matrixInversion ( const std::vector< TwoTimes > &  ntuple, unsigned  scrodID, unsigned  scrodChannel, double  meanTimeDifference, TH1F &  Hchi2, TH1F &  Hndf, TH1F &  HDeltaT  )

private

Method by matrix inversion.

Parameters

ntuple	ntuple data
scrodID	SCROD ID
scrodChannel	channel number within SCROD (0 - 127)
meanTimeDifference	average time difference [samples]
Hchi2	histogram to store normalized chi^2
Hndf	histogram to store degrees of freedom
HDeltaT	histogram to store fittet double pulse delay

Returns

true on success

Definition at line 398 of file TOPTimeBaseCalibratorModule.cc.

prepare() [1/11]

void prepare ( )

finalprivatevirtual

Replacement for initialize().

Register calibration dataobjects here as well

Reimplemented from CalibrationCollectorModule.

Definition at line 79 of file TOPAlignmentCollectorModule.cc.

prepare() [2/11]

void prepare ( )

finalprivatevirtual

Replacement for initialize().

Register calibration dataobjects here as well.

Reimplemented from CalibrationCollectorModule.

Definition at line 52 of file TOPAsicShiftsBS13dCollectorModule.cc.

prepare() [3/11]

void prepare ( )

finalprivatevirtual

Replacement for initialize().

Register calibration dataobjects here as well

Reimplemented from CalibrationCollectorModule.

Definition at line 39 of file TOPChannelMaskCollectorModule.cc.

prepare() [4/11]

void prepare ( )

finalprivatevirtual

Replacement for initialize().

Register calibration dataobjects here as well

Reimplemented from CalibrationCollectorModule.

Definition at line 49 of file TOPCommonT0BFCollectorModule.cc.

prepare() [5/11]

void prepare ( )

finalprivatevirtual

Replacement for initialize().

Register calibration dataobjects here as well

Reimplemented from CalibrationCollectorModule.

Definition at line 75 of file TOPCommonT0LLCollectorModule.cc.

prepare() [6/11]

void prepare ( )

finalprivatevirtual

Replacement for initialize().

Register calibration dataobjects here as well

Reimplemented from CalibrationCollectorModule.

Definition at line 50 of file TOPModuleT0DeltaTCollectorModule.cc.

prepare() [7/11]

void prepare ( )

finalprivatevirtual

Replacement for initialize().

Register calibration dataobjects here as well

Reimplemented from CalibrationCollectorModule.

Definition at line 75 of file TOPModuleT0LLCollectorModule.cc.

prepare() [8/11]

void prepare ( )

finalprivatevirtual

Replacement for initialize().

Register calibration dataobjects here as well

Reimplemented from CalibrationCollectorModule.

Definition at line 41 of file TOPOffsetCollectorModule.cc.

prepare() [9/11]

void prepare ( )

finalprivatevirtual

Replacement for initialize().

Register calibration dataobjects here as well

Reimplemented from CalibrationCollectorModule.

Definition at line 58 of file TOPPhotonYieldsCollectorModule.cc.

prepare() [10/11]

void prepare ( )

finalprivatevirtual

Replacement for initialize().

Register calibration dataobjects here as well

Reimplemented from CalibrationCollectorModule.

Definition at line 52 of file TOPPulseHeightCollectorModule.cc.

prepare() [11/11]

void prepare ( )

finalprivatevirtual

Replacement for initialize().

Register calibration dataobjects here as well

Reimplemented from CalibrationCollectorModule.

Definition at line 76 of file TOPValidationCollectorModule.cc.

printTheError()

bool printTheError ( )

private

Error messages suppression logic.

Returns

true to print the error message, false to suppress it

Definition at line 234 of file TOPUnpackerModule.cc.

saveAsHistogram() [1/4]

void saveAsHistogram ( const std::vector< double > &  vec, const std::string &  name, const std::string &  title, const std::string &  xTitle = "", const std::string &  yTitle = ""  ) const

private

Save vector to histogram and write it out.

Parameters

vec	vector of bin values
name	histogram name
title	histogram title
xTitle	x-axis title
yTitle	y-axis title

Definition at line 299 of file TOPDoublePulseGeneratorModule.cc.

saveAsHistogram() [2/4]

void saveAsHistogram ( const std::vector< double > &  vec, const std::string &  name, const std::string &  title, const std::string &  xTitle = "", const std::string &  yTitle = ""  ) const

private

Save vector to histogram and write it out.

Parameters

vec	vector of bin values
name	histogram name
title	histogram title
xTitle	x-axis title
yTitle	y-axis title

Definition at line 711 of file TOPTimeBaseCalibratorModule.cc.

saveAsHistogram() [3/4]

void saveAsHistogram ( const std::vector< double > &  vec, const std::vector< double > &  err, const std::string &  name, const std::string &  title, const std::string &  xTitle = "", const std::string &  yTitle = ""  ) const

private

Save vector and errors to histogram and write it out.

Parameters

vec	vector of bin values
err	vector of bin errors
name	histogram name
title	histogram title
xTitle	x-axis title
yTitle	y-axis title

Definition at line 730 of file TOPTimeBaseCalibratorModule.cc.

saveAsHistogram() [4/4]

void saveAsHistogram ( const TMatrixDSym &  M, const std::string &  name, const std::string &  title  ) const

private

Save matrix to histogram and write it out.

Parameters

M	matrix
name	histogram name
title	histogram title

Definition at line 750 of file TOPTimeBaseCalibratorModule.cc.

setFinder()

int setFinder ( TOP::Chi2MinimumFinder1D &  finder, const TOP::PDFConstructor &  reco, double  timeMin, double  timeMax  )

private

Sets finder object with chi2 values.

Parameters

finder	finder object
reco	reconstruction object
timeMin	lower edge of photon time window
timeMax	upper edge of photon time window

Returns

number of photons in the time window

Definition at line 623 of file TOPBunchFinderModule.cc.

startRun()

void startRun ( )

finalprivatevirtual

Replacement for beginRun().

Do anything you would normally do in beginRun here

Reimplemented from CalibrationCollectorModule.

Definition at line 167 of file TOPValidationCollectorModule.cc.

terminate() [1/23]

void terminate ( void  )

overridevirtual

Termination action.

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

Reimplemented from Module.

Definition at line 252 of file TOPAlignerModule.cc.

terminate() [2/23]

void terminate ( void  )

overridevirtual

Termination action.

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

Reimplemented from Module.

Definition at line 405 of file TOPBackgroundModule.cc.

terminate() [3/23]

void terminate ( void  )

overridevirtual

Termination action.

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

Reimplemented from Module.

Definition at line 557 of file TOPBunchFinderModule.cc.

terminate() [4/23]

void terminate ( void  )

overridevirtual

Termination action.

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

Reimplemented from Module.

Definition at line 259 of file TOPChannelT0CalibratorModule.cc.

terminate() [5/23]

void terminate ( void  )

overridevirtual

Termination action.

Output MC extraction

Reimplemented from Module.

Definition at line 93 of file TOPChannelT0MCModule.cc.

terminate() [6/23]

void terminate ( void  )

overridevirtual

Termination action.

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

Reimplemented from Module.

Definition at line 243 of file TOPCommonT0CalibratorModule.cc.

terminate() [7/23]

void terminate ( void  )

overridevirtual

Termination action.

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

Reimplemented from Module.

Definition at line 249 of file TOPCosmicT0FinderModule.cc.

terminate() [8/23]

void terminate ( void  )

overridevirtual

Termination action.

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

Reimplemented from Module.

Definition at line 154 of file TOPTriggerDigitizerModule.cc.

terminate() [9/23]

void terminate ( void  )

overridevirtual

Termination action.

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

Reimplemented from Module.

Definition at line 223 of file TOPDoublePulseGeneratorModule.cc.

terminate() [10/23]

void terminate ( void  )

overridevirtual

Termination action.

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

Reimplemented from HistoModule.

Definition at line 173 of file TOPGainEfficiencyCalculatorModule.cc.

terminate() [11/23]

void terminate ( void  )

overridevirtual

Termination action.

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

Reimplemented from HistoModule.

Definition at line 260 of file TOPLaserHitSelectorModule.cc.

terminate() [12/23]

void terminate ( void  )

overridevirtual

Termination action.

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

Reimplemented from HistoModule.

Definition at line 335 of file TOPInterimFENtupleModule.cc.

terminate() [13/23]

void terminate ( void  )

overridevirtual

Termination action.

Do fits , clean-up, close files, summarize statistics, etc.

Reimplemented from Module.

Definition at line 113 of file TOPLaserCalibratorModule.cc.

terminate() [14/23]

void terminate ( void  )

overridevirtual

Termination action.

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

Reimplemented from Module.

Definition at line 146 of file TOPMCTrackMakerModule.cc.

terminate() [15/23]

void terminate ( void  )

overridevirtual

Termination action.

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

Reimplemented from Module.

Definition at line 252 of file TOPModuleT0CalibratorModule.cc.

terminate() [16/23]

void terminate ( void  )

overridevirtual

Termination action.

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

Reimplemented from Module.

Definition at line 250 of file TOPNtupleModule.cc.

terminate() [17/23]

void terminate ( void  )

overridevirtual

Termination action.

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

Reimplemented from Module.

Definition at line 127 of file TOPPackerModule.cc.

terminate() [18/23]

void terminate ( void  )

overridevirtual

Termination action.

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

Reimplemented from HistoModule.

Definition at line 182 of file TOPPDFCheckerModule.cc.

terminate() [19/23]

void terminate ( void  )

overridevirtual

Termination action.

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

Reimplemented from Module.

Definition at line 481 of file TOPRawDigitConverterModule.cc.

terminate() [20/23]

void terminate ( void  )

overridevirtual

Termination action.

The historam writing takes place here

Reimplemented from HistoModule.

Definition at line 453 of file TOPTBCComparatorModule.cc.

terminate() [21/23]

void terminate ( void  )

overridevirtual

Termination action.

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

Reimplemented from Module.

Definition at line 220 of file TOPTimeBaseCalibratorModule.cc.

terminate() [22/23]

void terminate ( void  )

overridevirtual

Termination action.

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

Reimplemented from Module.

Definition at line 1254 of file TOPUnpackerModule.cc.

terminate() [23/23]

void terminate ( void  )

overridevirtual

Termination action.

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

Reimplemented from Module.

Definition at line 144 of file TOPWaveformFeatureExtractorModule.cc.

unpackHeadersInterimFEVer01()

bool unpackHeadersInterimFEVer01 ( const int *  buffer, int  bufferSize, bool  swapBytes  )

private

Tries to unpack raw data assuming it is in feature-extraction interim format.

Does not write out anything, just checks integrity.

Parameters

buffer	raw data buffer
bufferSize	buffer size
swapBytes	if true, swap bytes in buffer

Returns

true if buffer resembles interim format, false if not.

Definition at line 345 of file TOPUnpackerModule.cc.

unpackInterimFEVer01()

int unpackInterimFEVer01 ( const int *  buffer, int  bufferSize, bool  pedestalSubtracted  )

private

Unpack raw data given in feature-extraction interim format.

Parameters

buffer	raw data buffer
bufferSize	buffer size
pedestalSubtracted	false for version 2, true for version 3

Returns

number of words remaining in data buffer

Definition at line 415 of file TOPUnpackerModule.cc.

unpackProdDebug()

int unpackProdDebug ( const int *  buffer, int  bufferSize, TOP::RawDataType  dataFormat, bool  pedestalSubtracted  )

private

Unpack raw data given in production debugging format.

Parameters

buffer	raw data buffer
bufferSize	buffer size
dataFormat	data format
pedestalSubtracted	true, if pedestal is subtracted in waveforms

Returns

number of words remaining in data buffer

Definition at line 821 of file TOPUnpackerModule.cc.

unpackProductionDraft()

void unpackProductionDraft ( const int *  buffer, int  bufferSize  )

private

Unpack raw data given in a tentative production format (will vanish in future)

Parameters

buffer	raw data buffer
bufferSize	buffer size

Definition at line 243 of file TOPUnpackerModule.cc.

unpackType0Ver16()

void unpackType0Ver16 ( const int *  buffer, int  bufferSize  )

private

Unpack raw data given in feature-extraction production format.

Parameters

buffer	raw data buffer
bufferSize	buffer size

Definition at line 279 of file TOPUnpackerModule.cc.