CSCThrTurnOnFcn Class Reference

#include <CSCThrTurnOnFcn.h>

Inheritance diagram for CSCThrTurnOnFcn:

Public Member Functions
double	operator() (const std::vector< double > &) const override
	Provide the chi-squared function for the given data. More...
void	setData (const std::vector< float > &x, const std::vector< float > &y)
	Cache the current data, x and y. More...
void	setErrors (const std::vector< float > &er)
	Set the errors. More...
void	setNorm (float n)
	Set the norm (if needed) More...
double	Up () const override
	@ What? More...

Private Attributes
std::vector< float >	ery
float	norm
std::vector< float >	xdata
	data More...
std::vector< float >	ydata

Detailed Description

Model functional form for fitting AFEB turn-on threshold information from Muon Endcap CSC's. This version is for ROOT Minuit2. Based on CSCPulseHeightFcn as an example.

Definition at line 15 of file CSCThrTurnOnFcn.h.

Member Function Documentation

double CSCThrTurnOnFcn::operator() ( const std::vector< double > &  par ) const

override

Provide the chi-squared function for the given data.

Definition at line 8 of file CSCThrTurnOnFcn.cc.

References hltPixelTracks_cff::chi2, change_name::diff, ery, personalPlayback::fn, mps_fire::i, LogDebug, N, norm, findQualityFiles::size, x, xdata, y, and ydata.

Referenced by setNorm().

                                                                     {
  double x, y, er, fn;
  double N = norm;
  double chi2 = 0.;

  int size = xdata.size();
  for (int i = 0; i < size; ++i) {
    x = xdata[i];
    y = ydata[i];
    er = ery[i];

    fn = (x - par[0]) / (par[1] * 1.4142);
    fn = N * (1. - TMath::Erf(fn)) / 2.;

    double diff = y - fn;
    chi2 += diff * diff / (er * er);

    //    std::cout<<"CSC AFEB threshold fit "<<i+1<<" "<<x<<" "<<y<<" "
    //                   <<er<<" "<<fn<<" "<<chi2<<" "
    //                   <<par[0]<<" "<<par[1]<<"\n";

    LogDebug("CSC") << " AFEB threshold fit " << i + 1 << " " << x << " " << y << " " << er << " " << fn << " " << chi2
                    << " " << par[0] << " " << par[1] << "\n";
  }
  //  std::cout<<"Chi2 "<<chi2<<std::endl;
  return chi2;
}

void CSCThrTurnOnFcn::setData ( const std::vector< float > &  x, const std::vector< float > &  y  )

inline

Cache the current data, x and y.

Definition at line 25 of file CSCThrTurnOnFcn.h.

References mps_fire::i.

Referenced by CSCFitAFEBThr::ThresholdNoise().

                                                                     {
    for (unsigned int i = 0; i < x.size(); i++) {
      xdata.push_back(x[i]);
      ydata.push_back(y[i]);
    }
  };

void CSCThrTurnOnFcn::setErrors ( const std::vector< float > &  er )

inline

Set the errors.

Definition at line 33 of file CSCThrTurnOnFcn.h.

References mps_fire::i.

Referenced by CSCFitAFEBThr::ThresholdNoise().

                                             {
    for (unsigned int i = 0; i < er.size(); i++)
      ery.push_back(er[i]);
  };

void CSCThrTurnOnFcn::setNorm ( float  n )

inline

Set the norm (if needed)

Definition at line 39 of file CSCThrTurnOnFcn.h.

References dqmiodumpmetadata::n, and operator()().

Referenced by CSCFitAFEBThr::ThresholdNoise().

{ norm = n; };

double CSCThrTurnOnFcn::Up ( ) const

inlineoverride

@ What?

Definition at line 45 of file CSCThrTurnOnFcn.h.

{ return 1.; }

Member Data Documentation

std::vector<float> CSCThrTurnOnFcn::ery

private

Definition at line 20 of file CSCThrTurnOnFcn.h.

Referenced by operator()().

float CSCThrTurnOnFcn::norm

private

Definition at line 21 of file CSCThrTurnOnFcn.h.

Referenced by operator()().

std::vector<float> CSCThrTurnOnFcn::xdata

private

data

Definition at line 18 of file CSCThrTurnOnFcn.h.

Referenced by operator()().

std::vector<float> CSCThrTurnOnFcn::ydata

private

Definition at line 19 of file CSCThrTurnOnFcn.h.

Referenced by operator()().