FPGATrackSimInputUtils Namespace Reference

Typedefs
typedef std::bitset< NBITS >	ParentBitmask

Enumerations
enum	Bits {   TAU_PARENT_BIT, B_PARENT_BIT, PION_PARENT_BIT, PION_IMMEDIATE_PARENT_BIT,   NBITS }

Functions
const ParentBitmask	construct_truth_bitmap (HepMC::ConstGenParticlePtr p)

Typedef Documentation

ParentBitmask

typedef std::bitset<NBITS> FPGATrackSimInputUtils::ParentBitmask

Definition at line 16 of file FPGATrackSimInputUtils.h.

Enumeration Type Documentation

Bits

enum FPGATrackSimInputUtils::Bits

Enumerator
TAU_PARENT_BIT
B_PARENT_BIT
PION_PARENT_BIT
PION_IMMEDIATE_PARENT_BIT
NBITS

Definition at line 15 of file FPGATrackSimInputUtils.h.

{ TAU_PARENT_BIT , B_PARENT_BIT , PION_PARENT_BIT , PION_IMMEDIATE_PARENT_BIT , NBITS } Bits;

Function Documentation

construct_truth_bitmap()

const ParentBitmask FPGATrackSimInputUtils::construct_truth_bitmap

(

HepMC::ConstGenParticlePtr

p

)

Definition at line 10 of file FPGATrackSimInputUtils.cxx.

  {
    ParentBitmask result;
    result.reset();
    typedef std::pair<HepMC::ConstGenParticlePtr, unsigned int> Parent;
    std::vector<Parent> parents;
    parents.push_back(Parent(particle, 0));
    while (!parents.empty()) {
      HepMC::ConstGenParticlePtr p = parents.back().first;
      const unsigned int level = parents.back().second;
      if (std::abs(p->pdg_id()) == 15) { result.set(TAU_PARENT_BIT, 1); }
      if (std::abs(p->pdg_id()) == 5) { result.set(B_PARENT_BIT, 1); }
      if (std::abs(p->pdg_id()) == 211) { result.set(PION_PARENT_BIT, 1); }
      if (std::abs(p->pdg_id()) == 211 && level <= 1) { result.set(PION_IMMEDIATE_PARENT_BIT, 1); }
      if (result.count() == NBITS) { break; }
      parents.pop_back();
      if (!(p->production_vertex())) { continue; }
#ifdef HEPMC3
      for (HepMC::ConstGenParticlePtr i: p->production_vertex()->particles_in() ) {
        parents.push_back(Parent(i, level + 1));
      }
#else
      for (HepMC::GenVertex::particle_iterator i = p->production_vertex()->particles_begin(HepMC::parents), f = p->production_vertex()->particles_end(HepMC::parents); i != f; ++i) {
        parents.push_back(Parent(*i, level + 1));
      }
 
#endif
    }
    return result;
  }