MCMatchAlg.h

Go to the documentation of this file.

#ifndef RECOTOOL_MCMATCHALG_H
#define RECOTOOL_MCMATCHALG_H

#include <iostream>
//#include "DataFormat/simch.h"
//#include "DataFormat/hit.h"
//#include "LArUtil/Geometry.h"
#include "lardataobj/Simulation/SimChannel.h"
#include "lardataobj/RecoBase/Hit.h"
#include "larcore/Geometry/Geometry.h"
#include "MCBTAlg.h"

namespace btutil {
  class MCMatchAlg {
    
  public:
    
    MCMatchAlg();
    
    virtual ~MCMatchAlg(){};

    
    bool BuildMap(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(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(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