#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	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () 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 ( edm::ParameterSet const & conf )

explicit

Definition at line 31 of file HFRecoEcalCandidateProducer.cc.

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

                                                                                    :
   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>();
 
 } 

Member Function Documentation

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

override

Definition at line 59 of file HFRecoEcalCandidateProducer.cc.

References algo_, edm::Event::getByToken(), hfclustersHFEM_, hfclustersSC_, HFDBversion_, primaryVertexFilter_cfi::maxAbsZ, primaryVertexFilter_cfi::maxd0, eostools::move(), HFRecoEcalCandidateAlgo::produce(), edm::Handle< T >::product(), edm::Event::put(), electrons_cff::vertices, and vertices_.

                                                                                    {  
   
   
   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));
 
 }