reco::PFDisplacedVertexSeed Class Reference

Block of elements. More...

#include <PFDisplacedVertexSeed.h>

Public Member Functions
void	addElement (TrackBaseRef)
void	Dump (std::ostream &out=std::cout) const
	cout function More...
const std::vector< TrackBaseRef > &	elements () const
bool	isEmpty () const
	Check if it is a new Seed. More...
void	mergeWith (const PFDisplacedVertexSeed &displacedVertex)
	Merge two Seeds if their seed Points are close enough. More...
const double	nTracks () const
	PFDisplacedVertexSeed ()
	Default constructor. More...
void	reserveElements (size_t)
	Reserve space for elements. More...
const GlobalPoint &	seedPoint () const
const double	totalWeight () const
void	updateSeedPoint (const GlobalPoint &dcaPoint, const TrackBaseRef, const TrackBaseRef, double weight=1)

Private Attributes
std::vector< TrackBaseRef >	elements_
	---—— MEMBERS ----—— /// More...
GlobalPoint	seedPoint_
	Seed point which indicated the approximative position of the vertex. More...
float	totalWeight_

Friends
std::ostream &	operator<< (std::ostream &out, const PFDisplacedVertexSeed &co)

Detailed Description

Block of elements.

Author

Gouzevitch Maxime

Date

November 2009

A DisplacedVertexSeed is an intermediate format, usually not persistent, used by PFDisplacedVertexFinder to keep the information for vertex fit. This format is produced after parsing of PFDisplacedVertexCandidate which by construction may contain many seeds. This format contains:

• a set of Refs to tracks.
• a Point which indicated the approximative position of the vertex.

Definition at line 25 of file PFDisplacedVertexSeed.h.

Constructor & Destructor Documentation

PFDisplacedVertexSeed()

PFDisplacedVertexSeed::PFDisplacedVertexSeed

(

)

Default constructor.

Definition at line 9 of file PFDisplacedVertexSeed.cc.

: seedPoint_(GlobalPoint(0, 0, 0)), totalWeight_(0) {}

Member Function Documentation

addElement()

void PFDisplacedVertexSeed::addElement

(

TrackBaseRef

element

)

Add a track Reference to the current Seed If the track reference is already in the collection, it is ignored

Definition at line 11 of file PFDisplacedVertexSeed.cc.

References elements_, spr::find(), and eostools::move().

Referenced by updateSeedPoint().

                                                           {
  if (std::find(elements_.begin(), elements_.end(), element) == elements_.end()) {
    elements_.emplace_back(std::move(element));
  }
}

Dump()

void PFDisplacedVertexSeed::Dump

(

std::ostream &

out = std::cout

)

const

cout function

Definition at line 63 of file PFDisplacedVertexSeed.cc.

References elements_, MillePedeFileConverter_cfg::out, Pi, DiDispStaMuonMonitor_cfi::pt, seedPoint(), and mathSSE::sqrt().

Referenced by PFDisplacedVertexFinder::fitVertexFromSeed().

                                                   {
  if (!out)
    return;

  out << "\t--- DisplacedVertexSeed ---  " << endl;
  out << "\tnumber of elements: " << elements_.size() << endl;

  out << "\t Seed Point x = " << seedPoint().x() << "\t Seed Point y = " << seedPoint().y()
      << "\t Seed Point z = " << seedPoint().z() << 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 (auto const& ie : elements_) {
    math::XYZPoint Pi(ie.get()->innerPosition());
    math::XYZPoint Po(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 = ie->pt();

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

    out << "ie = " << ie.key() << " pt = "
        << pt
        //       <<" inn hit pos x = " << Pi.x() << " y = " << Pi.y() << " z = " << Pi.z()
        << " out hit pos x = " << Po.x() << " y = " << Po.y() << " z = " << Po.z() << endl;
  }

  out << endl;
}

elements()

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

inline

Returns

vector of unique references to tracks

Definition at line 48 of file PFDisplacedVertexSeed.h.

References elements_.

Referenced by PFDisplacedVertexFinder::fitVertexFromSeed(), and mergeWith().

{ return elements_; }

isEmpty()

bool reco::PFDisplacedVertexSeed::isEmpty ( ) const

inline

Check if it is a new Seed.

Definition at line 45 of file PFDisplacedVertexSeed.h.

References elements_.

Referenced by plotting.Plot::clone(), and updateSeedPoint().

{ return (elements_.empty()); }

mergeWith()

void PFDisplacedVertexSeed::mergeWith

(

const PFDisplacedVertexSeed &

displacedVertex

)

Merge two Seeds if their seed Points are close enough.

Definition at line 41 of file PFDisplacedVertexSeed.cc.

References MillePedeFileConverter_cfg::e, elements(), elements_, spr::find(), reserveElements(), seedPoint(), seedPoint_, totalWeight(), totalWeight_, mps_merge::weight, Basic3DVector< T >::x(), Basic3DVector< T >::y(), and Basic3DVector< T >::z().

                                                                                  {
  double weight = displacedVertex.totalWeight();
  const GlobalPoint& dcaPoint = displacedVertex.seedPoint();

  Basic3DVector<double> vertexSeedVector(seedPoint_);
  Basic3DVector<double> dcaVector(dcaPoint);

  dcaVector = (dcaVector * weight + vertexSeedVector * totalWeight_) / (totalWeight_ + weight);
  GlobalPoint P(dcaVector.x(), dcaVector.y(), dcaVector.z());
  totalWeight_ += weight;
  seedPoint_ = P;

  reserveElements(elements_.size() + displacedVertex.elements().size());
  auto const oldSize = elements_.size();
  //avoid checking elements we just added from displacedVertex.elements()
  for (auto const& e : displacedVertex.elements()) {
    if (std::find(elements_.begin(), elements_.begin() + oldSize, e) == elements_.begin() + oldSize) {
      elements_.emplace_back(e);
    }
  }
}

nTracks()

const double reco::PFDisplacedVertexSeed::nTracks ( ) const

inline

Definition at line 50 of file PFDisplacedVertexSeed.h.

References elements_.

{ return elements_.size(); }

reserveElements()

void PFDisplacedVertexSeed::reserveElements

(

size_t

newSize

)

Reserve space for elements.

Definition at line 17 of file PFDisplacedVertexSeed.cc.

References elements_.

Referenced by mergeWith(), and updateSeedPoint().

{ elements_.reserve(newSize); }

seedPoint()

const GlobalPoint& reco::PFDisplacedVertexSeed::seedPoint ( ) const

inline

Returns

the seedPoint for the vertex fitting

Definition at line 53 of file PFDisplacedVertexSeed.h.

References seedPoint_.

Referenced by Dump(), PFDisplacedVertexFinder::fitVertexFromSeed(), PFDisplacedVertexFinder::isCloseTo(), and mergeWith().

{ return seedPoint_; }

totalWeight()

const double reco::PFDisplacedVertexSeed::totalWeight ( ) const

inline

Returns

the total weight

Definition at line 56 of file PFDisplacedVertexSeed.h.

References totalWeight_.

Referenced by mergeWith().

{ return totalWeight_; }

updateSeedPoint()

void PFDisplacedVertexSeed::updateSeedPoint	(	const GlobalPoint &	dcaPoint,
		const TrackBaseRef	r1,
		const TrackBaseRef	r2,
		double	weight = 1
	)

Add a track Ref to the Seed and recalculate the seedPoint with a new dcaPoint A weight different from 1 may be assign to the new DCA point

Definition at line 19 of file PFDisplacedVertexSeed.cc.

References addElement(), elements_, isEmpty(), diffTwoXMLs::r1, diffTwoXMLs::r2, reserveElements(), seedPoint_, totalWeight_, mps_merge::weight, Basic3DVector< T >::x(), Basic3DVector< T >::y(), and Basic3DVector< T >::z().

                                                           {
  if (isEmpty()) {
    seedPoint_ = dcaPoint;
    totalWeight_ = weight;
  } else {
    Basic3DVector<double> vertexSeedVector(seedPoint_);
    Basic3DVector<double> dcaVector(dcaPoint);

    dcaVector = (dcaVector * weight + vertexSeedVector * totalWeight_) / (totalWeight_ + weight);
    GlobalPoint P(dcaVector.x(), dcaVector.y(), dcaVector.z());
    totalWeight_ += weight;
    seedPoint_ = P;
  }

  reserveElements(elements_.size() + 2);
  addElement(r1);
  addElement(r2);
}

Friends And Related Function Documentation

operator<<

std::ostream& operator<< ( std::ostream &  out, const PFDisplacedVertexSeed &  co  )

friend

Member Data Documentation

elements_

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

private

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

Set of tracks refs associated to the seed

Definition at line 67 of file PFDisplacedVertexSeed.h.

Referenced by addElement(), Dump(), elements(), isEmpty(), mergeWith(), nTracks(), reserveElements(), and updateSeedPoint().

seedPoint_

GlobalPoint reco::PFDisplacedVertexSeed::seedPoint_

private

Seed point which indicated the approximative position of the vertex.

Definition at line 69 of file PFDisplacedVertexSeed.h.

Referenced by mergeWith(), seedPoint(), and updateSeedPoint().

totalWeight_

float reco::PFDisplacedVertexSeed::totalWeight_

private

Total weight of the points used to calculate the seed point. Necessary for UpdateSeed Point function

Definition at line 72 of file PFDisplacedVertexSeed.h.

Referenced by mergeWith(), totalWeight(), and updateSeedPoint().