Inheritance diagram for InclusiveVertexFinder:

Public Member Functions
	InclusiveVertexFinder (const edm::ParameterSet &params)

virtual void	produce (edm::Event &event, const edm::EventSetup &es) override

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducer ()

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
	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

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::TrackRef &track) const

Private Attributes
std::auto_ptr < TracksClusteringFromDisplacedSeed >	clusterizer

double	maxLIP

unsigned int	maxNTracks

unsigned int	minHits

double	minPt

edm::EDGetTokenT< reco::BeamSpot >	token_beamSpot

edm::EDGetTokenT < reco::VertexCollection >	token_primaryVertex

edm::EDGetTokenT < reco::TrackCollection >	token_tracks

double	vertexMinAngleCosine

double	vertexMinDLen2DSig

double	vertexMinDLenSig

std::auto_ptr < VertexReconstructor >	vtxReco

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Static Public Member Functions inherited from edm::EDProducer
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

Definition at line 33 of file InclusiveVertexFinder.cc.

Constructor & Destructor Documentation

InclusiveVertexFinder::InclusiveVertexFinder ( const edm::ParameterSet & params )

Definition at line 59 of file InclusiveVertexFinder.cc.

References edm::ParameterSet::getParameter(), token_beamSpot, token_primaryVertex, and 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
     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<reco::TrackCollection>(params.getParameter<edm::InputTag>("tracks"));
     produces<reco::VertexCollection>();
     //produces<reco::VertexCollection>("multi");
 }

Member Function Documentation

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

private

void InclusiveVertexFinder::produce	(	edm::Event &	event,
		const edm::EventSetup &	es
	)

overridevirtual

Implements edm::EDProducer.

Definition at line 89 of file InclusiveVertexFinder.cc.

References funct::abs(), SiPixelRawToDigiRegional_cfi::beamSpot, clusterizer, gather_cfg::cout, reco::Vertex::covariance(), dir, Vector3DBase< T, FrameTag >::dot(), Measurement1D::error(), edm::EventSetup::get(), i, TransientVertex::isValid(), j, maxLIP, maxNTracks, reco::Vertex::position(), dt_dqm_sourceclient_common_cff::reco, reco::TransientTrack::setBeamSpot(), Measurement1D::significance(), lumiQTWidget::t, token_beamSpot, token_primaryVertex, token_tracks, trackFilter(), testEve_cfg::tracks, groupFilesInBlocks::tt, csvLumiCalc::unit, Vector3DBase< T, FrameTag >::unit(), Unknown, findQualityFiles::v, Measurement1D::value(), AdaptiveVertexFitter::vertex(), vertexMinAngleCosine, vertexMinDLen2DSig, vertexMinDLenSig, vtxReco, w(), reco::Vertex::xError(), reco::Vertex::yError(), and reco::Vertex::zError().

 {
     using namespace reco;
 
   double sigmacut = 3.0;
   double Tini = 256.;
   double ratio = 0.25;
   VertexDistance3D vdist;
   VertexDistanceXY vdist2d;
   MultiVertexFitter theMultiVertexFitter;
   AdaptiveVertexFitter theAdaptiveFitter(
                                             GeometricAnnealing(sigmacut, Tini, ratio),
                                             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<TrackCollection> tracks;
     event.getByToken(token_tracks, tracks);
 
     edm::ESHandle<TransientTrackBuilder> trackBuilder;
     es.get<TransientTrackRecord>().get("TransientTrackBuilder",
                                        trackBuilder);
 
 
         std::auto_ptr<VertexCollection> recoVertices(new VertexCollection);
         if(primaryVertices->size()!=0) {
      
     const reco::Vertex &pv = (*primaryVertices)[0];
         
     std::vector<TransientTrack> tts;
         //Fill transient track vector 
     for(TrackCollection::const_iterator track = tracks->begin();
         track != tracks->end(); ++track) {
         TrackRef ref(tracks, track - tracks->begin());
         if (!trackFilter(ref))
             continue;
                 if( std::abs(ref->dz(pv.position())) > maxLIP)
             continue;
         TransientTrack tt = trackBuilder->build(ref);
         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() == 0 || cluster->tracks.size() > maxNTracks ) 
              continue;
         
             cluster->tracks.push_back(cluster->seedingTrack); //add the seed to the list of tracks to fit
         std::vector<TransientVertex> vertices;
         vertices = vtxReco->vertices(cluster->tracks, bs);  // attempt with config given reconstructor
                 TransientVertex singleFitVertex;
                 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) {
 //          if(v->degreesOfFreedom() > 0.2)
                         {
                          Measurement1D dlen= vdist.distance(pv,*v);
                          Measurement1D dlen2= vdist2d.distance(pv,*v);
              reco::Vertex vv(*v);
 #ifdef VTXDEBUG
                        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) {
                                 reco::TrackRef t = i->trackBaseRef().castTo<reco::TrackRef>();
                                 float w = v->trackWeight(*i);
                                 if (w > 0.5) dir+=i->impactPointState().globalDirection();
 #ifdef VTXDEBUG
                                 std::cout << "\t[" << (*t).pt() << ": "
                                           << (*t).eta() << ", "
                                           << (*t).phi() << "], "
                                           << w << std::endl;
 #endif
                         }
                GlobalPoint ppv(pv.position().x(),pv.position().y(),pv.position().z());
                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);
 
 }

bool InclusiveVertexFinder::trackFilter ( const reco::TrackRef & track ) const

private

Definition at line 78 of file InclusiveVertexFinder.cc.

References minHits, and minPt.

Referenced by produce().

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

Member Data Documentation

std::auto_ptr<TracksClusteringFromDisplacedSeed> InclusiveVertexFinder::clusterizer

private

Definition at line 55 of file InclusiveVertexFinder.cc.

Referenced by produce().

double InclusiveVertexFinder::maxLIP

private

Definition at line 48 of file InclusiveVertexFinder.cc.

Referenced by produce().

unsigned int InclusiveVertexFinder::maxNTracks

private

Definition at line 47 of file InclusiveVertexFinder.cc.

Referenced by produce().

unsigned int InclusiveVertexFinder::minHits

private

Definition at line 46 of file InclusiveVertexFinder.cc.

Referenced by trackFilter().

double InclusiveVertexFinder::minPt

private

Definition at line 49 of file InclusiveVertexFinder.cc.

Referenced by trackFilter().

edm::EDGetTokenT<reco::BeamSpot> InclusiveVertexFinder::token_beamSpot

private

Definition at line 43 of file InclusiveVertexFinder.cc.

Referenced by InclusiveVertexFinder(), and produce().

edm::EDGetTokenT<reco::VertexCollection> InclusiveVertexFinder::token_primaryVertex

private

Definition at line 44 of file InclusiveVertexFinder.cc.

Referenced by InclusiveVertexFinder(), and produce().

edm::EDGetTokenT<reco::TrackCollection> InclusiveVertexFinder::token_tracks

private

Definition at line 45 of file InclusiveVertexFinder.cc.

Referenced by InclusiveVertexFinder(), and produce().

double InclusiveVertexFinder::vertexMinAngleCosine

private

Definition at line 50 of file InclusiveVertexFinder.cc.

Referenced by produce().

double InclusiveVertexFinder::vertexMinDLen2DSig

private

Definition at line 51 of file InclusiveVertexFinder.cc.

Referenced by produce().

double InclusiveVertexFinder::vertexMinDLenSig

private

Definition at line 52 of file InclusiveVertexFinder.cc.

Referenced by produce().

std::auto_ptr<VertexReconstructor> InclusiveVertexFinder::vtxReco

private

Definition at line 54 of file InclusiveVertexFinder.cc.

Referenced by produce().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation