ShowerRecoManager.cxx

Go to the documentation of this file.

#include "ShowerRecoManager.h"

#include "lardata/Utilities/GeometryUtilities.h"
#include "lardataobj/RecoBase/Shower.h"
#include "larreco/RecoAlg/CMTool/CMToolBase/CMTException.h"
#include "larreco/RecoAlg/CMTool/CMToolBase/CMatchBookKeeper.h"
#include "larreco/RecoAlg/CMTool/CMToolBase/CMatchManager.h"
#include "larreco/RecoAlg/ClusterRecoUtil/ClusterParamsAlg.h"
#include "larreco/ShowerFinder/ShowerReco3D/ShowerRecoAlgBase.h"

namespace showerreco {

  ShowerRecoManager::ShowerRecoManager(unsigned int const num_planes)
    : fShowerAlgo(nullptr), fMatchMgr(nullptr)
  {
    fMatch = true;
    fMatchMgr = new ::cmtool::CMatchManager(num_planes);
  }

  void ShowerRecoManager::Reset()
  {
    if (fShowerAlgo) fShowerAlgo->Reset();
    fMatchMgr->Reset();
  }

  ClusterAss_t ShowerRecoManager::Reconstruct(geo::GeometryCore const& geom,
                                              geo::WireReadoutGeom const& wireReadoutGeom,
                                              detinfo::DetectorClocksData const& clockData,
                                              detinfo::DetectorPropertiesData const& detProp,
                                              const std::vector<std::vector<util::PxHit>>& clusters,
                                              std::vector<::recob::Shower>& showers)
  {
    util::GeometryUtilities const gser{geom, wireReadoutGeom, clockData, detProp};
    showers.clear();
    fMatchMgr->SetClusters(gser, clusters);

    ClusterAss_t res_ass;
    // Run matching & retrieve matched cluster indices
    try {
      fMatchMgr->Process(gser);
    }
    catch (::cmtool::CMTException& e) {
      e.what();
      return res_ass;
    }
    res_ass = fMatchMgr->GetBookKeeper().GetResult();

    Process(geom, wireReadoutGeom, clockData, detProp, res_ass, showers);

    return res_ass;
  }

  void ShowerRecoManager::Reconstruct(geo::GeometryCore const& geom,
                                      geo::WireReadoutGeom const& wireReadoutGeom,
                                      detinfo::DetectorClocksData const& clockData,
                                      detinfo::DetectorPropertiesData const& detProp,
                                      const std::vector<std::vector<util::PxHit>>& clusters,
                                      const ClusterAss_t& ass,
                                      std::vector<::recob::Shower>& showers)
  {
    showers.clear();
    util::GeometryUtilities const gser{geom, wireReadoutGeom, clockData, detProp};
    fMatchMgr->SetClusters(gser, clusters);

    Process(geom, wireReadoutGeom, clockData, detProp, ass, showers);
  }

  void ShowerRecoManager::Process(geo::GeometryCore const& geom,
                                  geo::WireReadoutGeom const& wireReadoutGeom,
                                  detinfo::DetectorClocksData const& clockData,
                                  detinfo::DetectorPropertiesData const& detProp,
                                  const ClusterAss_t& ass,
                                  std::vector<::recob::Shower>& showers)
  {

    for (auto const& pair : ass) {
      std::vector<::cluster::ClusterParamsAlg> cpans;

      cpans.reserve(pair.size());

      for (auto const& index : pair)
        cpans.push_back(fMatchMgr->GetInputClusters()[index]);

      fShowerAlgo->AppendInputClusters(cpans);
    }

    // Run shower reco
    showers = fShowerAlgo->Reconstruct(geom, wireReadoutGeom, clockData, detProp);
  }

}