gar::rec::TPCCluster Class Reference

#include <TPCCluster.h>

Public Member Functions
	TPCCluster ()
	TPCCluster (const float sig, const float *pos, const float startT, const float endT, const float Time, const float RMS, const float *cov)
bool	operator== (const TPCCluster &rhs) const
bool	operator!= (const TPCCluster &rhs) const
gar::rec::IDNumber	getIDNumber () const
const float *	Position () const
float const &	Signal () const
float	StartTime () const
float	EndTime () const
float	Time () const
float	RMS () const
const float *	CovMatPacked () const
void	operator+= (gar::rec::TPCCluster const &h)

Private Attributes
gar::rec::IDNumber	fIDnumero
float	fSignal
	size of the signal for this TPCCluster (integral of ADC values) More...
float	fPosition [3]
	position of the TPCCluster More...
float	fTime
	time of TPCCluster charge arrival at the readout plane (ticks) More...
float	fStartTime
	start time of the TPCCluster (ticks) More...
float	fEndTime
	end time of the TPCCluster (ticks) More...
float	fRMS
	TPCCluster width calculated with RMS (in ticks) More...
float	fCovMat [6]
	packed covariance matrix, assuming symmetry. xx, xy, xz, yy, yz, zz More...

Static Private Attributes
static gar::rec::IDNumber const	FirstNumber = 100100000

Friends
std::ostream &	operator<< (std::ostream &o, gar::rec::TPCCluster const &h)

Detailed Description

Definition at line 16 of file TPCCluster.h.

Constructor & Destructor Documentation

gar::rec::TPCCluster::TPCCluster

(

)

Definition at line 15 of file TPCCluster.cxx.

    {
      IDNumberGen::create(FirstNumber);
      fIDnumero = IDNumberGen::create()->getNewOne();
      return;
    }

gar::rec::TPCCluster::TPCCluster	(	const float	sig,
		const float *	pos,
		const float	startT,
		const float	endT,
		const float	Time,
		const float	RMS,
		const float *	cov
	)

Definition at line 37 of file TPCCluster.cxx.

      : fSignal   (sig   )
      , fTime     (Time  )
      , fStartTime(startT)
      , fEndTime  (endT  )
      , fRMS      (RMS   )
    {
      IDNumberGen::create(FirstNumber);
      fIDnumero = IDNumberGen::create()->getNewOne();
      
      fPosition[0] = pos[0];
      fPosition[1] = pos[1];
      fPosition[2] = pos[2];
      for (size_t i=0; i<6; ++i)
        {
          fCovMat[i] = cov[i];
        }
      
      return;
    }

Member Function Documentation

const float * gar::rec::TPCCluster::CovMatPacked ( ) const

inline

Definition at line 74 of file TPCCluster.h.

{ return &fCovMat[0]; }

float gar::rec::TPCCluster::EndTime ( ) const

inline

Definition at line 71 of file TPCCluster.h.

{ return fEndTime;      }

gar::rec::IDNumber gar::rec::TPCCluster::getIDNumber

(

)

const

Definition at line 32 of file TPCCluster.cxx.

{return fIDnumero;}

bool gar::rec::TPCCluster::operator!=

(

const TPCCluster &

rhs

)

const

Definition at line 28 of file TPCCluster.cxx.

                                                           {
      return (this->fIDnumero != rhs.fIDnumero);
    }

void gar::rec::TPCCluster::operator+=

(

gar::rec::TPCCluster const &

h

)

Definition at line 68 of file TPCCluster.cxx.

    {
      
      // add h to this TPCCluster.  Set the new position, etc to be the weighted
      // average of the individual TPCClusters
      float totSig = fSignal + h.Signal();
      
      if(totSig == 0.){
        MF_LOG_WARNING("TPCCluster")
          << "attempting to add two TPCClusters and neithr has any signal, bail.";
        return;
      }
 
      for(size_t i = 0; i < 3; ++i)
        fPosition[i] = (fPosition[i] * fSignal + h.Position()[i] * h.Signal()) / totSig;

      // don't do a weighted average but just make a wider TPCCluster
      
      //fStartTime = (fStartTime * fSignal + h.StartTime() * h.Signal()) / totSig;
      //fEndTime   = (fEndTime   * fSignal + h.EndTime()   * h.Signal()) / totSig;
      
      fStartTime = TMath::Min(fStartTime,h.StartTime());
      fEndTime   = TMath::Max(fEndTime,h.EndTime());

      float avgtime = (fTime * fSignal + h.Time() * h.Signal()) / totSig;

      fRMS = TMath::Sqrt(  (fSignal*(TMath::Sq(fTime-avgtime)+TMath::Sq(fRMS))
                            + h.Signal()*(TMath::Sq(h.Time()-avgtime)+TMath::Sq(h.RMS())))/totSig );
                          
      fTime      = avgtime;

      fSignal += h.Signal();
      
      return;
    }

bool gar::rec::TPCCluster::operator==

(

const TPCCluster &

rhs

)

const

Definition at line 24 of file TPCCluster.cxx.

                                                           {
      return (this->fIDnumero == rhs.fIDnumero);
    }

const float * gar::rec::TPCCluster::Position ( ) const

inline

Definition at line 68 of file TPCCluster.h.

{ return &fPosition[0]; }

float gar::rec::TPCCluster::RMS ( ) const

inline

Definition at line 72 of file TPCCluster.h.

{ return fRMS;          }

float const & gar::rec::TPCCluster::Signal ( ) const

inline

Definition at line 69 of file TPCCluster.h.

{ return fSignal;       }

float gar::rec::TPCCluster::StartTime ( ) const

inline

Definition at line 70 of file TPCCluster.h.

{ return fStartTime;    }

float gar::rec::TPCCluster::Time ( ) const

inline

Definition at line 73 of file TPCCluster.h.

{ return fTime;         }    

Friends And Related Function Documentation

std::ostream& operator<< ( std::ostream &  o, gar::rec::TPCCluster const &  h  )

friend

Definition at line 106 of file TPCCluster.cxx.

    {
      const float* cov;
      cov = h.CovMatPacked();
      o << "TPCCluster "
        << "\n\tID number = "
        << h.getIDNumber()
        << "\n\tposition = ("
        << h.Position()[0]
        << ", "
        << h.Position()[1]
        << ", "
        << h.Position()[2]
        << ")"
        << "\n\tcovariance matrix = \n\t"
        << cov[0] << "\t" << cov[1] << "\t" << cov[2] << "\n"
        <<         "\t\t" << cov[3] << "\t" << cov[4] << "\n"
        <<                         "\t\t\t" << cov[5]
        << "\n\tsignal = "
        << h.Signal()
        << "\n\tdrift direction rms = "
        << h.RMS()
        << "\n\tstart time: "
        << h.StartTime()
        << "\n\tend time: "
        << h.EndTime();
      
      return o;
    }

Member Data Documentation

float gar::rec::TPCCluster::fCovMat[6]

private

packed covariance matrix, assuming symmetry. xx, xy, xz, yy, yz, zz

Definition at line 33 of file TPCCluster.h.

float gar::rec::TPCCluster::fEndTime

private

end time of the TPCCluster (ticks)

Definition at line 31 of file TPCCluster.h.

gar::rec::IDNumber gar::rec::TPCCluster::fIDnumero

private

Definition at line 25 of file TPCCluster.h.

gar::rec::IDNumber const gar::rec::TPCCluster::FirstNumber = 100100000

staticprivate

Definition at line 24 of file TPCCluster.h.

float gar::rec::TPCCluster::fPosition[3]

private

position of the TPCCluster

Definition at line 28 of file TPCCluster.h.

float gar::rec::TPCCluster::fRMS

private

TPCCluster width calculated with RMS (in ticks)

Definition at line 32 of file TPCCluster.h.

float gar::rec::TPCCluster::fSignal

private

size of the signal for this TPCCluster (integral of ADC values)

Definition at line 27 of file TPCCluster.h.

float gar::rec::TPCCluster::fStartTime

private

start time of the TPCCluster (ticks)

Definition at line 30 of file TPCCluster.h.

float gar::rec::TPCCluster::fTime

private

time of TPCCluster charge arrival at the readout plane (ticks)

Definition at line 29 of file TPCCluster.h.

---

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

• garsoft/ReconstructionDataProducts/TPCCluster.h
• garsoft/ReconstructionDataProducts/TPCCluster.cxx