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
	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 Types
using	RecHitCollectionType = edm::SortedCollection< RecHitType >

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

Private Attributes
const edm::ESGetToken < CaloGeometry, CaloGeometryRecord >	caloGeometryToken_
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<>
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

template<typename RecHitType>
class HLTRecHitInAllL1RegionsProducer< RecHitType >

Definition at line 101 of file HLTRecHitInAllL1RegionsProducer.cc.

Member Typedef Documentation

template<typename RecHitType >

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

private

Definition at line 102 of file HLTRecHitInAllL1RegionsProducer.cc.

Constructor & Destructor Documentation

template<typename RecHitType >

HLTRecHitInAllL1RegionsProducer< RecHitType >::HLTRecHitInAllL1RegionsProducer

(

const edm::ParameterSet &

ps

)

Definition at line 127 of file HLTRecHitInAllL1RegionsProducer.cc.

References DeDxTools::esConsumes().

    : caloGeometryToken_{esConsumes()} {
  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");
    // meh I was going to use a factory but it was going to be overly complex for my needs
    l1RegionData_.emplace_back(createL1RegionData(type, pset, consumesCollector()));
  }
  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 106 of file HLTRecHitInAllL1RegionsProducer.cc.

{}

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 253 of file HLTRecHitInAllL1RegionsProducer.cc.

References consumesColl, and 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 146 of file HLTRecHitInAllL1RegionsProducer.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), edm::ParameterSet::addParameter(), edm::ParameterSetDescription::addVPSet(), defaultModuleLabel(), submitPVResolutionJobs::desc, HLT_FULL_cff::InputTag, HLT_FULL_cff::l1InputRegions, HLT_FULL_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 210 of file HLTRecHitInAllL1RegionsProducer.cc.

References edmPickEvents::event, edm::EventSetup::getData(), CaloSubdetectorGeometry::getGeometry(), edm::HandleBase::id(), input, edm::HandleBase::isValid(), eostools::move(), mergeVDriftHistosByStation::name, HLT_FULL_cff::recHits, HLT_FULL_cff::region, GeneralSetup::setup(), and GeomDetEnumerators::subDetGeom.

                                                                                                     {
  // get the collection geometry:
  auto const& caloGeom = setup.getData(caloGeometryToken_);

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

  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 = caloGeom.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 >

const edm::ESGetToken<CaloGeometry, CaloGeometryRecord> HLTRecHitInAllL1RegionsProducer< RecHitType >::caloGeometryToken_

private

Definition at line 123 of file HLTRecHitInAllL1RegionsProducer.cc.

template<typename RecHitType >

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

private

Definition at line 116 of file HLTRecHitInAllL1RegionsProducer.cc.

template<typename RecHitType >

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

private

Definition at line 119 of file HLTRecHitInAllL1RegionsProducer.cc.

template<typename RecHitType >

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

private

Definition at line 118 of file HLTRecHitInAllL1RegionsProducer.cc.

template<typename RecHitType >

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

private

Definition at line 121 of file HLTRecHitInAllL1RegionsProducer.cc.