reco::tau Namespace Reference

Namespaces
	cone
	disc
	helpers
	qcuts
	xclean

Classes
class	AssociationMatchRefSelector
class	Combinatoric
class	CombinatoricGenerator
class	CombinatoricIterator
class	PATTauDiscriminationByMVAIsolationRun2
class	PFRecoTauChargedHadronBuilderPlugin
class	PFRecoTauChargedHadronFromPFCandidatePlugin
class	PFRecoTauChargedHadronFromTrackPlugin
class	PFRecoTauChargedHadronQualityPlugin
class	PFRecoTauChargedHadronStringQuality
class	PFRecoTauDiscriminationByMVAIsolationRun2
class	PFRecoTauEnergyAlgorithmPlugin
class	PFRecoTauMassPlugin
class	RecoTauBuilderCombinatoricPlugin
class	RecoTauBuilderConePlugin
class	RecoTauBuilderPlugin
class	RecoTauChargeCleanerPlugin
class	RecoTauChargedHadronMultiplicityCleanerPlugin
class	RecoTauCleanerPlugin
class	RecoTauConstructor
class	RecoTauDiscriminantCleanerPlugin
class	RecoTauDiscriminantFromDiscriminator
class	RecoTauDiscriminantFunctionPlugin
class	RecoTauDiscriminantInvariantWidth
class	RecoTauDiscriminantPlugin
class	RecoTauDiscriminantVectorFunctionPlugin
class	RecoTauDiscriminationBinnedIsolation
class	RecoTauDiscriminationBinnedIsolationImpl
class	RecoTauElectronRejectionPlugin
class	RecoTauEnergyRecoveryPlugin2
class	RecoTauEventHolderPlugin
class	RecoTauImpactParameterSignificancePlugin
class	RecoTauIsolationMasking
class	RecoTauLexicographicalRanking
class	RecoTauModifierPlugin
class	RecoTauNamedPlugin
class	RecoTauObjectEmbedder
class	RecoTauPhotonFilter
class	RecoTauPiZeroBuilderPlugin
class	RecoTauPiZeroCombinatoricPlugin
class	RecoTauPiZeroQualityPlugin
class	RecoTauPiZeroStringQuality
class	RecoTauPiZeroStripPlugin
class	RecoTauPiZeroStripPlugin2
class	RecoTauPiZeroStripPlugin3
class	RecoTauPiZeroTrivialPlugin
class	RecoTauQualityCuts
class	RecoTauRandomCleanerPlugin
class	RecoTauSoftTwoProngTausCleanerPlugin
class	RecoTauStringCleanerPlugin
class	RecoTauTagInfoWorkaroundModifer
class	RecoTauTwoProngFilter
class	RecoTauVertexAssociator
class	SortByDescendingPt
class	SortPFCandsDescendingPt

Typedefs
typedef tau::CombinatoricGenerator< ChargedHadronList >	ChargedHadronCombo
typedef std::vector< reco::PFRecoTauChargedHadron >	ChargedHadronList
typedef tau::CombinatoricGenerator< PiZeroList >	PiZeroCombo
typedef std::vector< RecoTauPiZero >	PiZeroList

Enumerations
enum	{   kOldDMwoLT, kOldDMwLT, kNewDMwoLT, kNewDMwLT,   kDBoldDMwLT, kDBnewDMwLT, kPWoldDMwLT, kPWnewDMwLT,   kDBoldDMwLTwGJ, kDBnewDMwLTwGJ }

Functions
template<typename RefVectorType , typename BaseView >
RefVectorType	castView (const edm::Handle< BaseView > &view)
	Convert a BaseView (View<T>) to a TRefVector. More...
unsigned int	chargedHadronsInDecayMode (PFTau::hadronicDecayMode mode)
	Reverse mapping of decay modes into multiplicities. More...
template<typename Container , class OverlapFunction >
Container	cleanOverlaps (const Container &dirty)
reco::Candidate::LorentzVector	compChargedHadronP4fromPThetaPhi (double, double, double)
reco::Candidate::LorentzVector	compChargedHadronP4fromPxPyPz (double, double, double)
std::string	discPluginName (const std::string &mvaName)
float	eratio (const reco::PFTau &tau)
	return ratio of energy in ECAL over sum of energy in ECAL and HCAL More...
float	eratio (const pat::Tau &tau)
std::pair< edm::ParameterSet, edm::ParameterSet >	factorizePUQCuts (const edm::ParameterSet &inputSet)
bool	fillIsoMVARun2Inputs (float *mvaInput, const pat::Tau &tau, int mvaOpt, const std::string &nameCharged, const std::string &nameNeutral, const std::string &namePu, const std::string &nameOutside, const std::string &nameFootprint)
template<typename Iterator >
std::vector< PFCandidatePtr >	filterPFCandidates (const Iterator &begin, const Iterator &end, int particleId, bool sort=true)
std::vector< PFCandidatePtr >	flattenPiZeros (const std::vector< RecoTauPiZero >::const_iterator &, const std::vector< RecoTauPiZero >::const_iterator &)
	Flatten a list of pi zeros into a list of there constituent PFCandidates. More...
std::vector< PFCandidatePtr >	flattenPiZeros (const std::vector< RecoTauPiZero > &)
PFTau::hadronicDecayMode	getDecayMode (const reco::GenJet *genJet)
float	lead_track_chi2 (const reco::PFTau &tau)
	return chi2 of the leading track ==> deprecated? <== More...
template<typename InputIterator >
InputIterator	leadPFCand (InputIterator begin, InputIterator end)
unsigned int	n_photons_total (const reco::PFTau &tau)
	return total number of pf photon candidates with pT>500 MeV, which are associated to signal More...
unsigned int	n_photons_total (const pat::Tau &tau)
std::vector< PFCandidatePtr >	pfCandidates (const PFJet &jet, int particleId, bool sort=true)
std::vector< PFCandidatePtr >	pfCandidates (const PFJet &jet, const std::vector< int > &particleIds, bool sort=true)
	Extract pfCandidates of a that match a list of particle Ids from a PFJet. More...
std::vector< PFCandidatePtr >	pfChargedCands (const PFJet &jet, bool sort=true)
	Extract all non-neutral candidates from a PFJet. More...
std::vector< PFCandidatePtr >	pfGammas (const PFJet &jet, bool sort=true)
	Extract all pfGammas from a PFJet. More...
unsigned int	piZerosInDecayMode (PFTau::hadronicDecayMode mode)
float	pt_weighted_deta_strip (const reco::PFTau &tau, int dm)
float	pt_weighted_deta_strip (const pat::Tau &tau, int dm)
float	pt_weighted_dphi_strip (const reco::PFTau &tau, int dm)
float	pt_weighted_dphi_strip (const pat::Tau &tau, int dm)
float	pt_weighted_dr_iso (const reco::PFTau &tau, int dm)
float	pt_weighted_dr_iso (const pat::Tau &tau, int dm)
float	pt_weighted_dr_signal (const reco::PFTau &tau, int dm)
float	pt_weighted_dr_signal (const pat::Tau &tau, int dm)
float	pt_weighted_dx (const reco::PFTau &tau, int mode=0, int var=0, int decaymode=-1)
float	pt_weighted_dx (const pat::Tau &tau, int mode=0, int var=0, int decaymode=-1)
template<typename T >
float	pt_weighted_dx_T (const typename T::Tau_t &tau, int mode, int var, int decaymode)
static const edm::ParameterSetDescriptionFillerPluginFactory::PMaker< edm::ParameterSetDescriptionFiller< PATTauDiscriminationByMVAIsolationRun2 > >	s_filler__LINE__ ("PATTauDiscriminationByMVAIsolationRun2")
static const edm::ParameterSetDescriptionFillerPluginFactory::PMaker< edm::ParameterSetDescriptionFiller< PFRecoTauDiscriminationByMVAIsolationRun2 > >	s_filler__LINE__ ("PFRecoTauDiscriminationByMVAIsolationRun2")
static const edm::MakerPluginFactory::PMaker< edm::WorkerMaker< PATTauDiscriminationByMVAIsolationRun2 > >	s_maker__LINE__ ("PATTauDiscriminationByMVAIsolationRun2")
static const edm::MakerPluginFactory::PMaker< edm::WorkerMaker< PFRecoTauDiscriminationByMVAIsolationRun2 > >	s_maker__LINE__ ("PFRecoTauDiscriminationByMVAIsolationRun2")
void	setChargedHadronP4 (reco::PFRecoTauChargedHadron &chargedHadron, double scaleFactor_neutralPFCands=1.0)
template<typename InputIterator >
int	sumPFCandCharge (InputIterator begin, InputIterator end)
	Sum the charge of a collection of PFCandidates. More...
template<typename InputIterator >
reco::Candidate::LorentzVector	sumPFCandP4 (InputIterator begin, InputIterator end)
template<typename InputIterator >
double	sumPFCandPt (InputIterator begin, InputIterator end)
	Sum the pT of a collection of PFCandidates. More...
template<typename InputIterator , typename FunctionPtr , typename ReturnType >
ReturnType	sumPFVector (InputIterator begin, InputIterator end, FunctionPtr func, ReturnType init)
	Sum the four vectors in a collection of PFCandidates. More...
template<typename InputIterator >
InputIterator	takeNElements (const InputIterator &begin, const InputIterator &end, size_t N)
PFTau::hadronicDecayMode	translateDecayMode (unsigned int nCharged, unsigned int nPiZero)
PFTau::hadronicDecayMode	translateGenDecayModeToReco (const std::string &genName)
	Convert a genTau decay mode string ('oneProng0Pi0') to the RECO enum. More...
std::string	translateRecoDecayModeToGen (PFTau::hadronicDecayMode decayMode)
	Convert a RECO enum decay mode to a string ('oneProng0Pi0') More...

Typedef Documentation

typedef tau::CombinatoricGenerator<ChargedHadronList> reco::tau::ChargedHadronCombo

Definition at line 24 of file RecoTauBuilderCombinatoricPlugin.cc.

typedef std::vector<reco::PFRecoTauChargedHadron> reco::tau::ChargedHadronList

Definition at line 23 of file RecoTauBuilderCombinatoricPlugin.cc.

typedef tau::CombinatoricGenerator<PiZeroList> reco::tau::PiZeroCombo

Definition at line 26 of file RecoTauBuilderCombinatoricPlugin.cc.

typedef std::vector< RecoTauPiZero > reco::tau::PiZeroList

Definition at line 25 of file RecoTauBuilderCombinatoricPlugin.cc.

Enumeration Type Documentation

anonymous enum

Enumerator
kOldDMwoLT
kOldDMwLT
kNewDMwoLT
kNewDMwLT
kDBoldDMwLT
kDBnewDMwLT
kPWoldDMwLT
kPWnewDMwLT
kDBoldDMwLTwGJ
kDBnewDMwLTwGJ

Definition at line 65 of file PFRecoTauClusterVariables.h.

       {kOldDMwoLT, kOldDMwLT, kNewDMwoLT, kNewDMwLT, kDBoldDMwLT, kDBnewDMwLT, kPWoldDMwLT, kPWnewDMwLT, 
        kDBoldDMwLTwGJ, kDBnewDMwLTwGJ};

Function Documentation

template<typename RefVectorType , typename BaseView >

RefVectorType reco::tau::castView

(

const edm::Handle< BaseView > &

view

)

Convert a BaseView (View<T>) to a TRefVector.

Definition at line 74 of file RecoTauCommonUtilities.h.

References mps_fire::i, and convertSQLitetoXML_cfg::output.

                                                        {
  typedef typename RefVectorType::value_type OutputRef;
  // Double check at compile time that the inheritance is okay.  It can still
  // fail at runtime if you pass it the wrong collection.
  BOOST_STATIC_ASSERT(
      (boost::is_base_of<typename BaseView::value_type,
                         typename RefVectorType::member_type>::value));
  RefVectorType output;
  size_t nElements = view->size();
  output.reserve(nElements);
  // Cast each of our Refs
  for (size_t i = 0; i < nElements; ++i) {
    output.push_back(view->refAt(i).template castTo<OutputRef>());
  }
  return output;
}

unsigned int reco::tau::chargedHadronsInDecayMode

(

PFTau::hadronicDecayMode

mode

)

Reverse mapping of decay modes into multiplicities.

Definition at line 45 of file PFTauDecayModeTools.cc.

References reco::PFTau::kOneProngNPiZero, and mode.

                                                                    {
   int modeAsInt = static_cast<int>(mode);
   return (modeAsInt / PFTau::kOneProngNPiZero) + 1;
}

template<typename Container , class OverlapFunction >

Container reco::tau::cleanOverlaps

(

const Container &

dirty

)

Definition at line 35 of file RecoTauCleaningTools.h.

References clean, and analyzePatCleaning_cfg::overlaps.

                                                {
  typedef typename Container::const_iterator Iterator;
  // Output container of clean objects
  Container clean;
  OverlapFunction overlapChecker;
  for (Iterator candidate = dirty.begin(); candidate != dirty.end();
       ++candidate) {
    // Check if this overlaps with a pizero already in the clean list
    bool overlaps = false;
    for (Iterator cleaned = clean.begin();
         cleaned != clean.end() && !overlaps; ++cleaned) {
      overlaps = overlapChecker(*candidate, *cleaned);
    }
    // If it didn't overlap with anything clean, add it to the clean list
    if (!overlaps)
      clean.insert(clean.end(), *candidate);
  }
  return clean;
}

reco::Candidate::LorentzVector reco::tau::compChargedHadronP4fromPThetaPhi	(	double	chargedHadronP,
		double	chargedHadronTheta,
		double	chargedHadronPhi
	)

Definition at line 58 of file pfRecoTauChargedHadronAuxFunctions.cc.

References compChargedHadronP4fromPxPyPz().

{
  double chargedHadronPx = chargedHadronP*TMath::Cos(chargedHadronPhi)*TMath::Sin(chargedHadronTheta);
  double chargedHadronPy = chargedHadronP*TMath::Sin(chargedHadronPhi)*TMath::Sin(chargedHadronTheta);
  double chargedHadronPz = chargedHadronP*TMath::Cos(chargedHadronTheta);    
  return compChargedHadronP4fromPxPyPz(chargedHadronPx, chargedHadronPy, chargedHadronPz);
}

reco::Candidate::LorentzVector reco::tau::compChargedHadronP4fromPxPyPz	(	double	chargedHadronPx,
		double	chargedHadronPy,
		double	chargedHadronPz
	)

Definition at line 50 of file pfRecoTauChargedHadronAuxFunctions.cc.

References mathSSE::sqrt().

Referenced by compChargedHadronP4fromPThetaPhi(), reco::tau::PFRecoTauEnergyAlgorithmPlugin::operator()(), and setChargedHadronP4().

{
  const double chargedPionMass = 0.13957; // GeV  
  double chargedHadronEn = sqrt(chargedHadronPx*chargedHadronPx + chargedHadronPy*chargedHadronPy + chargedHadronPz*chargedHadronPz + chargedPionMass*chargedPionMass);  
  reco::Candidate::LorentzVector chargedHadronP4(chargedHadronPx, chargedHadronPy, chargedHadronPz, chargedHadronEn);
  return chargedHadronP4;
}

std::string reco::tau::discPluginName ( const std::string &  mvaName )

inline

Definition at line 23 of file RecoTauDiscriminantPlugins.h.

References taus_updatedMVAIds_cff::mvaName.

Referenced by PFTauMVAInputDiscriminantTranslator::PFTauMVAInputDiscriminantTranslator().

                                                          {
  return "RecoTauDiscrimination" + mvaName;
}

float reco::tau::eratio

(

const reco::PFTau &

tau

)

return ratio of energy in ECAL over sum of energy in ECAL and HCAL

Definition at line 57 of file PFRecoTauClusterVariables.cc.

References reco::PFTau::signalPFCands(), and pileupDistInMC::total.

Referenced by reco::tau::PFRecoTauDiscriminationByMVAIsolationRun2::discriminate(), and fillIsoMVARun2Inputs().

                                     {
    float ecal_en_in_signal_pf_cands = 0;
    float hcal_en_in_signal_pf_cands = 0;
    for (const auto& signal_cand : tau.signalPFCands()) {
      ecal_en_in_signal_pf_cands += signal_cand->ecalEnergy();
      hcal_en_in_signal_pf_cands += signal_cand->hcalEnergy();
    }
    float total = ecal_en_in_signal_pf_cands + hcal_en_in_signal_pf_cands;
    if(total==0){ 
      return -1;
    }
    return ecal_en_in_signal_pf_cands/total;
  }

float reco::tau::eratio

(

const pat::Tau &

tau

)

Definition at line 70 of file PFRecoTauClusterVariables.cc.

References pat::Tau::ecalEnergy(), pat::Tau::hcalEnergy(), and pileupDistInMC::total.

                                  {
    float ecal_en_in_signal_cands = tau.ecalEnergy();
    float hcal_en_in_signal_cands = tau.hcalEnergy();
    float total = ecal_en_in_signal_cands + hcal_en_in_signal_cands;
    if(total == 0){ 
      return -1;
    }
    return ecal_en_in_signal_cands/total;
  }

std::pair< edm::ParameterSet, edm::ParameterSet > reco::tau::factorizePUQCuts

(

const edm::ParameterSet &

inputSet

)

Definition at line 317 of file RecoTauQualityCuts.cc.

References edm::ParameterSet::copyFrom(), TkAlMuonSelectors_cfi::cut, edm::ParameterSet::getParameterNames(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by PFRecoTauDiscriminationByIsolation::PFRecoTauDiscriminationByIsolation().

{
  edm::ParameterSet puCuts;
  edm::ParameterSet nonPUCuts;

  std::vector<std::string> inputNames = input.getParameterNames();
  BOOST_FOREACH( const std::string& cut, inputNames ) {
    if ( cut == "minTrackVertexWeight" || 
     cut == "maxDeltaZ"            ||
     cut == "maxDeltaZToLeadTrack" ) {
      puCuts.copyFrom(input, cut);
    } else {
      nonPUCuts.copyFrom(input, cut);
    }
  }
  return std::make_pair(puCuts, nonPUCuts);
}

bool reco::tau::fillIsoMVARun2Inputs	(	float *	mvaInput,
		const pat::Tau &	tau,
		int	mvaOpt,
		const std::string &	nameCharged,
		const std::string &	nameNeutral,
		const std::string &	namePu,
		const std::string &	nameOutside,
		const std::string &	nameFootprint
	)

Definition at line 158 of file PFRecoTauClusterVariables.cc.

References funct::abs(), tauProducer_cfi::chargedIsoPtSum, pat::Tau::decayMode(), pat::Tau::dxy(), pat::Tau::dxy_Sig(), MillePedeFileConverter_cfg::e, eratio(), reco::LeafCandidate::eta(), f, pat::Tau::flightLength(), pat::Tau::flightLengthSig(), tauProducer_cfi::footprintCorrection, pat::Tau::hasSecondaryVertex(), pat::Tau::ip3d(), pat::Tau::ip3d_Sig(), kDBnewDMwLT, kDBnewDMwLTwGJ, kDBoldDMwLT, kDBoldDMwLTwGJ, kNewDMwLT, kNewDMwoLT, kOldDMwLT, kOldDMwoLT, kPWnewDMwLT, kPWoldDMwLT, pat::Tau::leadingTrackNormChi2(), cmsBatch::log, SiStripPI::max, min(), n_photons_total(), tauProducer_cfi::neutralIsoPtSum, reco::LeafCandidate::p(), reco::LeafCandidate::p4(), tauProducer_cfi::photonPtSumOutsideSignalCone, funct::pow(), reco::LeafCandidate::pt(), pt_weighted_deta_strip(), pt_weighted_dphi_strip(), pt_weighted_dr_iso(), pt_weighted_dr_signal(), tauProducer_cfi::puCorrPtSum, mathSSE::sqrt(), and pat::Tau::tauID().

Referenced by reco::tau::PATTauDiscriminationByMVAIsolationRun2::discriminate().

  {
    int tauDecayMode = tau.decayMode();
    
  if ( ((mvaOpt == kOldDMwoLT || mvaOpt == kOldDMwLT || mvaOpt == kDBoldDMwLT || mvaOpt == kPWoldDMwLT || mvaOpt ==   kDBoldDMwLTwGJ)
          && (tauDecayMode == 0 || tauDecayMode == 1 || tauDecayMode == 2 || tauDecayMode == 10))
         ||
       ((mvaOpt == kNewDMwoLT || mvaOpt == kNewDMwLT || mvaOpt == kDBnewDMwLT || mvaOpt == kPWnewDMwLT || mvaOpt ==   kDBnewDMwLTwGJ)
        && (tauDecayMode == 0 || tauDecayMode == 1 || tauDecayMode == 2 || tauDecayMode == 5 || tauDecayMode == 6 ||   tauDecayMode == 10 || tauDecayMode == 11))
    ) {
    
      float chargedIsoPtSum = tau.tauID(nameCharged);
      float neutralIsoPtSum = tau.tauID(nameNeutral);
      float puCorrPtSum     = tau.tauID(namePu);
      float photonPtSumOutsideSignalCone = tau.tauID(nameOutside);
      float footprintCorrection = tau.tauID(nameFootprint);
        
      float decayDistX = tau.flightLength().x();
      float decayDistY = tau.flightLength().y();
      float decayDistZ = tau.flightLength().z();
      float decayDistMag = std::sqrt(decayDistX*decayDistX + decayDistY*decayDistY + decayDistZ*decayDistZ);
        
      // --- The following 5 variables differ slightly between AOD & MiniAOD
      //     because they are recomputed using packedCandidates saved in the tau
      float nPhoton = reco::tau::n_photons_total(tau);
      float ptWeightedDetaStrip = reco::tau::pt_weighted_deta_strip(tau, tauDecayMode);
      float ptWeightedDphiStrip = reco::tau::pt_weighted_dphi_strip(tau, tauDecayMode);
      float ptWeightedDrSignal = reco::tau::pt_weighted_dr_signal(tau, tauDecayMode);
      float ptWeightedDrIsolation = reco::tau::pt_weighted_dr_iso(tau, tauDecayMode);
      // ---
      float leadingTrackChi2 = tau.leadingTrackNormChi2();
      float eRatio = reco::tau::eratio(tau);
  
      // Difference between measured and maximally allowed Gottfried-Jackson angle
      float gjAngleDiff = -999;
      if ( tauDecayMode == 10 ) {
          double mTau = 1.77682;
          double mAOne = tau.p4().M();
          double pAOneMag = tau.p();
          double argumentThetaGJmax = (std::pow(mTau,2) - std::pow(mAOne,2) ) / ( 2 * mTau * pAOneMag );
        double argumentThetaGJmeasured = ( tau.p4().px() * decayDistX + tau.p4().py() * decayDistY + tau.p4().pz() * decayDistZ   ) / ( pAOneMag * decayDistMag );
          if ( std::abs(argumentThetaGJmax) <= 1. && std::abs(argumentThetaGJmeasured) <= 1. ) {
              double thetaGJmax = std::asin( argumentThetaGJmax );
              double thetaGJmeasured = std::acos( argumentThetaGJmeasured );
              gjAngleDiff = thetaGJmeasured - thetaGJmax;
          }
      }
        
      if ( mvaOpt == kOldDMwoLT || mvaOpt == kNewDMwoLT ) {
        mvaInput[0]  = std::log(std::max(1.f, (float)tau.pt()));
        mvaInput[1]  = std::abs((float)tau.eta());
        mvaInput[2]  = std::log(std::max(1.e-2f, chargedIsoPtSum));
        mvaInput[3]  = std::log(std::max(1.e-2f, neutralIsoPtSum - 0.125f*puCorrPtSum));
        mvaInput[4]  = std::log(std::max(1.e-2f, puCorrPtSum));
        mvaInput[5]  = tauDecayMode;
      } else if ( mvaOpt == kOldDMwLT || mvaOpt == kNewDMwLT  ) {
        mvaInput[0]  = std::log(std::max(1.f, (float)tau.pt()));
        mvaInput[1]  = std::abs((float)tau.eta());
        mvaInput[2]  = std::log(std::max(1.e-2f, chargedIsoPtSum));
        mvaInput[3]  = std::log(std::max(1.e-2f, neutralIsoPtSum - 0.125f*puCorrPtSum));
        mvaInput[4]  = std::log(std::max(1.e-2f, puCorrPtSum));
        mvaInput[5]  = tauDecayMode;
        mvaInput[6]  = std::copysign(+1.f, tau.dxy());
        mvaInput[7]  = std::sqrt(std::min(1.f, std::abs(tau.dxy())));
        mvaInput[8]  = std::min(10.f, std::abs(tau.dxy_Sig()));
        mvaInput[9]  = ( tau.hasSecondaryVertex() ) ? 1. : 0.;
        mvaInput[10] = std::sqrt(decayDistMag);
        mvaInput[11] = std::min(10.f, tau.flightLengthSig());
      } else if ( mvaOpt == kDBoldDMwLT || mvaOpt == kDBnewDMwLT ) {
        mvaInput[0]  = std::log(std::max(1.f, (float)tau.pt()));
        mvaInput[1]  = std::abs((float)tau.eta());
        mvaInput[2]  = std::log(std::max(1.e-2f, chargedIsoPtSum));
        mvaInput[3]  = std::log(std::max(1.e-2f, neutralIsoPtSum));
        mvaInput[4]  = std::log(std::max(1.e-2f, puCorrPtSum));
        mvaInput[5]  = std::log(std::max(1.e-2f, photonPtSumOutsideSignalCone));
        mvaInput[6]  = tauDecayMode;
        mvaInput[7]  = std::min(30.f, nPhoton);
        mvaInput[8]  = std::min(0.5f, ptWeightedDetaStrip);
        mvaInput[9]  = std::min(0.5f, ptWeightedDphiStrip);
        mvaInput[10] = std::min(0.5f, ptWeightedDrSignal);
        mvaInput[11] = std::min(0.5f, ptWeightedDrIsolation);
        mvaInput[12] = std::min(100.f, leadingTrackChi2);
        mvaInput[13] = std::min(1.f, eRatio);
        mvaInput[14]  = std::copysign(+1.f, tau.dxy());
        mvaInput[15]  = std::sqrt(std::min(1.f, std::abs(tau.dxy())));
        mvaInput[16]  = std::min(10.f, std::abs(tau.dxy_Sig()));
        mvaInput[17]  = std::copysign(+1.f, tau.ip3d());
        mvaInput[18]  = std::sqrt(std::min(1.f, std::abs(tau.ip3d())));
        mvaInput[19]  = std::min(10.f, std::abs(tau.ip3d_Sig()));
        mvaInput[20]  = ( tau.hasSecondaryVertex() ) ? 1. : 0.;
        mvaInput[21] = std::sqrt(decayDistMag);
        mvaInput[22] = std::min(10.f, tau.flightLengthSig());
      } else if ( mvaOpt == kPWoldDMwLT || mvaOpt == kPWnewDMwLT ) {
        mvaInput[0]  = std::log(std::max(1.f, (float)tau.pt()));
        mvaInput[1]  = std::abs((float)tau.eta());
        mvaInput[2]  = std::log(std::max(1.e-2f, chargedIsoPtSum));
        mvaInput[3]  = std::log(std::max(1.e-2f, neutralIsoPtSum));
        mvaInput[4]  = std::log(std::max(1.e-2f, footprintCorrection));
        mvaInput[5]  = std::log(std::max(1.e-2f, photonPtSumOutsideSignalCone));
        mvaInput[6]  = tauDecayMode;
        mvaInput[7]  = std::min(30.f, nPhoton);
        mvaInput[8]  = std::min(0.5f, ptWeightedDetaStrip);
        mvaInput[9]  = std::min(0.5f, ptWeightedDphiStrip);
        mvaInput[10] = std::min(0.5f, ptWeightedDrSignal);
        mvaInput[11] = std::min(0.5f, ptWeightedDrIsolation);
        mvaInput[12] = std::min(100.f, leadingTrackChi2);
        mvaInput[13] = std::min(1.f, eRatio);
        mvaInput[14]  = std::copysign(+1.f, tau.dxy());
        mvaInput[15]  = std::sqrt(std::min(1.f, std::abs(tau.dxy())));
        mvaInput[16]  = std::min(10.f, std::abs(tau.dxy_Sig()));
        mvaInput[17]  = std::copysign(+1.f, tau.ip3d());
        mvaInput[18]  = std::sqrt(std::min(1.f, std::abs(tau.ip3d())));
        mvaInput[19]  = std::min(10.f, std::abs(tau.ip3d_Sig()));
        mvaInput[20]  = ( tau.hasSecondaryVertex() ) ? 1. : 0.;
        mvaInput[21] = std::sqrt(decayDistMag);
        mvaInput[22] = std::min(10.f, tau.flightLengthSig());
      } else if ( mvaOpt == kDBoldDMwLTwGJ || mvaOpt == kDBnewDMwLTwGJ ) {
        mvaInput[0]  = std::log(std::max(1.f, (float)tau.pt()));
        mvaInput[1]  = std::abs((float)tau.eta());
        mvaInput[2]  = std::log(std::max(1.e-2f, chargedIsoPtSum));
        mvaInput[3]  = std::log(std::max(1.e-2f, neutralIsoPtSum));
        mvaInput[4]  = std::log(std::max(1.e-2f, puCorrPtSum));
        mvaInput[5]  = std::log(std::max(1.e-2f, photonPtSumOutsideSignalCone));
        mvaInput[6]  = tauDecayMode;
        mvaInput[7]  = std::min(30.f, nPhoton);
        mvaInput[8]  = std::min(0.5f, ptWeightedDetaStrip);
        mvaInput[9]  = std::min(0.5f, ptWeightedDphiStrip);
        mvaInput[10] = std::min(0.5f, ptWeightedDrSignal);
        mvaInput[11] = std::min(0.5f, ptWeightedDrIsolation);
        mvaInput[12] = std::min(1.f, eRatio);
        mvaInput[13]  = std::copysign(+1.f, tau.dxy());
        mvaInput[14]  = std::sqrt(std::min(1.f, std::abs(tau.dxy())));
        mvaInput[15]  = std::min(10.f, std::abs(tau.dxy_Sig()));
        mvaInput[16]  = std::copysign(+1.f, tau.ip3d());
        mvaInput[17]  = std::sqrt(std::min(1.f, std::abs(tau.ip3d())));
        mvaInput[18]  = std::min(10.f, std::abs(tau.ip3d_Sig()));
        mvaInput[19]  = ( tau.hasSecondaryVertex() ) ? 1. : 0.;
        mvaInput[20] = std::sqrt(decayDistMag);
        mvaInput[21] = std::min(10.f, tau.flightLengthSig());
        mvaInput[22] = std::max(-1.f, gjAngleDiff);
      }
  
      return true;
    }
    return false;
  }

template<typename Iterator >

std::vector<PFCandidatePtr> reco::tau::filterPFCandidates	(	const Iterator &	begin,
		const Iterator &	end,
		int	particleId,
		bool	sort = true
	)

Filter a collection of objects that are convertible to PFCandidatePtrs by PFCandidate ID

Definition at line 40 of file RecoTauCommonUtilities.h.

References end, flattenPiZeros(), metsig::jet, convertSQLitetoXML_cfg::output, pfCandidates(), pfChargedCands(), and pfGammas().

Referenced by pfCandidates(), and pfChargedCands().

                                    {
  std::vector<PFCandidatePtr> output;
  for(Iterator iter = begin; iter != end; ++iter) {
    reco::PFCandidatePtr ptr(*iter);
    if (ptr->particleId() == particleId)
      output.push_back(ptr);
  }
  if (sort) std::sort(output.begin(), output.end(), SortPFCandsDescendingPt());
  return output;
}

std::vector< PFCandidatePtr > reco::tau::flattenPiZeros	(	const std::vector< RecoTauPiZero >::const_iterator &	piZerosBegin,
		const std::vector< RecoTauPiZero >::const_iterator &	piZerosEnd
	)

Flatten a list of pi zeros into a list of there constituent PFCandidates.

Definition at line 16 of file RecoTauCommonUtilities.cc.

References convertSQLitetoXML_cfg::output.

Referenced by filterPFCandidates(), flattenPiZeros(), reco::tau::xclean::CrossCleanPtrs< PtrIter >::initialize(), and reco::tau::RecoTauPhotonFilter::operator()().

                                                                                                                                     {
  std::vector<PFCandidatePtr> output;

  for(std::vector<RecoTauPiZero>::const_iterator piZero = piZerosBegin;
      piZero != piZerosEnd; ++piZero) {
    for(size_t iDaughter = 0; iDaughter < piZero->numberOfDaughters();
        ++iDaughter) {
      output.push_back(PFCandidatePtr(piZero->daughterPtr(iDaughter)));
    }
  }
  return output;
}

std::vector< PFCandidatePtr > reco::tau::flattenPiZeros

(

const std::vector< RecoTauPiZero > &

piZeros

)

Definition at line 30 of file RecoTauCommonUtilities.cc.

References flattenPiZeros().

                                                        {
  return flattenPiZeros(piZeros.begin(), piZeros.end()); 
}

PFTau::hadronicDecayMode reco::tau::getDecayMode

(

const reco::GenJet *

genJet

)

Definition at line 80 of file PFTauDecayModeTools.cc.

References JetMCTagUtils::genTauDecayMode(), reco::PFTau::kNull, and translateGenDecayModeToReco().

Referenced by tautools::RecoTauDecayModeTruthMatchPlugin::operator()().

                                                              {
  if (!genJet)
    return reco::PFTau::kNull;
  return translateGenDecayModeToReco(JetMCTagUtils::genTauDecayMode(*genJet));
}

float reco::tau::lead_track_chi2

(

const reco::PFTau &

tau

)

return chi2 of the leading track ==> deprecated? <==

Definition at line 45 of file PFRecoTauClusterVariables.cc.

References objects.autophobj::float, edm::Ptr< T >::isNonnull(), edm::Ref< C, T, F >::isNonnull(), and reco::PFTau::leadPFChargedHadrCand().

Referenced by reco::tau::PFRecoTauDiscriminationByMVAIsolationRun2::discriminate().

                                              {
    float LeadingTracknormalizedChi2 = 0;
    const reco::PFCandidatePtr& leadingPFCharged = tau.leadPFChargedHadrCand() ;
    if (leadingPFCharged.isNonnull()) {
      reco::TrackRef tref = leadingPFCharged -> trackRef();
      if (tref.isNonnull()) {
        LeadingTracknormalizedChi2 = (float)(tref -> normalizedChi2());
      }
    }
    return LeadingTracknormalizedChi2;
  }

template<typename InputIterator >

InputIterator reco::tau::leadPFCand	(	InputIterator	begin,
		InputIterator	end
	)

Definition at line 132 of file RecoTauCommonUtilities.h.

References begin, and end.

Referenced by reco::tau::RecoTauConstructor::get(), reco::tau::RecoTauVertexAssociator::getLeadTrack(), pat::Tau::leadPFChargedHadrCandsignedSipt(), reco::tau::RecoTauBuilderConePlugin::operator()(), reco::operator<<(), and reco::PFTau::~PFTau().

                       {
    double max_pt = 0;
    InputIterator max_cand = begin;
    for(InputIterator cand = begin; cand != end; ++cand) {
      if( (*cand)->pt() > max_pt ) {
        max_pt = (*cand)->pt();
        max_cand = cand;
      }
    }
    return max_cand;
  }

unsigned int reco::tau::n_photons_total

(

const reco::PFTau &

tau

)

return total number of pf photon candidates with pT>500 MeV, which are associated to signal

Definition at line 133 of file PFRecoTauClusterVariables.cc.

References reco::PFTau::isolationPFGammaCands(), and reco::PFTau::signalPFGammaCands().

Referenced by reco::tau::PFRecoTauDiscriminationByMVAIsolationRun2::discriminate(), fillIsoMVARun2Inputs(), and pt_weighted_dr_iso().

                                                     {
    unsigned int n_photons = 0;
    for (auto& cand : tau.signalPFGammaCands()) {
      if (cand->pt() > 0.5)
        ++n_photons;
    }
    for (auto& cand : tau.isolationPFGammaCands()) {
      if (cand->pt() > 0.5)
        ++n_photons;
    }
    return n_photons;
  }

unsigned int reco::tau::n_photons_total

(

const pat::Tau &

tau

)

Definition at line 145 of file PFRecoTauClusterVariables.cc.

References pat::Tau::isolationGammaCands(), and pat::Tau::signalGammaCands().

                                                  {
    unsigned int n_photons = 0;
    for (auto& cand : tau.signalGammaCands()) {
      if (cand->pt() > 0.5) 
        ++n_photons;
    }
    for (auto& cand : tau.isolationGammaCands()) {  
      if (cand->pt() > 0.5)
    ++n_photons;
    }
    return n_photons;
  }  

std::vector< reco::PFCandidatePtr > reco::tau::pfCandidates	(	const PFJet &	jet,
		int	particleId,
		bool	sort = true
	)

Extract pfCandidates of a given particle Id from a PFJet. If sort is true, candidates will be sorted by descending PT

Definition at line 34 of file RecoTauCommonUtilities.cc.

References filterPFCandidates(), and reco::PFJet::getPFConstituents().

Referenced by filterPFCandidates(), reco::tau::RecoTauBuilderCombinatoricPlugin::operator()(), reco::tau::RecoTauBuilderConePlugin::operator()(), reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::operator()(), reco::tau::RecoTauPiZeroStripPlugin::operator()(), reco::tau::RecoTauPiZeroStripPlugin2::operator()(), reco::tau::RecoTauPiZeroStripPlugin3::operator()(), and pfGammas().

                               {
  PFCandPtrs pfCands = jet.getPFConstituents();
  PFCandPtrs selectedPFCands = filterPFCandidates(
      pfCands.begin(), pfCands.end(), particleId, sort);
  return selectedPFCands;
}

std::vector< reco::PFCandidatePtr > reco::tau::pfCandidates	(	const PFJet &	jet,
		const std::vector< int > &	particleIds,
		bool	sort = true
	)

Extract pfCandidates of a that match a list of particle Ids from a PFJet.

Definition at line 42 of file RecoTauCommonUtilities.cc.

References filterPFCandidates(), reco::PFJet::getPFConstituents(), and convertSQLitetoXML_cfg::output.

                                                  {
  PFCandPtrs&& pfCands = jet.getPFConstituents();
  PFCandPtrs output;
  // Get each desired candidate type, unsorted for now
  for(std::vector<int>::const_iterator particleId = particleIds.begin();
      particleId != particleIds.end(); ++particleId) {
    PFCandPtrs&& selectedPFCands = filterPFCandidates(pfCands.begin(), pfCands.end(), *particleId, false);
    output.insert(output.end(), selectedPFCands.begin(), selectedPFCands.end());
  }
  if (sort) std::sort(output.begin(), output.end(), SortPFCandsDescendingPt());
  return output;
}

std::vector< reco::PFCandidatePtr > reco::tau::pfChargedCands	(	const PFJet &	jet,
		bool	sort = true
	)

Extract all non-neutral candidates from a PFJet.

Definition at line 60 of file RecoTauCommonUtilities.cc.

References reco::PFCandidate::e, filterPFCandidates(), reco::PFJet::getPFConstituents(), reco::PFCandidate::h, reco::PFCandidate::mu, and convertSQLitetoXML_cfg::output.

Referenced by filterPFCandidates(), reco::tau::RecoTauVertexAssociator::getLeadTrack(), reco::tau::RecoTauBuilderCombinatoricPlugin::operator()(), and reco::tau::RecoTauBuilderConePlugin::operator()().

                                                            {
  PFCandPtrs&& pfCands = jet.getPFConstituents();
  PFCandPtrs output;
  PFCandPtrs&& mus = filterPFCandidates(pfCands.begin(), pfCands.end(), reco::PFCandidate::mu, false);
  PFCandPtrs&& es = filterPFCandidates(pfCands.begin(), pfCands.end(), reco::PFCandidate::e, false);
  PFCandPtrs&& chs = filterPFCandidates(pfCands.begin(), pfCands.end(), reco::PFCandidate::h, false);
  output.reserve(mus.size() + es.size() + chs.size());
  output.insert(output.end(), mus.begin(), mus.end());
  output.insert(output.end(), es.begin(), es.end());
  output.insert(output.end(), chs.begin(), chs.end());
  if (sort) std::sort(output.begin(), output.end(), SortPFCandsDescendingPt());
  return output;
}

std::vector< reco::PFCandidatePtr > reco::tau::pfGammas	(	const PFJet &	jet,
		bool	sort = true
	)

Extract all pfGammas from a PFJet.

Definition at line 56 of file RecoTauCommonUtilities.cc.

References reco::PFCandidate::gamma, and pfCandidates().

Referenced by filterPFCandidates(), reco::tau::RecoTauPiZeroTrivialPlugin::operator()(), reco::tau::RecoTauPiZeroCombinatoricPlugin::operator()(), reco::tau::RecoTauBuilderConePlugin::operator()(), and reco::tau::RecoTauConstructor::reserveTauChargedHadron().

                                                                        {
  return pfCandidates(jet, reco::PFCandidate::gamma, sort);
}

unsigned int reco::tau::piZerosInDecayMode

(

PFTau::hadronicDecayMode

mode

)

Definition at line 50 of file PFTauDecayModeTools.cc.

References reco::PFTau::kOneProngNPiZero, and mode.

                                                             {
   int modeAsInt = static_cast<int>(mode);
   return (modeAsInt % PFTau::kOneProngNPiZero);
}

float reco::tau::pt_weighted_deta_strip ( const reco::PFTau &  tau, int  dm  )

inline

return sum of pt weighted values of deta relative to tau candidate for all pf photon candidates, which are associated to signal

Definition at line 39 of file PFRecoTauClusterVariables.h.

References pt_weighted_dx().

Referenced by reco::tau::PFRecoTauDiscriminationByMVAIsolationRun2::discriminate(), and fillIsoMVARun2Inputs().

                                                                    {
    return pt_weighted_dx(tau, dm==10 ? 2 : 1, 1, dm);
  }

float reco::tau::pt_weighted_deta_strip ( const pat::Tau &  tau, int  dm  )

inline

Definition at line 42 of file PFRecoTauClusterVariables.h.

References pt_weighted_dx().

                                                                 {
    return pt_weighted_dx(tau, dm==10 ? 2 : 1, 1, dm);
  }

float reco::tau::pt_weighted_dphi_strip ( const reco::PFTau &  tau, int  dm  )

inline

return sum of pt weighted values of dphi relative to tau candidate for all pf photon candidates, which are associated to signal

Definition at line 47 of file PFRecoTauClusterVariables.h.

References pt_weighted_dx().

Referenced by reco::tau::PFRecoTauDiscriminationByMVAIsolationRun2::discriminate(), and fillIsoMVARun2Inputs().

                                                                    {
    return pt_weighted_dx(tau, dm==10 ? 2 : 1, 2, dm);
  }

float reco::tau::pt_weighted_dphi_strip ( const pat::Tau &  tau, int  dm  )

inline

Definition at line 50 of file PFRecoTauClusterVariables.h.

References pt_weighted_dx().

                                                                 {
    return pt_weighted_dx(tau, dm==10 ? 2 : 1, 2, dm);
  }  

float reco::tau::pt_weighted_dr_iso ( const reco::PFTau &  tau, int  dm  )

inline

return sum of pt weighted values of dr relative to tau candidate for all pf photon candidates, which are inside an isolation conde but not associated to signal

Definition at line 55 of file PFRecoTauClusterVariables.h.

References pt_weighted_dx().

Referenced by reco::tau::PFRecoTauDiscriminationByMVAIsolationRun2::discriminate(), and fillIsoMVARun2Inputs().

                                                                {
    return pt_weighted_dx(tau, 2, 0, dm);
  }

float reco::tau::pt_weighted_dr_iso ( const pat::Tau &  tau, int  dm  )

inline

Definition at line 58 of file PFRecoTauClusterVariables.h.

References n_photons_total(), pt_weighted_dx(), and metsig::tau.

                                                             {
    return pt_weighted_dx(tau, 2, 0, dm);
  }

float reco::tau::pt_weighted_dr_signal ( const reco::PFTau &  tau, int  dm  )

inline

return sum of pt weighted values of dr relative to tau candidate for all pf photon candidates, which are associated to signal

Definition at line 31 of file PFRecoTauClusterVariables.h.

References pt_weighted_dx().

Referenced by reco::tau::PFRecoTauDiscriminationByMVAIsolationRun2::discriminate(), and fillIsoMVARun2Inputs().

                                                                   {
    return pt_weighted_dx(tau, 0, 0, dm);
  }

float reco::tau::pt_weighted_dr_signal ( const pat::Tau &  tau, int  dm  )

inline

Definition at line 34 of file PFRecoTauClusterVariables.h.

References pt_weighted_dx().

                                                                {
    return pt_weighted_dx(tau, 0, 0, dm);
  }

float reco::tau::pt_weighted_dx	(	const reco::PFTau &	tau,
		int	mode = 0,
		int	var = 0,
		int	decaymode = -1
	)

return sum of pt weighted values of distance to tau candidate for all pf photon candidates, which are associated to signal; depending on var the distance is in 0=:dr, 1=:deta, 2=:dphi

Definition at line 126 of file PFRecoTauClusterVariables.cc.

References runTauDisplay::decaymode, ALCARECOPromptCalibProdSiPixelAli0T_cff::mode, metsig::tau, and JetChargeProducer_cfi::var.

Referenced by pt_weighted_deta_strip(), pt_weighted_dphi_strip(), pt_weighted_dr_iso(), and pt_weighted_dr_signal().

                                                                               {
    return pt_weighted_dx_T<PFTau_traits>(tau, mode, var, decaymode);
  }

float reco::tau::pt_weighted_dx	(	const pat::Tau &	tau,
		int	mode = 0,
		int	var = 0,
		int	decaymode = -1
	)

Definition at line 129 of file PFRecoTauClusterVariables.cc.

References runTauDisplay::decaymode, ALCARECOPromptCalibProdSiPixelAli0T_cff::mode, metsig::tau, and JetChargeProducer_cfi::var.

                                                                            {
    return pt_weighted_dx_T<PATTau_traits>(tau, mode, var, decaymode);
  }

template<typename T >

float reco::tau::pt_weighted_dx_T	(	const typename T::Tau_t &	tau,
		int	mode,
		int	var,
		int	decaymode
	)

return sum of pt weighted values of distance to tau candidate for all pf photon candidates, which are associated to signal; depending on var the distance is in 0=:dr, 1=:deta, 2=:dphi

Definition at line 82 of file PFRecoTauClusterVariables.cc.

References funct::abs(), egammaForCoreTracking_cff::cands, reco::deltaPhi(), reco::deltaR(), runTauDisplay::dr, RemoveAddSevLevel::flag, SiStripPI::max, min(), EnergyCorrector::pt, and metsig::tau.

                                                                                       {
    float sum_pt = 0.;
    float sum_dx_pt = 0.;
    float signalrad = std::max(0.05, std::min(0.1, 3./std::max(1., tau.pt())));
    int is3prong = (decaymode==10);
    const auto& cands = getGammas_T<T>(tau, mode < 2);
    for (const auto& cand : cands) {
      // only look at electrons/photons with pT > 0.5
      if (cand->pt() < 0.5){
        continue;
      }
      float dr = reco::deltaR(cand->eta(), cand->phi(), tau.eta(), tau.phi());
      float deta = std::abs(cand->eta() - tau.eta());
      float dphi = std::abs(reco::deltaPhi(cand->phi(), tau.phi()));
      float pt = cand->pt();
      bool flag = isInside(pt, deta, dphi);
      if(is3prong==0){
        if (mode == 2 || (mode == 0 && dr < signalrad) || (mode == 1 && dr > signalrad)) {
          sum_pt += pt;
          if (var == 0)
            sum_dx_pt += pt * dr;
          else if (var == 1)
            sum_dx_pt += pt * deta;
          else if (var == 2)
            sum_dx_pt += pt * dphi;
        }
      }
      else if(is3prong==1){
        if( (mode==2 && flag==false) || (mode==1 && flag==true) || mode==0){
          sum_pt += pt;
          if (var == 0)
            sum_dx_pt += pt * dr;
          else if (var == 1)
            sum_dx_pt += pt * deta;
          else if (var == 2)
            sum_dx_pt += pt * dphi;
        }
      }
    }
    if (sum_pt > 0.){
      return sum_dx_pt/sum_pt;
    }
    return 0.;
  }

static const edm::ParameterSetDescriptionFillerPluginFactory::PMaker<edm::ParameterSetDescriptionFiller< PATTauDiscriminationByMVAIsolationRun2 > > reco::tau::s_filler__LINE__ ( "PATTauDiscriminationByMVAIsolationRun2"  )

static

static const edm::ParameterSetDescriptionFillerPluginFactory::PMaker<edm::ParameterSetDescriptionFiller< PFRecoTauDiscriminationByMVAIsolationRun2 > > reco::tau::s_filler__LINE__ ( "PFRecoTauDiscriminationByMVAIsolationRun2"  )

static

static const edm::MakerPluginFactory ::PMaker< edm::WorkerMaker< PATTauDiscriminationByMVAIsolationRun2 > > reco::tau::s_maker__LINE__ ( "PATTauDiscriminationByMVAIsolationRun2"  )

static

static const edm::MakerPluginFactory ::PMaker< edm::WorkerMaker< PFRecoTauDiscriminationByMVAIsolationRun2 > > reco::tau::s_maker__LINE__ ( "PFRecoTauDiscriminationByMVAIsolationRun2"  )

static

void reco::tau::setChargedHadronP4	(	reco::PFRecoTauChargedHadron &	chargedHadron,
		double	scaleFactor_neutralPFCands = 1.0
	)

Definition at line 10 of file pfRecoTauChargedHadronAuxFunctions.cc.

References reco::PFRecoTauChargedHadron::algoIs(), compChargedHadronP4fromPxPyPz(), reco::PFRecoTauChargedHadron::getChargedPFCandidate(), reco::PFRecoTauChargedHadron::getNeutralPFCandidates(), reco::PFRecoTauChargedHadron::getTrack(), edm::Ptr< T >::isNonnull(), reco::PFRecoTauChargedHadron::kChargedPFCandidate, reco::PFRecoTauChargedHadron::kPFNeutralHadron, reco::PFRecoTauChargedHadron::kTrack, electrons_cff::ptRatio, reco::LeafCandidate::setP4(), and HiIsolationCommonParameters_cff::track.

Referenced by reco::tau::RecoTauConstructor::addTauChargedHadron(), reco::tau::PFRecoTauEnergyAlgorithmPlugin::operator()(), reco::tau::PFRecoTauChargedHadronFromPFCandidatePlugin::operator()(), reco::tau::PFRecoTauChargedHadronFromTrackPlugin::operator()(), and PFRecoTauChargedHadronProducer::produce().

{
  double chargedHadronP  = 0.;
  double chargedHadronPx = 0.;
  double chargedHadronPy = 0.;
  double chargedHadronPz = 0.;
  double SumNeutrals = 0.;
  if ( chargedHadron.algoIs(reco::PFRecoTauChargedHadron::kChargedPFCandidate) ||
       chargedHadron.algoIs(reco::PFRecoTauChargedHadron::kPFNeutralHadron)   ) {
    const reco::PFCandidatePtr& chargedPFCand = chargedHadron.getChargedPFCandidate();
    assert(chargedPFCand.isNonnull());
    chargedHadronP     += chargedPFCand->p();
    chargedHadronPx     = chargedPFCand->px();
    chargedHadronPy     = chargedPFCand->py();
    chargedHadronPz     = chargedPFCand->pz();
  } else if ( chargedHadron.algoIs(reco::PFRecoTauChargedHadron::kTrack) ) {
    const reco::PFRecoTauChargedHadron::TrackPtr& track = chargedHadron.getTrack();
    assert(track.isNonnull());
    chargedHadronP     += track->p();
    chargedHadronPx     = track->px();
    chargedHadronPy     = track->py();
    chargedHadronPz     = track->pz();
  } else assert(0);
  const std::vector<reco::PFCandidatePtr>& neutralPFCands = chargedHadron.getNeutralPFCandidates();
  for ( std::vector<reco::PFCandidatePtr>::const_iterator neutralPFCand = neutralPFCands.begin();
    neutralPFCand != neutralPFCands.end(); ++neutralPFCand ) {
    SumNeutrals += (*neutralPFCand)->p();
  }
  double noNeutrals=chargedHadronP;
  chargedHadronP+=scaleFactor_neutralPFCands*SumNeutrals;
  double ptRatio=chargedHadronP/noNeutrals;
  chargedHadronPx*=ptRatio;
  chargedHadronPy*=ptRatio;
  chargedHadronPz*=ptRatio;
  
      
  reco::Candidate::LorentzVector chargedHadronP4 = compChargedHadronP4fromPxPyPz(chargedHadronPx, chargedHadronPy, chargedHadronPz);
  chargedHadron.setP4(chargedHadronP4);
}

template<typename InputIterator >

int reco::tau::sumPFCandCharge	(	InputIterator	begin,
		InputIterator	end
	)

Sum the charge of a collection of PFCandidates.

Definition at line 127 of file RecoTauCommonUtilities.h.

References reco::LeafCandidate::charge(), and sumPFVector().

                       {
    return sumPFVector(begin, end, &PFCandidate::charge, 0);
  }

template<typename InputIterator >

reco::Candidate::LorentzVector reco::tau::sumPFCandP4	(	InputIterator	begin,
		InputIterator	end
	)

Definition at line 114 of file RecoTauCommonUtilities.h.

References reco::LeafCandidate::p4(), and sumPFVector().

                                            {
  return sumPFVector(begin, end, &PFCandidate::p4,
      reco::Candidate::LorentzVector());
}

template<typename InputIterator >

double reco::tau::sumPFCandPt	(	InputIterator	begin,
		InputIterator	end
	)

Sum the pT of a collection of PFCandidates.

Definition at line 121 of file RecoTauCommonUtilities.h.

References reco::LeafCandidate::pt(), and sumPFVector().

Referenced by reco::tau::RecoTauConstructor::get().

                       {
    return sumPFVector(begin, end, &PFCandidate::pt, 0.0);
  }

template<typename InputIterator , typename FunctionPtr , typename ReturnType >

ReturnType reco::tau::sumPFVector	(	InputIterator	begin,
		InputIterator	end,
		FunctionPtr	func,
		ReturnType	init
	)

Sum the four vectors in a collection of PFCandidates.

Definition at line 104 of file RecoTauCommonUtilities.h.

References end, patCandidates_cff::func, init, and convertSQLitetoXML_cfg::output.

Referenced by sumPFCandCharge(), sumPFCandP4(), and sumPFCandPt().

                                         {
    ReturnType output = init;
    for(InputIterator cand = begin; cand != end; ++cand) {
      //#define CALL_MEMBER_FN(object,ptrToMember)  ((object).*(ptrToMember))
      output += ((**cand).*(func))();
    }
    return output;
  }

template<typename InputIterator >

InputIterator reco::tau::takeNElements	(	const InputIterator &	begin,
		const InputIterator &	end,
		size_t	N
	)

Definition at line 96 of file RecoTauCommonUtilities.h.

References begin, and N.

Referenced by reco::tau::RecoTauPiZeroCombinatoricPlugin::operator()(), reco::tau::RecoTauBuilderCombinatoricPlugin::operator()(), reco::tau::RecoTauPiZeroStripPlugin::operator()(), reco::tau::RecoTauPiZeroStripPlugin2::operator()(), and reco::tau::RecoTauPiZeroStripPlugin3::operator()().

                                                                    {
  size_t input_size = end - begin;
  return (N > input_size) ? end : begin + N;
}

PFTau::hadronicDecayMode reco::tau::translateDecayMode	(	unsigned int	nCharged,
		unsigned int	nPiZero
	)

Definition at line 55 of file PFTauDecayModeTools.cc.

References reco::PFTau::kNull, reco::PFTau::kOneProngNPiZero, reco::PFTau::kRareDecayMode, RecoTauCombinatoricProducer_cfi::maxPiZeros, and hpstanc_transforms::nPiZeros.

Referenced by reco::tau::RecoTauObjectEmbedder< T >::operator()(), reco::tau::RecoTauDiscriminantInvariantWidth::RecoTauDiscriminantInvariantWidth(), and RecoTauMVATransform::RecoTauMVATransform().

                                                  {
   // If no tracks exist, this is definitely not a tau!
   if(!nCharged) return PFTau::kNull;
   // Find the maximum number of PiZeros our parameterization can hold
   const unsigned int maxPiZeros = PFTau::kOneProngNPiZero;
   // Determine our track index
   unsigned int trackIndex = (nCharged-1)*(maxPiZeros+1);
   // Check if we handle the given number of tracks
   if(trackIndex >= PFTau::kRareDecayMode) return PFTau::kRareDecayMode;

   nPiZeros = (nPiZeros <= maxPiZeros) ? nPiZeros : maxPiZeros;
   return static_cast<PFTau::hadronicDecayMode>(trackIndex + nPiZeros);
}

PFTau::hadronicDecayMode reco::tau::translateGenDecayModeToReco

(

const std::string &

genName

)

Convert a genTau decay mode string ('oneProng0Pi0') to the RECO enum.

Definition at line 70 of file PFTauDecayModeTools.cc.

Referenced by getDecayMode().

                           {
  return decayModeStringToId(name.c_str(), dmTranslatorMap);
}

std::string reco::tau::translateRecoDecayModeToGen

(

PFTau::hadronicDecayMode

decayMode

)

Convert a RECO enum decay mode to a string ('oneProng0Pi0')

Definition at line 75 of file PFTauDecayModeTools.cc.

                                                                        {
  return decayModeIdToString(decayMode, dmTranslatorMap);
}