HLTRecHitInAllL1RegionsProducer< RecHitType > Class Template Reference

#include <HLTRecHitInAllL1RegionsProducer.h>

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	hasAbilityToProduceInLumis () const final
bool	hasAbilityToProduceInRuns () 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 85 of file HLTRecHitInAllL1RegionsProducer.h.

Member Typedef Documentation

template<typename RecHitType >

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

private

Definition at line 87 of file HLTRecHitInAllL1RegionsProducer.h.

Constructor & Destructor Documentation

template<typename RecHitType >

HLTRecHitInAllL1RegionsProducer< RecHitType >::HLTRecHitInAllL1RegionsProducer

(

const edm::ParameterSet &

ps

)

Definition at line 112 of file HLTRecHitInAllL1RegionsProducer.h.

References edm::ParameterSet::getParameter(), 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 92 of file HLTRecHitInAllL1RegionsProducer.h.

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 248 of file HLTRecHitInAllL1RegionsProducer.h.

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 128 of file HLTRecHitInAllL1RegionsProducer.h.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), edm::ParameterSet::addParameter(), edm::ParameterSetDescription::addVPSet(), defaultModuleLabel(), 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 195 of file HLTRecHitInAllL1RegionsProducer.h.

References event(), edm::EventSetup::get(), CaloSubdetectorGeometry::getGeometry(), CaloGeometry::getSubdetectorGeometry(), edm::HandleBase::id(), input, edm::HandleBase::isValid(), eostools::move(), dataset::name, rpcPointValidation_cfi::recHit, 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 100 of file HLTRecHitInAllL1RegionsProducer.h.

template<typename RecHitType >

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

private

Definition at line 103 of file HLTRecHitInAllL1RegionsProducer.h.

template<typename RecHitType >

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

private

Definition at line 102 of file HLTRecHitInAllL1RegionsProducer.h.

template<typename RecHitType >

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

private

Definition at line 105 of file HLTRecHitInAllL1RegionsProducer.h.