#include <HiSuperClusterProducer.h>

Inheritance diagram for HiSuperClusterProducer:

Public Member Functions
virtual void	endJob ()

	HiSuperClusterProducer (const edm::ParameterSet &ps)

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

	~HiSuperClusterProducer () 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
bool	counterExceeded () const

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

void	outputValidationInfo (reco::SuperClusterCollection &superclusterCollection)

void	produceSuperclustersForECALPart (edm::Event &evt, const edm::EDGetTokenT< reco::BasicClusterCollection > &clustersToken, std::string superclusterColection)

Private Attributes
float	barrelBCEnergyThreshold_

float	barrelEtaSearchRoad_

float	barrelPhiSearchRoad_

std::string	barrelSuperclusterCollection_

HiBremRecoveryClusterAlgo *	bremAlgo_p

bool	doBarrel_

bool	doEndcaps_

edm::EDGetTokenT< reco::BasicClusterCollection >	ebClustersToken_

edm::EDGetTokenT< reco::BasicClusterCollection >	eeClustersToken_

float	endcapBCEnergyThreshold_

float	endcapEtaSearchRoad_

float	endcapPhiSearchRoad_

std::string	endcapSuperclusterCollection_

int	nEvt_

int	nMaxPrintout_

int	noSuperClusters

float	seedTransverseEnergyThreshold_

double	totalE

HiBremRecoveryClusterAlgo::VerbosityLevel	verbosity

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 19 of file HiSuperClusterProducer.h.

Constructor & Destructor Documentation

◆ HiSuperClusterProducer()

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

Definition at line 25 of file HiSuperClusterProducer.cc.

                                                                         {
   // The verbosity level
   std::string verbosityString = ps.getParameter<std::string>("VerbosityLevel");
   if (verbosityString == "DEBUG")
     verbosity = HiBremRecoveryClusterAlgo::pDEBUG;
   else if (verbosityString == "WARNING")
     verbosity = HiBremRecoveryClusterAlgo::pWARNING;
   else if (verbosityString == "INFO")
     verbosity = HiBremRecoveryClusterAlgo::pINFO;
   else
     verbosity = HiBremRecoveryClusterAlgo::pERROR;
  
   endcapSuperclusterCollection_ = ps.getParameter<std::string>("endcapSuperclusterCollection");
   barrelSuperclusterCollection_ = ps.getParameter<std::string>("barrelSuperclusterCollection");
  
   doBarrel_ = ps.getParameter<bool>("doBarrel");
   doEndcaps_ = ps.getParameter<bool>("doEndcaps");
  
   barrelEtaSearchRoad_ = ps.getParameter<double>("barrelEtaSearchRoad");
   barrelPhiSearchRoad_ = ps.getParameter<double>("barrelPhiSearchRoad");
   endcapEtaSearchRoad_ = ps.getParameter<double>("endcapEtaSearchRoad");
   endcapPhiSearchRoad_ = ps.getParameter<double>("endcapPhiSearchRoad");
   seedTransverseEnergyThreshold_ = ps.getParameter<double>("seedTransverseEnergyThreshold");
   barrelBCEnergyThreshold_ = ps.getParameter<double>("barrelBCEnergyThreshold");
   endcapBCEnergyThreshold_ = ps.getParameter<double>("endcapBCEnergyThreshold");
  
   if (verbosityString == "INFO") {
     std::cout << "Barrel BC Energy threshold = " << barrelBCEnergyThreshold_ << std::endl;
     std::cout << "Endcap BC Energy threshold = " << endcapBCEnergyThreshold_ << std::endl;
   }
  
   bremAlgo_p = new HiBremRecoveryClusterAlgo(barrelEtaSearchRoad_,
                                              barrelPhiSearchRoad_,
                                              endcapEtaSearchRoad_,
                                              endcapPhiSearchRoad_,
                                              seedTransverseEnergyThreshold_,
                                              barrelBCEnergyThreshold_,
                                              endcapBCEnergyThreshold_,
                                              verbosity);
  
   produces<reco::SuperClusterCollection>(endcapSuperclusterCollection_);
   produces<reco::SuperClusterCollection>(barrelSuperclusterCollection_);
  
   eeClustersToken_ = consumes<reco::BasicClusterCollection>(edm::InputTag(
       ps.getParameter<std::string>("endcapClusterProducer"), ps.getParameter<std::string>("endcapClusterCollection")));
   ebClustersToken_ = consumes<reco::BasicClusterCollection>(edm::InputTag(
       ps.getParameter<std::string>("barrelClusterProducer"), ps.getParameter<std::string>("barrelClusterCollection")));
  
   totalE = 0;
   noSuperClusters = 0;
   nEvt_ = 0;
 }

References barrelBCEnergyThreshold_, barrelEtaSearchRoad_, barrelPhiSearchRoad_, barrelSuperclusterCollection_, bremAlgo_p, gather_cfg::cout, doBarrel_, doEndcaps_, ebClustersToken_, eeClustersToken_, endcapBCEnergyThreshold_, endcapEtaSearchRoad_, endcapPhiSearchRoad_, endcapSuperclusterCollection_, edm::ParameterSet::getParameter(), nEvt_, noSuperClusters, HiBremRecoveryClusterAlgo::pDEBUG, HiBremRecoveryClusterAlgo::pERROR, HiBremRecoveryClusterAlgo::pINFO, HiBremRecoveryClusterAlgo::pWARNING, seedTransverseEnergyThreshold_, AlCaHLTBitMon_QueryRunRegistry::string, totalE, and verbosity.

◆ ~HiSuperClusterProducer()

HiSuperClusterProducer::~HiSuperClusterProducer ( )

override

Definition at line 78 of file HiSuperClusterProducer.cc.

78 { delete bremAlgo_p; }

References bremAlgo_p.

Member Function Documentation

◆ counterExceeded()

bool HiSuperClusterProducer::counterExceeded ( ) const

inlineprivate

Definition at line 66 of file HiSuperClusterProducer.h.

66 { return ((nEvt_ > nMaxPrintout_) || (nMaxPrintout_ < 0)); }

References nEvt_, and nMaxPrintout_.

◆ endJob()

void HiSuperClusterProducer::endJob ( void )

virtual

Definition at line 80 of file HiSuperClusterProducer.cc.

                                     {
   double averEnergy = 0.;
   std::ostringstream str;
   str << "HiSuperClusterProducer::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("HiSuperClusterProducerInfo") << str.str() << "\n";
 }

References noSuperClusters, str, and totalE.

◆ getClusterPtrVector()

void HiSuperClusterProducer::getClusterPtrVector	(	edm::Event &	evt,
		const edm::EDGetTokenT< reco::BasicClusterCollection > &	clustersToken,
		reco::CaloClusterPtrVector *	clusterPtrVector_p
	)

private

Definition at line 128 of file HiSuperClusterProducer.cc.

                                                                                                {
   edm::Handle<reco::BasicClusterCollection> bccHandle;
  
   evt.getByToken(clustersToken, bccHandle);
  
   if (!(bccHandle.isValid())) {
     edm::LogError("HiSuperClusterProducerError") << "could not get a handle on the BasicCluster Collection!";
     clusterPtrVector_p = nullptr;
   }
  
   const reco::BasicClusterCollection* clusterCollection_p = bccHandle.product();
   for (unsigned int i = 0; i < clusterCollection_p->size(); i++) {
     clusterPtrVector_p->push_back(reco::CaloClusterPtr(bccHandle, i));
   }
 }

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

Referenced by produceSuperclustersForECALPart().

◆ outputValidationInfo()

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

private

◆ produce()

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

override

Definition at line 93 of file HiSuperClusterProducer.cc.

                                                                            {
   if (doEndcaps_)
     produceSuperclustersForECALPart(evt, eeClustersToken_, endcapSuperclusterCollection_);
  
   if (doBarrel_)
     produceSuperclustersForECALPart(evt, ebClustersToken_, barrelSuperclusterCollection_);
  
   nEvt_++;
 }

References barrelSuperclusterCollection_, doBarrel_, doEndcaps_, ebClustersToken_, eeClustersToken_, endcapSuperclusterCollection_, nEvt_, and produceSuperclustersForECALPart().

◆ produceSuperclustersForECALPart()

void HiSuperClusterProducer::produceSuperclustersForECALPart	(	edm::Event &	evt,
		const edm::EDGetTokenT< reco::BasicClusterCollection > &	clustersToken,
		std::string	superclusterColection
	)

private

Definition at line 103 of file HiSuperClusterProducer.cc.

                                       {
   // get the cluster collection out and turn it to a BasicClusterRefVector:
   reco::CaloClusterPtrVector* clusterPtrVector_p = new reco::CaloClusterPtrVector;
   getClusterPtrVector(evt, clustersToken, clusterPtrVector_p);
  
   // run the brem recovery and get the SC collection
   auto superclusters_ap =
       std::make_unique<reco::SuperClusterCollection>(bremAlgo_p->makeSuperClusters(*clusterPtrVector_p));
  
   // count the total energy and the number of superclusters
   reco::SuperClusterCollection::iterator it;
   for (it = superclusters_ap->begin(); it != superclusters_ap->end(); it++) {
     totalE += it->energy();
     noSuperClusters++;
   }
  
   // put the SC collection in the event
   evt.put(std::move(superclusters_ap), superclusterCollection);
  
   delete clusterPtrVector_p;
 }