CSCThrTurnOnFcn Class Reference

#include <CSCThrTurnOnFcn.h>

Inheritance diagram for CSCThrTurnOnFcn:

Public Member Functions
double	operator() (std::span< const 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 19 of file CSCThrTurnOnFcn.h.

Member Function Documentation

operator()()

double CSCThrTurnOnFcn::operator() ( std::span< const double >  par ) const

override

Provide the chi-squared function for the given data.

Definition at line 9 of file CSCThrTurnOnFcn.cc.

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

                                                                  {
#else
double CSCThrTurnOnFcn::operator()(const std::vector<double>& par) const {
#endif
  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;
}

setData()

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

inline

Cache the current data, x and y.

Definition at line 29 of file CSCThrTurnOnFcn.h.

References mps_fire::i, x, xdata, y, and ydata.

Referenced by CSCFitAFEBThr::ThresholdNoise().

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

setErrors()

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

inline

Set the errors.

Definition at line 37 of file CSCThrTurnOnFcn.h.

References ery, and mps_fire::i.

Referenced by CSCFitAFEBThr::ThresholdNoise().

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

setNorm()

void CSCThrTurnOnFcn::setNorm ( float  n )

inline

Set the norm (if needed)

Definition at line 43 of file CSCThrTurnOnFcn.h.

References create_idmaps::n, and norm.

Referenced by CSCFitAFEBThr::ThresholdNoise().

{ norm = n; };

Up()

double CSCThrTurnOnFcn::Up ( ) const

inlineoverride

@ What?

Definition at line 53 of file CSCThrTurnOnFcn.h.

{ return 1.; }

Member Data Documentation

ery

std::vector<float> CSCThrTurnOnFcn::ery

private

Definition at line 24 of file CSCThrTurnOnFcn.h.

Referenced by operator()(), and setErrors().

norm

float CSCThrTurnOnFcn::norm

private

Definition at line 25 of file CSCThrTurnOnFcn.h.

Referenced by operator()(), and setNorm().

xdata

std::vector<float> CSCThrTurnOnFcn::xdata

private

data

Definition at line 22 of file CSCThrTurnOnFcn.h.

Referenced by operator()(), and setData().

ydata

std::vector<float> CSCThrTurnOnFcn::ydata

private

Definition at line 23 of file CSCThrTurnOnFcn.h.

Referenced by operator()(), and setData().