Namespace to encapsulate code needed for simulation and reconstrucion of the SVD. More...

Classes
class	DefaultWave
	A functor to provide a simple model of APV25 strip response. More...

class	EmpiricalDistributionFunction
	Empirical distribution function object is basic for mainpulation of probabilities. More...

class	GeoSVDCreator
	The creator for the SVD geometry of the Belle II detector. More...

class	NNWaveFitter
	The class uses a neural network to find a probability distribution of arrival times for a sextet of APX25 signal samples. More...

class	NNWaveFitTool
	The class holds arrays of bins and bin centers, and a wave generator object containing information on the waveform function. More...

class	NoiseMap
	Base Class to represent strip-dependent noise map. More...

class	RawCluster
	Class representing a raw cluster candidate during clustering of the SVD. More...

class	SensorInfo
	Specific implementation of SensorInfo for SVD Sensors which provides additional sensor specific information. More...

class	SimpleClusterCandidate
	Class representing a cluster candidate during simple clustering of the SVD. More...

struct	stripInCluster
	structure containing the relevant informations of eachstrip of the cluster More...

struct	StripInRawCluster
	structure containing the relevant informations of each strip of the raw cluster More...

class	SVDBackgroundModule
	SVD Background module. More...

class	SVDBeamBackHitFilterModule
	The SVDBeamBackHitFilter module. More...

class	SVDClusterCharge
	Abstract Class representing the SVD cluster charge. More...

class	SVDClusterizerDirectModule
	SVD Direct Clusterizer. More...

class	SVDClusterizerModule
	The SVD Clusterizer. More...

class	SVDClusterPosition
	Abstract Class representing the SVD cluster position. More...

class	SVDClusterTime
	Abstract Class representing the SVD cluster time. More...

class	SVDCoG3Time
	Derived Class representing the SVD cluster time computed with the CoG3 algorithm. More...

class	SVDCoG6Time
	Derived Class representing the SVD cluster time computed with the CoG6 algorithm. More...

class	SVDCoGOnlyPosition
	Derived Class representing the SVD cluster position computed with the CoGOnly algorithm. More...

class	SVDDatabaseTestModule
	Return the calibration results for the noises and the constants measured during the svd local runs (charge, ADC pulse, peaking time, pulse width). More...

class	SVDDigitizerModule
	The SVD Digitizer module. More...

class	SVDDQMDoseModule
	The SVD dose-monitoring DQM module. More...

class	SVDDQMInjectionModule
	The SVD Occupancy after Injection DQM module. More...

class	SVDELS3Charge
	Derived Class representing the SVD cluster charge computed with the ELS3 algorithm. More...

class	SVDELS3Time
	Derived Class representing the SVD cluster time computed with the ELS3 algorithm. More...

class	SVDMaxSampleCharge
	Derived Class representing the SVD cluster charge computed summing the max sample of each strip. More...

class	SVDMaxSumAlgorithm
	Class implementing the MaxSum algorithm. More...

class	SVDMissingAPVsClusterCreatorModule
	SVDMissingAPVsClusterCreatorModule: The SVD MissingAPVsClusterCreator. More...

class	SVDNNClusterizerModule
	SVD NN Clusterizer. More...

class	SVDNNShapeReconstructorModule
	The SVD NNShapeReconstructor. More...

class	SVDOldDefaultPosition
	Derived Class representing the SVD cluster position computed with the old algorithm (up to release-05). More...

class	SVDPackerModule
	SVDPackerModule: The SVD Raw Hits Creator. More...

class	SVDRecoChargeFactory
	Cluster Charge Factory Class. More...

class	SVDRecoDigitCreatorModule
	The SVD RecoDigit Creator. More...

class	SVDReconstructionBase
	Class to check whether the reconstruction algorithms are available or not. More...

class	SVDRecoPositionFactory
	Cluster Position Factory Class. More...

class	SVDRecoTimeFactory
	Cluster Time Factory Class. More...

class	SVDShaperDigitSorterModule
	The SVDShaperDigitSorter module. More...

class	SVDSimpleClusterizerModule
	SVDSimpleClusterizerModule: The SVD SimpleClusterizer. More...

class	SVDSumSamplesCharge
	Derived Class representing the SVD cluster charge computed summing the samples of each strip. More...

class	SVDUnpackerModule
	SVDUnpackerModule: The SVD Raw Hits Decoder. More...

class	SVDWaveform
	The SVD waveform class. More...

class	TauEncoder
	Encoder/decoder for neural network tau values. More...

class	WaveFitter
	Waveform fitter class. More...

class	WaveGenerator
	Waveform generator This is a functor to calculate APV samples from waveform. More...

Typedefs
typedef std::map< short int, SVDWaveform >	StripWaveforms
	Map of all channels' waveforms in one sensor side.

typedef std::pair< StripWaveforms, StripWaveforms >	SensorWaveforms
	Waveforms of u- and v- channels in one sensor.

typedef std::map< VxdID, SensorWaveforms >	Waveforms
	Map of all waveforms in all sensors.

typedef std::vector< double >	nnFitterBinData
	Vector of values defined for bins, such as bin times or bin probabilities.

typedef std::vector< double >	nnFitterBins
	Vector of bin edges, nnFitterBinData.size() + 1.

typedef VXD::SensitiveDetector< SVDSimHit, SVDTrueHit >	SensitiveDetector
	The SVD Sensitive Detector class.

typedef double	apvSampleBaseType
	Vector of input samples in float form.

typedef std::array< apvSampleBaseType, nAPVSamples >	apvSamples
	vector od apvSample BaseType objects

typedef std::function< double(double)>	WaveformShape
	WaveformShape type.

Functions
double	w_expo (double t)
	Gamma waveform shape, x.exp(-x) This is only historically useful.

double	w_poly3 (double t)
	Polynomial waveform shape, x.

double	w_betaprime (double t)
	Beta-prime waveform shape, x^alpha/(1+x)^beta.

double	w_adjacentU (double t)
	Adjacent-channel waveform U-side.

double	w_adjacentV (double t)
	Adjacent-channel waveform V-side.

template<typename T >
void	zeroSuppress (T &a, double thr)
	pass zero suppression

template<typename T >
bool	pass3Samples (const T &a, double thr)
	pass 3-samples

double	tau_raw2real (double raw_tau)
	Convert Hao's raw tau (integral, in latency units) to correct betaprime scale.

	TEST (NNTimeFitter, DISABLED_CompareNetworkCoefficient)
	Compare NN fitter parameter dump from Python with internal data representation in C++ to verify that the network coefficients were transferred without a significant loss of accuracy.

	TEST (NNTimeFitter, DISABLED_CompareFits)
	Read a sample of fits from the Python NN suite and check that the C++ NN fitter produces the same results.

	TEST (SVDModeByte, Constructor)
	Check that both constructors yield the same result.

	TEST (SVDModeByte, getters)
	Check getters and setters.

	TEST (SVDRecoDigit, ConstructEmpty)
	Check empty object creation and data getters.

	TEST (SVDRecoDigit, ConstructFromContainer)
	Check object creation and data getters using a stl container of probabilities.

	TEST (SVDRecoDigit, ToString)
	Check object creation and data getters using a stl container of probabilities.

	TEST (SVDRecoTimeBase, ConstructEmpty)
	Check empty object creation and data getters.

	TEST (SVDRecoTimeBase, ConstructFromContainer)
	Check object creation and data getters.

	TEST (SVDRecoTimeBase, toString)
	Check object printout.

	TEST (SVDShaperDigit, ConstructEmpty)
	Check empty object creation and data getters.

	TEST (SVDShaperDigit, ConstructFromContainer)
	Check object creation and data getters using a stl container of samples.

	TEST (SVDShaperDigit, ConstructFromCArray)
	Check object creation and data getters using a c-array of samples.

	TEST (SVDShaperDigit, ConstructDefaultTimeMode)
	Check standard object creation without FADC time and mode information.

	TEST (SVDShaperDigit, SampleTrimming)
	test sample trimming

	TEST (SVDShaperDigit, TrimmingFunction)
	test sample trimming

	TEST (SVDWaveform, Getters)
	Check object creation and simple object getters.

	TEST (SVDWaveform, Constructors)
	Check the different constructors.

	TEST (SVDWaveform, Assignment)
	Check assignment operators.

	TEST (SVDWaveform, Waveform)
	Check the waveform values and relations.

	TEST (SVDWaveform, toString)
	Check the waveform string representation.

	TEST (TNiel, output)
	Check correct functioning of the class.

Variables
geometry::CreatorFactory< GeoSVDCreator >	GeoSVDFactory ("SVDCreator")
	Register the creator.

const std::size_t	nAPVSamples = 6
	Number of APV samples.

const double	dt_APV = 31.44 * Unit::ns
	APV sampling time.

const apvSamples	apvTimeBase
	APV time base - times for the 6 signals.

Detailed Description

Namespace to encapsulate code needed for simulation and reconstrucion of the SVD.

Typedef Documentation

◆ apvSampleBaseType

typedef double apvSampleBaseType

Vector of input samples in float form.

Definition at line 38 of file SVDSimulationTools.h.

◆ apvSamples

typedef std::array<apvSampleBaseType, nAPVSamples> apvSamples

vector od apvSample BaseType objects

Definition at line 39 of file SVDSimulationTools.h.

◆ nnFitterBinData

typedef std::vector<double> nnFitterBinData

Vector of values defined for bins, such as bin times or bin probabilities.

Definition at line 30 of file NNWaveFitTool.h.

◆ nnFitterBins

typedef std::vector<double> nnFitterBins

Vector of bin edges, nnFitterBinData.size() + 1.

Definition at line 33 of file NNWaveFitTool.h.

◆ SensitiveDetector

typedef VXD::SensitiveDetector<SVDSimHit, SVDTrueHit> SensitiveDetector

The SVD Sensitive Detector class.

Definition at line 24 of file SensitiveDetector.h.

◆ SensorWaveforms

typedef std::pair<StripWaveforms, StripWaveforms> SensorWaveforms

Waveforms of u- and v- channels in one sensor.

Definition at line 44 of file SVDDigitizerModule.h.

◆ StripWaveforms

typedef std::map<short int, SVDWaveform> StripWaveforms

Map of all channels' waveforms in one sensor side.

Definition at line 41 of file SVDDigitizerModule.h.

◆ Waveforms

typedef std::map<VxdID, SensorWaveforms> Waveforms

Map of all waveforms in all sensors.

Definition at line 47 of file SVDDigitizerModule.h.

◆ WaveformShape

typedef std::function<double(double)> WaveformShape

WaveformShape type.

This is the type for a naked waveform function, giving a single value for a properly scaled arguemnt. The functions are scaled to 1.0 at mode, location (shift) 0.0, and scale 1.0.

Definition at line 54 of file SVDSimulationTools.h.

Function Documentation

◆ pass3Samples()

bool pass3Samples	(	const T &	a,
		double	thr
	)

inline

pass 3-samples

Definition at line 193 of file SVDSimulationTools.h.

    {
      return (search_n(a.begin(), a.end(), 3, thr, std::greater<double>()) != a.end());
    }

◆ tau_raw2real()

double tau_raw2real ( double raw_tau )

inline

Convert Hao's raw tau (integral, in latency units) to correct betaprime scale.

Includes scaling and fit adjustmenta

Definition at line 205 of file SVDSimulationTools.h.

205{ return 7.32313 * raw_tau; }

◆ TEST() [1/22]

TEST	(	NNTimeFitter	,
		DISABLED_CompareFits
	)

Read a sample of fits from the Python NN suite and check that the C++ NN fitter produces the same results.

This is the result of select_t0_fits.columns: Index([ 'test', 'amplitude', 't0', 'tau', 'sigma', 's1', 's2', 's3', 's4', 's5', 's6', 'normed_tau', 't0_bin', 'abin', 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 'a_fit', 'a_sigma', 'chi2_ndf', 't_fit', 't_sigma', 't_true', 'tau', 'a_true', 't_bin', 'a_bin', 't_diff','a_diff'], dtype='object') Columns labeled by numerals contain probabilities.

Definition at line 53 of file NNFitterTest.cc.

    {
      const size_t max_lines = 100; // maximum number of lines to be read
 
      // Create an instance of the NN fitter and the fitter tool.
      NNWaveFitter fitter("SVDTimeNet_6samples");
      auto fitTool = fitter.getFitTool();
      size_t nProbs = fitTool.getBinCenters().size();
 
      ifstream infile("svd/data/test_sample.csv");
 
      // Read the rows one by one and compare results
      string line;
      getline(infile, line);
 
      for (size_t i_line = 0; i_line < max_lines; i_line++) {
 
        getline(infile, line);
        if (line.size() < 10) break;
        istringstream sline(line);
 
        // Parse header. We want the dimennsion of the probability array.
        // not needed
        string cell;
        getline(sline, cell, ','); // index
        getline(sline, cell, ','); // test
 
        // true values. Read from the file, though not used.
        getline(sline, cell, ',');
        [[maybe_unused]] double true_amp = stod(cell);
        getline(sline, cell, ',');
        [[maybe_unused]] double true_t0 = stod(cell);
        getline(sline, cell, ',');
        [[maybe_unused]] double width = stod(cell);
        getline(sline, cell, ',');
        [[maybe_unused]] double noise = stod(cell);
 
        // normalized samples
        apvSamples normedSamples;
        for (size_t iSample = 0; iSample < nAPVSamples; ++iSample) {
          getline(sline, cell, ',');
          normedSamples[iSample] = stod(cell);
        }
 
        // not needed
        getline(sline, cell, ',');
        getline(sline, cell, ',');
        getline(sline, cell, ',');
 
        // probabilities
        nnFitterBinData ProbsPy(nProbs);
        for (size_t iSample = 0; iSample < nProbs; ++iSample) {
          getline(sline, cell, ',');
          ProbsPy[iSample] = stod(cell);
        }
 
        // fit results
        getline(sline, cell, ',');
        double fitPy_amp = stod(cell);
        getline(sline, cell, ',');
        double fitPy_ampSigma = stod(cell);
        getline(sline, cell, ',');
        [[maybe_unused]] double fitPy_chi2 = stod(cell);
        getline(sline, cell, ',');
        double fitPy_t0 = stod(cell);
        getline(sline, cell, ',');
        double fitPy_t0Sigma = stod(cell);
 
        // now do the Cpp fit
        const shared_ptr<nnFitterBinData> ProbsCpp = fitter.getFit(normedSamples, width);
        for (size_t iBin = 0; iBin < nProbs; ++iBin)
          EXPECT_NEAR((*ProbsCpp)[iBin], ProbsPy[iBin], 5.0e-3);
 
        double t0_cpp, t0_err_cpp;
        tie(t0_cpp, t0_err_cpp) = fitTool.getTimeShift(*ProbsCpp);
        EXPECT_NEAR(t0_cpp, fitPy_t0, 5);
        EXPECT_NEAR(t0_err_cpp, fitPy_t0Sigma, 2);
 
        double amp_cpp, amp_err_cpp, chi2_cpp;
        tie(amp_cpp, amp_err_cpp, chi2_cpp) = fitTool.getAmplitudeChi2(normedSamples, t0_cpp, width);
        EXPECT_NEAR(amp_cpp, fitPy_amp, 1.0);
        EXPECT_NEAR(amp_err_cpp, fitPy_ampSigma, 0.1);
        // FIXME: This is calculated slightly differently in Python, and it shows.
        // EXPECT_NEAR(chi2_cpp, fitPy_chi2, 1);
      }
    }

◆ TEST() [2/22]

TEST	(	NNTimeFitter	,
		DISABLED_CompareNetworkCoefficient
	)

Compare NN fitter parameter dump from Python with internal data representation in C++ to verify that the network coefficients were transferred without a significant loss of accuracy.

FIXME (Oct 12,2017: This test fails on buildbot, but ran successfully on all machines I tried.

Definition at line 35 of file NNFitterTest.cc.

    {
      // Create an instance of the NN fitter
      NNWaveFitter fitter("svd/data/SVDTimeNet.xml");
      EXPECT_TRUE(fitter.checkCoefficients("svd/data/classifier.txt", 1.0e-6));
    }

◆ TEST() [3/22]

TEST	(	SVDModeByte	,
		Constructor
	)

Check that both constructors yield the same result.

Definition at line 23 of file SVDModeByte.cc.

    {
      // Default construct
      SVDModeByte s0;
      SVDModeByte s_def(151);
      SVDModeByte s_en(SVDRunType::zero_suppressed, SVDEventType::global_run,
                       SVDDAQModeType::daq_6samples, uint8_t(7));
      SVDModeByte s_num(2, 0, 2, 7);
 
      EXPECT_EQ(s0, s_def);
      EXPECT_EQ(s_def, s_en);
      EXPECT_EQ(s_en, s_num);
      // Empty
      EXPECT_EQ(uint8_t(151), s0.getID());
      EXPECT_EQ((string)s0, "0-suppr/global/6 samples/???");
      // Malformed
      s0.setTriggerBin(3);
      EXPECT_EQ((string)s0, "0-suppr/global/6 samples/3");
      s0.setDAQMode(1);
      EXPECT_EQ((string)s0, "0-suppr/global/3 samples/3");
    }

◆ TEST() [4/22]

TEST	(	SVDModeByte	,
		getters
	)

Check getters and setters.

Definition at line 46 of file SVDModeByte.cc.

    {
      // Construct
      SVDModeByte s(3, 1, 1, 3);
      EXPECT_EQ(s.getTriggerBin(), 3);
      EXPECT_EQ(s.getDAQMode(), SVDDAQModeType::daq_3samples);
      EXPECT_EQ(s.getEventType(), SVDEventType::local_run);
      EXPECT_EQ(s.getRunType(), SVDRunType::zero_suppressed_timefit);
      // Default
      SVDModeByte s0;
      EXPECT_EQ(s0.getTriggerBin(), 7);
      EXPECT_EQ(s0.getDAQMode(), SVDDAQModeType::daq_6samples);
    }

◆ TEST() [5/22]

TEST	(	SVDRecoDigit	,
		ConstructEmpty
	)

Check empty object creation and data getters.

Definition at line 26 of file SVDRecoDigit.cc.

    {
      // Create an empty 6-digit
      SVDRecoDigit digit;
      // Test getters
      EXPECT_EQ(VxdID{0}, digit.getSensorID());
      EXPECT_TRUE(digit.isUStrip());
      EXPECT_EQ(0, digit.getCellID());
      EXPECT_EQ(0, digit.getAmplitude());
      EXPECT_EQ(10, digit.getAmplitudeError());
      EXPECT_EQ(0, digit.getTime());
      EXPECT_EQ(100, digit.getTimeError());
      const SVDRecoDigit::OutputProbArray pdf = digit.getProbabilities();
      for (auto prob : pdf)
        EXPECT_EQ(SVDRecoDigit::OutputProbType(1.0), prob);
      EXPECT_EQ(100.0, digit.getChi2Ndf());
    }

◆ TEST() [6/22]

TEST	(	SVDRecoDigit	,
		ConstructFromContainer
	)

Check object creation and data getters using a stl container of probabilities.

Definition at line 47 of file SVDRecoDigit.cc.

    {
      // Create an arbitrary recodigit
      VxdID sensorID(3, 4, 1);
      short int cellID = 132;
      float init_charge = 23456;
      float init_chargeErr = 1234;
      float init_time = -16;
      float init_timeErr = 3.21;
      float init_chi2 = 2.34;
 
      std::vector<float> init_probs({0.0, 0.01, 0.10, 0.79, 0.06, 0.04});
      SVDRecoDigit digit(sensorID, false, cellID, init_charge, init_chargeErr, init_time,
                         init_timeErr, init_probs, init_chi2);
      // Test getters
      EXPECT_EQ(sensorID, digit.getSensorID());
      EXPECT_FALSE(digit.isUStrip());
      EXPECT_EQ(cellID, digit.getCellID());
      EXPECT_EQ(init_charge, digit.getCharge());
      EXPECT_EQ(init_chargeErr, digit.getAmplitudeError());
      EXPECT_EQ(init_time, digit.getTime());
      EXPECT_EQ(init_timeErr, digit.getTimeError());
      const SVDRecoDigit::OutputProbArray pdf = digit.getProbabilities();
      for (size_t ip = 0; ip < pdf.size(); ++ip)
        EXPECT_LE(fabs(static_cast<SVDRecoDigit::OutputProbType>(init_probs[ip]) - pdf[ip]), 1.0 / SVDRecoDigit::storedProbArrayNorm);
      // Also, norm must be 1
      float outputNorm = std::accumulate(pdf.begin(), pdf.end(), 0.0);
      EXPECT_LE(fabs(1.0 - outputNorm), 1.0e-6);
      EXPECT_EQ(init_chi2, digit.getChi2Ndf());
 
    }

◆ TEST() [7/22]

TEST	(	SVDRecoDigit	,
		ToString
	)

Check object creation and data getters using a stl container of probabilities.

Definition at line 81 of file SVDRecoDigit.cc.

    {
      // Create an arbitrary recodigit
      VxdID sensorID(3, 4, 1);
      short int cellID = 132;
      float init_charge = 23456;
      float init_chargeErr = 1234;
      float init_time = -16;
      float init_timeErr = 3.21;
      float init_chi2 = 2.34;
 
      std::vector<float> init_probs({0.0, 0.01, 0.10, 0.79, 0.06, 0.04});
      SVDRecoDigit digit(sensorID, false, cellID, init_charge, init_chargeErr, init_time,
                         init_timeErr, init_probs, init_chi2);
      // Test getters
      std::string
      digitString("VXDID : 25632 = 3.4.1 strip: V-132 Amplitude: 23456 +/- 1234 Time: -16 +/- 3.21\n probabilities: 0 655 6553 51772 3932 2621 Chi2/ndf: 2.34\n");
      EXPECT_EQ(digitString, digit.toString());
    }

◆ TEST() [8/22]

TEST	(	SVDRecoTimeBase	,
		ConstructEmpty
	)

Check empty object creation and data getters.

Definition at line 25 of file SVDRecoTimeBase.cc.

    {
      // Create an empty time base
      SVDRecoTimeBase timeBase;
      // Test getters
      EXPECT_EQ(VxdID{0}, timeBase.getSensorID());
      EXPECT_TRUE(timeBase.isUStrip());
      const SVDRecoTimeBase::BinEdgesArray bins = timeBase.getBins();
      EXPECT_EQ(0, timeBase.getNBins());
      for (auto  bin : bins)
        EXPECT_EQ(SVDRecoTimeBase::BinnedDataType(0.0), bin);
    }

◆ TEST() [9/22]

TEST	(	SVDRecoTimeBase	,
		ConstructFromContainer
	)

Check object creation and data getters.

Definition at line 41 of file SVDRecoTimeBase.cc.

    {
      // Create an arbitrary recotimeBase
      VxdID sensorID(3, 4, 1);
      std::vector<float> init_bins({ -5.0, -3.0, -1.0, 1.0, 3.0, 5.0});
      SVDRecoTimeBase timeBase(sensorID, false, init_bins);
      // Test getters
      EXPECT_EQ(sensorID, timeBase.getSensorID());
      EXPECT_FALSE(timeBase.isUStrip());
      EXPECT_EQ(init_bins.size() - 1, timeBase.getNBins());
      const SVDRecoTimeBase::BinEdgesArray& bins = timeBase.getBins();
      for (size_t ib = 0; ib < bins.size(); ++ib)
        EXPECT_EQ(init_bins[ib], bins[ib]);
    }

◆ TEST() [10/22]

TEST	(	SVDRecoTimeBase	,
		toString
	)

Check object printout.

Definition at line 58 of file SVDRecoTimeBase.cc.

    {
      // Create an arbitrary recotimeBase
      VxdID sensorID(3, 4, 1);
      std::vector<float> init_bins({ -5.0, -3.0, -1.0, 1.0, 3.0, 5.0});
      SVDRecoTimeBase timeBase(sensorID, false, init_bins);
      // Test toString method
      EXPECT_EQ("VXDID : 25632 = 3.4.1 side: V bins: -5 -3 -1 1 3 5 ", timeBase.toString());
    }

◆ TEST() [11/22]

TEST	(	SVDShaperDigit	,
		ConstructDefaultTimeMode
	)

Check standard object creation without FADC time and mode information.

Definition at line 90 of file SVDShaperDigit.cc.

    {
      // Create an arbitrary 6-digit
      VxdID sensorID(3, 4, 1);
      short int cellID = 132;
 
      double init_samples[6] = {0, 5, 10, 9, 6, 5};
      SVDShaperDigit digit(sensorID, false, cellID, init_samples);
      // Test that time and error are set correctly.
      EXPECT_EQ(0.0, digit.getFADCTime());
    }

◆ TEST() [12/22]

TEST	(	SVDShaperDigit	,
		ConstructEmpty
	)

Check empty object creation and data getters.

Definition at line 26 of file SVDShaperDigit.cc.

    {
      // Create an empty 6-digit
      SVDShaperDigit digit;
      // Test getters
      EXPECT_EQ(VxdID{0}, digit.getSensorID());
      EXPECT_TRUE(digit.isUStrip());
      EXPECT_EQ(0, digit.getCellID());
      const SVDShaperDigit::APVFloatSamples& samples = digit.getSamples();
      for (auto sample : samples)
        EXPECT_EQ(SVDShaperDigit::APVFloatSampleType(0), sample);
      EXPECT_EQ(0.0, digit.getFADCTime());
    }

◆ TEST() [13/22]

TEST	(	SVDShaperDigit	,
		ConstructFromCArray
	)

Check object creation and data getters using a c-array of samples.

Definition at line 66 of file SVDShaperDigit.cc.

    {
      // Create an arbitrary 6-digit
      VxdID sensorID(3, 4, 1);
      short int cellID = 132;
      char digitFADCTime(-16);
 
      // floats must work, too.
      float init_samples[SVDShaperDigit::c_nAPVSamples] = {0, 5, 10, 9, 6, 5};
      SVDShaperDigit digit(
        sensorID, false, cellID, init_samples, digitFADCTime);
      // Test getters
      EXPECT_EQ(sensorID, digit.getSensorID());
      EXPECT_FALSE(digit.isUStrip());
      EXPECT_EQ(cellID, digit.getCellID());
      const SVDShaperDigit::APVFloatSamples& samples = digit.getSamples();
      for (size_t isample = 0; isample < SVDShaperDigit::c_nAPVSamples; ++isample)
        EXPECT_EQ(static_cast<SVDShaperDigit::APVFloatSampleType>(init_samples[isample]),
                  samples[isample]);
      EXPECT_EQ(static_cast<float>(digitFADCTime), digit.getFADCTime());
    }

◆ TEST() [14/22]

TEST	(	SVDShaperDigit	,
		ConstructFromContainer
	)

Check object creation and data getters using a stl container of samples.

Definition at line 43 of file SVDShaperDigit.cc.

    {
      // Create an arbitrary 6-digit
      VxdID sensorID(3, 4, 1);
      short int cellID = 132;
      char digitFADCTime(-16);
 
      std::vector<int> init_samples({0, 5, 10, 9, 6, 5});
      SVDShaperDigit digit(sensorID, false, cellID, init_samples, digitFADCTime);
      // Test getters
      EXPECT_EQ(sensorID, digit.getSensorID());
      EXPECT_FALSE(digit.isUStrip());
      EXPECT_EQ(cellID, digit.getCellID());
      const SVDShaperDigit::APVFloatSamples& samples = digit.getSamples();
      for (size_t isample = 0; isample < SVDShaperDigit::c_nAPVSamples; ++isample)
        EXPECT_EQ(static_cast<SVDShaperDigit::APVFloatSampleType>(init_samples[isample]),
                  samples[isample]);
      EXPECT_EQ(static_cast<float>(digitFADCTime), digit.getFADCTime());
    }

◆ TEST() [15/22]

TEST	(	SVDShaperDigit	,
		SampleTrimming
	)

test sample trimming

Definition at line 103 of file SVDShaperDigit.cc.

    {
      // Create an arbitrary 6-digit
      VxdID sensorID(3, 4, 1);
      short int cellID = 132;
 
      auto sampleMin =
        static_cast<SVDShaperDigit::APVFloatSampleType>(std::numeric_limits<SVDShaperDigit::APVRawSampleType>::lowest());
      auto sampleMax =
        static_cast<SVDShaperDigit::APVFloatSampleType>(std::numeric_limits<SVDShaperDigit::APVRawSampleType>::max());
 
      std::vector<float> init_samples({0, -10, 255, 256, 500, -1});
      SVDShaperDigit digit(sensorID, false, cellID, init_samples);
      const auto& samples = digit.getSamples();
      for (size_t isample = 0; isample < SVDShaperDigit::c_nAPVSamples; ++isample) {
        auto trimmedSample = std::max(sampleMin, std::min(sampleMax, init_samples[isample]));
        EXPECT_EQ(trimmedSample, samples[isample]);
      }
    }

◆ TEST() [16/22]

TEST	(	SVDShaperDigit	,
		TrimmingFunction
	)

test sample trimming

Definition at line 124 of file SVDShaperDigit.cc.

    {
 
      // Create an arbitrary 6-digit
      auto sampleMin =
        static_cast<SVDShaperDigit::APVFloatSampleType>(std::numeric_limits<SVDShaperDigit::APVRawSampleType>::lowest());
      auto sampleMax =
        static_cast<SVDShaperDigit::APVFloatSampleType>(std::numeric_limits<SVDShaperDigit::APVRawSampleType>::max());
 
      std::vector<float> init_samples({0, -10, 255, 256, 500, -1});
      for (size_t isample = 0; isample < SVDShaperDigit::c_nAPVSamples; ++isample) {
        auto trimmedSample = std::max(sampleMin, std::min(sampleMax, init_samples[isample]));
        EXPECT_EQ(trimmedSample, SVDShaperDigit::trimToSampleRange(init_samples[isample]));
      }
    }

◆ TEST() [17/22]

TEST	(	SVDWaveform	,
		Assignment
	)

Check assignment operators.

Definition at line 94 of file SVDWaveform.cc.

    {
      // Create an artificial waveform and load it with data.
      vector<double> charges;
      charges.push_back(10);
      charges.push_back(11);
      charges.push_back(30);
      vector<double> times;
      for (int i = 0; i < 3; ++i) times.push_back(i * 15.0e-9);
      vector<RelationElement::index_type> mcParticles;
      for (int i = 0; i < 3; ++i) mcParticles.push_back(3 * i);
      vector<RelationElement::index_type> trueHits;
      for (int i = 0; i < 3; ++i) trueHits.push_back(7 * i);
      const float tau = 25.0e-9;
      SVDWaveform waveform;
      for (int i = 0; i < 3; ++i) waveform.add(times.at(i), charges.at(i), tau, mcParticles.at(i), trueHits.at(i));
      // Test the assignment operator
      SVDWaveform waveform2 = waveform;
      double totalcharge = 0;
      for (int i = 0; i < 3; ++i) totalcharge += charges.at(i);
      EXPECT_EQ(waveform2.getCharge(), totalcharge);
      //Test SVDWaveform::getElementaryWaveforms()
      SVDWaveform::elementary_waveform_list functions = waveform2.getElementaryWaveforms();
      for (int i = 0; i < 3; ++i) {
        EXPECT_EQ(functions.at(i).m_initTime, times.at(i));
        EXPECT_EQ(functions.at(i).m_charge, charges.at(i));
        EXPECT_EQ(functions.at(i).m_tau, tau);
        EXPECT_EQ(functions.at(i).m_particle, mcParticles.at(i));
        EXPECT_EQ(functions.at(i).m_truehit, trueHits.at(i));
      }
    }

◆ TEST() [18/22]

TEST	(	SVDWaveform	,
		Constructors
	)

Check the different constructors.

Definition at line 59 of file SVDWaveform.cc.

    {
      // Create an artificial waveform and load it with data.
      vector<double> charges;
      charges.push_back(10);
      charges.push_back(11);
      charges.push_back(30);
      vector<double> times;
      for (int i = 0; i < 3; ++i) times.push_back(i * 15.0e-9);
      vector<RelationElement::index_type> mcParticles;
      for (int i = 0; i < 3; ++i) mcParticles.push_back(3 * i);
      vector<RelationElement::index_type> trueHits;
      for (int i = 0; i < 3; ++i) trueHits.push_back(7 * i);
      const float tau = 25.0e-9;
      SVDWaveform waveform;
      for (int i = 0; i < 3; ++i) waveform.add(times.at(i), charges.at(i), tau, mcParticles.at(i), trueHits.at(i));
      // Test the copy constructor
      SVDWaveform waveform2(waveform);
      double totalcharge = 0;
      for (int i = 0; i < 3; ++i) totalcharge += charges.at(i);
      EXPECT_EQ(waveform2.getCharge(), totalcharge);
      //Test SVDWaveform::getElementaryWaveforms()
      SVDWaveform::elementary_waveform_list functions = waveform2.getElementaryWaveforms();
      size_t n = functions.size();
      EXPECT_EQ(n, 3);
      for (int i = 0; i < 3; ++i) {
        EXPECT_EQ(functions.at(i).m_initTime, times.at(i));
        EXPECT_EQ(functions.at(i).m_charge, charges.at(i));
        EXPECT_EQ(functions.at(i).m_tau, tau);
        EXPECT_EQ(functions.at(i).m_particle, mcParticles.at(i));
        EXPECT_EQ(functions.at(i).m_truehit, trueHits.at(i));
      }
    }

◆ TEST() [19/22]

TEST	(	SVDWaveform	,
		Getters
	)

Check object creation and simple object getters.

Definition at line 27 of file SVDWaveform.cc.

    {
      // First create an artificial waveform and load it with data.
      vector<double> charges;
      charges.push_back(10);
      charges.push_back(11);
      charges.push_back(30);
      vector<double> times;
      for (int i = 0; i < 3; ++i) times.push_back(i * 15.0e-9);
      vector<RelationElement::index_type> mcParticles;
      for (int i = 0; i < 3; ++i) mcParticles.push_back(3 * i);
      vector<RelationElement::index_type> trueHits;
      for (int i = 0; i < 3; ++i) trueHits.push_back(7 * i);
      const float tau = 25.0e-9;
      SVDWaveform waveform;
      for (int i = 0; i < 3; ++i) waveform.add(times.at(i), charges.at(i), tau, mcParticles.at(i), trueHits.at(i));
      // Test SVDWaveform::getCharge()
      double totalcharge = 0;
      for (int i = 0; i < 3; ++i) totalcharge += charges.at(i);
      EXPECT_EQ(waveform.getCharge(), totalcharge);
      //Test SVDWaveform::getElementaryWaveforms()
      SVDWaveform::elementary_waveform_list functions = waveform.getElementaryWaveforms();
      for (int i = 0; i < 3; ++i) {
        EXPECT_EQ(functions.at(i).m_initTime, times.at(i));
        EXPECT_EQ(functions.at(i).m_charge, charges.at(i));
        EXPECT_EQ(functions.at(i).m_tau, tau);
        EXPECT_EQ(functions.at(i).m_particle, mcParticles.at(i));
        EXPECT_EQ(functions.at(i).m_truehit, trueHits.at(i));
      }
    }

◆ TEST() [20/22]

TEST	(	SVDWaveform	,
		toString
	)

Check the waveform string representation.

Definition at line 156 of file SVDWaveform.cc.

    {
      // Create an artificial waveform and load it with data.
      vector<double> charges;
      charges.push_back(10);
      charges.push_back(11);
      charges.push_back(30);
      vector<double> times;
      for (int i = 0; i < 3; ++i) times.push_back(i * 15.0e-9);
      vector<RelationElement::index_type> mcParticles;
      for (int i = 0; i < 3; ++i) mcParticles.push_back(3 * i);
      vector<RelationElement::index_type> trueHits;
      for (int i = 0; i < 3; ++i) trueHits.push_back(7 * i);
      const float tau = 50.0e-9; // 25 ns shaping time
      SVDWaveform waveform;
      for (int i = 0; i < 3; ++i) waveform.add(times.at(i), charges.at(i), tau, mcParticles.at(i), trueHits.at(i));
      std::ostringstream os;
      for (int i = 0; i < 3; ++i)
        os << i + 1 << '\t' << times.at(i) << '\t' << charges.at(i) << '\t' << tau << std::endl;
      EXPECT_EQ(waveform.toString(), os.str());
    }

◆ TEST() [21/22]

TEST	(	SVDWaveform	,
		Waveform
	)

Check the waveform values and relations.

Definition at line 127 of file SVDWaveform.cc.

    {
      // Create an artificial waveform and load it with data.
      vector<double> charges;
      charges.push_back(10);
      charges.push_back(11);
      charges.push_back(30);
      vector<double> times;
      for (int i = 0; i < 3; ++i) times.push_back(i * 15.0e-9);
      vector<RelationElement::index_type> mcParticles;
      for (int i = 0; i < 3; ++i) mcParticles.push_back(3 * i);
      vector<RelationElement::index_type> trueHits;
      for (int i = 0; i < 3; ++i) trueHits.push_back(7 * i);
      const float tau = 25.0e-9; // 25 ns shaping time
      SVDWaveform waveform;
      for (int i = 0; i < 3; ++i) waveform.add(times.at(i), charges.at(i), tau, mcParticles.at(i), trueHits.at(i));
      // Test operator()
      double time = 40.0e-9; // 40 ns
      // expected function value
      double expected_value = 0;
      for (SVDWaveform::ElementaryWaveform wave : waveform.getElementaryWaveforms()) {
        double value = waveform.waveform(time, wave);
        expected_value += value;
      }
      double function_value = waveform(time);
      EXPECT_EQ(function_value, expected_value);
    }

◆ TEST() [22/22]

TEST	(	TNiel	,
		output
	)

Check correct functioning of the class.

Definition at line 25 of file niel_fun.cc.

    {
      // First create a niel function and load it with data.
      TNiel niel("vxd/tests/niel_test.csv");
      // Value from within the range
      EXPECT_EQ(niel.getNielFactor(3.0), 3.126E-2);
      // Underflow value
      EXPECT_EQ(niel.getNielFactor(0.0), 1.543E-2);
      // Overflow value
      EXPECT_EQ(niel.getNielFactor(20.0), 5.160E-2);
    }

◆ w_adjacentU()

double w_adjacentU ( double t )

inline

Adjacent-channel waveform U-side.

The exact parametrisation is to be determined.

Parameters

t	Properly scaled time, (t - t0)/tau.

Returns: Waveform value (TBD).

Definition at line 99 of file SVDSimulationTools.h.

    {
 
      double p[8];
 
      // waveform parameters
      p[0] = -7.379;
      p[1] = -0.3594;
      p[2] =  0.07564;
      p[3] =  7.391;
      p[4] = -0.08004;
      p[5] = -2.139e-5;
      p[6] = -1.781e-7;
      p[7] =  0.01827;
 
      if (t < 0.0)
        return 0.0;
      else
        return (p[0] + p[1] * t) * exp(-p[2] * t) + (p[3] + p[4] * t + p[5] * pow(t, 2) + p[6] * pow(t, 3)) * exp(-p[7] * t);
    }

◆ w_adjacentV()

double w_adjacentV ( double t )

inline

Adjacent-channel waveform V-side.

The exact parametrisation is to be determined.

Parameters

t	Properly scaled time, (t - t0)/tau.

Returns: Waveform value (TBD).

Definition at line 125 of file SVDSimulationTools.h.

    {
 
      double p[8];
 
      // waveform parameters
      p[0] =  0.2043;
      p[1] = -0.0197;
      p[2] =  0.01448;
      p[3] = -0.2076;
      p[4] =  0.0387;
      p[5] =  0.01055;
      p[6] = -4.03e-6;
      p[7] =  0.07193;
 
      if (t < 0.0)
        return 0.0;
      else
        return (p[0] + p[1] * t) * exp(-p[2] * t) + (p[3] + p[4] * t + p[5] * pow(t, 2) + p[6] * pow(t, 3)) * exp(-p[7] * t);
    }

◆ w_betaprime()

double w_betaprime ( double t )

inline

Beta-prime waveform shape, x^alpha/(1+x)^beta.

This is the function for general use.

Parameters

t	Properly scaled time, (t - t0)/tau.

Returns: Waveform value 149.012 * t^2 / (1+t)^10 for t > 0, else 0

Definition at line 86 of file SVDSimulationTools.h.

    {
      if (t < 0.0)
        return 0.0;
      else
        return 149.012 * pow(t, 2) * pow(1.0 + t, -10);
    }

◆ w_expo()

double w_expo ( double t )

inline

Gamma waveform shape, x.exp(-x) This is only historically useful.

Use beta-prime instead.

Parameters

t	properly scaled time, (t-t0)/tau

Returns: Waveform value t * exp(1 - t) for t > 0, else 0

Definition at line 61 of file SVDSimulationTools.h.

    {
      if (t < 0.0) return 0.0;
      else return t * exp(1.0 - t);
    }

◆ w_poly3()

double w_poly3 ( double t )

inline

Polynomial waveform shape, x.

(1-x)^2 This is used in signal calibration in SVD hardware tests. Do not use for signal simulation.

Parameters

t	properly scaled time, (t-t0)/tau

Returns: Waveform value 27/4 * t * (1 - t)^2 for t in <0,1>, else 0

Definition at line 73 of file SVDSimulationTools.h.

    {
      if (t < 0.0 || t > 1.0)
        return 0.0;
      else
        return 6.75 * t * (1.0 - t) * (1.0 - t);
    }

◆ zeroSuppress()

void zeroSuppress	(	T &	a,
		double	thr
	)

inline

pass zero suppression

Definition at line 182 of file SVDSimulationTools.h.

    {
      std::replace_if(a.begin(), a.end(), std::bind(std::less<double>(), std::placeholders::_1, thr), 0.0);
    }

Variable Documentation

◆ apvTimeBase

const apvSamples apvTimeBase

Initial value:

=

Belle2::SVD::dt_APV

const double dt_APV

APV sampling time.

Definition: SVDSimulationTools.h:35

APV time base - times for the 6 signals.

Definition at line 42 of file SVDSimulationTools.h.

◆ dt_APV

const double dt_APV = 31.44 * Unit::ns

APV sampling time.

Definition at line 35 of file SVDSimulationTools.h.

◆ nAPVSamples

const std::size_t nAPVSamples = 6

Number of APV samples.

Definition at line 32 of file SVDSimulationTools.h.

Classes

Typedefs

Functions

Variables

Detailed Description

Typedef Documentation

◆ apvSampleBaseType

◆ apvSamples

◆ nnFitterBinData

◆ nnFitterBins

◆ SensitiveDetector

◆ SensorWaveforms

◆ StripWaveforms

◆ Waveforms

◆ WaveformShape

Function Documentation

◆ pass3Samples()

◆ tau_raw2real()

◆ TEST() [1/22]

◆ TEST() [2/22]

◆ TEST() [3/22]

◆ TEST() [4/22]

◆ TEST() [5/22]

◆ TEST() [6/22]

◆ TEST() [7/22]

◆ TEST() [8/22]

◆ TEST() [9/22]

◆ TEST() [10/22]

◆ TEST() [11/22]

◆ TEST() [12/22]

◆ TEST() [13/22]

◆ TEST() [14/22]

◆ TEST() [15/22]

◆ TEST() [16/22]

◆ TEST() [17/22]

◆ TEST() [18/22]

◆ TEST() [19/22]

◆ TEST() [20/22]

◆ TEST() [21/22]

◆ TEST() [22/22]

◆ w_adjacentU()

◆ w_adjacentV()

◆ w_betaprime()

◆ w_expo()

◆ w_poly3()

◆ zeroSuppress()

Variable Documentation

◆ apvTimeBase

◆ dt_APV

◆ nAPVSamples