StripCompactDigiSimLinks.cc

Go to the documentation of this file.

#include "SimDataFormats/TrackerDigiSimLink/interface/StripCompactDigiSimLinks.h"

#include <algorithm>

StripCompactDigiSimLinks::Links StripCompactDigiSimLinks::getLinks(const StripCompactDigiSimLinks::key_type &key) const {
  std::vector<TrackRecord>::const_iterator last = trackRecords_.end();
  std::vector<TrackRecord>::const_iterator match = std::lower_bound(trackRecords_.begin(), last, key);
  if ((match != last) && (*match == key)) {
    // std::vector<TrackRecord>::const_iterator next = match+1;
    unsigned int end = (match + 1 == last ? hitRecords_.size() : (match + 1)->start);
    return Links(hitRecords_.begin() + match->start, hitRecords_.begin() + end);
  } else {
    return Links(hitRecords_.end(), hitRecords_.end());
  }
}

StripCompactDigiSimLinks::StripCompactDigiSimLinks(const StripCompactDigiSimLinks::Filler &filler) {
  trackRecords_.reserve(filler.keySize());
  hitRecords_.reserve(filler.dataSize());
  for (auto const &pair : filler.storage()) {
    trackRecords_.push_back(TrackRecord(pair.first, hitRecords_.size()));
    hitRecords_.insert(hitRecords_.end(), pair.second.begin(), pair.second.end());
  }
}

StripCompactDigiSimLinks::Filler::~Filler() {}

void StripCompactDigiSimLinks::Filler::insert(const StripCompactDigiSimLinks::key_type &key,
                                              const StripCompactDigiSimLinks::HitRecord &record) {
  Filler::Map::iterator it = storage_.find(key);
  if (it == storage_.end()) {
    storage_.insert(std::make_pair(key, std::vector<HitRecord>(1, record)));
    num_keys_++;
    num_values_++;
  } else {
    it->second.push_back(record);
    num_values_++;
  }
}

std::map<uint32_t, std::vector<StripCompactDigiSimLinks::RevLink> > StripCompactDigiSimLinks::makeReverseMap() const {
  std::map<uint32_t, std::vector<StripCompactDigiSimLinks::RevLink> > ret;
  typedef std::vector<TrackRecord>::const_iterator trk_it;
  typedef std::vector<HitRecord>::const_iterator hit_it;
  hit_it hstart = hitRecords_.begin(), ith = hstart;
  for (trk_it itt = trackRecords_.begin(), endt = trackRecords_.end(); itt != endt; ++itt) {
    hit_it edh = (itt + 1 != endt ? hstart + (itt + 1)->start : hitRecords_.end());
    for (; ith < edh; ++ith) {
      ret[ith->detId].push_back(RevLink(*itt, *ith));
    }
  }
  return ret;
}