reco::tau::qcuts Namespace Reference

Functions
bool	minPackedCandVertexWeight (const pat::PackedCandidate &pCand, const reco::VertexRef *pv, double cut)

Function Documentation

bool reco::tau::qcuts::minPackedCandVertexWeight	(	const pat::PackedCandidate &	pCand,
		const reco::VertexRef *	pv,
		double	cut
	)

Definition at line 55 of file RecoTauQualityCuts.cc.

References pat::PackedCandidate::eta(), edm::Ref< C, T, F >::isNonnull(), edm::Ref< C, T, F >::isNull(), edm::Ref< C, T, F >::key(), LogDebug, pat::PackedCandidate::phi(), pat::PackedCandidate::pt(), pat::PackedCandidate::pvAssociationQuality(), quality, pat::PackedCandidate::UsedInFitLoose, pat::PackedCandidate::UsedInFitTight, pat::PackedCandidate::vertexRef(), and histoStyle::weight.

Referenced by reco::tau::RecoTauQualityCuts::filterChargedCand().

                                                                                                         {
      if (pv->isNull()) {
        edm::LogError("QCutsNoPrimaryVertex") << "Primary vertex Ref in "
                                              << "RecoTauQualityCuts is invalid. - minPackedCandVertexWeight";
        return false;
      }
      //there is some low granular information on track weight in the vertex available with packed cands
      double weight = -9.9;
      if (pCand.vertexRef().isNonnull() && pCand.vertexRef().key() == pv->key()) {
        int quality = pCand.pvAssociationQuality();
        if (quality == pat::PackedCandidate::UsedInFitTight)
          weight = 0.6;  //0.6 as proxy for weight above 0.5
        else if (quality == pat::PackedCandidate::UsedInFitLoose)
          weight = 0.1;  //0.6 as proxy for weight below 0.5
      }
      LogDebug("TauQCuts") << " packedCand: Pt = " << pCand.pt() << ", eta = " << pCand.eta()
                           << ", phi = " << pCand.phi();
      LogDebug("TauQCuts") << " vertex: x = " << (*pv)->position().x() << ", y = " << (*pv)->position().y()
                           << ", z = " << (*pv)->position().z();
      LogDebug("TauQCuts") << "--> trackWeight from packedCand = " << weight << " (cut = " << cut << ")";
      return (weight >= cut);
    }