hit::RFFHitFinderAlg Class Reference

#include <RFFHitFinderAlg.h>

Public Member Functions
	RFFHitFinderAlg (fhicl::ParameterSet const &)
void	SetFitterParamsVectors (geo::Geometry const &)
void	Run (std::vector< recob::Wire > const &, std::vector< recob::Hit > &, geo::Geometry const &)

Private Member Functions
void	SetFitterParams (unsigned int)
void	EmplaceHit (std::vector< recob::Hit > &, recob::Wire const &, float const &, raw::TDCtick_t const &, raw::TDCtick_t const &, geo::SigType_t const &, geo::WireID const &)

Private Attributes
const float	SQRT_TWO_PI = 2.506628
std::vector< float >	fMatchThresholdVec
std::vector< unsigned int >	fMergeMultiplicityVec
std::vector< float >	fAmpThresholdVec
RFFHitFitter	fFitter

Detailed Description

Definition at line 31 of file RFFHitFinderAlg.h.

Constructor & Destructor Documentation

hit::RFFHitFinderAlg::RFFHitFinderAlg

(

fhicl::ParameterSet const &

p

)

Title: RFFHitFinderAlg Class Author: Wes Ketchum (wketc.nosp@m.hum@.nosp@m.lanl..nosp@m.gov)

Description: Class that runs the RFF HitFinder. Implements an RFFHitFitter, and takes the result and stores it in recob::Hit objects.

Input: recob::Wire Output: recob::Hit

Definition at line 21 of file RFFHitFinderAlg.cxx.

{
    fMatchThresholdVec = p.get< std::vector<float> >("MeanMatchThreshold");
    fMergeMultiplicityVec = p.get< std::vector<unsigned int> >("MinMergeMultiplicity");
    fAmpThresholdVec = p.get< std::vector<float> >("AmplitudeThreshold",std::vector<float>(1,0.0));
}

Member Function Documentation

void hit::RFFHitFinderAlg::EmplaceHit ( std::vector< recob::Hit > &  hitVector, recob::Wire const &  wire, float const &  summedADCTotal, raw::TDCtick_t const &  startTick, raw::TDCtick_t const &  endTick, geo::SigType_t const &  sigtype, geo::WireID const &  wireID  )

private

Definition at line 87 of file RFFHitFinderAlg.cxx.

{

    float totalArea = 0.0;
    std::vector<float> areaVector(fFitter.NHits());
    std::vector<float> areaErrorVector(fFitter.NHits());
    std::vector<float> areaFracVector(fFitter.NHits());

    for(size_t ihit=0; ihit < fFitter.NHits(); ihit++){
        areaVector[ihit] = fFitter.AmplitudeVector()[ihit]*fFitter.SigmaVector()[ihit]*SQRT_TWO_PI;
        areaErrorVector[ihit] =
            SQRT_TWO_PI*std::sqrt(fFitter.AmplitudeVector()[ihit]*fFitter.SigmaErrorVector()[ihit]*fFitter.AmplitudeVector()[ihit]*fFitter.SigmaErrorVector()[ihit] +
                            fFitter.AmplitudeErrorVector()[ihit]*fFitter.SigmaVector()[ihit]*fFitter.AmplitudeErrorVector()[ihit]*fFitter.SigmaVector()[ihit]);
        totalArea += areaVector[ihit];
    }

    for(size_t ihit=0; ihit < fFitter.NHits(); ihit++)
    {
        areaFracVector[ihit] = areaVector[ihit]/totalArea;

        hitVector.emplace_back(wire.Channel(),
                               startTick,
                               endTick,
                               fFitter.MeanVector()[ihit]+(float)startTick,
                               fFitter.MeanErrorVector()[ihit],
                               fFitter.SigmaVector()[ihit],
                               fFitter.AmplitudeVector()[ihit],
                               fFitter.AmplitudeErrorVector()[ihit],
                               summedADCTotal*areaFracVector[ihit],
                               areaVector[ihit],
                               areaErrorVector[ihit],
                               fFitter.NHits(),
                               ihit,
                               -999.,
                               -999,
                               wire.View(),
                               sigtype,
                               wireID);
    }

}

void hit::RFFHitFinderAlg::Run	(	std::vector< recob::Wire > const &	wireVector,
		std::vector< recob::Hit > &	hitVector,
		geo::Geometry const &	geo
	)

Definition at line 60 of file RFFHitFinderAlg.cxx.

{
    hitVector.reserve(wireVector.size());
    for(auto const& wire : wireVector)
    {
        geo::SigType_t const& sigtype = geo.SignalType(wire.Channel());
        geo::WireID const& wireID = geo.ChannelToWire(wire.Channel()).at(0);

        SetFitterParams(wire.View());

        for(auto const& roi : wire.SignalROI().get_ranges())
        {
            fFitter.RunFitter(roi.data());

            const float summedADCTotal = std::accumulate(roi.data().begin(),roi.data().end(),0.0);
            const raw::TDCtick_t startTick = roi.begin_index();
            const raw::TDCtick_t endTick = roi.begin_index()+roi.size();

            EmplaceHit(hitVector,wire,summedADCTotal,startTick,endTick,sigtype,wireID);
        }//end loop over ROIs on wire

    }//end loop over wires

}

void hit::RFFHitFinderAlg::SetFitterParams ( unsigned int  p )

private

Definition at line 55 of file RFFHitFinderAlg.cxx.

{
    fFitter.SetFitterParams(fMatchThresholdVec[p],fMergeMultiplicityVec[p],fAmpThresholdVec[p]);
}

void hit::RFFHitFinderAlg::SetFitterParamsVectors

(

geo::Geometry const &

geo

)

Definition at line 28 of file RFFHitFinderAlg.cxx.

{
    const unsigned int n_planes = geo.Nplanes();

    //If size zero, throw. If size one, assume same for all planes.
    //If size > 1 but < n_planes, throw. If size = n_plane, good.

    if(fMatchThresholdVec.size()==0 ||
       fMergeMultiplicityVec.size()==0 ||
       fAmpThresholdVec.size()==0)
        throw std::runtime_error("Error in RFFHitFinderAlg: Configured with zero planes.");

    if( (fMatchThresholdVec.size()>1 && fMatchThresholdVec.size()<n_planes) ||
        (fMergeMultiplicityVec.size()>1 && fMergeMultiplicityVec.size()<n_planes) ||
        (fAmpThresholdVec.size()>1 && fAmpThresholdVec.size()<n_planes) )
        throw std::runtime_error("Error in RFFHitFinderAlg: Configured with incorrect n_planes.");

    if(fMatchThresholdVec.size()==1)
        fMatchThresholdVec.resize(n_planes,fMatchThresholdVec[0]);

    if(fMergeMultiplicityVec.size()==1)
        fMergeMultiplicityVec.resize(n_planes,fMergeMultiplicityVec[0]);

    if(fAmpThresholdVec.size()==1)
        fAmpThresholdVec.resize(n_planes,fAmpThresholdVec[0]);
}

Member Data Documentation

std::vector<float> hit::RFFHitFinderAlg::fAmpThresholdVec

private

Definition at line 47 of file RFFHitFinderAlg.h.

RFFHitFitter hit::RFFHitFinderAlg::fFitter

private

Definition at line 58 of file RFFHitFinderAlg.h.

std::vector<float> hit::RFFHitFinderAlg::fMatchThresholdVec

private

Definition at line 45 of file RFFHitFinderAlg.h.

std::vector<unsigned int> hit::RFFHitFinderAlg::fMergeMultiplicityVec

private

Definition at line 46 of file RFFHitFinderAlg.h.

const float hit::RFFHitFinderAlg::SQRT_TWO_PI = 2.506628

private

Definition at line 33 of file RFFHitFinderAlg.h.

---

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

• larreco/larreco/HitFinder/RFFHitFinderAlg.h
• larreco/larreco/HitFinder/RFFHitFinderAlg.cxx