L1HPSPFTauBuilder Class Reference

#include <L1HPSPFTauBuilder.h>

Public Member Functions
void	addL1PFCandidates (const std::vector< l1t::PFCandidateRef > &l1PFCands)
void	buildL1PFTau ()
l1t::HPSPFTau	getL1PFTau () const
	L1HPSPFTauBuilder (const edm::ParameterSet &cfg)
void	reset ()
void	setL1PFCandProductID (const edm::ProductID &l1PFCandProductID)
void	setL1PFTauSeed (const l1t::PFCandidateRef &l1PFCandSeed)
void	setL1PFTauSeed (const reco::CaloJetRef &l1JetSeed, const std::vector< l1t::PFCandidateRef > &l1PFCands)
void	setVertex (const l1t::TkPrimaryVertexRef &primaryVertex)
	~L1HPSPFTauBuilder ()=default

Private Member Functions
l1t::PFCandidateRefVector	convertToRefVector (const std::vector< l1t::PFCandidateRef > &l1PFCands)
bool	isWithinIsolationCone (const l1t::PFCandidate &l1PFCand)
bool	isWithinSignalCone (const l1t::PFCandidate &l1PFCand)
bool	isWithinStrip (const l1t::PFCandidate &l1PFCand)

Private Attributes
bool	debug_
bool	isJetSeeded_
std::vector< l1t::PFCandidateRef >	isoAllL1PFCandidates_
std::vector< l1t::PFCandidateRef >	isoChargedHadrons_
std::vector< l1t::PFCandidateRef >	isoElectrons_
double	isolationConeSize2_
double	isolationConeSize_
std::vector< L1HPSPFTauQualityCut >	isolationQualityCutsDzCutDisabled_
std::vector< L1HPSPFTauQualityCut >	isolationQualityCutsDzCutEnabledPileup_
std::vector< L1HPSPFTauQualityCut >	isolationQualityCutsDzCutEnabledPrimary_
std::vector< l1t::PFCandidateRef >	isoMuons_
std::vector< l1t::PFCandidateRef >	isoNeutralHadrons_
std::vector< l1t::PFCandidateRef >	isoPhotons_
bool	isPFCandSeeded_
reco::CaloJetRef	l1JetSeed_
edm::ProductID	l1PFCandProductID_
l1t::PFCandidateRef	l1PFCandSeed_
l1t::HPSPFTau	l1PFTau_
double	l1PFTauSeedEta_
double	l1PFTauSeedPhi_
double	l1PFTauSeedZVtx_
double	maxSignalConeSize_
double	minSignalConeSize_
l1t::TkPrimaryVertexRef	primaryVertex_
std::vector< l1t::PFCandidateRef >	signalAllL1PFCandidates_
std::vector< l1t::PFCandidateRef >	signalChargedHadrons_
double	signalConeSize2_
double	signalConeSize_
reco::FormulaEvaluator	signalConeSizeFormula_
std::vector< l1t::PFCandidateRef >	signalElectrons_
std::vector< l1t::PFCandidateRef >	signalMuons_
std::vector< l1t::PFCandidateRef >	signalNeutralHadrons_
std::vector< l1t::PFCandidateRef >	signalPhotons_
std::vector< L1HPSPFTauQualityCut >	signalQualityCutsDzCutDisabled_
std::vector< L1HPSPFTauQualityCut >	signalQualityCutsDzCutEnabledPrimary_
std::vector< l1t::PFCandidateRef >	stripAllL1PFCandidates_
std::vector< l1t::PFCandidateRef >	stripElectrons_
reco::Particle::LorentzVector	stripP4_
std::vector< l1t::PFCandidateRef >	stripPhotons_
double	stripSizeEta_
double	stripSizePhi_
std::vector< l1t::PFCandidateRef >	sumAllL1PFCandidates_
double	sumAllL1PFCandidatesPt_
std::vector< l1t::PFCandidateRef >	sumChargedHadrons_
double	sumChargedIsoPileup_
std::vector< l1t::PFCandidateRef >	sumElectrons_
std::vector< l1t::PFCandidateRef >	sumMuons_
std::vector< l1t::PFCandidateRef >	sumNeutralHadrons_
std::vector< l1t::PFCandidateRef >	sumPhotons_
bool	useStrips_

Detailed Description

Definition at line 14 of file L1HPSPFTauBuilder.h.

Constructor & Destructor Documentation

L1HPSPFTauBuilder::L1HPSPFTauBuilder

(

const edm::ParameterSet &

cfg

)

Definition at line 10 of file L1HPSPFTauBuilder.cc.

References cms::cuda::assert(), gather_cfg::cout, debug_, edm::ParameterSet::getParameter(), isolationConeSize2_, isolationConeSize_, isolationQualityCutsDzCutDisabled_, isolationQualityCutsDzCutEnabledPileup_, isolationQualityCutsDzCutEnabledPrimary_, maxSignalConeSize_, minSignalConeSize_, readL1PFTauQualityCuts(), signalQualityCutsDzCutDisabled_, and signalQualityCutsDzCutEnabledPrimary_.

    : signalConeSizeFormula_(
          std::regex_replace(cfg.getParameter<std::string>("signalConeSize"), std::regex("pt"), "x")),
      minSignalConeSize_(cfg.getParameter<double>("minSignalConeSize")),
      maxSignalConeSize_(cfg.getParameter<double>("maxSignalConeSize")),
      useStrips_(cfg.getParameter<bool>("useStrips")),
      stripSizeEta_(cfg.getParameter<double>("stripSizeEta")),
      stripSizePhi_(cfg.getParameter<double>("stripSizePhi")),
      isolationConeSize_(cfg.getParameter<double>("isolationConeSize")),
      debug_(cfg.getUntrackedParameter<bool>("debug", false)) {
  assert(maxSignalConeSize_ >= minSignalConeSize_);

  isolationConeSize2_ = isolationConeSize_ * isolationConeSize_;

  if (debug_) {
    std::cout << "setting Quality cuts for signal PFCands:" << std::endl;
  }
  edm::ParameterSet cfg_signalQualityCuts = cfg.getParameter<edm::ParameterSet>("signalQualityCuts");
  signalQualityCutsDzCutDisabled_ = readL1PFTauQualityCuts(cfg_signalQualityCuts, "disabled", debug_);
  signalQualityCutsDzCutEnabledPrimary_ = readL1PFTauQualityCuts(cfg_signalQualityCuts, "enabled_primary", debug_);
  if (debug_) {
    std::cout << "setting Quality cuts for isolation PFCands:" << std::endl;
  }
  edm::ParameterSet cfg_isolationQualityCuts = cfg.getParameter<edm::ParameterSet>("isolationQualityCuts");
  isolationQualityCutsDzCutDisabled_ = readL1PFTauQualityCuts(cfg_isolationQualityCuts, "disabled", debug_);
  isolationQualityCutsDzCutEnabledPrimary_ =
      readL1PFTauQualityCuts(cfg_isolationQualityCuts, "enabled_primary", debug_);
  isolationQualityCutsDzCutEnabledPileup_ = readL1PFTauQualityCuts(cfg_isolationQualityCuts, "enabled_pileup", debug_);
}

L1HPSPFTauBuilder::~L1HPSPFTauBuilder ( )

default

Member Function Documentation

void L1HPSPFTauBuilder::addL1PFCandidates

(

const std::vector< l1t::PFCandidateRef > &

l1PFCands

)

Definition at line 183 of file L1HPSPFTauBuilder.cc.

References cms::cuda::assert(), l1t::PFCandidate::ChargedHadron, gather_cfg::cout, debug_, reco::deltaR(), l1t::PFCandidate::Electron, reco::FormulaEvaluator::evaluate(), isJetSeeded_, isoAllL1PFCandidates_, isoChargedHadrons_, isoElectrons_, isolationQualityCutsDzCutEnabledPileup_, isolationQualityCutsDzCutEnabledPrimary_, isoMuons_, isoNeutralHadrons_, isoPhotons_, isPFCandSeeded_, isSelected(), isWithinIsolationCone(), isWithinSignalCone(), isWithinStrip(), l1PFTauSeedEta_, l1PFTauSeedPhi_, l1PFTauSeedZVtx_, maxSignalConeSize_, minSignalConeSize_, l1t::PFCandidate::Muon, l1t::PFCandidate::NeutralHadron, l1t::PFCandidate::Photon, primaryVertex_, printPFCand(), signalAllL1PFCandidates_, signalChargedHadrons_, signalConeSize2_, signalConeSize_, signalConeSizeFormula_, signalElectrons_, signalMuons_, signalNeutralHadrons_, signalPhotons_, signalQualityCutsDzCutEnabledPrimary_, stripAllL1PFCandidates_, stripElectrons_, stripP4_, stripPhotons_, sumAllL1PFCandidates_, sumAllL1PFCandidatesPt_, sumChargedHadrons_, sumChargedIsoPileup_, sumElectrons_, sumMuons_, sumNeutralHadrons_, sumPhotons_, and useStrips_.

Referenced by HPSPFTauProducer::produce().

                                                                                       {
  if (debug_) {
    std::cout << "<L1HPSPFTauBuilder::addL1PFCandidates>:" << std::endl;
  }

  // do not build tau candidates for which no reference z-position exists,
  // as in this case charged PFCands originating from the primary (hard-scatter) interaction
  // cannot be distinguished from charged PFCands originating from pileup
  if (!(isPFCandSeeded_ || isJetSeeded_))
    return;

  for (const auto& l1PFCand : l1PFCands) {
    if (!isWithinIsolationCone(*l1PFCand))
      continue;
    sumAllL1PFCandidates_.push_back(l1PFCand);
    if (l1PFCand->id() == l1t::PFCandidate::ChargedHadron) {
      sumChargedHadrons_.push_back(l1PFCand);
    } else if (l1PFCand->id() == l1t::PFCandidate::Electron) {
      sumElectrons_.push_back(l1PFCand);
    } else if (l1PFCand->id() == l1t::PFCandidate::NeutralHadron) {
      sumNeutralHadrons_.push_back(l1PFCand);
    } else if (l1PFCand->id() == l1t::PFCandidate::Photon) {
      sumPhotons_.push_back(l1PFCand);
    } else if (l1PFCand->id() == l1t::PFCandidate::Muon) {
      sumMuons_.push_back(l1PFCand);
    }
  }

  for (const auto& l1PFCand : sumAllL1PFCandidates_) {
    sumAllL1PFCandidatesPt_ += l1PFCand->pt();
  }
  std::vector<double> emptyV;
  std::vector<double> sumAllL1PFCandidatesPt(1);
  sumAllL1PFCandidatesPt[0] = sumAllL1PFCandidatesPt_;

  signalConeSize_ = signalConeSizeFormula_.evaluate(sumAllL1PFCandidatesPt, emptyV);

  if (signalConeSize_ < minSignalConeSize_)
    signalConeSize_ = minSignalConeSize_;
  if (signalConeSize_ > maxSignalConeSize_)
    signalConeSize_ = maxSignalConeSize_;
  signalConeSize2_ = signalConeSize_ * signalConeSize_;

  for (const auto& l1PFCand : sumAllL1PFCandidates_) {
    if (debug_) {
      printPFCand(std::cout, *l1PFCand, primaryVertex_);
    }

    bool isSignalPFCand = false;
    bool isStripPFCand = false;
    bool isElectron_or_Photon =
        l1PFCand->id() == l1t::PFCandidate::Electron || l1PFCand->id() == l1t::PFCandidate::Photon;
    bool isChargedHadron = l1PFCand->id() == l1t::PFCandidate::ChargedHadron;
    if (isWithinSignalCone(*l1PFCand) && !(isChargedHadron && signalChargedHadrons_.size() > 3)) {
      isSignalPFCand = true;
    }
    if (isElectron_or_Photon && isWithinStrip(*l1PFCand)) {
      if (useStrips_) {
        isSignalPFCand = true;
      }
      isStripPFCand = true;
    }
    bool passesSignalQualityCuts = isSelected(signalQualityCutsDzCutEnabledPrimary_, *l1PFCand, l1PFTauSeedZVtx_);
    if (isSignalPFCand && passesSignalQualityCuts) {
      signalAllL1PFCandidates_.push_back(l1PFCand);
      if (l1PFCand->id() == l1t::PFCandidate::ChargedHadron) {
        signalChargedHadrons_.push_back(l1PFCand);
      } else if (l1PFCand->id() == l1t::PFCandidate::Electron) {
        signalElectrons_.push_back(l1PFCand);
      } else if (l1PFCand->id() == l1t::PFCandidate::NeutralHadron) {
        signalNeutralHadrons_.push_back(l1PFCand);
      } else if (l1PFCand->id() == l1t::PFCandidate::Photon) {
        signalPhotons_.push_back(l1PFCand);
      } else if (l1PFCand->id() == l1t::PFCandidate::Muon) {
        signalMuons_.push_back(l1PFCand);
      }
    }
    if (isStripPFCand && passesSignalQualityCuts) {
      stripAllL1PFCandidates_.push_back(l1PFCand);
      if (l1PFCand->id() == l1t::PFCandidate::Electron) {
        stripElectrons_.push_back(l1PFCand);
        stripP4_ += l1PFCand->p4();
      } else if (l1PFCand->id() == l1t::PFCandidate::Photon) {
        stripPhotons_.push_back(l1PFCand);
        stripP4_ += l1PFCand->p4();
      } else
        assert(0);
    }

    bool isIsolationPFCand = isWithinIsolationCone(*l1PFCand) && !isSignalPFCand;
    bool passesIsolationQualityCuts = isSelected(isolationQualityCutsDzCutEnabledPrimary_, *l1PFCand, l1PFTauSeedZVtx_);
    if (isIsolationPFCand && passesIsolationQualityCuts) {
      isoAllL1PFCandidates_.push_back(l1PFCand);
      if (l1PFCand->id() == l1t::PFCandidate::ChargedHadron) {
        isoChargedHadrons_.push_back(l1PFCand);
      } else if (l1PFCand->id() == l1t::PFCandidate::Electron) {
        isoElectrons_.push_back(l1PFCand);
      } else if (l1PFCand->id() == l1t::PFCandidate::NeutralHadron) {
        isoNeutralHadrons_.push_back(l1PFCand);
      } else if (l1PFCand->id() == l1t::PFCandidate::Photon) {
        isoPhotons_.push_back(l1PFCand);
      } else if (l1PFCand->id() == l1t::PFCandidate::Muon) {
        isoMuons_.push_back(l1PFCand);
      }
    }

    if (debug_) {
      std::cout << "dR = " << reco::deltaR(l1PFCand->eta(), l1PFCand->phi(), l1PFTauSeedEta_, l1PFTauSeedPhi_) << ":"
                << " isSignalPFCand = " << isSignalPFCand << ", isStripPFCand = " << isStripPFCand
                << " (passesSignalQualityCuts = " << passesSignalQualityCuts << "),"
                << " isIsolationPFCand = " << isIsolationPFCand
                << " (passesIsolationQualityCuts = " << passesIsolationQualityCuts << ")" << std::endl;
    }
  }

  for (const auto& l1PFCand : l1PFCands) {
    if (!isWithinIsolationCone(*l1PFCand))
      continue;

    if (l1PFCand->charge() != 0 && isSelected(isolationQualityCutsDzCutEnabledPileup_, *l1PFCand, l1PFTauSeedZVtx_)) {
      sumChargedIsoPileup_ += l1PFCand->pt();
    }
  }
}

void L1HPSPFTauBuilder::buildL1PFTau

(

)

Definition at line 339 of file L1HPSPFTauBuilder.cc.

References l1t::PFCandidate::ChargedHadron, convertToRefVector(), l1t::PFCandidate::Electron, edm::Ref< C, T, F >::isNull(), isoAllL1PFCandidates_, isoChargedHadrons_, isoElectrons_, isolationConeSize_, isoMuons_, isoNeutralHadrons_, isoPhotons_, l1t::HPSPFTau::kOneProng0Pi0, l1t::HPSPFTau::kOneProng1Pi0, l1t::HPSPFTau::kThreeProng0Pi0, l1t::HPSPFTau::kThreeProng1Pi0, l1JetSeed_, l1PFCandSeed_, l1PFTau_, l1t::HPSPFTau::leadChargedPFCand(), l1t::PFCandidate::Photon, primaryVertex_, l1t::HPSPFTau::setIsoAllL1PFCandidates(), l1t::HPSPFTau::setIsoChargedHadrons(), l1t::HPSPFTau::setIsoElectrons(), l1t::HPSPFTau::setisolationConeSize(), l1t::HPSPFTau::setIsoMuons(), l1t::HPSPFTau::setIsoNeutralHadrons(), l1t::HPSPFTau::setIsoPhotons(), l1t::HPSPFTau::setLeadChargedPFCand(), reco::LeafCandidate::setP4(), l1t::HPSPFTau::setPassLooseIso(), l1t::HPSPFTau::setPassMediumIso(), l1t::HPSPFTau::setPassTightIso(), l1t::HPSPFTau::setPassVLooseIso(), l1t::HPSPFTau::setPrimaryVertex(), l1t::HPSPFTau::setSeedChargedPFCand(), l1t::HPSPFTau::setSeedJet(), l1t::HPSPFTau::setSignalAllL1PFCandidates(), l1t::HPSPFTau::setSignalChargedHadrons(), l1t::HPSPFTau::setSignalConeSize(), l1t::HPSPFTau::setSignalElectrons(), l1t::HPSPFTau::setSignalMuons(), l1t::HPSPFTau::setSignalNeutralHadrons(), l1t::HPSPFTau::setSignalPhotons(), l1t::HPSPFTau::setStripAllL1PFCandidates(), l1t::HPSPFTau::setStripElectrons(), l1t::HPSPFTau::setStripP4(), l1t::HPSPFTau::setStripPhotons(), l1t::HPSPFTau::setSumAllL1PFCandidates(), l1t::HPSPFTau::setSumAllL1PFCandidatesPt(), l1t::HPSPFTau::setSumChargedHadrons(), l1t::HPSPFTau::setSumChargedIso(), l1t::HPSPFTau::setSumChargedIsoPileup(), l1t::HPSPFTau::setSumCombinedIso(), l1t::HPSPFTau::setSumElectrons(), l1t::HPSPFTau::setSumMuons(), l1t::HPSPFTau::setSumNeutralHadrons(), l1t::HPSPFTau::setSumNeutralIso(), l1t::HPSPFTau::setSumPhotons(), l1t::HPSPFTau::setTauType(), signalAllL1PFCandidates_, l1t::HPSPFTau::signalChargedHadrons(), signalChargedHadrons_, signalConeSize_, signalElectrons_, signalMuons_, signalNeutralHadrons_, signalPhotons_, edm::RefVector< C, T, F >::size(), stripAllL1PFCandidates_, stripElectrons_, stripP4_, stripPhotons_, sumAllL1PFCandidates_, sumAllL1PFCandidatesPt_, sumChargedHadrons_, l1t::HPSPFTau::sumChargedIso(), sumChargedIsoPileup_, sumElectrons_, sumMuons_, sumNeutralHadrons_, and sumPhotons_.

Referenced by HPSPFTauProducer::produce().

                                     {
  reco::Particle::LorentzVector l1PFTau_p4;
  for (const auto& l1PFCand : signalAllL1PFCandidates_) {
    if (l1PFCand->id() == l1t::PFCandidate::ChargedHadron || l1PFCand->id() == l1t::PFCandidate::Electron ||
        l1PFCand->id() == l1t::PFCandidate::Photon) {
      l1PFTau_p4 += l1PFCand->p4();
      if (l1PFCand->charge() != 0 &&
          (l1PFTau_.leadChargedPFCand().isNull() || l1PFCand->pt() > l1PFTau_.leadChargedPFCand()->pt())) {
        l1PFTau_.setLeadChargedPFCand(l1PFCand);
      }
    }
  }
  l1PFTau_.setP4(l1PFTau_p4);

  l1PFTau_.setSeedChargedPFCand(l1PFCandSeed_);
  l1PFTau_.setSeedJet(l1JetSeed_);

  l1PFTau_.setSignalAllL1PFCandidates(convertToRefVector(signalAllL1PFCandidates_));
  l1PFTau_.setSignalChargedHadrons(convertToRefVector(signalChargedHadrons_));
  l1PFTau_.setSignalElectrons(convertToRefVector(signalElectrons_));
  l1PFTau_.setSignalNeutralHadrons(convertToRefVector(signalNeutralHadrons_));
  l1PFTau_.setSignalPhotons(convertToRefVector(signalPhotons_));
  l1PFTau_.setSignalMuons(convertToRefVector(signalMuons_));

  l1PFTau_.setStripAllL1PFCandidates(convertToRefVector(stripAllL1PFCandidates_));
  l1PFTau_.setStripElectrons(convertToRefVector(stripElectrons_));
  l1PFTau_.setStripPhotons(convertToRefVector(stripPhotons_));

  l1PFTau_.setIsoAllL1PFCandidates(convertToRefVector(isoAllL1PFCandidates_));
  l1PFTau_.setIsoChargedHadrons(convertToRefVector(isoChargedHadrons_));
  l1PFTau_.setIsoElectrons(convertToRefVector(isoElectrons_));
  l1PFTau_.setIsoNeutralHadrons(convertToRefVector(isoNeutralHadrons_));
  l1PFTau_.setIsoPhotons(convertToRefVector(isoPhotons_));
  l1PFTau_.setIsoMuons(convertToRefVector(isoMuons_));

  l1PFTau_.setSumAllL1PFCandidates(convertToRefVector(sumAllL1PFCandidates_));
  l1PFTau_.setSumChargedHadrons(convertToRefVector(sumChargedHadrons_));
  l1PFTau_.setSumElectrons(convertToRefVector(sumElectrons_));
  l1PFTau_.setSumNeutralHadrons(convertToRefVector(sumNeutralHadrons_));
  l1PFTau_.setSumPhotons(convertToRefVector(sumPhotons_));
  l1PFTau_.setSumMuons(convertToRefVector(sumMuons_));

  l1PFTau_.setPrimaryVertex(primaryVertex_);

  if (l1PFTau_.signalChargedHadrons().size() > 1) {
    if (stripP4_.pt() < 5.)
      l1PFTau_.setTauType(l1t::HPSPFTau::kThreeProng0Pi0);
    else
      l1PFTau_.setTauType(l1t::HPSPFTau::kThreeProng1Pi0);
  } else {
    if (stripP4_.pt() < 5.)
      l1PFTau_.setTauType(l1t::HPSPFTau::kOneProng0Pi0);
    else
      l1PFTau_.setTauType(l1t::HPSPFTau::kOneProng1Pi0);
  }

  l1PFTau_.setStripP4(stripP4_);

  l1PFTau_.setSumAllL1PFCandidatesPt(sumAllL1PFCandidatesPt_);
  l1PFTau_.setSignalConeSize(signalConeSize_);
  l1PFTau_.setisolationConeSize(isolationConeSize_);

  double sumChargedIso = 0.;
  double sumNeutralIso = 0.;
  for (const auto& l1PFCand : isoAllL1PFCandidates_) {
    if (l1PFCand->charge() != 0) {
      sumChargedIso += l1PFCand->pt();
    } else if (l1PFCand->id() == l1t::PFCandidate::Photon) {
      sumNeutralIso += l1PFCand->pt();
    }
  }
  l1PFTau_.setSumChargedIso(sumChargedIso);
  l1PFTau_.setSumNeutralIso(sumNeutralIso);
  const double weightNeutralIso = 1.;
  const double offsetNeutralIso = 0.;
  l1PFTau_.setSumCombinedIso(sumChargedIso + weightNeutralIso * (sumNeutralIso - offsetNeutralIso));
  l1PFTau_.setSumChargedIsoPileup(sumChargedIsoPileup_);

  if (l1PFTau_.sumChargedIso() < 20.0) {
    l1PFTau_.setPassVLooseIso(true);
  }
  if (l1PFTau_.sumChargedIso() < 10.0) {
    l1PFTau_.setPassLooseIso(true);
  }
  if (l1PFTau_.sumChargedIso() < 5.0) {
    l1PFTau_.setPassMediumIso(true);
  }
  if (l1PFTau_.sumChargedIso() < 2.5) {
    l1PFTau_.setPassTightIso(true);
  }
}

l1t::PFCandidateRefVector L1HPSPFTauBuilder::convertToRefVector ( const std::vector< l1t::PFCandidateRef > &  l1PFCands )

private

Definition at line 431 of file L1HPSPFTauBuilder.cc.

References l1PFCandProductID_, and edm::RefVector< C, T, F >::push_back().

Referenced by buildL1PFTau().

                                                                                                           {
  l1t::PFCandidateRefVector l1PFCandsRefVector(l1PFCandProductID_);
  for (const auto& l1PFCand : l1PFCands) {
    l1PFCandsRefVector.push_back(l1PFCand);
  }
  return l1PFCandsRefVector;
}

l1t::HPSPFTau L1HPSPFTauBuilder::getL1PFTau ( ) const

inline

Definition at line 28 of file L1HPSPFTauBuilder.h.

References l1PFTau_.

Referenced by HPSPFTauProducer::produce().

{ return l1PFTau_; }

bool L1HPSPFTauBuilder::isWithinIsolationCone ( const l1t::PFCandidate &  l1PFCand )

private

Definition at line 330 of file L1HPSPFTauBuilder.cc.

References spr::deltaEta, srCondWrite_cfg::deltaPhi, reco::LeafCandidate::eta(), isolationConeSize2_, l1PFTauSeedEta_, l1PFTauSeedPhi_, and reco::LeafCandidate::phi().

Referenced by addL1PFCandidates().

                                                                            {
  double deltaEta = l1PFCand.eta() - l1PFTauSeedEta_;
  double deltaPhi = l1PFCand.phi() - l1PFTauSeedPhi_;
  if ((deltaEta * deltaEta + deltaPhi * deltaPhi) < isolationConeSize2_)
    return true;
  else
    return false;
}

bool L1HPSPFTauBuilder::isWithinSignalCone ( const l1t::PFCandidate &  l1PFCand )

private

Definition at line 310 of file L1HPSPFTauBuilder.cc.

References spr::deltaEta, srCondWrite_cfg::deltaPhi, reco::LeafCandidate::eta(), isJetSeeded_, isPFCandSeeded_, l1PFTauSeedEta_, l1PFTauSeedPhi_, reco::LeafCandidate::phi(), and signalConeSize2_.

Referenced by addL1PFCandidates().

                                                                         {
  if (isPFCandSeeded_ || isJetSeeded_) {
    double deltaEta = l1PFCand.eta() - l1PFTauSeedEta_;
    double deltaPhi = l1PFCand.phi() - l1PFTauSeedPhi_;
    if ((deltaEta * deltaEta + deltaPhi * deltaPhi) < signalConeSize2_)
      return true;
  }
  return false;
}

bool L1HPSPFTauBuilder::isWithinStrip ( const l1t::PFCandidate &  l1PFCand )

private

Definition at line 320 of file L1HPSPFTauBuilder.cc.

References spr::deltaEta, srCondWrite_cfg::deltaPhi, reco::LeafCandidate::eta(), isJetSeeded_, isPFCandSeeded_, l1PFTauSeedEta_, l1PFTauSeedPhi_, reco::LeafCandidate::phi(), stripSizeEta_, and stripSizePhi_.

Referenced by addL1PFCandidates().

                                                                    {
  if (isPFCandSeeded_ || isJetSeeded_) {
    double deltaEta = l1PFCand.eta() - l1PFTauSeedEta_;
    double deltaPhi = l1PFCand.phi() - l1PFTauSeedPhi_;
    if (std::fabs(deltaEta) < stripSizeEta_ && std::fabs(deltaPhi) < stripSizePhi_)
      return true;
  }
  return false;
}

void L1HPSPFTauBuilder::reset

(

void

)

Definition at line 40 of file L1HPSPFTauBuilder.cc.

References isJetSeeded_, isoAllL1PFCandidates_, isoChargedHadrons_, isoElectrons_, isoMuons_, isoNeutralHadrons_, isoPhotons_, isPFCandSeeded_, l1JetSeed_, l1PFCandProductID_, l1PFCandSeed_, l1PFTau_, l1PFTauSeedEta_, l1PFTauSeedPhi_, l1PFTauSeedZVtx_, primaryVertex_, signalAllL1PFCandidates_, signalChargedHadrons_, signalConeSize2_, signalConeSize_, signalElectrons_, signalMuons_, signalNeutralHadrons_, signalPhotons_, stripAllL1PFCandidates_, stripElectrons_, stripP4_, stripPhotons_, sumAllL1PFCandidates_, sumAllL1PFCandidatesPt_, sumChargedHadrons_, sumChargedIsoPileup_, sumElectrons_, sumMuons_, sumNeutralHadrons_, and sumPhotons_.

Referenced by HPSPFTauProducer::produce().

                              {
  signalConeSize_ = 0.;
  signalConeSize2_ = 0.;

  l1PFCandProductID_ = edm::ProductID();
  isPFCandSeeded_ = false;
  l1PFCandSeed_ = l1t::PFCandidateRef();
  isJetSeeded_ = false;
  l1JetSeed_ = reco::CaloJetRef();
  l1PFTauSeedEta_ = 0.;
  l1PFTauSeedPhi_ = 0.;
  l1PFTauSeedZVtx_ = 0.;
  sumAllL1PFCandidatesPt_ = 0.;
  primaryVertex_ = l1t::TkPrimaryVertexRef();
  l1PFTau_ = l1t::HPSPFTau();

  stripP4_ = reco::Particle::LorentzVector(0., 0., 0., 0.);

  signalAllL1PFCandidates_.clear();
  signalChargedHadrons_.clear();
  signalElectrons_.clear();
  signalNeutralHadrons_.clear();
  signalPhotons_.clear();
  signalMuons_.clear();

  stripAllL1PFCandidates_.clear();
  stripElectrons_.clear();
  stripPhotons_.clear();

  isoAllL1PFCandidates_.clear();
  isoChargedHadrons_.clear();
  isoElectrons_.clear();
  isoNeutralHadrons_.clear();
  isoPhotons_.clear();
  isoMuons_.clear();

  sumAllL1PFCandidates_.clear();
  sumChargedHadrons_.clear();
  sumElectrons_.clear();
  sumNeutralHadrons_.clear();
  sumPhotons_.clear();
  sumMuons_.clear();

  sumChargedIsoPileup_ = 0.;
}

void L1HPSPFTauBuilder::setL1PFCandProductID

(

const edm::ProductID &

l1PFCandProductID

)

Definition at line 86 of file L1HPSPFTauBuilder.cc.

References l1PFCandProductID_.

Referenced by HPSPFTauProducer::produce().

                                                                                  {
  l1PFCandProductID_ = l1PFCandProductID;
}

void L1HPSPFTauBuilder::setL1PFTauSeed

(

const l1t::PFCandidateRef &

l1PFCandSeed

)

Definition at line 92 of file L1HPSPFTauBuilder.cc.

References gather_cfg::cout, debug_, edm::Ref< C, T, F >::isNonnull(), isPFCandSeeded_, l1PFCandSeed_, l1PFTauSeedEta_, l1PFTauSeedPhi_, l1PFTauSeedZVtx_, primaryVertex_, and printPFCand().

Referenced by HPSPFTauProducer::produce().

                                                                            {
  if (debug_) {
    std::cout << "<L1HPSPFTauBuilder::setL1PFTauSeed>:" << std::endl;
    std::cout << "seeding HPSPFTau with ChargedPFCand:";
    printPFCand(std::cout, *l1PFCandSeed, primaryVertex_);
  }

  l1PFCandSeed_ = l1PFCandSeed;
  l1PFTauSeedEta_ = l1PFCandSeed->eta();
  l1PFTauSeedPhi_ = l1PFCandSeed->phi();
  if (l1PFCandSeed->charge() != 0 && l1PFCandSeed->pfTrack().isNonnull()) {
    l1PFTauSeedZVtx_ = l1PFCandSeed->pfTrack()->vertex().z();
    isPFCandSeeded_ = true;
  }
}

void L1HPSPFTauBuilder::setL1PFTauSeed	(	const reco::CaloJetRef &	l1JetSeed,
		const std::vector< l1t::PFCandidateRef > &	l1PFCands
	)

Definition at line 147 of file L1HPSPFTauBuilder.cc.

References l1t::PFCandidate::ChargedHadron, gather_cfg::cout, debug_, reco::deltaR(), l1t::PFCandidate::Electron, isJetSeeded_, l1JetSeed_, l1PFTauSeedEta_, l1PFTauSeedPhi_, l1PFTauSeedZVtx_, l1t::PFCandidate::Muon, and l1t::PFCandidate::Photon.

                                                                                      {
  if (debug_) {
    std::cout << "<L1HPSPFTauBuilder::setL1PFTauSeed>:" << std::endl;
    std::cout << "seeding HPSPFTau with Jet:";
    std::cout << " pT = " << l1JetSeed->pt() << ", eta = " << l1JetSeed->eta() << ", phi = " << l1JetSeed->phi()
              << std::endl;
  }

  l1JetSeed_ = l1JetSeed;
  reco::Candidate::LorentzVector l1PFTauSeed_p4;
  float l1PFTauSeedZVtx = 0.;
  bool l1PFTauSeed_hasVtx = false;
  float max_chargedPFCand_pt = -1.;
  for (const auto& l1PFCand : l1PFCands) {
    double dR = reco::deltaR(l1PFCand->eta(), l1PFCand->phi(), l1JetSeed->eta(), l1JetSeed->phi());
    if (dR > 0.4)
      continue;
    if (l1PFCand->id() == l1t::PFCandidate::ChargedHadron || l1PFCand->id() == l1t::PFCandidate::Electron ||
        l1PFCand->id() == l1t::PFCandidate::Photon || l1PFCand->id() == l1t::PFCandidate::Muon) {
      l1PFTauSeed_p4 += l1PFCand->p4();
      if (l1PFCand->charge() != 0 && l1PFCand->pfTrack().isNonnull() && l1PFCand->pt() > max_chargedPFCand_pt) {
        l1PFTauSeedZVtx = l1PFCand->pfTrack()->vertex().z();
        l1PFTauSeed_hasVtx = true;
        max_chargedPFCand_pt = l1PFCand->pt();
      }
    }
  }
  if (l1PFTauSeed_p4.pt() > 1. && l1PFTauSeed_hasVtx) {
    l1PFTauSeedEta_ = l1PFTauSeed_p4.eta();
    l1PFTauSeedPhi_ = l1PFTauSeed_p4.phi();
    l1PFTauSeedZVtx_ = l1PFTauSeedZVtx;
    isJetSeeded_ = true;
  }
}

void L1HPSPFTauBuilder::setVertex

(

const l1t::TkPrimaryVertexRef &

primaryVertex

)

Definition at line 90 of file L1HPSPFTauBuilder.cc.

References beam_dqm_sourceclient-live_cfg::primaryVertex, and primaryVertex_.

Referenced by HPSPFTauProducer::produce().

{ primaryVertex_ = primaryVertex; }

Member Data Documentation

bool L1HPSPFTauBuilder::debug_

private

Definition at line 97 of file L1HPSPFTauBuilder.h.

Referenced by addL1PFCandidates(), L1HPSPFTauBuilder(), and setL1PFTauSeed().

bool L1HPSPFTauBuilder::isJetSeeded_

private

Definition at line 59 of file L1HPSPFTauBuilder.h.

Referenced by addL1PFCandidates(), isWithinSignalCone(), isWithinStrip(), reset(), and setL1PFTauSeed().

std::vector<l1t::PFCandidateRef> L1HPSPFTauBuilder::isoAllL1PFCandidates_

private

Definition at line 81 of file L1HPSPFTauBuilder.h.

Referenced by addL1PFCandidates(), buildL1PFTau(), and reset().

std::vector<l1t::PFCandidateRef> L1HPSPFTauBuilder::isoChargedHadrons_

private

Definition at line 82 of file L1HPSPFTauBuilder.h.

Referenced by addL1PFCandidates(), buildL1PFTau(), and reset().

std::vector<l1t::PFCandidateRef> L1HPSPFTauBuilder::isoElectrons_

private

Definition at line 83 of file L1HPSPFTauBuilder.h.

Referenced by addL1PFCandidates(), buildL1PFTau(), and reset().

double L1HPSPFTauBuilder::isolationConeSize2_

private

Definition at line 49 of file L1HPSPFTauBuilder.h.

Referenced by isWithinIsolationCone(), and L1HPSPFTauBuilder().

double L1HPSPFTauBuilder::isolationConeSize_

private

Definition at line 48 of file L1HPSPFTauBuilder.h.

Referenced by buildL1PFTau(), and L1HPSPFTauBuilder().

std::vector<L1HPSPFTauQualityCut> L1HPSPFTauBuilder::isolationQualityCutsDzCutDisabled_

private

Definition at line 53 of file L1HPSPFTauBuilder.h.

Referenced by L1HPSPFTauBuilder().

std::vector<L1HPSPFTauQualityCut> L1HPSPFTauBuilder::isolationQualityCutsDzCutEnabledPileup_

private

Definition at line 55 of file L1HPSPFTauBuilder.h.

Referenced by addL1PFCandidates(), and L1HPSPFTauBuilder().

std::vector<L1HPSPFTauQualityCut> L1HPSPFTauBuilder::isolationQualityCutsDzCutEnabledPrimary_

private

Definition at line 54 of file L1HPSPFTauBuilder.h.

Referenced by addL1PFCandidates(), and L1HPSPFTauBuilder().

std::vector<l1t::PFCandidateRef> L1HPSPFTauBuilder::isoMuons_

private

Definition at line 86 of file L1HPSPFTauBuilder.h.

Referenced by addL1PFCandidates(), buildL1PFTau(), and reset().

std::vector<l1t::PFCandidateRef> L1HPSPFTauBuilder::isoNeutralHadrons_

private

Definition at line 84 of file L1HPSPFTauBuilder.h.

Referenced by addL1PFCandidates(), buildL1PFTau(), and reset().

std::vector<l1t::PFCandidateRef> L1HPSPFTauBuilder::isoPhotons_

private

Definition at line 85 of file L1HPSPFTauBuilder.h.

Referenced by addL1PFCandidates(), buildL1PFTau(), and reset().

bool L1HPSPFTauBuilder::isPFCandSeeded_

private

Definition at line 57 of file L1HPSPFTauBuilder.h.

Referenced by addL1PFCandidates(), isWithinSignalCone(), isWithinStrip(), reset(), and setL1PFTauSeed().

reco::CaloJetRef L1HPSPFTauBuilder::l1JetSeed_

private

Definition at line 60 of file L1HPSPFTauBuilder.h.

Referenced by buildL1PFTau(), reset(), and setL1PFTauSeed().

edm::ProductID L1HPSPFTauBuilder::l1PFCandProductID_

private

Definition at line 56 of file L1HPSPFTauBuilder.h.

Referenced by convertToRefVector(), reset(), and setL1PFCandProductID().

l1t::PFCandidateRef L1HPSPFTauBuilder::l1PFCandSeed_

private

Definition at line 58 of file L1HPSPFTauBuilder.h.

Referenced by buildL1PFTau(), reset(), and setL1PFTauSeed().

l1t::HPSPFTau L1HPSPFTauBuilder::l1PFTau_

private

Definition at line 66 of file L1HPSPFTauBuilder.h.

Referenced by buildL1PFTau(), getL1PFTau(), and reset().

double L1HPSPFTauBuilder::l1PFTauSeedEta_

private

Definition at line 61 of file L1HPSPFTauBuilder.h.

Referenced by addL1PFCandidates(), isWithinIsolationCone(), isWithinSignalCone(), isWithinStrip(), reset(), and setL1PFTauSeed().

double L1HPSPFTauBuilder::l1PFTauSeedPhi_

private

Definition at line 62 of file L1HPSPFTauBuilder.h.

Referenced by addL1PFCandidates(), isWithinIsolationCone(), isWithinSignalCone(), isWithinStrip(), reset(), and setL1PFTauSeed().

double L1HPSPFTauBuilder::l1PFTauSeedZVtx_

private

Definition at line 63 of file L1HPSPFTauBuilder.h.

Referenced by addL1PFCandidates(), reset(), and setL1PFTauSeed().

double L1HPSPFTauBuilder::maxSignalConeSize_

private

Definition at line 42 of file L1HPSPFTauBuilder.h.

Referenced by addL1PFCandidates(), and L1HPSPFTauBuilder().

double L1HPSPFTauBuilder::minSignalConeSize_

private

Definition at line 41 of file L1HPSPFTauBuilder.h.

Referenced by addL1PFCandidates(), and L1HPSPFTauBuilder().

l1t::TkPrimaryVertexRef L1HPSPFTauBuilder::primaryVertex_

private

Definition at line 65 of file L1HPSPFTauBuilder.h.

Referenced by addL1PFCandidates(), buildL1PFTau(), reset(), setL1PFTauSeed(), and setVertex().

std::vector<l1t::PFCandidateRef> L1HPSPFTauBuilder::signalAllL1PFCandidates_

private

Definition at line 70 of file L1HPSPFTauBuilder.h.

Referenced by addL1PFCandidates(), buildL1PFTau(), and reset().

std::vector<l1t::PFCandidateRef> L1HPSPFTauBuilder::signalChargedHadrons_

private

Definition at line 71 of file L1HPSPFTauBuilder.h.

Referenced by addL1PFCandidates(), buildL1PFTau(), and reset().

double L1HPSPFTauBuilder::signalConeSize2_

private

Definition at line 40 of file L1HPSPFTauBuilder.h.

Referenced by addL1PFCandidates(), isWithinSignalCone(), and reset().

double L1HPSPFTauBuilder::signalConeSize_

private

Definition at line 39 of file L1HPSPFTauBuilder.h.

Referenced by addL1PFCandidates(), buildL1PFTau(), and reset().

reco::FormulaEvaluator L1HPSPFTauBuilder::signalConeSizeFormula_

private

Definition at line 37 of file L1HPSPFTauBuilder.h.

Referenced by addL1PFCandidates().

std::vector<l1t::PFCandidateRef> L1HPSPFTauBuilder::signalElectrons_

private

Definition at line 72 of file L1HPSPFTauBuilder.h.

Referenced by addL1PFCandidates(), buildL1PFTau(), and reset().

std::vector<l1t::PFCandidateRef> L1HPSPFTauBuilder::signalMuons_

private

Definition at line 75 of file L1HPSPFTauBuilder.h.

Referenced by addL1PFCandidates(), buildL1PFTau(), and reset().

std::vector<l1t::PFCandidateRef> L1HPSPFTauBuilder::signalNeutralHadrons_

private

Definition at line 73 of file L1HPSPFTauBuilder.h.

Referenced by addL1PFCandidates(), buildL1PFTau(), and reset().

std::vector<l1t::PFCandidateRef> L1HPSPFTauBuilder::signalPhotons_

private

Definition at line 74 of file L1HPSPFTauBuilder.h.

Referenced by addL1PFCandidates(), buildL1PFTau(), and reset().

std::vector<L1HPSPFTauQualityCut> L1HPSPFTauBuilder::signalQualityCutsDzCutDisabled_

private

Definition at line 51 of file L1HPSPFTauBuilder.h.

Referenced by L1HPSPFTauBuilder().

std::vector<L1HPSPFTauQualityCut> L1HPSPFTauBuilder::signalQualityCutsDzCutEnabledPrimary_

private

Definition at line 52 of file L1HPSPFTauBuilder.h.

Referenced by addL1PFCandidates(), and L1HPSPFTauBuilder().

std::vector<l1t::PFCandidateRef> L1HPSPFTauBuilder::stripAllL1PFCandidates_

private

Definition at line 77 of file L1HPSPFTauBuilder.h.

Referenced by addL1PFCandidates(), buildL1PFTau(), and reset().

std::vector<l1t::PFCandidateRef> L1HPSPFTauBuilder::stripElectrons_

private

Definition at line 78 of file L1HPSPFTauBuilder.h.

Referenced by addL1PFCandidates(), buildL1PFTau(), and reset().

reco::Particle::LorentzVector L1HPSPFTauBuilder::stripP4_

private

Definition at line 68 of file L1HPSPFTauBuilder.h.

Referenced by addL1PFCandidates(), buildL1PFTau(), and reset().

std::vector<l1t::PFCandidateRef> L1HPSPFTauBuilder::stripPhotons_

private

Definition at line 79 of file L1HPSPFTauBuilder.h.

Referenced by addL1PFCandidates(), buildL1PFTau(), and reset().

double L1HPSPFTauBuilder::stripSizeEta_

private

Definition at line 45 of file L1HPSPFTauBuilder.h.

Referenced by isWithinStrip().

double L1HPSPFTauBuilder::stripSizePhi_

private

Definition at line 46 of file L1HPSPFTauBuilder.h.

Referenced by isWithinStrip().

std::vector<l1t::PFCandidateRef> L1HPSPFTauBuilder::sumAllL1PFCandidates_

private

Definition at line 88 of file L1HPSPFTauBuilder.h.

Referenced by addL1PFCandidates(), buildL1PFTau(), and reset().

double L1HPSPFTauBuilder::sumAllL1PFCandidatesPt_

private

Definition at line 64 of file L1HPSPFTauBuilder.h.

Referenced by addL1PFCandidates(), buildL1PFTau(), and reset().

std::vector<l1t::PFCandidateRef> L1HPSPFTauBuilder::sumChargedHadrons_

private

Definition at line 89 of file L1HPSPFTauBuilder.h.

Referenced by addL1PFCandidates(), buildL1PFTau(), and reset().

double L1HPSPFTauBuilder::sumChargedIsoPileup_

private

Definition at line 95 of file L1HPSPFTauBuilder.h.

Referenced by addL1PFCandidates(), buildL1PFTau(), and reset().

std::vector<l1t::PFCandidateRef> L1HPSPFTauBuilder::sumElectrons_

private

Definition at line 90 of file L1HPSPFTauBuilder.h.

Referenced by addL1PFCandidates(), buildL1PFTau(), and reset().

std::vector<l1t::PFCandidateRef> L1HPSPFTauBuilder::sumMuons_

private

Definition at line 93 of file L1HPSPFTauBuilder.h.

Referenced by addL1PFCandidates(), buildL1PFTau(), and reset().

std::vector<l1t::PFCandidateRef> L1HPSPFTauBuilder::sumNeutralHadrons_

private

Definition at line 91 of file L1HPSPFTauBuilder.h.

Referenced by addL1PFCandidates(), buildL1PFTau(), and reset().

std::vector<l1t::PFCandidateRef> L1HPSPFTauBuilder::sumPhotons_

private

Definition at line 92 of file L1HPSPFTauBuilder.h.

Referenced by addL1PFCandidates(), buildL1PFTau(), and reset().

bool L1HPSPFTauBuilder::useStrips_

private

Definition at line 44 of file L1HPSPFTauBuilder.h.

Referenced by addL1PFCandidates().