CSCThrTurnOnFcn Class Reference

#include <CSCThrTurnOnFcn.h>

Inheritance diagram for CSCThrTurnOnFcn:

Public Member Functions
virtual double	operator() (const std::vector< double > &) const
	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...
virtual double	Up () const
	@ 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 18 of file CSCThrTurnOnFcn.h.

Member Function Documentation

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

virtual

Provide the chi-squared function for the given data.

Definition at line 9 of file CSCThrTurnOnFcn.cc.

References beam_dqm_sourceclient-live_cfg::chi2, mps_update::diff, personalPlayback::fn, i, LogDebug, N, findQualityFiles::size, x(), and detailsBasic3DVector::y.

                                     {
  
  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 29 of file CSCThrTurnOnFcn.h.

References i, xdata, and ydata.

                                                                       { 
    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 36 of file CSCThrTurnOnFcn.h.

References ery, and i.

                                               {
   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 41 of file CSCThrTurnOnFcn.h.

References gen::n, and norm.

{norm=n;};

virtual double CSCThrTurnOnFcn::Up ( ) const

inlinevirtual

@ What?

Definition at line 47 of file CSCThrTurnOnFcn.h.

{return 1.;}

Member Data Documentation

std::vector<float> CSCThrTurnOnFcn::ery

private

Definition at line 24 of file CSCThrTurnOnFcn.h.

Referenced by setErrors().

float CSCThrTurnOnFcn::norm

private

Definition at line 25 of file CSCThrTurnOnFcn.h.

Referenced by setNorm().

std::vector<float> CSCThrTurnOnFcn::xdata

private

data

Definition at line 22 of file CSCThrTurnOnFcn.h.

Referenced by setData().

std::vector<float> CSCThrTurnOnFcn::ydata

private

Definition at line 23 of file CSCThrTurnOnFcn.h.

Referenced by setData().