TopQuarkAnalysis/TopObjectResolutions/interface/Electron.h - master - githubmirror/cmssw

TopQuarkAnalysis/TopObjectResolutions/interface/Electron.h
235 строк8 KB
Shahzad Malik Muzaffar
[ANALYSIS] Changes suggested by new llvm18 clang-format
05 сен 2024, 17:04
05 сен 2024, 17:04a467cd9
#ifndef TopObjetcResolutionsElectron_h
#define TopObjetcResolutionsElectron_h

#include <cmath>

namespace res {
  class HelperElectron {
  public:
    HelperElectron() {}
    ~HelperElectron() {}

    inline double pt(double pt, double eta);
    inline double eta(double pt, double eta);
    inline double phi(double pt, double eta);
    inline double a(double pt, double eta);
    inline double b(double pt, double eta);
    inline double c(double pt, double eta);
    inline double d(double pt, double eta);
    inline double et(double pt, double eta);
    inline double theta(double pt, double eta);
  };
}  // namespace res

inline double res::HelperElectron::pt(double pt, double eta) {
  double res = 0.2 * sqrt(pt);
  if (fabs(eta) <= 0.8)
    res += -0.28 * fabs(eta) + 0.54;
  else if (0.8 < fabs(eta) && fabs(eta) <= 1.4)
    res += 1.52 * fabs(eta) - 1.07;
  else
    res += -0.158 * eta * eta + 0.97;
  return res;
}

inline double res::HelperElectron::a(double pt, double eta) {
  double res = 1000.;
  if (fabs(eta) < 0.17)
    res = 0.0114228 + 0.1135 * exp(-(0.111787 * pt));
  else if (fabs(eta) < 0.35)
    res = 0.010297 + 0.0163219 * exp(-(0.0248655 * pt));
  else if (fabs(eta) < 0.5)
    res = 0.0081634 + 0.0230718 * exp(-(0.0250963 * pt));
  else if (fabs(eta) < 0.7)
    res = 0.0122627 + 0.0228103 * exp(-(0.0420525 * pt));
  else if (fabs(eta) < 0.9)
    res = 0.0124118 + 0.0384051 * exp(-(0.0545988 * pt));
  else if (fabs(eta) < 1.15)
    res = 0.0133398 + 0.0307574 * exp(-(0.0316605 * pt));
  else if (fabs(eta) < 1.3)
    res = 0.0144956 + 0.0355239 * exp(-(0.0273916 * pt));
  else if (fabs(eta) < 1.6)
    res = -13.9017 + 13.9315 * exp(-(7.01823e-06 * pt));
  else if (fabs(eta) < 1.9)
    res = 0.0106309 + 0.0230149 * exp(-(0.00792017 * pt));
  else if (fabs(eta) < 2.5)
    res = 0.0187822 + 0.119922 * exp(-(0.140598 * pt));
  return res;
}

inline double res::HelperElectron::b(double pt, double eta) {
  double res = 1000.;
  if (fabs(eta) < 0.17)
    res = -36.5147 + 36.516 * exp(-(-9.36847e-06 * pt));
  else if (fabs(eta) < 0.35)
    res = -0.132294 + 0.134289 * exp(-(-0.0021212 * pt));
  else if (fabs(eta) < 0.5)
    res = -96.9689 + 96.9704 * exp(-(-2.83158e-06 * pt));
  else if (fabs(eta) < 0.7)
    res = -122.755 + 122.757 * exp(-(-1.90859e-06 * pt));
  else if (fabs(eta) < 0.9)
    res = -226.455 + 226.457 * exp(-(-8.91395e-07 * pt));
  else if (fabs(eta) < 1.15)
    res = -249.279 + 249.281 * exp(-(-7.57645e-07 * pt));
  else if (fabs(eta) < 1.3)
    res = -16.7465 + 16.7481 * exp(-(-1.17848e-05 * pt));
  else if (fabs(eta) < 1.6)
    res = -128.535 + 128.537 * exp(-(-1.52236e-06 * pt));
  else if (fabs(eta) < 1.9)
    res = -66.1731 + 66.1762 * exp(-(-2.49121e-06 * pt));
  else if (fabs(eta) < 2.5)
    res = -0.0199509 + 0.0237796 * exp(-(-0.00591733 * pt));
  return res;
}

inline double res::HelperElectron::c(double pt, double eta) {
  double res = 1000.;
  if (fabs(eta) < 0.17)
    res = -68.6101 + 68.6161 * exp(-(-7.88064e-07 * pt));
  else if (fabs(eta) < 0.35)
    res = -130.219 + 130.224 * exp(-(-8.39696e-07 * pt));
  else if (fabs(eta) < 0.5)
    res = -262.825 + 262.83 * exp(-(-4.26008e-07 * pt));
  else if (fabs(eta) < 0.7)
    res = -47.898 + 47.9041 * exp(-(-2.07568e-06 * pt));
  else if (fabs(eta) < 0.9)
    res = -178.795 + 178.799 * exp(-(-5.78263e-07 * pt));
  else if (fabs(eta) < 1.15)
    res = -249.26 + 249.267 * exp(-(-3.16408e-07 * pt));
  else if (fabs(eta) < 1.3)
    res = -173.603 + 173.606 * exp(-(-1.58982e-06 * pt));
  else if (fabs(eta) < 1.6)
    res = -98.726 + 98.7326 * exp(-(-2.67151e-06 * pt));
  else if (fabs(eta) < 1.9)
    res = -82.6028 + 82.6122 * exp(-(-2.79483e-06 * pt));
  else if (fabs(eta) < 2.5)
    res = -119.94 + 119.95 * exp(-(-1.69882e-06 * pt));
  return res;
}

inline double res::HelperElectron::d(double pt, double eta) {
  double res = 1000.;
  if (fabs(eta) < 0.17)
    res = 0.0114197 + 0.081753 * exp(-(0.0969625 * pt));
  else if (fabs(eta) < 0.35)
    res = 0.0097562 + 0.0137737 * exp(-(0.0187112 * pt));
  else if (fabs(eta) < 0.5)
    res = 0.00844899 + 0.021595 * exp(-(0.0241679 * pt));
  else if (fabs(eta) < 0.7)
    res = 0.0124758 + 0.0243678 * exp(-(0.0459914 * pt));
  else if (fabs(eta) < 0.9)
    res = 0.0117518 + 0.0357218 * exp(-(0.0463595 * pt));
  else if (fabs(eta) < 1.15)
    res = 0.0133048 + 0.0292511 * exp(-(0.0305812 * pt));
  else if (fabs(eta) < 1.3)
    res = 0.0144949 + 0.0337369 * exp(-(0.0269872 * pt));
  else if (fabs(eta) < 1.6)
    res = 0.0137634 + 0.0343427 * exp(-(0.0208025 * pt));
  else if (fabs(eta) < 1.9)
    res = 0.00572644 + 0.0270719 * exp(-(0.00536847 * pt));
  else if (fabs(eta) < 2.5)
    res = 0.0189177 + 0.169591 * exp(-(0.152597 * pt));
  return res;
}

inline double res::HelperElectron::theta(double pt, double eta) {
  double res = 1000.;
  if (fabs(eta) < 0.17)
    res = 0.000282805 + 0.000157786 * exp(-(0.0343273 * pt));
  else if (fabs(eta) < 0.35)
    res = 0.000184362 + 4.34076e-05 * exp(-(-0.0131909 * pt));
  else if (fabs(eta) < 0.5)
    res = 0.000249332 + 5.83114e-05 * exp(-(0.0508729 * pt));
  else if (fabs(eta) < 0.7)
    res = -6.56357e-05 + 0.000325051 * exp(-(0.00177319 * pt));
  else if (fabs(eta) < 0.9)
    res = 0.000182277 + 0.000125324 * exp(-(0.0581923 * pt));
  else if (fabs(eta) < 1.15)
    res = 0.000140771 + 0.000407914 * exp(-(0.0971668 * pt));
  else if (fabs(eta) < 1.3)
    res = 0.000125551 + 0.001266 * exp(-(0.180176 * pt));
  else if (fabs(eta) < 1.6)
    res = 0.000107631 + 101920 * exp(-(1.17024 * pt));
  else if (fabs(eta) < 1.9)
    res = 8.33927e-05 + 158936 * exp(-(1.20127 * pt));
  else if (fabs(eta) < 2.5)
    res = 6.55271e-05 + 0.12459 * exp(-(0.437044 * pt));
  return res;
}

inline double res::HelperElectron::phi(double pt, double eta) {
  double res = 1000.;
  if (fabs(eta) < 0.17)
    res = 0.000175676 + 0.000471783 * exp(-(0.0383161 * pt));
  else if (fabs(eta) < 0.35)
    res = 0.000202185 + 0.00048635 * exp(-(0.0373331 * pt));
  else if (fabs(eta) < 0.5)
    res = 0.000150868 + 0.000444216 * exp(-(0.0268835 * pt));
  else if (fabs(eta) < 0.7)
    res = 0.000243624 + 0.00182347 * exp(-(0.0850746 * pt));
  else if (fabs(eta) < 0.9)
    res = 0.000254463 + 0.000431233 * exp(-(0.0446507 * pt));
  else if (fabs(eta) < 1.15)
    res = 0.000309592 + 0.000918965 * exp(-(0.0555677 * pt));
  else if (fabs(eta) < 1.3)
    res = 0.000502204 + 0.000277996 * exp(-(0.076721 * pt));
  else if (fabs(eta) < 1.6)
    res = 0.000361181 + 0.000655126 * exp(-(0.0238519 * pt));
  else if (fabs(eta) < 1.9)
    res = 0.000321587 + 0.00155721 * exp(-(0.0337709 * pt));
  else if (fabs(eta) < 2.5)
    res = 0.000819101 + 0.00205336 * exp(-(0.0992806 * pt));
  return res;
}

inline double res::HelperElectron::et(double pt, double eta) {
  double res = 1000.;
  if (fabs(eta) < 0.17)
    res = 0.326238 + 0.00760789 * pt;
  else if (fabs(eta) < 0.35)
    res = 0.40493 + 0.00659958 * pt;
  else if (fabs(eta) < 0.5)
    res = 0.369785 + 0.00690331 * pt;
  else if (fabs(eta) < 0.7)
    res = 0.437539 + 0.00703808 * pt;
  else if (fabs(eta) < 0.9)
    res = 0.456138 + 0.0078252 * pt;
  else if (fabs(eta) < 1.15)
    res = 0.518685 + 0.00907836 * pt;
  else if (fabs(eta) < 1.3)
    res = 0.733672 + 0.00953255 * pt;
  else if (fabs(eta) < 1.6)
    res = 1.02678 + 0.0116056 * pt;
  else if (fabs(eta) < 1.9)
    res = 0.948368 + 0.00977619 * pt;
  else if (fabs(eta) < 2.5)
    res = 0.418302 + 0.0127816 * pt;
  return res;
}

inline double res::HelperElectron::eta(double pt, double eta) {
  double res = 1000.;
  if (fabs(eta) < 0.17)
    res = 0.000266154 + 0.000104322 * exp(-(0.0140464 * pt));
  else if (fabs(eta) < 0.35)
    res = -0.251539 + 0.251791 * exp(-(-7.37147e-07 * pt));
  else if (fabs(eta) < 0.5)
    res = 0.000290074 + 1.54664e-12 * exp(-(-0.115541 * pt));
  else if (fabs(eta) < 0.7)
    res = 8.37182e-05 + 0.000233453 * exp(-(0.00602386 * pt));
  else if (fabs(eta) < 0.9)
    res = 0.000229422 + 0.000114253 * exp(-(0.0188935 * pt));
  else if (fabs(eta) < 1.15)
    res = 0.000191525 + 0.000404238 * exp(-(0.0554545 * pt));
  else if (fabs(eta) < 1.3)
    res = 0.000195461 + 1.43699e-07 * exp(-(-0.0315088 * pt));
  else if (fabs(eta) < 1.6)
    res = 0.000223422 + 2.05169e+07 * exp(-(1.41408 * pt));
  else if (fabs(eta) < 1.9)
    res = -0.867114 + 0.867336 * exp(-(4.08639e-07 * pt));
  else if (fabs(eta) < 2.5)
    res = -0.866567 + 0.866746 * exp(-(-8.62871e-07 * pt));
  return res;
}

#endif