#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

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () const final

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

Constructor & Destructor Documentation

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

Definition at line 25 of file HiSuperClusterProducer.cc.

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

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

HiSuperClusterProducer::~HiSuperClusterProducer ( )

override

Definition at line 78 of file HiSuperClusterProducer.cc.

References bremAlgo_p.

78 { delete bremAlgo_p; }

HiSuperClusterProducer::bremAlgo_p

HiBremRecoveryClusterAlgo * bremAlgo_p

Definition: HiSuperClusterProducer.h:51

Member Function Documentation

bool HiSuperClusterProducer::counterExceeded ( ) const

inlineprivate

Definition at line 66 of file HiSuperClusterProducer.h.

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

HiSuperClusterProducer::nMaxPrintout_

int nMaxPrintout_

Definition: HiSuperClusterProducer.h:29

HiSuperClusterProducer::nEvt_

int nEvt_

Definition: HiSuperClusterProducer.h:30

void HiSuperClusterProducer::endJob ( void )

virtual

Definition at line 80 of file HiSuperClusterProducer.cc.

References noSuperClusters, str, and totalE.

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

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.

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

Referenced by produceSuperclustersForECALPart().

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

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

private

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

override

Definition at line 93 of file HiSuperClusterProducer.cc.

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

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

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

private

Definition at line 103 of file HiSuperClusterProducer.cc.

References bremAlgo_p, getClusterPtrVector(), HiBremRecoveryClusterAlgo::makeSuperClusters(), eostools::move(), noSuperClusters, edm::Event::put(), and totalE.

Referenced by produce().

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