reco::PFDisplacedVertexCandidate Class Reference

A block of tracks linked together. More...

#include <PFDisplacedVertexCandidate.h>

Classes
struct	VertexLink

Public Types
typedef std::map< float, std::pair< int, int > >	DistMap
typedef std::vector< float >	DistVector
typedef std::map< unsigned int, VertexLink >	VertexLinkData
enum	VertexLinkTest { LINKTEST_DCA, LINKTEST_DUMMY, LINKTEST_ALL }

Public Member Functions
void	addElement (const TrackBaseRef)
	add a track Reference to the current Candidate More...
void	associatedElements (const unsigned i, const VertexLinkData &vertexLinkData, std::multimap< float, unsigned > &sortedAssociates, const VertexLinkTest test=LINKTEST_DCA) const
	associate 2 elements More...
const GlobalPoint	dcaPoint (unsigned ie1, unsigned ie2) const
DistVector	distVector () const
void	Dump (std::ostream &out=std::cout) const
	cout function More...
const std::vector< TrackBaseRef > &	elements () const
bool	isValid () const
	A Vertex Candidate is valid if it has at least two tracks. More...
unsigned	nTracks () const
	PFDisplacedVertexCandidate ()
	Default constructor. More...
DistMap	r2Map () const
	--—— Provide useful information --—— /// More...
DistVector	r2Vector () const
void	setLink (unsigned i1, unsigned i2, const float dist, const GlobalPoint &dcaPoint, const VertexLinkTest test=LINKTEST_DCA)
const TrackBaseRef &	tref (unsigned ie) const
const VertexLinkData &	vertexLinkData () const

Private Member Functions
const float	dist (unsigned ie1, unsigned ie2) const
	--—— Internal tools --—— /// More...
bool	matrix2vector (unsigned i, unsigned j, unsigned &index) const
bool	testLink (unsigned ie1, unsigned ie2) const
	test if a link between two tracks is valid: value_link =! -1 More...
unsigned	vertexLinkDataSize () const
	--—— Storage of the information --—— /// More...

Private Attributes
std::vector< TrackBaseRef >	elements_
	--—— MEMBERS --—— /// More...
VertexLinkData	vertexLinkData_
	map of links between tracks More...

Friends
std::ostream &	operator<< (std::ostream &, const PFDisplacedVertexCandidate &)
	cout function More...

Detailed Description

A block of tracks linked together.

Author

Gouzevitch Maxime

Date

November 2009

A DisplacedVertexCandidate is a format produced by the DisplacedVertexCandidateFinder to hold a Collection of Refs to the tracks linked together by DCA. It contains:

• a set of reco::Tracks candidates for a common displaced vertex.
• a set of links between these reco::Tracks.
• a Seed point calculated as average from DCA Points. ie. points Pm in the middle between two points P1 and P2 belonging to tracks 1 and 2: P1–Pm–P2

Definition at line 28 of file PFDisplacedVertexCandidate.h.

Member Typedef Documentation

DistMap

typedef std::map<float, std::pair<int, int> > reco::PFDisplacedVertexCandidate::DistMap

A type to provide the information about the position of DCA Points or values of DCA.

Definition at line 49 of file PFDisplacedVertexCandidate.h.

DistVector

typedef std::vector<float> reco::PFDisplacedVertexCandidate::DistVector

Definition at line 50 of file PFDisplacedVertexCandidate.h.

VertexLinkData

typedef std::map<unsigned int, VertexLink> reco::PFDisplacedVertexCandidate::VertexLinkData

Definition at line 45 of file PFDisplacedVertexCandidate.h.

Member Enumeration Documentation

VertexLinkTest

enum reco::PFDisplacedVertexCandidate::VertexLinkTest

Test used for the track linkind. For the moment only DCA is used, but other are possibles like distance between inner hits.

Enumerator
LINKTEST_DCA
LINKTEST_DUMMY
LINKTEST_ALL

Definition at line 43 of file PFDisplacedVertexCandidate.h.

{ LINKTEST_DCA, LINKTEST_DUMMY, LINKTEST_ALL };

Constructor & Destructor Documentation

PFDisplacedVertexCandidate()

PFDisplacedVertexCandidate::PFDisplacedVertexCandidate

(

)

Default constructor.

Definition at line 9 of file PFDisplacedVertexCandidate.cc.

{}

Member Function Documentation

addElement()

void PFDisplacedVertexCandidate::addElement

(

const TrackBaseRef

element

)

add a track Reference to the current Candidate

Definition at line 11 of file PFDisplacedVertexCandidate.cc.

Referenced by PFDisplacedVertexCandidateFinder::associate().

{ elements_.push_back(element); }

associatedElements()

void PFDisplacedVertexCandidate::associatedElements	(	const unsigned	i,
		const VertexLinkData &	vertexLinkData,
		std::multimap< float, unsigned > &	sortedAssociates,
		const VertexLinkTest	test = LINKTEST_DCA
	)		const

associate 2 elements

Definition at line 39 of file PFDisplacedVertexCandidate.cc.

References mps_fire::i, ALPAKA_ACCELERATOR_NAMESPACE::caPixelDoublets::if(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, and builder_last_value_cfg::test.

                                                                                     {
  sortedAssociates.clear();

  // i is too large
  if (i > elements_.size())
    return;
  // assert(i>=0); // i >= 0, since i is unsigned
  for (unsigned ie = 0; ie < elements_.size(); ie++) {
    // considered element itself
    if (ie == i) {
      continue;
    }

    // Order the elements by increasing distance !

    unsigned index = 0;
    if (!matrix2vector(i, ie, index))
      continue;

    float c2 = -1;
    VertexLinkData::const_iterator it = vertexLinkData.find(index);
    if (it != vertexLinkData.end() && (((1 << test) & it->second.test_) != 0 || (test == LINKTEST_ALL)))
      c2 = it->second.distance_;

    // not associated
    if (c2 < 0) {
      continue;
    }

    sortedAssociates.insert(pair<float, unsigned>(c2, ie));
  }
}

dcaPoint()

const GlobalPoint PFDisplacedVertexCandidate::dcaPoint	(	unsigned	ie1,
		unsigned	ie2
	)		const

Returns

DCA point between two tracks

Definition at line 126 of file PFDisplacedVertexCandidate.cc.

References ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it.

Referenced by PFDisplacedVertexFinder::findSeedsFromCandidate().

                                                                                       {
  GlobalPoint dcaPoint(1e10, 1e10, 1e10);

  unsigned index = 0;
  if (!matrix2vector(ie1, ie2, index))
    return dcaPoint;
  VertexLinkData::const_iterator it = vertexLinkData_.find(index);
  if (it != vertexLinkData_.end())
    dcaPoint = it->second.dcaPoint_;

  return dcaPoint;
}

dist()

const float PFDisplacedVertexCandidate::dist ( unsigned  ie1, unsigned  ie2  ) const

private

--—— Internal tools --—— ///

Returns

distance of link between two tracks

Definition at line 148 of file PFDisplacedVertexCandidate.cc.

References ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it.

                                                                             {
  float dist = -1;

  unsigned index = 0;
  if (!matrix2vector(ie1, ie2, index))
    return dist;
  VertexLinkData::const_iterator it = vertexLinkData_.find(index);
  if (it != vertexLinkData_.end())
    dist = it->second.distance_;

  return dist;
}

distVector()

PFDisplacedVertexCandidate::DistVector PFDisplacedVertexCandidate::distVector

(

)

const

Returns

the vector of DCA useful for DCA

Definition at line 111 of file PFDisplacedVertexCandidate.cc.

References ztail::d.

                                                                                  {
  DistVector distVector;

  for (unsigned ie1 = 0; ie1 < elements_.size(); ie1++)
    for (unsigned ie2 = ie1 + 1; ie2 < elements_.size(); ie2++) {
      float d = dist(ie1, ie2);
      if (d < -0.5)
        continue;

      distVector.push_back(d);
    }

  return distVector;
}

Dump()

void PFDisplacedVertexCandidate::Dump

(

std::ostream &

out = std::cout

)

const

cout function

Definition at line 192 of file PFDisplacedVertexCandidate.cc.

References bookConverter::elements, MillePedeFileConverter_cfg::out, Pi, DiDispStaMuonMonitor_cfi::pt, and mathSSE::sqrt().

                                                        {
  if (!out)
    return;

  const vector<TrackBaseRef>& elements = elements_;
  out << "\t--- DisplacedVertexCandidate ---  " << endl;
  out << "\tnumber of elements: " << elements.size() << endl;

  // Build element label (string) : elid from type, layer and occurence number
  // use stringstream instead of sprintf to concatenate string and integer into string
  for (unsigned ie = 0; ie < elements.size(); ie++) {
    math::XYZPoint Pi(elements[ie].get()->innerPosition());
    math::XYZPoint Po(elements[ie].get()->outerPosition());

    float innermost_radius = sqrt(Pi.x() * Pi.x() + Pi.y() * Pi.y() + Pi.z() * Pi.z());
    float outermost_radius = sqrt(Po.x() * Po.x() + Po.y() * Po.y() + Po.z() * Po.z());
    float innermost_rho = sqrt(Pi.x() * Pi.x() + Pi.y() * Pi.y());
    float outermost_rho = sqrt(Po.x() * Po.x() + Po.y() * Po.y());

    double pt = elements[ie]->pt();

    out << "ie = " << elements[ie].key() << " pt = " << pt << " innermost hit radius = " << innermost_radius
        << " rho = " << innermost_rho << " outermost hit radius = " << outermost_radius << " rho = " << outermost_rho
        << endl;
  }

  out << endl;
}

elements()

const std::vector<TrackBaseRef>& reco::PFDisplacedVertexCandidate::elements ( ) const

inline

Returns

the vector of Refs to tracks

Definition at line 95 of file PFDisplacedVertexCandidate.h.

References elements_.

Referenced by PFDisplacedVertexCandidateFinder::associate(), PFDisplacedVertexCandidateFinder::findDisplacedVertexCandidates(), and PFDisplacedVertexCandidateFinder::packLinks().

{ return elements_; }

isValid()

bool reco::PFDisplacedVertexCandidate::isValid ( void  ) const

inline

A Vertex Candidate is valid if it has at least two tracks.

Definition at line 89 of file PFDisplacedVertexCandidate.h.

References elements_.

Referenced by ntupleDataFormat._Object::_checkIsValid(), PFDisplacedVertexCandidateFinder::findDisplacedVertexCandidates(), and core.AutoHandle.AutoHandle::ReallyLoad().

{ return elements_.size() > 1; }

matrix2vector()

bool PFDisplacedVertexCandidate::matrix2vector ( unsigned  i, unsigned  j, unsigned &  index  ) const

private

makes the correspondance between a 2d element matrix and the 1D vector which is the most compact way to store the matrix

Definition at line 172 of file PFDisplacedVertexCandidate.cc.

References combine::missing, and edm::swap().

                                                                                                      {
  unsigned size = elements_.size();
  if (iindex == jindex || iindex >= size || jindex >= size) {
    return false;
  }

  if (iindex > jindex)
    swap(iindex, jindex);

  index = jindex - iindex - 1;

  if (iindex > 0) {
    index += iindex * size;
    unsigned missing = iindex * (iindex + 1) / 2;
    index -= missing;
  }

  return true;
}

nTracks()

unsigned reco::PFDisplacedVertexCandidate::nTracks ( ) const

inline

Returns

the number of tracks associated to the candidate

Definition at line 98 of file PFDisplacedVertexCandidate.h.

References elements_.

{ return elements_.size(); }

r2Map()

PFDisplacedVertexCandidate::DistMap PFDisplacedVertexCandidate::r2Map

(

)

const

--—— Provide useful information --—— ///

Returns

the map of Radius^2 to DCA Points

Definition at line 77 of file PFDisplacedVertexCandidate.cc.

References diffTwoXMLs::r2.

Referenced by PFDisplacedVertexFinder::findSeedsFromCandidate().

                                                                          {
  DistMap r2Map;

  for (unsigned ie1 = 0; ie1 < elements_.size(); ie1++)
    for (unsigned ie2 = ie1 + 1; ie2 < elements_.size(); ie2++) {
      GlobalPoint P = dcaPoint(ie1, ie2);
      if (P.x() > 1e9)
        continue;

      float r2 = P.x() * P.x() + P.y() * P.y() + P.z() * P.z();

      r2Map.insert(pair<float, pair<int, int> >(r2, pair<int, int>(ie1, ie2)));
    }

  return r2Map;
}

r2Vector()

PFDisplacedVertexCandidate::DistVector PFDisplacedVertexCandidate::r2Vector

(

)

const

Returns

the vector of Radius^2 to DCA Points useful for FWLite

Definition at line 94 of file PFDisplacedVertexCandidate.cc.

References diffTwoXMLs::r2.

                                                                                {
  DistVector r2Vector;

  for (unsigned ie1 = 0; ie1 < elements_.size(); ie1++)
    for (unsigned ie2 = ie1 + 1; ie2 < elements_.size(); ie2++) {
      GlobalPoint P = dcaPoint(ie1, ie2);
      if (P.x() > 1e9)
        continue;

      float r2 = P.x() * P.x() + P.y() * P.y() + P.z() * P.z();

      r2Vector.push_back(r2);
    }

  return r2Vector;
}

setLink()

void PFDisplacedVertexCandidate::setLink	(	unsigned	i1,
		unsigned	i2,
		const float	dist,
		const GlobalPoint &	dcaPoint,
		const VertexLinkTest	test = LINKTEST_DCA
	)

set a link between elements of indices i1 and i2, of "distance" dist the link is set in the linkData vector provided as an argument. As indicated by the 'const' statement, 'this' is not modified.

Definition at line 13 of file PFDisplacedVertexCandidate.cc.

References cms::cuda::assert(), testProducerWithPsetDescEmpty_cfi::i1, testProducerWithPsetDescEmpty_cfi::i2, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, MainPageGenerator::l, convertSQLiteXML::ok, and builder_last_value_cfg::test.

Referenced by PFDisplacedVertexCandidateFinder::packLinks().

                                                                                                                    {
  assert(test < LINKTEST_ALL);

  unsigned index = 0;
  bool ok = matrix2vector(i1, i2, index);

  if (ok) {
    //ignore the  -1, -1 pair
    if (dist > -0.5) {
      VertexLink& l = vertexLinkData_[index];
      l.distance_ = dist;
      l.dcaPoint_ = dcaPoint;
      l.test_ |= (1 << test);
    } else  //delete if existing
    {
      VertexLinkData::iterator it = vertexLinkData_.find(index);
      if (it != vertexLinkData_.end())
        vertexLinkData_.erase(it);
    }

  } else {
    assert(0);
  }
}

testLink()

bool PFDisplacedVertexCandidate::testLink ( unsigned  ie1, unsigned  ie2  ) const

private

test if a link between two tracks is valid: value_link =! -1

Definition at line 141 of file PFDisplacedVertexCandidate.cc.

References ztail::d.

                                                                          {
  float d = dist(ie1, ie2);
  if (d < -0.5)
    return false;
  return true;
}

tref()

const TrackBaseRef& reco::PFDisplacedVertexCandidate::tref ( unsigned  ie ) const

inline

Returns

the reference to a given tracks

Definition at line 92 of file PFDisplacedVertexCandidate.h.

References elements_.

Referenced by PFDisplacedVertexFinder::findSeedsFromCandidate().

{ return elements_[ie]; }

vertexLinkData()

const VertexLinkData& reco::PFDisplacedVertexCandidate::vertexLinkData ( ) const

inline

Returns

the map of link data

Definition at line 101 of file PFDisplacedVertexCandidate.h.

References vertexLinkData_.

{ return vertexLinkData_; }

vertexLinkDataSize()

unsigned PFDisplacedVertexCandidate::vertexLinkDataSize ( ) const

private

--—— Storage of the information --—— ///

Those are the tools from PFBlockAlgo

Returns

size of linkData_, calculated from the number of elements

Definition at line 163 of file PFDisplacedVertexCandidate.cc.

References create_idmaps::n.

                                                              {
  unsigned n = elements_.size();

  // number of possible undirected links between n elements.
  // reflective links impossible.

  return n * (n - 1) / 2;
}

Friends And Related Function Documentation

operator<<

std::ostream& operator<< ( std::ostream &  , const PFDisplacedVertexCandidate &    )

friend

cout function

Member Data Documentation

elements_

std::vector<TrackBaseRef> reco::PFDisplacedVertexCandidate::elements_

private

--—— MEMBERS --—— ///

vector of refs to the associated tracks

Definition at line 131 of file PFDisplacedVertexCandidate.h.

Referenced by elements(), isValid(), nTracks(), and tref().

vertexLinkData_

VertexLinkData reco::PFDisplacedVertexCandidate::vertexLinkData_

private

map of links between tracks

Definition at line 134 of file PFDisplacedVertexCandidate.h.

Referenced by vertexLinkData().