beamlinereco::CFDHitFinder< T > Class Template Reference

#include "/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-01-21/BeamlineReco/CFDHitFinderAlg.h"

Public Member Functions
	CFDHitFinder ()
virtual	~CFDHitFinder ()=default
void	SetWaveform (std::vector< uint16_t > &waveform, unsigned int channelNo, uint32_t timestamp)
void	SetParams (const std::map< CFDParams, double > &paramSet)
void	SetDiscriminationType (const std::string type)
void	SetInterpolationType (const std::string type)
void	SetFilterType (const std::string type)
const std::string	GetDiscriminationType () const
const std::string	GetInterpolationType () const
const std::string	GetFilterType () const
const std::map< T, hit_t< T > > &	GetHitCollection () const
const T &	GetPedestal () const
const T &	GetNoiseSigma () const
const T &	GetDiscriminationThreshold () const
virtual void	Go ()

Private Member Functions
void	SetChannel (unsigned int channelNo)
void	SetTimestamp (uint32_t timestamp)
void	SetParam (CFDParams param, T value)
virtual void	FindPedestal ()
virtual void	FindRawHitLogic ()
virtual void	FindCFDHits ()
virtual bool	BackwardFindingOfHitStart (size_t hitPeakTimeIndex, T hitPeakValue, T &hitStartTimeIndex, T &hitRiseTimeInIndex)
virtual bool	ForwardFindingOfHitFallTime (size_t hitPeakTimeIndex, T &hitFallTimeInIndex)
virtual T	IntegrateWaveformInADC (size_t hitStartTimeIndex, size_t hitEndTimeIndex)
virtual T	FindHitWidth (size_t hitPeakTimeIndex, T hitPeakValueADC)
virtual T	FindPeakValue (size_t hitPeakTimeIndex, T &hitPeakValue)
virtual void	Reset ()

Private Attributes
unsigned int	_ChannelNumber
uint32_t	_Timestamp
std::string	_DiscriminationType
std::string	_InterpolationType
std::string	_FilterType
std::vector< T >	_WaveformADCNanosec
std::vector< T >	_NonfilterWaveformADCNanosec
std::vector< bool >	_RawHitLogicNanosec
T	_PedestalInADC
T	_NoiseSigmaInADC
T	_DiscriminationThresholdInADC
std::map< CFDParams, T >	_CFDParamSet
std::map< T, hit_t< T > >	_HitCollection

Detailed Description

template<class T>
class beamlinereco::CFDHitFinder< T >

Definition at line 71 of file CFDHitFinderAlg.h.

Constructor & Destructor Documentation

template<class T>

beamlinereco::CFDHitFinder< T >::CFDHitFinder ( )

inline

Definition at line 73 of file CFDHitFinderAlg.h.

                       {
            _HitCollection.clear();
            _RawHitLogicNanosec.clear();
            _PedestalInADC = 0;
            _NoiseSigmaInADC = 0;
            _DiscriminationThresholdInADC = 0;
            _WaveformADCNanosec.clear();
            _ChannelNumber = 993;
            _Timestamp = 0;
            _NonfilterWaveformADCNanosec.clear();
            _CFDParamSet.clear();
        }

template<class T>

virtual beamlinereco::CFDHitFinder< T >::~CFDHitFinder ( )

virtualdefault

Member Function Documentation

template<class T>

bool beamlinereco::CFDHitFinder< T >::BackwardFindingOfHitStart ( size_t  hitPeakTimeIndex, T  hitPeakValue, T &  hitStartTimeIndex, T &  hitRiseTimeInIndex  )

privatevirtual

Definition at line 452 of file CFDHitFinderAlg.h.

References art::errors::Configuration, allTimeWatchdog::endl, beamlinereco::kDiscriminationThreshold, PandAna.Demos.pi0_spectra::p0, plot_validation_datamc::p1, T, confusionMatrixTree::t, getGoodRuns4SAM::t1, t2, submit_syst::x, and submit_syst::y.

                                                                                                                                                {
    bool isThisHitContainedFromTheRisingEdge = true;

    if (_DiscriminationType == "CFD") {
      _DiscriminationThresholdInADC = _PedestalInADC - _CFDParamSet[kDiscriminationThreshold] * (_PedestalInADC - hitPeakValue);
    } else if (_DiscriminationType == "LE") {
      _DiscriminationThresholdInADC = _PedestalInADC - _CFDParamSet[kDiscriminationThreshold] * _NoiseSigmaInADC;
    } else {
      throw art::Exception(art::errors::Configuration)
        << "Unknown DiscriminationType: " << _DiscriminationType << "." << std::endl;
    }

    size_t tmp_index = hitPeakTimeIndex;
    while (_WaveformADCNanosec.at(tmp_index) < _DiscriminationThresholdInADC) {
        if (tmp_index == 2 && _WaveformADCNanosec.at(tmp_index - 1) < _DiscriminationThresholdInADC) {
            isThisHitContainedFromTheRisingEdge = false;
            break;
        }
        tmp_index = tmp_index - 1;
    }

    if ((_InterpolationType == "LinearNearest") || (tmp_index < 2)) {
      // Draw a line between two points
      T V1 = _WaveformADCNanosec.at(tmp_index);
      T V2 = _WaveformADCNanosec.at(tmp_index - 1);
      size_t t1 = tmp_index;
      size_t t2 = tmp_index - 1;
      T slope = (V2 - V1) / ((T)t2 - (T)t1);
      hitStartTimeIndex = (_DiscriminationThresholdInADC - V1) / slope + (T)t1;
    }
    if ((_InterpolationType == "LinearFit") && (tmp_index >= 2)) {
      // Fit the surrounding 4
      T x[4] = {
        (T)(tmp_index - 2),
        (T)(tmp_index - 1),
        (T)(tmp_index),
        (T)(tmp_index + 1)
      };
      T y[4] = {
        _WaveformADCNanosec.at(tmp_index - 2),
        _WaveformADCNanosec.at(tmp_index - 1),
        _WaveformADCNanosec.at(tmp_index),
        _WaveformADCNanosec.at(tmp_index + 1)
      };
      TGraph *gfit = new TGraph(4, x, y);
      TFitResultPtr frp = gfit->Fit("pol1", "SQNF");
      double p0 = frp->Parameter(0);
      double p1 = frp->Parameter(1);
      T V = p0 + p1 * (T)tmp_index;
      T t = (T)tmp_index;
      hitStartTimeIndex = (T)t + (_DiscriminationThresholdInADC - V) / p1;
    }

    if (!isThisHitContainedFromTheRisingEdge) {
        hitRiseTimeInIndex = (T) hitPeakTimeIndex - (T) tmp_index;
    } else {
        tmp_index = hitPeakTimeIndex;
        while (_WaveformADCNanosec.at(tmp_index) <= _PedestalInADC) {
            if (tmp_index == 2 && _WaveformADCNanosec.at(tmp_index - 1) < _PedestalInADC) {
                isThisHitContainedFromTheRisingEdge = false;
                break;
            }
            tmp_index = tmp_index - 1;
        }
        hitRiseTimeInIndex = (T) hitPeakTimeIndex - (T) tmp_index;
    }

    return isThisHitContainedFromTheRisingEdge;
}

template<class T >

void beamlinereco::CFDHitFinder< T >::FindCFDHits ( )

privatevirtual

Definition at line 330 of file CFDHitFinderAlg.h.

References beamlinereco::hit_t< T >::AmplitudeInADC, beamlinereco::hit_t< T >::AmplitudeInMiliVolt, beamlinereco::hit_t< T >::DigitizerChannel, beamlinereco::hit_t< T >::FallTimeInNanoSec, MECModelEnuComparisons::i, beamlinereco::hit_t< T >::IntegratedChargeInADCNanoSec, beamlinereco::hit_t< T >::IntegratedChargeInADCTimeTicks, beamlinereco::hit_t< T >::IsContained, calib::j, beamlinereco::kADCDynamicRange, beamlinereco::kADCNBits, beamlinereco::kRawHitFinderTicksFromEnd, beamlinereco::kTimeSamplingInterval, beamlinereco::hit_t< T >::RiseTimeInNanoSec, T, beamlinereco::hit_t< T >::Timestamp, beamlinereco::hit_t< T >::TPeakInNanoSec, beamlinereco::hit_t< T >::TStartInNanoSec, and beamlinereco::hit_t< T >::WidthInNanoSec.

                                              {
    std::vector<T> hitPeakValueVector; hitPeakValueVector.clear();
    std::vector<size_t> hitPeakTimeIndexVector; hitPeakTimeIndexVector.clear();
    std::vector<T> hitPeakTimeInNanoSecVector; hitPeakTimeInNanoSecVector.clear();
    std::vector<T> hitStartTimeIndexVector; hitStartTimeIndexVector.clear();
    std::vector<T> hitRiseTimeInIndexVector; hitRiseTimeInIndexVector.clear();
    std::vector<T> hitFallTimeInIndexVector; hitFallTimeInIndexVector.clear();
    std::vector<T> hitWidthValueVector; hitWidthValueVector.clear();
    std::vector<bool> isHitContainedVector; isHitContainedVector.clear();

    bool logicFallingPart = true;
    for (size_t i = 1; i < _RawHitLogicNanosec.size() - 1; i++) {
        T hitPeakValue = 1E10;
        size_t hitPeakTimeIndex = 0;
        if (_RawHitLogicNanosec.at(i) && !(_RawHitLogicNanosec.at(i - 1)) && logicFallingPart) {
            logicFallingPart = false;
            for (size_t j = i; j < _RawHitLogicNanosec.size() - 1; j++) {
                hitPeakValue = hitPeakValue > _WaveformADCNanosec.at(j) ? (T) _WaveformADCNanosec.at(j) : hitPeakValue;
                hitPeakTimeIndex = hitPeakValue == _WaveformADCNanosec.at(j) ? j : hitPeakTimeIndex;
                if (((_RawHitLogicNanosec.at(j) && !(_RawHitLogicNanosec.at(j + 1))) || ((_RawHitLogicNanosec.size() - j) < _CFDParamSet[kRawHitFinderTicksFromEnd])) && !logicFallingPart) {
                    logicFallingPart = true;
                    T hitStartTimeIndex = (T)hitPeakTimeIndex;
                    T hitRiseTimeInIndex = (T)hitPeakTimeIndex;
                    T hitFallTimeInIndex = (T)hitPeakTimeIndex;
                    T hitPeakTimeInNanoSec = FindPeakValue(hitPeakTimeIndex, hitPeakValue);
                    T hitWidthInTicks = FindHitWidth(hitPeakTimeIndex, hitPeakValue);
                    bool containedFromRisingEdge = BackwardFindingOfHitStart(hitPeakTimeIndex, hitPeakValue, hitStartTimeIndex, hitRiseTimeInIndex);
                    bool containedFromFallingEdge = ForwardFindingOfHitFallTime(hitPeakTimeIndex, hitFallTimeInIndex);

                    hitPeakValueVector.push_back(hitPeakValue);
                    hitPeakTimeIndexVector.push_back(hitPeakTimeIndex);
                    hitPeakTimeInNanoSecVector.push_back(hitPeakTimeInNanoSec);
                    hitStartTimeIndexVector.push_back(hitStartTimeIndex);
                    hitRiseTimeInIndexVector.push_back(hitRiseTimeInIndex);
                    hitFallTimeInIndexVector.push_back(hitFallTimeInIndex);
                    hitWidthValueVector.push_back(hitWidthInTicks);
                    isHitContainedVector.push_back(containedFromRisingEdge && containedFromFallingEdge);
                    break;
                }
            }
        }
    }

    for (size_t i = 0; i < hitStartTimeIndexVector.size(); i++) {
        hit_t<double> aHit;
        aHit.Timestamp                       = _Timestamp;
        aHit.DigitizerChannel                = _ChannelNumber;
        aHit.TStartInNanoSec                 = hitStartTimeIndexVector.at(i) * _CFDParamSet[kTimeSamplingInterval];
        aHit.TPeakInNanoSec                  = hitPeakTimeInNanoSecVector.at(i) * _CFDParamSet[kTimeSamplingInterval];
        aHit.AmplitudeInMiliVolt             = hitPeakValueVector.at(i) * (_CFDParamSet[kADCDynamicRange] / (TMath::Power(2, _CFDParamSet[kADCNBits]) - 1));
        aHit.AmplitudeInADC                  = hitPeakValueVector.at(i);
        aHit.WidthInNanoSec                  = hitWidthValueVector.at(i) * _CFDParamSet[kTimeSamplingInterval];
        aHit.RiseTimeInNanoSec               = hitRiseTimeInIndexVector.at(i) * _CFDParamSet[kTimeSamplingInterval];
        aHit.FallTimeInNanoSec                   = hitFallTimeInIndexVector.at(i) * _CFDParamSet[kTimeSamplingInterval];
        aHit.IsContained                     = isHitContainedVector.at(i);
        aHit.IntegratedChargeInADCTimeTicks  = IntegrateWaveformInADC(hitPeakTimeIndexVector.at(i) - hitRiseTimeInIndexVector.at(i),
                                                                      hitPeakTimeIndexVector.at(i) + hitFallTimeInIndexVector.at(i));
        aHit.IntegratedChargeInADCNanoSec    = aHit.IntegratedChargeInADCTimeTicks * _CFDParamSet[kTimeSamplingInterval];

        _HitCollection[hitStartTimeIndexVector.at(i)] = aHit;
    }
}

template<class T>

T beamlinereco::CFDHitFinder< T >::FindHitWidth ( size_t  hitPeakTimeIndex, T  hitPeakValueADC  )

privatevirtual

Definition at line 539 of file CFDHitFinderAlg.h.

References MECModelEnuComparisons::i, beamlinereco::kNSamplingPoints, T, getGoodRuns4SAM::t1, t2, y1, and submit_syst::y2.

                                                                                      {
    // Calculate FWHM

    // Find half-max
    T halfMax = _PedestalInADC - (_PedestalInADC - hitPeakValueADC)/2.;

    T halfMaxLeft = 9993;  // left point of half-max in ticks
    T halfMaxRight = 9993; // right point of half-max in ticks

    // Scan backwards
    for (size_t i=hitPeakTimeIndex; i>0; i--) {
        if (_WaveformADCNanosec.at(i) == halfMax) { // exactly equal
            halfMaxLeft = (T)i;
            break;
        } else if (_WaveformADCNanosec.at(i) > halfMax) {
            // Simple two-point linear interpolation
            T y1 = _WaveformADCNanosec.at(i);
            T y2 = _WaveformADCNanosec.at(i + 1);
            size_t t1 = i;
            size_t t2 = i + 1;
            T slope = (y2 - y1) / ((T)t2 - (T)t1);
            halfMaxLeft = (halfMax - y1) / slope + (T)t1;
            break;
        }
    }

    // Scan forwards
    for (size_t i=hitPeakTimeIndex; i<_CFDParamSet[kNSamplingPoints]-1; i++) {
        if (_WaveformADCNanosec.at(i) == halfMax) { // exactly equal
            halfMaxRight = (T)i;
            break;
        } else if (_WaveformADCNanosec.at(i) > halfMax) {
            // Simple two-point linear interpolation
            T y1 = _WaveformADCNanosec.at(i - 1);
            T y2 = _WaveformADCNanosec.at(i);
            size_t t1 = i - 1;
            size_t t2 = i;
            T slope = (y2 - y1) / ((T)t2 - (T)t1);
            halfMaxRight = (halfMax - y1) / slope + (T)t1;
            break;
        }
    }

    T widthInTicks = halfMaxRight - halfMaxLeft;
    return widthInTicks;
}

template<class T>

T beamlinereco::CFDHitFinder< T >::FindPeakValue ( size_t  hitPeakTimeIndex, T &  hitPeakValue  )

privatevirtual

Definition at line 394 of file CFDHitFinderAlg.h.

References PandAna.Demos.pi0_spectra::p0, plot_validation_datamc::p1, plot_validation_datamc::p2, make_associated_cosmic_defs::p3, std::sqrt(), T, submit_syst::x, and submit_syst::y.

                                                                                     {
    // Fit a third-order polynomial
    // Translate x to near zero to help out the fitter
    T x[7] = { -3, -2, -1, 0, 1, 2, 3 };
    T y[7] = {
        _WaveformADCNanosec.at(hitPeakTimeIndex - 3),
        _WaveformADCNanosec.at(hitPeakTimeIndex - 2),
        _WaveformADCNanosec.at(hitPeakTimeIndex - 1),
        _WaveformADCNanosec.at(hitPeakTimeIndex),
        _WaveformADCNanosec.at(hitPeakTimeIndex + 1),
        _WaveformADCNanosec.at(hitPeakTimeIndex + 2),
        _WaveformADCNanosec.at(hitPeakTimeIndex + 3)
    };

    TGraph *gfit = new TGraph(7, x, y);

    // Make sure fitter is Minuit2 -- works better than Minut here
    TVirtualFitter::SetDefaultFitter("Minuit2");
    TFitResultPtr frp = gfit->Fit("pol3", "SFNQ");
    T p0 = frp->Parameter(0);
    T p1 = frp->Parameter(1);
    T p2 = frp->Parameter(2);
    T p3 = frp->Parameter(3);

    T disc = p2*p2 - 3.*p1*p3;
    if (disc < 0) { // Discriminant negative => no extrema
        hitPeakValue = _WaveformADCNanosec.at(hitPeakTimeIndex);
        return hitPeakTimeIndex;
    }

    // Positive root corresponds to the minimum for a third-order polynomial
    T hitPeakTime = (-p2 + sqrt(disc))/(3.*p3);
    hitPeakValue = p0 + p1*hitPeakTime + p2*hitPeakTime*hitPeakTime + p3*hitPeakTime*hitPeakTime*hitPeakTime;
    hitPeakTime += (T)hitPeakTimeIndex; // translate back to actual time
    return hitPeakTime;
}

template<class T >

void beamlinereco::CFDHitFinder< T >::FindPedestal ( )

privatevirtual

Definition at line 242 of file CFDHitFinderAlg.h.

References ana::Sqrt(), and T.

                                               {
/* 2nd way
    TH1D* pedHist = new TH1D("", "", 10000, 0, 5000);
    for (size_t i = 700; i < 1024; i++) {
        pedHist->Fill(_NonfilterWaveformADCNanosec.at(i));
    }

    _PedestalInADC   = (T) pedHist->GetMean();
    _NoiseSigmaInADC = (T) pedHist->GetStdDev();

    delete (pedHist);
*/

/* 3rd way*/
    T NoiseBand = 1E9;
    T OldPedestal = 1E9;
    T Pedestal  = 0;

    while (TMath::Abs(Pedestal - OldPedestal) >= 2) {
        T average = 0;
        T stdDev = 0;
        size_t countSamples = 0;
        for (typename std::vector<T>::iterator itr = _NonfilterWaveformADCNanosec.begin(); itr != _NonfilterWaveformADCNanosec.end(); itr++) {
            if ((*itr <= Pedestal + NoiseBand) && (*itr >= Pedestal - NoiseBand)) {
                average += *itr;
                countSamples++;
            }
        }
        average = average / (T)countSamples;
        for (typename std::vector<T>::iterator itr = _NonfilterWaveformADCNanosec.begin(); itr != _NonfilterWaveformADCNanosec.end(); itr++) {
            if ((*itr <= Pedestal + NoiseBand) && (*itr >= Pedestal - NoiseBand)) {
                stdDev += TMath::Power((*itr) - average, 2);
            }
        }

        stdDev = TMath::Sqrt(stdDev / (T)countSamples);
        OldPedestal = Pedestal;
        Pedestal = average;
        NoiseBand = stdDev;
    }

    _PedestalInADC = Pedestal;
    _NoiseSigmaInADC = NoiseBand;
/**/
}

template<class T >

void beamlinereco::CFDHitFinder< T >::FindRawHitLogic ( )

privatevirtual

Definition at line 289 of file CFDHitFinderAlg.h.

References MECModelEnuComparisons::i, calib::j, beamlinereco::kConsecutiveHitSeperationDurationInTicks, beamlinereco::kRawHitFinderThresholdInNoiseSigma, beamlinereco::kShortRawHitIgnoringDurationInTicks, and T.

                                                  {
    T rawHitThreshold = (_PedestalInADC) - ((_NoiseSigmaInADC) * _CFDParamSet[kRawHitFinderThresholdInNoiseSigma]);
    for (size_t i = 0; i < _WaveformADCNanosec.size(); i++) {
        _RawHitLogicNanosec.push_back(_WaveformADCNanosec.at(i) < rawHitThreshold);
    }

    for (size_t i = 1; i < _RawHitLogicNanosec.size() - 1; i++) {
        if ((not _RawHitLogicNanosec.at(i - 1)) and _RawHitLogicNanosec.at(i)) {
            size_t countDuration = 0;
            for (size_t j = i; j < _RawHitLogicNanosec.size() - 1; j++) {
                if (_RawHitLogicNanosec.at(j)) countDuration++;
                if (not(_RawHitLogicNanosec.at(j + 1))) break;
            }

            if (countDuration < _CFDParamSet[kShortRawHitIgnoringDurationInTicks]) {
                for (size_t j = i; j < i + countDuration; j++) {
                    _RawHitLogicNanosec.at(j) = false;
                }
            }
        }
    }

    for (size_t i = 1; i < _RawHitLogicNanosec.size() - 1; i++) {
        if (_RawHitLogicNanosec.at(i - 1) && not(_RawHitLogicNanosec.at(i))) {
            size_t countDuration = 0;
            for (size_t j = i; j < _RawHitLogicNanosec.size() - 1; j++) {
                if (not(_RawHitLogicNanosec.at(j))) countDuration++;
                if (_RawHitLogicNanosec.at(j + 1)) break;
            }

            if (countDuration < _CFDParamSet[kConsecutiveHitSeperationDurationInTicks]) {
                for (size_t j = i; j < i + countDuration; j++) {
                    _RawHitLogicNanosec.at(j) = true;
                }
            }
        }
    }

}

template<class T>

bool beamlinereco::CFDHitFinder< T >::ForwardFindingOfHitFallTime ( size_t  hitPeakTimeIndex, T &  hitFallTimeInIndex  )

privatevirtual

Definition at line 523 of file CFDHitFinderAlg.h.

References beamlinereco::kNSamplingPoints, and T.

                                                                                                            {
    bool isThisHitContainedFromTheFallingEdge = true;
    size_t tmp_index = hitPeakTimeIndex;
    while (_WaveformADCNanosec.at(tmp_index) < _PedestalInADC) {
        if (tmp_index == _CFDParamSet[kNSamplingPoints] - 2 && _WaveformADCNanosec.at(tmp_index + 1) < _PedestalInADC) {
            isThisHitContainedFromTheFallingEdge = false;
            break;
        }
        tmp_index = tmp_index + 1;
    }

    hitFallTimeInIndex = (T)tmp_index - (T)hitPeakTimeIndex;
    return isThisHitContainedFromTheFallingEdge;
}

template<class T>

const T& beamlinereco::CFDHitFinder< T >::GetDiscriminationThreshold ( ) const

inline

Definition at line 100 of file CFDHitFinderAlg.h.

{ return _DiscriminationThresholdInADC; }

template<class T>

const std::string beamlinereco::CFDHitFinder< T >::GetDiscriminationType ( ) const

inline

Definition at line 94 of file CFDHitFinderAlg.h.

{ return _DiscriminationType; }

template<class T>

const std::string beamlinereco::CFDHitFinder< T >::GetFilterType ( ) const

inline

Definition at line 96 of file CFDHitFinderAlg.h.

{ return _FilterType; }

template<class T>

const std::map<T, hit_t<T> >& beamlinereco::CFDHitFinder< T >::GetHitCollection ( ) const

inline

Definition at line 97 of file CFDHitFinderAlg.h.

Referenced by beamlinereco::CFDHitReco::produce().

{ return _HitCollection; }

template<class T>

const std::string beamlinereco::CFDHitFinder< T >::GetInterpolationType ( ) const

inline

Definition at line 95 of file CFDHitFinderAlg.h.

{ return _InterpolationType; }

template<class T>

const T& beamlinereco::CFDHitFinder< T >::GetNoiseSigma ( ) const

inline

Definition at line 99 of file CFDHitFinderAlg.h.

{ return _NoiseSigmaInADC; }

template<class T>

const T& beamlinereco::CFDHitFinder< T >::GetPedestal ( ) const

inline

Definition at line 98 of file CFDHitFinderAlg.h.

Referenced by beamlinereco::CFDHitReco::produce().

{ return _PedestalInADC; }

template<class T >

void beamlinereco::CFDHitFinder< T >::Go ( )

virtual

Definition at line 587 of file CFDHitFinderAlg.h.

Referenced by beamlinereco::CFDHitReco::produce().

                                     {
    if (_WaveformADCNanosec.size() == 0) return;

    FindPedestal(); // Get pedestal and noise band sigma
    FindRawHitLogic(); // From pedestal and noise band, get rawHitThreshold and RawHitLogic vector
    FindCFDHits(); // From RawHitLogic vector, get hitPeakValue and hitPeakTimeIndex. From hitPeakValue and pedestal, get cfdHitThresholdValue. Go back from hitPeakTimeIndex to the point when
}

template<class T >

T beamlinereco::CFDHitFinder< T >::IntegrateWaveformInADC ( size_t  hitStartTimeIndex, size_t  hitEndTimeIndex  )

privatevirtual

Definition at line 432 of file CFDHitFinderAlg.h.

References MECModelEnuComparisons::i, beamlinereco::kIntergratedWindowFixed, beamlinereco::kIntergratedWindowLowerLimitIndex, beamlinereco::kIntergratedWindowUpperLimitIndex, beamlinereco::kNSamplingPoints, and T.

                                                                                                      {

    T integration = 0;

    if (!_CFDParamSet[kIntergratedWindowFixed]) {
        hitStartTimeIndex = hitStartTimeIndex < 0 ? 0 : hitStartTimeIndex;
        hitEndTimeIndex   = hitEndTimeIndex > (_CFDParamSet[kNSamplingPoints] - 1) ? _CFDParamSet[kNSamplingPoints] - 1 : hitEndTimeIndex;
        for (size_t i = hitStartTimeIndex; i <= hitEndTimeIndex; i++) {
            integration = integration + (_WaveformADCNanosec.at(i) - _PedestalInADC);
        }
    } else {
        for (size_t i = _CFDParamSet[kIntergratedWindowLowerLimitIndex]; i <= _CFDParamSet[kIntergratedWindowUpperLimitIndex]; i++) {
            integration = integration + (_WaveformADCNanosec.at(i) - _PedestalInADC);
        }
    }

    return integration;
}

template<class T>

virtual void beamlinereco::CFDHitFinder< T >::Reset ( )

inlineprivatevirtual

Definition at line 136 of file CFDHitFinderAlg.h.

                                    {
            _HitCollection.clear();
            _RawHitLogicNanosec.clear();
            _PedestalInADC = 0;
            _NoiseSigmaInADC = 0;
            _DiscriminationThresholdInADC = 0;
            _WaveformADCNanosec.clear();
            _ChannelNumber = 993;
            _Timestamp = 0;
            _NonfilterWaveformADCNanosec.clear();
        }

template<class T >

void beamlinereco::CFDHitFinder< T >::SetChannel ( unsigned int  channelNo )

private

Definition at line 217 of file CFDHitFinderAlg.h.

                                                                   {
    _ChannelNumber = channelNo;
}

template<class T >

void beamlinereco::CFDHitFinder< T >::SetDiscriminationType

(

const std::string

type

)

Definition at line 232 of file CFDHitFinderAlg.h.

                                                                            {
        _DiscriminationType = type;
}

template<class T >

void beamlinereco::CFDHitFinder< T >::SetFilterType

(

const std::string

type

)

Definition at line 237 of file CFDHitFinderAlg.h.

                                                                    {
  _FilterType = type;
}

template<class T >

void beamlinereco::CFDHitFinder< T >::SetInterpolationType

(

const std::string

type

)

Definition at line 227 of file CFDHitFinderAlg.h.

                                                                           {
        _InterpolationType = type;
}

template<class T>

void beamlinereco::CFDHitFinder< T >::SetParam ( CFDParams  param, T  value  )

private

Definition at line 205 of file CFDHitFinderAlg.h.

References ParseFluxesFile::param.

                                                                   {
    _CFDParamSet[param] = value;
}

template<class T >

void beamlinereco::CFDHitFinder< T >::SetParams

(

const std::map< CFDParams, double > &

paramSet

)

Definition at line 210 of file CFDHitFinderAlg.h.

Referenced by beamlinereco::CFDHitReco::reconfigure().

                                                                                     {
    for (auto & itr : paramSet) {
        SetParam(itr.first, itr.second);
    }
}

template<class T >

void beamlinereco::CFDHitFinder< T >::SetTimestamp ( uint32_t  timestamp )

private

Definition at line 222 of file CFDHitFinderAlg.h.

                                                          {
    _Timestamp = ts;
}

template<class T >

void beamlinereco::CFDHitFinder< T >::SetWaveform	(	std::vector< uint16_t > &	waveform,
		unsigned int	channelNo,
		uint32_t	timestamp
	)

Definition at line 154 of file CFDHitFinderAlg.h.

References art::errors::Configuration, allTimeWatchdog::endl, beamlinereco::Kalman::GetFilteredValue(), MECModelEnuComparisons::i, K, beamlinereco::kGSFilterDegree, beamlinereco::kGSFilterWindow, beamlinereco::kIsWaveformNegativePolarity, beamlinereco::kKalmanFilterGain, beamlinereco::kKalmanFilterMeasurementNoiseCovariance, beamlinereco::kKalmanFilterProcessNoiseCovariance, beamlinereco::kNSamplingPoints, P, r(), slidt::Reset(), beamlinereco::SGSmoothing::Smooth(), and T.

Referenced by beamlinereco::CFDHitReco::produce().

                                                                                                                       {
  Reset();

  SetChannel(channelNo);
  SetTimestamp(timestamp);

  double polaritySignFactor = _CFDParamSet[kIsWaveformNegativePolarity] ? 1. : -1.;

  T* nonfilterwaveform = (T*) malloc(_CFDParamSet[kNSamplingPoints] * sizeof(T));
  T* filterwaveform    = (T*) malloc(_CFDParamSet[kNSamplingPoints] * sizeof(T));

  for (auto itr = waveform.begin(); itr != waveform.end(); itr++) {
    _NonfilterWaveformADCNanosec.push_back((T)(*itr) * polaritySignFactor);
    nonfilterwaveform[itr - waveform.begin()] = (T)(*itr) * polaritySignFactor;
  }

  if (_FilterType == "GS" || _FilterType == "gs") { // Savitzky-Golay smoothing
    SGSmoothing::Smooth(_CFDParamSet[kNSamplingPoints],
                        nonfilterwaveform,
                        filterwaveform,
                        _CFDParamSet[kGSFilterWindow],
                        _CFDParamSet[kGSFilterDegree]);
  } else if (_FilterType == "Kalman" || _FilterType == "kalman" || _FilterType == "KALMAN") { // Kalman filter
    float q = _CFDParamSet[kKalmanFilterProcessNoiseCovariance];
    float r = _CFDParamSet[kKalmanFilterMeasurementNoiseCovariance];
    float P = .05;
    float K = _CFDParamSet[kKalmanFilterGain];
    Kalman filter = Kalman(q,r,P,K);
    for (unsigned int i=0; i<_CFDParamSet[kNSamplingPoints]; i++) {
      filterwaveform[i] = filter.GetFilteredValue(nonfilterwaveform[i]);
    }
  } else if (_FilterType == "" || _FilterType == "None" || _FilterType == "none" || _FilterType == "NONE") {
    // no smoothing
    for (auto itr = waveform.begin(); itr != waveform.end(); itr++) {
      filterwaveform[itr - waveform.begin()] = (T)(*itr) * polaritySignFactor;
    }
  } else {
    // unknown -- throw exception
    throw art::Exception(art::errors::Configuration)
      << "Unknown raw hit finder filter type: " << _FilterType << "." << std::endl;
  }

  for (size_t i = 0; i < _CFDParamSet[kNSamplingPoints]; i++) {
    _WaveformADCNanosec.push_back(*(filterwaveform + i));
  }

  free(nonfilterwaveform);
  free(filterwaveform);
}

Member Data Documentation

template<class T>

std::map<CFDParams, T> beamlinereco::CFDHitFinder< T >::_CFDParamSet

private

Definition at line 123 of file CFDHitFinderAlg.h.

template<class T>

unsigned int beamlinereco::CFDHitFinder< T >::_ChannelNumber

private

Definition at line 111 of file CFDHitFinderAlg.h.

template<class T>

T beamlinereco::CFDHitFinder< T >::_DiscriminationThresholdInADC

private

Definition at line 121 of file CFDHitFinderAlg.h.

template<class T>

std::string beamlinereco::CFDHitFinder< T >::_DiscriminationType

private

Definition at line 113 of file CFDHitFinderAlg.h.

template<class T>

std::string beamlinereco::CFDHitFinder< T >::_FilterType

private

Definition at line 115 of file CFDHitFinderAlg.h.

template<class T>

std::map<T, hit_t<T> > beamlinereco::CFDHitFinder< T >::_HitCollection

private

Definition at line 124 of file CFDHitFinderAlg.h.

template<class T>

std::string beamlinereco::CFDHitFinder< T >::_InterpolationType

private

Definition at line 114 of file CFDHitFinderAlg.h.

template<class T>

T beamlinereco::CFDHitFinder< T >::_NoiseSigmaInADC

private

Definition at line 120 of file CFDHitFinderAlg.h.

template<class T>

std::vector<T> beamlinereco::CFDHitFinder< T >::_NonfilterWaveformADCNanosec

private

Definition at line 117 of file CFDHitFinderAlg.h.

template<class T>

T beamlinereco::CFDHitFinder< T >::_PedestalInADC

private

Definition at line 119 of file CFDHitFinderAlg.h.

template<class T>

std::vector<bool> beamlinereco::CFDHitFinder< T >::_RawHitLogicNanosec

private

Definition at line 118 of file CFDHitFinderAlg.h.

template<class T>

uint32_t beamlinereco::CFDHitFinder< T >::_Timestamp

private

Definition at line 112 of file CFDHitFinderAlg.h.

template<class T>

std::vector<T> beamlinereco::CFDHitFinder< T >::_WaveformADCNanosec

private

Definition at line 116 of file CFDHitFinderAlg.h.

---

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

• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-01-21/BeamlineReco/CFDHitFinderAlg.h