SVDELS3Charge Class Reference

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

#include <SVDELS3Charge.h>

Inheritance diagram for SVDELS3Charge:

Collaboration diagram for SVDELS3Charge:

Public Member Functions
virtual	~SVDELS3Charge ()
	virtual destructor
void	computeClusterCharge (Belle2::SVD::RawCluster &rawCluster, double &charge, double &SNR, double &seedCharge) override
	compute the cluster charge, charge error and SNR with ELS3
void	applyMaxSampleCharge (const Belle2::SVD::RawCluster &rawCluster, double &charge, double &SNR, double &seedCharge)
	MaxSample Charge Algorithm.
void	applySumSamplesCharge (const Belle2::SVD::RawCluster &rawCluster, double &charge, double &SNR, double &seedCharge)
	SumSamples Charge Algorithm.
void	applyELS3Charge (const Belle2::SVD::RawCluster &rawCluster, double &charge, double &SNR, double &seedCharge)
	ELS3 Charge Algorithm.

Protected Attributes
DBObjPtr< HardwareClockSettings >	m_hwClock
	Hardware Clocks.
double	m_apvClockPeriod = 1. / m_hwClock->getClockFrequency(Const::EDetector::SVD, "sampling")
	APV clock period.
SVDPulseShapeCalibrations	m_PulseShapeCal
	SVDPulseShaper calibration wrapper.
SVDNoiseCalibrations	m_NoiseCal
	SVDNoise calibration wrapper.

Detailed Description

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

ELS3: Least-Squares method with the Exponentially decaying waveform.

Definition at line 26 of file SVDELS3Charge.h.

Constructor & Destructor Documentation

~SVDELS3Charge()

virtual ~SVDELS3Charge ( )

inlinevirtual

virtual destructor

Definition at line 33 of file SVDELS3Charge.h.

{};

Member Function Documentation

applyELS3Charge()

void applyELS3Charge ( const Belle2::SVD::RawCluster &  rawCluster, double &  charge, double &  SNR, double &  seedCharge  )

inherited

ELS3 Charge Algorithm.

Definition at line 95 of file SVDClusterCharge.cc.

    {
 
      // ISSUES:
      // 1. samples always in electrons for charge
      // 2. hardcoded ELS3 tau
      // 3. seed charge is maxSample
 
      float   m_ELS3tau = 55;
 
      //get seed charge
      int seedCellID = rawCluster.getStripsInRawCluster().at(rawCluster.getSeedInternalIndex()).cellID;
 
      seedCharge = m_PulseShapeCal.getChargeFromADC(rawCluster.getSensorID(), rawCluster.isUSide(), seedCellID,
                                                    rawCluster.getSeedMaxSample());
 
      //initialize noise and charge
      double noise = 0;
      charge = 0;
 
      //take the MaxSum 3 samples
      SVDMaxSumAlgorithm maxSum = SVDMaxSumAlgorithm(rawCluster.getClsSamples(false));
      std::vector<float> selectedSamples = maxSum.getSelectedSamples();
 
      auto begin = selectedSamples.begin();
 
      const double E = std::exp(- m_apvClockPeriod / m_ELS3tau);
      const double E2 = E * E;
      const double E3 = E * E * E;
      const double E4 = E * E * E * E;
      double a0 = (*begin);
      double a1 = (*(begin + 1));
      double a2 = (*(begin + 2));
 
      //compute raw time
      const double w = (a0 -  E2 * a2) / (2 * a0 + E * a1);
      auto rawtime_num  = 2 * E4 + w * E2;
      auto rawtime_den =  1 - E4 - w * (2 + E2);
      float rawtime = - m_apvClockPeriod * rawtime_num / rawtime_den;
 
      a0 = m_PulseShapeCal.getChargeFromADC(rawCluster.getSensorID(), rawCluster.isUSide(), rawCluster.getStripsInRawCluster()[0].cellID,
                                            a0);
      a1 = m_PulseShapeCal.getChargeFromADC(rawCluster.getSensorID(), rawCluster.isUSide(), rawCluster.getStripsInRawCluster()[0].cellID,
                                            a1);
      a2 = m_PulseShapeCal.getChargeFromADC(rawCluster.getSensorID(), rawCluster.isUSide(), rawCluster.getStripsInRawCluster()[0].cellID,
                                            a2);
 
 
      double num = (1. / E - E3) * a1 + (2 + E2) * a2 - (1 + 2 * E2) * a0;
      double den =  m_apvClockPeriod / m_ELS3tau * std::exp(1 + rawtime / m_ELS3tau) * (1 + 4 * E2 + E4);
 
      charge = num / den;
 
      //compute Noise
      std::vector<Belle2::SVD::StripInRawCluster> strips = rawCluster.getStripsInRawCluster();
 
      for (int i = 0; i < (int)strips.size(); i++) {
 
        Belle2::SVD::StripInRawCluster strip = strips.at(i);
 
        double tmp_noise = m_PulseShapeCal.getChargeFromADC(rawCluster.getSensorID(), rawCluster.isUSide(), strip.cellID, strip.noise);
        noise += tmp_noise * tmp_noise;
      }
 
      SNR = charge / sqrt(noise);
 
    }

applyMaxSampleCharge()

void applyMaxSampleCharge ( const Belle2::SVD::RawCluster &  rawCluster, double &  charge, double &  SNR, double &  seedCharge  )

inherited

MaxSample Charge Algorithm.

Definition at line 24 of file SVDClusterCharge.cc.

    {
 
      //get seed charge
      int seedCellID = rawCluster.getStripsInRawCluster().at(rawCluster.getSeedInternalIndex()).cellID;
 
      seedCharge = m_PulseShapeCal.getChargeFromADC(rawCluster.getSensorID(), rawCluster.isUSide(), seedCellID,
                                                    rawCluster.getSeedMaxSample());
 
      std::vector<Belle2::SVD::StripInRawCluster> strips = rawCluster.getStripsInRawCluster();
 
      //initialize noise and charge
      double noise = 0;
      charge = 0;
 
 
      for (int i = 0; i < (int)strips.size(); i++) {
 
        Belle2::SVD::StripInRawCluster strip = strips.at(i);
 
        double rawCharge = *std::max_element(begin(strip.samples), end(strip.samples));
 
        // calibrate (ADC -> electrons)
        double stripCharge = m_PulseShapeCal.getChargeFromADC(rawCluster.getSensorID(), rawCluster.isUSide(), strip.cellID, rawCharge);
 
        double tmp_noise = m_PulseShapeCal.getChargeFromADC(rawCluster.getSensorID(), rawCluster.isUSide(), strip.cellID, strip.noise);
        noise += tmp_noise * tmp_noise;
 
        charge += stripCharge;
      }
 
      SNR = charge / sqrt(noise);
    }

applySumSamplesCharge()

void applySumSamplesCharge ( const Belle2::SVD::RawCluster &  rawCluster, double &  charge, double &  SNR, double &  seedCharge  )

inherited

SumSamples Charge Algorithm.

Definition at line 59 of file SVDClusterCharge.cc.

    {
 
      std::vector<Belle2::SVD::StripInRawCluster> strips = rawCluster.getStripsInRawCluster();
 
      //initialize noise and charge
      double noise = 0;
      charge = 0;
      seedCharge = 0;
 
      for (int i = 0; i < (int)strips.size(); i++) {
 
        Belle2::SVD::StripInRawCluster strip = strips.at(i);
 
        double rawCharge = 0;
        for (auto sample : strip.samples)
          rawCharge += sample;
 
        // calibrate (ADC -> electrons)
        double stripCharge = m_PulseShapeCal.getChargeFromADC(rawCluster.getSensorID(), rawCluster.isUSide(), strip.cellID, rawCharge);
 
        if (stripCharge > seedCharge)
          seedCharge = stripCharge;
 
        double tmp_noise = m_PulseShapeCal.getChargeFromADC(rawCluster.getSensorID(), rawCluster.isUSide(), strip.cellID, strip.noise);
        noise += tmp_noise * tmp_noise;
 
        charge += stripCharge;
      }
 
      SNR = charge / sqrt(noise);
 
    }

computeClusterCharge()

void computeClusterCharge ( Belle2::SVD::RawCluster &  rawCluster, double &  charge, double &  SNR, double &  seedCharge  )

overridevirtual

compute the cluster charge, charge error and SNR with ELS3

Implements SVDClusterCharge.

Definition at line 25 of file SVDELS3Charge.cc.

    {
      applyELS3Charge(rawCluster, charge, SNR, seedCharge);
    }

Member Data Documentation

m_apvClockPeriod

double m_apvClockPeriod = 1. / m_hwClock->getClockFrequency(Const::EDetector::SVD, "sampling")

protectedinherited

APV clock period.

Definition at line 59 of file SVDClusterCharge.h.

m_hwClock

DBObjPtr<HardwareClockSettings> m_hwClock

protectedinherited

Hardware Clocks.

Definition at line 56 of file SVDClusterCharge.h.

m_NoiseCal

SVDNoiseCalibrations m_NoiseCal

protectedinherited

SVDNoise calibration wrapper.

Definition at line 65 of file SVDClusterCharge.h.

m_PulseShapeCal

SVDPulseShapeCalibrations m_PulseShapeCal

protectedinherited

SVDPulseShaper calibration wrapper.

Definition at line 62 of file SVDClusterCharge.h.

---

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

• svd/reconstruction/include/SVDELS3Charge.h
• svd/reconstruction/src/SVDELS3Charge.cc