Lepjets_Event_Lep.cc

Go to the documentation of this file.

//
//
// File: src/Lepjets_Event_Lep.cc
// Purpose: Represent a `lepton' in a Lepjets_Event.
// Created: Jul, 2000, sss, based on run 1 mass analysis code.
//
// CMSSW File      : src/Lepjets_Event_Lep.cc
// Original Author : Scott Stuart Snyder <snyder@bnl.gov> for D0
// Imported to CMSSW by Haryo Sumowidagdo <Suharyo.Sumowidagdo@cern.ch>
//

#include "TopQuarkAnalysis/TopHitFit/interface/Lepjets_Event_Lep.h"

namespace hitfit {

  Lepjets_Event_Lep::Lepjets_Event_Lep(const Fourvec& p, int type, const Vector_Resolution& res)
      //
      // Purpose: Constructor.
      //
      // Inputs:
      //   p -           The 4-momentum.
      //   type -        The type code.
      //   res -         The resolution.
      //
      : _p(p), _type(type), _res(res) {}

  Fourvec& Lepjets_Event_Lep::p()
  //
  // Purpose: Access the 4-momentum.
  //
  // Returns:
  //   The 4-momentum.
  //
  {
    return _p;
  }

  const Fourvec& Lepjets_Event_Lep::p() const
  //
  // Purpose: Access the 4-momentum.
  //
  // Returns:
  //   The 4-momentum.
  //
  {
    return _p;
  }

  int& Lepjets_Event_Lep::type()
  //
  // Purpose: Access the type code.
  //
  // Returns:
  //   The type code.
  //
  {
    return _type;
  }

  int Lepjets_Event_Lep::type() const
  //
  // Purpose: Access the type code.
  //
  // Returns:
  //   The type code.
  //
  {
    return _type;
  }

  const Vector_Resolution& Lepjets_Event_Lep::res() const
  //
  // Purpose: Access the resolutions.
  //
  // Returns:
  //   The resolutions.
  //
  {
    return _res;
  }

  Vector_Resolution& Lepjets_Event_Lep::res()
  //
  // Purpose: Access the resolutions.
  //
  // Returns:
  //   The resolutions.
  //
  {
    return _res;
  }

  double Lepjets_Event_Lep::p_sigma() const
  //
  // Purpose: Return the momentum (or 1/p) resolution for this object.
  //
  // Returns:
  //   The momentum (or 1/p) resolution for this object.
  //
  {
    return _res.p_sigma(_p);
  }

  double Lepjets_Event_Lep::eta_sigma() const
  //
  // Purpose: Return the eta resolution for this object.
  //
  // Returns:
  //   The eta resolution for this object.
  //
  {
    return _res.eta_sigma(_p);
  }

  double Lepjets_Event_Lep::phi_sigma() const
  //
  // Purpose: Return the phi resolution for this object.
  //
  // Returns:
  //   The phi resolution for this object.
  //
  {
    return _res.phi_sigma(_p);
  }

  void Lepjets_Event_Lep::smear(CLHEP::HepRandomEngine& engine, bool smear_dir /*= false*/)
  //
  // Purpose: Smear this object according to its resolutions.
  //
  // Inputs:
  //   engine -      The underlying RNG.
  //   smear_dir -   If false, smear the momentum only.
  //
  {
    _res.smear(_p, engine, smear_dir);
  }

  std::ostream& Lepjets_Event_Lep::dump(std::ostream& s, bool full /*= false*/) const
  //
  // Purpose: Dump out this object.
  //
  // Inputs:
  //   s -           The stream to which to write.
  //   full -        If true, dump the resolutions too.
  //
  // Returns:
  //   The stream S.
  //
  {
    s << "[" << _type << "] " << _p << "; pt: " << _p.perp() << ", eta: " << _p.eta() << ", phi: " << _p.phi();
    if (full) {
      s << "\n    " << _res;
    }
    return s;
  }

  std::ostream& operator<<(std::ostream& s, const Lepjets_Event_Lep& l)
  //
  // Purpose: Dump out this object.
  //
  // Inputs:
  //   s -           The stream to which to write.
  //   l -           The object to dump.
  //
  // Returns:
  //   The stream S.
  //
  {
    return l.dump(s);
  }

  bool Lepjets_Event_Lep::operator<(const Lepjets_Event_Lep& x) const
  //
  // Purpose: Sort objects by pt.
  //
  // Retruns:
  //   True if this object's pt is less than that of x.
  {
    return p().perp2() < x.p().perp2();
  }

}  // namespace hitfit