reco::VertexFilter Class Reference

#include <VertexFilter.h>

Public Member Functions
bool	operator() (const reco::Vertex &pv, const TemplatedSecondaryVertex< reco::Vertex > &sv, const GlobalVector &direction) const
bool	operator() (const reco::Vertex &pv, const TemplatedSecondaryVertex< reco::VertexCompositePtrCandidate > &sv, const GlobalVector &direction) const
	VertexFilter (const edm::ParameterSet &params)
	~VertexFilter ()

Private Attributes
double	distSig2dMax
double	distSig2dMin
double	distSig3dMax
double	distSig3dMin
double	distVal2dMax
double	distVal2dMin
double	distVal3dMax
double	distVal3dMin
double	fracPV
double	massMax
double	maxDeltaRToJetAxis
double	minTrackWeight
unsigned int	multiplicityMin
bool	useTrackWeights
V0Filter	v0Filter

Detailed Description

Definition at line 25 of file VertexFilter.h.

Constructor & Destructor Documentation

VertexFilter::VertexFilter

(

const edm::ParameterSet &

params

)

Definition at line 22 of file VertexFilter.cc.

                                                        :
    useTrackWeights(params.getParameter<bool>("useTrackWeights")),
    minTrackWeight(params.getParameter<double>("minimumTrackWeight")),
    massMax(params.getParameter<double>("massMax")),
    fracPV(params.getParameter<double>("fracPV")),
    multiplicityMin(params.getParameter<unsigned int>("multiplicityMin")),
    distVal2dMin(params.getParameter<double>("distVal2dMin")),
    distVal2dMax(params.getParameter<double>("distVal2dMax")),
    distVal3dMin(params.getParameter<double>("distVal3dMin")),
    distVal3dMax(params.getParameter<double>("distVal3dMax")),
    distSig2dMin(params.getParameter<double>("distSig2dMin")),
    distSig2dMax(params.getParameter<double>("distSig2dMax")),
    distSig3dMin(params.getParameter<double>("distSig3dMin")),
    distSig3dMax(params.getParameter<double>("distSig3dMax")),
    maxDeltaRToJetAxis(params.getParameter<double>("maxDeltaRToJetAxis")),
    v0Filter(params.getParameter<edm::ParameterSet>("v0Filter"))
{
}

reco::VertexFilter::~VertexFilter ( )

inline

Definition at line 28 of file VertexFilter.h.

{}

Member Function Documentation

bool VertexFilter::operator()	(	const reco::Vertex &	pv,
		const TemplatedSecondaryVertex< reco::Vertex > &	sv,
		const GlobalVector &	direction
	)		const

Definition at line 41 of file VertexFilter.cc.

References funct::abs(), vertexTools::computeSharedTracks(), deltaR(), reco::TemplatedSecondaryVertex< SV >::dist2d(), reco::TemplatedSecondaryVertex< SV >::dist3d(), distSig2dMax, distSig2dMin, distSig3dMax, distSig3dMin, distVal2dMax, distVal2dMin, distVal3dMax, distVal3dMin, Measurement1D::error(), fracPV, massMax, maxDeltaRToJetAxis, minTrackWeight, multiplicityMin, reco::Vertex::position(), Measurement1D::significance(), useTrackWeights, v0Filter, Measurement1D::value(), reco::TrackKinematics::vectorSum(), and reco::TrackKinematics::weightedVectorSum().

{
    std::vector<TrackRef> svTracks;
    for(std::vector<TrackBaseRef>::const_iterator iter = sv.tracks_begin();
        iter != sv.tracks_end(); iter++)
        if (sv.trackWeight(*iter) >= minTrackWeight)
            svTracks.push_back(iter->castTo<TrackRef>());

    // minimum number of tracks at vertex

    if (svTracks.size() < multiplicityMin)
        return false;

    // invalid errors

    if (sv.dist2d().error() < 0 || sv.dist3d().error() < 0)
        return false;

    // flight distance limits (value and significance, 2d and 3d)

    if (sv.dist2d().value()        < distVal2dMin ||
        sv.dist2d().value()        > distVal2dMax ||
        sv.dist3d().value()        < distVal3dMin ||
        sv.dist3d().value()        > distVal3dMax ||
        sv.dist2d().significance() < distSig2dMin ||
        sv.dist2d().significance() > distSig2dMax ||
        sv.dist3d().significance() < distSig3dMin ||
        sv.dist3d().significance() > distSig3dMax)
        return false;

    // SV direction filter

    if (Geom::deltaR(sv.position() - pv.position(),
                     (maxDeltaRToJetAxis > 0) ? direction : -direction)
                        > std::abs(maxDeltaRToJetAxis))
        return false;

    // compute fourvector sum of tracks as vertex and cut on inv. mass

    TrackKinematics kin(sv);

    double mass = useTrackWeights ? kin.weightedVectorSum().M()
                                  : kin.vectorSum().M();

    if (mass > massMax)
        return false;

    // find shared tracks between PV and SV

    if (fracPV < 1.0) {
        double fractionSharedTracks =
            vertexTools::computeSharedTracks(pv, svTracks, minTrackWeight);
        if (fractionSharedTracks > fracPV)
            return false;
    }

    // check for V0 vertex

    if (sv.hasRefittedTracks())
        return v0Filter(sv.refittedTracks());
    else
        return v0Filter(svTracks);
}

bool VertexFilter::operator()	(	const reco::Vertex &	pv,
		const TemplatedSecondaryVertex< reco::VertexCompositePtrCandidate > &	sv,
		const GlobalVector &	direction
	)		const

Definition at line 106 of file VertexFilter.cc.

References funct::abs(), vertexTools::computeSharedTracks(), deltaR(), reco::TemplatedSecondaryVertex< SV >::dist2d(), reco::TemplatedSecondaryVertex< SV >::dist3d(), distSig2dMax, distSig2dMin, distSig3dMax, distSig3dMin, distVal2dMax, distVal2dMin, distVal3dMax, distVal3dMin, Measurement1D::error(), fracPV, massMax, maxDeltaRToJetAxis, minTrackWeight, multiplicityMin, reco::Vertex::position(), Measurement1D::significance(), useTrackWeights, v0Filter, Measurement1D::value(), reco::TrackKinematics::vectorSum(), and reco::TrackKinematics::weightedVectorSum().

                                                                 {

    const std::vector<CandidatePtr> & svTracks = sv.daughterPtrVector();

    // minimum number of tracks at vertex

    if (svTracks.size() < multiplicityMin)
        return false;

    // invalid errors

    if (sv.dist2d().error() < 0 || sv.dist3d().error() < 0)
        return false;

    // flight distance limits (value and significance, 2d and 3d)

    if (sv.dist2d().value()        < distVal2dMin ||
            sv.dist2d().value()        > distVal2dMax ||
            sv.dist3d().value()        < distVal3dMin ||
            sv.dist3d().value()        > distVal3dMax ||
            sv.dist2d().significance() < distSig2dMin ||
            sv.dist2d().significance() > distSig2dMax ||
            sv.dist3d().significance() < distSig3dMin ||
            sv.dist3d().significance() > distSig3dMax)
        return false;

    // SV direction filter

    if (Geom::deltaR(sv.vertex() - pv.position(),
                (maxDeltaRToJetAxis > 0) ? direction : -direction)
            > std::abs(maxDeltaRToJetAxis))
        return false;
    // compute fourvector sum of tracks as vertex and cut on inv. mass

    TrackKinematics kin(sv);

    double mass = useTrackWeights ? kin.weightedVectorSum().M()
        : kin.vectorSum().M();

    if (mass > massMax)
        return false;

    // find shared tracks between PV and SV

    if (fracPV < 1.0) {
        double fractionSharedTracks =
            vertexTools::computeSharedTracks(pv, svTracks, minTrackWeight);
        if (fractionSharedTracks  > fracPV)
            return false;
    }

    // check for V0 vertex

    return v0Filter(svTracks);
}

Member Data Documentation

double reco::VertexFilter::distSig2dMax

private

Definition at line 48 of file VertexFilter.h.

Referenced by operator()().

double reco::VertexFilter::distSig2dMin

private

Definition at line 47 of file VertexFilter.h.

Referenced by operator()().

double reco::VertexFilter::distSig3dMax

private

Definition at line 50 of file VertexFilter.h.

Referenced by operator()().

double reco::VertexFilter::distSig3dMin

private

Definition at line 49 of file VertexFilter.h.

Referenced by operator()().

double reco::VertexFilter::distVal2dMax

private

Definition at line 43 of file VertexFilter.h.

Referenced by operator()().

double reco::VertexFilter::distVal2dMin

private

Definition at line 42 of file VertexFilter.h.

Referenced by operator()().

double reco::VertexFilter::distVal3dMax

private

Definition at line 45 of file VertexFilter.h.

Referenced by operator()().

double reco::VertexFilter::distVal3dMin

private

Definition at line 44 of file VertexFilter.h.

Referenced by operator()().

double reco::VertexFilter::fracPV

private

Definition at line 39 of file VertexFilter.h.

Referenced by operator()().

double reco::VertexFilter::massMax

private

Definition at line 38 of file VertexFilter.h.

Referenced by operator()().

double reco::VertexFilter::maxDeltaRToJetAxis

private

Definition at line 52 of file VertexFilter.h.

Referenced by operator()().

double reco::VertexFilter::minTrackWeight

private

Definition at line 37 of file VertexFilter.h.

Referenced by operator()().

unsigned int reco::VertexFilter::multiplicityMin

private

Definition at line 40 of file VertexFilter.h.

Referenced by operator()().

bool reco::VertexFilter::useTrackWeights

private

Definition at line 36 of file VertexFilter.h.

Referenced by operator()().

V0Filter reco::VertexFilter::v0Filter

private

Definition at line 53 of file VertexFilter.h.

Referenced by operator()().