sim::MCShowerRecoPart Class Reference

#include "MCShowerRecoPart.h"

Public Member Functions
	MCShowerRecoPart (fhicl::ParameterSet const &pset)
	Default constructor with fhicl parameters. More...
void	ConstructShower (const MCRecoPart &part_v)
	Main function to read-in data and fill variables in this algorithm to reconstruct MC shower. More...
const std::vector< unsigned int > &	ShowerDaughters (const unsigned int shower_id) const
const std::vector< unsigned int >	ShowerMothers () const
int	ShowerIndex (const unsigned int part_index) const

Static Public Attributes
static const unsigned int	kINVALID_UINT = std::numeric_limits<unsigned int>::max()
static const int	kINVALID_INT = std::numeric_limits<int>::max()

Protected Attributes
bool	_debug_mode
	lots of stdout stream More...
std::vector< int >	_shower_id
	Track index to shower index map. More...
std::map< unsigned int, unsigned int >	_shower_index
	Shower Primary Index ID => Shower Index Map. More...
std::vector< std::vector< unsigned int > >	_shower_daughters
	Shower time-ordered daughters. More...

Detailed Description

Definition at line 21 of file MCShowerRecoPart.h.

Constructor & Destructor Documentation

sim::MCShowerRecoPart::MCShowerRecoPart ( fhicl::ParameterSet const &  pset )

explicit

Default constructor with fhicl parameters.

Definition at line 23 of file MCShowerRecoPart.cxx.

References _debug_mode, and fhicl::ParameterSet::get().

  {
    _debug_mode = pset.get<bool>("DebugMode");
  }

Member Function Documentation

void sim::MCShowerRecoPart::ConstructShower

(

const MCRecoPart &

part_v

)

Main function to read-in data and fill variables in this algorithm to reconstruct MC shower.

Definition at line 30 of file MCShowerRecoPart.cxx.

References _debug_mode, _shower_daughters, _shower_id, _shower_index, and sim::MCRecoPart::_track_index.

Referenced by sim::MCShowerRecoAlg::Reconstruct().

  {
    _shower_id.clear();
    _shower_id.resize(part_v.size(), -1);
    _shower_index.clear();
    _shower_daughters.clear();

    if (!part_v.size()) return;

    // Construct MCShower
    std::vector<std::multimap<double, unsigned int>> daughter_map;
    for (size_t i = 0; i < part_v.size(); ++i) {

      auto const& mcp = part_v[i];

      int candidate_mom_index = -1;
      if (mcp._pdgcode == 22 || mcp._pdgcode == 11 || mcp._pdgcode == -11) candidate_mom_index = i;

      unsigned int mom_track = mcp._mother;
      auto mom_iter = part_v._track_index.find(mom_track);
      while (mom_iter != part_v._track_index.end()) {

        unsigned int mom_index = (*mom_iter).second;

        if (part_v.at(mom_index)._pdgcode == 22 || part_v.at(mom_index)._pdgcode == 11 ||
            part_v.at(mom_index)._pdgcode == -11)

          candidate_mom_index = mom_index;

        mom_iter = part_v._track_index.find(part_v.at(mom_index)._mother);
      }

      if (candidate_mom_index >= 0) {

        auto candidate_mom_iter = _shower_index.find(candidate_mom_index);
        if (candidate_mom_iter == _shower_index.end()) {
          _shower_index.insert(
            std::make_pair((unsigned int)candidate_mom_index, (unsigned int)_shower_index.size()));
          daughter_map.push_back(std::multimap<double, unsigned int>());
        }
        unsigned int shower_index = (*_shower_index.find(candidate_mom_index)).second;
        daughter_map.at(shower_index)
          .insert(std::make_pair((double)(mcp._start_vtx[3]), (unsigned int)i));
        _shower_id.at(i) = shower_index;
      }
      else if (_debug_mode) {

        std::cout << "Found a particle that does not belong to a shower!" << std::endl
                  << Form(" PDGID: %d ... Track %d @ (%g,%g,%g,%g) with (%g,%g,%g,%g)",
                          mcp._pdgcode,
                          mcp._track_id,
                          mcp._start_vtx[0],
                          mcp._start_vtx[1],
                          mcp._start_vtx[2],
                          mcp._start_vtx[3],
                          mcp._start_mom[0],
                          mcp._start_mom[1],
                          mcp._start_mom[2],
                          mcp._start_mom[3])
                  << std::endl
                  << std::endl;
      }
    }

    if (_debug_mode)
      std::cout << Form("Found %zu MCShowers....", _shower_index.size()) << std::endl;

    _shower_daughters.resize(_shower_index.size(), std::vector<unsigned int>());
    for (auto const& mom : _shower_index) {

      _shower_daughters.at(mom.second).reserve(daughter_map.at(mom.second).size());
      for (auto const& part_index : daughter_map.at(mom.second))

        _shower_daughters.at(mom.second).push_back(part_index.second);

      auto const& mcp = part_v.at(mom.first);
      if (_debug_mode)
        std::cout
          << Form("PDGID: %d ... Track %d @ (%g,%g,%g,%g) with (%g,%g,%g,%g) ... %zu daughters!",
                  mcp._pdgcode,
                  mcp._track_id,
                  mcp._start_vtx[0],
                  mcp._start_vtx[1],
                  mcp._start_vtx[2],
                  mcp._start_vtx[3],
                  mcp._start_mom[0],
                  mcp._start_mom[1],
                  mcp._start_mom[2],
                  mcp._start_mom[3],
                  _shower_daughters.at(mom.second).size())
          << std::endl;
    }

    if (_debug_mode) std::cout << std::endl;
  }

const std::vector<unsigned int>& sim::MCShowerRecoPart::ShowerDaughters ( const unsigned int  shower_id ) const

inline

Returns a list ot daughter particle index numbers for the specified shower with the shower index number as an input

Definition at line 39 of file MCShowerRecoPart.h.

Referenced by sim::MCShowerRecoAlg::Reconstruct().

    {
      if (shower_id >= _shower_daughters.size())
        throw cet::exception(__FUNCTION__) << "Invalid shower index!";
      return _shower_daughters.at(shower_id);
    }

int sim::MCShowerRecoPart::ShowerIndex ( const unsigned int  part_index ) const

inline

Take particle index number and returns shower index number to which this particle belongs. Returns -1 if a particle does not belong to any MC shower. Returns kINVALID_INT if input is invalid.

Definition at line 64 of file MCShowerRecoPart.h.

    {
      if (_shower_id.size() <= part_index) return kINVALID_INT;
      return _shower_id.at(part_index);
    }

const std::vector<unsigned int> sim::MCShowerRecoPart::ShowerMothers ( ) const

inline

Returns a list of shower-mother's particle index. Order respects shower index number.

Definition at line 49 of file MCShowerRecoPart.h.

Referenced by sim::MCShowerRecoAlg::Reconstruct().

    {
      std::vector<unsigned int> mothers(_shower_index.size(), 0);
      for (auto mother_iter = _shower_index.begin(); mother_iter != _shower_index.end();
           ++mother_iter)
        mothers.at((*mother_iter).second) = (*mother_iter).first;
      return mothers;
    }

Member Data Documentation

bool sim::MCShowerRecoPart::_debug_mode

protected

lots of stdout stream

Definition at line 72 of file MCShowerRecoPart.h.

Referenced by ConstructShower(), and MCShowerRecoPart().

std::vector<std::vector<unsigned int> > sim::MCShowerRecoPart::_shower_daughters

protected

Shower time-ordered daughters.

Definition at line 84 of file MCShowerRecoPart.h.

Referenced by ConstructShower().

std::vector<int> sim::MCShowerRecoPart::_shower_id

protected

Track index to shower index map.

Definition at line 75 of file MCShowerRecoPart.h.

Referenced by ConstructShower().

std::map<unsigned int, unsigned int> sim::MCShowerRecoPart::_shower_index

protected

Shower Primary Index ID => Shower Index Map.

Definition at line 81 of file MCShowerRecoPart.h.

Referenced by ConstructShower().

const int sim::MCShowerRecoPart::kINVALID_INT = std::numeric_limits<int>::max()

static

Definition at line 25 of file MCShowerRecoPart.h.

const unsigned int sim::MCShowerRecoPart::kINVALID_UINT = std::numeric_limits<unsigned int>::max()

static

Definition at line 24 of file MCShowerRecoPart.h.

---

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

• larsim/v09_43_00/source/larsim/MCSTReco/MCShowerRecoPart.h
• larsim/v09_43_00/source/larsim/MCSTReco/MCShowerRecoPart.cxx