#include <PFDisplacedVertexCandidateFinder.h>

Public Types
typedef reco::PFDisplacedVertexCandidateCollection::const_iterator	IBC

typedef std::list< reco::TrackBaseRef >::iterator	IE

typedef std::list< reco::TrackBaseRef >::const_iterator	IEC

typedef std::vector< bool >	Mask
	Mask used to spot if a track is free or not. More...

Public Member Functions
void	findDisplacedVertexCandidates ()
	-----— Main function which find vertices -----— /// More...

	PFDisplacedVertexCandidateFinder ()

void	setDebug (bool debug)
	sets debug printout flag More...

void	setInput (const edm::Handle< reco::TrackCollection > &trackh, const MagneticField *magField)

void	setParameters (double dcaCut, double primaryVertexCut, double dcaPInnerHitCut, const edm::ParameterSet &ps_trk)
	---—— Set different algo parameters —— /// More...

void	setPrimaryVertex (edm::Handle< reco::VertexCollection > mainVertexHandle, edm::Handle< reco::BeamSpot > beamSpotHandle)

std::unique_ptr< reco::PFDisplacedVertexCandidateCollection >	transferVertexCandidates ()

const std::unique_ptr< reco::PFDisplacedVertexCandidateCollection > &	vertexCandidates () const

	~PFDisplacedVertexCandidateFinder ()

Private Member Functions
IE	associate (IE next, IE last, reco::PFDisplacedVertexCandidate &tempVertexCandidate)
	--—— Different steps of the finder algorithm --—— /// More...

GlobalTrajectoryParameters	getGlobalTrajectoryParameters (const reco::Track *) const
	--—— TOOLS --—— // More...

bool	goodPtResolution (const reco::TrackBaseRef &trackref) const
	Quality Criterion on the Pt resolution to select a Track. More...

bool	isSelected (const reco::TrackBaseRef &trackref)

void	link (const reco::TrackBaseRef &el1, const reco::TrackBaseRef &el2, double &dist, GlobalPoint &P, reco::PFDisplacedVertexCandidate::VertexLinkTest &linktest)
	Check whether 2 elements are linked and fill the link parameters. More...

void	packLinks (reco::PFDisplacedVertexCandidate &vertexCandidate)

Private Attributes
double	dcaCut_
	–— Algo parameters for the vertex finder -— /// More...

double	dcaPInnerHitCut2_
	Maximum distance between the DCA Point and the inner hit of the track. More...

bool	debug_
	if true, debug printouts activated More...

double	dxy_

std::list< reco::TrackBaseRef >	eventTracks_
	The track refs. More...

std::vector< GlobalTrajectoryParameters >	eventTrackTrajectories_

const MagneticField *	magField_

double	nChi2_max_

double	primaryVertexCut2_
	Do not reconstruct vertices wich are too close to the beam pipe. More...

double	pt_min_

double	pt_min_prim_

math::XYZPoint	pvtx_

TwoTrackMinimumDistance	theMinimum_

Mask	trackMask_

std::unique_ptr< reco::PFDisplacedVertexCandidateCollection >	vertexCandidates_
	-----— Members -----— /// More...

unsigned	vertexCandidatesSize_

Friends
std::ostream &	operator<< (std::ostream &, const PFDisplacedVertexCandidateFinder &)

Detailed Description

Definition at line 24 of file PFDisplacedVertexCandidateFinder.h.

Member Typedef Documentation

typedef reco::PFDisplacedVertexCandidateCollection::const_iterator PFDisplacedVertexCandidateFinder::IBC

Definition at line 38 of file PFDisplacedVertexCandidateFinder.h.

typedef std::list< reco::TrackBaseRef >::iterator PFDisplacedVertexCandidateFinder::IE

Definition at line 36 of file PFDisplacedVertexCandidateFinder.h.

typedef std::list< reco::TrackBaseRef >::const_iterator PFDisplacedVertexCandidateFinder::IEC

Definition at line 37 of file PFDisplacedVertexCandidateFinder.h.

typedef std::vector<bool> PFDisplacedVertexCandidateFinder::Mask

Mask used to spot if a track is free or not.

Definition at line 34 of file PFDisplacedVertexCandidateFinder.h.

Constructor & Destructor Documentation

PFDisplacedVertexCandidateFinder::PFDisplacedVertexCandidateFinder ( )

Definition at line 15 of file PFDisplacedVertexCandidateFinder.cc.

References TwoTrackMinimumDistance::SlowMode, and theMinimum_.

                                                                    : 
   vertexCandidates_( new PFDisplacedVertexCandidateCollection ),
   dcaCut_(1000),
   primaryVertexCut2_(0.0),
   dcaPInnerHitCut2_(1000.0),
   vertexCandidatesSize_(50),
   debug_(false) {
 
   TwoTrackMinimumDistance theMinimum(TwoTrackMinimumDistance::SlowMode);
   theMinimum_ = theMinimum;
 }

PFDisplacedVertexCandidateFinder::~PFDisplacedVertexCandidateFinder ( )

Definition at line 28 of file PFDisplacedVertexCandidateFinder.cc.

References gather_cfg::cout, debug_, and eventTracks_.

                                                                     {
 
 #ifdef PFLOW_DEBUG
   if(debug_)
     cout<<"~PFDisplacedVertexCandidateFinder - number of remaining elements: "
     <<eventTracks_.size()<<endl;
 #endif
   
 }

Member Function Documentation

PFDisplacedVertexCandidateFinder::IE PFDisplacedVertexCandidateFinder::associate	(	IE	next,
		IE	last,
		reco::PFDisplacedVertexCandidate &	tempVertexCandidate
	)

private

--—— Different steps of the finder algorithm --—— ///

Recursive procedure to associate tracks together

Definition at line 129 of file PFDisplacedVertexCandidateFinder.cc.

References reco::PFDisplacedVertexCandidate::addElement(), gather_cfg::cout, debug_, eventTracks_, mps_fire::i, link(), GetRecoTauVFromDQM_MC_cff::next, trackMask_, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by findDisplacedVertexCandidates(), and vertexCandidates().

                                                                          {
 
     
 #ifdef PFLOW_DEBUG
   if(debug_ ) cout<<"== Start the association procedure =="<<endl;
 #endif
 
   if( last!= eventTracks_.end() ) {
     double dist = -1;
     GlobalPoint P(0,0,0);
     PFDisplacedVertexCandidate::VertexLinkTest linktest;
     link( (*last), (*next), dist, P, linktest); 
 
 
     if(dist<-0.5) {
 #ifdef PFLOW_DEBUG
       if(debug_ ) cout<<"link failed"<<endl;
 #endif
       return ++next; // association failed
     }
     else {
       // add next element to the current pflowblock
       tempVertexCandidate.addElement( (*next));  
       trackMask_[(*next).key()] = false;
 #ifdef PFLOW_DEBUG   
       if(debug_ ) 
     cout<<"link parameters "
         << " *next = " << (*next).key()
         << " *last = " << (*last).key()
         << "  dist = " << dist
         << " P.x = " << P.x() 
         << " P.y = " << P.y()
         << " P.z = " << P.z()
         << endl;
 #endif
     }
   }
   else {
     // add next element to this eflowblock
 #ifdef PFLOW_DEBUG   
     if(debug_ ) cout<<"adding to block element "
             << (*next).key()
             <<endl;
 #endif
     tempVertexCandidate.addElement( (*next));  
     trackMask_[(*next).key()] = false;
 
 
   }
 
   // recursive call: associate next and other unused elements 
 #ifdef PFLOW_DEBUG  
   if(debug_ ) {
     for(unsigned i=0; i<trackMask_.size(); i++) 
       cout << " Mask[" << i << "] = " << trackMask_[i];
     cout << "" << endl;
   }
 #endif
 
   for(IE ie = eventTracks_.begin(); 
       ie != eventTracks_.end();) {
     
     if( ie == last || ie == next ) { ++ie; continue;} 
 
     // *ie already included to a block
     if( !trackMask_[(*ie).key()] ) { ++ie; continue;}
 
 #ifdef PFLOW_DEBUG      
     if(debug_ ) cout<<"calling associate "
             << (*next).key()
             <<" & "
             << (*ie).key()
             <<endl;
 #endif
     ie = associate(next, ie, tempVertexCandidate);
     
   }       
 
 #ifdef PFLOW_DEBUG
   if(debug_ ) {
     cout<<"**** removing element "<<endl;
   }
 #endif
 
   IE iteratorToNextFreeElement = eventTracks_.erase( next );
 
 #ifdef PFLOW_DEBUG
   if(debug_ ) cout<< "== End the association procedure ==" <<endl;
 #endif
 
   return iteratorToNextFreeElement;
 }

void PFDisplacedVertexCandidateFinder::findDisplacedVertexCandidates ( )

-----— Main function which find vertices -----— ///

Definition at line 94 of file PFDisplacedVertexCandidateFinder.cc.

References associate(), gather_cfg::cout, debug_, eventTracks_, reco::PFDisplacedVertexCandidate::isValid(), packLinks(), vertexCandidates_, and vertexCandidatesSize_.

Referenced by vertexCandidates().

                                                                 {
 
   if (debug_) cout << "========= Start Finding Displaced Vertex Candidates =========" << endl;
   
   // The vertexCandidates have not been passed to the event, and need to be cleared
   if(vertexCandidates_.get() )vertexCandidates_->clear();
   else 
     vertexCandidates_.reset( new PFDisplacedVertexCandidateCollection );
 
   vertexCandidates_->reserve(vertexCandidatesSize_);
   for(IE ie = eventTracks_.begin(); 
       ie != eventTracks_.end();) {
 
     // Run the recursive procedure to find all tracks link together 
     // In one blob called Candidate
 
     PFDisplacedVertexCandidate tempVertexCandidate;
 
     ie = associate( eventTracks_.end() , ie, tempVertexCandidate);
 
     // Build remaining links in current block
     if(tempVertexCandidate.isValid()) {
       packLinks( tempVertexCandidate );
       vertexCandidates_->push_back( tempVertexCandidate );
     }
   }       
 
   if(debug_) cout << "========= End Finding Displaced Vertex Candidates =========" << endl;
 
 }

GlobalTrajectoryParameters PFDisplacedVertexCandidateFinder::getGlobalTrajectoryParameters ( const reco::Track * track ) const

private

--—— TOOLS --—— //

Allows to calculate the helix aproximation for a given track which may be then extrapolated to any point.

Definition at line 341 of file PFDisplacedVertexCandidateFinder.cc.

References reco::TrackBase::charge(), magField_, reco::TrackBase::momentum(), position, reco::TrackBase::vx(), reco::TrackBase::vy(), and reco::TrackBase::vz().

Referenced by packLinks(), setInput(), and vertexCandidates().

 {
 
   const GlobalPoint position(track->vx(),
                  track->vy(),
                  track->vz());
  
   const GlobalVector momentum(track->momentum().x(),
                   track->momentum().y(),
                   track->momentum().z());
 
   GlobalTrajectoryParameters gtp(position,momentum,
                  track->charge(),magField_);
 
   return gtp;
 }

bool PFDisplacedVertexCandidateFinder::goodPtResolution ( const reco::TrackBaseRef & trackref ) const

private

Quality Criterion on the Pt resolution to select a Track.

Definition at line 361 of file PFDisplacedVertexCandidateFinder.cc.

References gather_cfg::cout, debug_, PVValHelper::dxy, reco::TrackBase::dxy(), dxy_, nChi2_max_, reco::TrackBase::normalizedChi2(), EnergyCorrector::pt, reco::TrackBase::pt(), pt_min_, pt_min_prim_, reco::TrackBase::ptError(), and pvtx_.

Referenced by isSelected(), and vertexCandidates().

                                                                                       {
 
   double nChi2 = trackref->normalizedChi2(); 
   double pt = trackref->pt();
   double dpt = trackref->ptError();
   double dxy = trackref->dxy(pvtx_);
 
   double pt_error = dpt/pt*100;
 
   if (debug_) cout << " PFDisplacedVertexFinder: PFrecTrack->Track Pt= "
            << pt << " dPt/Pt = " << pt_error << "% nChi2 = " << nChi2 << endl;
   if (nChi2 > nChi2_max_ || pt < pt_min_){
 
     if (debug_) cout << " PFBlockAlgo: skip badly measured or low pt track"
              << " nChi2_cut = " << 5 
              << " pt_cut = " << 0.2 << endl;
     return false;
   }
   //  cout << "dxy = " << dxy << endl; 
   if (fabs(dxy) < dxy_ && pt < pt_min_prim_) return false;
   //  if (fabs(dxy) < 0.2 && pt < 0.8) return false;
 
   
 
   return true;
 }

bool PFDisplacedVertexCandidateFinder::isSelected ( const reco::TrackBaseRef & trackref )

inlineprivate

A function which gather the information if a track is available for vertexing

Definition at line 110 of file PFDisplacedVertexCandidateFinder.h.

References goodPtResolution(), and operator<<.

Referenced by setInput().

111 { return goodPtResolution(trackref);}

PFDisplacedVertexCandidateFinder::goodPtResolution

bool goodPtResolution(const reco::TrackBaseRef &trackref) const

Quality Criterion on the Pt resolution to select a Track.

Definition: PFDisplacedVertexCandidateFinder.cc:361

void PFDisplacedVertexCandidateFinder::link	(	const reco::TrackBaseRef &	el1,
		const reco::TrackBaseRef &	el2,
		double &	dist,
		GlobalPoint &	P,
		reco::PFDisplacedVertexCandidate::VertexLinkTest &	linktest
	)

private

Check whether 2 elements are linked and fill the link parameters.

Definition at line 226 of file PFDisplacedVertexCandidateFinder.cc.

References TwoTrackMinimumDistance::calculate(), TwoTrackMinimumDistance::crossingPoint(), dcaCut_, TwoTrackMinimumDistance::distance(), reco::TrackBase::eta(), eventTrackTrajectories_, edm::RefToBase< T >::key(), reco::TrackBase::phi(), primaryVertexCut2_, reco::TrackBase::pt(), theMinimum_, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by associate(), packLinks(), and vertexCandidates().

                         {
  
   if ( fabs(el1->eta()-el2->eta()) > 1) {dist = -1; return;}
   if ( el1->pt()>2 && el2->pt()>2 && fabs(el1->phi()-el2->phi()) > 1) {dist = -1; return;}
 
   const GlobalTrajectoryParameters& gt1 = eventTrackTrajectories_[el1.key()];
   const GlobalTrajectoryParameters& gt2 = eventTrackTrajectories_[el2.key()];
   
 
   // Closest approach
   theMinimum_.calculate(gt1,gt2);
 
   // Fill the parameters
   dist    = theMinimum_.distance();
   P       = theMinimum_.crossingPoint();
   
   vertexLinkTest = PFDisplacedVertexCandidate::LINKTEST_DCA; //rechit by default 
 
 
   // Check if the link test fails
   if (dist > dcaCut_) {dist = -1; return;}
 
   // Check if the closses approach point is too close to the primary vertex/beam pipe 
   double rho2 = P.x()*P.x()+P.y()*P.y();
 
   if ( rho2 < primaryVertexCut2_) {dist = -1; return;}
 
   // Check if the closses approach point is out of range of the tracker 
 
   double tob_rho_limit2 = 10000;
   double tec_z_limit = 270;
 
   if ( rho2 > tob_rho_limit2) {dist = -1; return;}
   if ( fabs(P.z()) > tec_z_limit) {dist = -1; return;}
 
   return;
 
 
   /*
   // Check if the inner hit of one of the tracks is not too far away from the vertex
   double dx = el1->innerPosition().x()-P.x();
   double dy = el1->innerPosition().y()-P.y();
   double dz = el1->innerPosition().z()-P.z();
   double dist2 = dx*dx+dy*dy+dz*dz; 
 
   if (dist2 > dcaPInnerHitCut2_) {
     dist = -1; 
     if (debug_) cout << "track " << el1.key() << " dist to vertex " << sqrt(dist2) << endl; 
     return;
   }
 
   dx = el2->innerPosition().x()-P.x();
   dy = el2->innerPosition().y()-P.y();
   dz = el2->innerPosition().z()-P.z();
   dist2 = dx*dx+dy*dy+dz*dz; 
 
   if (dist2 > dcaPInnerHitCut2_) {
     dist = -1; 
     if (debug_) cout << "track " << el2.key() << " dist to vertex " << sqrt(dist2) << endl; 
     return;
   }
   */
   
 }

void PFDisplacedVertexCandidateFinder::packLinks ( reco::PFDisplacedVertexCandidate & vertexCandidate )

private

Compute missing links in the displacedVertexCandidates (the recursive procedure does not build all links)

Definition at line 299 of file PFDisplacedVertexCandidateFinder.cc.

References gather_cfg::cout, debug_, reco::PFDisplacedVertexCandidate::elements(), getGlobalTrajectoryParameters(), link(), and reco::PFDisplacedVertexCandidate::setLink().

Referenced by findDisplacedVertexCandidates(), and vertexCandidates().

                                                                                         {
   
   
   const vector < TrackBaseRef >& els = vertexCandidate.elements();
   
   //First Loop: update all link data
   for( unsigned i1=0; i1<els.size(); i1++ ) {
     for( unsigned i2=i1+1; i2<els.size(); i2++ ) {
       
       // no reflexive link
       if( i1==i2 ) continue;
       
       double dist = -1;
       GlobalPoint P(0,0,0); 
       PFDisplacedVertexCandidate::VertexLinkTest linktest; 
 
       link( els[i1], els[i2], dist, P, linktest);
 
 
 #ifdef PFLOW_DEBUG
       if( debug_ )
     cout << "Setting link between elements " << i1 << " key " << els[i1].key() 
          << " and " << i2 << " key " << els[i2].key() 
          << " of dist =" << dist << " computed from link test "
          << linktest << endl;
 #endif
 
       if(dist >-0.5) vertexCandidate.setLink( i1, i2, dist, P, linktest );
     }
   }
 
 }

void PFDisplacedVertexCandidateFinder::setDebug ( bool debug )

inline

sets debug printout flag

Definition at line 56 of file PFDisplacedVertexCandidateFinder.h.

References debug, debug_, and setInput().

56 {debug_ = debug;}

debug

#define debug

Definition: HDRShower.cc:19

PFDisplacedVertexCandidateFinder::debug_

bool debug_

if true, debug printouts activated

Definition: PFDisplacedVertexCandidateFinder.h:163

void PFDisplacedVertexCandidateFinder::setInput	(	const edm::Handle< reco::TrackCollection > &	trackh,
		const MagneticField *	magField
	)

Set the imput collection of tracks and calculate their trajectory parameters the Global Trajectory Parameters

Definition at line 60 of file PFDisplacedVertexCandidateFinder.cc.

References eventTracks_, eventTrackTrajectories_, edm::Ref< C, T, F >::get(), getGlobalTrajectoryParameters(), mps_fire::i, isSelected(), edm::HandleBase::isValid(), magField_, and trackMask_.

Referenced by setDebug().

                                                       {
 
   magField_ = magField;
   trackMask_.clear();
   trackMask_.resize(trackh->size());
   eventTrackTrajectories_.clear();
   eventTrackTrajectories_.resize(trackh->size());
 
   for(unsigned i=0;i<trackMask_.size(); i++) trackMask_[i] = true; 
 
   eventTracks_.clear();
   if(trackh.isValid()) {
     for(unsigned i=0;i<trackh->size(); i++) {
 
       TrackRef tref( trackh, i); 
       TrackBaseRef tbref(tref);
 
       if( !isSelected( tbref ) ) {
     trackMask_[i] = false;
     continue;      
       }
       
       const Track* trk = tref.get();
       eventTracks_.push_back( tbref );
       eventTrackTrajectories_[i] = getGlobalTrajectoryParameters(trk);
     }
   }
 }

void PFDisplacedVertexCandidateFinder::setParameters	(	double	dcaCut,
		double	primaryVertexCut,
		double	dcaPInnerHitCut,
		const edm::ParameterSet &	ps_trk
	)

inline

---—— Set different algo parameters —— ///

Sets algo parameters for the vertex candidate finder

Definition at line 44 of file PFDisplacedVertexCandidateFinder.h.

References particleFlowDisplacedVertexCandidate_cfi::dcaCut, dcaCut_, particleFlowDisplacedVertexCandidate_cfi::dcaPInnerHitCut, dcaPInnerHitCut2_, dxy_, edm::ParameterSet::getParameter(), nChi2_max_, particleFlowDisplacedVertex_cfi::primaryVertexCut, primaryVertexCut2_, pt_min_, and pt_min_prim_.

                                             {
     dcaCut_ = dcaCut;
     primaryVertexCut2_ = primaryVertexCut*primaryVertexCut;
     dcaPInnerHitCut2_ = dcaPInnerHitCut*dcaPInnerHitCut;
     nChi2_max_ = ps_trk.getParameter<double>("nChi2_max");
     pt_min_ = ps_trk.getParameter<double>("pt_min");
     pt_min_prim_ = ps_trk.getParameter<double>("pt_min_prim");
     dxy_ = ps_trk.getParameter<double>("dxy");
   }

void PFDisplacedVertexCandidateFinder::setPrimaryVertex	(	edm::Handle< reco::VertexCollection >	mainVertexHandle,
		edm::Handle< reco::BeamSpot >	beamSpotHandle
	)

Definition at line 39 of file PFDisplacedVertexCandidateFinder.cc.

References ecalDrivenElectronSeedsParameters_cff::beamSpot, edm::HandleBase::isValid(), pvtx_, reco::BeamSpot::x0(), reco::BeamSpot::y0(), and reco::BeamSpot::z0().

Referenced by vertexCandidates().

                                                                     {
 
   const math::XYZPoint beamSpot = beamSpotHandle.isValid() ? 
     math::XYZPoint(beamSpotHandle->x0(), beamSpotHandle->y0(), beamSpotHandle->z0()) : 
     math::XYZPoint(0, 0, 0);
 
   // The primary vertex is taken from the refitted list, 
   // if does not exist from the average offline beam spot position  
   // if does not exist (0,0,0) is used
   pvtx_ = mainVertexHandle.isValid() ? 
     math::XYZPoint(mainVertexHandle->begin()->x(), 
            mainVertexHandle->begin()->y(), 
            mainVertexHandle->begin()->z()) :
     beamSpot;
 }

std::unique_ptr< reco::PFDisplacedVertexCandidateCollection > PFDisplacedVertexCandidateFinder::transferVertexCandidates ( )

inline

Returns: auto_ptr to collection of DisplacedVertexCandidates

Definition at line 65 of file PFDisplacedVertexCandidateFinder.h.

References eostools::move(), and vertexCandidates_.

65 {return std::move(vertexCandidates_);}

PFDisplacedVertexCandidateFinder::vertexCandidates_

std::unique_ptr< reco::PFDisplacedVertexCandidateCollection > vertexCandidates_

-----— Members -----— ///

Definition: PFDisplacedVertexCandidateFinder.h:119

eostools.move

def move(src, dest)

Definition: eostools.py:511

const std::unique_ptr< reco::PFDisplacedVertexCandidateCollection >& PFDisplacedVertexCandidateFinder::vertexCandidates ( ) const

inline

Definition at line 67 of file PFDisplacedVertexCandidateFinder.h.

References associate(), findDisplacedVertexCandidates(), getGlobalTrajectoryParameters(), goodPtResolution(), plotBeamSpotDB::last, link(), eostools::move(), GetRecoTauVFromDQM_MC_cff::next, packLinks(), setPrimaryVertex(), and vertexCandidates_.

Referenced by operator<<().

68 {return std::move(vertexCandidates_);}

PFDisplacedVertexCandidateFinder::vertexCandidates_

std::unique_ptr< reco::PFDisplacedVertexCandidateCollection > vertexCandidates_

-----— Members -----— ///

Definition: PFDisplacedVertexCandidateFinder.h:119

eostools.move

def move(src, dest)

Definition: eostools.py:511

Friends And Related Function Documentation

std::ostream& operator<<	(	std::ostream &	out,
		const PFDisplacedVertexCandidateFinder &	a
	)

friend

Definition at line 398 of file PFDisplacedVertexCandidateFinder.cc.

Referenced by isSelected().

                                                                                 {
   if(! out) return out;
   
   out<<"====== Particle Flow Block Algorithm ======= ";
   out<<endl;
   out<<"number of unassociated elements : "<<a.eventTracks_.size()<<endl;
   out<<endl;
   
   out << " Tracks selection based on " << std::endl;
   out << " pvtx_ = " << a.pvtx_ << std::endl;
   out << " fabs(dxy) < " << a.dxy_ << " and pt < "<< a.pt_min_prim_ << std::endl;
   out << " nChi2 < " << a.nChi2_max_ << " and pt < "<< a.pt_min_ << std::endl;
 
   out<<endl;
 
 
   for(PFDisplacedVertexCandidateFinder::IEC ie = a.eventTracks_.begin(); 
       ie != a.eventTracks_.end(); ie++) {
 
     double pt = (*ie).get()->pt(); 
 
     math::XYZPoint Pi = (*ie).get()->innerPosition(); 
     math::XYZPoint Po = (*ie).get()->outerPosition(); 
 
     double innermost_radius = sqrt(Pi.x()*Pi.x() + Pi.y()*Pi.y() + Pi.z()*Pi.z());
     double outermost_radius = sqrt(Po.x()*Po.x() + Po.y()*Po.y() + Po.z()*Po.z());
     double innermost_rho = sqrt(Pi.x()*Pi.x() + Pi.y()*Pi.y());
     double outermost_rho = sqrt(Po.x()*Po.x() + Po.y()*Po.y());
     
     out<<"ie = " << (*ie).key() 
        <<" pt = " << pt
        <<" innermost hit radius = " << innermost_radius << " rho = " << innermost_rho
        <<" outermost hit radius = " << outermost_radius << " rho = " << outermost_rho
        <<endl;
   }
 
 
   const std::unique_ptr< reco::PFDisplacedVertexCandidateCollection >& vertexCandidates
     = std::move(a.vertexCandidates()); 
     
   if(!vertexCandidates.get() ) {
     out<<"vertexCandidates already transfered"<<endl;
   }
   else {
     out<<"number of vertexCandidates : "<<vertexCandidates->size()<<endl;
     out<<endl;
  
     
     for(PFDisplacedVertexCandidateFinder::IBC ib=vertexCandidates->begin(); 
     ib != vertexCandidates->end(); ib++)
       ib->Dump();
 
 
     
   }
  
   return out;
 }

Member Data Documentation

double PFDisplacedVertexCandidateFinder::dcaCut_

private

–— Algo parameters for the vertex finder -— ///

Distance of minimal approach below which two tracks are considered as linked together

Definition at line 138 of file PFDisplacedVertexCandidateFinder.h.

Referenced by link(), and setParameters().

double PFDisplacedVertexCandidateFinder::dcaPInnerHitCut2_

private

Maximum distance between the DCA Point and the inner hit of the track.

Definition at line 142 of file PFDisplacedVertexCandidateFinder.h.

Referenced by setParameters().

bool PFDisplacedVertexCandidateFinder::debug_

private

if true, debug printouts activated

Definition at line 163 of file PFDisplacedVertexCandidateFinder.h.

Referenced by associate(), findDisplacedVertexCandidates(), goodPtResolution(), packLinks(), setDebug(), and ~PFDisplacedVertexCandidateFinder().

double PFDisplacedVertexCandidateFinder::dxy_

private

Definition at line 151 of file PFDisplacedVertexCandidateFinder.h.

Referenced by goodPtResolution(), operator<<(), and setParameters().

std::list< reco::TrackBaseRef > PFDisplacedVertexCandidateFinder::eventTracks_

private

The track refs.

Definition at line 123 of file PFDisplacedVertexCandidateFinder.h.

Referenced by associate(), findDisplacedVertexCandidates(), operator<<(), setInput(), and ~PFDisplacedVertexCandidateFinder().

std::vector< GlobalTrajectoryParameters > PFDisplacedVertexCandidateFinder::eventTrackTrajectories_

private

The Trajectories vector allow to calculate snd to store only once the track trajectory parameters

Definition at line 132 of file PFDisplacedVertexCandidateFinder.h.

Referenced by link(), and setInput().

const MagneticField* PFDisplacedVertexCandidateFinder::magField_

private

Definition at line 166 of file PFDisplacedVertexCandidateFinder.h.

Referenced by getGlobalTrajectoryParameters(), and setInput().

double PFDisplacedVertexCandidateFinder::nChi2_max_

private

Tracks preselection to reduce the combinatorics in PFDisplacedVertexCandidates this cuts are repeated then in a smarter way in the PFDisplacedVertexFinder be sure you are consistent between them

Definition at line 147 of file PFDisplacedVertexCandidateFinder.h.

Referenced by goodPtResolution(), operator<<(), and setParameters().

double PFDisplacedVertexCandidateFinder::primaryVertexCut2_

private

Do not reconstruct vertices wich are too close to the beam pipe.

Definition at line 140 of file PFDisplacedVertexCandidateFinder.h.

Referenced by link(), and setParameters().

double PFDisplacedVertexCandidateFinder::pt_min_

private

Definition at line 148 of file PFDisplacedVertexCandidateFinder.h.

Referenced by goodPtResolution(), operator<<(), and setParameters().

double PFDisplacedVertexCandidateFinder::pt_min_prim_