#include <HBHEIsolatedNoiseReflagger.h>

Inheritance diagram for HBHEIsolatedNoiseReflagger:

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

	~HBHEIsolatedNoiseReflagger () override

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

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

Private Member Functions
void	DumpHBHEHitMap (std::vector< HBHEHitMap > &i) const

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

Private Attributes
edm::ESGetToken< CaloTowerConstituentsMap, CaloGeometryRecord >	ctcmToken_

bool	debug_

edm::ESGetToken< EcalChannelStatus, EcalChannelStatusRcd >	ecalChStatusToken_

edm::ESGetToken< EcalSeverityLevelAlgo, EcalSeverityLevelAlgoRcd >	ecalSevToken_

edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	geoToken_

edm::ESGetToken< HcalChannelQuality, HcalChannelQualityRcd >	hcalChStatusToken_

edm::ESGetToken< HcalSeverityLevelComputer, HcalSeverityLevelComputerRcd >	hcalSevToken_

const HcalFrontEndMap *	hfemap

edm::ESGetToken< HcalFrontEndMap, HcalFrontEndMapRcd >	hfemapToken_

double	LooseDiHitEne_

double	LooseEcalIsol_

double	LooseHcalIsol_

double	LooseHPDEne1_

double	LooseHPDEne2_

int	LooseHPDHits1_

int	LooseHPDHits2_

double	LooseMonoHitEne_

double	LooseRBXEne1_

double	LooseRBXEne2_

int	LooseRBXHits1_

int	LooseRBXHits2_

double	LooseTrackIsol_

ObjectValidator	objvalidator_

double	RBXEneThreshold_

double	TightDiHitEne_

double	TightEcalIsol_

double	TightHcalIsol_

double	TightHPDEne1_

double	TightHPDEne2_

int	TightHPDHits1_

int	TightHPDHits2_

double	TightMonoHitEne_

double	TightRBXEne1_

double	TightRBXEne2_

int	TightRBXHits1_

int	TightRBXHits2_

double	TightTrackIsol_

edm::EDGetTokenT< EcalRecHitCollection >	tok_EB_

edm::EDGetTokenT< EcalRecHitCollection >	tok_EE_

edm::EDGetTokenT< HBHERecHitCollection >	tok_hbhe_

edm::EDGetTokenT< std::vector< reco::TrackExtrapolation > >	tok_trackExt_

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
using	CacheTypes = CacheContexts< T... >

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T... >

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Detailed Description

Definition at line 31 of file HBHEIsolatedNoiseReflagger.h.

Constructor & Destructor Documentation

◆ HBHEIsolatedNoiseReflagger()

HBHEIsolatedNoiseReflagger::HBHEIsolatedNoiseReflagger ( const edm::ParameterSet & iConfig )

explicit

Definition at line 23 of file HBHEIsolatedNoiseReflagger.cc.

References ctcmToken_, ecalChStatusToken_, ecalSevToken_, deDxTools::esConsumes(), geoToken_, edm::ParameterSet::getParameter(), hcalChStatusToken_, hcalSevToken_, hfemapToken_, tok_EB_, tok_EE_, tok_hbhe_, and tok_trackExt_.

     :
 
       LooseHcalIsol_(iConfig.getParameter<double>("LooseHcalIsol")),
       LooseEcalIsol_(iConfig.getParameter<double>("LooseEcalIsol")),
       LooseTrackIsol_(iConfig.getParameter<double>("LooseTrackIsol")),
       TightHcalIsol_(iConfig.getParameter<double>("TightHcalIsol")),
       TightEcalIsol_(iConfig.getParameter<double>("TightEcalIsol")),
       TightTrackIsol_(iConfig.getParameter<double>("TightTrackIsol")),
 
       LooseRBXEne1_(iConfig.getParameter<double>("LooseRBXEne1")),
       LooseRBXEne2_(iConfig.getParameter<double>("LooseRBXEne2")),
       LooseRBXHits1_(iConfig.getParameter<int>("LooseRBXHits1")),
       LooseRBXHits2_(iConfig.getParameter<int>("LooseRBXHits2")),
       TightRBXEne1_(iConfig.getParameter<double>("TightRBXEne1")),
       TightRBXEne2_(iConfig.getParameter<double>("TightRBXEne2")),
       TightRBXHits1_(iConfig.getParameter<int>("TightRBXHits1")),
       TightRBXHits2_(iConfig.getParameter<int>("TightRBXHits2")),
 
       LooseHPDEne1_(iConfig.getParameter<double>("LooseHPDEne1")),
       LooseHPDEne2_(iConfig.getParameter<double>("LooseHPDEne2")),
       LooseHPDHits1_(iConfig.getParameter<int>("LooseHPDHits1")),
       LooseHPDHits2_(iConfig.getParameter<int>("LooseHPDHits2")),
       TightHPDEne1_(iConfig.getParameter<double>("TightHPDEne1")),
       TightHPDEne2_(iConfig.getParameter<double>("TightHPDEne2")),
       TightHPDHits1_(iConfig.getParameter<int>("TightHPDHits1")),
       TightHPDHits2_(iConfig.getParameter<int>("TightHPDHits2")),
 
       LooseDiHitEne_(iConfig.getParameter<double>("LooseDiHitEne")),
       TightDiHitEne_(iConfig.getParameter<double>("TightDiHitEne")),
       LooseMonoHitEne_(iConfig.getParameter<double>("LooseMonoHitEne")),
       TightMonoHitEne_(iConfig.getParameter<double>("TightMonoHitEne")),
 
       RBXEneThreshold_(iConfig.getParameter<double>("RBXEneThreshold")),
 
       debug_(iConfig.getUntrackedParameter<bool>("debug", true)),
       objvalidator_(iConfig) {
   tok_hbhe_ = consumes<HBHERecHitCollection>(iConfig.getParameter<edm::InputTag>("hbheInput"));
   tok_EB_ = consumes<EcalRecHitCollection>(iConfig.getParameter<edm::InputTag>("ebInput"));
   tok_EE_ = consumes<EcalRecHitCollection>(iConfig.getParameter<edm::InputTag>("eeInput"));
   tok_trackExt_ =
       consumes<std::vector<reco::TrackExtrapolation> >(iConfig.getParameter<edm::InputTag>("trackExtrapolationInput"));
 
   // ES tokens
   ecalChStatusToken_ = esConsumes<EcalChannelStatus, EcalChannelStatusRcd>();
   hcalChStatusToken_ = esConsumes<HcalChannelQuality, HcalChannelQualityRcd>(edm::ESInputTag("", "withTopo"));
   hcalSevToken_ = esConsumes<HcalSeverityLevelComputer, HcalSeverityLevelComputerRcd>();
   ecalSevToken_ = esConsumes<EcalSeverityLevelAlgo, EcalSeverityLevelAlgoRcd>();
   ctcmToken_ = esConsumes<CaloTowerConstituentsMap, CaloGeometryRecord>();
   hfemapToken_ = esConsumes<HcalFrontEndMap, HcalFrontEndMapRcd>();
   geoToken_ = esConsumes();
 
   produces<HBHERecHitCollection>();
 }

◆ ~HBHEIsolatedNoiseReflagger()

HBHEIsolatedNoiseReflagger::~HBHEIsolatedNoiseReflagger ( )

override

Definition at line 78 of file HBHEIsolatedNoiseReflagger.cc.

78 {}

Member Function Documentation

◆ DumpHBHEHitMap()

void HBHEIsolatedNoiseReflagger::DumpHBHEHitMap ( std::vector< HBHEHitMap > & i ) const

private

Definition at line 255 of file HBHEIsolatedNoiseReflagger.cc.

References hfemap, mps_fire::i, hit::id, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, HcalFrontEndMap::lookupRBX(), and HcalFrontEndMap::lookupRMIndex().

Referenced by produce().

                                                                               {
   for (std::vector<HBHEHitMap>::const_iterator it = i.begin(); it != i.end(); ++it) {
     edm::LogInfo("HBHEIsolatedNoiseReflagger")
         << "hit energy=" << it->hitEnergy() << "; # hits=" << it->nHits()
         << "; hcal energy same=" << it->hcalEnergySameTowers() << "; ecal energy same=" << it->ecalEnergySameTowers()
         << "; track energy same=" << it->trackEnergySameTowers()
         << "; neighbor hcal energy=" << it->hcalEnergyNeighborTowers() << std::endl;
     edm::LogInfo("HBHEIsolatedNoiseReflagger") << "hits:" << std::endl;
     for (HBHEHitMap::hitmap_const_iterator it2 = it->beginHits(); it2 != it->endHits(); ++it2) {
       const HBHERecHit* hit = it2->first;
       edm::LogInfo("HBHEIsolatedNoiseReflagger")
           << "RBX #=" << hfemap->lookupRBX(hit->id()) << "; HPD #=" << hfemap->lookupRMIndex(hit->id()) << "; "
           << (*hit) << std::endl;
     }
   }
   return;
 }

◆ produce()

void HBHEIsolatedNoiseReflagger::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	evSetup
	)

overrideprivate

Definition at line 80 of file HBHEIsolatedNoiseReflagger.cc.

References edm::SortedCollection< T, SORT >::begin(), ctcmToken_, debug_, DumpHBHEHitMap(), ecalChStatusToken_, ecalSevToken_, edm::SortedCollection< T, SORT >::end(), geoToken_, edm::EventSetup::getData(), HcalCaloFlagLabels::HBHEIsolatedNoise, hcalChStatusToken_, hcalSevToken_, hfemap, hfemapToken_, mps_fire::i, iEvent, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, LooseDiHitEne_, LooseEcalIsol_, LooseHcalIsol_, LooseHPDEne1_, LooseHPDEne2_, LooseHPDHits1_, LooseHPDHits2_, LooseMonoHitEne_, LooseRBXEne1_, LooseRBXEne2_, LooseRBXHits1_, LooseRBXHits2_, LooseTrackIsol_, eostools::move(), TrackingDataMCValidation_Standalone_cff::nhits, objvalidator_, RBXEneThreshold_, ObjectValidator::setEBRecHitCollection(), ObjectValidator::setEcalChannelStatus(), ObjectValidator::setEcalSeverityLevelAlgo(), ObjectValidator::setEERecHitCollection(), CaloRecHit::setFlagField(), ObjectValidator::setHcalChannelQuality(), ObjectValidator::setHcalSeverityLevelComputer(), TightDiHitEne_, TightEcalIsol_, TightHcalIsol_, TightHPDEne1_, TightHPDEne2_, TightHPDHits1_, TightHPDHits2_, TightMonoHitEne_, TightRBXEne1_, TightRBXEne2_, TightRBXHits1_, TightRBXHits2_, TightTrackIsol_, tok_EB_, tok_EE_, tok_hbhe_, and tok_trackExt_.

                                                                                        {
   // get the ECAL channel status map
   const EcalChannelStatus* dbEcalChStatus = &evSetup.getData(ecalChStatusToken_);
 
   // get the HCAL channel status map
   const HcalChannelQuality* dbHcalChStatus = &evSetup.getData(hcalChStatusToken_);
 
   // get the severity level computers
   const HcalSeverityLevelComputer* hcalSevLvlComputer = &evSetup.getData(hcalSevToken_);
   const EcalSeverityLevelAlgo* ecalSevLvlAlgo = &evSetup.getData(ecalSevToken_);
 
   // get the calotower mappings
   const CaloTowerConstituentsMap& ctcm = evSetup.getData(ctcmToken_);
 
   // get hcal frontend map
   const HcalFrontEndMap* hfemap = &evSetup.getData(hfemapToken_);
 
   // get the HB/HE hits
   edm::Handle<HBHERecHitCollection> hbhehits_h;
   iEvent.getByToken(tok_hbhe_, hbhehits_h);
 
   // get the ECAL hits
   edm::Handle<EcalRecHitCollection> ebhits_h;
   iEvent.getByToken(tok_EB_, ebhits_h);
   edm::Handle<EcalRecHitCollection> eehits_h;
   iEvent.getByToken(tok_EE_, eehits_h);
 
   // get the tracks
   edm::Handle<std::vector<reco::TrackExtrapolation> > trackextraps_h;
   iEvent.getByToken(tok_trackExt_, trackextraps_h);
 
   // set the status maps and severity level computers for the hit validator
   objvalidator_.setHcalChannelQuality(dbHcalChStatus);
   objvalidator_.setEcalChannelStatus(dbEcalChStatus);
   objvalidator_.setHcalSeverityLevelComputer(hcalSevLvlComputer);
   objvalidator_.setEcalSeverityLevelAlgo(ecalSevLvlAlgo);
   objvalidator_.setEBRecHitCollection(&(*ebhits_h));
   objvalidator_.setEERecHitCollection(&(*eehits_h));
 
   // organizer the hits
   PhysicsTowerOrganizer pto(
       hbhehits_h, ebhits_h, eehits_h, trackextraps_h, objvalidator_, ctcm, evSetup.getData(geoToken_));
   HBHEHitMapOrganizer organizer(hbhehits_h, objvalidator_, pto, hfemap);
 
   // get the rbxs, hpds, dihits, and monohits
   std::vector<HBHEHitMap> rbxs;
   std::vector<HBHEHitMap> hpds;
   std::vector<HBHEHitMap> dihits;
   std::vector<HBHEHitMap> monohits;
   organizer.getRBXs(rbxs, LooseRBXEne1_ < TightRBXEne1_ ? LooseRBXEne1_ : TightRBXEne1_);
   organizer.getHPDs(hpds, LooseHPDEne1_ < TightHPDEne1_ ? LooseHPDEne1_ : TightHPDEne1_);
   organizer.getDiHits(dihits, LooseDiHitEne_ < TightDiHitEne_ ? LooseDiHitEne_ : TightDiHitEne_);
   organizer.getMonoHits(monohits, LooseMonoHitEne_ < TightMonoHitEne_ ? LooseMonoHitEne_ : TightMonoHitEne_);
 
   if (debug_ && (!rbxs.empty() || !hpds.empty() || !dihits.empty() || !monohits.empty())) {
     edm::LogInfo("HBHEIsolatedNoiseReflagger") << "RBXs:" << std::endl;
     DumpHBHEHitMap(rbxs);
     edm::LogInfo("HBHEIsolatedNoiseReflagger") << "\nHPDs:" << std::endl;
     DumpHBHEHitMap(hpds);
     edm::LogInfo("HBHEIsolatedNoiseReflagger") << "\nDiHits:" << std::endl;
     DumpHBHEHitMap(dihits);
     edm::LogInfo("HBHEIsolatedNoiseReflagger") << "\nMonoHits:" << std::endl;
     DumpHBHEHitMap(monohits);
   }
 
   //  bool result=true;
 
   // determine which hits are noisy
   std::set<const HBHERecHit*> noisehits;
   for (int i = 0; i < static_cast<int>(rbxs.size()); i++) {
     int nhits = rbxs[i].nHits();
     double ene = rbxs[i].hitEnergy();
     double trkfide = rbxs[i].hitEnergyTrackFiducial();
     double isolhcale = rbxs[i].hcalEnergySameTowers() + rbxs[i].hcalEnergyNeighborTowers();
     double isolecale = rbxs[i].ecalEnergySameTowers();
     double isoltrke = rbxs[i].trackEnergySameTowers() + rbxs[i].trackEnergyNeighborTowers();
     //
     // RBX mistag reduction
     bool isLooseIso = false;
     bool isTightIso = false;
     if (ene > RBXEneThreshold_ && ene > 0) {  // New absolute iso-cut for high energy RBX clusters
       if (isolhcale < LooseHcalIsol_ * RBXEneThreshold_ && isolecale < LooseEcalIsol_ * RBXEneThreshold_ &&
           isoltrke < LooseTrackIsol_ * RBXEneThreshold_)
         isLooseIso = true;
       if (isolhcale < TightHcalIsol_ * RBXEneThreshold_ && isolecale < TightEcalIsol_ * RBXEneThreshold_ &&
           isoltrke < TightTrackIsol_ * RBXEneThreshold_)
         isTightIso = true;
     }
     if (ene <= RBXEneThreshold_ && ene > 0) {  // Old relative iso-cut for low energy RBX clusters
       if (isolhcale / ene < LooseHcalIsol_ && isolecale / ene < LooseEcalIsol_ && isoltrke / ene < LooseTrackIsol_)
         isLooseIso = true;
       if (isolhcale / ene < TightHcalIsol_ && isolecale / ene < TightEcalIsol_ && isoltrke / ene < TightTrackIsol_)
         isTightIso = true;
     }
     //
     if ((isLooseIso && ((trkfide > LooseRBXEne1_ && nhits >= LooseRBXHits1_) ||
                         (trkfide > LooseRBXEne2_ && nhits >= LooseRBXHits2_))) ||
         (isTightIso && ((trkfide > TightRBXEne1_ && nhits >= TightRBXHits1_) ||
                         (trkfide > TightRBXEne2_ && nhits >= TightRBXHits2_)))) {
       for (HBHEHitMap::hitmap_const_iterator it = rbxs[i].beginHits(); it != rbxs[i].endHits(); ++it)
         noisehits.insert(it->first);
       //      result=false;
     }
   }
 
   for (int i = 0; i < static_cast<int>(hpds.size()); i++) {
     int nhits = hpds[i].nHits();
     double ene = hpds[i].hitEnergy();
     double trkfide = hpds[i].hitEnergyTrackFiducial();
     double isolhcale = hpds[i].hcalEnergySameTowers() + hpds[i].hcalEnergyNeighborTowers();
     double isolecale = hpds[i].ecalEnergySameTowers();
     double isoltrke = hpds[i].trackEnergySameTowers() + hpds[i].trackEnergyNeighborTowers();
     if ((ene > 0 && isolhcale / ene < LooseHcalIsol_ && isolecale / ene < LooseEcalIsol_ &&
          isoltrke / ene < LooseTrackIsol_ &&
          ((trkfide > LooseHPDEne1_ && nhits >= LooseHPDHits1_) ||
           (trkfide > LooseHPDEne2_ && nhits >= LooseHPDHits2_))) ||
         (ene > 0 && isolhcale / ene < TightHcalIsol_ && isolecale / ene < TightEcalIsol_ &&
          isoltrke / ene < TightTrackIsol_ &&
          ((trkfide > TightHPDEne1_ && nhits >= TightHPDHits1_) ||
           (trkfide > TightHPDEne2_ && nhits >= TightHPDHits2_)))) {
       for (HBHEHitMap::hitmap_const_iterator it = hpds[i].beginHits(); it != hpds[i].endHits(); ++it)
         noisehits.insert(it->first);
       //      result=false;
     }
   }
 
   for (int i = 0; i < static_cast<int>(dihits.size()); i++) {
     double ene = dihits[i].hitEnergy();
     double trkfide = dihits[i].hitEnergyTrackFiducial();
     double isolhcale = dihits[i].hcalEnergySameTowers() + dihits[i].hcalEnergyNeighborTowers();
     double isolecale = dihits[i].ecalEnergySameTowers();
     double isoltrke = dihits[i].trackEnergySameTowers() + dihits[i].trackEnergyNeighborTowers();
     if ((ene > 0 && isolhcale / ene < LooseHcalIsol_ && isolecale / ene < LooseEcalIsol_ &&
          isoltrke / ene < LooseTrackIsol_ && trkfide > 0.99 * ene && trkfide > LooseDiHitEne_) ||
         (ene > 0 && isolhcale / ene < TightHcalIsol_ && isolecale / ene < TightEcalIsol_ &&
          isoltrke / ene < TightTrackIsol_ && ene > TightDiHitEne_)) {
       for (HBHEHitMap::hitmap_const_iterator it = dihits[i].beginHits(); it != dihits[i].endHits(); ++it)
         noisehits.insert(it->first);
       //      result=false;
     }
   }
 
   for (int i = 0; i < static_cast<int>(monohits.size()); i++) {
     double ene = monohits[i].hitEnergy();
     double trkfide = monohits[i].hitEnergyTrackFiducial();
     double isolhcale = monohits[i].hcalEnergySameTowers() + monohits[i].hcalEnergyNeighborTowers();
     double isolecale = monohits[i].ecalEnergySameTowers();
     double isoltrke = monohits[i].trackEnergySameTowers() + monohits[i].trackEnergyNeighborTowers();
     if ((ene > 0 && isolhcale / ene < LooseHcalIsol_ && isolecale / ene < LooseEcalIsol_ &&
          isoltrke / ene < LooseTrackIsol_ && trkfide > 0.99 * ene && trkfide > LooseMonoHitEne_) ||
         (ene > 0 && isolhcale / ene < TightHcalIsol_ && isolecale / ene < TightEcalIsol_ &&
          isoltrke / ene < TightTrackIsol_ && ene > TightMonoHitEne_)) {
       for (HBHEHitMap::hitmap_const_iterator it = monohits[i].beginHits(); it != monohits[i].endHits(); ++it)
         noisehits.insert(it->first);
       //      result=false;
     }
   }
 
   // prepare the output HBHE RecHit collection
   auto pOut = std::make_unique<HBHERecHitCollection>();
   // loop over rechits, and set the new bit you wish to use
   for (HBHERecHitCollection::const_iterator it = hbhehits_h->begin(); it != hbhehits_h->end(); ++it) {
     const HBHERecHit* hit = &(*it);
     HBHERecHit newhit(*hit);
     if (noisehits.end() != noisehits.find(hit)) {
       newhit.setFlagField(1, HcalCaloFlagLabels::HBHEIsolatedNoise);
     }
     pOut->push_back(newhit);
   }
 
   iEvent.put(std::move(pOut));
 
   return;
 }