Inheritance diagram for PFBadHcalPseudoClusterProducer:

Public Member Functions
	PFBadHcalPseudoClusterProducer (const edm::ParameterSet &)

	~PFBadHcalPseudoClusterProducer () 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
virtual void	init (const EventSetup &c)

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

Private Attributes
std::vector< reco::PFCluster >	badAreasC_

std::vector< reco::PFRecHit >	badAreasRH_

unsigned long long	cacheId_geom_

unsigned long long	cacheId_quality_

bool	debug_

bool	enabled_

edm::ESHandle< CaloGeometry >	hGeom_

edm::ESHandle< HcalChannelQuality >	hQuality_

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 41 of file PFBadHcalPseudoClusterProducer.cc.

Constructor & Destructor Documentation

◆ PFBadHcalPseudoClusterProducer()

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

explicit

Definition at line 62 of file PFBadHcalPseudoClusterProducer.cc.

     : enabled_(ps.getParameter<bool>("enable")),
       debug_(ps.getUntrackedParameter<bool>("debug", false)),
       cacheId_quality_(0),
       cacheId_geom_(0) {
   produces<std::vector<reco::PFCluster>>();
   produces<std::vector<reco::PFRecHit>>("hits");
 }

◆ ~PFBadHcalPseudoClusterProducer()

PFBadHcalPseudoClusterProducer::~PFBadHcalPseudoClusterProducer ( )

override

Definition at line 71 of file PFBadHcalPseudoClusterProducer.cc.

71 {}

Member Function Documentation

◆ fillDescriptions()

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

static

Definition at line 165 of file PFBadHcalPseudoClusterProducer.cc.

                                                                                                 {
   edm::ParameterSetDescription desc;
   desc.add<bool>("enable", false)
       ->setComment("activate the module (if false, it doesn't check the DB and produces an empty collection)");
   desc.addUntracked<bool>("debug", false);
   descriptions.add("particleFlowBadHcalPseudoCluster", desc);
 }

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), and edm::ParameterSetDescription::addUntracked().

◆ init()

void PFBadHcalPseudoClusterProducer::init ( const EventSetup & c )

privatevirtual

Definition at line 73 of file PFBadHcalPseudoClusterProducer.cc.

                                                                   {
   badAreasC_.clear();
   badAreasRH_.clear();
  
   edm::ESHandle<HcalChannelQuality> hQuality_;
   iSetup.get<HcalChannelQualityRcd>().get("withTopo", hQuality_);
   const HcalChannelQuality& chanquality = *hQuality_;
  
   edm::ESHandle<CaloGeometry> hGeom_;
   iSetup.get<CaloGeometryRecord>().get(hGeom_);
   const CaloGeometry& caloGeom = *hGeom_;
   const CaloSubdetectorGeometry* hbGeom = caloGeom.getSubdetectorGeometry(DetId::Hcal, HcalBarrel);
   const CaloSubdetectorGeometry* heGeom = caloGeom.getSubdetectorGeometry(DetId::Hcal, HcalEndcap);
  
   int statusMask = ((1 << HcalChannelStatus::HcalCellOff) | (1 << HcalChannelStatus::HcalCellMask) |
                     (1 << HcalChannelStatus::HcalCellDead));
   // histogram the number of bad depths at each ieta, iphi
   std::map<std::pair<int, int>, int> good, bads;
   std::map<std::pair<int, int>, std::pair<int, HcalSubdetector>> minDepths;
   for (const DetId& i : chanquality.getAllChannels()) {
     if (i.det() != DetId::Hcal)
       continue;  // not an hcal cell
     HcalDetId id = HcalDetId(i);
     if (id.subdet() != HcalBarrel && id.subdet() != HcalEndcap)
       continue;  // we don't deal with HO and HF here
     int status = chanquality.getValues(id)->getValue();
     auto tower = std::make_pair(id.ieta(), id.iphi());
     if (status & statusMask) {
       bads[tower]++;
       if (debug_)
         std::cout << "Channel " << i() << " (subdet " << id.subdet() << ", zside " << id.zside() << ", ieta "
                   << id.ieta() << ", iphi " << id.iphi() << " depth " << id.depth() << " has status " << status
                   << "  masked " << (status & statusMask) << std::endl;
     } else {
       good[tower]++;
     }
     auto& minD = minDepths[tower];
     if (minD.second == HcalEmpty || minD.first > id.depth()) {
       minD.first = id.depth();
       minD.second = id.subdet();
     }
   }
  
   const float dummyEnergy = 1e-5;  // non-zero, but small (even if it shouldn't ever be used)
   for (const auto& rec : bads) {
     int ieta = rec.first.first, iphi = rec.first.second, nbad = rec.second, ngood = good[rec.first];
     auto minDepth = minDepths[rec.first];
     bool barrel = minDepth.second == HcalBarrel;
     HcalDetId id(minDepth.second, ieta, iphi, minDepth.first);
     bool isBad = (nbad > 0 && nbad >= ngood);
     if (debug_)
       std::cout << "At ieta " << id.ieta() << ", iphi " << id.iphi() << " I have " << nbad << " bad depths, " << ngood
                 << " good depths. First depth is in " << (barrel ? "HB" : "HE") << " depth " << minDepth.first << "; "
                 << (isBad ? " MARK BAD" : " ignore") << std::endl;
     if (!isBad)
       continue;
     PFLayer::Layer layer = (barrel ? PFLayer::HCAL_BARREL1 : PFLayer::HCAL_ENDCAP);
     // make a PF RecHit
     std::shared_ptr<const CaloCellGeometry> thisCell = (barrel ? hbGeom : heGeom)->getGeometry(id);
     const GlobalPoint& pos = thisCell->getPosition();
     badAreasRH_.emplace_back(thisCell, id(), layer, dummyEnergy);
     // make a PF cluster (but can't add the rechit, as for that I need an edm::Ref)
     badAreasC_.emplace_back(layer, dummyEnergy, pos.x(), pos.y(), pos.z());
     badAreasC_.back().setFlags(reco::CaloCluster::badHcalMarker);
   }
  
   cacheId_quality_ = iSetup.get<HcalChannelQualityRcd>().cacheIdentifier();
   cacheId_geom_ = iSetup.get<CaloGeometryRecord>().cacheIdentifier();
 }

References badAreasC_, badAreasRH_, reco::CaloCluster::badHcalMarker, Reference_intrackfit_cff::barrel, cacheId_geom_, cacheId_quality_, gather_cfg::cout, debug_, MillePedeFileConverter_cfg::e, edm::EventSetup::get(), get, HcalCondObjectContainer< Item >::getAllChannels(), ecaldqm::getGeometry(), CaloGeometry::getSubdetectorGeometry(), HcalChannelStatus::getValue(), HcalCondObjectContainer< Item >::getValues(), DetId::Hcal, PFLayer::HCAL_BARREL1, PFLayer::HCAL_ENDCAP, HcalBarrel, HcalChannelStatus::HcalCellDead, HcalChannelStatus::HcalCellMask, HcalChannelStatus::HcalCellOff, HcalEmpty, HcalEndcap, hGeom_, hQuality_, mps_fire::i, triggerObjects_cff::id, LEDCalibrationChannels::ieta, LEDCalibrationChannels::iphi, mps_update::status, and hgcalTowerProducer_cfi::tower.

Referenced by produce().

◆ produce()

void PFBadHcalPseudoClusterProducer::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprivate

Definition at line 144 of file PFBadHcalPseudoClusterProducer.cc.

                                                                                           {
   if (enabled_) {
     if (cacheId_quality_ != iSetup.get<HcalChannelQualityRcd>().cacheIdentifier() ||
         cacheId_geom_ != iSetup.get<CaloGeometryRecord>().cacheIdentifier()) {
       init(iSetup);
     }
   }
  
   auto outRH = std::make_unique<std::vector<reco::PFRecHit>>(badAreasRH_);
   auto rhHandle = iEvent.put(std::move(outRH), "hits");
  
   auto outC = std::make_unique<std::vector<reco::PFCluster>>(badAreasC_);
  
   // now we go set references
   for (unsigned int i = 0, n = rhHandle->size(); i < n; ++i) {
     (*outC)[i].addRecHitFraction(reco::PFRecHitFraction(reco::PFRecHitRef(rhHandle, i), 1.0));
   }
  
   iEvent.put(std::move(outC));
 }