TemplatedInclusiveVertexFinder< InputContainer, VTX > Class Template Reference

#include <InclusiveVertexFinder.h>

Inheritance diagram for TemplatedInclusiveVertexFinder< InputContainer, VTX >:

Public Types
typedef std::vector< VTX >	Product
typedef InputContainer::value_type	TRK
Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T...>	CacheTypes
typedef CacheTypes::GlobalCache	GlobalCache
typedef AbilityChecker< T...>	HasAbility
typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache
typedef LuminosityBlockContextT < LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext
typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache
typedef CacheTypes::RunCache	RunCache
typedef RunContextT< RunCache, GlobalCache >	RunContext
typedef CacheTypes::RunSummaryCache	RunSummaryCache
Public Types inherited from edm::stream::EDProducerBase
typedef EDProducerAdaptorBase	ModuleType
Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Member Functions
virtual void	produce (edm::Event &event, const edm::EventSetup &es) override
	TemplatedInclusiveVertexFinder (const edm::ParameterSet &params)
Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default
Public Member Functions inherited from edm::stream::EDProducerBase
	EDProducerBase ()
ModuleDescription const &	moduleDescription () const
virtual	~EDProducerBase ()
Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...
virtual	~ProducerBase ()
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
	EDConsumerBase ()
ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Private Member Functions
std::pair< std::vector < reco::TransientTrack > , GlobalPoint >	nearTracks (const reco::TransientTrack &seed, const std::vector< reco::TransientTrack > &tracks, const reco::Vertex &primaryVertex) const
bool	trackFilter (const reco::Track &track) const

Private Attributes
std::auto_ptr < TracksClusteringFromDisplacedSeed >	clusterizer
double	fitterRatio
double	fitterSigmacut
double	fitterTini
double	maxLIP
unsigned int	maxNTracks
unsigned int	minHits
double	minPt
edm::EDGetTokenT< reco::BeamSpot >	token_beamSpot
edm::EDGetTokenT < reco::VertexCollection >	token_primaryVertex
edm::EDGetTokenT< InputContainer >	token_tracks
bool	useVertexFitter
bool	useVertexReco
double	vertexMinAngleCosine
double	vertexMinDLen2DSig
double	vertexMinDLenSig
std::auto_ptr < VertexReconstructor >	vtxReco

Additional Inherited Members
Static Public Member Functions inherited from edm::stream::EDProducerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
Protected Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

template<class InputContainer, class VTX>
class TemplatedInclusiveVertexFinder< InputContainer, VTX >

Definition at line 36 of file InclusiveVertexFinder.h.

Member Typedef Documentation

template<class InputContainer , class VTX >

typedef std::vector<VTX> TemplatedInclusiveVertexFinder< InputContainer, VTX >::Product

Definition at line 38 of file InclusiveVertexFinder.h.

template<class InputContainer , class VTX >

typedef InputContainer::value_type TemplatedInclusiveVertexFinder< InputContainer, VTX >::TRK

Definition at line 39 of file InclusiveVertexFinder.h.

Constructor & Destructor Documentation

template<class InputContainer , class VTX >

TemplatedInclusiveVertexFinder< InputContainer, VTX >::TemplatedInclusiveVertexFinder

(

const edm::ParameterSet &

params

)

Definition at line 68 of file InclusiveVertexFinder.h.

References edm::ParameterSet::getParameter(), TemplatedInclusiveVertexFinder< InputContainer, VTX >::token_beamSpot, TemplatedInclusiveVertexFinder< InputContainer, VTX >::token_primaryVertex, and TemplatedInclusiveVertexFinder< InputContainer, VTX >::token_tracks.

                                                                                                                :
    minHits(params.getParameter<unsigned int>("minHits")),
    maxNTracks(params.getParameter<unsigned int>("maxNTracks")),
        maxLIP(params.getParameter<double>("maximumLongitudinalImpactParameter")),
    minPt(params.getParameter<double>("minPt")), //0.8
        vertexMinAngleCosine(params.getParameter<double>("vertexMinAngleCosine")), //0.98
        vertexMinDLen2DSig(params.getParameter<double>("vertexMinDLen2DSig")), //2.5
        vertexMinDLenSig(params.getParameter<double>("vertexMinDLenSig")), //0.5
        fitterSigmacut(params.getParameter<double>("fitterSigmacut")),
        fitterTini(params.getParameter<double>("fitterTini")),
        fitterRatio(params.getParameter<double>("fitterRatio")),
    useVertexFitter(params.getParameter<bool>("useDirectVertexFitter")),
    useVertexReco(params.getParameter<bool>("useVertexReco")),
    vtxReco(new ConfigurableVertexReconstructor(params.getParameter<edm::ParameterSet>("vertexReco"))),
        clusterizer(new TracksClusteringFromDisplacedSeed(params.getParameter<edm::ParameterSet>("clusterizer")))

{
    token_beamSpot = consumes<reco::BeamSpot>(params.getParameter<edm::InputTag>("beamSpot"));
    token_primaryVertex = consumes<reco::VertexCollection>(params.getParameter<edm::InputTag>("primaryVertices"));
    token_tracks = consumes<InputContainer>(params.getParameter<edm::InputTag>("tracks"));
    produces<Product>();
    //produces<reco::VertexCollection>("multi");
}

Member Function Documentation

template<class InputContainer , class VTX >

std::pair<std::vector<reco::TransientTrack>,GlobalPoint> TemplatedInclusiveVertexFinder< InputContainer, VTX >::nearTracks ( const reco::TransientTrack &  seed, const std::vector< reco::TransientTrack > &  tracks, const reco::Vertex &  primaryVertex  ) const

private

template<class InputContainer , class VTX >

void TemplatedInclusiveVertexFinder< InputContainer, VTX >::produce ( edm::Event &  event, const edm::EventSetup &  es  )

overridevirtual

Implements edm::stream::EDProducerBase.

Definition at line 103 of file InclusiveVertexFinder.h.

References funct::abs(), SiPixelRawToDigiRegional_cfi::beamSpot, tthelpers::buildTT(), HLT_25ns14e33_v1_cff::clusterizer, HLT_25ns14e33_v1_cff::clusters, gather_cfg::cout, reco::Vertex::covariance(), dir, Vector3DBase< T, FrameTag >::dot(), reco::TrackBase::dz(), Measurement1D::error(), HLT_25ns14e33_v1_cff::fitterRatio, HLT_25ns14e33_v1_cff::fitterSigmacut, HLT_25ns14e33_v1_cff::fitterTini, edm::EventSetup::get(), i, TransientVertex::isValid(), reco::TransientTrack::isValid(), j, HLT_25ns14e33_v1_cff::maxNTracks, reco::Vertex::position(), HLT_25ns14e33_v1_cff::primaryVertices, MetAnalyzer::pv(), dt_dqm_sourceclient_common_cff::reco, reco::TransientTrack::setBeamSpot(), Measurement1D::significance(), reco::TransientTrack::track(), valSkim_cff::trackFilter, testEve_cfg::tracks, groupFilesInBlocks::tt, csvLumiCalc::unit, Vector3DBase< T, FrameTag >::unit(), reco::BeamSpot::Unknown, HLT_25ns14e33_v1_cff::useVertexReco, findQualityFiles::v, Measurement1D::value(), AdaptiveVertexFitter::vertex(), HLT_25ns14e33_v1_cff::vertexMinAngleCosine, HLT_25ns14e33_v1_cff::vertexMinDLen2DSig, HLT_25ns14e33_v1_cff::vertexMinDLenSig, HLT_25ns14e33_v1_cff::vertices, and w.

{
    using namespace reco;

  VertexDistance3D vdist;
  VertexDistanceXY vdist2d;
  MultiVertexFitter theMultiVertexFitter;
  AdaptiveVertexFitter theAdaptiveFitter(
                                            GeometricAnnealing(fitterSigmacut, fitterTini, fitterRatio),
                                            DefaultLinearizationPointFinder(),
                                            KalmanVertexUpdator<5>(),
                                            KalmanVertexTrackCompatibilityEstimator<5>(),
                                            KalmanVertexSmoother() );


    edm::Handle<BeamSpot> beamSpot;
    event.getByToken(token_beamSpot,beamSpot);

    edm::Handle<VertexCollection> primaryVertices;
    event.getByToken(token_primaryVertex, primaryVertices);

    edm::Handle<InputContainer> tracks;
    event.getByToken(token_tracks, tracks);

    edm::ESHandle<TransientTrackBuilder> trackBuilder;
    es.get<TransientTrackRecord>().get("TransientTrackBuilder",
                                       trackBuilder);


        std::auto_ptr<Product> recoVertices(new Product);
        if(primaryVertices->size()!=0) {
     
    const reco::Vertex &pv = (*primaryVertices)[0];
    GlobalPoint ppv(pv.position().x(),pv.position().y(),pv.position().z());
        
    std::vector<TransientTrack> tts;
        //Fill transient track vector 
    for(typename InputContainer::const_iterator track = tracks->begin();
        track != tracks->end(); ++track) {
//TransientTrack tt = trackBuilder->build(ref);
//TrackRef ref(tracks, track - tracks->begin());
        TransientTrack tt(tthelpers::buildTT(tracks,trackBuilder, track - tracks->begin()));
        if(!tt.isValid()) continue;
        if (!trackFilter(tt.track()))
            continue;
                if( std::abs(tt.track().dz(pv.position())) > maxLIP)
            continue;
        tt.setBeamSpot(*beamSpot);
        tts.push_back(tt);
    }
        std::vector<TracksClusteringFromDisplacedSeed::Cluster> clusters = clusterizer->clusters(pv,tts);

        //Create BS object from PV to feed in the AVR
    BeamSpot::CovarianceMatrix cov;
    for(unsigned int i = 0; i < 7; i++) {
        for(unsigned int j = 0; j < 7; j++) {
            if (i < 3 && j < 3)
                cov(i, j) = pv.covariance(i, j);
            else
                cov(i, j) = 0.0;
        }
    }
    BeamSpot bs(pv.position(), 0.0, 0.0, 0.0, 0.0, cov, BeamSpot::Unknown);


        int i=0;
#ifdef VTXDEBUG

    std::cout <<  "CLUSTERS " << clusters.size() << std::endl; 
#endif

    for(std::vector<TracksClusteringFromDisplacedSeed::Cluster>::iterator cluster = clusters.begin();
        cluster != clusters.end(); ++cluster,++i)
        {
                if(cluster->tracks.size() < 2 || cluster->tracks.size() > maxNTracks ) 
             continue;
        std::vector<TransientVertex> vertices;
        if(useVertexReco) {
            vertices = vtxReco->vertices(cluster->tracks, bs);  // attempt with config given reconstructor
        }
                TransientVertex singleFitVertex;
        if(useVertexFitter) {
            singleFitVertex = theAdaptiveFitter.vertex(cluster->tracks,cluster->seedPoint); //attempt with direct fitting
            if(singleFitVertex.isValid())
                vertices.push_back(singleFitVertex);
        }
        for(std::vector<TransientVertex>::const_iterator v = vertices.begin();
            v != vertices.end(); ++v) {
            Measurement1D dlen= vdist.distance(pv,*v);
            Measurement1D dlen2= vdist2d.distance(pv,*v);
#ifdef VTXDEBUG
            VTX vv(*v);
            std::cout << "V chi2/n: " << v->normalisedChiSquared() << " ndof: " <<v->degreesOfFreedom() ;
            std::cout << " dlen: " << dlen.value() << " error: " << dlen.error() << " signif: " << dlen.significance();
            std::cout << " dlen2: " << dlen2.value() << " error2: " << dlen2.error() << " signif2: " << dlen2.significance();
            std::cout << " pos: " << vv.position() << " error: " <<vv.xError() << " " << vv.yError() << " " << vv.zError() << std::endl;
#endif
            GlobalVector dir;  
            std::vector<reco::TransientTrack> ts = v->originalTracks();
            for(std::vector<reco::TransientTrack>::const_iterator i = ts.begin();
                    i != ts.end(); ++i) {
                float w = v->trackWeight(*i);
                if (w > 0.5) dir+=i->impactPointState().globalDirection();
#ifdef VTXDEBUG
                std::cout << "\t[" << (*i).track().pt() << ": "
                    << (*i).track().eta() << ", "
                    << (*i).track().phi() << "], "
                    << w << std::endl;
#endif
                        }
            GlobalPoint sv((*v).position().x(),(*v).position().y(),(*v).position().z());
            float vscal = dir.unit().dot((sv-ppv).unit()) ;
            //                        std::cout << "Vscal: " <<  vscal << std::endl;
            if(dlen.significance() > vertexMinDLenSig  && vscal > vertexMinAngleCosine &&  v->normalisedChiSquared() < 10 && dlen2.significance() > vertexMinDLen2DSig)
            {    
                recoVertices->push_back(*v);
#ifdef VTXDEBUG
                            std::cout << "ADDED" << std::endl;
#endif
                        }
                      
                   }
        }
#ifdef VTXDEBUG

        std::cout <<  "Final put  " << recoVertices->size() << std::endl;
#endif  
        }
 
    event.put(recoVertices);

}

template<class InputContainer , class VTX >

bool TemplatedInclusiveVertexFinder< InputContainer, VTX >::trackFilter ( const reco::Track &  track ) const

private

Definition at line 92 of file InclusiveVertexFinder.h.

References reco::TrackBase::hitPattern(), OccupancyClient_cfi::minHits, HLT_25ns14e33_v1_cff::minPt, reco::HitPattern::numberOfValidHits(), and reco::TrackBase::pt().

{
    if (track.hitPattern().numberOfValidHits() < (int)minHits)
        return false;
    if (track.pt() < minPt )
        return false;
 
    return true;
}

Member Data Documentation

template<class InputContainer , class VTX >

std::auto_ptr<TracksClusteringFromDisplacedSeed> TemplatedInclusiveVertexFinder< InputContainer, VTX >::clusterizer

private

Definition at line 64 of file InclusiveVertexFinder.h.

template<class InputContainer , class VTX >

double TemplatedInclusiveVertexFinder< InputContainer, VTX >::fitterRatio

private

Definition at line 60 of file InclusiveVertexFinder.h.

template<class InputContainer , class VTX >

double TemplatedInclusiveVertexFinder< InputContainer, VTX >::fitterSigmacut

private

Definition at line 58 of file InclusiveVertexFinder.h.

template<class InputContainer , class VTX >

double TemplatedInclusiveVertexFinder< InputContainer, VTX >::fitterTini

private

Definition at line 59 of file InclusiveVertexFinder.h.

template<class InputContainer , class VTX >

double TemplatedInclusiveVertexFinder< InputContainer, VTX >::maxLIP

private

Definition at line 53 of file InclusiveVertexFinder.h.

template<class InputContainer , class VTX >

unsigned int TemplatedInclusiveVertexFinder< InputContainer, VTX >::maxNTracks

private

Definition at line 52 of file InclusiveVertexFinder.h.

template<class InputContainer , class VTX >

unsigned int TemplatedInclusiveVertexFinder< InputContainer, VTX >::minHits

private

Definition at line 51 of file InclusiveVertexFinder.h.

template<class InputContainer , class VTX >

double TemplatedInclusiveVertexFinder< InputContainer, VTX >::minPt

private

Definition at line 54 of file InclusiveVertexFinder.h.

template<class InputContainer , class VTX >

edm::EDGetTokenT<reco::BeamSpot> TemplatedInclusiveVertexFinder< InputContainer, VTX >::token_beamSpot

private

Definition at line 48 of file InclusiveVertexFinder.h.

Referenced by TemplatedInclusiveVertexFinder< InputContainer, VTX >::TemplatedInclusiveVertexFinder().

template<class InputContainer , class VTX >

edm::EDGetTokenT<reco::VertexCollection> TemplatedInclusiveVertexFinder< InputContainer, VTX >::token_primaryVertex

private

Definition at line 49 of file InclusiveVertexFinder.h.

Referenced by TemplatedInclusiveVertexFinder< InputContainer, VTX >::TemplatedInclusiveVertexFinder().

template<class InputContainer , class VTX >

edm::EDGetTokenT<InputContainer> TemplatedInclusiveVertexFinder< InputContainer, VTX >::token_tracks

private

Definition at line 50 of file InclusiveVertexFinder.h.

Referenced by TemplatedInclusiveVertexFinder< InputContainer, VTX >::TemplatedInclusiveVertexFinder().

template<class InputContainer , class VTX >

bool TemplatedInclusiveVertexFinder< InputContainer, VTX >::useVertexFitter

private

Definition at line 61 of file InclusiveVertexFinder.h.

template<class InputContainer , class VTX >

bool TemplatedInclusiveVertexFinder< InputContainer, VTX >::useVertexReco

private

Definition at line 62 of file InclusiveVertexFinder.h.

template<class InputContainer , class VTX >

double TemplatedInclusiveVertexFinder< InputContainer, VTX >::vertexMinAngleCosine

private

Definition at line 55 of file InclusiveVertexFinder.h.

template<class InputContainer , class VTX >

double TemplatedInclusiveVertexFinder< InputContainer, VTX >::vertexMinDLen2DSig

private

Definition at line 56 of file InclusiveVertexFinder.h.

template<class InputContainer , class VTX >

double TemplatedInclusiveVertexFinder< InputContainer, VTX >::vertexMinDLenSig

private

Definition at line 57 of file InclusiveVertexFinder.h.

template<class InputContainer , class VTX >

std::auto_ptr<VertexReconstructor> TemplatedInclusiveVertexFinder< InputContainer, VTX >::vtxReco

private

Definition at line 63 of file InclusiveVertexFinder.h.