MultiHitFromChi2EDProducer Class Reference

Inheritance diagram for MultiHitFromChi2EDProducer:

Public Member Functions
	MultiHitFromChi2EDProducer (const edm::ParameterSet &iConfig)
void	produce (edm::Event &iEvent, const edm::EventSetup &iSetup) override
	~MultiHitFromChi2EDProducer () override=default
Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default
	EDProducer (const EDProducer &)=delete
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final
const EDProducer &	operator= (const EDProducer &)=delete

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Private Attributes
edm::EDGetTokenT< IntermediateHitDoublets >	doubletToken_
MultiHitGeneratorFromChi2	generator_
edm::RunningAverage	localRA_

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
using	CacheTypes = CacheContexts< T... >
using	GlobalCache = typename CacheTypes::GlobalCache
using	HasAbility = AbilityChecker< T... >
using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
using	RunCache = typename CacheTypes::RunCache
using	RunContext = RunContextT< RunCache, GlobalCache >
using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Detailed Description

Definition at line 17 of file MultiHitFromChi2EDProducer.cc.

Constructor & Destructor Documentation

MultiHitFromChi2EDProducer()

MultiHitFromChi2EDProducer::MultiHitFromChi2EDProducer

(

const edm::ParameterSet &

iConfig

)

Definition at line 34 of file MultiHitFromChi2EDProducer.cc.

    : doubletToken_(consumes<IntermediateHitDoublets>(iConfig.getParameter<edm::InputTag>("doublets"))),
      generator_(iConfig, consumesCollector()) {
  produces<RegionsSeedingHitSets>();
  produces<edm::OwnVector<BaseTrackerRecHit> >();
}

~MultiHitFromChi2EDProducer()

MultiHitFromChi2EDProducer::~MultiHitFromChi2EDProducer ( )

overridedefault

Member Function Documentation

fillDescriptions()

void MultiHitFromChi2EDProducer::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions )

static

Definition at line 41 of file MultiHitFromChi2EDProducer.cc.

References edm::ConfigurationDescriptions::add(), submitPVResolutionJobs::desc, MultiHitGeneratorFromChi2::fillDescriptions(), MultiHitGeneratorFromChi2::fillDescriptionsLabel(), ProducerED_cfi::InputTag, label, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                            {
  edm::ParameterSetDescription desc;

  desc.add<edm::InputTag>("doublets", edm::InputTag("hitPairEDProducer"));

  MultiHitGeneratorFromChi2::fillDescriptions(desc);

  auto label = MultiHitGeneratorFromChi2::fillDescriptionsLabel() + std::string("EDProducerDefault");
  descriptions.add(label, desc);
}

produce()

void MultiHitFromChi2EDProducer::produce ( edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

override

Definition at line 52 of file MultiHitFromChi2EDProducer.cc.

References a, b, LayerHitMapCache::clear(), doubletToken_, RegionsSeedingHitSets::dummyFiller(), Exception, JetChargeProducer_cfi::exp, LayerHitMapCache::extend(), newFWLiteAna::found, generator_, MultiHitGeneratorFromChi2::hitSets(), iEvent, MultiHitGeneratorFromChi2::initES(), MainPageGenerator::l, LayerTriplets::layers(), localRA_, LogDebug, LogTrace, eostools::move(), SeedingLayerSetsHits::numberOfLayersInSet(), nano_mu_digi_cff::region, HLT_2023v12_cff::seedingHitSets, OrderedMultiHits::size(), edm::RunningAverage::update(), and edm::RunningAverage::upper().

                                                                                      {
  edm::Handle<IntermediateHitDoublets> hdoublets;
  iEvent.getByToken(doubletToken_, hdoublets);
  const auto& regionDoublets = *hdoublets;

  const SeedingLayerSetsHits& seedingLayerHits = regionDoublets.seedingLayerHits();
  if (seedingLayerHits.numberOfLayersInSet() < 3) {
    throw cms::Exception("LogicError")
        << "MultiHitFromChi2EDProducer expects SeedingLayerSetsHits::numberOfLayersInSet() to be >= 3, got "
        << seedingLayerHits.numberOfLayersInSet()
        << ". This is likely caused by a configuration error of this module, HitPairEDProducer, or "
           "SeedingLayersEDProducer.";
  }

  auto seedingHitSets = std::make_unique<RegionsSeedingHitSets>();
  if (regionDoublets.empty()) {
    iEvent.put(std::move(seedingHitSets));
    return;
  }
  seedingHitSets->reserve(regionDoublets.regionSize(), localRA_.upper());
  generator_.initES(iSetup);

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

  OrderedMultiHits multihits;
  multihits.reserve(localRA_.upper());
  std::vector<std::unique_ptr<BaseTrackerRecHit> > refittedHitStorage;
  refittedHitStorage.reserve(localRA_.upper() * 2);

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

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

    auto seedingHitSetsFiller = RegionsSeedingHitSets::dummyFiller();
    seedingHitSetsFiller = seedingHitSets->beginRegion(&region);

    hitCache.clear();
    hitCache.extend(regionLayerPairs.layerHitMapCache());

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

    for (const auto& layerPair : regionLayerPairs) {
      LogTrace("MultiHitFromChi2EDProducer")
          << "  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;

      generator_.hitSets(region, multihits, layerPair.doublets(), thirdLayers, hitCache, refittedHitStorage);

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

      for (const SeedingHitSet& hitSet : multihits) {
        seedingHitSetsFiller.emplace_back(hitSet);
      }
      multihits.clear();
    }
  }
  localRA_.update(seedingHitSets->size());

  auto storage = std::make_unique<edm::OwnVector<BaseTrackerRecHit> >();
  storage->reserve(refittedHitStorage.size());
  for (auto& ptr : refittedHitStorage)
    storage->push_back(ptr.release());

  seedingHitSets->shrink_to_fit();
  storage->shrink_to_fit();
  iEvent.put(std::move(seedingHitSets));
  iEvent.put(std::move(storage));
}

Member Data Documentation

doubletToken_

edm::EDGetTokenT<IntermediateHitDoublets> MultiHitFromChi2EDProducer::doubletToken_

private

Definition at line 27 of file MultiHitFromChi2EDProducer.cc.

Referenced by produce().

generator_

MultiHitGeneratorFromChi2 MultiHitFromChi2EDProducer::generator_

private

Definition at line 31 of file MultiHitFromChi2EDProducer.cc.

Referenced by produce().

localRA_

edm::RunningAverage MultiHitFromChi2EDProducer::localRA_

private

Definition at line 29 of file MultiHitFromChi2EDProducer.cc.

Referenced by produce().