Inheritance diagram for InclusiveVertexFinder:

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

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

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

virtual	~EDProducer ()

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

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

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

virtual	~ProducerBase ()

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
edm::InputTag	beamSpotCollection

std::auto_ptr < TracksClusteringFromDisplacedSeed >	clusterizer

double	fitterRatio

double	fitterSigmacut

double	fitterTini

double	maxLIP

unsigned int	maxNTracks

unsigned int	minHits

double	minPt

edm::InputTag	primaryVertexCollection

edm::InputTag	trackCollection

bool	useVertexFitter

bool	useVertexReco

double	vertexMinAngleCosine

double	vertexMinDLen2DSig

double	vertexMinDLenSig

std::auto_ptr < VertexReconstructor >	vtxReco

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

typedef WorkerT< EDProducer >	WorkerType

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::EDProducer
CurrentProcessingContext const *	currentContext () const

Protected Member Functions inherited from edm::ProducerBase
template<class TProducer , class TMethod >
void	callWhenNewProductsRegistered (TProducer *iProd, TMethod iMethod)

Detailed Description

Definition at line 33 of file InclusiveVertexFinder.cc.

Constructor & Destructor Documentation

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

Definition at line 63 of file InclusiveVertexFinder.cc.

                                                                           :
     beamSpotCollection(params.getParameter<edm::InputTag>("beamSpot")),
     primaryVertexCollection(params.getParameter<edm::InputTag>("primaryVertices")),
     trackCollection(params.getParameter<edm::InputTag>("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")))
 
 {
     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
	)

virtual

Implements edm::EDProducer.

Definition at line 98 of file InclusiveVertexFinder.cc.

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

 {
     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.getByLabel(beamSpotCollection, beamSpot);
 
     edm::Handle<VertexCollection> primaryVertices;
     event.getByLabel(primaryVertexCollection, primaryVertices);
 
     edm::Handle<TrackCollection> tracks;
     event.getByLabel(trackCollection, 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() < 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) {
 //          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 87 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

edm::InputTag InclusiveVertexFinder::beamSpotCollection

private

Definition at line 43 of file InclusiveVertexFinder.cc.

Referenced by produce().

std::auto_ptr<TracksClusteringFromDisplacedSeed> InclusiveVertexFinder::clusterizer

private

Definition at line 59 of file InclusiveVertexFinder.cc.

Referenced by produce().

double InclusiveVertexFinder::fitterRatio

private

Definition at line 55 of file InclusiveVertexFinder.cc.

Referenced by produce().

double InclusiveVertexFinder::fitterSigmacut

private

Definition at line 53 of file InclusiveVertexFinder.cc.

Referenced by produce().

double InclusiveVertexFinder::fitterTini

private

Definition at line 54 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::InputTag InclusiveVertexFinder::primaryVertexCollection

private

Definition at line 44 of file InclusiveVertexFinder.cc.

Referenced by produce().

edm::InputTag InclusiveVertexFinder::trackCollection

private

Definition at line 45 of file InclusiveVertexFinder.cc.

Referenced by produce().

bool InclusiveVertexFinder::useVertexFitter

private

Definition at line 56 of file InclusiveVertexFinder.cc.

Referenced by produce().

bool InclusiveVertexFinder::useVertexReco

private

Definition at line 57 of file InclusiveVertexFinder.cc.

Referenced by 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 58 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