SigmaPtDiff Class Reference

Returns ( sigmaPt/Pt(data)^2 - sigmaPt/Pt(MC)^2 ) More...

#include <SigmaPtDiff.h>

Public Member Functions
double	etaByPoints (const double &eta)
	SigmaPtDiff ()
double	squaredDiff (const double &eta)

Public Attributes
SigmaPt	sigmaPt

Detailed Description

Returns ( sigmaPt/Pt(data)^2 - sigmaPt/Pt(MC)^2 )

Definition at line 68 of file SigmaPtDiff.h.

Constructor & Destructor Documentation

SigmaPtDiff::SigmaPtDiff ( )

inline

Definition at line 71 of file SigmaPtDiff.h.

References benchmark_cfg::errors, Parameters::parameters, SigmaPt::setParErr(), and sigmaPt.

  {
    std::vector<double> parameters;
    std::vector<double> errors;
    parameters.push_back(1.66);
    parameters.push_back(0.021);
    parameters.push_back(0.);
    parameters.push_back(0.);
    parameters.push_back(0.);
    parameters.push_back(0.0058);
    parameters.push_back(0.);
    parameters.push_back(0.03);
    parameters.push_back(1.8);
    parameters.push_back(0.);
    parameters.push_back(0.);
    parameters.push_back(0.);
    parameters.push_back(0.);
    parameters.push_back(0.);
    parameters.push_back(0.);
    errors.push_back(0.09);
    errors.push_back(0.002);
    errors.push_back(0.);
    errors.push_back(0.);
    errors.push_back(0.);
    errors.push_back(0.0009);
    errors.push_back(0.);
    errors.push_back(0.03);
    errors.push_back(0.3);
    errors.push_back(0.);
    errors.push_back(0.);
    errors.push_back(0.);
    errors.push_back(0.);
    errors.push_back(0.);
    errors.push_back(0.);

    sigmaPt.setParErr( parameters, errors );
  }

Member Function Documentation

double SigmaPtDiff::etaByPoints ( const double &  eta )

inline

Definition at line 108 of file SigmaPtDiff.h.

Referenced by squaredDiff().

  {
    if(      eta <= -2.2 ) return 0.0233989;
    else if( eta <= -2.0 ) return 0.0197057;
    else if( eta <= -1.8 ) return 0.014693;
    else if( eta <= -1.6 ) return 0.0146727;
    else if( eta <= -1.4 ) return 0.0141323;
    else if( eta <= -1.2 ) return 0.0159712;
    else if( eta <= -1.0 ) return 0.0117224;
    else if( eta <= -0.8 ) return 0.010726;
    else if( eta <= -0.6 ) return 0.0104777;
    else if( eta <= -0.4 ) return 0.00814458;
    else if( eta <= -0.2 ) return 0.00632501;
    else if( eta <=  0.0 ) return 0.00644172;
    else if( eta <=  0.2 ) return 0.00772645;
    else if( eta <=  0.4 ) return 0.010103;
    else if( eta <=  0.6 ) return 0.0099275;
    else if( eta <=  0.8 ) return 0.0100309;
    else if( eta <=  1.0 ) return 0.0125116;
    else if( eta <=  1.2 ) return 0.0147211;
    else if( eta <=  1.4 ) return 0.0151623;
    else if( eta <=  1.6 ) return 0.015259;
    else if( eta <=  1.8 ) return 0.014499;
    else if( eta <=  2.0 ) return 0.0165215;
    else if( eta <=  2.2 ) return 0.0212348;
    return 0.0227285;
  }

double SigmaPtDiff::squaredDiff ( const double &  eta )

inline

Definition at line 136 of file SigmaPtDiff.h.

References etaByPoints(), SigmaPt::sigma(), and sigmaPt.

Referenced by ErrorsPropagationAnalyzer::fillHistograms(), and smearFunctionType7::smear().

  {
    double sigmaPtPlus = sigmaPt(eta) + sigmaPt.sigma(eta);
    double sigmaPtMinus = sigmaPt(eta) - sigmaPt.sigma(eta);
    if( fabs(sigmaPtPlus*sigmaPtPlus - etaByPoints(eta)*etaByPoints(eta)) > fabs(sigmaPtMinus*sigmaPtMinus - etaByPoints(eta)*etaByPoints(eta)) ) {
      return( fabs(sigmaPtPlus*sigmaPtPlus - etaByPoints(eta)*etaByPoints(eta)) );
    }
    return( fabs(sigmaPtMinus*sigmaPtMinus - etaByPoints(eta)*etaByPoints(eta)) );
  }

Member Data Documentation

SigmaPt SigmaPtDiff::sigmaPt

Definition at line 145 of file SigmaPtDiff.h.

Referenced by SigmaPtDiff(), and squaredDiff().