PartToHitTable.cxx

Go to the documentation of this file.

#include <map>
#include "Simulation/PartToHitTable.h"
#include "Simulation/ParticleList.h"
#include "Simulation/FLSHitList.h"
using namespace sim;

PartToHitTable::PartToHitTable(const ParticleList& plist,
                               const FLSHitList& hitlist) :
  fHitLists(plist.fParticles.size())
{
  unsigned register int i;
  //
  // The sizes can be large, so try to do this as efficiently as
  // possible. The following makes only two passes over the lists
  //
  // First build a look up table to map trackIds to an index inside
  // the particle list
  std::map<int,int> trk2indx;
  for (i=0; i<plist.fParticles.size(); ++i) {
    trk2indx[plist.fParticles[i].fTrackId] = i;
  }
  
  // Now, use that map to build the lists
  for (i=0; i<hitlist.fHits.size(); ++i) {
    int indx = trk2indx[hitlist.fHits[i].fTrackId];
    fHitLists[indx].push_back(&hitlist.fHits[i]);
  }
}
  
//......................................................................

const PartToHitTable::HitList& 
PartToHitTable::operator[](unsigned int i) const 
{
  assert(i<fHitLists.size());
  return fHitLists[i];
}

---