KHitGroup.cxx

Go to the documentation of this file.

#include "lardata/RecoObjects/KHitGroup.h"
#include "cetlib_except/exception.h"

namespace trkf {

  KHitGroup::KHitGroup(bool has_path, double path) : fPlane(-1), fHasPath(has_path), fPath(path) {}

  KHitGroup::~KHitGroup() {}

  void KHitGroup::addHit(const std::shared_ptr<const KHitBase>& hit)
  {
    // Make sure that the measurement pointer is not null (throw exception if null).

    if (hit.get() == 0) throw cet::exception("KHitGroup") << "Attempt to add null measurement.\n";

    // If the stored common surface pointer has not yet been initialized,
    // initialize it now.
    // Otherwise, make sure that the new measurement surface matches
    // the common surface (pointer to same surface object).
    // Throw exception if the new surface doesn't match.

    if (hit->getMeasPlane() < 0)
      throw cet::exception("KHitGroup")
        << __func__ << ": invalid hit plane " << hit->getMeasPlane() << "\n";

    if (fSurf.get() == 0) {
      fSurf = hit->getMeasSurface();
      fPlane = hit->getMeasPlane();
    }
    else {
      if (fSurf.get() != hit->getMeasSurface().get())
        throw cet::exception("KHitGroup") << "Attempt to add non-matching measurement.\n";
      if (hit->getMeasPlane() != fPlane) {
        throw cet::exception("KHitGroup") << __func__ << ": hit plane mismatch, "
                                          << hit->getMeasPlane() << " vs. " << fPlane << "\n";
      }
      if (fPlane < 0)
        throw cet::exception("KHitGroup") << __func__ << ": invalid plane " << fPlane << "\n";
    }

    // Everything OK.  Add the measurement.

    fHits.push_back(hit);
  }

  bool KHitGroup::operator==(const KHitGroup& obj) const
  {
    bool result = (!fHasPath && !obj.fHasPath) || (fHasPath && obj.fHasPath && fPath == obj.fPath);
    return result;
  }

  bool KHitGroup::operator<(const KHitGroup& obj) const
  {
    bool result = false;
    if (fHasPath != obj.fHasPath)
      throw cet::exception("KHitGroup") << "Attempt to compare incomparable objects.\n";
    if (fHasPath) result = fPath < obj.fPath;
    return result;
  }

} // end namespace trkf