#include <CaloExtractorByAssociator.h>

Inheritance diagram for muonisolation::CaloExtractorByAssociator:

Public Member Functions
	CaloExtractorByAssociator ()
	constructors More...

	CaloExtractorByAssociator (const edm::ParameterSet &par, edm::ConsumesCollector &&iC)

reco::IsoDeposit	deposit (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::Track &track) const override
	no-op: by design of this extractor the deposits are pulled out all at a time More...

std::vector< reco::IsoDeposit >	deposits (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::Track &track) const override
	return deposits for 3 calorimeter subdetectors (ecal, hcal, ho) – in this order More...

void	fillVetos (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::TrackCollection &tracks) override
	allows to set extra vetoes (in addition to the muon) – no-op at this point More...

	~CaloExtractorByAssociator () override
	destructor More...

Public Member Functions inherited from reco::isodeposit::IsoDepositExtractor
virtual reco::IsoDeposit	deposit (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::TrackBaseRef &track) const
	make single IsoDeposit based on trackRef as input More...

virtual reco::IsoDeposit	deposit (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::Candidate &track) const

virtual reco::IsoDeposit	deposit (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::CandidateBaseRef &track) const
	make single IsoDeposit based on a CandidateBaseRef as input More...

virtual std::vector< reco::IsoDeposit >	deposits (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::TrackBaseRef &track) const

virtual std::vector< reco::IsoDeposit >	deposits (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::Candidate &track) const

virtual std::vector< reco::IsoDeposit >	deposits (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::CandidateBaseRef &track) const

virtual void	initEvent (const edm::Event &, const edm::EventSetup &)
	perform some constly initialization More...

virtual	~IsoDepositExtractor ()
	Destructor. More...

Private Member Functions
double	noiseEcal (const CaloTower &tower) const
	Determine noise for HCAL and ECAL (take some defaults for the time being) More...

double	noiseHcal (const CaloTower &tower) const

double	noiseHOcal (const CaloTower &tower) const

double	noiseRecHit (const DetId &detId) const

Private Attributes
edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	bFieldToken_

edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	caloGeomToken_

edm::ESGetToken< EcalPFRecHitThresholds, EcalPFRecHitThresholdsRcd >	ecalPFRechitThresholdsToken_

bool	ecalRecHitThresh_

EcalPFRecHitThresholds *	ecalThresholds = nullptr

edm::ESGetToken< HcalChannelQuality, HcalChannelQualityRcd >	hcalChannelQualityToken_

HcalChannelQuality *	hcalChStatus_

HcalPFCuts *	hcalCuts = nullptr

bool	hcalCutsFromDB_

edm::ESGetToken< HcalPFCuts, HcalPFCutsRcd >	hcalCutsToken_

HcalSeverityLevelComputer *	hcalSevLvlComputer_

edm::ESGetToken< HcalSeverityLevelComputer, HcalSeverityLevelComputerRcd >	hcalSevLvlComputerToken_

HcalTopology *	hcalTopology_

edm::ESGetToken< HcalTopology, HcalRecNumberingRecord >	hcalTopologyToken_

TrackDetectorAssociator *	theAssociator

TrackAssociatorParameters *	theAssociatorParameters
	associator, its' parameters and the propagator More...

bool	theCenterConeOnCalIntersection

std::vector< std::string >	theDepositInstanceLabels
	multiple deposits: labels – expect 3 labels beginning with "e", "h", "ho" More...

std::string	theDepositLabel
	Label of deposit – suggest to set to "" (all info is in collection name anyways) More...

double	theDR_Max
	max cone size in which towers are considered More...

double	theDR_Veto_E
	cone sizes inside which the Et (towers) are not counted More...

double	theDR_Veto_H

double	theDR_Veto_HO

int	theMaxSeverityHB

int	theMaxSeverityHE

double	theNoise_EB

double	theNoise_EE

double	theNoise_HB

double	theNoise_HE

double	theNoise_HO

double	theNoiseTow_EB

double	theNoiseTow_EE

bool	thePrintTimeReport
	flag to turn on/off printing of a time report More...

std::string	thePropagatorName
	propagator name to feed into the track associator More...

MuonServiceProxy *	theService
	the event setup proxy, it takes care the services update More...

double	theThreshold_E

double	theThreshold_H

double	theThreshold_HO

bool	theUseEcalRecHitsFlag
	use towers or rec hits More...

bool	theUseHcalRecHitsFlag

bool	theUseHORecHitsFlag

std::vector< DetId >	theVetoCollection
	Vector of calo Ids to veto – not used. More...

Detailed Description

Definition at line 49 of file CaloExtractorByAssociator.h.

Constructor & Destructor Documentation

◆ CaloExtractorByAssociator() [1/2]

muonisolation::CaloExtractorByAssociator::CaloExtractorByAssociator ( )

inline

constructors

Definition at line 52 of file CaloExtractorByAssociator.h.

52 {};

◆ CaloExtractorByAssociator() [2/2]

CaloExtractorByAssociator::CaloExtractorByAssociator	(	const edm::ParameterSet &	par,
		edm::ConsumesCollector &&	iC
	)

Definition at line 40 of file CaloExtractorByAssociator.cc.

References caloGeomToken_, ecalPFRechitThresholdsToken_, ecalRecHitThresh_, edm::ParameterSet::getParameter(), hcalChannelQualityToken_, hcalCutsFromDB_, hcalCutsToken_, hcalSevLvlComputerToken_, hcalTopologyToken_, TrackAssociatorParameters::loadParameters(), MuonServiceProxy_cff::MuonServiceProxy, theAssociator, theAssociatorParameters, theService, and HLT_2024v14_cff::TrackAssociatorParameters.

     : theUseEcalRecHitsFlag(par.getParameter<bool>("UseEcalRecHitsFlag")),
       theUseHcalRecHitsFlag(par.getParameter<bool>("UseHcalRecHitsFlag")),
       theUseHORecHitsFlag(par.getParameter<bool>("UseHORecHitsFlag")),
       theDepositLabel(par.getUntrackedParameter<string>("DepositLabel")),
       theDepositInstanceLabels(par.getParameter<std::vector<std::string> >("DepositInstanceLabels")),
       thePropagatorName(par.getParameter<std::string>("PropagatorName")),
       theThreshold_E(par.getParameter<double>("Threshold_E")),
       theThreshold_H(par.getParameter<double>("Threshold_H")),
       theThreshold_HO(par.getParameter<double>("Threshold_HO")),
       theMaxSeverityHB(par.getParameter<int>("MaxSeverityHB")),
       theMaxSeverityHE(par.getParameter<int>("MaxSeverityHE")),
       theDR_Veto_E(par.getParameter<double>("DR_Veto_E")),
       theDR_Veto_H(par.getParameter<double>("DR_Veto_H")),
       theDR_Veto_HO(par.getParameter<double>("DR_Veto_HO")),
       theCenterConeOnCalIntersection(par.getParameter<bool>("CenterConeOnCalIntersection")),
       theDR_Max(par.getParameter<double>("DR_Max")),
       theNoise_EB(par.getParameter<double>("Noise_EB")),
       theNoise_EE(par.getParameter<double>("Noise_EE")),
       theNoise_HB(par.getParameter<double>("Noise_HB")),
       theNoise_HE(par.getParameter<double>("Noise_HE")),
       theNoise_HO(par.getParameter<double>("Noise_HO")),
       theNoiseTow_EB(par.getParameter<double>("NoiseTow_EB")),
       theNoiseTow_EE(par.getParameter<double>("NoiseTow_EE")),
       theService(nullptr),
       theAssociator(nullptr),
       bFieldToken_(iC.esConsumes()),
       thePrintTimeReport(par.getUntrackedParameter<bool>("PrintTimeReport")) {
   ParameterSet serviceParameters = par.getParameter<ParameterSet>("ServiceParameters");
   theService = new MuonServiceProxy(serviceParameters, edm::ConsumesCollector(iC));
 
   //theAssociatorParameters = new TrackAssociatorParameters(par.getParameter<edm::ParameterSet>("TrackAssociatorParameters"), iC);
   theAssociatorParameters = new TrackAssociatorParameters();
   theAssociatorParameters->loadParameters(par.getParameter<edm::ParameterSet>("TrackAssociatorParameters"), iC);
   theAssociator = new TrackDetectorAssociator();
 
   ecalRecHitThresh_ = par.getParameter<bool>("EcalRecHitThresh");
   hcalCutsFromDB_ = par.getParameter<bool>("HcalCutsFromDB");
 
   caloGeomToken_ = iC.esConsumes();
   ecalPFRechitThresholdsToken_ = iC.esConsumes();
   hcalCutsToken_ = iC.esConsumes();
   hcalTopologyToken_ = iC.esConsumes();
   hcalChannelQualityToken_ = iC.esConsumes(edm::ESInputTag("", "withTopo"));
   hcalSevLvlComputerToken_ = iC.esConsumes();
 }

◆ ~CaloExtractorByAssociator()

CaloExtractorByAssociator::~CaloExtractorByAssociator ( )

override

destructor

Definition at line 87 of file CaloExtractorByAssociator.cc.

References theAssociator, theAssociatorParameters, and theService.

                                                       {
   if (theAssociatorParameters)
     delete theAssociatorParameters;
   if (theService)
     delete theService;
   if (theAssociator)
     delete theAssociator;
 }

Member Function Documentation

◆ deposit()

IsoDeposit CaloExtractorByAssociator::deposit	(	const edm::Event &	ev,
		const edm::EventSetup &	evSetup,
		const reco::Track &	track
	)		const

overridevirtual

no-op: by design of this extractor the deposits are pulled out all at a time

Implements reco::isodeposit::IsoDepositExtractor.

Definition at line 104 of file CaloExtractorByAssociator.cc.

                                                                        {
   IsoDeposit::Direction muonDir(muon.eta(), muon.phi());
   IsoDeposit dep(muonDir);
 
   //   LogWarning("CaloExtractorByAssociator")
   //     <<"single deposit is not an option here\n"
   //     <<"use ::deposits --> extract all and reweight as necessary";
 
   return dep;
 }

◆ deposits()

std::vector< IsoDeposit > CaloExtractorByAssociator::deposits	(	const edm::Event &	ev,
		const edm::EventSetup &	evSetup,
		const reco::Track &	track
	)		const

overridevirtual

return deposits for 3 calorimeter subdetectors (ecal, hcal, ho) – in this order

Make separate deposits: for ECAL, HCAL, HO.

check configuration consistency could've been made at construction stage (fix later?)

this should be (eventually) set to the eta-phi of the crossing point of a straight line tangent to a muon at IP and the calorimeter

each deposit type veto is at the point of intersect with that detector

first check if the hit is inside the veto cone by dR-alone

and now pitch those in the crossed list

first check if the hit is inside the veto cone by dR-alone

and now pitch those in the crossed list

first check if the hit is inside the veto cone by dR-alone

and now pitch those in the crossed list

use calo towers

first check if the tower is inside the veto cone by dR-alone

and now pitch those in the crossed list

add the Et of the tower to deposits if it's not a vetoed; put into muonEnergy otherwise

Reimplemented from reco::isodeposit::IsoDepositExtractor.

Definition at line 118 of file CaloExtractorByAssociator.cc.

                                                                                      {
   theService->update(eventSetup);
   theAssociator->setPropagator(&*(theService->propagator(thePropagatorName)));
 
   const EcalPFRecHitThresholds* ecalThresholds = &eventSetup.getData(ecalPFRechitThresholdsToken_);
   const HcalPFCuts* hcalCuts = &eventSetup.getData(hcalCutsToken_);
   const HcalTopology* hcalTopology_ = &eventSetup.getData(hcalTopologyToken_);
   const HcalChannelQuality* hcalChStatus_ = &eventSetup.getData(hcalChannelQualityToken_);
   const HcalSeverityLevelComputer* hcalSevLvlComputer_ = &eventSetup.getData(hcalSevLvlComputerToken_);
 
   if (theDepositInstanceLabels.size() != 3) {
     LogError("MuonIsolation") << "Configuration is inconsistent: Need 3 deposit instance labels";
   }
   if (!(theDepositInstanceLabels[0].compare(0, 1, std::string("e")) == 0) ||
       !(theDepositInstanceLabels[1].compare(0, 1, std::string("h")) == 0) ||
       !(theDepositInstanceLabels[2].compare(0, 2, std::string("ho")) == 0)) {
     LogWarning("MuonIsolation")
         << "Deposit instance labels do not look like  (e*, h*, ho*):"
         << "proceed at your own risk. The extractor interprets lab0=from ecal; lab1=from hcal; lab2=from ho";
   }
 
   typedef IsoDeposit::Veto Veto;
   IsoDeposit::Direction muonDir(muon.eta(), muon.phi());
 
   IsoDeposit depEcal(muonDir);
   IsoDeposit depHcal(muonDir);
   IsoDeposit depHOcal(muonDir);
 
   auto const& bField = eventSetup.getData(bFieldToken_);
 
   reco::TransientTrack tMuon(muon, &bField);
   FreeTrajectoryState iFTS = tMuon.initialFreeState();
   TrackDetMatchInfo mInfo = theAssociator->associate(event, eventSetup, iFTS, *theAssociatorParameters);
 
   depEcal.setVeto(
       Veto(reco::isodeposit::Direction(mInfo.trkGlobPosAtEcal.eta(), mInfo.trkGlobPosAtEcal.phi()), theDR_Veto_E));
   depHcal.setVeto(
       Veto(reco::isodeposit::Direction(mInfo.trkGlobPosAtHcal.eta(), mInfo.trkGlobPosAtHcal.phi()), theDR_Veto_H));
   depHOcal.setVeto(
       Veto(reco::isodeposit::Direction(mInfo.trkGlobPosAtHO.eta(), mInfo.trkGlobPosAtHO.phi()), theDR_Veto_HO));
 
   if (theCenterConeOnCalIntersection) {
     reco::isodeposit::Direction dirTmp = depEcal.veto().vetoDir;
     double dRtmp = depEcal.veto().dR;
     depEcal = IsoDeposit(dirTmp);
     depEcal.setVeto(Veto(dirTmp, dRtmp));
 
     dirTmp = depHcal.veto().vetoDir;
     dRtmp = depHcal.veto().dR;
     depHcal = IsoDeposit(dirTmp);
     depHcal.setVeto(Veto(dirTmp, dRtmp));
 
     dirTmp = depHOcal.veto().vetoDir;
     dRtmp = depHOcal.veto().dR;
     depHOcal = IsoDeposit(dirTmp);
     depHOcal.setVeto(Veto(dirTmp, dRtmp));
   }
 
   if (theUseEcalRecHitsFlag) {
     //Ecal
     auto const& caloGeom = eventSetup.getData(caloGeomToken_);
     std::vector<const EcalRecHit*>::const_iterator eHitCI = mInfo.ecalRecHits.begin();
     for (; eHitCI != mInfo.ecalRecHits.end(); ++eHitCI) {
       const EcalRecHit* eHitCPtr = *eHitCI;
       GlobalPoint eHitPos = caloGeom.getPosition(eHitCPtr->detid());
       double deltaR2 = reco::deltaR2(muon, eHitPos);
       double cosTheta = 1. / cosh(eHitPos.eta());
       double energy = eHitCPtr->energy();
       double et = energy * cosTheta;
       if (deltaR2 > std::max(dRMax_CandDep * dRMax_CandDep, theDR_Max * theDR_Max))
         continue;
 
       if (ecalThresholds != nullptr) {  // use thresholds from rechit
         float rhThres = (ecalThresholds != nullptr) ? (*ecalThresholds)[eHitCPtr->detid()] : 0.f;
         if (energy <= rhThres)
           continue;
       } else {  // use thresholds from config
         if (et <= theThreshold_E || energy <= 3 * noiseRecHit(eHitCPtr->detid()))
           continue;
       }
 
       bool vetoHit = false;
       if (deltaR2 < std::pow(theDR_Veto_E, 2)) {
         LogDebug("RecoMuon|CaloExtractorByAssociator") << " >>> Veto ECAL hit: Calo deltaR2= " << deltaR2;
         LogDebug("RecoMuon|CaloExtractorByAssociator")
             << " >>> Calo eta phi ethcal: " << eHitPos.eta() << " " << eHitPos.phi() << " " << et;
         vetoHit = true;
       }
       if (!vetoHit) {
         for (unsigned int iH = 0; iH < mInfo.crossedEcalIds.size() && !vetoHit; ++iH) {
           if (mInfo.crossedEcalIds[iH].rawId() == eHitCPtr->detid().rawId())
             vetoHit = true;
         }
       }
 
       //check theDR_Max only here to keep vetoHits being added to the veto energy
       if (deltaR2 > std::pow(theDR_Max, 2) && !vetoHit)
         continue;
 
       if (vetoHit) {
         depEcal.addCandEnergy(et);
       } else {
         depEcal.addDeposit(reco::isodeposit::Direction(eHitPos.eta(), eHitPos.phi()), et);
       }
     }
   }
 
   if (theUseHcalRecHitsFlag) {
     //Hcal
     auto const& caloGeom = eventSetup.getData(caloGeomToken_);
     std::vector<const HBHERecHit*>::const_iterator hHitCI = mInfo.hcalRecHits.begin();
     for (; hHitCI != mInfo.hcalRecHits.end(); ++hHitCI) {
       const HBHERecHit* hHitCPtr = *hHitCI;
       GlobalPoint hHitPos = caloGeom.getPosition(hHitCPtr->detid());
       double deltaR2 = reco::deltaR2(muon, hHitPos);
       double cosTheta = 1. / cosh(hHitPos.eta());
       double energy = hHitCPtr->energy();
       double et = energy * cosTheta;
       if (deltaR2 > std::max(dRMax_CandDep * dRMax_CandDep, theDR_Max * theDR_Max))
         continue;
 
       // check Hcal Cuts from DB
       if (hcalCuts != nullptr) {
         const HcalPFCut* item = hcalCuts->getValues(hHitCPtr->id().rawId());
         if (energy <= item->noiseThreshold())
           continue;
       } else {
         if (et <= theThreshold_H || energy <= 3 * noiseRecHit(hHitCPtr->detid()))
           continue;
       }
 
       const HcalDetId hid(hHitCPtr->detid());
       DetId did = hcalTopology_->idFront(hid);
       const uint32_t flag = hHitCPtr->flags();
       const uint32_t dbflag = hcalChStatus_->getValues(did)->getValue();
       bool recovered = hcalSevLvlComputer_->recoveredRecHit(did, flag);
       int severity = hcalSevLvlComputer_->getSeverityLevel(did, flag, dbflag);
 
       const bool goodHB = hid.subdet() == HcalBarrel and (severity <= theMaxSeverityHB or recovered);
       const bool goodHE = hid.subdet() == HcalEndcap and (severity <= theMaxSeverityHE or recovered);
       if (!goodHB and !goodHE)
         continue;
 
       bool vetoHit = false;
       if (deltaR2 < std::pow(theDR_Veto_H, 2)) {
         LogDebug("RecoMuon|CaloExtractorByAssociator") << " >>> Veto HBHE hit: Calo deltaR2= " << deltaR2;
         LogDebug("RecoMuon|CaloExtractorByAssociator")
             << " >>> Calo eta phi ethcal: " << hHitPos.eta() << " " << hHitPos.phi() << " " << et;
         vetoHit = true;
       }
       if (!vetoHit) {
         for (unsigned int iH = 0; iH < mInfo.crossedHcalIds.size() && !vetoHit; ++iH) {
           if (mInfo.crossedHcalIds[iH].rawId() == hHitCPtr->detid().rawId())
             vetoHit = true;
         }
       }
 
       //check theDR_Max only here to keep vetoHits being added to the veto energy
       if (deltaR2 > std::pow(theDR_Max, 2) && !vetoHit)
         continue;
 
       if (vetoHit) {
         depHcal.addCandEnergy(et);
       } else {
         depHcal.addDeposit(reco::isodeposit::Direction(hHitPos.eta(), hHitPos.phi()), et);
       }
     }
   }
 
   if (theUseHORecHitsFlag) {
     //HOcal
     auto const& caloGeom = eventSetup.getData(caloGeomToken_);
     std::vector<const HORecHit*>::const_iterator hoHitCI = mInfo.hoRecHits.begin();
     for (; hoHitCI != mInfo.hoRecHits.end(); ++hoHitCI) {
       const HORecHit* hoHitCPtr = *hoHitCI;
       GlobalPoint hoHitPos = caloGeom.getPosition(hoHitCPtr->detid());
       double deltaR2 = reco::deltaR2(muon, hoHitPos);
       double cosTheta = 1. / cosh(hoHitPos.eta());
       double energy = hoHitCPtr->energy();
       double et = energy * cosTheta;
       if (deltaR2 > std::max(dRMax_CandDep * dRMax_CandDep, theDR_Max * theDR_Max) ||
           !(et > theThreshold_HO && energy > 3 * noiseRecHit(hoHitCPtr->detid())))
         continue;
 
       bool vetoHit = false;
       if (deltaR2 < std::pow(theDR_Veto_HO, 2)) {
         LogDebug("RecoMuon|CaloExtractorByAssociator") << " >>> Veto HO hit: Calo deltaR2= " << deltaR2;
         LogDebug("RecoMuon|CaloExtractorByAssociator")
             << " >>> Calo eta phi ethcal: " << hoHitPos.eta() << " " << hoHitPos.phi() << " " << et;
         vetoHit = true;
       }
       if (!vetoHit) {
         for (unsigned int iH = 0; iH < mInfo.crossedHOIds.size() && !vetoHit; ++iH) {
           if (mInfo.crossedHOIds[iH].rawId() == hoHitCPtr->detid().rawId())
             vetoHit = true;
         }
       }
 
       //check theDR_Max only here to keep vetoHits being added to the veto energy
       if (deltaR2 > std::pow(theDR_Max, 2) && !vetoHit)
         continue;
 
       if (vetoHit) {
         depHOcal.addCandEnergy(et);
       } else {
         depHOcal.addDeposit(reco::isodeposit::Direction(hoHitPos.eta(), hoHitPos.phi()), et);
       }
     }
   }
 
   if (!theUseEcalRecHitsFlag or !theUseHcalRecHitsFlag or !theUseHORecHitsFlag) {
     std::vector<const CaloTower*>::const_iterator calCI = mInfo.towers.begin();
     for (; calCI != mInfo.towers.end(); ++calCI) {
       const CaloTower* calCPtr = *calCI;
       double deltaR2 = reco::deltaR2(muon, *calCPtr);
       if (deltaR2 > std::max(dRMax_CandDep * dRMax_CandDep, theDR_Max * theDR_Max))
         continue;
 
       //even more copy-pasting .. need to refactor
       double etecal = calCPtr->emEt();
       double eecal = calCPtr->emEnergy();
       bool doEcal = etecal > theThreshold_E && eecal > 3 * noiseEcal(*calCPtr);
       double ethcal = calCPtr->hadEt();
       double ehcal = calCPtr->hadEnergy();
       bool doHcal = ethcal > theThreshold_H && ehcal > 3 * noiseHcal(*calCPtr);
       double ethocal = calCPtr->outerEt();
       double ehocal = calCPtr->outerEnergy();
       bool doHOcal = ethocal > theThreshold_HO && ehocal > 3 * noiseHOcal(*calCPtr);
       if ((!doEcal) && (!doHcal) && (!doHcal))
         continue;
 
       bool vetoTowerEcal = false;
       double deltar2Ecal = reco::deltaR2(mInfo.trkGlobPosAtEcal, *calCPtr);
       if (deltar2Ecal < std::pow(theDR_Veto_E, 2)) {
         LogDebug("RecoMuon|CaloExtractorByAssociator") << " >>> Veto ecal tower: Calo deltaR= " << deltar2Ecal;
         LogDebug("RecoMuon|CaloExtractorByAssociator")
             << " >>> Calo eta phi ethcal: " << calCPtr->eta() << " " << calCPtr->phi() << " " << ethcal;
         vetoTowerEcal = true;
       }
       bool vetoTowerHcal = false;
       double deltar2Hcal = reco::deltaR2(mInfo.trkGlobPosAtHcal, *calCPtr);
       if (deltar2Hcal < std::pow(theDR_Veto_H, 2)) {
         LogDebug("RecoMuon|CaloExtractorByAssociator") << " >>> Veto hcal tower: Calo deltaR= " << deltar2Hcal;
         LogDebug("RecoMuon|CaloExtractorByAssociator")
             << " >>> Calo eta phi ethcal: " << calCPtr->eta() << " " << calCPtr->phi() << " " << ethcal;
         vetoTowerHcal = true;
       }
       bool vetoTowerHOCal = false;
       double deltar2HOcal = reco::deltaR2(mInfo.trkGlobPosAtHO, *calCPtr);
       if (deltar2HOcal < std::pow(theDR_Veto_HO, 2)) {
         LogDebug("RecoMuon|CaloExtractorByAssociator") << " >>> Veto HO tower: Calo deltaR= " << deltar2HOcal;
         LogDebug("RecoMuon|CaloExtractorByAssociator")
             << " >>> Calo eta phi ethcal: " << calCPtr->eta() << " " << calCPtr->phi() << " " << ethcal;
         vetoTowerHOCal = true;
       }
 
       if (!(vetoTowerHOCal && vetoTowerHcal && vetoTowerEcal)) {
         for (unsigned int iH = 0; iH < mInfo.crossedTowerIds.size(); ++iH) {
           if (mInfo.crossedTowerIds[iH].rawId() == calCPtr->id().rawId()) {
             vetoTowerEcal = true;
             vetoTowerHcal = true;
             vetoTowerHOCal = true;
             break;
           }
         }
       }
 
       if (deltaR2 > std::pow(theDR_Max, 2) && !(vetoTowerEcal || vetoTowerHcal || vetoTowerHOCal))
         continue;
 
       reco::isodeposit::Direction towerDir(calCPtr->eta(), calCPtr->phi());
       if (doEcal and !theUseEcalRecHitsFlag) {
         if (vetoTowerEcal)
           depEcal.addCandEnergy(etecal);
         else if (deltaR2 <= std::pow(theDR_Max, 2))
           depEcal.addDeposit(towerDir, etecal);
       }
       if (doHcal and !theUseHcalRecHitsFlag) {
         if (vetoTowerHcal)
           depHcal.addCandEnergy(ethcal);
         else if (deltaR2 <= std::pow(theDR_Max, 2))
           depHcal.addDeposit(towerDir, ethcal);
       }
       if (doHOcal and !theUseHORecHitsFlag) {
         if (vetoTowerHOCal)
           depHOcal.addCandEnergy(ethocal);
         else if (deltaR2 <= std::pow(theDR_Max, 2))
           depHOcal.addDeposit(towerDir, ethocal);
       }
     }
   }
 
   std::vector<IsoDeposit> resultDeps;
   resultDeps.push_back(depEcal);
   resultDeps.push_back(depHcal);
   resultDeps.push_back(depHOcal);
 
   return resultDeps;
 }

◆ fillVetos()

void CaloExtractorByAssociator::fillVetos	(	const edm::Event &	ev,
		const edm::EventSetup &	evSetup,
		const reco::TrackCollection &	tracks
	)

overridevirtual

allows to set extra vetoes (in addition to the muon) – no-op at this point

Implements reco::isodeposit::IsoDepositExtractor.

Definition at line 96 of file CaloExtractorByAssociator.cc.

                                                                         {
   //   LogWarning("CaloExtractorByAssociator")
   //     <<"fillVetos does nothing now: IsoDeposit provides enough functionality\n"
   //     <<"to remove a deposit at/around given (eta, phi)";
 }

◆ noiseEcal()

double CaloExtractorByAssociator::noiseEcal ( const CaloTower & tower ) const

private

Determine noise for HCAL and ECAL (take some defaults for the time being)

Definition at line 437 of file CaloExtractorByAssociator.cc.

References PVValHelper::eta, hgchebackDigitizer_cfi::noise, theNoiseTow_EB, theNoiseTow_EE, and l1tHGCalTowerProducer_cfi::tower.

Referenced by deposits().

                                                                         {
   double noise = theNoiseTow_EB;
   double eta = tower.eta();
   if (fabs(eta) > 1.479)
     noise = theNoiseTow_EE;
   return noise;
 }

◆ noiseHcal()

double CaloExtractorByAssociator::noiseHcal ( const CaloTower & tower ) const

private

Definition at line 445 of file CaloExtractorByAssociator.cc.

References hgchebackDigitizer_cfi::noise, theNoise_HB, theNoise_HE, and l1tHGCalTowerProducer_cfi::tower.

Referenced by deposits().

                                                                         {
   double noise = fabs(tower.eta()) > 1.479 ? theNoise_HE : theNoise_HB;
   return noise;
 }

◆ noiseHOcal()

double CaloExtractorByAssociator::noiseHOcal ( const CaloTower & tower ) const

private

Definition at line 450 of file CaloExtractorByAssociator.cc.

References hgchebackDigitizer_cfi::noise, and theNoise_HO.

Referenced by deposits().

                                                                          {
   double noise = theNoise_HO;
   return noise;
 }

◆ noiseRecHit()

double CaloExtractorByAssociator::noiseRecHit ( const DetId & detId ) const

private

Definition at line 455 of file CaloExtractorByAssociator.cc.

References hcalRecHitTable_cff::detId, DetId::Ecal, EcalBarrel, EcalEndcap, DetId::Hcal, HcalBarrel, HcalEndcap, HcalOuter, hgchebackDigitizer_cfi::noise, theNoise_EB, theNoise_EE, theNoise_HB, theNoise_HE, and theNoise_HO.

Referenced by deposits().

                                                                       {
   double noise = 100;
   DetId::Detector det = detId.det();
   if (det == DetId::Ecal) {
     EcalSubdetector subDet = (EcalSubdetector)(detId.subdetId());
     if (subDet == EcalBarrel) {
       noise = theNoise_EB;
     } else if (subDet == EcalEndcap) {
       noise = theNoise_EE;
     }
   } else if (det == DetId::Hcal) {
     HcalSubdetector subDet = (HcalSubdetector)(detId.subdetId());
     if (subDet == HcalBarrel) {
       noise = theNoise_HB;
     } else if (subDet == HcalEndcap) {
       noise = theNoise_HE;
     } else if (subDet == HcalOuter) {
       noise = theNoise_HO;
     }
   }
   return noise;
 }

Member Data Documentation

◆ bFieldToken_

edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> muonisolation::CaloExtractorByAssociator::bFieldToken_

private

Definition at line 122 of file CaloExtractorByAssociator.h.

Referenced by deposits().

◆ caloGeomToken_

edm::ESGetToken<CaloGeometry, CaloGeometryRecord> muonisolation::CaloExtractorByAssociator::caloGeomToken_

private

Definition at line 123 of file CaloExtractorByAssociator.h.

Referenced by CaloExtractorByAssociator(), and deposits().

◆ ecalPFRechitThresholdsToken_

edm::ESGetToken<EcalPFRecHitThresholds, EcalPFRecHitThresholdsRcd> muonisolation::CaloExtractorByAssociator::ecalPFRechitThresholdsToken_

private

Definition at line 126 of file CaloExtractorByAssociator.h.

Referenced by CaloExtractorByAssociator(), and deposits().

◆ ecalRecHitThresh_

bool muonisolation::CaloExtractorByAssociator::ecalRecHitThresh_

private

Definition at line 127 of file CaloExtractorByAssociator.h.

Referenced by CaloExtractorByAssociator().

◆ ecalThresholds

EcalPFRecHitThresholds* muonisolation::CaloExtractorByAssociator::ecalThresholds = nullptr

private

Definition at line 128 of file CaloExtractorByAssociator.h.

Referenced by deposits().

◆ hcalChannelQualityToken_

edm::ESGetToken<HcalChannelQuality, HcalChannelQualityRcd> muonisolation::CaloExtractorByAssociator::hcalChannelQualityToken_

private

Definition at line 136 of file CaloExtractorByAssociator.h.

Referenced by CaloExtractorByAssociator(), and deposits().

◆ hcalChStatus_

HcalChannelQuality* muonisolation::CaloExtractorByAssociator::hcalChStatus_

private

Definition at line 139 of file CaloExtractorByAssociator.h.

Referenced by deposits().

◆ hcalCuts

HcalPFCuts* muonisolation::CaloExtractorByAssociator::hcalCuts = nullptr

private

Definition at line 133 of file CaloExtractorByAssociator.h.

Referenced by deposits().

◆ hcalCutsFromDB_

bool muonisolation::CaloExtractorByAssociator::hcalCutsFromDB_

private

Definition at line 132 of file CaloExtractorByAssociator.h.

Referenced by CaloExtractorByAssociator().

◆ hcalCutsToken_

edm::ESGetToken<HcalPFCuts, HcalPFCutsRcd> muonisolation::CaloExtractorByAssociator::hcalCutsToken_

private

Definition at line 131 of file CaloExtractorByAssociator.h.

Referenced by CaloExtractorByAssociator(), and deposits().

◆ hcalSevLvlComputer_

HcalSeverityLevelComputer* muonisolation::CaloExtractorByAssociator::hcalSevLvlComputer_

private

Definition at line 140 of file CaloExtractorByAssociator.h.

Referenced by deposits().

◆ hcalSevLvlComputerToken_

edm::ESGetToken<HcalSeverityLevelComputer, HcalSeverityLevelComputerRcd> muonisolation::CaloExtractorByAssociator::hcalSevLvlComputerToken_

private

Definition at line 137 of file CaloExtractorByAssociator.h.

Referenced by CaloExtractorByAssociator(), and deposits().

◆ hcalTopology_

HcalTopology* muonisolation::CaloExtractorByAssociator::hcalTopology_

private

Definition at line 138 of file CaloExtractorByAssociator.h.

Referenced by deposits().

◆ hcalTopologyToken_

edm::ESGetToken<HcalTopology, HcalRecNumberingRecord> muonisolation::CaloExtractorByAssociator::hcalTopologyToken_

private

Definition at line 135 of file CaloExtractorByAssociator.h.

Referenced by CaloExtractorByAssociator(), and deposits().

◆ theAssociator

TrackDetectorAssociator* muonisolation::CaloExtractorByAssociator::theAssociator

private

Definition at line 120 of file CaloExtractorByAssociator.h.

Referenced by CaloExtractorByAssociator(), deposits(), and ~CaloExtractorByAssociator().

◆ theAssociatorParameters

TrackAssociatorParameters* muonisolation::CaloExtractorByAssociator::theAssociatorParameters

private

associator, its' parameters and the propagator

Definition at line 119 of file CaloExtractorByAssociator.h.

Referenced by CaloExtractorByAssociator(), deposits(), and ~CaloExtractorByAssociator().

◆ theCenterConeOnCalIntersection

bool muonisolation::CaloExtractorByAssociator::theCenterConeOnCalIntersection

private

centers the cone on the veto direction – makes more sense for very displaced tracks like in cosmics

Definition at line 98 of file CaloExtractorByAssociator.h.

Referenced by deposits().

◆ theDepositInstanceLabels

std::vector<std::string> muonisolation::CaloExtractorByAssociator::theDepositInstanceLabels

private

multiple deposits: labels – expect 3 labels beginning with "e", "h", "ho"

Definition at line 79 of file CaloExtractorByAssociator.h.

Referenced by deposits().

◆ theDepositLabel

std::string muonisolation::CaloExtractorByAssociator::theDepositLabel

private

Label of deposit – suggest to set to "" (all info is in collection name anyways)

Definition at line 76 of file CaloExtractorByAssociator.h.

◆ theDR_Max

double muonisolation::CaloExtractorByAssociator::theDR_Max

private

max cone size in which towers are considered

Definition at line 100 of file CaloExtractorByAssociator.h.

Referenced by deposits().

◆ theDR_Veto_E

double muonisolation::CaloExtractorByAssociator::theDR_Veto_E

private

cone sizes inside which the Et (towers) are not counted

Definition at line 93 of file CaloExtractorByAssociator.h.

Referenced by deposits().

◆ theDR_Veto_H

double muonisolation::CaloExtractorByAssociator::theDR_Veto_H

private

Definition at line 94 of file CaloExtractorByAssociator.h.

Referenced by deposits().

◆ theDR_Veto_HO

double muonisolation::CaloExtractorByAssociator::theDR_Veto_HO

private

Definition at line 95 of file CaloExtractorByAssociator.h.

Referenced by deposits().

◆ theMaxSeverityHB

int muonisolation::CaloExtractorByAssociator::theMaxSeverityHB

private

Definition at line 89 of file CaloExtractorByAssociator.h.

Referenced by deposits().

◆ theMaxSeverityHE

int muonisolation::CaloExtractorByAssociator::theMaxSeverityHE

private

Definition at line 90 of file CaloExtractorByAssociator.h.

Referenced by deposits().

◆ theNoise_EB

double muonisolation::CaloExtractorByAssociator::theNoise_EB

private

the noise "sigmas" for a hit or tower to be considered consider if Energy > 3.*sigma

Definition at line 104 of file CaloExtractorByAssociator.h.

Referenced by noiseRecHit().

◆ theNoise_EE

double muonisolation::CaloExtractorByAssociator::theNoise_EE

private

Definition at line 105 of file CaloExtractorByAssociator.h.

Referenced by noiseRecHit().

◆ theNoise_HB

double muonisolation::CaloExtractorByAssociator::theNoise_HB

private

Definition at line 106 of file CaloExtractorByAssociator.h.

Referenced by noiseHcal(), and noiseRecHit().

◆ theNoise_HE

double muonisolation::CaloExtractorByAssociator::theNoise_HE

private

Definition at line 107 of file CaloExtractorByAssociator.h.

Referenced by noiseHcal(), and noiseRecHit().

◆ theNoise_HO

double muonisolation::CaloExtractorByAssociator::theNoise_HO

private

Definition at line 108 of file CaloExtractorByAssociator.h.

Referenced by noiseHOcal(), and noiseRecHit().

◆ theNoiseTow_EB

double muonisolation::CaloExtractorByAssociator::theNoiseTow_EB

private

Definition at line 109 of file CaloExtractorByAssociator.h.

Referenced by noiseEcal().

◆ theNoiseTow_EE

double muonisolation::CaloExtractorByAssociator::theNoiseTow_EE

private

Definition at line 110 of file CaloExtractorByAssociator.h.

Referenced by noiseEcal().

◆ thePrintTimeReport

bool muonisolation::CaloExtractorByAssociator::thePrintTimeReport

private

flag to turn on/off printing of a time report

Definition at line 143 of file CaloExtractorByAssociator.h.

◆ thePropagatorName

std::string muonisolation::CaloExtractorByAssociator::thePropagatorName

private

propagator name to feed into the track associator

Definition at line 82 of file CaloExtractorByAssociator.h.

Referenced by deposits().

◆ theService

MuonServiceProxy* muonisolation::CaloExtractorByAssociator::theService

private

the event setup proxy, it takes care the services update

Definition at line 116 of file CaloExtractorByAssociator.h.

Referenced by CaloExtractorByAssociator(), deposits(), and ~CaloExtractorByAssociator().

◆ theThreshold_E

double muonisolation::CaloExtractorByAssociator::theThreshold_E

private

Cone cuts and thresholds min values of Et to be included in deposits

Definition at line 86 of file CaloExtractorByAssociator.h.

Referenced by deposits().

◆ theThreshold_H

double muonisolation::CaloExtractorByAssociator::theThreshold_H

private

Definition at line 87 of file CaloExtractorByAssociator.h.

Referenced by deposits().

◆ theThreshold_HO

double muonisolation::CaloExtractorByAssociator::theThreshold_HO

private

Definition at line 88 of file CaloExtractorByAssociator.h.

Referenced by deposits().

◆ theUseEcalRecHitsFlag

bool muonisolation::CaloExtractorByAssociator::theUseEcalRecHitsFlag

private

use towers or rec hits

Definition at line 71 of file CaloExtractorByAssociator.h.

Referenced by deposits().

◆ theUseHcalRecHitsFlag

bool muonisolation::CaloExtractorByAssociator::theUseHcalRecHitsFlag

private

Definition at line 72 of file CaloExtractorByAssociator.h.

Referenced by deposits().

◆ theUseHORecHitsFlag

bool muonisolation::CaloExtractorByAssociator::theUseHORecHitsFlag

private

Definition at line 73 of file CaloExtractorByAssociator.h.

Referenced by deposits().

◆ theVetoCollection

std::vector<DetId> muonisolation::CaloExtractorByAssociator::theVetoCollection

private

Vector of calo Ids to veto – not used.

Definition at line 113 of file CaloExtractorByAssociator.h.

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ CaloExtractorByAssociator() [1/2]

◆ CaloExtractorByAssociator() [2/2]

◆ ~CaloExtractorByAssociator()

Member Function Documentation

◆ deposit()

◆ deposits()

◆ fillVetos()

◆ noiseEcal()

◆ noiseHcal()

◆ noiseHOcal()

◆ noiseRecHit()

Member Data Documentation

◆ bFieldToken_

◆ caloGeomToken_

◆ ecalPFRechitThresholdsToken_

◆ ecalRecHitThresh_

◆ ecalThresholds

◆ hcalChannelQualityToken_

◆ hcalChStatus_

◆ hcalCuts

◆ hcalCutsFromDB_

◆ hcalCutsToken_

◆ hcalSevLvlComputer_

◆ hcalSevLvlComputerToken_

◆ hcalTopology_

◆ hcalTopologyToken_

◆ theAssociator

◆ theAssociatorParameters

◆ theCenterConeOnCalIntersection

◆ theDepositInstanceLabels

◆ theDepositLabel

◆ theDR_Max

◆ theDR_Veto_E

◆ theDR_Veto_H

◆ theDR_Veto_HO

◆ theMaxSeverityHB

◆ theMaxSeverityHE

◆ theNoise_EB

◆ theNoise_EE

◆ theNoise_HB

◆ theNoise_HE

◆ theNoise_HO

◆ theNoiseTow_EB

◆ theNoiseTow_EE

◆ thePrintTimeReport

◆ thePropagatorName

◆ theService

◆ theThreshold_E

◆ theThreshold_H

◆ theThreshold_HO

◆ theUseEcalRecHitsFlag

◆ theUseHcalRecHitsFlag

◆ theUseHORecHitsFlag

◆ theVetoCollection