gar::rec::Hit Class Reference

#include <Hit.h>

Public Member Functions
	Hit ()
	Hit (unsigned int chan, float sig, float *pos, float startT, float endT, float Time, float RMS)
const float *	Position () const
float const &	Signal () const
unsigned int	Channel () const
float	StartTime () const
float	EndTime () const
float	Time () const
float	RMS () const
void	operator+= (gar::rec::Hit const &h)
bool	operator< (gar::rec::Hit const &h) const
bool	operator== (gar::rec::Hit const &h) const

Private Attributes
unsigned int	fChannel
	channel recording this hit More...
float	fSignal
	size of the signal for this hit (integral of ADC values) More...
float	fPosition [3]
	position of the hit More...
float	fTime
	time of hit charge arrival at the readout plane (ticks) More...
float	fStartTime
	start time of the hit (ticks) More...
float	fEndTime
	end time of the hit (ticks) More...
float	fRMS
	Hit width calculated with RMS (in ticks) More...

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

Detailed Description

Definition at line 15 of file Hit.h.

Constructor & Destructor Documentation

gar::rec::Hit::Hit

(

)

Definition at line 15 of file Hit.cxx.

    {
      return;
    }

gar::rec::Hit::Hit	(	unsigned int	chan,
		float	sig,
		float *	pos,
		float	startT,
		float	endT,
		float	Time,
		float	RMS
	)

Definition at line 21 of file Hit.cxx.

    : fChannel  (chan  )
    , fSignal   (sig   )
    , fTime     (Time  )
    , fStartTime(startT)
    , fEndTime  (endT  )
    , fRMS      (RMS   )
    {
      
      fPosition[0] = pos[0];
      fPosition[1] = pos[1];
      fPosition[2] = pos[2];
      
      return;
    }

Member Function Documentation

unsigned int gar::rec::Hit::Channel ( ) const

inline

Definition at line 64 of file Hit.h.

{ return fChannel;      }

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

inline

Definition at line 68 of file Hit.h.

{ return fEndTime;      }

void gar::rec::Hit::operator+=

(

gar::rec::Hit const &

h

)

Definition at line 47 of file Hit.cxx.

    {
      // only add hits on the same channel
      if(fChannel != h.Channel()){
        MF_LOG_WARNING("Hit")
        << "attempting to add hits from different channels, that will not work";
        return;
      }
      
      // add h to this hit.  Set the new position, etc to be the weighted
      // average of the individual hits
      float totSig = fSignal + h.Signal();
      
      if(totSig == 0.){
        MF_LOG_WARNING("Hit")
        << "attempting to add two hits 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 hit
      
      //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::Hit::operator<

(

gar::rec::Hit const &

h

)

const

Definition at line 91 of file Hit.cxx.

    {
      if (Channel()   < h.Channel()  ) return true;
      if (Channel()   > h.Channel()  ) return false;
      if (StartTime() < h.StartTime()) return true;
      if (StartTime() > h.StartTime()) return false;
      if (EndTime()   < h.EndTime()  ) return true;
      if (EndTime()   > h.EndTime()  ) return false;
      // I dunno wot den
      //return true; //original
      return false;
      
    }

bool gar::rec::Hit::operator==

(

gar::rec::Hit const &

h

)

const

Definition at line 106 of file Hit.cxx.

    {
        //if(Position()!=h.Position())
        //    return false;
        if(Channel()!=h.Channel())
            return false;
        //if(Signal()!=h.Signal())
        //    return false;
        if(StartTime()!=h.StartTime())
            return false;
        if(EndTime()!=h.EndTime())
            return false;
        //if(RMS()!=h.RMS())
        //    return false;
        //if(Time()!=h.Time())
        //    return false;

        return true;

    }

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

inline

Definition at line 65 of file Hit.h.

{ return &fPosition[0]; }

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

inline

Definition at line 69 of file Hit.h.

{ return fRMS;          }

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

inline

Definition at line 66 of file Hit.h.

{ return fSignal;       }

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

inline

Definition at line 67 of file Hit.h.

{ return fStartTime;    }

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

inline

Definition at line 70 of file Hit.h.

{ return fTime;         }    

Friends And Related Function Documentation

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

friend

Definition at line 128 of file Hit.cxx.

    {
      
      o << "Hit on channel "
      << h.Channel()
      << "\n\tposition = ("
      << h.Position()[0]
      << ", "
      << h.Position()[1]
      << ", "
      << h.Position()[2]
      << ")"
      << "\n\tsignal = "
      << h.Signal()
      << "\n\tstart time: "
      << h.StartTime()
      << " end time: "
      << h.EndTime();
      
      return o;
    }

Member Data Documentation

unsigned int gar::rec::Hit::fChannel

private

channel recording this hit

Definition at line 24 of file Hit.h.

float gar::rec::Hit::fEndTime

private

end time of the hit (ticks)

Definition at line 29 of file Hit.h.

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

private

position of the hit

Definition at line 26 of file Hit.h.

float gar::rec::Hit::fRMS

private

Hit width calculated with RMS (in ticks)

Definition at line 30 of file Hit.h.

float gar::rec::Hit::fSignal

private

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

Definition at line 25 of file Hit.h.

float gar::rec::Hit::fStartTime

private

start time of the hit (ticks)

Definition at line 28 of file Hit.h.

float gar::rec::Hit::fTime

private

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

Definition at line 27 of file Hit.h.

---

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

• garsoft/ReconstructionDataProducts/Hit.h
• garsoft/ReconstructionDataProducts/Hit.cxx