reco::Cluster3D Class Reference

#include <Cluster3D.h>

Public Member Functions
	Cluster3D ()
	Cluster3D (unsigned statusBits, const PrincipalComponents &pcaResults, float totalCharge, const float *startPosition, const float *endPosition, int idx)
unsigned	getStatusBits () const
const PrincipalComponents &	getPcaResults () const
float	getTotalCharge () const
const float *	getStartPosition () const
const float *	getEndPosition () const
int	getClusterIdx () const
void	setStatusBit (unsigned bits) const
void	clearStatusBits (unsigned bits) const
Cluster3D	operator+ (Cluster3D)

Private Attributes
unsigned	m_statusBits
	Default constructor. More...
PrincipalComponents	m_pcaResults
	Output of the prinicipal componenets analysis. More...
float	m_totalCharge
	Total charge in the cluster. More...
float	m_startPosition [3]
	"start" position for cluster (world coordinates) More...
float	m_endPosition [3]
	"end" position for cluster More...
int	m_clusterIdx
	ID for this cluster. More...

Friends
std::ostream &	operator<< (std::ostream &o, const Cluster3D &c)
bool	operator< (const Cluster3D &a, const Cluster3D &b)

Detailed Description

Definition at line 259 of file Cluster3D.h.

Constructor & Destructor Documentation

reco::Cluster3D::Cluster3D

(

)

Definition at line 243 of file Cluster3D.cxx.

                     : m_statusBits(0),
                         m_pcaResults(PrincipalComponents()),
                         m_totalCharge(0.),
                         m_startPosition{0.,0.,0.},
                                 m_endPosition{0.,0.,0.},
                                 m_clusterIdx(0)
{}

reco::Cluster3D::Cluster3D	(	unsigned	statusBits,
		const PrincipalComponents &	pcaResults,
		float	totalCharge,
		const float *	startPosition,
		const float *	endPosition,
		int	idx
	)

Definition at line 251 of file Cluster3D.cxx.

                                                     :
        m_statusBits(statusBits),
        m_pcaResults(pcaResults),
        m_totalCharge(totalCharge),
        m_startPosition{startPosition[0],startPosition[1],startPosition[2]},
        m_endPosition{endPosition[0],endPosition[1],endPosition[2]},
        m_clusterIdx(idx)
       {}

Member Function Documentation

void reco::Cluster3D::clearStatusBits ( unsigned  bits ) const

inline

Definition at line 290 of file Cluster3D.h.

{m_statusBits &= ~bits;}

int reco::Cluster3D::getClusterIdx ( ) const

inline

Definition at line 287 of file Cluster3D.h.

{return m_clusterIdx;}

const float* reco::Cluster3D::getEndPosition ( ) const

inline

Definition at line 286 of file Cluster3D.h.

{return m_endPosition;}

const PrincipalComponents& reco::Cluster3D::getPcaResults ( ) const

inline

Definition at line 283 of file Cluster3D.h.

{return m_pcaResults;}

const float* reco::Cluster3D::getStartPosition ( ) const

inline

Definition at line 285 of file Cluster3D.h.

{return m_startPosition;}

unsigned reco::Cluster3D::getStatusBits ( ) const

inline

Definition at line 282 of file Cluster3D.h.

{return m_statusBits;}

float reco::Cluster3D::getTotalCharge ( ) const

inline

Definition at line 284 of file Cluster3D.h.

{return m_totalCharge;}

Cluster3D reco::Cluster3D::operator+

(

Cluster3D

a

)

Definition at line 268 of file Cluster3D.cxx.

{
/*
    // throw exception if the clusters are not from the same plane
    if( a.View() != this->View() )
      throw cet::exception("Cluster+operator") << "Attempting to sum clusters from "
                 << "different views is not allowed\n";

    // check the start and end positions - for now the
    // smallest wire number means start position, largest means end position
    std::vector<float> astart(a.StartPos());
    std::vector<float> aend  (a.EndPos()  );
    std::vector<float> start(StartPos());
    std::vector<float> end  (EndPos()  );
    std::vector<float> sigstart(SigmaStartPos());
    std::vector<float> sigend  (SigmaEndPos()  );

    if(astart[0] < fStartPos[0]){
      start = astart;
      sigstart = a.SigmaStartPos();
    }

    if(aend[0] > fEndPos[0]){
      end = aend;
      sigend = a.SigmaEndPos();
    }

    //take weighted mean in obtaining average slope and differential charge,
    //based on total charge each cluster
    float dtdw = ((this->Charge()*dTdW()) + (a.Charge()*a.dTdW()))/(this->Charge() + a.Charge());
    float dqdw = ((this->Charge()*dQdW()) + (a.Charge()*a.dQdW()))/(this->Charge() + a.Charge());

    //hits.sort();//sort the PtrVector to organize Hits of new Cluster
    float sigdtdw = TMath::Max(SigmadTdW(), a.SigmadTdW());
    float sigdqdw = TMath::Max(SigmadQdW(), a.SigmadQdW());

    Cluster sum(//hits,
    start[0], sigstart[0],
    start[1], sigstart[1],
    end[0],   sigend[0],
    end[1],   sigend[1],
    dtdw, sigdtdw,
    dqdw, sigdqdw,
    this->Charge() + a.Charge(),
    this->View(),
    ID());
*/
    //return sum;
    return a;
}

void reco::Cluster3D::setStatusBit ( unsigned  bits ) const

inline

Definition at line 289 of file Cluster3D.h.

{m_statusBits |= bits;}

Friends And Related Function Documentation

bool operator< ( const Cluster3D &  a, const Cluster3D &  b  )

friend

Definition at line 342 of file Cluster3D.cxx.

{
/*
    if(a.View() != b.View())
      return a.View() < b.View();
    if(a.ID() != b. ID())
      return a.ID() < b.ID();
    if(a.StartPos()[0] != b.StartPos()[0])
      return a.StartPos()[0] < b.StartPos()[0];
    if(a.EndPos()[0] != b.EndPos()[0])
      return a.EndPos()[0] < b.EndPos()[0];
*/
    if (a.getStartPosition()[2] < b.getStartPosition()[2]) return true;

    return false; //They are equal
}

std::ostream& operator<< ( std::ostream &  o, const Cluster3D &  c  )

friend

Definition at line 322 of file Cluster3D.cxx.

{
    o << std::setiosflags(std::ios::fixed) << std::setprecision(2);
    o << "Cluster ID "    << std::setw(5)  << std::right << c.getClusterIdx();
//      << " : View = "     << std::setw(3)  << std::right << c.View()
//      << " StartWire = "  << std::setw(7)  << std::right << c.StartPos()[0]
//      << " EndWire = "    << std::setw(7)  << std::right << c.EndPos()[0]
//      << " StartTime = "  << std::setw(9)  << std::right << c.StartPos()[1]
//      << " EndTime = "    << std::setw(9)  << std::right << c.EndPos()[1]
//      << " dTdW = "       << std::setw(9)  << std::right << c.dTdW()
//      << " dQdW = "       << std::setw(9)  << std::right << c.dQdW()
//      << " Charge = "     << std::setw(10) << std::right << c.Charge();

    return o;
  }

Member Data Documentation

int reco::Cluster3D::m_clusterIdx

private

ID for this cluster.

Definition at line 272 of file Cluster3D.h.

float reco::Cluster3D::m_endPosition[3]

private

"end" position for cluster

Definition at line 271 of file Cluster3D.h.

PrincipalComponents reco::Cluster3D::m_pcaResults

private

Output of the prinicipal componenets analysis.

Definition at line 268 of file Cluster3D.h.

float reco::Cluster3D::m_startPosition[3]

private

"start" position for cluster (world coordinates)

Definition at line 270 of file Cluster3D.h.

unsigned reco::Cluster3D::m_statusBits

mutableprivate

Default constructor.

Status bits for the cluster

Definition at line 267 of file Cluster3D.h.

float reco::Cluster3D::m_totalCharge

private

Total charge in the cluster.

Definition at line 269 of file Cluster3D.h.

---

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

• larreco/larreco/RecoAlg/Cluster3DAlgs/Cluster3D.h
• larreco/larreco/RecoAlg/Cluster3DAlgs/Cluster3D.cxx