pmtana::RiseTimeGaussFit Class Reference

Inheritance diagram for pmtana::RiseTimeGaussFit:

Classes
struct	Config

Public Member Functions
	RiseTimeGaussFit (art::ToolConfigTable< Config > const &config)
double	RiseTime (const pmtana::Waveform_t &wf_pulse, const pmtana::PedestalMean_t &ped_pulse, bool _positive) const override
std::size_t	findFirstMax (const std::vector< double > &arr, double threshold) const

Private Attributes
double	fMinAmp
double	fInitSigma
double	fTolerance
int	fNbins

Detailed Description

Definition at line 24 of file RiseTimeGaussFit_tool.cc.

Constructor & Destructor Documentation

pmtana::RiseTimeGaussFit::RiseTimeGaussFit ( art::ToolConfigTable< Config > const &  config )

explicit

Definition at line 53 of file RiseTimeGaussFit_tool.cc.

References fInitSigma, fNbins, and fTolerance.

    : fMinAmp{config().MinAmp()}
    , fInitSigma{config().InitSigma()}
    , fTolerance{config().Tolerance()}
    , fNbins{config().Nbins()}
  {}

Member Function Documentation

std::size_t pmtana::RiseTimeGaussFit::findFirstMax	(	const std::vector< double > &	arr,
		double	threshold
	)		const

Linear search, O(N), 1st peak should be close to the start of the vector for scintillation LAr Signals. Returns the position of the first local max

Definition at line 106 of file RiseTimeGaussFit_tool.cc.

References DEFINE_ART_CLASS_TOOL, and n.

Referenced by RiseTime().

  { 
    double max = arr[0];
    for (std::size_t i = 1, n = arr.size(); i < n; ++i) {
      if (arr[i] >= max)
        max = arr[i];

      else if (max < threshold) //didn't pass the threshold, keep searching
        max = arr[i];

      else
        return i - 1;
    }

    //No local maxima found.
    mf::LogInfo("RiseTimeGaussFit") << "No local max found above fixed threshold: " << threshold;
    return 0;
  }

double pmtana::RiseTimeGaussFit::RiseTime ( const pmtana::Waveform_t &  wf_pulse, const pmtana::PedestalMean_t &  ped_pulse, bool  _positive  ) const

overridevirtual

Implements pmtana::RiseTimeCalculatorBase.

Definition at line 60 of file RiseTimeGaussFit_tool.cc.

References util::abs(), f, findFirstMax(), fInitSigma, fMinAmp, fNbins, and fTolerance.

  {

    // Pedestal-subtracted pulse
    std::vector<double> wf_aux(ped_pulse);
    if (_positive) {
      for (size_t ix = 0; ix < wf_aux.size(); ix++) {
        wf_aux[ix] = ((double)wf_pulse[ix]) - wf_aux[ix];
      }
    }
    else {
      for (size_t ix = 0; ix < wf_aux.size(); ix++) {
        wf_aux[ix] = wf_aux[ix] - ((double)wf_pulse[ix]);
      }
    }

    // Find first local maximum
    size_t first_max = findFirstMax(wf_aux, fMinAmp);

    // Create & fill th1
    TH1F* h_aux = new TH1F("aux", "aux", wf_aux.size(), -0.5, wf_aux.size() - 0.5);
    for (long unsigned int j = 0; j < wf_aux.size(); j++)
      h_aux->SetBinContent(j + 1, wf_aux[j]);

    // Function & initial values for the fit, ROOT Gaussian:  [p0]*e**(-0.5 (x-[p1])**2 / [p2]**2)
    TF1* f = new TF1("f", "gaus", double(first_max) - fNbins, double(first_max) + fNbins);
    f->SetParameters(wf_aux[first_max], first_max, fInitSigma);
    // Fit
    h_aux->Fit(f, "q", "SAME", first_max - fNbins, first_max + fNbins);
    double t_fit = f->GetParameter(1);
    double peak_time;
    if (
      std::abs(t_fit - first_max) <
      fTolerance) { //check fit is close in distance to the original max, use max peak as time otherwise
      peak_time = t_fit;
    }
    else {
      peak_time = first_max;
      mf::LogInfo("RiseTimeGaussFit") << "No good fit found, keeping 1st max bin instead";
    }

    return peak_time;
  }

Member Data Documentation

double pmtana::RiseTimeGaussFit::fInitSigma

private

Definition at line 48 of file RiseTimeGaussFit_tool.cc.

Referenced by RiseTime(), and RiseTimeGaussFit().

double pmtana::RiseTimeGaussFit::fMinAmp

private

Definition at line 47 of file RiseTimeGaussFit_tool.cc.

Referenced by RiseTime().

int pmtana::RiseTimeGaussFit::fNbins

private

Definition at line 50 of file RiseTimeGaussFit_tool.cc.

Referenced by RiseTime(), and RiseTimeGaussFit().

double pmtana::RiseTimeGaussFit::fTolerance

private

Definition at line 49 of file RiseTimeGaussFit_tool.cc.

Referenced by RiseTime(), and RiseTimeGaussFit().

---

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

• larana/v09_15_05/source/larana/OpticalDetector/OpHitFinder/RiseTimeTools/RiseTimeGaussFit_tool.cc