btutil::MCBTAlg Class Reference

#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 (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

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 51 of file MCBTAlg.h.

Constructor & Destructor Documentation

btutil::MCBTAlg::MCBTAlg ( )

inline

Definition at line 55 of file MCBTAlg.h.

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

{}

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

Definition at line 8 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 159 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

(

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 95 of file MCBTAlg.cxx.

References _event_info, _num_parts, btutil::WireRange_t::ch, btutil::WireRange_t::end, btutil::WireRange_t::start, and detinfo::DetectorClocks::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];

    const detinfo::DetectorClocks* ts = lar::providerFrom<detinfo::DetectorClocksService>();
    //auto ts = ::larutil::TimeService::GetME();

    auto itlow = ch_info.lower_bound((unsigned int)(ts->TPCTick2TDC(hit.start)));
    auto itup  = ch_info.upper_bound((unsigned int)(ts->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

(

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 136 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(h);
      for(size_t i=0; i<res.size(); ++i) res[i] += tmp_res[i];
    }
    return res;
  }

std::vector< double > MCBTAlg::MCQFrac

(

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 123 of file MCBTAlg.cxx.

References MCQ().

Referenced by MCBTDemo::analyze().

  { 
    auto res = MCQ(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

(

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 146 of file MCBTAlg.cxx.

References MCQ().

  {
    auto res = MCQ(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 88 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 109 of file MCBTAlg.h.

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

{ return _num_parts-1; }

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

protected

Definition at line 43 of file MCBTAlg.cxx.

References _event_info, _num_parts, _sum_mcq, _trkid_to_index, geo::GeometryCore::ChannelToWire(), btutil::kINVALID_INDEX, and geo::GeometryCore::Nplanes().

Referenced by Reset().

  {

    art::ServiceHandle<geo::Geometry> geo;
    //auto geo = ::larutil::Geometry::GetME();
    _sum_mcq.resize(geo->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 = geo->ChannelToWire(ch)[0].Plane;
      //size_t plane = geo->ChannelToPlane(ch);

      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 165 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 176 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 14 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 29 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);
  }

Member Data Documentation

std::vector< ::btutil::ch_info_t> btutil::MCBTAlg::_event_info

protected

Definition at line 119 of file MCBTAlg.h.

Referenced by MCQ(), ProcessSimChannel(), and Reset().

size_t btutil::MCBTAlg::_num_parts

protected

Definition at line 122 of file MCBTAlg.h.

Referenced by MCQ(), ProcessSimChannel(), Register(), and Reset().

std::vector<std::vector<double> > btutil::MCBTAlg::_sum_mcq

protected

Definition at line 121 of file MCBTAlg.h.

Referenced by MCQSum(), ProcessSimChannel(), and Reset().

std::vector<size_t> btutil::MCBTAlg::_trkid_to_index

protected

Definition at line 120 of file MCBTAlg.h.

Referenced by Index(), ProcessSimChannel(), Register(), and Reset().

---

The documentation for this class was generated from the following files:

• /cvmfs/larsoft.opensciencegrid.org/products/larreco/v07_10_02/source/larreco/MCComp/MCBTAlg.h
• /cvmfs/larsoft.opensciencegrid.org/products/larreco/v07_10_02/source/larreco/MCComp/MCBTAlg.cxx