DTResidualFitter Class Reference

#include <DTResidualFitter.h>

Public Member Functions
	DTResidualFitter (bool debug=false)
DTResidualFitResult	fitResiduals (TH1F &histo, int nSigmas=1)
	~DTResidualFitter ()

Private Attributes
bool	debug_

Detailed Description

Definition at line 21 of file DTResidualFitter.h.

Constructor & Destructor Documentation

DTResidualFitter()

DTResidualFitter::DTResidualFitter

(

bool

debug = false

)

Definition at line 15 of file DTResidualFitter.cc.

: debug_(debug) {}

~DTResidualFitter()

DTResidualFitter::~DTResidualFitter

(

)

Definition at line 17 of file DTResidualFitter.cc.

{}

Member Function Documentation

fitResiduals()

DTResidualFitResult DTResidualFitter::fitResiduals	(	TH1F &	histo,
		int	nSigmas = 1
	)

Definition at line 19 of file DTResidualFitter.cc.

References debug_, timingPdfMaker::histo, and fileinputsource_cfi::option.

Referenced by dtCalibration::DTVDriftSegment::compute().

                                                                           {
  TString option("R");
  if (!debug_)
    option += "Q";

  float under = histo.GetBinContent(0) / histo.GetEntries();
  float over = histo.GetBinContent(histo.GetNbinsX() + 1) / histo.GetEntries();
  float minFit = histo.GetMean() - histo.GetRMS();
  float maxFit = histo.GetMean() + histo.GetRMS();

  if ((under > 0.1) || (over > 0.1))
    edm::LogError("DTResidualFitter") << "WARNING in histogram: " << histo.GetName() << "\n"
                                      << "             entries: " << histo.GetEntries() << "\n"
                                      << "           underflow: " << under * 100. << "% \n"
                                      << "            overflow: " << over * 100. << "%";

  TString funcName = TString(histo.GetName()) + "_gaus";
  TF1* fitFunc = new TF1(funcName, "gaus", minFit, maxFit);

  histo.Fit(fitFunc, option);

  minFit = fitFunc->GetParameter(1) - nSigmas * fitFunc->GetParameter(2);
  maxFit = fitFunc->GetParameter(1) + nSigmas * fitFunc->GetParameter(2);
  fitFunc->SetRange(minFit, maxFit);
  histo.Fit(fitFunc, option);

  return DTResidualFitResult(
      fitFunc->GetParameter(1), fitFunc->GetParError(1), fitFunc->GetParameter(2), fitFunc->GetParError(2));
}

Member Data Documentation

debug_

bool DTResidualFitter::debug_

private

Definition at line 29 of file DTResidualFitter.h.

Referenced by fitResiduals().