#include "Alignment/MuonAlignmentAlgorithms/interface/MuonResidualsFitter.h"
#include <fstream>
#include <set>
#include "TMath.h"
#include "TH1.h"
#include "TF1.h"
#include "TRobustEstimator.h"

Functions
double	MuonResidualsFitter_compute_log_convolution (double toversigma, double gammaoversigma, double max, double step, double power)

Double_t	MuonResidualsFitter_GaussPowerTails_TF1 (Double_t xvec, Double_t par)

double	MuonResidualsFitter_integrate_pureGaussian (double low, double high, double center, double sigma)

double	MuonResidualsFitter_logGaussPowerTails (double residual, double center, double sigma)

double	MuonResidualsFitter_logPowerLawTails (double residual, double center, double sigma, double gamma)

double	MuonResidualsFitter_logPureGaussian (double residual, double center, double sigma)

double	MuonResidualsFitter_logPureGaussian2D (double x, double y, double x0, double y0, double sx, double sy, double r)

double	MuonResidualsFitter_logROOTVoigt (double residual, double center, double sigma, double gamma)

Double_t	MuonResidualsFitter_powerLawTails_TF1 (Double_t xvec, Double_t par)

Double_t	MuonResidualsFitter_pureGaussian_TF1 (Double_t xvec, Double_t par)

Double_t	MuonResidualsFitter_ROOTVoigt_TF1 (Double_t xvec, Double_t par)

Variables
const double	MuonResidualsFitter_gsbinsize = 0.01

double	MuonResidualsFitter_lookup_table [MuonResidualsFitter_numgsbins][MuonResidualsFitter_numtsbins]

const int	MuonResidualsFitter_numgsbins = 500

const int	MuonResidualsFitter_numtsbins = 500

bool	MuonResidualsFitter_table_initialized = false

static TMinuit *	MuonResidualsFitter_TMinuit

const double	MuonResidualsFitter_tsbinsize = 0.1

Function Documentation

◆ MuonResidualsFitter_compute_log_convolution()

double MuonResidualsFitter_compute_log_convolution	(	double	toversigma,
		double	gammaoversigma,
		double	max,
		double	step,
		double	power
	)

Definition at line 54 of file MuonResidualsFitter.cc.

References JetChargeProducer_cfi::exp, dqm-mbProfile::log, M_PI, SiStripPI::max, funct::pow(), cms::alpakatools::detail::power(), and mathSSE::sqrt().

Referenced by MuonResidualsFitter::initialize_table().

                                                                                      {
   if (gammaoversigma == 0.)
     return (-toversigma * toversigma / 2.) - log(sqrt(2 * M_PI));
 
   double sum = 0.;
   double uplus = 0.;
   double integrandplus_last = 0.;
   double integrandminus_last = 0.;
   for (double inc = 0.; uplus < max; inc += step) {
     double uplus_last = uplus;
     uplus = pow(inc, power);
 
     double integrandplus = exp(-pow(uplus - toversigma, 2) / 2.) / (uplus * uplus / gammaoversigma + gammaoversigma) /
                            M_PI / sqrt(2. * M_PI);
     double integrandminus = exp(-pow(-uplus - toversigma, 2) / 2.) / (uplus * uplus / gammaoversigma + gammaoversigma) /
                             M_PI / sqrt(2. * M_PI);
 
     sum += integrandplus * (uplus - uplus_last);
     sum += 0.5 * fabs(integrandplus_last - integrandplus) * (uplus - uplus_last);
 
     sum += integrandminus * (uplus - uplus_last);
     sum += 0.5 * fabs(integrandminus_last - integrandminus) * (uplus - uplus_last);
 
     integrandplus_last = integrandplus;
     integrandminus_last = integrandminus;
   }
   return log(sum);
 }

◆ MuonResidualsFitter_GaussPowerTails_TF1()

Double_t MuonResidualsFitter_GaussPowerTails_TF1	(	Double_t *	xvec,
		Double_t *	par
	)

Definition at line 142 of file MuonResidualsFitter.cc.

References JetChargeProducer_cfi::exp, and MuonResidualsFitter_logGaussPowerTails().

Referenced by MuonResiduals1DOFFitter::plot(), MuonResidualsPositionFitter::plot(), MuonResidualsAngleFitter::plot(), MuonResidualsBfieldAngleFitter::plot(), MuonResiduals6DOFrphiFitter::plot(), MuonResiduals5DOFFitter::plot(), and MuonResiduals6DOFFitter::plot().

                                                                                 {
   return par[0] * exp(MuonResidualsFitter_logGaussPowerTails(xvec[0], par[1], par[2]));
 }

◆ MuonResidualsFitter_integrate_pureGaussian()

double MuonResidualsFitter_integrate_pureGaussian	(	double	low,
		double	high,
		double	center,
		double	sigma
	)

Definition at line 146 of file MuonResidualsFitter.cc.

References JetChargeProducer_cfi::exp, LaserClient_cfi::high, LaserClient_cfi::low, and mathSSE::sqrt().

                                                                                                         {
   static const double isqr2 = 1. / sqrt(2.);
   return (erf((high + center) * isqr2 / sigma) - erf((low + center) * isqr2 / sigma)) * exp(0.5 / sigma / sigma) * 0.5;
 }

◆ MuonResidualsFitter_logGaussPowerTails()

double MuonResidualsFitter_logGaussPowerTails	(	double	residual,
		double	center,
		double	sigma
	)

Definition at line 128 of file MuonResidualsFitter.cc.

References a, JetChargeProducer_cfi::exp, dqm-mbProfile::log, visualization-live-secondInstance_cfg::m, M_PI, dqmiodumpmetadata::n, funct::pow(), alignCSCRings::s, mathSSE::sqrt(), and x.

Referenced by MuonResiduals1DOFFitter_FCN(), MuonResiduals5DOFFitter_FCN(), MuonResiduals6DOFFitter_FCN(), MuonResiduals6DOFrphiFitter_FCN(), MuonResidualsAngleFitter_FCN(), MuonResidualsBfieldAngleFitter_FCN(), MuonResidualsFitter_GaussPowerTails_TF1(), and MuonResidualsPositionFitter_FCN().

                                                                                             {
   double x = residual - center;
   double s = fabs(sigma);
   double m = 2. * s;
   static const double e2 = exp(-2.), sqerf = sqrt(2. * M_PI) * erf(sqrt(2.));
   double a = pow(m, 4) * e2;
   double n = sqerf * s + 2. * a * pow(m, -3) / 3.;
 
   if (fabs(x) < m)
     return -x * x / (2. * s * s) - log(n);
   else
     return log(a) - 4. * log(fabs(x)) - log(n);
 }

◆ MuonResidualsFitter_logPowerLawTails()

double MuonResidualsFitter_logPowerLawTails	(	double	residual,
		double	center,
		double	sigma,
		double	gamma
	)

Definition at line 85 of file MuonResidualsFitter.cc.

References reco::ceil(), CustomPhysics_cfi::gamma, createfilelist::int, dqm-mbProfile::log, M_PI, MuonResidualsFitter_gsbinsize, MuonResidualsFitter_lookup_table, MuonResidualsFitter_numgsbins, MuonResidualsFitter_numtsbins, MuonResidualsFitter_tsbinsize, and heppy_batch::val.

Referenced by MuonResiduals1DOFFitter_FCN(), MuonResiduals5DOFFitter_FCN(), MuonResiduals6DOFFitter_FCN(), MuonResiduals6DOFrphiFitter_FCN(), MuonResidualsAngleFitter_FCN(), MuonResidualsBfieldAngleFitter_FCN(), MuonResidualsFitter_powerLawTails_TF1(), and MuonResidualsPositionFitter_FCN().

                                                                                                         {
   sigma = fabs(sigma);
   double gammaoversigma = fabs(gamma / sigma);
   double toversigma = fabs((residual - center) / sigma);
 
   int gsbin0 = int(floor(gammaoversigma / MuonResidualsFitter_gsbinsize));
   int gsbin1 = int(ceil(gammaoversigma / MuonResidualsFitter_gsbinsize));
   int tsbin0 = int(floor(toversigma / MuonResidualsFitter_tsbinsize));
   int tsbin1 = int(ceil(toversigma / MuonResidualsFitter_tsbinsize));
 
   bool gsisbad = (gsbin0 >= MuonResidualsFitter_numgsbins || gsbin1 >= MuonResidualsFitter_numgsbins);
   bool tsisbad = (tsbin0 >= MuonResidualsFitter_numtsbins || tsbin1 >= MuonResidualsFitter_numtsbins);
 
   if (gsisbad || tsisbad) {
     return log(gammaoversigma / M_PI) - log(toversigma * toversigma + gammaoversigma * gammaoversigma) - log(sigma);
   } else {
     double val00 = MuonResidualsFitter_lookup_table[gsbin0][tsbin0];
     double val01 = MuonResidualsFitter_lookup_table[gsbin0][tsbin1];
     double val10 = MuonResidualsFitter_lookup_table[gsbin1][tsbin0];
     double val11 = MuonResidualsFitter_lookup_table[gsbin1][tsbin1];
 
     double val0 = val00 + ((toversigma / MuonResidualsFitter_tsbinsize) - tsbin0) * (val01 - val00);
     double val1 = val10 + ((toversigma / MuonResidualsFitter_tsbinsize) - tsbin0) * (val11 - val10);
 
     double val = val0 + ((gammaoversigma / MuonResidualsFitter_gsbinsize) - gsbin0) * (val1 - val0);
 
     return val - log(sigma);
   }
 }

◆ MuonResidualsFitter_logPureGaussian()

double MuonResidualsFitter_logPureGaussian	(	double	residual,
		double	center,
		double	sigma
	)

Definition at line 31 of file MuonResidualsFitter.cc.

References dqm-mbProfile::log, M_PI, and funct::pow().

Referenced by MuonResiduals1DOFFitter_FCN(), MuonResiduals5DOFFitter_FCN(), MuonResiduals6DOFFitter_FCN(), MuonResiduals6DOFrphiFitter_FCN(), MuonResidualsAngleFitter_FCN(), MuonResidualsBfieldAngleFitter_FCN(), MuonResidualsFitter_pureGaussian_TF1(), and MuonResidualsPositionFitter_FCN().

                                                                                          {
   sigma = fabs(sigma);
   static const double cgaus = 0.5 * log(2. * M_PI);
   return (-pow(residual - center, 2) * 0.5 / sigma / sigma) - cgaus - log(sigma);
 }

◆ MuonResidualsFitter_logPureGaussian2D()

double MuonResidualsFitter_logPureGaussian2D	(	double	x,
		double	y,
		double	x0,
		double	y0,
		double	sx,
		double	sy,
		double	r
	)

Definition at line 43 of file MuonResidualsFitter.cc.

References PVValHelper::dx, PVValHelper::dy, dqm-mbProfile::log, M_PI, funct::pow(), alignCSCRings::r, mathSSE::sqrt(), fftjetcommon_cfi::sx, fftjetcommon_cfi::sy, x, and detailsBasic3DVector::y.

Referenced by MuonResiduals5DOFFitter_FCN(), MuonResiduals6DOFFitter_FCN(), and MuonResiduals6DOFrphiFitter_FCN().

                                                                                                                        {
   sx = fabs(sx);
   sy = fabs(sy);
   static const double c2gaus = log(2. * M_PI);
   double one_r2 = 1. - r * r;
   double dx = x - x0;
   double dy = y - y0;
   return -0.5 / one_r2 * (pow(dx / sx, 2) + pow(dy / sy, 2) - 2. * r * dx / sx * dy / sy) - c2gaus -
          log(sx * sy * sqrt(one_r2));
 }

◆ MuonResidualsFitter_logROOTVoigt()

double MuonResidualsFitter_logROOTVoigt	(	double	residual,
		double	center,
		double	sigma,
		double	gamma
	)

Definition at line 120 of file MuonResidualsFitter.cc.

References CustomPhysics_cfi::gamma, and dqm-mbProfile::log.

Referenced by MuonResiduals1DOFFitter_FCN(), MuonResiduals5DOFFitter_FCN(), MuonResiduals6DOFFitter_FCN(), MuonResiduals6DOFrphiFitter_FCN(), MuonResidualsAngleFitter_FCN(), MuonResidualsBfieldAngleFitter_FCN(), and MuonResidualsPositionFitter_FCN().

                                                                                                     {
   return log(TMath::Voigt(residual - center, fabs(sigma), fabs(gamma) * 2.));
 }

◆ MuonResidualsFitter_powerLawTails_TF1()

Double_t MuonResidualsFitter_powerLawTails_TF1	(	Double_t *	xvec,
		Double_t *	par
	)

Definition at line 116 of file MuonResidualsFitter.cc.

References JetChargeProducer_cfi::exp, and MuonResidualsFitter_logPowerLawTails().

Referenced by MuonResiduals1DOFFitter::plot(), MuonResidualsPositionFitter::plot(), MuonResidualsAngleFitter::plot(), MuonResidualsBfieldAngleFitter::plot(), MuonResiduals6DOFrphiFitter::plot(), MuonResiduals5DOFFitter::plot(), and MuonResiduals6DOFFitter::plot().

                                                                               {
   return par[0] * exp(MuonResidualsFitter_logPowerLawTails(xvec[0], par[1], par[2], par[3]));
 }

◆ MuonResidualsFitter_pureGaussian_TF1()

Double_t MuonResidualsFitter_pureGaussian_TF1	(	Double_t *	xvec,
		Double_t *	par
	)

Definition at line 38 of file MuonResidualsFitter.cc.

References JetChargeProducer_cfi::exp, and MuonResidualsFitter_logPureGaussian().

Referenced by MuonResiduals1DOFFitter::plot(), MuonResidualsPositionFitter::plot(), MuonResidualsAngleFitter::plot(), MuonResidualsBfieldAngleFitter::plot(), MuonResiduals6DOFrphiFitter::plot(), MuonResiduals5DOFFitter::plot(), and MuonResiduals6DOFFitter::plot().

                                                                              {
   return par[0] * exp(MuonResidualsFitter_logPureGaussian(xvec[0], par[1], par[2]));
 }

◆ MuonResidualsFitter_ROOTVoigt_TF1()

Double_t MuonResidualsFitter_ROOTVoigt_TF1	(	Double_t *	xvec,
		Double_t *	par
	)

Definition at line 124 of file MuonResidualsFitter.cc.

Referenced by MuonResiduals1DOFFitter::plot(), MuonResidualsPositionFitter::plot(), MuonResidualsAngleFitter::plot(), MuonResidualsBfieldAngleFitter::plot(), MuonResiduals6DOFrphiFitter::plot(), MuonResiduals5DOFFitter::plot(), and MuonResiduals6DOFFitter::plot().

                                                                           {
   return par[0] * TMath::Voigt(xvec[0] - par[1], fabs(par[2]), fabs(par[3]) * 2.);
 }