MCMatchAlg.h

Go to the documentation of this file.

#ifndef RECOTOOL_MCMATCHALG_H
#define RECOTOOL_MCMATCHALG_H

#include "canvas/Persistency/Common/Ptr.h"

#include "MCBTAlg.h"
#include "lardataobj/RecoBase/Hit.h"
#include "lardataobj/Simulation/SimChannel.h"
namespace detinfo {
  class DetectorClocksData;
}

#include <cstddef>
#include <utility>
#include <vector>

namespace btutil {
  class MCMatchAlg {

  public:
    MCMatchAlg();

    virtual ~MCMatchAlg(){};

    bool BuildMap(detinfo::DetectorClocksData const& clockData,
                  const std::vector<unsigned int>& g4_trackid_v,
                  const std::vector<sim::SimChannel>& simch_v,
                  const std::vector<std::vector<art::Ptr<recob::Hit>>>& cluster_v);

    bool BuildMap(detinfo::DetectorClocksData const& clockData,
                  const std::vector<std::vector<unsigned int>>& g4_trackid_v,
                  const std::vector<sim::SimChannel>& simch_v,
                  const std::vector<std::vector<art::Ptr<recob::Hit>>>& cluster_v);

    const std::vector<int>& BestClusters(const size_t mcshower_index) const;

    double ClusterCorrectness(const size_t cluster_index, const size_t mcshower_index) const;

    std::pair<size_t, double> ShowerCorrectness(
      const std::vector<unsigned int> cluster_indices) const;

    std::pair<double, double> ClusterEP(const size_t cluster_index,
                                        const size_t mcshower_index) const;

    std::pair<double, double> BestClusterEP(const size_t mcshower_index,
                                            const size_t plane_id) const;

    const MCBTAlg& BTAlg() const { return fBTAlgo; }

  protected:
    bool BuildMap(detinfo::DetectorClocksData const& clockData,
                  const std::vector<std::vector<art::Ptr<recob::Hit>>>& cluster_v);

    MCBTAlg fBTAlgo;

    std::vector<size_t> _view_to_plane;

    std::vector<std::vector<double>> _summed_mcq;
    std::vector<std::vector<double>> _cluster_mcq_v;

    std::vector<unsigned char> _cluster_plane_id;
    std::vector<std::vector<int>> _bmatch_id;
  };
}
#endif
 // end of doxygen group