#include <Multi5x5ClusterProducer.h>

Inheritance diagram for Multi5x5ClusterProducer:

Public Member Functions
	Multi5x5ClusterProducer (const edm::ParameterSet &ps)

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

	~Multi5x5ClusterProducer () 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
void	clusterizeECALPart (edm::Event &evt, const edm::EventSetup &es, const edm::EDGetTokenT< EcalRecHitCollection > &token, const std::string &clusterCollection, const reco::CaloID::Detectors detector)

bool	counterExceeded () const

const EcalRecHitCollection *	getCollection (edm::Event &evt, const edm::EDGetTokenT< EcalRecHitCollection > &token)

void	outputValidationInfo (reco::CaloClusterPtrVector &clusterPtrVector)

Private Attributes
std::string	barrelClusterCollection_

edm::EDGetTokenT< EcalRecHitCollection >	barrelHitToken_

edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	caloGeometryToken_

bool	doBarrel_

bool	doEndcap_

std::string	endcapClusterCollection_

edm::EDGetTokenT< EcalRecHitCollection >	endcapHitToken_

Multi5x5ClusterAlgo *	island_p

int	nEvt_

int	nMaxPrintout_

PositionCalc	posCalculator_

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 27 of file Multi5x5ClusterProducer.h.

Constructor & Destructor Documentation

◆ Multi5x5ClusterProducer()

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

Definition at line 35 of file Multi5x5ClusterProducer.cc.

                                                                           {
   // Parameters to identify the hit collections
   barrelHitToken_ = consumes<EcalRecHitCollection>(ps.getParameter<edm::InputTag>("barrelHitTag"));
  
   endcapHitToken_ = consumes<EcalRecHitCollection>(ps.getParameter<edm::InputTag>("endcapHitTag"));
  
   //EventSetup Token for CaloGeometry
   caloGeometryToken_ = esConsumes<CaloGeometry, CaloGeometryRecord>();
  
   // should cluster algo be run in barrel and endcap?
   doEndcap_ = ps.getParameter<bool>("doEndcap");
   doBarrel_ = ps.getParameter<bool>("doBarrel");
  
   // The names of the produced cluster collections
   barrelClusterCollection_ = ps.getParameter<std::string>("barrelClusterCollection");
   endcapClusterCollection_ = ps.getParameter<std::string>("endcapClusterCollection");
  
   // Island algorithm parameters
   double barrelSeedThreshold = ps.getParameter<double>("IslandBarrelSeedThr");
   double endcapSeedThreshold = ps.getParameter<double>("IslandEndcapSeedThr");
  
   const std::vector<std::string> flagnames = ps.getParameter<std::vector<std::string> >("RecHitFlagToBeExcluded");
  
   const std::vector<int> v_chstatus = StringToEnumValue<EcalRecHit::Flags>(flagnames);
  
   // Parameters for the position calculation:
   edm::ParameterSet posCalcParameters = ps.getParameter<edm::ParameterSet>("posCalcParameters");
   posCalculator_ = PositionCalc(posCalcParameters);
  
   // Produces a collection of barrel and a collection of endcap clusters
   produces<reco::BasicClusterCollection>(endcapClusterCollection_);
   produces<reco::BasicClusterCollection>(barrelClusterCollection_);
  
   bool reassignSeedCrysToClusterItSeeds = false;
   if (ps.exists("reassignSeedCrysToClusterItSeeds"))
     reassignSeedCrysToClusterItSeeds = ps.getParameter<bool>("reassignSeedCrysToClusterItSeeds");
  
   island_p = new Multi5x5ClusterAlgo(
       barrelSeedThreshold, endcapSeedThreshold, v_chstatus, posCalculator_, reassignSeedCrysToClusterItSeeds);
  
   nEvt_ = 0;
 }

References barrelClusterCollection_, barrelHitToken_, caloGeometryToken_, doBarrel_, doEndcap_, endcapClusterCollection_, endcapHitToken_, edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), island_p, nEvt_, MonitorAlCaEcalPi0_cfi::posCalcParameters, posCalculator_, multi5x5BasicClusters_cfi::reassignSeedCrysToClusterItSeeds, and AlCaHLTBitMon_QueryRunRegistry::string.

◆ ~Multi5x5ClusterProducer()

Multi5x5ClusterProducer::~Multi5x5ClusterProducer ( )

override

Definition at line 78 of file Multi5x5ClusterProducer.cc.

78 { delete island_p; }

References island_p.

Member Function Documentation

◆ clusterizeECALPart()

void Multi5x5ClusterProducer::clusterizeECALPart	(	edm::Event &	evt,
		const edm::EventSetup &	es,
		const edm::EDGetTokenT< EcalRecHitCollection > &	token,
		const std::string &	clusterCollection,
		const reco::CaloID::Detectors	detector
	)

private

Definition at line 98 of file Multi5x5ClusterProducer.cc.

                                                                                      {
   // get the hit collection from the event:
   const EcalRecHitCollection* hitCollection_p = getCollection(evt, token);
  
   // get the geometry and topology from the event setup:
   edm::ESHandle<CaloGeometry> geoHandle = es.getHandle(caloGeometryToken_);
  
   const CaloSubdetectorGeometry* geometry_p;
   std::unique_ptr<CaloSubdetectorTopology> topology_p;
  
   if (detector == reco::CaloID::DET_ECAL_BARREL) {
     geometry_p = geoHandle->getSubdetectorGeometry(DetId::Ecal, EcalBarrel);
     topology_p = std::make_unique<EcalBarrelTopology>(*geoHandle);
   } else {
     geometry_p = geoHandle->getSubdetectorGeometry(DetId::Ecal, EcalEndcap);
     topology_p = std::make_unique<EcalEndcapTopology>(*geoHandle);
   }
  
   const CaloSubdetectorGeometry* geometryES_p;
   geometryES_p = geoHandle->getSubdetectorGeometry(DetId::Ecal, EcalPreshower);
  
   // Run the clusterization algorithm:
   reco::BasicClusterCollection clusters;
   clusters = island_p->makeClusters(hitCollection_p, geometry_p, topology_p.get(), geometryES_p, detector);
  
   // create a unique_ptr to a BasicClusterCollection, copy the barrel clusters into it and put in the Event:
   auto clusters_p = std::make_unique<reco::BasicClusterCollection>();
   clusters_p->assign(clusters.begin(), clusters.end());
   edm::OrphanHandle<reco::BasicClusterCollection> bccHandle;
   if (detector == reco::CaloID::DET_ECAL_BARREL)
     bccHandle = evt.put(std::move(clusters_p), barrelClusterCollection_);
   else
     bccHandle = evt.put(std::move(clusters_p), endcapClusterCollection_);
 }

References barrelClusterCollection_, caloGeometryToken_, bsc_activity_cfg::clusters, reco::CaloID::DET_ECAL_BARREL, hgcalTestNeighbor_cfi::detector, DetId::Ecal, EcalBarrel, EcalEndcap, EcalPreshower, endcapClusterCollection_, getCollection(), edm::EventSetup::getHandle(), CaloGeometry::getSubdetectorGeometry(), island_p, Multi5x5ClusterAlgo::makeClusters(), eostools::move(), edm::Event::put(), and unpackBuffers-CaloStage2::token.

Referenced by produce().

◆ counterExceeded()

bool Multi5x5ClusterProducer::counterExceeded ( ) const

inlineprivate

Definition at line 53 of file Multi5x5ClusterProducer.h.

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

References nEvt_, and nMaxPrintout_.

◆ getCollection()

const EcalRecHitCollection * Multi5x5ClusterProducer::getCollection	(	edm::Event &	evt,
		const edm::EDGetTokenT< EcalRecHitCollection > &	token
	)

private

Definition at line 91 of file Multi5x5ClusterProducer.cc.

                                                                       {
   edm::Handle<EcalRecHitCollection> rhcHandle;
   evt.getByToken(token, rhcHandle);
   return rhcHandle.product();
 }

References edm::Event::getByToken(), edm::Handle< T >::product(), and unpackBuffers-CaloStage2::token.

Referenced by clusterizeECALPart().

◆ outputValidationInfo()

void Multi5x5ClusterProducer::outputValidationInfo ( reco::CaloClusterPtrVector & clusterPtrVector )

private

◆ produce()

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

override

Definition at line 80 of file Multi5x5ClusterProducer.cc.

                                                                             {
   if (doEndcap_) {
     clusterizeECALPart(evt, es, endcapHitToken_, endcapClusterCollection_, reco::CaloID::DET_ECAL_ENDCAP);
   }
   if (doBarrel_) {
     clusterizeECALPart(evt, es, barrelHitToken_, barrelClusterCollection_, reco::CaloID::DET_ECAL_BARREL);
   }
  
   nEvt_++;
 }