gar::rec::alg::AlgCluster< T > Class Template Reference

#include <NNClusters.h>

Public Member Functions
	AlgCluster (float xendcap)
gar::rec::Cluster *	operator() (GenericCluster< T > *c)

Protected Member Functions
	AlgCluster ()

Protected Attributes
float	_xendcap

Detailed Description

template<class T>
class gar::rec::alg::AlgCluster< T >

Helper class that creates an gar::rec::Cluster from a generic cluster with hit types that have a Position() and a Energy() method.

Definition at line 333 of file NNClusters.h.

Constructor & Destructor Documentation

template<class T>

gar::rec::alg::AlgCluster< T >::AlgCluster ( float  xendcap )

inline

Definition at line 337 of file NNClusters.h.

: _xendcap( xendcap ) {}

template<class T>

gar::rec::alg::AlgCluster< T >::AlgCluster ( )

inlineprotected

Definition at line 440 of file NNClusters.h.

{} ;

Member Function Documentation

template<class T>

gar::rec::Cluster* gar::rec::alg::AlgCluster< T >::operator() ( GenericCluster< T > *  c )

inline

Definition at line 339 of file NNClusters.h.

                                                                {

          gar::rec::Cluster* clu = new gar::rec::Cluster();

          unsigned n = c->size() ;
          unsigned i = 0 ;
          float rmin = 9999.;
          float rmax = 0.;
          std::vector<float> tmin;
          std::vector<float> tmax;

          std::vector<float> a, t, x, y, z;
          a.resize(n);
          t.resize(n);
          x.resize(n);
          y.resize(n);
          z.resize(n);

          for( typename GenericCluster<T>::iterator hi = c->begin(); hi != c->end() ; hi++) {

            T* hit = (*hi)->first ;

            a[i] = hit->Energy() ;
            t[i] = hit->Time().first ;
            x[i] = hit->Position()[0] ;
            y[i] = hit->Position()[1] ;
            z[i] = hit->Position()[2] ;

            clu->addHit( hit , a[i] ) ;

            //get time of first layer and last layer to compute the difference
            float r = std::sqrt(y[i]*y[i] + z[i]*z[i]);
            if( std::fabs(x[i]) < _xendcap ){
                //case in barrel
                if(rmin >= r)
                {
                    rmin = r;
                    tmin.push_back(t[i]);
                }
                if(rmax <= r)
                {
                    rmax = r;
                    tmax.push_back(t[i]);
                }
            }
            else{
                //case endcap
                if( rmin >= std::fabs(x[i]) )
                {
                    rmin = std::fabs(x[i]);
                    tmin.push_back(t[i]);
                }
                if( rmax <= std::fabs(x[i]) )
                {
                    rmax = std::fabs(x[i]);
                    tmax.push_back(t[i]);
                }
            }

            ++i ;
          }

          //get minimum time of the first layer
          auto time_first = std::min_element( tmin.begin(), tmin.end() );
          assert(time_first != tmin.end());
          //get minimum time of the last layer
          auto time_last = std::min_element( tmax.begin(), tmax.end() );
          assert(time_last != tmax.end());

          util::ClusterShapes cs(n, a, t, x, y, z) ;

          clu->setEnergy( cs.getTotalAmplitude()  ) ;
          clu->setTime( cs.getAverageTime(), *time_first - *time_last ) ;
          clu->setPosition( cs.getCenterOfGravity() ) ;

          // direction of cluster's PCA
          float* d = cs.getEigenVecInertia() ;
          clu->setEigenVectors( d );

          //Main eigenvector
          CLHEP::Hep3Vector v( d[0], d[1], d[2] ) ;

          clu->setITheta( v.theta() )  ;
          clu->setIPhi( v.phi() ) ;

          float param[6] ;

          // param[0] = cs.getElipsoid_r1() ;
          param[0] = cs.getElipsoid_r_forw();
          param[1] = cs.getElipsoid_r_back();
          param[2] = cs.getElipsoid_r2() ;
          param[3] = cs.getElipsoid_r3() ;
          param[4] = cs.getElipsoid_vol() ;
          param[5] = cs.getWidth() ;

          clu->setShape( param ) ;

          return clu;
        }

Member Data Documentation

template<class T>

float gar::rec::alg::AlgCluster< T >::_xendcap

protected

Definition at line 440 of file NNClusters.h.

---

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

• garsoft/RecoAlg/NNClusters.h