EtaDepResElement.cc

Go to the documentation of this file.

//
//
//
// File   : src/EtaDepResolution.cc
// Author : Haryo Sumowidagdo <Suharyo.Sumowidagdo@cern.ch>
// Purpose: Hold on to an element eta-dependent resolution.
//
 
#include <sstream>
#include <stdexcept>
#include "TopQuarkAnalysis/TopHitFit/interface/EtaDepResElement.h"
 
namespace hitfit {
 
  EtaDepResElement::EtaDepResElement(double eta1, double eta2, const Vector_Resolution& res) : _Vector_Resolution(res) {
    SetEta(eta1, eta2);
  }
 
  EtaDepResElement::EtaDepResElement(double eta1, double eta2, std::string res) : _Vector_Resolution(res) {
    SetEta(eta1, eta2);
  }
 
  EtaDepResElement::EtaDepResElement(double eta1,
                                     double eta2,
                                     const Resolution& p_res,
                                     const Resolution& eta_res,
                                     const Resolution& phi_res,
                                     bool use_et)
      : _Vector_Resolution(p_res, eta_res, phi_res, use_et) {
    SetEta(eta1, eta2);
  }
 
  EtaDepResElement::~EtaDepResElement() {}
 
  void EtaDepResElement::SetEta(double eta1, double eta2) {
    if (fabs(eta1 - eta2) < (1.0 / double(InverseEtaPrecision))) {
      throw std::runtime_error("EtaDepResElement::equal EtaMin and EtaMax");
    }
 
    if (eta1 < eta2) {
      _EtaMin = eta1;
      _EtaMax = eta2;
    } else {
      _EtaMin = eta2;
      _EtaMax = eta1;
    }
  }
 
  bool operator<(const EtaDepResElement& a, const EtaDepResElement& b) {
    if (a.IsOverlap(b)) {
      return false;
    }
    return !(a._EtaMax > b._EtaMin);
  }
 
  const double EtaDepResElement::EtaMin() const { return _EtaMin; }
 
  const double EtaDepResElement::EtaMax() const { return _EtaMax; }
 
  bool EtaDepResElement::IsOverlap(const EtaDepResElement& e) const {
    return (IsInInterval(e._EtaMin) || IsInInterval(e._EtaMax));
  }
 
  bool EtaDepResElement::IsNotOverlap(const EtaDepResElement& e) const { return !(IsOverlap(e)); }
 
  bool EtaDepResElement::IsInInterval(const double& eta) const { return ((_EtaMin < eta) && (eta < _EtaMax)); }
 
  bool EtaDepResElement::IsOnEdge(const double& eta) const {
    bool nearEtaMin = fabs(eta - _EtaMin) < (1.0 / double(InverseEtaPrecision));
    bool nearEtaMax = fabs(eta - _EtaMax) < (1.0 / double(InverseEtaPrecision));
    return nearEtaMin || nearEtaMax;
  }
 
  bool EtaDepResElement::IsOnEdge(const EtaDepResElement& e) const {
    return (e.IsOnEdge(_EtaMin) || e.IsOnEdge(_EtaMax));
  }
 
  const Vector_Resolution EtaDepResElement::GetResolution() const { return _Vector_Resolution; }
 
  std::ostream& operator<<(std::ostream& s, const EtaDepResElement& e) {
    s << "(" << e._EtaMin << " to " << e._EtaMax << ")"
      << " / " << e.GetResolution();
    return s;
  }
 
}  // namespace hitfit