#include <GapClusterizerInZ.h>

Inheritance diagram for GapClusterizerInZ:

Public Member Functions
std::vector< std::vector < reco::TransientTrack > >	clusterize (const std::vector< reco::TransientTrack > &tracks) const

	GapClusterizerInZ (const edm::ParameterSet &conf)

std::vector< TransientVertex >	vertices (const std::vector< reco::TransientTrack > &tracks) const

float	zSeparation () const

	~GapClusterizerInZ ()

Public Member Functions inherited from TrackClusterizerInZ
	TrackClusterizerInZ ()

	TrackClusterizerInZ (const edm::ParameterSet &conf)

virtual	~TrackClusterizerInZ ()

Private Attributes
bool	verbose_

float	zSep

Detailed Description

Description: separates event tracks into clusters along the beam line

Definition at line 17 of file GapClusterizerInZ.h.

Constructor & Destructor Documentation

GapClusterizerInZ::GapClusterizerInZ ( const edm::ParameterSet & conf )

Definition at line 21 of file GapClusterizerInZ.cc.

References gather_cfg::cout, edm::ParameterSet::getParameter(), and edm::ParameterSet::getUntrackedParameter().

 {
   // some defaults to avoid uninitialized variables
   verbose_= conf.getUntrackedParameter<bool>("verbose", false);
   zSep = conf.getParameter<double>("zSeparation");
   if(verbose_) {std::cout << "TrackClusterizerInZ:  algorithm=gap, zSeparation="<< zSep << std::endl;}
 }

GapClusterizerInZ::~GapClusterizerInZ ( )

inline

Definition at line 31 of file GapClusterizerInZ.h.

31 {};

Member Function Documentation

vector< vector< reco::TransientTrack > > GapClusterizerInZ::clusterize ( const std::vector< reco::TransientTrack > & tracks ) const

virtual

Implements TrackClusterizerInZ.

Definition at line 40 of file GapClusterizerInZ.cc.

References abs, and testEve_cfg::tracks.

 {
   
   vector<reco::TransientTrack> tks = tracks; // copy to be sorted
 
   vector< vector<reco::TransientTrack> > clusters;
   if (tks.empty()) return clusters;
 
   // sort in increasing order of z
   stable_sort(tks.begin(), tks.end(), recTrackLessZ);
 
   // init first cluster
   vector<reco::TransientTrack>::const_iterator it = tks.begin();
   vector <reco::TransientTrack> currentCluster; currentCluster.push_back(*it);
 
   it++;
   for ( ; it != tks.end(); it++) {
 
     double zPrev = currentCluster.back().stateAtBeamLine().trackStateAtPCA().position().z();
     double zCurr = (*it).stateAtBeamLine().trackStateAtPCA().position().z();
 
     if ( abs(zCurr - zPrev) < zSeparation() ) {
       // close enough ? cluster together
       currentCluster.push_back(*it);
     }
     else {
       // store current cluster, start new one
       clusters.push_back(currentCluster);
       currentCluster.clear();
       currentCluster.push_back(*it);
     }
   }
 
   // store last cluster
   clusters.push_back(currentCluster);
 
   return clusters;
 
 }

std::vector< TransientVertex > GapClusterizerInZ::vertices ( const std::vector< reco::TransientTrack > & tracks ) const

float GapClusterizerInZ::zSeparation ( ) const

Definition at line 31 of file GapClusterizerInZ.cc.

 {
   return zSep;
 }

Member Data Documentation

bool GapClusterizerInZ::verbose_

private

Definition at line 35 of file GapClusterizerInZ.h.

float GapClusterizerInZ::zSep