Inheritance diagram for EgammaHLTNxNClusterProducer:

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

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

	~EgammaHLTNxNClusterProducer () 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

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Private Member Functions
bool	checkStatusOfEcalRecHit (const EcalChannelStatus &channelStatus, const EcalRecHit &rh)

void	makeNxNClusters (edm::Event &evt, const edm::EventSetup &es, const EcalRecHitCollection *hits, const reco::CaloID::Detectors detector)

Private Attributes
const std::string	barrelClusterCollection_

const edm::EDGetTokenT< EcalRecHitCollection >	barrelHitProducer_

const int	clusEtaSize_

const int	clusPhiSize_

const double	clusSeedThr_

const double	clusSeedThrEndCap_

const int	debug_

const bool	doBarrel_

const bool	doEndcaps_

const std::string	endcapClusterCollection_

const edm::EDGetTokenT< EcalRecHitCollection >	endcapHitProducer_

const int	flagLevelRecHitsToUse_

const int	maxNumberofClusters_

const int	maxNumberofSeeds_

PositionCalc	posCalculator_

const int	statusLevelRecHitsToUse_

const bool	useDBStatus_

const bool	useRecoFlag_

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

Description: simple NxN ( 3x3 etc) clustering ,( for low energy photon reconstrution, currently used for pi0/eta HLT path)

Implementation: <Notes on="" implementation>="">

Definition at line 46 of file EgammaHLTNxNClusterProducer.cc.

Constructor & Destructor Documentation

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

Definition at line 88 of file EgammaHLTNxNClusterProducer.cc.

References barrelClusterCollection_, and endcapClusterCollection_.

     : doBarrel_(ps.getParameter<bool>("doBarrel")),
       doEndcaps_(ps.getParameter<bool>("doEndcaps")),
       barrelHitProducer_(consumes<EcalRecHitCollection>(ps.getParameter<edm::InputTag>("barrelHitProducer"))),
       endcapHitProducer_(consumes<EcalRecHitCollection>(ps.getParameter<edm::InputTag>("endcapHitProducer"))),
       clusEtaSize_(ps.getParameter<int>("clusEtaSize")),
       clusPhiSize_(ps.getParameter<int>("clusPhiSize")),
       barrelClusterCollection_(ps.getParameter<std::string>("barrelClusterCollection")),
       endcapClusterCollection_(ps.getParameter<std::string>("endcapClusterCollection")),
       clusSeedThr_(ps.getParameter<double>("clusSeedThr")),
       clusSeedThrEndCap_(ps.getParameter<double>("clusSeedThrEndCap")),
       useRecoFlag_(ps.getParameter<bool>("useRecoFlag")),
       flagLevelRecHitsToUse_(ps.getParameter<int>("flagLevelRecHitsToUse")),
       useDBStatus_(ps.getParameter<bool>("useDBStatus")),
       statusLevelRecHitsToUse_(ps.getParameter<int>("statusLevelRecHitsToUse")),
       maxNumberofSeeds_(ps.getParameter<int>("maxNumberofSeeds")),
       maxNumberofClusters_(ps.getParameter<int>("maxNumberofClusters")),
       debug_(ps.getParameter<int>("debugLevel")),
       posCalculator_(PositionCalc(ps.getParameter<edm::ParameterSet>("posCalcParameters"))) {
   produces<reco::BasicClusterCollection>(barrelClusterCollection_);
   produces<reco::BasicClusterCollection>(endcapClusterCollection_);
 }

EgammaHLTNxNClusterProducer::~EgammaHLTNxNClusterProducer ( )

override

Definition at line 111 of file EgammaHLTNxNClusterProducer.cc.

111 {}

Member Function Documentation

bool EgammaHLTNxNClusterProducer::checkStatusOfEcalRecHit	(	const EcalChannelStatus &	channelStatus,
		const EcalRecHit &	rh
	)

private

from recoFlag()

good

good || PoorCalib

good || PoorCalib || LeadingEdgeRecovered || kNeighboursRecovered,

Definition at line 176 of file EgammaHLTNxNClusterProducer.cc.

References RemoveAddSevLevel::flag, flagLevelRecHitsToUse_, EcalRecHit::id(), createfilelist::int, DetId::rawId(), EcalRecHit::recoFlag(), mps_update::status, statusLevelRecHitsToUse_, useDBStatus_, and useRecoFlag_.

Referenced by makeNxNClusters().

                                                                                 {
   if (useRecoFlag_) {  
     int flag = rh.recoFlag();
     if (flagLevelRecHitsToUse_ == 0) {  
       if (flag != 0)
         return false;
     } else if (flagLevelRecHitsToUse_ == 1) {  
       if (flag != 0 && flag != 4)
         return false;
     } else if (flagLevelRecHitsToUse_ == 2) {  
       if (flag != 0 && flag != 4 && flag != 6 && flag != 7)
         return false;
     }
   }
   if (useDBStatus_) {  
     int status = int(channelStatus[rh.id().rawId()].getStatusCode());
     if (status > statusLevelRecHitsToUse_)
       return false;
   }
 
   return true;
 }

void EgammaHLTNxNClusterProducer::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 113 of file EgammaHLTNxNClusterProducer.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                              {
   edm::ParameterSetDescription desc;
   desc.add<bool>(("doBarrel"), true);
   desc.add<bool>(("doEndcaps"), true);
   desc.add<edm::InputTag>(("barrelHitProducer"), edm::InputTag("hltEcalRegionalPi0EtaRecHit", "EcalRecHitsEB"));
   desc.add<edm::InputTag>(("endcapHitProducer"), edm::InputTag("hltEcalRegionalPi0EtaRecHit", "EcalRecHitsEE"));
   desc.add<int>(("clusEtaSize"), 3);
   desc.add<int>(("clusPhiSize"), 3);
   desc.add<std::string>(("barrelClusterCollection"), "Simple3x3ClustersBarrel");
   desc.add<std::string>(("endcapClusterCollection"), "Simple3x3ClustersEndcap");
   desc.add<double>(("clusSeedThr"), 0.5);
   desc.add<double>(("clusSeedThrEndCap"), 1.0);
   desc.add<bool>(("useRecoFlag"), false);
   desc.add<int>(("flagLevelRecHitsToUse"), 1);
   desc.add<bool>(("useDBStatus"), true);
   desc.add<int>(("statusLevelRecHitsToUse"), 1);
 
   edm::ParameterSetDescription posCalcPSET;
   posCalcPSET.add<double>("T0_barl", 7.4);
   posCalcPSET.add<double>("T0_endc", 3.1);
   posCalcPSET.add<double>("T0_endcPresh", 1.2);
   posCalcPSET.add<double>("W0", 4.2);
   posCalcPSET.add<double>("X0", 0.89);
   posCalcPSET.add<bool>("LogWeighted", true);
   desc.add<edm::ParameterSetDescription>("posCalcParameters", posCalcPSET);
 
   desc.add<int>(("maxNumberofSeeds"), 1000);
   desc.add<int>(("maxNumberofClusters"), 200);
   desc.add<int>(("debugLevel"), 0);
   descriptions.add(("hltEgammaHLTNxNClusterProducer"), desc);
 }

void EgammaHLTNxNClusterProducer::makeNxNClusters	(	edm::Event &	evt,
		const edm::EventSetup &	es,
		const EcalRecHitCollection *	hits,
		const reco::CaloID::Detectors	detector
	)

private

get status from DB

if too much clusters made, then return 0 also

Definition at line 200 of file EgammaHLTNxNClusterProducer.cc.

References barrelClusterCollection_, edm::SortedCollection< T, SORT >::begin(), PositionCalc::Calculate_Location(), checkStatusOfEcalRecHit(), clusEtaSize_, clusPhiSize_, clusSeedThr_, clusSeedThrEndCap_, bsc_activity_cfg::clusters, debug_, DEFINE_FWK_MODULE, reco::CaloID::DET_ECAL_BARREL, DetId::Ecal, EcalBarrel, EcalEndcap, EcalPreshower, edm::SortedCollection< T, SORT >::end(), endcapClusterCollection_, HCALHighEnergyHPDFilter_cfi::energy, edm::EventID::event(), spr::find(), edm::SortedCollection< T, SORT >::find(), edm::EventSetup::get(), CaloGeometry::getSubdetectorGeometry(), edm::EventBase::id(), reco::CaloCluster::island, LogDebug, maxNumberofClusters_, maxNumberofSeeds_, eostools::move(), posCalculator_, edm::Event::put(), edm::EventID::run(), useDBStatus_, x, and y.

Referenced by produce().

                                                                                       {
   edm::ESHandle<EcalChannelStatus> csHandle;
   if (useDBStatus_)
     es.get<EcalChannelStatusRcd>().get(csHandle);
   const EcalChannelStatus &channelStatus = *csHandle;
 
   std::vector<EcalRecHit> seeds;
 
   double clusterSeedThreshold;
   if (detector == reco::CaloID::DET_ECAL_BARREL) {
     clusterSeedThreshold = clusSeedThr_;
   } else {
     clusterSeedThreshold = clusSeedThrEndCap_;
   }
 
   for (EcalRecHitCollection::const_iterator itt = hits->begin(); itt != hits->end(); itt++) {
     double energy = itt->energy();
     if (!checkStatusOfEcalRecHit(channelStatus, *itt))
       continue;
     if (energy > clusterSeedThreshold)
       seeds.push_back(*itt);
 
     if (int(seeds.size()) > maxNumberofSeeds_) {  //too many seeds, like beam splash events
       seeds.clear();                              
       break;
     }
   }
 
   // get the geometry and topology from the event setup:
   edm::ESHandle<CaloGeometry> geoHandle;
   es.get<CaloGeometryRecord>().get(geoHandle);
 
   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);
 
   std::vector<reco::BasicCluster> clusters;
   std::vector<DetId> usedXtals;
 
   // sort seed according to Energy
   sort(seeds.begin(), seeds.end(), [](auto const &x, auto const &y) { return (x.energy() > y.energy()); });
 
   for (std::vector<EcalRecHit>::iterator itseed = seeds.begin(); itseed != seeds.end(); itseed++) {
     DetId seed_id = itseed->id();
 
     std::vector<DetId>::iterator itdet = find(usedXtals.begin(), usedXtals.end(), seed_id);
     if (itdet != usedXtals.end())
       continue;
 
     std::vector<DetId> clus_v = topology_p->getWindow(seed_id, clusEtaSize_, clusPhiSize_);
     std::vector<std::pair<DetId, float> > clus_used;
 
     float clus_energy = 0;
 
     for (std::vector<DetId>::iterator det = clus_v.begin(); det != clus_v.end(); det++) {
       DetId detid = *det;
 
       //not yet used
       std::vector<DetId>::iterator itdet = find(usedXtals.begin(), usedXtals.end(), detid);
       if (itdet != usedXtals.end())
         continue;
       //inside the collection
       EcalRecHitCollection::const_iterator hit = hits->find(detid);
       if (hit == hits->end())
         continue;
 
       if (!checkStatusOfEcalRecHit(channelStatus, *hit))
         continue;
 
       usedXtals.push_back(detid);
       clus_used.push_back(std::pair<DetId, float>(detid, 1.));
       clus_energy += hit->energy();
 
     }  
 
     if (clus_energy <= 0)
       continue;
 
     math::XYZPoint clus_pos = posCalculator_.Calculate_Location(clus_used, hits, geometry_p, geometryES_p);
 
     if (debug_ >= 2)
       LogDebug("") << "nxn_cluster in run " << evt.id().run() << " event " << evt.id().event()
                    << " energy: " << clus_energy << " eta: " << clus_pos.Eta() << " phi: " << clus_pos.Phi()
                    << " nRecHits: " << clus_used.size() << std::endl;
 
     clusters.push_back(reco::BasicCluster(
         clus_energy, clus_pos, reco::CaloID(detector), clus_used, reco::CaloCluster::island, seed_id));
     if (int(clusters.size()) > maxNumberofClusters_) {  
       clusters.clear();
       break;
     }
   }
 
   //Create empty output collections
   auto clusters_p = std::make_unique<reco::BasicClusterCollection>();
   clusters_p->assign(clusters.begin(), clusters.end());
   if (detector == reco::CaloID::DET_ECAL_BARREL) {
     if (debug_ >= 1)
       LogDebug("") << "nxnclusterProducer: " << clusters_p->size() << " made in barrel" << std::endl;
     evt.put(std::move(clusters_p), barrelClusterCollection_);
   } else {
     if (debug_ >= 1)
       LogDebug("") << "nxnclusterProducer: " << clusters_p->size() << " made in endcap" << std::endl;
     evt.put(std::move(clusters_p), endcapClusterCollection_);
   }
 }

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

override

Definition at line 145 of file EgammaHLTNxNClusterProducer.cc.

References barrelHitProducer_, debug_, reco::CaloID::DET_ECAL_BARREL, reco::CaloID::DET_ECAL_ENDCAP, doBarrel_, doEndcaps_, endcapHitProducer_, edm::EventID::event(), edm::Event::getByToken(), edm::EventBase::id(), edm::HandleBase::isValid(), LogDebug, makeNxNClusters(), edm::Handle< T >::product(), edm::EventID::run(), and edm::SortedCollection< T, SORT >::size().

                                                                                 {
   if (doBarrel_) {
     edm::Handle<EcalRecHitCollection> barrelRecHitsHandle;
     evt.getByToken(barrelHitProducer_, barrelRecHitsHandle);
     if (!barrelRecHitsHandle.isValid()) {
       LogDebug("") << "EgammaHLTNxNClusterProducer Error! can't get product eb hit!" << std::endl;
     }
 
     const EcalRecHitCollection *hits_eb = barrelRecHitsHandle.product();
     if (debug_ >= 2)
       LogDebug("") << "EgammaHLTNxNClusterProducer nEBrechits: " << evt.id().run() << " event " << evt.id().event()
                    << " " << hits_eb->size() << std::endl;
 
     makeNxNClusters(evt, es, hits_eb, reco::CaloID::DET_ECAL_BARREL);
   }
 
   if (doEndcaps_) {
     edm::Handle<EcalRecHitCollection> endcapRecHitsHandle;
     evt.getByToken(endcapHitProducer_, endcapRecHitsHandle);
     if (!endcapRecHitsHandle.isValid()) {
       LogDebug("") << "EgammaHLTNxNClusterProducer Error! can't get product ee hit!" << std::endl;
     }
 
     const EcalRecHitCollection *hits_ee = endcapRecHitsHandle.product();
     if (debug_ >= 2)
       LogDebug("") << "EgammaHLTNxNClusterProducer nEErechits: " << evt.id().run() << " event " << evt.id().event()
                    << " " << hits_ee->size() << std::endl;
     makeNxNClusters(evt, es, hits_ee, reco::CaloID::DET_ECAL_ENDCAP);
   }
 }