Trk::Z0PVTrackCompatibilityEstimator Class Reference

Z0PVTrackCompatibilityEstimator.h. More...

#include <Z0PVTrackCompatibilityEstimator.h>

Inheritance diagram for Trk::Z0PVTrackCompatibilityEstimator:

Collaboration diagram for Trk::Z0PVTrackCompatibilityEstimator:

Public Member Functions
	Z0PVTrackCompatibilityEstimator (const std::string &type, const std::string &name, const IInterface *parent)
virtual	~Z0PVTrackCompatibilityEstimator ()
virtual StatusCode	initialize ()
virtual StatusCode	finalize ()
bool	isCompatible (const xAOD::TrackParticle *track, const xAOD::VertexContainer *pvContainer, unsigned int pvIndex=0) const
bool	isCompatible (const TrackParameters *track, const xAOD::VertexContainer *pvContainer, unsigned int pvIndex=0) const

Private Attributes
ToolHandle< Trk::ITrackToVertexIPEstimator >	m_ipEstimator
double	m_z0SignificanceMax

Detailed Description

Z0PVTrackCompatibilityEstimator.h.

Assigns tracks to the primary vertex based on their z0 significance. Tracks that are close to a pileup vertex are not assigned to the primary vertex.

ruwie.nosp@m.del@.nosp@m.physi.nosp@m.k.un.nosp@m.i-bon.nosp@m.n.de

---

Changes:

David Shope david.nosp@m..ric.nosp@m.hard..nosp@m.shop.nosp@m.e@cer.nosp@m.n.ch (2016-04-19)

EDM Migration to xAOD - Replace Trk::VxContainer with xAOD::VertexContainer

Definition at line 44 of file Z0PVTrackCompatibilityEstimator.h.

Constructor & Destructor Documentation

Z0PVTrackCompatibilityEstimator()

Trk::Z0PVTrackCompatibilityEstimator::Z0PVTrackCompatibilityEstimator	(	const std::string &	type,
		const std::string &	name,
		const IInterface *	parent
	)

Definition at line 20 of file Z0PVTrackComaptibilityEstimator.cxx.

                                                             :
    base_class(type, name, parent),
    m_ipEstimator("Trk::TrackToVertexIPEstimator"),
    m_z0SignificanceMax(30.)
  {
    declareProperty("IPEstimator", m_ipEstimator, "Impact parameter estimator tool");
    declareProperty("Z0SignificanceMax", m_z0SignificanceMax);
  }

~Z0PVTrackCompatibilityEstimator()

Trk::Z0PVTrackCompatibilityEstimator::~Z0PVTrackCompatibilityEstimator ( )

virtualdefault

Member Function Documentation

finalize()

StatusCode Trk::Z0PVTrackCompatibilityEstimator::finalize ( )

virtual

Definition at line 48 of file Z0PVTrackComaptibilityEstimator.cxx.

                                                      {
 
    ATH_MSG_INFO( "Finalizing Z0PVTrackCompatibilityEstimator" );
 
    return StatusCode::SUCCESS; 
  } 

initialize()

StatusCode Trk::Z0PVTrackCompatibilityEstimator::initialize ( )

virtual

Definition at line 34 of file Z0PVTrackComaptibilityEstimator.cxx.

                                                        {
    
    ATH_MSG_INFO( "Initializing Z0PVTrackCompatibilityEstimator" );
 
    StatusCode sc = m_ipEstimator.retrieve();
    if( !sc.isSuccess() ){
      msg(MSG::FATAL) << "Unable to retrieve impact parameter estimator tool" << endmsg;
      return sc;
    }
 
    return StatusCode::SUCCESS; 
  }

isCompatible() [1/2]

bool Trk::Z0PVTrackCompatibilityEstimator::isCompatible	(	const TrackParameters *	track,
		const xAOD::VertexContainer *	pvContainer,
		unsigned int	pvIndex = 0
	)		const

Definition at line 63 of file Z0PVTrackComaptibilityEstimator.cxx.

                                                                                                         {
    
    // TODO: do we want to doRemoval here? (I have it as false for now) 
    std::unique_ptr<ImpactParametersAndSigma> ip =  m_ipEstimator->estimate(track, track, (*primaryVertexContainer)[index], false);
    
    bool pileup = false;
  
    double sumptPV = 0;
    
    /*
    for(std::vector<Trk::VxTrackAtVertex*>::const_iterator track = 
      (*primaryVertexContainer)[index]->vxTrackAtVertex()->begin();
    track != (*primaryVertexContainer)[index]->vxTrackAtVertex()->end(); ++track){
      
      sumptPV += (*track)->initialPerigee()->momentum().perp();
    }
    */
 
    // TODO: perhaps calculate sumptPV using VxTracksAtVertex in xAOD::Vertex? -David S.
    for(unsigned int i = 0; i < (*primaryVertexContainer)[index]->nTrackParticles(); ++i){
 
      sumptPV += (*primaryVertexContainer)[index]->trackParticle(i)->perigeeParameters().momentum().perp();
    }
 
    for(unsigned int i = 0; i < primaryVertexContainer->size(); i++){
      
      double sumpt = 0.;
      
      if( i != index &&
      (*primaryVertexContainer)[i]->vertexType() != xAOD::VxType::NoVtx ){
        /*
    for(std::vector<Trk::VxTrackAtVertex*>::const_iterator track = 
          (*primaryVertexContainer)[i]->vxTrackAtVertex()->begin();
        track != (*primaryVertexContainer)[i]->vxTrackAtVertex()->end(); ++track){
    
      sumpt += (*track)->initialPerigee()->momentum().perp();
    }
    */
 
        // TODO: perhaps calculate sumptPV using VxTracksAtVertex in xAOD::Vertex? -David S.
        for(unsigned int itrk = 0; itrk < (*primaryVertexContainer)[i]->nTrackParticles(); ++itrk){
 
          sumpt += (*primaryVertexContainer)[i]->trackParticle(itrk)->perigeeParameters().momentum().perp();
        }
    double vxReach = 
      std::abs((*primaryVertexContainer)[i]->position().z()
           - (*primaryVertexContainer)[index]->position().z()) * sumpt / (sumpt + sumptPV);
    
    if(std::abs(ip->IPz0 + (*primaryVertexContainer)[index]->position().z() -
            (*primaryVertexContainer)[i]->position().z()) < vxReach){
      
      pileup = true;
    }
      }
    }
  
    return std::abs(ip->IPz0) / ip->sigmaz0 < m_z0SignificanceMax && !pileup;
  }

isCompatible() [2/2]

bool Trk::Z0PVTrackCompatibilityEstimator::isCompatible	(	const xAOD::TrackParticle *	track,
		const xAOD::VertexContainer *	pvContainer,
		unsigned int	pvIndex = 0
	)		const

Definition at line 57 of file Z0PVTrackComaptibilityEstimator.cxx.

                                                                                                                                 {
 
    return isCompatible( &(track->perigeeParameters()), primaryVertexContainer, index );
  }

Member Data Documentation

m_ipEstimator

ToolHandle< Trk::ITrackToVertexIPEstimator > Trk::Z0PVTrackCompatibilityEstimator::m_ipEstimator

private

Definition at line 60 of file Z0PVTrackCompatibilityEstimator.h.

m_z0SignificanceMax

double Trk::Z0PVTrackCompatibilityEstimator::m_z0SignificanceMax

private

Definition at line 62 of file Z0PVTrackCompatibilityEstimator.h.

---

The documentation for this class was generated from the following files:

• Z0PVTrackCompatibilityEstimator.h
• Z0PVTrackComaptibilityEstimator.cxx