#include <SymmetryFit.h>

Static Public Member Functions
static TF1 *	fitfunction ()

static std::string	name (std::string base)

static std::vector< double >	pol2_from_pol2 (TH1 *hist)

static std::vector< double >	pol2_from_pol3 (TH1 *hist)

static TH1 *	symmetryChi2 (std::string, const std::vector< TH1 * > &, const std::pair< unsigned, unsigned >)

static TH1 *	symmetryChi2 (const TH1 *, const std::pair< unsigned, unsigned >)

Private Member Functions
float	chi2 (std::pair< unsigned, unsigned >)

float	chi2_element (std::pair< unsigned, unsigned >)

std::vector< std::pair< unsigned, unsigned > >	continuousRanges () const

void	fillchi2 ()

std::pair< unsigned, unsigned >	findUsableMinMax () const

int	fit ()

void	makeChi2Histogram ()

SymmetryFit	operator+= (const SymmetryFit &R)

	SymmetryFit (const TH1 *, const std::pair< unsigned, unsigned >)

Private Attributes
TH1 *	chi2_

const unsigned	minDF_

const std::pair< unsigned, unsigned >	minmaxUsable_

unsigned	ndf_

const std::pair< unsigned, unsigned >	range_

const TH1 *	symm_candidate_

Detailed Description

Definition at line 8 of file SymmetryFit.h.

Constructor & Destructor Documentation

SymmetryFit::SymmetryFit	(	const TH1 *	h,
		const std::pair< unsigned, unsigned >	r
	)

private

Definition at line 38 of file SymmetryFit.cc.

References fillchi2(), and makeChi2Histogram().

Referenced by name().

     : symm_candidate_(h),
       minDF_(r.second - r.first),
       range_(r),
       minmaxUsable_(findUsableMinMax()),
       ndf_(minmaxUsable_.first < minmaxUsable_.second ? minmaxUsable_.second - minmaxUsable_.first : 0),
       chi2_(nullptr) {
   makeChi2Histogram();
   fillchi2();
 }

Member Function Documentation

float SymmetryFit::chi2 ( std::pair< unsigned, unsigned > point )

private

Definition at line 106 of file SymmetryFit.cc.

References chi2_element(), mps_fire::i, minmaxUsable_, ndf_, and XX.

Referenced by fillchi2(), and name().

                                                          {
   point.first -= minmaxUsable_.first;
   point.second += minmaxUsable_.first;
   float XX = 0;
   unsigned i = ndf_;
   while (i-- > 0) {
     XX += chi2_element(point);
     point.first--;
     point.second++;
   }
   return XX;
 }

float SymmetryFit::chi2_element ( std::pair< unsigned, unsigned > range )

private

Definition at line 119 of file SymmetryFit.cc.

References StorageManager_cfg::e1, funct::pow(), symm_candidate_, and testProducerWithPsetDescEmpty_cfi::y1.

Referenced by chi2(), and name().

                                                                  {
   float y0(symm_candidate_->GetBinContent(range.first)), y1(symm_candidate_->GetBinContent(range.second)),
       e0(symm_candidate_->GetBinError(range.first)), e1(symm_candidate_->GetBinError(range.second));
   assert(e0 && e1);
 
   return pow(y0 - y1, 2) / (e0 * e0 + e1 * e1);
 }

std::vector< std::pair< unsigned, unsigned > > SymmetryFit::continuousRanges ( ) const

private

Definition at line 81 of file SymmetryFit.cc.

References mps_fire::i, SiStripMonitorCluster_cfi::Nbins, FastTimerService_cff::range, diffTwoXMLs::ranges, and symm_candidate_.

Referenced by findUsableMinMax(), and name().

                                                                             {
   std::vector<std::pair<unsigned, unsigned> > ranges;
   const unsigned Nbins = symm_candidate_->GetNbinsX();
   unsigned i = 0;
   while (++i <= Nbins) {
     if (symm_candidate_->GetBinError(i)) {
       std::pair<unsigned, unsigned> range(i, i + 1);
       while (++i <= Nbins && symm_candidate_->GetBinError(i))
         range.second++;
       ranges.push_back(range);
     }
   }
   return ranges;
 }

void SymmetryFit::fillchi2 ( )

private

Definition at line 96 of file SymmetryFit.cc.

References chi2(), chi2_, mps_fire::i, dttmaxenums::L, minDF_, ndf_, dttmaxenums::R, and range_.

Referenced by name(), and SymmetryFit().

                            {
   if (ndf_ < minDF_)
     return;
 
   for (int i = 1; i <= chi2_->GetNbinsX(); ++i) {
     const unsigned L(range_.first - 1 + (i - 1) / 2), R(range_.first + i / 2);
     chi2_->SetBinContent(i, chi2(std::make_pair(L, R)));
   }
 }

std::pair< unsigned, unsigned > SymmetryFit::findUsableMinMax ( ) const

private

Definition at line 58 of file SymmetryFit.cc.

References generateEDF::cont, continuousRanges(), dqmdumpme::first, dttmaxenums::L, SiStripPI::max, min(), dttmaxenums::R, range_, and edm::second().

Referenced by name().

                                                                 {
   unsigned NEAR(0), FAR(0);
   const std::vector<std::pair<unsigned, unsigned> > cont = continuousRanges();
   for (unsigned L = 0; L < cont.size(); L++) {
     if (cont[L].first > range_.first)
       continue;
     for (unsigned R = L; R < cont.size(); R++) {
       if (cont[R].second < range_.second)
         continue;
 
       const unsigned far = std::min(range_.first - cont[L].first, cont[R].second - range_.second);
       const unsigned near = std::max(range_.second < cont[L].second ? 0 : range_.second - cont[L].second,
                                      cont[R].first < range_.first ? 0 : cont[R].first - range_.first);
 
       if ((far > near) && (far - near) > (FAR - NEAR)) {
         FAR = far;
         NEAR = near;
       }
     }
   }
   return std::make_pair(NEAR, FAR);
 }

int SymmetryFit::fit ( )

private

Definition at line 127 of file SymmetryFit.cc.

References chi2_, f, fitfunction(), ndf_, AlCaHLTBitMon_ParallelJobs::p, and pol2_from_pol3().

Referenced by trackingPlots.Iteration::modules(), name(), and symmetryChi2().

                      {
   std::vector<double> p = pol2_from_pol3(chi2_);
   if (p.empty() || p[0] < chi2_->GetBinCenter(1) || p[0] > chi2_->GetBinCenter(chi2_->GetNbinsX()))
     return 7;
 
   std::unique_ptr<TF1> f(fitfunction());
   f->SetParameter(0, p[0]);
   f->SetParLimits(0, p[0], p[0]);
   f->SetParameter(1, p[1]);
   f->SetParLimits(1, p[1], p[1]);
   f->SetParameter(2, p[2]);
   f->SetParLimits(2, p[2], p[2]);
   f->SetParameter(3, ndf_);
   f->SetParLimits(3, ndf_, ndf_);  //Fixed
   chi2_->Fit(f.get(), "WQ");
   return 0;
 }

TF1 * SymmetryFit::fitfunction ( )

static

Definition at line 145 of file SymmetryFit.cc.

References f.

Referenced by fit(), and name().

                               {
   TF1* f = new TF1("SymmetryFit", "((x-[0])/[1])**2+[2]+0*[3]");
   f->SetParName(0, "middle");
   f->SetParName(1, "uncertainty");
   f->SetParName(2, "chi2");
   f->SetParName(3, "NDF");
   return f;
 }

void SymmetryFit::makeChi2Histogram ( )

private

Definition at line 49 of file SymmetryFit.cc.

References chi2_, LaserClient_cfi::low, name(), SiStripMonitorCluster_cfi::Nbins, range_, AlCaHLTBitMon_QueryRunRegistry::string, symm_candidate_, and up.

Referenced by name(), and SymmetryFit().

                                     {
   std::string XXname = name(symm_candidate_->GetName());
   unsigned Nbins = 2 * (range_.second - range_.first) + 1;
   double binwidth = symm_candidate_->GetBinWidth(1);
   double low = symm_candidate_->GetBinCenter(range_.first) - 3 * binwidth / 4;
   double up = symm_candidate_->GetBinCenter(range_.second - 1) + 3 * binwidth / 4;
   chi2_ = new TH1F(XXname.c_str(), "", Nbins, low, up);
 }

static std::string SymmetryFit::name ( std::string base )

inlinestatic

Definition at line 12 of file SymmetryFit.h.

References chi2(), chi2_element(), continuousRanges(), fillchi2(), findUsableMinMax(), fit(), fitfunction(), compare::hist, makeChi2Histogram(), pol2_from_pol2(), pol2_from_pol3(), and SymmetryFit().

Referenced by config.CFG::__str__(), validation.Sample::digest(), VIDSelectorBase.VIDSelectorBase::initialize(), makeChi2Histogram(), LA_Filler_Fitter::method(), and Vispa.Views.PropertyView.Property::valueChanged().

12 { return base + "_symmchi2"; }

newFWLiteAna.base

base

Main Program

Definition: newFWLiteAna.py:92

SymmetryFit SymmetryFit::operator+= ( const SymmetryFit & R )

inlineprivate

Definition at line 26 of file SymmetryFit.h.

References chi2_, and ndf_.

                                                {
     ndf_ += R.ndf_;
     chi2_->Add(R.chi2_);
     return *this;
   }

std::vector< double > SymmetryFit::pol2_from_pol2 ( TH1 * hist )

static

Definition at line 154 of file SymmetryFit.cc.

References TrackCollections2monitor_cff::func, AlCaHLTBitMon_ParallelJobs::p, mathSSE::sqrt(), mps_update::status, and findQualityFiles::v.

Referenced by name().

                                                        {
   std::vector<double> v;
 
   //Need our own copy for thread safety
   TF1 func("mypol2", "pol2");
   int status = hist->Fit(&func, "WQ");
   if (!status) {
     std::vector<double> p;
     p.push_back(func.GetParameter(0));
     p.push_back(func.GetParameter(1));
     p.push_back(func.GetParameter(2));
     if (p[2] > 0) {
       v.push_back(-0.5 * p[1] / p[2]);
       v.push_back(1. / sqrt(p[2]));
       v.push_back(p[0] - 0.25 * p[1] * p[1] / p[2]);
     }
   }
   return v;
 }

std::vector< double > SymmetryFit::pol2_from_pol3 ( TH1 * hist )

static

Definition at line 174 of file SymmetryFit.cc.

References TrackCollections2monitor_cff::func, AlCaHLTBitMon_ParallelJobs::p, funct::pow(), mathSSE::sqrt(), mps_update::status, and findQualityFiles::v.

Referenced by fit(), and name().

                                                        {
   std::vector<double> v;
 
   auto func = std::make_unique<TF1>("mypol3", "pol3");
   int status = hist->Fit(func.get(), "WQ");
   if (!status) {
     std::vector<double> p;
     p.push_back(func->GetParameter(0));
     p.push_back(func->GetParameter(1));
     p.push_back(func->GetParameter(2));
     p.push_back(func->GetParameter(3));
     double radical = p[2] * p[2] - 3 * p[1] * p[3];
     if (radical > 0) {
       double x0 = (-p[2] + sqrt(radical)) / (3 * p[3]);
       v.push_back(x0);
       v.push_back(pow(radical, -0.25));
       v.push_back(p[0] + p[1] * x0 + p[2] * x0 * x0 + p[3] * x0 * x0 * x0);
     }
   }
   return v;
 }

TH1 * SymmetryFit::symmetryChi2	(	std::string	basename,
		const std::vector< TH1 * > &	candidates,
		const std::pair< unsigned, unsigned >	range
	)

static

Definition at line 6 of file SymmetryFit.cc.

References chi2_, fit(), and mps_update::status.

Referenced by LA_Filler_Fitter::make_and_fit_symmchi2().

                                                                         {
   TH1* fake = (TH1*)(candidates[0]->Clone(basename.c_str()));
   fake->Reset();
   SymmetryFit combined(fake, range);
   delete fake;
 
   for (auto const* candidate : candidates) {
     SymmetryFit sf(candidate, range);
     combined += sf;
     delete sf.chi2_;
   }
 
   int status = combined.fit();
   if (status) {
     delete combined.chi2_;
     combined.chi2_ = nullptr;
   }
   return combined.chi2_;
 }

TH1 * SymmetryFit::symmetryChi2	(	const TH1 *	candidate,
		const std::pair< unsigned, unsigned >	range
	)

static

Definition at line 28 of file SymmetryFit.cc.

References chi2_, fit(), and mps_update::status.

                                                                                             {
   SymmetryFit sf(candidate, range);
   int status = sf.fit();
   if (status) {
     delete sf.chi2_;
     sf.chi2_ = nullptr;
   }
   return sf.chi2_;
 }