#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 Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector < edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, std::unordered_multimap< std::string, edm::ProductResolverIndex > const &iIndicies, std::string const &moduleLabel)

virtual	~ProducerBase () noexcept(false)

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector < ProductResolverIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

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::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_FULL_cff::clusterizer, HLT_FULL_cff::clusters, gather_cfg::cout, reco::Vertex::covariance(), dir, Vector3DBase< T, FrameTag >::dot(), reco::TrackBase::dz(), Measurement1D::error(), HLT_FULL_cff::fitterRatio, HLT_FULL_cff::fitterSigmacut, HLT_FULL_cff::fitterTini, edm::EventSetup::get(), i, reco::TransientTrack::isValid(), TransientVertex::isValid(), j, HLT_FULL_cff::maxNTracks, eostools::move(), reco::Vertex::position(), HLT_FULL_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_FULL_cff::useVertexReco, findQualityFiles::v, Measurement1D::value(), AdaptiveVertexFitter::vertex(), HLT_FULL_cff::vertexMinAngleCosine, HLT_FULL_cff::vertexMinDLen2DSig, HLT_FULL_cff::vertexMinDLenSig, HLT_FULL_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);
 
 
         auto recoVertices = std::make_unique<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(std::move(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_FULL_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.

Public Types

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

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

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation

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