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 23 of file DTResidualFitter.h.

Constructor & Destructor Documentation

DTResidualFitter::DTResidualFitter

(

bool

debug = false

)

Definition at line 15 of file DTResidualFitter.cc.

:debug_(debug) {}

DTResidualFitter::~DTResidualFitter

(

)

Definition at line 17 of file DTResidualFitter.cc.

{}

Member Function Documentation

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

Definition at line 19 of file DTResidualFitter.cc.

References debug_, and TSGForRoadSearch_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

bool DTResidualFitter::debug_

private

Definition at line 31 of file DTResidualFitter.h.

Referenced by fitResiduals().