#include "MCBTAlg.h"

Public Member Functions
	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

Protected Member Functions
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)

Protected Attributes
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

Detailed Description

Definition at line 54 of file MCBTAlg.h.

Constructor & Destructor Documentation

btutil::MCBTAlg::MCBTAlg ( )

inline

Definition at line 57 of file MCBTAlg.h.

References Reset(), and lar::dump::vector().

57 {}

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

Definition at line 16 of file MCBTAlg.cxx.

References Reset().

   {
     Reset(g4_trackid_v, simch_v);
   }

Member Function Documentation

size_t MCBTAlg::Index ( const unsigned int g4_track_id ) const

Definition at line 164 of file MCBTAlg.cxx.

References _trkid_to_index, and btutil::kINVALID_INDEX.

   {
     if (g4_track_id >= _trkid_to_index.size()) return kINVALID_INDEX;
     return _trkid_to_index[g4_track_id];
   }

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

Relate Hit => MCShower/MCTrack (called MCX). Returns a vector of double w/ length = # of relevant MCX + 1. Each entry is # drifted electrons from each relevant MCX. The last element contains a sum of drifted electrons that do not belong to any of relevant MCX.

Definition at line 97 of file MCBTAlg.cxx.

References _event_info, _num_parts, btutil::WireRange_t::ch, btutil::WireRange_t::end, btutil::WireRange_t::start, and detinfo::DetectorClocksData::TPCTick2TDC().

Referenced by MCBTDemo::analyze(), ShowerQuality::analyze(), btutil::MCMatchAlg::BuildMap(), MCQ(), and MCQFrac().

   {
     std::vector<double> res(_num_parts, 0);
 
     if (_event_info.size() <= hit.ch) return res;
 
     auto const& ch_info = _event_info[hit.ch];
 
     auto itlow = ch_info.lower_bound((unsigned int)(clockData.TPCTick2TDC(hit.start)));
     auto itup = ch_info.upper_bound((unsigned int)(clockData.TPCTick2TDC(hit.end)) + 1);
 
     while (itlow != ch_info.end() && itlow != itup) {
 
       auto const& edep_info = (*itlow).second;
 
       for (size_t part_index = 0; part_index < _num_parts; ++part_index)
 
         res[part_index] += edep_info[part_index];
 
       ++itlow;
     }
     return res;
   }

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

Relate Cluster => MCX. Returns a vector of double w/ length = # of relevant MCXs + 1. Each entry is # drifted electrons from each relevant MCX. The last element contains a sum of drifted electrons that do not belong to any of relevant MCX.

Definition at line 137 of file MCBTAlg.cxx.

References _num_parts, and MCQ().

   {
     std::vector<double> res(_num_parts, 0);
     for (auto const& h : hit_v) {
       auto tmp_res = MCQ(clockData, h);
       for (size_t i = 0; i < res.size(); ++i)
         res[i] += tmp_res[i];
     }
     return res;
   }

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

Relate Hit => MCX. Returns a vector of double w/ length = # of relevant MCXs + 1. Each entry is a fraction of # drifted electrons within the specified time range from each relevant MCX. The last element contains a sum of drifted electrons that do not belong to any of relevant MCX.

Definition at line 122 of file MCBTAlg.cxx.

References MCQ(), and sum.

Referenced by MCBTDemo::analyze().

   {
     auto res = MCQ(clockData, hit);
     if (!res.size()) return res;
 
     double sum = 0;
     for (auto const& v : res)
       sum += v;
     for (size_t i = 0; i < (res.size() - 1); ++i)
       res[i] /= (sum - (*res.rbegin()));
     (*res.rbegin()) /= sum;
     return res;
   }

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

Relate Cluster => MCX. Returns a vector of double w/ length = # of relevant MCXs + 1. Each entry is a fraction of # drifted electrons within the specified time range from each relevant MCX. The last element contains a sum of drifted electrons that do not belong to any of relevant MCX.

Definition at line 149 of file MCBTAlg.cxx.

References MCQ(), and sum.

   {
     auto res = MCQ(clockData, hit_v);
     if (!res.size()) return res;
 
     double sum = 0;
     for (auto const& v : res)
       sum += v;
     for (size_t i = 0; i < (res.size() - 1); ++i)
       res[i] /= (sum - (*res.rbegin()));
     (*res.rbegin()) /= sum;
     return res;
   }

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

Returns MC charge sum per MCX for a specified plane

Definition at line 90 of file MCBTAlg.cxx.

References _sum_mcq.

Referenced by MCBTDemo::analyze().

   {
     if (plane_id > _sum_mcq.size())
       throw MCBTException(Form("Invalid plane requested: %zu", plane_id));
     return _sum_mcq[plane_id];
   }

size_t btutil::MCBTAlg::NumParts ( ) const

inline

Definition at line 115 of file MCBTAlg.h.

Referenced by btutil::MCMatchAlg::BuildMap().

115 { return _num_parts - 1; }

btutil::MCBTAlg::_num_parts

size_t _num_parts

Definition: MCBTAlg.h:127

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

protected

Definition at line 50 of file MCBTAlg.cxx.

References _event_info, _num_parts, _sum_mcq, _trkid_to_index, Get, and btutil::kINVALID_INDEX.

Referenced by Reset().

   {
     auto const& wireReadoutGeom = art::ServiceHandle<geo::WireReadout const>()->Get();
     _sum_mcq.resize(wireReadoutGeom.Nplanes(), std::vector<double>(_num_parts, 0));
 
     for (auto const& sch : simch_v) {
 
       auto const ch = sch.Channel();
       if (_event_info.size() <= ch) _event_info.resize(ch + 1);
 
       auto& ch_info = _event_info[ch];
 
       size_t plane = wireReadoutGeom.ChannelToWire(ch)[0].Plane;
 
       for (auto const& time_ide : sch.TDCIDEMap()) {
 
         auto const& time = time_ide.first;
         auto const& ide_v = time_ide.second;
 
         auto& edep_info = ch_info[time];
 
         if (!edep_info.size()) edep_info.resize(_num_parts, 0);
 
         for (auto const& ide : ide_v) {
 
           size_t index = kINVALID_INDEX;
           if (ide.trackID < (int)(_trkid_to_index.size())) { index = _trkid_to_index[ide.trackID]; }
           if (_num_parts <= index) {
             (*edep_info.rbegin()) += ide.numElectrons;
             (*(_sum_mcq[plane]).rbegin()) += ide.numElectrons;
           }
           else {
             edep_info[index] += ide.numElectrons;
             _sum_mcq[plane][index] += ide.numElectrons;
           }
         }
       }
     }
   }

void MCBTAlg::Register ( const unsigned int & g4_track_id )

protected

Definition at line 170 of file MCBTAlg.cxx.

References _num_parts, _trkid_to_index, and btutil::kINVALID_INDEX.

Referenced by Reset().

   {
     if (_trkid_to_index.size() <= g4_track_id)
       _trkid_to_index.resize(g4_track_id + 1, kINVALID_INDEX);
 
     if (_trkid_to_index[g4_track_id] == kINVALID_INDEX) {
       _trkid_to_index[g4_track_id] = _num_parts;
       ++_num_parts;
     }
   }

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

protected

Definition at line 181 of file MCBTAlg.cxx.

References _num_parts, _trkid_to_index, and btutil::kINVALID_INDEX.

   {
     unsigned int max_id = 0;
     for (auto const& id : track_id_v)
       if (max_id < id) max_id = id;
     if (_trkid_to_index.size() <= max_id) _trkid_to_index.resize(max_id + 1, kINVALID_INDEX);
 
     for (auto const& id : track_id_v) {
 
       if (_trkid_to_index[id] == kINVALID_INDEX)
 
         _trkid_to_index[id] = _num_parts;
 
       else
 
         throw MCBTException(Form("Doubly used TrackID: %d", id));
     }
     ++_num_parts;
   }

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

Definition at line 22 of file MCBTAlg.cxx.

References _event_info, _num_parts, _sum_mcq, _trkid_to_index, ProcessSimChannel(), and Register().

Referenced by btutil::MCMatchAlg::BuildMap(), and MCBTAlg().

   {
     _num_parts = 0;
     _sum_mcq.clear();
     _trkid_to_index.clear();
     _event_info.clear();
     //
     for (auto const& id : g4_trackid_v)
       Register(id);
     _num_parts++;
     ProcessSimChannel(simch_v);
   }

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

Definition at line 36 of file MCBTAlg.cxx.

References _event_info, _num_parts, _sum_mcq, _trkid_to_index, ProcessSimChannel(), and Register().

   {
     _num_parts = 0;
     _sum_mcq.clear();
     _trkid_to_index.clear();
     _event_info.clear();
     //
     for (auto const& id : g4_trackid_v)
       Register(id);
     _num_parts++;
     ProcessSimChannel(simch_v);
   }