icarus_tool::Filter Class Reference

Inheritance diagram for icarus_tool::Filter:

Public Member Functions
	Filter (const fhicl::ParameterSet &pset)
	~Filter ()
void	configure (const fhicl::ParameterSet &pset) override
void	setResponse (size_t numBins, double correct3D, double timeScaleFctr) override
void	outputHistograms (art::TFileDirectory &) const override
size_t	getPlane () const override
const icarusutil::FrequencyVec &	getResponseVec () const override

Private Attributes
size_t	fPlane
std::vector< double >	fParameters
std::string	fFunctionString
double	fFilterWidthCorrectionFactor
icarusutil::FrequencyVec	fFilterVec

Additional Inherited Members
Private Member Functions inherited from icarus_tool::IFilter
virtual	~IFilter () noexcept=default

Detailed Description

Definition at line 23 of file Filter_tool.cc.

Constructor & Destructor Documentation

icarus_tool::Filter::Filter ( const fhicl::ParameterSet &  pset )

explicit

Definition at line 50 of file Filter_tool.cc.

{
    configure(pset);
}

icarus_tool::Filter::~Filter

(

)

Definition at line 55 of file Filter_tool.cc.

{
    return;
}

Member Function Documentation

void icarus_tool::Filter::configure ( const fhicl::ParameterSet &  pset )

overridevirtual

Implements icarus_tool::IFilter.

Definition at line 60 of file Filter_tool.cc.

{
    // Start by recovering the parameters
    fPlane                       = pset.get<size_t>("Plane");
    fParameters                  = pset.get<std::vector<double>>("FilterParametersVec");
    fFunctionString              = pset.get<std::string>("FilterFunction");
    fFilterWidthCorrectionFactor = pset.get<double>("FilterWidthCorrectionFactor");
    
    return;
}

size_t icarus_tool::Filter::getPlane ( ) const

inlineoverridevirtual

Implements icarus_tool::IFilter.

Definition at line 34 of file Filter_tool.cc.

{return fPlane;}

const icarusutil::FrequencyVec& icarus_tool::Filter::getResponseVec ( ) const

inlineoverridevirtual

Implements icarus_tool::IFilter.

Definition at line 35 of file Filter_tool.cc.

{return fFilterVec;}

void icarus_tool::Filter::outputHistograms ( art::TFileDirectory &  histDir ) const

overridevirtual

Implements icarus_tool::IFilter.

Definition at line 124 of file Filter_tool.cc.

{
    // It is assumed that the input TFileDirectory has been set up to group histograms into a common
    // folder at the calling routine's level. Here we create one more level of indirection to keep
    // histograms made by this tool separate.
    std::string dirName = "FilterPlane_" + std::to_string(fPlane);
    
    art::TFileDirectory dir = histDir.mkdir(dirName.c_str());
    
    auto const clockData = art::ServiceHandle<detinfo::DetectorClocksService const>()->DataForJob();
    double      samplingRate = 1.e-3 * sampling_rate(clockData); // Note sampling rate is in ns, convert to us
    double      numBins      = fFilterVec.size();
    double      maxFreq      = 1.e3 / samplingRate;      // Max frequency in MHz
    double      minFreq      = maxFreq / numBins;
    std::string histName     = "FilterPlane_" + std::to_string(fPlane);
    TProfile*   hist         = dir.make<TProfile>(histName.c_str(), "Filter;Frequency(kHz)", numBins/2, minFreq, 0.5*maxFreq);
    
    for(int bin = 0; bin < numBins/2; bin++)
    {
        double freq = bin * minFreq;
        
        hist->Fill(freq, std::abs(fFilterVec.at(bin)), 1.);
    }
    
    
    return;
}

void icarus_tool::Filter::setResponse ( size_t  numBins, double  correct3D, double  timeScaleFctr  )

overridevirtual

Implements icarus_tool::IFilter.

Definition at line 71 of file Filter_tool.cc.

{
    // Note that here we are working in frequency space, not in the time domain...
    // FIXME: Assumes that the job-level clock data is sufficient.
    auto const clockData = art::ServiceHandle<detinfo::DetectorClocksService const>()->DataForJob();
    double      samplingRate = 1.e-3 * sampling_rate(clockData); // Note sampling rate is in ns, convert to us
    double      maxFreq      = 1.e3 / (2. * samplingRate);      // highest frequency in cycles/us (MHz)
    double      freqRes      = maxFreq / double(numBins/2);     // frequency resolution in cycles/us
    
    std::string funcName = "tempFilter";
    TF1 function(funcName.c_str(),fFunctionString.c_str());

    // Set the range on the function, probably not really necessary
    function.SetRange(0, maxFreq);

    // now to scale the filter function!
    // only scale params 1,2 &3
    double timeFactor = 1. / (timeScaleFctr * correct3D * fFilterWidthCorrectionFactor);
    size_t paramIdx(0);

    for(const auto& parameter : fParameters) function.SetParameter(paramIdx++, timeFactor * parameter);
    
    // Don't assume that the filter vec has not already been initialized...
    // Note we are setting up the FFT to be the "full" so add a bin for folding over
    fFilterVec.resize(numBins,icarusutil::ComplexVal(0.,0.));
    
    // Keep track of the peak value
    double peakVal(std::numeric_limits<double>::min());

    size_t nyquistBin = numBins/2 + 1;

    // Now ready to set the response vector
    for(size_t bin = 0; bin < nyquistBin; bin++)
    {
        // This takes a sampling rate in ns -> gives a frequency in cycles/us
        double freq = bin * freqRes;
        double f    = function.Eval(freq);

        peakVal = std::max(peakVal, f);
        
        fFilterVec[bin] = icarusutil::ComplexVal(f, 0.);
    }

    // Since we are returning the "full" FFT... folder over the first half
    for(size_t bin = nyquistBin; bin < numBins; bin++)
        fFilterVec[bin] = std::conj(fFilterVec[numBins - bin]);
    
    // "Normalize" to peak value
    for(auto& filterValue : fFilterVec) filterValue = filterValue / peakVal;

    return;
}

Member Data Documentation

icarusutil::FrequencyVec icarus_tool::Filter::fFilterVec

private

Definition at line 45 of file Filter_tool.cc.

double icarus_tool::Filter::fFilterWidthCorrectionFactor

private

Definition at line 42 of file Filter_tool.cc.

std::string icarus_tool::Filter::fFunctionString

private

Definition at line 41 of file Filter_tool.cc.

std::vector<double> icarus_tool::Filter::fParameters

private

Definition at line 40 of file Filter_tool.cc.

size_t icarus_tool::Filter::fPlane

private

Definition at line 39 of file Filter_tool.cc.

---

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

• Filter_tool.cc