MCBTAlg.h

Go to the documentation of this file.

#ifndef RECOTOOL_MCBTALG_H
#define RECOTOOL_MCBTALG_H

#include "lardataobj/Simulation/SimChannel.h"

#include <map>
#include <vector>

namespace detinfo {
  class DetectorClocksData;
}

namespace btutil {

  struct WireRange_t {
    unsigned int ch;
    double start, end;
    WireRange_t()
    {
      ch = std::numeric_limits<unsigned int>::max();
      start = end = std::numeric_limits<double>::max();
    }
    WireRange_t(unsigned int c, double s, double e)
    {
      ch = c;
      start = s;
      end = e;
    }
  };

  typedef std::vector<double> edep_info_t; // vector of energy deposition

  typedef std::map<unsigned int, ::btutil::edep_info_t>
    ch_info_t; // vector of time (index) for each edep (value)

  class MCBTAlg {

  public:
    MCBTAlg() {}

    MCBTAlg(const std::vector<unsigned int>& g4_trackid_v,
            const std::vector<sim::SimChannel>& simch_v);

    void Reset(const std::vector<unsigned int>& g4_trackid_v,
               const std::vector<sim::SimChannel>& simch_v);

    void Reset(const std::vector<std::vector<unsigned int>>& g4_trackid_v,
               const std::vector<sim::SimChannel>& simch_v);

    const std::vector<double>& MCQSum(const size_t plane_id) const;

    std::vector<double> MCQ(detinfo::DetectorClocksData const& clockData,
                            const WireRange_t& hit) const;

    std::vector<double> MCQFrac(detinfo::DetectorClocksData const& clockData,
                                const WireRange_t& hit) const;

    std::vector<double> MCQ(detinfo::DetectorClocksData const& clockData,
                            const std::vector<btutil::WireRange_t>& hit_v) const;

    std::vector<double> MCQFrac(detinfo::DetectorClocksData const& clockData,
                                const std::vector<btutil::WireRange_t>& hit_v) const;

    size_t Index(const unsigned int g4_track_id) const;

    size_t NumParts() const { return _num_parts - 1; }

  protected:
    void Register(const unsigned int& g4_track_id);

    void Register(const std::vector<unsigned int>& g4_track_id);

    void ProcessSimChannel(const std::vector<sim::SimChannel>& simch_v);

    std::vector<::btutil::ch_info_t> _event_info;
    std::vector<size_t> _trkid_to_index;
    std::vector<std::vector<double>> _sum_mcq;
    size_t _num_parts;
  };
}
#endif
 // end of doxygen group