All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
Public Member Functions | Static Public Member Functions | Private Attributes | List of all members
reco_tool::PeakFitterICARUS Class Reference
Inheritance diagram for reco_tool::PeakFitterICARUS:
reco_tool::IPeakFitter

Public Member Functions

 PeakFitterICARUS (const fhicl::ParameterSet &pset)
 
void configure (const fhicl::ParameterSet &pset) override
 
void findPeakParameters (const std::vector< float > &, const ICandidateHitFinder::HitCandidateVec &, PeakParamsVec &, double &, int &) const override
 

Static Public Member Functions

static Double_t fitf (Double_t *x, Double_t *par)
 

Private Attributes

double fMinWidth
 minimum initial width for ICARUS fit More...
 
double fMaxWidthMult
 multiplier for max width for ICARUS fit More...
 
double fPeakRange
 set range limits for peak center More...
 
double fAmpRange
 set range limit for peak amplitude More...
 
TH1F fHistogram
 
TF1 fFit
 Cache of fit functions (one so far). More...
 
const geo::GeometryCore * fGeometry = lar::providerFrom<geo::Geometry>()
 
- Private Attributes inherited from reco_tool::IPeakFitter
float peakCenterError
 
float peakSigma
 
float peakSigmaError
 
float peakAmplitude
 
float peakAmplitudeError
 
float peakTauLeft
 
float peakTauLeftError
 
float peakTauRight
 
float peakTauRightError
 
float peakBaseline
 
float peakBaselineError
 

Additional Inherited Members

- Private Types inherited from reco_tool::IPeakFitter
using PeakFitParams_t = struct PeakFitParams{float peakCenter
 
using PeakParamsVec = std::vector< PeakFitParams_t >
 
using PeakParamsVec = std::vector< PeakFitParams_t >
 
- Private Member Functions inherited from reco_tool::IPeakFitter
virtual ~IPeakFitter () noexcept=default
 
virtual ~IPeakFitter ()=default
 

Detailed Description

Definition at line 24 of file PeakFitterICARUS_tool.cc.

Constructor & Destructor Documentation

reco_tool::PeakFitterICARUS::PeakFitterICARUS ( const fhicl::ParameterSet &  pset)
explicit

Definition at line 54 of file PeakFitterICARUS_tool.cc.

55  : fFit("ICARUSfunc", fitf, 0.0, 1.0, 5) // 5 parameters, fake range
56 {
57  configure(pset);
58 }
reco_tool::PeakFitterICARUS::configure
void configure(const fhicl::ParameterSet &pset) override
Definition: PeakFitterICARUS_tool.cc:60
reco_tool::PeakFitterICARUS::fitf
static Double_t fitf(Double_t *x, Double_t *par)
Definition: PeakFitterICARUS_tool.cc:188
reco_tool::PeakFitterICARUS::fFit
TF1 fFit
Cache of fit functions (one so far).
Definition: PeakFitterICARUS_tool.cc:47

Member Function Documentation

void reco_tool::PeakFitterICARUS::configure ( const fhicl::ParameterSet &  pset)
overridevirtual

Implements reco_tool::IPeakFitter.

Definition at line 60 of file PeakFitterICARUS_tool.cc.

61 {
62  // Start by recovering the parameters
63  fMinWidth = pset.get<double>("MinWidth", 0.5);
64  fMaxWidthMult = pset.get<double>("MaxWidthMult", 3.);
65  fPeakRange = pset.get<double>("PeakRangeFact", 2.);
66  fAmpRange = pset.get<double>("PeakAmpRange", 2.);
67 
68  fHistogram = TH1F("PeakFitterHitSignal","",500,0.,500.);
69 
70  fHistogram.Sumw2();
71 
72  return;
73 }
reco_tool::PeakFitterICARUS::fAmpRange
double fAmpRange
set range limit for peak amplitude
Definition: PeakFitterICARUS_tool.cc:44
reco_tool::PeakFitterICARUS::fMinWidth
double fMinWidth
minimum initial width for ICARUS fit
Definition: PeakFitterICARUS_tool.cc:41
reco_tool::PeakFitterICARUS::fPeakRange
double fPeakRange
set range limits for peak center
Definition: PeakFitterICARUS_tool.cc:43
reco_tool::PeakFitterICARUS::fMaxWidthMult
double fMaxWidthMult
multiplier for max width for ICARUS fit
Definition: PeakFitterICARUS_tool.cc:42
reco_tool::PeakFitterICARUS::fHistogram
TH1F fHistogram
Definition: PeakFitterICARUS_tool.cc:46
void reco_tool::PeakFitterICARUS::findPeakParameters ( const std::vector< float > &  roiSignalVec,
const ICandidateHitFinder::HitCandidateVec &  hitCandidateVec,
PeakParamsVec &  peakParamsVec,
double &  chi2PerNDF,
int &  NDF 
) const
overridevirtual

Implements reco_tool::IPeakFitter.

Definition at line 76 of file PeakFitterICARUS_tool.cc.

81 {
82  // The following is a translation of the original FitICARUSs function in the original
83  // GausHitFinder module originally authored by Jonathan Asaadi
84  //
85  // *** NOTE: this algorithm assumes the reference time for input hit candidates is to
86  // the first tick of the input waveform (ie 0)
87  //
88  if (hitCandidateVec.empty()) return;
89 
90  // in case of a fit failure, set the chi-square to infinity
91  chi2PerNDF = std::numeric_limits<double>::infinity();
92 
93  int startTime = hitCandidateVec.front().startTick;
94  int endTime = hitCandidateVec.back().stopTick;
95  int roiSize = endTime - startTime;
96 
97  // Check to see if we need a bigger histogram for fitting
98  if (roiSize > fHistogram.GetNbinsX())
99  {
100  std::string histName = "PeakFitterHitSignal_" + std::to_string(roiSize);
101  fHistogram = TH1F(histName.c_str(),"",roiSize,0.,roiSize);
102  fHistogram.Sumw2();
103  }
104 
105  for(int idx = 0; idx < roiSize; idx++) fHistogram.SetBinContent(idx+1,roiSignalVec.at(startTime+idx));
106 
107  // for(int idx = 0; idx < roiSize; idx++)
108  // std::cout << " bin " << idx << " fHistogram " << fHistogram.GetBinContent(idx+1) << std::endl;
109 
110 
111  // ### Setting the parameters for the ICARUS Fit ###
112  //int parIdx(0);
113  for(auto& candidateHit : hitCandidateVec)
114  {
115  double peakMean = candidateHit.hitCenter - float(startTime);
116  double peakWidth = candidateHit.hitSigma;
117  // std::cout << " peakWidth " << peakWidth << std::endl;
118  // std::cout << " hitcenter " << candidateHit.hitCenter << " starttime " << startTime <<std::endl;
119 
120 
121  double amplitude = candidateHit.hitHeight;
122  double meanLowLim = std::max(peakMean - fPeakRange * peakWidth, 0.);
123  double meanHiLim = std::min(peakMean + fPeakRange * peakWidth, double(roiSize));
124  // std::cout << " amplitude " << amplitude << std::endl;
125  // std::cout << " peakMean " << peakMean << std::endl;
126 
127  fFit.SetParameter(0,0);
128  fFit.SetParameter(1, amplitude);
129  fFit.SetParameter(2, peakMean);
130  fFit.SetParameter(3,peakWidth/2);
131  fFit.SetParameter(4,peakWidth/2);
132 
133  fFit.SetParLimits(0, -5, 5);
134  fFit.SetParLimits(1, 0.1 * amplitude, 10. * amplitude);
135  fFit.SetParLimits(2, meanLowLim,meanHiLim);
136  fFit.SetParLimits(3, std::max(fMinWidth, 0.1 * peakWidth), fMaxWidthMult * peakWidth);
137  fFit.SetParLimits(4, std::max(fMinWidth, 0.1 * peakWidth), fMaxWidthMult * peakWidth);
138 
139  // std::cout << " peakWidth limits " << std::max(fMinWidth, 0.1 * peakWidth) <<" " << fMaxWidthMult * peakWidth << std::endl;
140  // std::cout << " peakMean limits " << meanLowLim <<" " << meanHiLim << std::endl;
141 
142  }
143 
144  int fitResult(-1);
145 
146  // the range of the fit does not matter since we specify the fitting range
147  // explicitly (we do NOT use option "R")
148  try
149  { fitResult = fHistogram.Fit(&fFit,"QNWB","", 0., roiSize);}
150  catch(...)
151  {mf::LogWarning("GausHitFinder") << "Fitter failed finding a hit";}
152 
153  if(fitResult!=0)
154 // std::cout << " fit cannot converge " << std::endl;
155  // ##################################################
156  // ### Getting the fitted parameters from the fit ###
157  // ##################################################
158  NDF = roiSize-5;
159  chi2PerNDF = (fFit.GetChisquare() / NDF);
160 
161  // std::cout << " chi2 " << fFit.GetChisquare() << std::endl;
162  // std::cout << " ndf " << NDF << std::endl;
163 
164  // std::cout << " chi2ndf " << chi2PerNDF<< std::endl;
165  //parIdx = 0;
166  for(size_t idx = 0; idx < hitCandidateVec.size(); idx++)
167  {
168  PeakFitParams_t peakParams;
169 
170  peakParams.peakAmplitude = fFit.GetParameter(1);
171  peakParams.peakAmplitudeError = fFit.GetParError(1);
172  peakParams.peakCenter = fFit.GetParameter(2) + float(startTime);
173  peakParams.peakCenterError = fFit.GetParError(2);
174  //std::cout << " rising time " << fFit.GetParameter(3) << " falling time " <<fFit.GetParameter(4) << std::endl;
175  peakParams.peakTauLeft = fFit.GetParameter(3);
176  peakParams.peakTauLeftError = fFit.GetParError(3);
177  peakParams.peakTauRight = fFit.GetParameter(4);
178  peakParams.peakTauRightError = fFit.GetParError(4);
179  peakParams.peakBaseline = fFit.GetParameter(0);
180  peakParams.peakBaselineError = fFit.GetParError(0);
181 
182  peakParamsVec.emplace_back(peakParams);
183 
184  }
185 
186 }
reco_tool::PeakFitterICARUS::fMinWidth
double fMinWidth
minimum initial width for ICARUS fit
Definition: PeakFitterICARUS_tool.cc:41
reco_tool::IPeakFitter::PeakFitParams_t
struct PeakFitParams{float peakCenter PeakFitParams_t
Definition: icaruscode/icaruscode/TPC/SignalProcessing/HitFinder/HitFinderTools/IPeakFitter.h:31
reco_tool::PeakFitterICARUS::fPeakRange
double fPeakRange
set range limits for peak center
Definition: PeakFitterICARUS_tool.cc:43
reco_tool::PeakFitterICARUS::fMaxWidthMult
double fMaxWidthMult
multiplier for max width for ICARUS fit
Definition: PeakFitterICARUS_tool.cc:42
icarus::trigger::to_string
std::string to_string(WindowPattern const &pattern)
Definition: WindowPattern.h:115
reco_tool::PeakFitterICARUS::fHistogram
TH1F fHistogram
Definition: PeakFitterICARUS_tool.cc:46
reco_tool::PeakFitterICARUS::fFit
TF1 fFit
Cache of fit functions (one so far).
Definition: PeakFitterICARUS_tool.cc:47
Double_t reco_tool::PeakFitterICARUS::fitf ( Double_t *  x,
Double_t *  par 
)
static

Definition at line 188 of file PeakFitterICARUS_tool.cc.

189  {
190  // Double_t arg = 0;
191 
192  Double_t fitval = par[0]+par[1]*TMath::Exp(-(x[0]-par[2])/par[3])/(1+TMath::Exp(-(x[0]-par[3])/par[4]));
193  return fitval;
194  }
x
process_name opflash particleana ie x
Definition: run_opflash_electron.fcl:90

Member Data Documentation

double reco_tool::PeakFitterICARUS::fAmpRange
private

set range limit for peak amplitude

Definition at line 44 of file PeakFitterICARUS_tool.cc.

TF1 reco_tool::PeakFitterICARUS::fFit
mutableprivate

Cache of fit functions (one so far).

Definition at line 47 of file PeakFitterICARUS_tool.cc.

const geo::GeometryCore* reco_tool::PeakFitterICARUS::fGeometry = lar::providerFrom<geo::Geometry>()
private

Definition at line 49 of file PeakFitterICARUS_tool.cc.

TH1F reco_tool::PeakFitterICARUS::fHistogram
mutableprivate

Definition at line 46 of file PeakFitterICARUS_tool.cc.

double reco_tool::PeakFitterICARUS::fMaxWidthMult
private

multiplier for max width for ICARUS fit

Definition at line 42 of file PeakFitterICARUS_tool.cc.

double reco_tool::PeakFitterICARUS::fMinWidth
private

minimum initial width for ICARUS fit

Definition at line 41 of file PeakFitterICARUS_tool.cc.

double reco_tool::PeakFitterICARUS::fPeakRange
private

set range limits for peak center

Definition at line 43 of file PeakFitterICARUS_tool.cc.

The documentation for this class was generated from the following file:
Generated by   doxygen 1.8.5