hitfit::Top_Decaykin Class Reference

A class to hold functions to calculate kinematic quantities of interest in events. This class has no state, only static member functions. More...

#include <Top_Decaykin.h>

Static Public Member Functions
static double	cos_theta_star (const Fourvec &fermion, const Fourvec &W, const Fourvec &top)
	Calculate in top quark decay. More...
static double	cos_theta_star (const Lepjets_Event &ev)
	Calculate the lepton in top quark leptonic decay. More...
static double	cos_theta_star_hadt (const Lepjets_Event &ev)
	Calculate the hadronic in top quark leptonic decay. As there is no information on the weak isospin component of the fermion, the absolute value of will be returned (the solutions for up-type and down-type fermions will differ only in sign but not in magnitude). More...
static double	cos_theta_star_lept (const Lepjets_Event &ev)
	Calculate the lepton in top quark leptonic decay. More...
static std::ostream &	dump_ev (std::ostream &s, const Lepjets_Event &ev)
	Print the kinematic information for an event. More...
static Fourvec	hadt (const Lepjets_Event &ev)
	Sum up the appropriate four-momenta to find the hadronic top quark. More...
static Fourvec	hadw (const Lepjets_Event &ev)
	Sum up the appropriate four-momenta to find the hadronic boson. More...
static Fourvec	hadw1 (const Lepjets_Event &ev)
	Return the hadronic boson jet which have higher . More...
static Fourvec	hadw2 (const Lepjets_Event &ev)
	Return the hadronic boson jet which have lower . More...
static Fourvec	lept (const Lepjets_Event &ev)
	Sum up the appropriate four-momenta to find the leptonic top quark. More...
static Fourvec	lepw (const Lepjets_Event &ev)
	Sum up the appropriate four-momenta to find the leptonic boson. More...
static bool	solve_nu (const Lepjets_Event &ev, double wmass, double &nuz1, double &nuz2)
	Solve for the longitudinal momentum that makes the leptonic -boson to have a certain value of mass. Returns TRUE if there were real solutions. Returns FALSE if there were only complex solutions. In case of complex solutions, the real components of the solutions are given. More...
static bool	solve_nu (const Lepjets_Event &ev, double wmass, double &re_nuz1, double &im_nuz1, double &re_nuz2, double &im_nuz2)
	Solve for the longitudinal momentum that makes the leptonic -boson to have a certain value of mass. The complex component of the solutions are also given. Returns TRUE if there were real solutions. Returns FALSE if there were only complex solutions. In case of real solutions, the first solution is the one which have smaller absolute value. In case of imaginary solutions (which are complex conjugate of each other), the first solution is the one which have imaginary component in the lower half of the complex plane, i.e., the one which have negative imaginary component). More...
static bool	solve_nu_tmass (const Lepjets_Event &ev, double tmass, double &nuz1, double &nuz2)
	Solve for the neutrino longitudinal momentum that makes the leptonic top have a certain value of mass. Returns TRUE if there were real solutions. Returns FALSE if there were only complex solutions. In case of complex solutions, the real components of the solutions are given. More...
static bool	solve_nu_tmass (const Lepjets_Event &ev, double tmass, double &re_nuz1, double &im_nuz1, double &re_nuz2, double &im_nuz2)
	Solve for the neutrino longitudinal momentum that makes the leptonic top have a certain value of mass. The complex component of the solutions are also given. Returns TRUE if there were real solutions. Returns FALSE if there were only complex solutions. In case of real solutions, the first solution is the one which have smaller absolute value. In case of imaginary solutions (which are complex conjugate of each other), the first solution is the one which have imaginary component in the lower half of the complex plane, i.e., the one which have negative imaginary component). More...

Detailed Description

A class to hold functions to calculate kinematic quantities of interest in events. This class has no state, only static member functions.

Definition at line 51 of file Top_Decaykin.h.

Member Function Documentation

cos_theta_star() [1/2]

double hitfit::Top_Decaykin::cos_theta_star ( const Fourvec &  fermion, const Fourvec &  W, const Fourvec &  top  )

static

Calculate in top quark decay.

Parameters

fermion	The four-momentum of fermion from boson from top decay.
W	The four-momentum of boson from top decay.
top	The four-momentum of top.

Definition at line 394 of file Top_Decaykin.cc.

Referenced by cos_theta_star(), cos_theta_star_hadt(), and cos_theta_star_lept().

  {
    if (W.isLightlike() || W.isSpacelike()) {
      return 100.0;
    }

    CLHEP::HepBoost BoostWCM(W.findBoostToCM());

    CLHEP::Hep3Vector boost_v3fermion = BoostWCM(fermion).vect();
    CLHEP::Hep3Vector boost_v3top = BoostWCM(top).vect();

    double costhetastar = boost_v3fermion.cosTheta(-boost_v3top);

    return costhetastar;
  }

cos_theta_star() [2/2]

double hitfit::Top_Decaykin::cos_theta_star ( const Lepjets_Event &  ev )

static

Calculate the lepton in top quark leptonic decay.

Parameters

ev	The event to solve.

Definition at line 420 of file Top_Decaykin.cc.

References cos_theta_star(), makeMEIFBenchmarkPlots::ev, lept(), HLT_2022v15_cff::leptons, and lepw().

  {
    return cos_theta_star(leptons(ev), lepw(ev), lept(ev));
  }

cos_theta_star_hadt()

double hitfit::Top_Decaykin::cos_theta_star_hadt ( const Lepjets_Event &  ev )

static

Calculate the hadronic in top quark leptonic decay. As there is no information on the weak isospin component of the fermion, the absolute value of will be returned (the solutions for up-type and down-type fermions will differ only in sign but not in magnitude).

Parameters

ev	The event to solve.

Definition at line 446 of file Top_Decaykin.cc.

References cos_theta_star(), makeMEIFBenchmarkPlots::ev, hadt(), hadw(), and hadw1().

  {
    return fabs(cos_theta_star(hadw1(ev), hadw(ev), hadt(ev)));
  }

cos_theta_star_lept()

double hitfit::Top_Decaykin::cos_theta_star_lept ( const Lepjets_Event &  ev )

static

Calculate the lepton in top quark leptonic decay.

Parameters

ev	The event to solve.

Definition at line 433 of file Top_Decaykin.cc.

References cos_theta_star(), and makeMEIFBenchmarkPlots::ev.

  {
    return cos_theta_star(ev);
  }

dump_ev()

ostream & hitfit::Top_Decaykin::dump_ev ( std::ostream &  s, const Lepjets_Event &  ev  )

static

Print the kinematic information for an event.

Parameters

s	The stream of which to write.
ev	The event to be printed.

Definition at line 367 of file Top_Decaykin.cc.

References makeMEIFBenchmarkPlots::ev, hadt(), hadw(), lept(), lepw(), AlCaHLTBitMon_ParallelJobs::p, and alignCSCRings::s.

Referenced by hitfit::Top_Fit::fit_one_perm().

  {
    s << ev;
    Fourvec p;

    p = lepw(ev);
    s << "lepw " << p << " " << p.m() << "\n";
    p = lept(ev);
    s << "lept " << p << " " << p.m() << "\n";
    p = hadw(ev);
    s << "hadw " << p << " " << p.m() << "\n";
    p = hadt(ev);
    s << "hadt " << p << " " << p.m() << "\n";

    return s;
  }

hadt()

Fourvec hitfit::Top_Decaykin::hadt ( const Lepjets_Event &  ev )

static

Sum up the appropriate four-momenta to find the hadronic top quark.

Parameters

ev	The event.

Definition at line 339 of file Top_Decaykin.cc.

References makeMEIFBenchmarkPlots::ev, hitfit::hadb_label, and hadw().

Referenced by cos_theta_star_hadt(), dump_ev(), and hitfit::Top_Fit::fit_one_perm().

  {
    return (ev.sum(hadb_label) + hadw(ev));
  }

hadw()

Fourvec hitfit::Top_Decaykin::hadw ( const Lepjets_Event &  ev )

static

Sum up the appropriate four-momenta to find the hadronic boson.

Parameters

ev	The event.

Definition at line 282 of file Top_Decaykin.cc.

References makeMEIFBenchmarkPlots::ev, hitfit::hadw1_label, and hitfit::hadw2_label.

Referenced by cos_theta_star_hadt(), dump_ev(), hitfit::Top_Fit::fit_one_perm(), and hadt().

  {
    return (ev.sum(hadw1_label) + ev.sum(hadw2_label));
  }

hadw1()

Fourvec hitfit::Top_Decaykin::hadw1 ( const Lepjets_Event &  ev )

static

Return the hadronic boson jet which have higher .

Parameters

ev	The event.

Definition at line 296 of file Top_Decaykin.cc.

References makeMEIFBenchmarkPlots::ev, and hitfit::hadw1_label.

Referenced by cos_theta_star_hadt().

  {
    return ev.sum(hadw1_label);
  }

hadw2()

Fourvec hitfit::Top_Decaykin::hadw2 ( const Lepjets_Event &  ev )

static

Return the hadronic boson jet which have lower .

Parameters

ev	The event.

Definition at line 310 of file Top_Decaykin.cc.

References makeMEIFBenchmarkPlots::ev, and hitfit::hadw2_label.

  {
    return ev.sum(hadw2_label);
  }

lept()

Fourvec hitfit::Top_Decaykin::lept ( const Lepjets_Event &  ev )

static

Sum up the appropriate four-momenta to find the leptonic top quark.

Parameters

ev	The event.

Definition at line 353 of file Top_Decaykin.cc.

References makeMEIFBenchmarkPlots::ev, hitfit::lepb_label, and lepw().

Referenced by cos_theta_star(), dump_ev(), and hitfit::Top_Fit::fit_one_perm().

  {
    return (ev.sum(lepb_label) + lepw(ev));
  }

lepw()

Fourvec hitfit::Top_Decaykin::lepw ( const Lepjets_Event &  ev )

static

Sum up the appropriate four-momenta to find the leptonic boson.

Parameters

ev	The event.

Definition at line 325 of file Top_Decaykin.cc.

References makeMEIFBenchmarkPlots::ev, and HLT_2022v15_cff::leptons.

Referenced by cos_theta_star(), dump_ev(), and lept().

  {
    return (leptons(ev) + ev.met());
  }

solve_nu() [1/2]

bool hitfit::Top_Decaykin::solve_nu ( const Lepjets_Event &  ev, double  wmass, double &  nuz1, double &  nuz2  )

static

Solve for the longitudinal momentum that makes the leptonic -boson to have a certain value of mass. Returns TRUE if there were real solutions. Returns FALSE if there were only complex solutions. In case of complex solutions, the real components of the solutions are given.

Parameters

ev	Input: The event to solve.
wmass	Input: The desired mass of the boson in GeV.
nuz1	Output: First solution (smaller absolute value).
nuz2	Output: Second solution.

Definition at line 177 of file Top_Decaykin.cc.

References a, funct::abs(), b, c, ztail::d, makeMEIFBenchmarkPlots::ev, HLT_2022v15_cff::leptons, mathSSE::sqrt(), edm::swap(), and x.

Referenced by hitfit::Top_Fit::fit_one_perm().

  {
    bool discrim_flag = true;

    const Fourvec& vnu = ev.met();
    Fourvec vlep = leptons(ev);

    double x = vlep.x() * vnu.x() + vlep.y() * vnu.y() + wmass * wmass / 2;
    double a = vlep.z() * vlep.z() - vlep.e() * vlep.e();
    double b = 2 * x * vlep.z();
    double c = x * x - vnu.perp2() * vlep.e() * vlep.e();

    double d = b * b - 4 * a * c;
    if (d < 0) {
      d = 0;
      discrim_flag = false;
    }

    nuz1 = (-b + sqrt(d)) / 2 / a;
    nuz2 = (-b - sqrt(d)) / 2 / a;
    if (abs(nuz1) > abs(nuz2))
      swap(nuz1, nuz2);

    return discrim_flag;
  }

solve_nu() [2/2]

bool hitfit::Top_Decaykin::solve_nu ( const Lepjets_Event &  ev, double  wmass, double &  re_nuz1, double &  im_nuz1, double &  re_nuz2, double &  im_nuz2  )

static

Solve for the longitudinal momentum that makes the leptonic -boson to have a certain value of mass. The complex component of the solutions are also given. Returns TRUE if there were real solutions. Returns FALSE if there were only complex solutions. In case of real solutions, the first solution is the one which have smaller absolute value. In case of imaginary solutions (which are complex conjugate of each other), the first solution is the one which have imaginary component in the lower half of the complex plane, i.e., the one which have negative imaginary component).

Parameters

ev	Input: The event to solve.
wmass	Input: The desired mass of the boson in GeV.
re_nuz1	Output: Real component of the first solution.
im_nuz1	Output: Imaginary component of the first solution.
re_nuz2	Output: Real component of the second solution.
im_nuz2	Output: Imaginary component of the second solution.

Definition at line 220 of file Top_Decaykin.cc.

References a, b, c, ztail::d, makeMEIFBenchmarkPlots::ev, HLT_2022v15_cff::leptons, mathSSE::sqrt(), edm::swap(), and x.

  {
    bool discrim_flag = true;

    const Fourvec& vnu = ev.met();
    Fourvec vlep = leptons(ev);

    double x = vlep.x() * vnu.x() + vlep.y() * vnu.y() + wmass * wmass / 2;
    double a = vlep.z() * vlep.z() - vlep.e() * vlep.e();
    double b = 2 * x * vlep.z();
    double c = x * x - vnu.perp2() * vlep.e() * vlep.e();

    double d = b * b - 4 * a * c;
    if (d < 0) {
      discrim_flag = false;
    }

    if (discrim_flag) {
      re_nuz1 = (-b + sqrt(d)) / 2 / a;
      im_nuz1 = 0.0;
      re_nuz2 = (-b - sqrt(d)) / 2 / a;
      im_nuz2 = 0.0;
      if (fabs(re_nuz1) > fabs(re_nuz2)) {
        swap(re_nuz1, re_nuz2);
      }

    } else {
      // Take absolute value of the imaginary component of nuz, in case
      // a is negative, before multiplying by +1 or -1 to get the upper-half
      // or lower-half imaginary value.

      re_nuz1 = -b / 2 / a;
      im_nuz1 = -fabs(sqrt(-d) / a);
      re_nuz2 = -b / 2 / a;
      im_nuz2 = fabs(sqrt(-d) / a);
    }

    return discrim_flag;
  }

solve_nu_tmass() [1/2]

bool hitfit::Top_Decaykin::solve_nu_tmass ( const Lepjets_Event &  ev, double  tmass, double &  nuz1, double &  nuz2  )

static

Solve for the neutrino longitudinal momentum that makes the leptonic top have a certain value of mass. Returns TRUE if there were real solutions. Returns FALSE if there were only complex solutions. In case of complex solutions, the real components of the solutions are given.

Parameters

ev	Input:The event to solve.
tmass	Input: The desired value of top quark mass in GeV.
nuz1	Output: The first solution (smaller absolute value).
nuz2	Output: The second solution.

Definition at line 72 of file Top_Decaykin.cc.

References a, funct::abs(), alpha, b, c, ztail::d, createTree::dd, makeMEIFBenchmarkPlots::ev, hitfit::lepb_label, HLT_2022v15_cff::leptons, mathSSE::sqrt(), and edm::swap().

Referenced by hitfit::Top_Fit::fit_one_perm().

  {
    bool discrim_flag = true;

    const Fourvec& vnu = ev.met();
    Fourvec cprime = leptons(ev) + ev.sum(lepb_label);
    double alpha1 = tmass * tmass - cprime.m2();
    double a = 2 * 4 * (cprime.z() * cprime.z() - cprime.e() * cprime.e());
    double alpha = alpha1 + 2 * (cprime.x() * vnu.x() + cprime.y() * vnu.y());
    double b = 4 * alpha * cprime.z();
    double c = alpha * alpha - 4 * cprime.e() * cprime.e() * vnu.vect().perp2();
    double d = b * b - 2 * a * c;
    if (d < 0) {
      discrim_flag = false;
      d = 0;
    }

    double dd = sqrt(d);
    nuz1 = (-b + dd) / a;
    nuz2 = (-b - dd) / a;
    if (abs(nuz1) > abs(nuz2))
      swap(nuz1, nuz2);

    return discrim_flag;
  }

solve_nu_tmass() [2/2]

bool hitfit::Top_Decaykin::solve_nu_tmass ( const Lepjets_Event &  ev, double  tmass, double &  re_nuz1, double &  im_nuz1, double &  re_nuz2, double &  im_nuz2  )

static

Solve for the neutrino longitudinal momentum that makes the leptonic top have a certain value of mass. The complex component of the solutions are also given. Returns TRUE if there were real solutions. Returns FALSE if there were only complex solutions. In case of real solutions, the first solution is the one which have smaller absolute value. In case of imaginary solutions (which are complex conjugate of each other), the first solution is the one which have imaginary component in the lower half of the complex plane, i.e., the one which have negative imaginary component).

Parameters

ev	Input:The event to solve.
tmass	Input: The desired value of top quark mass in GeV.
re_nuz1	Output: Real component of the first solution.
im_nuz1	Output: Imaginary component of the first solution.
re_nuz2	Output: Real component of the second solution.
im_nuz2	Output: Imaginary component of the second solution.

Definition at line 115 of file Top_Decaykin.cc.

References a, funct::abs(), alpha, b, c, ztail::d, makeMEIFBenchmarkPlots::ev, hitfit::lepb_label, HLT_2022v15_cff::leptons, mathSSE::sqrt(), and edm::swap().

  {
    bool discrim_flag = true;

    const Fourvec& vnu = ev.met();
    Fourvec cprime = leptons(ev) + ev.sum(lepb_label);
    double alpha1 = tmass * tmass - cprime.m2();
    double a = 2 * 4 * (cprime.z() * cprime.z() - cprime.e() * cprime.e());
    // Haryo's note: Here a is equivalent to '2a' in the quadratic
    // equation ax^2 + bx + c = 0
    double alpha = alpha1 + 2 * (cprime.x() * vnu.x() + cprime.y() * vnu.y());
    double b = 4 * alpha * cprime.z();
    double c = alpha * alpha - 4 * cprime.e() * cprime.e() * vnu.vect().perp2();
    double d = b * b - 2 * a * c;
    if (d < 0) {
      discrim_flag = false;
    }

    if (discrim_flag) {
      re_nuz1 = (-b + sqrt(d)) / a;
      im_nuz1 = 0;
      re_nuz2 = (-b - sqrt(d)) / a;
      im_nuz2 = 0;
      if (abs(re_nuz1) > abs(re_nuz2)) {
        swap(re_nuz1, re_nuz2);
      }

    } else {
      // Take absolute value of the imaginary component of nuz, in case
      // a is negative, before multiplying by +1 or -1 to get the upper-half
      // or lower-half imaginary value.
      re_nuz1 = -b / a;
      im_nuz1 = -fabs(sqrt(-d) / a);
      re_nuz2 = -b / a;
      im_nuz2 = fabs(sqrt(-d) / a);
    }

    return discrim_flag;
  }