reco::PFDisplacedVertexSeed Class Reference

Block of elements. More...

#include <PFDisplacedVertexSeed.h>

Classes
struct	Compare

Public Types
typedef std::set< reco::TrackBaseRef >::iterator	IEset
	-----— Useful Types -----— /// More...

Public Member Functions
void	addElement (TrackBaseRef)
	Add a track Reference to the current Seed. More...
void	Dump (std::ostream &out=std::cout) const
	cout function More...
const std::set< TrackBaseRef, Compare > &	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...
const GlobalPoint &	seedPoint () const
const double	totalWeight () const
void	updateSeedPoint (const GlobalPoint &dcaPoint, const TrackBaseRef, const TrackBaseRef, double weight=1)

Private Attributes
std::set< TrackBaseRef, Compare >	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 28 of file PFDisplacedVertexSeed.h.

Member Typedef Documentation

typedef std::set< reco::TrackBaseRef >::iterator reco::PFDisplacedVertexSeed::IEset

-----— Useful Types -----— ///

Definition at line 34 of file PFDisplacedVertexSeed.h.

Constructor & Destructor Documentation

PFDisplacedVertexSeed::PFDisplacedVertexSeed

(

)

Default constructor.

Definition at line 10 of file PFDisplacedVertexSeed.cc.

Referenced by reco::PFDisplacedVertexSeed::Compare::operator()().

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

Member Function Documentation

void PFDisplacedVertexSeed::addElement

(

TrackBaseRef

)

Add a track Reference to the current Seed.

Definition at line 16 of file PFDisplacedVertexSeed.cc.

References elements_.

Referenced by mergeWith(), reco::PFDisplacedVertexSeed::Compare::operator()(), and updateSeedPoint().

                                                           {
  elements_.insert( element ); 
}

void PFDisplacedVertexSeed::Dump

(

std::ostream &

out = std::cout

)

const

cout function

Definition at line 71 of file PFDisplacedVertexSeed.cc.

References elements_, Pi, EnergyCorrector::pt, seedPoint(), and mathSSE::sqrt().

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

                                                     {
  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(IEset ie = elements_.begin(); ie !=  elements_.end(); ie++){

    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;


}

const std::set< TrackBaseRef, Compare >& reco::PFDisplacedVertexSeed::elements ( ) const

inline

Returns

set of references to tracks

Definition at line 64 of file PFDisplacedVertexSeed.h.

References elements_.

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

      {return elements_;}

bool reco::PFDisplacedVertexSeed::isEmpty ( ) const

inline

Check if it is a new Seed.

Definition at line 61 of file PFDisplacedVertexSeed.h.

References elements_.

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

{return (elements_.size() == 0);}

void PFDisplacedVertexSeed::mergeWith

(

const PFDisplacedVertexSeed &

displacedVertex

)

Merge two Seeds if their seed Points are close enough.

Definition at line 46 of file PFDisplacedVertexSeed.cc.

References addElement(), elements(), seedPoint(), seedPoint_, totalWeight(), totalWeight_, mps_merge::weight, Basic3DVector< T >::x(), Basic3DVector< T >::y(), and Basic3DVector< T >::z().

Referenced by reco::PFDisplacedVertexSeed::Compare::operator()().

                                                                                 {

  
  double weight = displacedVertex.totalWeight();
  set<TrackBaseRef, Compare> newElements= displacedVertex.elements();
  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;



  for (  set<TrackBaseRef, Compare>::const_iterator il = newElements.begin(); il != newElements.end(); il++)
    addElement(*il);



}

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

inline

Definition at line 67 of file PFDisplacedVertexSeed.h.

References elements_.

{return elements_.size();}

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

inline

Returns

the seedPoint for the vertex fitting

Definition at line 70 of file PFDisplacedVertexSeed.h.

References seedPoint_.

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

{return seedPoint_;}

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

inline

Returns

the total weight

Definition at line 73 of file PFDisplacedVertexSeed.h.

References gather_cfg::cout, Dump(), operator<<, MillePedeFileConverter_cfg::out, and totalWeight_.

Referenced by mergeWith().

{return totalWeight_;}

void PFDisplacedVertexSeed::updateSeedPoint	(	const GlobalPoint &	dcaPoint,
		const TrackBaseRef	,
		const TrackBaseRef	,
		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 21 of file PFDisplacedVertexSeed.cc.

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

Referenced by reco::PFDisplacedVertexSeed::Compare::operator()().

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

  }

  addElement(r1); 
  addElement(r2);

}

Friends And Related Function Documentation

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

friend

Referenced by totalWeight().

Member Data Documentation

std::set< TrackBaseRef , Compare > reco::PFDisplacedVertexSeed::elements_

private

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

Set of tracks refs associated to the seed

Definition at line 87 of file PFDisplacedVertexSeed.h.

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

GlobalPoint reco::PFDisplacedVertexSeed::seedPoint_

private

Seed point which indicated the approximative position of the vertex.

Definition at line 89 of file PFDisplacedVertexSeed.h.

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

float reco::PFDisplacedVertexSeed::totalWeight_

private

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

Definition at line 92 of file PFDisplacedVertexSeed.h.

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