hitTripletEDProducerT::Impl< T_HitTripletGenerator, T_SeedingHitSets, T_IntermediateHitTriplets > Class Template Reference

#include <HitTripletEDProducerT.h>

Inheritance diagram for hitTripletEDProducerT::Impl< T_HitTripletGenerator, T_SeedingHitSets, T_IntermediateHitTriplets >:

Public Member Functions
	Impl (const edm::ParameterSet &iConfig, edm::ConsumesCollector &iC)
void	produce (const IntermediateHitDoublets &regionDoublets, edm::Event &iEvent, const edm::EventSetup &iSetup) override
void	produces (edm::ProducesCollector producesCollector) const override
	~Impl () override=default
Public Member Functions inherited from hitTripletEDProducerT::ImplGeneratorBase< T_HitTripletGenerator >
	ImplGeneratorBase (const edm::ParameterSet &iConfig, edm::ConsumesCollector &iC)
	~ImplGeneratorBase () override=default
Public Member Functions inherited from hitTripletEDProducerT::ImplBase
	ImplBase ()=default
virtual	~ImplBase ()=default

Additional Inherited Members
Protected Attributes inherited from hitTripletEDProducerT::ImplGeneratorBase< T_HitTripletGenerator >
T_HitTripletGenerator	generator_
Protected Attributes inherited from hitTripletEDProducerT::ImplBase
edm::RunningAverage	localRA_

Detailed Description

template<typename T_HitTripletGenerator, typename T_SeedingHitSets, typename T_IntermediateHitTriplets>
class hitTripletEDProducerT::Impl< T_HitTripletGenerator, T_SeedingHitSets, T_IntermediateHitTriplets >

Definition at line 71 of file HitTripletEDProducerT.h.

Constructor & Destructor Documentation

Impl()

template<typename T_HitTripletGenerator , typename T_SeedingHitSets , typename T_IntermediateHitTriplets >

hitTripletEDProducerT::Impl< T_HitTripletGenerator, T_SeedingHitSets, T_IntermediateHitTriplets >::Impl ( const edm::ParameterSet &  iConfig, edm::ConsumesCollector &  iC  )

inline

Definition at line 73 of file HitTripletEDProducerT.h.

        : ImplGeneratorBase<T_HitTripletGenerator>(iConfig, iC) {}

~Impl()

template<typename T_HitTripletGenerator , typename T_SeedingHitSets , typename T_IntermediateHitTriplets >

hitTripletEDProducerT::Impl< T_HitTripletGenerator, T_SeedingHitSets, T_IntermediateHitTriplets >::~Impl ( )

overridedefault

Member Function Documentation

produce()

template<typename T_HitTripletGenerator , typename T_SeedingHitSets , typename T_IntermediateHitTriplets >

void hitTripletEDProducerT::Impl< T_HitTripletGenerator, T_SeedingHitSets, T_IntermediateHitTriplets >::produce ( const IntermediateHitDoublets &  regionDoublets, edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

inlineoverridevirtual

Implements hitTripletEDProducerT::ImplBase.

Definition at line 82 of file HitTripletEDProducerT.h.

References a, b, IntermediateHitDoublets::empty(), Exception, JetChargeProducer_cfi::exp, LayerHitMapCache::extend(), newFWLiteAna::found, hitTripletEDProducerT::ImplGeneratorBase< T_HitTripletGenerator >::generator_, iEvent, MainPageGenerator::l, IntermediateHitDoublets::layerPairsSize(), LayerTriplets::layers(), hitTripletEDProducerT::ImplBase::localRA_, LogDebug, LogTrace, nano_mu_digi_cff::region, IntermediateHitDoublets::regionSize(), IntermediateHitDoublets::seedingLayerHits(), OrderedHitTriplets::size(), edm::RunningAverage::update(), and edm::RunningAverage::upper().

                                                       {
      const SeedingLayerSetsHits& seedingLayerHits = regionDoublets.seedingLayerHits();

      auto seedingHitSetsProducer = T_SeedingHitSets();
      auto intermediateHitTripletsProducer = T_IntermediateHitTriplets(&seedingLayerHits);

      if (regionDoublets.empty()) {
        seedingHitSetsProducer.putEmpty(iEvent);
        intermediateHitTripletsProducer.putEmpty(iEvent);
        return;
      }

      seedingHitSetsProducer.reserve(regionDoublets.regionSize(), this->localRA_.upper());
      intermediateHitTripletsProducer.reserve(regionDoublets.regionSize(), this->localRA_.upper());

      // match-making of pair and triplet layers
      std::vector<LayerTriplets::LayerSetAndLayers> trilayers = LayerTriplets::layers(seedingLayerHits);

      OrderedHitTriplets triplets;
      triplets.reserve(this->localRA_.upper());
      size_t triplets_total = 0;

      LogDebug("HitTripletEDProducer") << "Creating triplets for " << regionDoublets.regionSize() << " regions, and "
                                       << trilayers.size() << " pair+3rd layers from "
                                       << regionDoublets.layerPairsSize() << " layer pairs";

      for (const auto& regionLayerPairs : regionDoublets) {
        const TrackingRegion& region = regionLayerPairs.region();

        auto hitCachePtr_filler_shs = seedingHitSetsProducer.beginRegion(&region, nullptr);
        auto hitCachePtr_filler_iht =
            intermediateHitTripletsProducer.beginRegion(&region, std::get<0>(hitCachePtr_filler_shs));
        auto hitCachePtr = std::get<0>(hitCachePtr_filler_iht);

        LayerHitMapCache& hitCache = *hitCachePtr;
        hitCache.extend(regionLayerPairs.layerHitMapCache());

        LogTrace("HitTripletEDProducer") << " starting region";

        for (const auto& layerPair : regionLayerPairs) {
          LogTrace("HitTripletEDProducer")
              << "  starting layer pair " << layerPair.innerLayerIndex() << "," << layerPair.outerLayerIndex();

          auto found = std::find_if(trilayers.begin(), trilayers.end(), [&](const LayerTriplets::LayerSetAndLayers& a) {
            return a.first[0].index() == layerPair.innerLayerIndex() &&
                   a.first[1].index() == layerPair.outerLayerIndex();
          });
          if (found == trilayers.end()) {
            auto exp = cms::Exception("LogicError") << "Did not find the layer pair from vector<pair+third layers>. "
                                                       "This is a sign of some internal inconsistency\n";
            exp << "I was looking for layer pair " << layerPair.innerLayerIndex() << "," << layerPair.outerLayerIndex()
                << ". Triplets have the following pairs:\n";
            for (const auto& a : trilayers) {
              exp << " " << a.first[0].index() << "," << a.first[1].index() << ": 3rd layers";
              for (const auto& b : a.second) {
                exp << " " << b.index();
              }
              exp << "\n";
            }
            throw exp;
          }
          const auto& thirdLayers = found->second;

          this->generator_.hitTriplets(region,
                                       triplets,
                                       iEvent,
                                       iSetup,
                                       layerPair.doublets(),
                                       thirdLayers,
                                       intermediateHitTripletsProducer.tripletLastLayerIndexVector(),
                                       hitCache);

#ifdef EDM_ML_DEBUG
          LogTrace("HitTripletEDProducer")
              << "  created " << triplets.size() << " triplets for layer pair " << layerPair.innerLayerIndex() << ","
              << layerPair.outerLayerIndex() << " and 3rd layers";
          for (const auto& l : thirdLayers) {
            LogTrace("HitTripletEDProducer") << "   " << l.index();
          }
#endif

          triplets_total += triplets.size();
          seedingHitSetsProducer.fill(std::get<1>(hitCachePtr_filler_shs), triplets);
          intermediateHitTripletsProducer.fill(
              std::get<1>(hitCachePtr_filler_iht), layerPair.layerPair(), thirdLayers, triplets);

          triplets.clear();
        }
      }
      this->localRA_.update(triplets_total);

      seedingHitSetsProducer.put(iEvent);
      intermediateHitTripletsProducer.put(iEvent);
    }

produces()

template<typename T_HitTripletGenerator , typename T_SeedingHitSets , typename T_IntermediateHitTriplets >

void hitTripletEDProducerT::Impl< T_HitTripletGenerator, T_SeedingHitSets, T_IntermediateHitTriplets >::produces ( edm::ProducesCollector  producesCollector ) const

inlineoverridevirtual

Implements hitTripletEDProducerT::ImplBase.

Definition at line 77 of file HitTripletEDProducerT.h.

                                                                         {
      T_SeedingHitSets::produces(producesCollector);
      T_IntermediateHitTriplets::produces(producesCollector);
    };