MCBTAlg.h

Go to the documentation of this file.

#ifndef RECOTOOL_MCBTALG_H
#define RECOTOOL_MCBTALG_H

#include <iostream>
#include <vector>
#include <map>
#include <set>
//#include "DataFormat/simch.h"
//#include "LArUtil/TimeService.h"
//#include "LArUtil/Geometry.h"
#include "lardataobj/Simulation/SimChannel.h"
#include "lardata/DetectorInfoServices/DetectorClocksService.h"
#include "larcore/Geometry/Geometry.h"
#include "MCBTAlgConstants.h"
#include "MCBTException.h"
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(const WireRange_t& hit) const;

    std::vector<double> MCQFrac(const WireRange_t& hit) const;

    std::vector<double> MCQ(const std::vector<btutil::WireRange_t>& hit_v) const;

    std::vector<double> MCQFrac(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