HLTRecHitInAllL1RegionsProducer< RecHitType > Class Template Reference

Inheritance diagram for HLTRecHitInAllL1RegionsProducer< RecHitType >:

Public Member Functions
	HLTRecHitInAllL1RegionsProducer (const edm::ParameterSet &ps)
void	produce (edm::Event &, const edm::EventSetup &) override
	~HLTRecHitInAllL1RegionsProducer () override
Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndRuns () const final

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

Private Types
using	RecHitCollectionType = edm::SortedCollection< RecHitType >

Private Member Functions
L1RegionDataBase *	createL1RegionData (const std::string &, const edm::ParameterSet &, edm::ConsumesCollector &&)

Private Attributes
std::vector< std::unique_ptr< L1RegionDataBase > >	l1RegionData_
std::vector< std::string >	productLabels_
std::vector< edm::InputTag >	recHitLabels_
std::vector< edm::EDGetTokenT< RecHitCollectionType > >	recHitTokens_

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

Detailed Description

template<typename RecHitType>
class HLTRecHitInAllL1RegionsProducer< RecHitType >

Definition at line 94 of file HLTRecHitInAllL1RegionsProducer.cc.

Member Typedef Documentation

template<typename RecHitType >

using HLTRecHitInAllL1RegionsProducer< RecHitType >::RecHitCollectionType = edm::SortedCollection<RecHitType>

private

Definition at line 95 of file HLTRecHitInAllL1RegionsProducer.cc.

Constructor & Destructor Documentation

template<typename RecHitType >

HLTRecHitInAllL1RegionsProducer< RecHitType >::HLTRecHitInAllL1RegionsProducer

(

const edm::ParameterSet &

ps

)

Definition at line 118 of file HLTRecHitInAllL1RegionsProducer.cc.

References edm::ParameterSet::getParameter(), HLT_2018_cff::l1InputRegions, muonDTDigis_cfi::pset, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                                        {
  const std::vector<edm::ParameterSet> l1InputRegions =
      para.getParameter<std::vector<edm::ParameterSet>>("l1InputRegions");
  for (auto& pset : l1InputRegions) {
    const std::string type = pset.getParameter<std::string>("type");
    l1RegionData_.emplace_back(createL1RegionData(
        type,
        pset,
        consumesCollector()));  //meh I was going to use a factory but it was going to be overly complex for my needs
  }
  recHitLabels_ = para.getParameter<std::vector<edm::InputTag>>("recHitLabels");
  productLabels_ = para.getParameter<std::vector<std::string>>("productLabels");

  for (unsigned int collNr = 0; collNr < recHitLabels_.size(); collNr++) {
    recHitTokens_.push_back(consumes<RecHitCollectionType>(recHitLabels_[collNr]));
    produces<RecHitCollectionType>(productLabels_[collNr]);
  }
}

template<typename RecHitType >

HLTRecHitInAllL1RegionsProducer< RecHitType >::~HLTRecHitInAllL1RegionsProducer ( )

inlineoverride

Definition at line 99 of file HLTRecHitInAllL1RegionsProducer.cc.

References fillDescriptions(), and AlCaHLTBitMon_QueryRunRegistry::string.

{}

Member Function Documentation

template<typename RecHitType >

L1RegionDataBase * HLTRecHitInAllL1RegionsProducer< RecHitType >::createL1RegionData ( const std::string &  type, const edm::ParameterSet &  para, edm::ConsumesCollector &&  consumesColl  )

private

Definition at line 249 of file HLTRecHitInAllL1RegionsProducer.cc.

References Exception.

                                                                                             {
  if (type == "L1EmParticle") {
    return new L1RegionData<l1extra::L1EmParticleCollection>(para, consumesColl);
  } else if (type == "L1JetParticle") {
    return new L1RegionData<l1extra::L1JetParticleCollection>(para, consumesColl);
  } else if (type == "L1MuonParticle") {
    return new L1RegionData<l1extra::L1MuonParticleCollection>(para, consumesColl);
  } else if (type == "EGamma") {
    return new L1RegionData<l1t::EGammaBxCollection>(para, consumesColl);
  } else if (type == "Jet") {
    return new L1RegionData<l1t::JetBxCollection>(para, consumesColl);
  } else if (type == "Muon") {
    return new L1RegionData<l1t::MuonBxCollection>(para, consumesColl);
  } else if (type == "Tau") {
    return new L1RegionData<l1t::TauBxCollection>(para, consumesColl);
  } else {
    //this is a major issue and could lead to rather subtle efficiency losses, so if its incorrectly configured, we're aborting the job!
    throw cms::Exception("InvalidConfig")
        << " type " << type
        << " is not recognised, this means the rec-hit you think you are keeping may not be and you should fix this "
           "error as it can lead to hard to find efficiency loses"
        << std::endl;
  }
}

template<typename RecHitType >

void HLTRecHitInAllL1RegionsProducer< RecHitType >::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions )

static

Definition at line 137 of file HLTRecHitInAllL1RegionsProducer.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), edm::ParameterSet::addParameter(), edm::ParameterSetDescription::addVPSet(), defaultModuleLabel(), HLT_2018_cff::InputTag, JetMonitor_cfi::jetPSet, HLT_2018_cff::l1InputRegions, HLT_2018_cff::recHitLabels, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                                             {
  edm::ParameterSetDescription desc;
  std::vector<std::string> productTags;
  productTags.push_back("EcalRegionalRecHitsEB");
  productTags.push_back("EcalRegionalRecHitsEE");
  desc.add<std::vector<std::string>>("productLabels", productTags);
  std::vector<edm::InputTag> recHitLabels;
  recHitLabels.push_back(edm::InputTag("hltEcalRegionalEgammaRecHit:EcalRecHitsEB"));
  recHitLabels.push_back(edm::InputTag("hltEcalRegionalEgammaRecHit:EcalRecHitsEE"));
  recHitLabels.push_back(edm::InputTag("hltESRegionalEgammaRecHit:EcalRecHitsES"));
  desc.add<std::vector<edm::InputTag>>("recHitLabels", recHitLabels);
  std::vector<edm::ParameterSet> l1InputRegions;

  edm::ParameterSet emIsoPSet;
  emIsoPSet.addParameter<std::string>("type", "L1EmParticle");
  emIsoPSet.addParameter<double>("minEt", 5);
  emIsoPSet.addParameter<double>("maxEt", 999);
  emIsoPSet.addParameter<double>("regionEtaMargin", 0.14);
  emIsoPSet.addParameter<double>("regionPhiMargin", 0.4);
  emIsoPSet.addParameter<edm::InputTag>("inputColl", edm::InputTag("hltL1extraParticles:NonIsolated"));
  l1InputRegions.push_back(emIsoPSet);
  edm::ParameterSet emNonIsoPSet;
  emNonIsoPSet.addParameter<std::string>("type", "L1EmParticle");
  emNonIsoPSet.addParameter<double>("minEt", 5);
  emNonIsoPSet.addParameter<double>("maxEt", 999);
  emNonIsoPSet.addParameter<double>("regionEtaMargin", 0.14);
  emNonIsoPSet.addParameter<double>("regionPhiMargin", 0.4);
  emNonIsoPSet.addParameter<edm::InputTag>("inputColl", edm::InputTag("hltL1extraParticles:Isolated"));
  l1InputRegions.push_back(emNonIsoPSet);

  // Why no Central Jets here? They are present in the python config, e.g. OnLine_HLT_GRun.py
  // SHarper: because these are the default parameters designed to reproduce the original (no jets) behaviour
  //
  edm::ParameterSet egPSet;
  egPSet.addParameter<std::string>("type", "EGamma");
  egPSet.addParameter<double>("minEt", 5);
  egPSet.addParameter<double>("maxEt", 999);
  egPSet.addParameter<double>("regionEtaMargin", 0.4);
  egPSet.addParameter<double>("regionPhiMargin", 0.5);
  egPSet.addParameter<edm::InputTag>("inputColl", edm::InputTag("hltCaloStage2Digis"));
  l1InputRegions.push_back(egPSet);

  edm::ParameterSet jetPSet;
  jetPSet.addParameter<std::string>("type", "EGamma");
  jetPSet.addParameter<double>("minEt", 200);
  jetPSet.addParameter<double>("maxEt", 999);
  jetPSet.addParameter<double>("regionEtaMargin", 0.4);
  jetPSet.addParameter<double>("regionPhiMargin", 0.5);
  jetPSet.addParameter<edm::InputTag>("inputColl", edm::InputTag("hltCaloStage2Digis"));
  l1InputRegions.push_back(jetPSet);

  edm::ParameterSetDescription l1InputRegionDesc;
  l1InputRegionDesc.add<std::string>("type");
  l1InputRegionDesc.add<double>("minEt");
  l1InputRegionDesc.add<double>("maxEt");
  l1InputRegionDesc.add<double>("regionEtaMargin");
  l1InputRegionDesc.add<double>("regionPhiMargin");
  l1InputRegionDesc.add<edm::InputTag>("inputColl");
  desc.addVPSet("l1InputRegions", l1InputRegionDesc, l1InputRegions);

  descriptions.add(defaultModuleLabel<HLTRecHitInAllL1RegionsProducer<RecHitType>>(), desc);
}

template<typename RecHitType >

void HLTRecHitInAllL1RegionsProducer< RecHitType >::produce ( edm::Event &  event, const edm::EventSetup &  setup  )

override

Definition at line 201 of file HLTRecHitInAllL1RegionsProducer.cc.

References event(), edm::EventSetup::get(), CaloSubdetectorGeometry::getGeometry(), CaloGeometry::getSubdetectorGeometry(), edm::HandleBase::id(), input, edm::HandleBase::isValid(), eostools::move(), Skims_PA_cff::name, rpcPointValidation_cfi::recHit, FastTrackerRecHitMaskProducer_cfi::recHits, HLT_2018_cff::region, and GeomDetEnumerators::subDetGeom.

                                                                                                     {
  // get the collection geometry:
  edm::ESHandle<CaloGeometry> caloGeomHandle;
  setup.get<CaloGeometryRecord>().get(caloGeomHandle);

  // Get the CaloGeometry
  edm::ESHandle<L1CaloGeometry> l1CaloGeom;
  setup.get<L1CaloGeometryRecord>().get(l1CaloGeom);

  std::vector<RectangularEtaPhiRegion> regions;
  std::for_each(l1RegionData_.begin(),
                l1RegionData_.end(),
                [&event, &regions, l1CaloGeom](const std::unique_ptr<L1RegionDataBase>& input) {
                  input->getEtaPhiRegions(event, regions, *l1CaloGeom);
                });

  for (size_t recHitCollNr = 0; recHitCollNr < recHitTokens_.size(); recHitCollNr++) {
    edm::Handle<RecHitCollectionType> recHits;
    event.getByToken(recHitTokens_[recHitCollNr], recHits);

    if (!(recHits.isValid())) {
      edm::LogError("ProductNotFound") << "could not get a handle on the " << typeid(RecHitCollectionType).name()
                                       << " named " << recHitLabels_[recHitCollNr].encode() << std::endl;
      continue;
    }

    auto filteredRecHits = std::make_unique<RecHitCollectionType>();

    if (!recHits->empty()) {
      const CaloSubdetectorGeometry* subDetGeom = caloGeomHandle->getSubdetectorGeometry(recHits->front().id());
      if (!regions.empty()) {
        for (const RecHitType& recHit : *recHits) {
          auto this_cell = subDetGeom->getGeometry(recHit.id());
          for (const auto& region : regions) {
            if (region.inRegion(this_cell->etaPos(), this_cell->phiPos())) {
              filteredRecHits->push_back(recHit);
              break;
            }
          }
        }
      }  //end check of empty regions
    }    //end check of empty rec-hits
    //   std::cout <<"putting fileter coll in "<<filteredRecHits->size()<<std::endl;
    event.put(std::move(filteredRecHits), productLabels_[recHitCollNr]);
  }  //end loop over all rec hit collections
}

Member Data Documentation

template<typename RecHitType >

std::vector<std::unique_ptr<L1RegionDataBase> > HLTRecHitInAllL1RegionsProducer< RecHitType >::l1RegionData_

private

Definition at line 109 of file HLTRecHitInAllL1RegionsProducer.cc.

template<typename RecHitType >

std::vector<std::string> HLTRecHitInAllL1RegionsProducer< RecHitType >::productLabels_

private

Definition at line 112 of file HLTRecHitInAllL1RegionsProducer.cc.

template<typename RecHitType >

std::vector<edm::InputTag> HLTRecHitInAllL1RegionsProducer< RecHitType >::recHitLabels_

private

Definition at line 111 of file HLTRecHitInAllL1RegionsProducer.cc.

template<typename RecHitType >

std::vector<edm::EDGetTokenT<RecHitCollectionType> > HLTRecHitInAllL1RegionsProducer< RecHitType >::recHitTokens_

private

Definition at line 114 of file HLTRecHitInAllL1RegionsProducer.cc.