---

SMS Class Reference

Class to compute the SMS location estimator The SMS estimator is the mean value of a set of observations with Small Median of Squared distances. More...

#include <RecoVertex/VertexTools/interface/SMS.h>

List of all members.

Public Types
enum	SMSType { None = 0, Interpolate = 1, Iterate = 2, Weighted = 4 }
Public Member Functions
GlobalPoint	location (const vector< pair< GlobalPoint, float > > &) const
GlobalPoint	location (const vector< GlobalPoint > &) const
	SMS (SMSType tp=(SMSType)(Interpolate|Iterate|Weighted), float q=0.5)
	Constructor.
Private Attributes
float	theRatio
SMSType	theType

---

Detailed Description

Class to compute the SMS location estimator The SMS estimator is the mean value of a set of observations with Small Median of Squared distances.

Definition at line 18 of file SMS.h.

---

Member Enumeration Documentation

enum SMS::SMSType

Enumerator:

None
Interpolate
Iterate
Weighted

Definition at line 21 of file SMS.h.

               { None        = 0,
                 Interpolate = 1,
                 Iterate     = 2,
                 Weighted    = 4 };

---

Constructor & Destructor Documentation

SMS::SMS	(	SMSType	tp = (SMSType) (Interpolate | Iterate | Weighted),
		float	q = 0.5
	)

Constructor.

Parameters:

	tp	What specific kind of SMS algorithm do you want?
	q	What fraction of data points are considered for the "next step"?

Definition at line 73 of file SMS.cc.

: theType(tp) , theRatio(q) {}

---

Member Function Documentation

GlobalPoint SMS::location

(

const vector< pair< GlobalPoint, float > > &

)

const

GlobalPoint SMS::location

(

const vector< GlobalPoint > &

data

)

const

Definition at line 76 of file SMS.cc.

References average(), GenMuonPlsPt100GeV_cfg::cout, funct::D, lat::endl(), i, int, Interpolate, Iterate, j, mergeAndRegister_online::pairs, python::multivaluedict::sort(), theRatio, theType, tmp, and Weighted.

Referenced by SmsModeFinder3d::operator()().

{
  if ( theType & Weighted )
  {
    cout << "[SMS] warning: Weighted SMS was asked for, but data are "
         << "weightless!" << endl;
  };
  int nobs=data.size();
  int nq=(int) ceil( theRatio*nobs);
  // cout << "nobs= " << nobs << "  nq= " << nq << endl;

  // Compute distances
  vector<MyPair> pairs;

  for ( vector< GlobalPoint >::const_iterator i=data.begin(); i!=data.end() ; ++i )
  {
    vector < float > D;
    // Compute squared distances to all points
    for ( vector< GlobalPoint >::const_iterator j=data.begin(); j!=data.end() ; ++j )
    { D.push_back ( (*j - *i).mag2() ); }
    // Find q-quantile in each row of the distance matrix
    sort( D.begin(), D.end() );
    MyPair tmp (  D[nq-1], &(*i) );
    pairs.push_back ( tmp );
  };

  // Sort pairs by first element
  sort( pairs.begin(), pairs.end(), Sorter() );
  if ( !(theType & SMS::Interpolate) &&
       !(theType & SMS::Iterate) )
  {
    // we dont interpolate, we dont iterate, so we can stop right here.
    // cout << "No interpolation, no iteration" << endl;
    return *(pairs.begin()->second);
  };

  // we dont iterate, or we dont have anything to iterate (anymore?)
  // so we stop here

  // cout << "nobs= " << nobs << "  nq= " << nq << endl;
  if (!(theType & SMS::Iterate) || nq<=2)
    return average ( pairs, nq );

  // we iterate (recursively)

  vector < GlobalPoint > data1;
  vector<MyPair>::iterator j;

  for ( j=pairs.begin(); j-pairs.begin()<nq; ++j)
     data1.push_back(*(j->second));

  return this->location( data1 );

}

---

Member Data Documentation

float SMS::theRatio [private]

Definition at line 38 of file SMS.h.

Referenced by location().

SMSType SMS::theType [private]

Definition at line 37 of file SMS.h.

Referenced by location().

---

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

• RecoVertex/VertexTools/interface/SMS.h
• RecoVertex/VertexTools/src/SMS.cc

---