PFDisplacedVertexCandidateFinder Class Reference

#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 &crossing_point, 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_
double	qoverpError_max_
TwoTrackMinimumDistance	theMinimum_
edm::soa::PtEtaPhiTable	track_table_
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 26 of file PFDisplacedVertexCandidateFinder.h.

Member Typedef Documentation

IBC

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

Definition at line 37 of file PFDisplacedVertexCandidateFinder.h.

IE

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

Definition at line 35 of file PFDisplacedVertexCandidateFinder.h.

IEC

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

Definition at line 36 of file PFDisplacedVertexCandidateFinder.h.

Mask

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

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

Definition at line 33 of file PFDisplacedVertexCandidateFinder.h.

Constructor & Destructor Documentation

PFDisplacedVertexCandidateFinder()

PFDisplacedVertexCandidateFinder::PFDisplacedVertexCandidateFinder

(

)

Definition at line 18 of file PFDisplacedVertexCandidateFinder.cc.

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

~PFDisplacedVertexCandidateFinder()

PFDisplacedVertexCandidateFinder::~PFDisplacedVertexCandidateFinder

(

)

Definition at line 28 of file PFDisplacedVertexCandidateFinder.cc.

References eventTracks_, and LogDebug.

                                                                    {
  LogDebug("PFDisplacedVertexCandidateFinder")
      << "~PFDisplacedVertexCandidateFinder - number of remaining elements: " << eventTracks_.size();
}

Member Function Documentation

associate()

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 111 of file PFDisplacedVertexCandidateFinder.cc.

References reco::PFDisplacedVertexCandidate::addElement(), reco::PFDisplacedVertexCandidate::elements(), eventTracks_, mps_fire::i, dqmdumpme::last, link(), LogDebug, LogTrace, GetRecoTauVFromDQM_MC_cff::next, trackMask_, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by findDisplacedVertexCandidates().

                                                                       {
  LogDebug("PFDisplacedVertexCandidateFinder") << "== Start the association procedure ==";

  if (last != eventTracks_.end()) {
    double dist = -1;
    GlobalPoint crossingPoint(0, 0, 0);
    PFDisplacedVertexCandidate::VertexLinkTest linktest;
    link(*last, *next, dist, crossingPoint, linktest);

    if (dist < -0.5) {
      LogDebug("PFDisplacedVertexCandidateFinder") << "link failed";

      return ++next;  // association failed
    } else {
      // add next element to the current pflowblock
      LogDebug("PFDisplacedVertexCandidateFinder")
          << "adding element " << (*next).key() << " to PFDisplacedVertexCandidate with "
          << tempVertexCandidate.elements().size() << " elements";
      tempVertexCandidate.addElement((*next));
      trackMask_[(*next).key()] = false;

      LogDebug("PFDisplacedVertexCandidateFinder")
          << "link parameters "
          << " *next = " << (*next).key() << " *last = " << (*last).key() << "  dist = " << crossingPoint
          << " P.x = " << crossingPoint.x() << " P.y = " << crossingPoint.y() << " P.z = " << crossingPoint.z();
    }
  } else {
    // add next element to this eflowblock
    LogDebug("PFDisplacedVertexCandidateFinder") << "adding to block element " << (*next).key();
    tempVertexCandidate.addElement((*next));
    trackMask_[(*next).key()] = false;
  }

  // recursive call: associate next and other unused elements
#ifdef EDM_ML_DEBUG
  for (unsigned i = 0; i < trackMask_.size(); i++)
    LogTrace("PFDisplacedVertexCandidateFinder") << " Mask[" << i << "] = " << trackMask_[i];
#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;
    }

    LogDebug("PFDisplacedVertexCandidateFinder") << "calling associate " << (*next).key() << " & " << (*ie).key();
    ie = associate(next, ie, tempVertexCandidate);
  }

  LogDebug("PFDisplacedVertexCandidateFinder") << "**** removing element ";

  IE iteratorToNextFreeElement = eventTracks_.erase(next);

  LogDebug("PFDisplacedVertexCandidateFinder") << "== End the association procedure ==";

  return iteratorToNextFreeElement;
}

findDisplacedVertexCandidates()

void PFDisplacedVertexCandidateFinder::findDisplacedVertexCandidates

(

)

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

Definition at line 82 of file PFDisplacedVertexCandidateFinder.cc.

References associate(), reco::PFDisplacedVertexCandidate::elements(), eventTracks_, reco::PFDisplacedVertexCandidate::isValid(), LogDebug, packLinks(), vertexCandidates_, and vertexCandidatesSize_.

Referenced by PFDisplacedVertexCandidateProducer::produce().

                                                                     {
  LogDebug("PFDisplacedVertexCandidateFinder") << "========= Start Finding Displaced Vertex Candidates =========";

  // The vertexCandidates have not been passed to the event, and need to be cleared
  if (vertexCandidates_.get())
    vertexCandidates_->clear();
  else
    vertexCandidates_ = std::make_unique<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);
    LogDebug("PFDisplacedVertexCandidateFinder") << "elements=" << tempVertexCandidate.elements().size();

    // Build remaining links in current block
    if (tempVertexCandidate.isValid()) {
      packLinks(tempVertexCandidate);
      vertexCandidates_->push_back(tempVertexCandidate);
    }
  }

  LogDebug("PFDisplacedVertexCandidateFinder") << "========= End Finding Displaced Vertex Candidates =========";
}

getGlobalTrajectoryParameters()

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 299 of file PFDisplacedVertexCandidateFinder.cc.

References magField_, position, and HLT_2024v14_cff::track.

Referenced by setInput().

                                                                                                                   {
  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;
}

goodPtResolution()

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

private

Quality Criterion on the Pt resolution to select a Track.

Definition at line 310 of file PFDisplacedVertexCandidateFinder.cc.

References funct::abs(), PVValHelper::dxy, reco::TrackBase::dxy(), dxy_, LogDebug, nChi2_max_, reco::TrackBase::normalizedChi2(), DiDispStaMuonMonitor_cfi::pt, reco::TrackBase::pt(), pt_min_prim_, reco::TrackBase::ptError(), pvtx_, run3scouting_cff::qoverpError, reco::TrackBase::qoverpError(), and qoverpError_max_.

Referenced by isSelected().

                                                                                          {
  double nChi2 = trackref->normalizedChi2();
  double pt = trackref->pt();
  double dpt = trackref->ptError();
  double dxy = trackref->dxy(pvtx_);

  double pt_error = dpt / pt * 100;
  double qoverpError = trackref->qoverpError();

  LogDebug("PFDisplacedVertexCandidateFinder")
      << " PFDisplacedVertexFinder: PFrecTrack->Track Pt= " << pt << " dPt/Pt = " << pt_error << "% nChi2 = " << nChi2;
  if (nChi2 > nChi2_max_ || pt < pt_min_ || qoverpError > qoverpError_max_) {
    LogDebug("PFDisplacedVertexCandidateFinder") << " PFBlockAlgo: skip badly measured or low pt track"
                                                 << " nChi2_cut = " << 5 << " pt_cut = " << 0.2;
    return false;
  }
  //  cout << "dxy = " << dxy << endl;
  if (std::abs(dxy) < dxy_ && pt < pt_min_prim_)
    return false;
  //  if (fabs(dxy) < 0.2 && pt < 0.8) return false;

  return true;
}

isSelected()

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

inlineprivate

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

Definition at line 104 of file PFDisplacedVertexCandidateFinder.h.

References goodPtResolution().

Referenced by setInput().

{ return goodPtResolution(trackref); }

link()

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

private

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

Definition at line 176 of file PFDisplacedVertexCandidateFinder.cc.

References funct::abs(), TwoTrackMinimumDistance::calculate(), TwoTrackMinimumDistance::crossingPoint(), dcaCut_, SiPixelRawToDigiRegional_cfi::deltaPhi, TwoTrackMinimumDistance::distance(), HLT_2024v14_cff::eta1, HLT_2024v14_cff::eta2, eventTrackTrajectories_, edm::RefToBase< T >::key(), LogDebug, primaryVertexCut2_, HLT_2024v14_cff::pt1, HLT_2024v14_cff::pt2, theMinimum_, track_table_, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by associate(), and packLinks().

                                                                                                      {
  using namespace edm::soa::col;
  const auto iel1 = el1.key();
  const auto iel2 = el2.key();
  const auto pt1 = track_table_.get<Pt>(iel1);
  const auto pt2 = track_table_.get<Pt>(iel2);
  const auto eta1 = track_table_.get<Eta>(iel1);
  const auto eta2 = track_table_.get<Eta>(iel2);
  const auto phi1 = track_table_.get<Phi>(iel1);
  const auto phi2 = track_table_.get<Phi>(iel2);

  if (std::abs(eta1 - eta2) > 1) {
    dist = -1;
    return;
  }
  if (pt1 > 2 && pt2 > 2 && std::abs(::deltaPhi(phi1, phi2)) > 1) {
    dist = -1;
    return;
  }

  const GlobalTrajectoryParameters& gt1 = eventTrackTrajectories_[iel1];
  const GlobalTrajectoryParameters& gt2 = eventTrackTrajectories_[iel2];

  // Closest approach
  theMinimum_.calculate(gt1, gt2);

  // Fill the parameters
  dist = theMinimum_.distance();
  LogDebug("PFDisplacedVertexCandidateFinder") << "link iel1=" << iel1 << " iel2=" << iel2 << " dist=" << dist;
  crossingPoint = 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 = crossingPoint.x() * crossingPoint.x() + crossingPoint.y() * crossingPoint.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 (std::abs(crossingPoint.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;
  }
  */
}

packLinks()

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

private

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

Definition at line 271 of file PFDisplacedVertexCandidateFinder.cc.

References reco::PFDisplacedVertexCandidate::elements(), testProducerWithPsetDescEmpty_cfi::i1, testProducerWithPsetDescEmpty_cfi::i2, link(), LogDebug, and reco::PFDisplacedVertexCandidate::setLink().

Referenced by findDisplacedVertexCandidates().

                                                                                            {
  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 crossingPoint(0, 0, 0);
      PFDisplacedVertexCandidate::VertexLinkTest linktest;

      link(els[i1], els[i2], dist, crossingPoint, linktest);
      LogDebug("PFDisplacedVertexCandidateFinder")
          << "Setting link between elements " << i1 << " key " << els[i1].key() << " and " << i2 << " key "
          << els[i2].key() << " of dist =" << dist << " computed from link test " << linktest;

      if (dist > -0.5)
        vertexCandidate.setLink(i1, i2, dist, crossingPoint, linktest);
    }
  }
}

setDebug()

void PFDisplacedVertexCandidateFinder::setDebug ( bool  debug )

inline

sets debug printout flag

Definition at line 54 of file PFDisplacedVertexCandidateFinder.h.

References debug, and debug_.

Referenced by PFDisplacedVertexCandidateProducer::PFDisplacedVertexCandidateProducer().

{ debug_ = debug; }

setInput()

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 50 of file PFDisplacedVertexCandidateFinder.cc.

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

Referenced by PFDisplacedVertexCandidateProducer::produce().

                                                                               {
  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);
    }
  }
  track_table_ = edm::soa::PtEtaPhiTable{*trackh};
}

setParameters()

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 42 of file PFDisplacedVertexCandidateFinder.h.

References dcaCut_, dcaPInnerHitCut2_, dxy_, edm::ParameterSet::getParameter(), nChi2_max_, primaryVertexCut2_, pt_min_, pt_min_prim_, and qoverpError_max_.

Referenced by PFDisplacedVertexCandidateProducer::PFDisplacedVertexCandidateProducer().

                                                                                                                    {
    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");
    qoverpError_max_ = ps_trk.getParameter<double>("qoverpError_max");
  }

setPrimaryVertex()

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

Definition at line 33 of file PFDisplacedVertexCandidateFinder.cc.

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

Referenced by PFDisplacedVertexCandidateProducer::produce().

                                                                                                {
  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;
}

transferVertexCandidates()

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

inline

Returns

auto_ptr to collection of DisplacedVertexCandidates

Definition at line 61 of file PFDisplacedVertexCandidateFinder.h.

References eostools::move(), and vertexCandidates_.

Referenced by PFDisplacedVertexCandidateProducer::produce().

                                                                                     {
    return std::move(vertexCandidates_);
  }

vertexCandidates()

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

inline

Definition at line 65 of file PFDisplacedVertexCandidateFinder.h.

References eostools::move(), and vertexCandidates_.

Referenced by operator<<().

                                                                                          {
    return std::move(vertexCandidates_);
  }

Friends And Related Function Documentation

operator<<

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

friend

Definition at line 334 of file PFDisplacedVertexCandidateFinder.cc.

                                                                                {
  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 << " std::abs(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

dcaCut_

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 127 of file PFDisplacedVertexCandidateFinder.h.

Referenced by link(), and setParameters().

dcaPInnerHitCut2_

double PFDisplacedVertexCandidateFinder::dcaPInnerHitCut2_

private

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

Definition at line 131 of file PFDisplacedVertexCandidateFinder.h.

Referenced by setParameters().

debug_

bool PFDisplacedVertexCandidateFinder::debug_

private

if true, debug printouts activated

Definition at line 152 of file PFDisplacedVertexCandidateFinder.h.

Referenced by setDebug().

dxy_

double PFDisplacedVertexCandidateFinder::dxy_

private

Definition at line 140 of file PFDisplacedVertexCandidateFinder.h.

Referenced by goodPtResolution(), and setParameters().

eventTracks_

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

private

The track refs.

Definition at line 113 of file PFDisplacedVertexCandidateFinder.h.

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

eventTrackTrajectories_

std::vector<GlobalTrajectoryParameters> PFDisplacedVertexCandidateFinder::eventTrackTrajectories_

private

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

Definition at line 121 of file PFDisplacedVertexCandidateFinder.h.

Referenced by link(), and setInput().

magField_

const MagneticField* PFDisplacedVertexCandidateFinder::magField_

private

Definition at line 155 of file PFDisplacedVertexCandidateFinder.h.

Referenced by getGlobalTrajectoryParameters(), and setInput().

nChi2_max_

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 136 of file PFDisplacedVertexCandidateFinder.h.

Referenced by goodPtResolution(), and setParameters().

primaryVertexCut2_

double PFDisplacedVertexCandidateFinder::primaryVertexCut2_

private

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

Definition at line 129 of file PFDisplacedVertexCandidateFinder.h.

Referenced by link(), and setParameters().

pt_min_

double PFDisplacedVertexCandidateFinder::pt_min_

private

Definition at line 137 of file PFDisplacedVertexCandidateFinder.h.

Referenced by setParameters().

pt_min_prim_

double PFDisplacedVertexCandidateFinder::pt_min_prim_

private

Definition at line 139 of file PFDisplacedVertexCandidateFinder.h.

Referenced by goodPtResolution(), and setParameters().

pvtx_

math::XYZPoint PFDisplacedVertexCandidateFinder::pvtx_

private

Definition at line 149 of file PFDisplacedVertexCandidateFinder.h.

Referenced by goodPtResolution(), and setPrimaryVertex().

qoverpError_max_

double PFDisplacedVertexCandidateFinder::qoverpError_max_

private

Definition at line 141 of file PFDisplacedVertexCandidateFinder.h.

Referenced by goodPtResolution(), and setParameters().

theMinimum_

TwoTrackMinimumDistance PFDisplacedVertexCandidateFinder::theMinimum_

private

Definition at line 147 of file PFDisplacedVertexCandidateFinder.h.

Referenced by link().

track_table_

edm::soa::PtEtaPhiTable PFDisplacedVertexCandidateFinder::track_table_

private

Definition at line 157 of file PFDisplacedVertexCandidateFinder.h.

Referenced by link(), and setInput().

trackMask_

Mask PFDisplacedVertexCandidateFinder::trackMask_

private

The trackMask allows to keep the information on the tracks which are still free and those which are already used or disabled.

Definition at line 118 of file PFDisplacedVertexCandidateFinder.h.

Referenced by associate(), and setInput().

vertexCandidates_

std::unique_ptr<reco::PFDisplacedVertexCandidateCollection> PFDisplacedVertexCandidateFinder::vertexCandidates_

private

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

Definition at line 110 of file PFDisplacedVertexCandidateFinder.h.

Referenced by findDisplacedVertexCandidates(), transferVertexCandidates(), and vertexCandidates().

vertexCandidatesSize_

unsigned PFDisplacedVertexCandidateFinder::vertexCandidatesSize_

private

Max number of expected vertexCandidates in the event Used to allocate the memory and avoid multiple copy

Definition at line 144 of file PFDisplacedVertexCandidateFinder.h.

Referenced by findDisplacedVertexCandidates().