CombinedHitPairGenerator.cc

Go to the documentation of this file.

#include "RecoTracker/TkHitPairs/interface/CombinedHitPairGenerator.h"
#include "RecoTracker/TkHitPairs/interface/HitPairGeneratorFromLayerPair.h"
#include "FWCore/Framework/interface/ConsumesCollector.h"
#include "FWCore/Framework/interface/Event.h"
#include "DataFormats/Common/interface/Handle.h"

CombinedHitPairGenerator::CombinedHitPairGenerator(const edm::ParameterSet& cfg, edm::ConsumesCollector& iC)
    : theSeedingLayerToken(iC.consumes<SeedingLayerSetsHits>(cfg.getParameter<edm::InputTag>("SeedingLayers"))) {
  theMaxElement = cfg.getParameter<unsigned int>("maxElement");
  theGenerator = std::make_unique<HitPairGeneratorFromLayerPair>(0, 1, &theLayerCache, theMaxElement);
}

CombinedHitPairGenerator::CombinedHitPairGenerator(const CombinedHitPairGenerator& cb)
    : HitPairGenerator(cb),
      theSeedingLayerToken(cb.theSeedingLayerToken),
      theGenerator(std::make_unique<HitPairGeneratorFromLayerPair>(0, 1, &theLayerCache, cb.theMaxElement)) {
  theMaxElement = cb.theMaxElement;
}

CombinedHitPairGenerator::~CombinedHitPairGenerator() {}

void CombinedHitPairGenerator::hitPairs(const TrackingRegion& region,
                                        OrderedHitPairs& result,
                                        const edm::Event& ev,
                                        const edm::EventSetup& es) {
  edm::Handle<SeedingLayerSetsHits> hlayers;
  ev.getByToken(theSeedingLayerToken, hlayers);
  const SeedingLayerSetsHits& layers = *hlayers;
  if (layers.numberOfLayersInSet() != 2)
    throw cms::Exception("Configuration")
        << "CombinedHitPairGenerator expects SeedingLayerSetsHits::numberOfLayersInSet() to be 2, got "
        << layers.numberOfLayersInSet();

  for (SeedingLayerSetsHits::SeedingLayerSet layerSet : layers) {
    theGenerator->hitPairs(region, result, ev, es, layerSet);
  }

  theLayerCache.clear();

  LogDebug("CombinedHitPairGenerator") << " total number of pairs provided back CHPG : " << result.size();
}