Inheritance diagram for l1tpf::PFClusterProducerFromL1EGClusters:

Public Member Functions
	PFClusterProducerFromL1EGClusters (const edm::ParameterSet &)

	~PFClusterProducerFromL1EGClusters () 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

Private Member Functions
void	produce (edm::Event &, const edm::EventSetup &) override

Private Attributes
l1tpf::corrector	corrector_

std::vector< double > const	etaBounds_

double	etCut_

std::vector< unsigned int > const	maxClustersEtaPhi_

std::vector< double > const	phiBounds_

l1tpf::ParametricResolution	resol_

edm::EDGetTokenT< BXVector< l1t::EGamma > >	src_

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 14 of file PFClusterProducerFromL1EGClusters.cc.

Constructor & Destructor Documentation

◆ PFClusterProducerFromL1EGClusters()

l1tpf::PFClusterProducerFromL1EGClusters::PFClusterProducerFromL1EGClusters ( const edm::ParameterSet & iConfig )

explicit

Definition at line 33 of file PFClusterProducerFromL1EGClusters.cc.

References etaBounds_, Exception, maxClustersEtaPhi_, and phiBounds_.

     : src_(consumes<BXVector<l1t::EGamma>>(iConfig.getParameter<edm::InputTag>("src"))),
       etCut_(iConfig.getParameter<double>("etMin")),
       etaBounds_(iConfig.getParameter<std::vector<double>>("etaBounds")),
       phiBounds_(iConfig.getParameter<std::vector<double>>("phiBounds")),
       maxClustersEtaPhi_(iConfig.getParameter<std::vector<unsigned int>>("maxClustersEtaPhi")),
       corrector_(iConfig.getParameter<std::string>("corrector"), -1),
       resol_(iConfig.getParameter<edm::ParameterSet>("resol")) {
   produces<l1t::PFClusterCollection>("all");
   produces<l1t::PFClusterCollection>("selected");
   if ((etaBounds_.size() - 1) * (phiBounds_.size() - 1) != maxClustersEtaPhi_.size()) {
     throw cms::Exception("Configuration")
         << "Size mismatch between eta/phi bounds and max clusters: " << (etaBounds_.size() - 1) << " x "
         << (phiBounds_.size() - 1) << " != " << maxClustersEtaPhi_.size() << "\n";
   }
   if (!std::is_sorted(etaBounds_.begin(), etaBounds_.end())) {
     throw cms::Exception("Configuration") << "etaBounds is not sorted\n";
   }
   if (!std::is_sorted(phiBounds_.begin(), phiBounds_.end())) {
     throw cms::Exception("Configuration") << "phiBounds is not sorted\n";
   }
 }

◆ ~PFClusterProducerFromL1EGClusters()

l1tpf::PFClusterProducerFromL1EGClusters::~PFClusterProducerFromL1EGClusters ( )

inlineoverride

Definition at line 17 of file PFClusterProducerFromL1EGClusters.cc.

17 {}

Member Function Documentation

◆ produce()

void l1tpf::PFClusterProducerFromL1EGClusters::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &
	)

overrideprivate

Definition at line 56 of file PFClusterProducerFromL1EGClusters.cc.

References funct::abs(), bsc_activity_cfg::clusters, PVValHelper::eta, l1tpf_calo::GridSelector::fill(), l1trig_cff::hwQual, iEvent, cuy::ii, dqmdumpme::indices, eostools::move(), MillePedeFileConverter_cfg::out, DiDispStaMuonMonitor_cfi::pt, l1tpf_calo::GridSelector::returnSorted(), and l1t::PFCluster::setPtError().

                                                                                             {
   std::unique_ptr<l1t::PFClusterCollection> out(new l1t::PFClusterCollection());
   std::unique_ptr<l1t::PFClusterCollection> out_sel(new l1t::PFClusterCollection());
   edm::Handle<BXVector<l1t::EGamma>> clusters;
   iEvent.getByToken(src_, clusters);
 
   l1tpf_calo::GridSelector selector = l1tpf_calo::GridSelector(etaBounds_, phiBounds_, maxClustersEtaPhi_);
 
   unsigned int index = 0;
   for (auto it = clusters->begin(), ed = clusters->end(); it != ed; ++it, ++index) {
     if (it->pt() <= etCut_)
       continue;
 
     l1t::PFCluster cluster(
         it->pt(), it->eta(), it->phi(), /*hOverE=*/0., /*isEM=*/true);  // it->hovere() seems to return random values
     if (corrector_.valid())
       corrector_.correctPt(cluster);
     cluster.setPtError(resol_(cluster.pt(), std::abs(cluster.eta())));
     cluster.setHwQual(it->hwQual());
     out->push_back(cluster);
     out->back().addConstituent(edm::Ptr<l1t::L1Candidate>(clusters, index));
     selector.fill(cluster.pt(), cluster.eta(), cluster.phi(), index);
   }
   std::vector<unsigned int> indices = selector.returnSorted();
   for (unsigned int ii = 0; ii < indices.size(); ii++) {
     unsigned int theIndex = indices[ii];
     l1t::PFCluster cluster((clusters->begin() + theIndex)->pt(),
                            (clusters->begin() + theIndex)->eta(),
                            (clusters->begin() + theIndex)->phi(),
                            /*hOverE=*/0.,
                            /*isEM=*/true);  // it->hovere() seems to return random values
     if (corrector_.valid())
       corrector_.correctPt(cluster);
     cluster.setPtError(resol_(cluster.pt(), std::abs(cluster.eta())));
     cluster.setHwQual((clusters->begin() + theIndex)->hwQual());
     out_sel->push_back(cluster);
     out_sel->back().addConstituent(edm::Ptr<l1t::L1Candidate>(clusters, theIndex));
   }
 
   iEvent.put(std::move(out), "all");
   iEvent.put(std::move(out_sel), "selected");
 }