KHitContainerWireLine.cxx

Go to the documentation of this file.

#include "lardata/RecoObjects/KHitContainerWireLine.h"

#include "art/Framework/Services/Registry/ServiceHandle.h"
#include "cetlib_except/exception.h"

#include "larcore/Geometry/Geometry.h"
#include "lardata/RecoObjects/KHitWireLine.h"

namespace trkf {

  void KHitContainerWireLine::fill(const detinfo::DetectorPropertiesData& detProp,
                                   const art::PtrVector<recob::Hit>& hits,
                                   int only_plane)
  {
    // Get services.

    art::ServiceHandle<geo::Geometry const> geom;

    // Loop over hits.

    for (art::PtrVector<recob::Hit>::const_iterator ihit = hits.begin(); ihit != hits.end();
         ++ihit) {
      const recob::Hit& hit = **ihit;

      // Extract the wireid from the Hit.
      geo::WireID hitWireID = hit.WireID();

      uint32_t channel = hit.Channel();

      // Choose plane.
      if (only_plane >= 0 && hitWireID.Plane != (unsigned int)(only_plane)) continue;

      // Make a new KHitGroup for each hit.

      getUnsorted().push_back(KHitGroup());
      KHitGroup* pgr = &(getUnsorted().back());
      if (!pgr) {
        throw cet::exception("KHitContainerWireLine")
          << __func__ << ": no group map for channel " << channel << "\n";
      }

      pgr->addHit(std::make_shared<KHitWireLine>(detProp, *ihit, pgr->getSurface()));
    }
  }

} // end namespace trkf