Inheritance diagram for PFBadHcalPseudoClusterProducer:

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

	~PFBadHcalPseudoClusterProducer () override

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () 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 42 of file PFBadHcalPseudoClusterProducer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 64 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 ( )

override

Definition at line 74 of file PFBadHcalPseudoClusterProducer.cc.

75 {

76 }

Member Function Documentation

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

static

Definition at line 164 of file PFBadHcalPseudoClusterProducer.cc.

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

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

void PFBadHcalPseudoClusterProducer::init ( const EventSetup & c )

privatevirtual

Definition at line 78 of file PFBadHcalPseudoClusterProducer.cc.

References badAreasC_, badAreasRH_, reco::CaloCluster::badHcalMarker, Reference_intrackfit_cff::barrel, cacheId_geom_, cacheId_quality_, gather_cfg::cout, debug_, MillePedeFileConverter_cfg::e, edm::EventSetup::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, mps_update::status, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by produce().

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

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

overrideprivate

Definition at line 141 of file PFBadHcalPseudoClusterProducer.cc.

References badAreasC_, badAreasRH_, cacheId_geom_, cacheId_quality_, edm::eventsetup::EventSetupRecord::cacheIdentifier(), enabled_, edm::EventSetup::get(), mps_fire::i, init(), eostools::move(), gen::n, and edm::Event::put().

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