#include <HFRecoEcalCandidateProducer.h>

Inheritance diagram for HFRecoEcalCandidateProducer:

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

void	produce (edm::Event &e, edm::EventSetup const &iSetup) override

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

Private Attributes
HFRecoEcalCandidateAlgo	algo_

double	Cut2D_

std::vector< double >	defaultDB_

double	defaultSlope2D_

bool	doPU_

edm::EDGetToken	hfclustersHFEM_

edm::EDGetToken	hfclustersSC_

std::vector< double >	HFDBvector_

int	HFDBversion_

reco::HFValueStruct	hfvars_

edm::EDGetToken	vertices_

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

Definition at line 24 of file HFRecoEcalCandidateProducer.h.

Constructor & Destructor Documentation

◆ HFRecoEcalCandidateProducer()

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

explicit

Definition at line 31 of file HFRecoEcalCandidateProducer.cc.

     : defaultDB_(std::vector<double>()),
       hfclustersSC_(consumes<reco::SuperClusterCollection>(conf.getParameter<edm::InputTag>("hfclusters"))),
       hfclustersHFEM_(
           consumes<reco::HFEMClusterShapeAssociationCollection>(conf.getParameter<edm::InputTag>("hfclusters"))),
  
       HFDBversion_(conf.existsAs<int>("HFDBversion") ? conf.getParameter<int>("HFDBversion") : 99),  //do nothing
       HFDBvector_(conf.existsAs<std::vector<double> >("HFDBvector")
                       ? conf.getParameter<std::vector<double> >("HFDBvector")
                       : defaultDB_),
       doPU_(false),
       Cut2D_(conf.getParameter<double>("intercept2DCut")),
       defaultSlope2D_(
           (Cut2D_ <= 0.83)
               ? (0.475)
               : ((Cut2D_ > 0.83 && Cut2D_ <= 0.9) ? (0.275) : (0.2))),  //fix for hlt unable to add slope variable now
       hfvars_(HFDBversion_, HFDBvector_),
       algo_(conf.existsAs<bool>("Correct") ? conf.getParameter<bool>("Correct") : true,
             conf.getParameter<double>("e9e25Cut"),
             conf.getParameter<double>("intercept2DCut"),
             conf.existsAs<double>("intercept2DSlope") ? conf.getParameter<double>("intercept2DSlope") : defaultSlope2D_,
             conf.getParameter<std::vector<double> >("e1e9Cut"),
             conf.getParameter<std::vector<double> >("eCOREe9Cut"),
             conf.getParameter<std::vector<double> >("eSeLCut"),
             hfvars_) {
   if (conf.existsAs<edm::InputTag>("VertexCollection")) {
     vertices_ = consumes<reco::VertexCollection>(conf.getParameter<edm::InputTag>("VertexCollection"));
   } else
     vertices_ = consumes<reco::VertexCollection>(edm::InputTag("offlinePrimaryVertices"));
  
   produces<reco::RecoEcalCandidateCollection>();
 }

References edm::ParameterSet::existsAs(), edm::ParameterSet::getParameter(), and vertices_.

Member Function Documentation

◆ produce()

void HFRecoEcalCandidateProducer::produce	(	edm::Event &	e,
		edm::EventSetup const &	iSetup
	)

override

Definition at line 64 of file HFRecoEcalCandidateProducer.cc.

                                                                                   {
   edm::Handle<reco::SuperClusterCollection> super_clus;
   edm::Handle<reco::HFEMClusterShapeAssociationCollection> hf_assoc;
  
   e.getByToken(hfclustersSC_, super_clus);
   e.getByToken(hfclustersHFEM_, hf_assoc);
  
   int nvertex = 0;
   if (HFDBversion_ != 99) {
     edm::Handle<reco::VertexCollection> pvHandle;
     e.getByToken(vertices_, pvHandle);
     const reco::VertexCollection& vertices = *pvHandle.product();
     static const int minNDOF = 4;
     static const double maxAbsZ = 15.0;
     static const double maxd0 = 2.0;
  
     //count verticies
  
     for (reco::VertexCollection::const_iterator vit = vertices.begin(); vit != vertices.end(); ++vit) {
       if (vit->ndof() > minNDOF && ((maxAbsZ <= 0) || fabs(vit->z()) <= maxAbsZ) &&
           ((maxd0 <= 0) || fabs(vit->position().rho()) <= maxd0))
         nvertex++;
     }
   } else {
     nvertex = 1;
   }
  
   // create return data
   auto retdata1 = std::make_unique<reco::RecoEcalCandidateCollection>();
  
   algo_.produce(super_clus, *hf_assoc, *retdata1, nvertex);
  
   e.put(std::move(retdata1));
 }