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 20 of file MCShowerRecoPart.h.

Constructor & Destructor Documentation

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

explicit

Default constructor with fhicl parameters.

Definition at line 21 of file MCShowerRecoPart.cxx.

  {
    _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 28 of file MCShowerRecoPart.cxx.

  {
    _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.

    {
      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 62 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 48 of file MCShowerRecoPart.h.

    {
      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 71 of file MCShowerRecoPart.h.

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

protected

Shower time-ordered daughters.

Definition at line 83 of file MCShowerRecoPart.h.

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

protected

Track index to shower index map.

Definition at line 74 of file MCShowerRecoPart.h.

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

protected

Shower Primary Index ID => Shower Index Map.

Definition at line 80 of file MCShowerRecoPart.h.

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

static

Definition at line 24 of file MCShowerRecoPart.h.

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

static

Definition at line 23 of file MCShowerRecoPart.h.

---

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

• MCShowerRecoPart.h
• MCShowerRecoPart.cxx