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 61 of file SigmaPtDiff.h.

Constructor & Destructor Documentation

SigmaPtDiff()

SigmaPtDiff::SigmaPtDiff ( )

inline

Definition at line 63 of file SigmaPtDiff.h.

References debug_messages_cfi::errors, 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

etaByPoints()

double SigmaPtDiff::etaByPoints ( const double &  eta )

inline

Definition at line 99 of file SigmaPtDiff.h.

References PVValHelper::eta.

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;
  }

squaredDiff()

double SigmaPtDiff::squaredDiff ( const double &  eta )

inline

Definition at line 149 of file SigmaPtDiff.h.

References PVValHelper::eta, 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

SigmaPt SigmaPtDiff::sigmaPt

Definition at line 158 of file SigmaPtDiff.h.

Referenced by SigmaPtDiff(), and squaredDiff().