#include <Multi5x5SuperClusterProducer.h>

Inheritance diagram for Multi5x5SuperClusterProducer:

Public Member Functions
void	endStream () override

	Multi5x5SuperClusterProducer (const edm::ParameterSet &ps)

void	produce (edm::Event &, const edm::EventSetup &) 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 Member Functions
reco::CaloClusterPtrVector	getClusterPtrVector (edm::Event &evt, const edm::EDGetTokenT< reco::BasicClusterCollection > &clustersToken) const

void	outputValidationInfo (reco::SuperClusterCollection &superclusterCollection)

void	produceSuperclustersForECALPart (edm::Event &evt, const edm::EDGetTokenT< reco::BasicClusterCollection > &clustersToken, const edm::EDPutTokenT< reco::SuperClusterCollection > &putToken)

Private Attributes
const float	barrelEtaSearchRoad_

const float	barrelPhiSearchRoad_

edm::EDPutTokenT< reco::SuperClusterCollection >	barrelPutToken_

std::unique_ptr< Multi5x5BremRecoveryClusterAlgo >	bremAlgo_p

const bool	doBarrel_

const bool	doEndcaps_

edm::EDGetTokenT< reco::BasicClusterCollection >	ebClustersToken_

edm::EDGetTokenT< reco::BasicClusterCollection >	eeClustersToken_

const float	endcapEtaSearchRoad_

const float	endcapPhiSearchRoad_

edm::EDPutTokenT< reco::SuperClusterCollection >	endcapPutToken_

int	noSuperClusters

const float	seedTransverseEnergyThreshold_

double	totalE

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 18 of file Multi5x5SuperClusterProducer.h.

Constructor & Destructor Documentation

◆ Multi5x5SuperClusterProducer()

Multi5x5SuperClusterProducer::Multi5x5SuperClusterProducer ( const edm::ParameterSet & ps )

Definition at line 23 of file Multi5x5SuperClusterProducer.cc.

     : barrelEtaSearchRoad_{static_cast<float>(ps.getParameter<double>("barrelEtaSearchRoad"))},
       barrelPhiSearchRoad_{static_cast<float>(ps.getParameter<double>("barrelPhiSearchRoad"))},
       endcapEtaSearchRoad_{static_cast<float>(ps.getParameter<double>("endcapEtaSearchRoad"))},
       endcapPhiSearchRoad_{static_cast<float>(ps.getParameter<double>("endcapPhiSearchRoad"))},
       seedTransverseEnergyThreshold_{static_cast<float>(ps.getParameter<double>("seedTransverseEnergyThreshold"))},
       doBarrel_{ps.getParameter<bool>("doBarrel")},
       doEndcaps_{ps.getParameter<bool>("doEndcaps")},
       totalE{0.},
       noSuperClusters{0} {
   if (doEndcaps_) {
     eeClustersToken_ = consumes<reco::BasicClusterCollection>(ps.getParameter<edm::InputTag>("endcapClusterTag"));
   }
   if (doBarrel_) {
     ebClustersToken_ = consumes<reco::BasicClusterCollection>(ps.getParameter<edm::InputTag>("barrelClusterTag"));
   }
  
   const edm::ParameterSet bremRecoveryPset = ps.getParameter<edm::ParameterSet>("bremRecoveryPset");
   bool dynamicPhiRoad = ps.getParameter<bool>("dynamicPhiRoad");
  
   bremAlgo_p = std::make_unique<Multi5x5BremRecoveryClusterAlgo>(bremRecoveryPset,
                                                                  barrelEtaSearchRoad_,
                                                                  barrelPhiSearchRoad_,
                                                                  endcapEtaSearchRoad_,
                                                                  endcapPhiSearchRoad_,
                                                                  dynamicPhiRoad,
                                                                  seedTransverseEnergyThreshold_);
  
   if (doEndcaps_) {
     endcapPutToken_ =
         produces<reco::SuperClusterCollection>(ps.getParameter<std::string>("endcapSuperclusterCollection"));
   }
   if (doBarrel_) {
     barrelPutToken_ =
         produces<reco::SuperClusterCollection>(ps.getParameter<std::string>("barrelSuperclusterCollection"));
   }
 }

References edm::ParameterSet::getParameter().

Member Function Documentation

◆ endStream()

void Multi5x5SuperClusterProducer::endStream ( )

override

Definition at line 61 of file Multi5x5SuperClusterProducer.cc.

                                              {
   double averEnergy = 0.;
   std::ostringstream str;
   str << "Multi5x5SuperClusterProducer::endJob()\n"
       << "  total # reconstructed super clusters: " << noSuperClusters << "\n"
       << "  total energy of all clusters: " << totalE << "\n";
   if (noSuperClusters > 0) {
     averEnergy = totalE / noSuperClusters;
     str << "  average SuperCluster energy = " << averEnergy << "\n";
   }
   edm::LogInfo("Multi5x5SuperClusterProducerInfo") << str.str() << "\n";
 }

References noSuperClusters, str, and totalE.

◆ getClusterPtrVector()

reco::CaloClusterPtrVector Multi5x5SuperClusterProducer::getClusterPtrVector	(	edm::Event &	evt,
		const edm::EDGetTokenT< reco::BasicClusterCollection > &	clustersToken
	)		const

private

Definition at line 102 of file Multi5x5SuperClusterProducer.cc.

                                                                                           {
   reco::CaloClusterPtrVector clusterPtrVector_p;
   edm::Handle<reco::BasicClusterCollection> bccHandle;
   evt.getByToken(clustersToken, bccHandle);
  
   const reco::BasicClusterCollection* clusterCollection_p = bccHandle.product();
   clusterPtrVector_p.reserve(clusterCollection_p->size());
   for (unsigned int i = 0; i < clusterCollection_p->size(); i++) {
     clusterPtrVector_p.push_back(reco::CaloClusterPtr(bccHandle, i));
   }
   return clusterPtrVector_p;
 }

References edm::Event::getByToken(), mps_fire::i, edm::Handle< T >::product(), edm::PtrVector< T >::push_back(), and edm::PtrVectorBase::reserve().

Referenced by produceSuperclustersForECALPart().

◆ outputValidationInfo()

void Multi5x5SuperClusterProducer::outputValidationInfo ( reco::SuperClusterCollection & superclusterCollection )

private

◆ produce()

void Multi5x5SuperClusterProducer::produce	(	edm::Event &	evt,
		const edm::EventSetup &	es
	)

override

Definition at line 74 of file Multi5x5SuperClusterProducer.cc.

                                                                                  {
   if (doEndcaps_)
     produceSuperclustersForECALPart(evt, eeClustersToken_, endcapPutToken_);
  
   if (doBarrel_)
     produceSuperclustersForECALPart(evt, ebClustersToken_, barrelPutToken_);
 }

References barrelPutToken_, doBarrel_, doEndcaps_, ebClustersToken_, eeClustersToken_, endcapPutToken_, and produceSuperclustersForECALPart().

◆ produceSuperclustersForECALPart()

void Multi5x5SuperClusterProducer::produceSuperclustersForECALPart	(	edm::Event &	evt,
		const edm::EDGetTokenT< reco::BasicClusterCollection > &	clustersToken,
		const edm::EDPutTokenT< reco::SuperClusterCollection > &	putToken
	)

private

Definition at line 82 of file Multi5x5SuperClusterProducer.cc.

                                                                 {
   // get the cluster collection out and turn it to a BasicClusterRefVector:
   reco::CaloClusterPtrVector clusterPtrVector_p = getClusterPtrVector(evt, clustersToken);
  
   // run the brem recovery and get the SC collection
   reco::SuperClusterCollection superclusters_ap(bremAlgo_p->makeSuperClusters(clusterPtrVector_p));
  
   // count the total energy and the number of superclusters
   for (auto const& sc : superclusters_ap) {
     totalE += sc.energy();
     noSuperClusters++;
   }
  
   // put the SC collection in the event
   evt.emplace(putToken, std::move(superclusters_ap));
 }