#include <Cluster1DCleaner.h>

Public Member Functions
	Cluster1DCleaner (const float zoffset, bool useErr)

std::vector< Cluster1D< T > >	clusters (const std::vector< Cluster1D< T > > &)

std::vector< Cluster1D< T > >	discardedCluster1Ds () const

Private Member Functions
float	average (const std::vector< Cluster1D< T > > &)

void	cleanCluster1Ds (const std::vector< Cluster1D< T > > &)

Private Attributes
std::vector< Cluster1D< T > >	theCleanedCluster1Ds

std::vector< Cluster1D< T > >	theDiscardedCluster1Ds

bool	theUseError

float	theZOffSet

Detailed Description

template<class T>
class Cluster1DCleaner< T >

Definition at line 13 of file Cluster1DCleaner.h.

Constructor & Destructor Documentation

template<class T>

Cluster1DCleaner< T >::Cluster1DCleaner	(	const float	zoffset,
		bool	useErr
	)

inline

Definition at line 15 of file Cluster1DCleaner.h.

References Cluster1DCleaner< T >::clusters(), Cluster1DCleaner< T >::theCleanedCluster1Ds, and Cluster1DCleaner< T >::theDiscardedCluster1Ds.

                                                      : theZOffSet(zoffset), theUseError(useErr) {
     theCleanedCluster1Ds.clear();
     theDiscardedCluster1Ds.clear();
   }

Member Function Documentation

template<class T >

float Cluster1DCleaner< T >::average ( const std::vector< Cluster1D< T > > & clust )

private

Definition at line 70 of file Cluster1DCleaner.h.

References runTheMatrix::err, and mathSSE::sqrt().

Referenced by average.Average::__str__(), Cluster1DCleaner< T >::cleanCluster1Ds(), and Cluster1DCleaner< T >::discardedCluster1Ds().

                                                                         {
   float ave = clust.front().position().value();
   float err = clust.front().position().error();
   for (typename std::vector<Cluster1D<T> >::const_iterator ic = (clust.begin()) + 1; ic != clust.end(); ic++) {
     float oldave = ave;
     float olderr = err;
     ave = (oldave / olderr / olderr + ic->position().value() / ic->position().error() / ic->position().error()) /
           (1. / olderr / olderr + 1. / ic->position().error() / ic->position().error());
     err = sqrt(olderr * olderr + ic->position().error() * ic->position().error());
   }
   return ave;
 }

template<class T >

void Cluster1DCleaner< T >::cleanCluster1Ds ( const std::vector< Cluster1D< T > > & _clust )

private

Definition at line 49 of file Cluster1DCleaner.h.

References Cluster1DCleaner< T >::average(), listHistos::discr, Cluster1DCleaner< T >::theCleanedCluster1Ds, Cluster1DCleaner< T >::theDiscardedCluster1Ds, Cluster1DCleaner< T >::theUseError, and Cluster1DCleaner< T >::theZOffSet.

Referenced by Cluster1DCleaner< T >::clusters(), and Cluster1DCleaner< T >::discardedCluster1Ds().

                                                                                 {
   std::vector<Cluster1D<T> > clust = _clust;
   theCleanedCluster1Ds.clear();
   theDiscardedCluster1Ds.clear();
   if (clust.size() == 0)
     return;
   float oldPos = average(clust);
   for (typename std::vector<Cluster1D<T> >::const_iterator ic = clust.begin(); ic != clust.end(); ic++) {
     float discr = theUseError ? fabs(((*ic).position().value() - oldPos) / (*ic).position().error())
                               : fabs(((*ic).position().value() - oldPos));
     if (discr < theZOffSet) {
       theCleanedCluster1Ds.push_back(*ic);
     } else {
       theDiscardedCluster1Ds.push_back(*ic);
     }
   }
   return;
 }

template<class T >

std::vector< Cluster1D< T > > Cluster1DCleaner< T >::clusters ( const std::vector< Cluster1D< T > > & _clust )

Definition at line 42 of file Cluster1DCleaner.h.

References Cluster1DCleaner< T >::cleanCluster1Ds(), and Cluster1DCleaner< T >::theCleanedCluster1Ds.

Referenced by Cluster1DCleaner< T >::Cluster1DCleaner().

                                                                                              {
   std::vector<Cluster1D<T> > clust = _clust;
   cleanCluster1Ds(clust);
   return theCleanedCluster1Ds;
 }