HitsStandard_tool.cc

Go to the documentation of this file.

// Framework Includes
#include "art/Utilities/ToolMacros.h"
#include "messagefacility/MessageLogger/MessageLogger.h"

// LArSoft includes
#include "lardataobj/RecoBase/Hit.h"
#include "larreco/SpacePointSolver/HitReaders/IHitReader.h"

// std includes
#include <cmath>
#include <ostream>
#include <vector>

//------------------------------------------------------------------------------------------------------------------------------------------
// implementation follows

namespace reco3d {

  class HitsStandard : virtual public IHitReader {
  public:
    explicit HitsStandard(const fhicl::ParameterSet&);

    ~HitsStandard();

    void configure(fhicl::ParameterSet const& pset) override;

    bool readHits(const std::vector<art::Ptr<recob::Hit>>&,           // input hits
                  std::vector<art::Ptr<recob::Hit>>&,                 // output hits plane 0
                  std::vector<art::Ptr<recob::Hit>>&,                 // output hits plane 1
                  std::vector<art::Ptr<recob::Hit>>&) const override; // output hits plane 2
  };

  HitsStandard::HitsStandard(fhicl::ParameterSet const& pset)
  {
    this->configure(pset);
  }

  //------------------------------------------------------------------------------------------------------------------------------------------

  HitsStandard::~HitsStandard() {}

  //------------------------------------------------------------------------------------------------------------------------------------------

  void HitsStandard::configure(fhicl::ParameterSet const&)
  {
    //    m_enableMonitoring   = pset.get<bool>  ("EnableMonitoring",  true  );

    return;
  }

  bool HitsStandard::readHits(const std::vector<art::Ptr<recob::Hit>>& inputHits, // input hits
                              std::vector<art::Ptr<recob::Hit>>& xhits,       // output hits plane 0
                              std::vector<art::Ptr<recob::Hit>>& uhits,       // output hits plane 1
                              std::vector<art::Ptr<recob::Hit>>& vhits) const // output hits plane 2
  {

    bool is2view = false;

    for (auto& hit : inputHits) {
      if (hit->Integral() < 0 || std::isnan(hit->Integral()) || std::isinf(hit->Integral())) {
        mf::LogWarning("HitsStandard")
          << "WARNING: bad recob::Hit::Integral() = " << hit->Integral() << ". Skipping."
          << std::endl;
        continue;
      }

      if (hit->SignalType() == geo::kCollection) {
        // For DualPhase, both view are collection. Arbitrarily map V to the main
        // "X" view. For Argoneut and Lariat, collection=V is also the right
        // convention.
        if (hit->View() == geo::kZ) { xhits.push_back(hit); }
        if (hit->View() == geo::kV) {
          xhits.push_back(hit);
          is2view = true;
        }
        if (hit->View() == geo::kU || hit->View() == geo::kY) {
          uhits.push_back(hit);
          is2view = true;
        }
      }
      else {
        if (hit->View() == geo::kU) uhits.push_back(hit);
        if (hit->View() == geo::kV) vhits.push_back(hit);
      }
    } // end for hit

    mf::LogDebug("HitsStandard") << ">>>>> Reading hits done" << std::endl;

    return is2view;
  }

  DEFINE_ART_CLASS_TOOL(HitsStandard)
} // namespace lar_cluster3d