Inheritance diagram for CorrectedECALPFClusterProducer:

Public Member Functions
	CorrectedECALPFClusterProducer (const edm::ParameterSet &conf)

void	produce (edm::Event &e, const edm::EventSetup &es) 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 Attributes
std::unique_ptr< PFClusterEMEnergyCorrector >	corrector_

edm::EDGetTokenT< reco::PFClusterCollection >	inputECAL_

edm::EDGetTokenT< reco::PFClusterCollection >	inputPS_

const double	minimumPSEnergy_

const bool	skipPS_

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 44 of file CorrectedECALPFClusterProducer.cc.

Constructor & Destructor Documentation

◆ CorrectedECALPFClusterProducer()

CorrectedECALPFClusterProducer::CorrectedECALPFClusterProducer ( const edm::ParameterSet & conf )

inline

Definition at line 46 of file CorrectedECALPFClusterProducer.cc.

       : minimumPSEnergy_(conf.getParameter<double>("minimumPSEnergy")), skipPS_(conf.getParameter<bool>("skipPS")) {
     const edm::InputTag& inputECAL = conf.getParameter<edm::InputTag>("inputECAL");
     inputECAL_ = consumes<reco::PFClusterCollection>(inputECAL);
  
     const edm::InputTag& inputPS = conf.getParameter<edm::InputTag>("inputPS");
     if (!skipPS_)
       inputPS_ = consumes<reco::PFClusterCollection>(inputPS);
  
     const edm::ParameterSet& corConf = conf.getParameterSet("energyCorrector");
     corrector_ = std::make_unique<PFClusterEMEnergyCorrector>(corConf, consumesCollector());
  
     produces<reco::PFCluster::EEtoPSAssociation>();
     produces<reco::PFClusterCollection>();
   }

References corrector_, edm::ParameterSet::getParameter(), edm::ParameterSet::getParameterSet(), HLT_FULL_cff::inputECAL, inputECAL_, HLT_FULL_cff::inputPS, inputPS_, and skipPS_.

Member Function Documentation

◆ fillDescriptions()

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

static

Definition at line 133 of file CorrectedECALPFClusterProducer.cc.

                                                                                                 {
   edm::ParameterSetDescription desc;
   desc.add<double>("minimumPSEnergy", 0.0);
   desc.ifValue(
       edm::ParameterDescription<bool>("skipPS", false, true),
       true >> (edm::ParameterDescription<edm::InputTag>("inputPS", edm::InputTag(""), true)) or
           false >> (edm::ParameterDescription<edm::InputTag>("inputPS", edm::InputTag("particleFlowClusterPS"), true)));
   {
     edm::ParameterSetDescription psd0;
     psd0.add<bool>("applyCrackCorrections", false);
     psd0.add<bool>("applyMVACorrections", false);
     psd0.add<bool>("srfAwareCorrection", false);
     psd0.add<bool>("setEnergyUncertainty", false);
     psd0.add<bool>("autoDetectBunchSpacing", true);
     psd0.add<int>("bunchSpacing", 25);
     psd0.add<double>("maxPtForMVAEvaluation", -99.);
     psd0.add<edm::InputTag>("recHitsEBLabel", edm::InputTag("ecalRecHit", "EcalRecHitsEB"));
     psd0.add<edm::InputTag>("recHitsEELabel", edm::InputTag("ecalRecHit", "EcalRecHitsEE"));
     psd0.add<edm::InputTag>("ebSrFlagLabel", edm::InputTag("ecalDigis"));
     psd0.add<edm::InputTag>("eeSrFlagLabel", edm::InputTag("ecalDigis"));
     desc.add<edm::ParameterSetDescription>("energyCorrector", psd0);
   }
   desc.add<edm::InputTag>("inputECAL", edm::InputTag("particleFlowClusterECALUncorrected"));
   descriptions.add("particleFlowClusterECAL", desc);
 }

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), submitPVResolutionJobs::desc, HLT_FULL_cff::InputTag, and or.

◆ produce()

void CorrectedECALPFClusterProducer::produce	(	edm::Event &	e,
		const edm::EventSetup &	es
	)

override

Definition at line 75 of file CorrectedECALPFClusterProducer.cc.

                                                                                  {
   auto clusters_out = std::make_unique<reco::PFClusterCollection>();
   auto association_out = std::make_unique<reco::PFCluster::EEtoPSAssociation>();
  
   edm::Handle<reco::PFClusterCollection> handleECAL;
   e.getByToken(inputECAL_, handleECAL);
   edm::Handle<reco::PFClusterCollection> handlePS;
   if (!skipPS_)
     e.getByToken(inputPS_, handlePS);
  
   auto const& ecals = *handleECAL;
  
   clusters_out->reserve(ecals.size());
   association_out->reserve(ecals.size());
   clusters_out->insert(clusters_out->end(), ecals.begin(), ecals.end());
  
   //build the EE->PS association
   if (!skipPS_) {
     auto const& pss = *handlePS;
     for (unsigned i = 0; i < pss.size(); ++i) {
       switch (pss[i].layer()) {  // just in case this isn't the ES...
         case PFLayer::PS1:
         case PFLayer::PS2:
           break;
         default:
           continue;
       }
       if (pss[i].energy() < minimumPSEnergy_)
         continue;
       int eematch = -1;
       auto min_dist = std::numeric_limits<double>::max();
       for (size_t ic = 0; ic < ecals.size(); ++ic) {
         if (ecals[ic].layer() != PFLayer::ECAL_ENDCAP)
           continue;
         auto dist = testPreshowerDistance(ecals[ic], pss[i]);
         if (dist == -1.0)
           dist = std::numeric_limits<double>::max();
         if (dist < min_dist) {
           eematch = ic;
           min_dist = dist;
         }
       }  // loop on EE clusters
       if (eematch >= 0) {
         edm::Ptr<reco::PFCluster> psclus(handlePS, i);
         association_out->push_back(std::make_pair(eematch, psclus));
       }
     }
   }
   std::sort(association_out->begin(), association_out->end(), sortByKey);
  
   corrector_->correctEnergies(e, es, *association_out, *clusters_out);
  
   association_out->shrink_to_fit();
  
   e.put(std::move(association_out));
   e.put(std::move(clusters_out));
 }