Algorithm collection class computing cluster parameters. More...

#include <StandardClusterParamsAlg.h>

Inheritance diagram for cluster::StandardClusterParamsAlg:

Public Types
using	Measure_t = ClusterParamsAlgBase::Measure_t

Public Types inherited from cluster::ClusterParamsAlgBase
using	Measure_t = details::Measure_t< float >
	Type used to return values with errors. More...

Public Member Functions
	StandardClusterParamsAlg ()
	Constructor. More...

virtual	~StandardClusterParamsAlg ()=default
	Destructor. More...

void	Clear () override
	Restores the class to post-configuration, pre-initialization state. More...

void	SetHits (util::GeometryUtilities const &gser, std::vector< recob::Hit const * > const &hits) override
	Sets the list of input hits. More...

void	SetHits (util::GeometryUtilities const &gser, std::vector< recob::Hit > const &hits) override
	Sets the list of input hits. More...

void	SetVerbose (int level=1) override
	Set the verbosity level. More...

size_t	NHits () override
	Returns the number of hits in the cluster. More...

float	MultipleHitDensity () override
	Fraction of wires in the cluster with more than one hit. More...

float	Width (util::GeometryUtilities const &gser) override
	Computes the width of the cluster. More...

size_t	NInputHits () const
	Returns the number of input hits. More...


Measure_t	StartCharge (util::GeometryUtilities const &gser) override
	Computes the charge on the first and last wire of the track. More...

Measure_t	EndCharge (util::GeometryUtilities const &gser) override


Measure_t	StartAngle () override
	Computes the angle of the cluster. More...

Measure_t	EndAngle () override


Measure_t	StartOpeningAngle () override
	Computes the opening angle at the start or end of the cluster. More...

Measure_t	EndOpeningAngle () override

Cluster charge
Measure_t	Integral () override
	Computes the total charge of the cluster from Hit::Integral() More...

Measure_t	IntegralStdDev () override
	Computes the standard deviation on the charge of the cluster hits. More...

Measure_t	SummedADC () override
	Computes the total charge of the cluster from Hit::SummedADC() More...

Measure_t	SummedADCStdDev () override
	Computes the standard deviation on the charge of the cluster hits. More...

Public Member Functions inherited from cluster::ClusterParamsAlgBase
virtual	~ClusterParamsAlgBase ()=default

Protected Attributes
ClusterParamsAlg	algo
	the actual algorithm class More...

Protected Attributes inherited from cluster::ClusterParamsAlgBase
int	verbose = 0
	verbosity level: 0 is normal, negative is even quieter More...

Additional Inherited Members
Static Protected Member Functions inherited from cluster::ClusterParamsAlgBase
static std::logic_error	NotImplemented (std::string function_name)

Detailed Description

Algorithm collection class computing cluster parameters.

See Also: ClusterParamsAlg

This class wraps ClusterParamsAlg class, designed in the context of shower reconstruction, to expose a standard ClusterParamsBaseAlg interface.

Definition at line 34 of file StandardClusterParamsAlg.h.

Member Typedef Documentation

using cluster::StandardClusterParamsAlg::Measure_t = ClusterParamsAlgBase::Measure_t

Definition at line 36 of file StandardClusterParamsAlg.h.

Constructor & Destructor Documentation

cluster::StandardClusterParamsAlg::StandardClusterParamsAlg ( )

Constructor.

Definition at line 20 of file StandardClusterParamsAlg.cxx.

 {
   SetVerbose(0);
 }

virtual cluster::StandardClusterParamsAlg::~StandardClusterParamsAlg ( )

virtualdefault

Destructor.

Member Function Documentation

void cluster::StandardClusterParamsAlg::Clear ( )

overridevirtual

Restores the class to post-configuration, pre-initialization state.

Reimplemented from cluster::ClusterParamsAlgBase.

Definition at line 35 of file StandardClusterParamsAlg.cxx.

 {
   algo.Initialize();
 }

cluster::StandardClusterParamsAlg::Measure_t cluster::StandardClusterParamsAlg::EndAngle ( )

overridevirtual

Reimplemented from cluster::ClusterParamsAlgBase.

Definition at line 78 of file StandardClusterParamsAlg.cxx.

 {
   return StartAngle(); // Ummm...this doesn't look right. FIXME
 }

cluster::StandardClusterParamsAlg::Measure_t cluster::StandardClusterParamsAlg::EndCharge ( util::GeometryUtilities const & gser )

overridevirtual

Reimplemented from cluster::ClusterParamsAlgBase.

Definition at line 60 of file StandardClusterParamsAlg.cxx.

 {
   if (NInputHits() == 0) return {0.F};
   return {(float)algo.EndCharge(gser)};
 }

cluster::StandardClusterParamsAlg::Measure_t cluster::StandardClusterParamsAlg::EndOpeningAngle ( )

overridevirtual

Reimplemented from cluster::ClusterParamsAlgBase.

Definition at line 95 of file StandardClusterParamsAlg.cxx.

 {
   if (NInputHits() < 3) return {0.F};
 
   algo.RefineDirection();
   return {(float)algo.GetParams().closing_angle_charge_wgt};
 }

cluster::StandardClusterParamsAlg::Measure_t cluster::StandardClusterParamsAlg::Integral ( )

overridevirtual

Computes the total charge of the cluster from Hit::Integral()

Returns: total charge of the cluster, in ADC count units

See Also: IntegralStdDev(), SummedADC()

ClusterParamsAlg computes the sum from all hits.

Reimplemented from cluster::ClusterParamsAlgBase.

Definition at line 105 of file StandardClusterParamsAlg.cxx.

 {
   if (NInputHits() == 0) return {0.F};
 
   algo.GetAverages();
   return {(float)algo.GetParams().sum_charge};
 }

cluster::StandardClusterParamsAlg::Measure_t cluster::StandardClusterParamsAlg::IntegralStdDev ( )

overridevirtual

Computes the standard deviation on the charge of the cluster hits.

Returns: the standard deviation of charge of hits, in ADC count units

See Also: Integral()

ClusterParamsAlg computes the standard deviation of the sample of charges from all hits. Hit charge is obtained by recob::Hit::Integral().

Reimplemented from cluster::ClusterParamsAlgBase.

Definition at line 115 of file StandardClusterParamsAlg.cxx.

 {
   if (NInputHits() < 2) return {0.F};
 
   algo.GetAverages();
   return {(float)algo.GetParams().rms_charge};
 }

float cluster::StandardClusterParamsAlg::MultipleHitDensity ( )

overridevirtual

Fraction of wires in the cluster with more than one hit.

Returns: fraction of wires with more than one hit, or 0 if no wires

Returns a quantity defined as NMultiHitWires / NWires, where NWires is the number of wires hosting at least one hit of this cluster, and NMultiHitWires is the number of wires which have more than just one hit.

Reimplemented from cluster::ClusterParamsAlgBase.

Definition at line 157 of file StandardClusterParamsAlg.cxx.

 {
   if (NInputHits() < 2) return 0.0F;
 
   algo.GetAverages();
   return algo.GetParams().N_Wires ? algo.GetParams().multi_hit_wires / algo.GetParams().N_Wires :
                                     0.;
 }

size_t cluster::StandardClusterParamsAlg::NHits ( )

overridevirtual

Returns the number of hits in the cluster.

Reimplemented from cluster::ClusterParamsAlgBase.

Definition at line 147 of file StandardClusterParamsAlg.cxx.

 {
   if (NInputHits() < 2) return NInputHits();
 
   algo.GetAverages();
   return (size_t)algo.GetParams().N_Hits;
 }

size_t cluster::StandardClusterParamsAlg::NInputHits ( ) const

Returns the number of input hits.

Definition at line 178 of file StandardClusterParamsAlg.cxx.

 {
   return algo.GetNHits();
 }

void cluster::StandardClusterParamsAlg::SetHits	(	util::GeometryUtilities const &	gser,
		std::vector< recob::Hit const * > const &	hits
	)

overridevirtual

Sets the list of input hits.

Parameters

hits	list of hits

Exceptions

undefined in case of error, this method can throw (anything)

Hits are translated into our own internal format. The original hits are not used afterward, and their distruction will not affect this object. This method calls Clear() at the beginning (although the protocol does not requires it).

Implements cluster::ClusterParamsAlgBase.

Definition at line 42 of file StandardClusterParamsAlg.cxx.

 {
   Clear();
   util::PxHitConverter pxhitconverter{gser};
   algo.SetHits(pxhitconverter.ToPxHitVector(hits));
 }

void cluster::StandardClusterParamsAlg::SetHits	(	util::GeometryUtilities const &	gser,
		std::vector< recob::Hit > const &	hits
	)

inlineoverridevirtual

Sets the list of input hits.

Parameters

hits	list of hits (hits will not be modified)

Exceptions

undefined in case of error, this method can throw (anything)

See Also: ClusterParamsAlgBase::SetHits(std::vector<recob::Hit> const&)

Reimplemented from cluster::ClusterParamsAlgBase.

Definition at line 68 of file StandardClusterParamsAlg.h.

     {
       ClusterParamsAlgBase::SetHits(gser, hits);
     }

void cluster::StandardClusterParamsAlg::SetVerbose ( int level = 1 )

overridevirtual

Set the verbosity level.

Reimplemented from cluster::ClusterParamsAlgBase.

Definition at line 27 of file StandardClusterParamsAlg.cxx.

 {
   ClusterParamsAlgBase::SetVerbose(level);
   algo.SetVerbose(level > 0);
 }

cluster::StandardClusterParamsAlg::Measure_t cluster::StandardClusterParamsAlg::StartAngle ( )

overridevirtual

Computes the angle of the cluster.

Returns: angle of the cluster in the wire x time space, in radians

Uses the coordinates from the hits, weighted by charge (Hit::Integral()) to compute a slope in the homogenized wire x time space. The homogenized space has both wires and ticks converted into a distance (by using detector parameters: wire pitch and drift velocity).

The angle is in the $[ -\pi, \pi ]$ range, with 0 corresponding to a cluster parallel to the wire plane and $\pi$ to a cluster orthogonal to the wire plane, going farther from it.

Note: Both the methods return the same value.

Reimplemented from cluster::ClusterParamsAlgBase.

Definition at line 68 of file StandardClusterParamsAlg.cxx.

 {
   if (NInputHits() < 2) return {0.F};
 
   algo.GetRoughAxis();
   return {(float)util::DegreesToRadians(algo.GetParams().cluster_angle_2d)};
 }

cluster::StandardClusterParamsAlg::Measure_t cluster::StandardClusterParamsAlg::StartCharge ( util::GeometryUtilities const & gser )

overridevirtual

Computes the charge on the first and last wire of the track.

Returns: the charge in ADC counts, with uncertainty

The implementation in ClusterParamsAlg provides an estimation of the charge collected in the first or last 1 cm of the cluster, using a linear fit on the deposited charge to reduce fluctuations.

Reimplemented from cluster::ClusterParamsAlgBase.

Definition at line 52 of file StandardClusterParamsAlg.cxx.

 {
   if (NInputHits() == 0) return {0.F};
   return {(float)algo.StartCharge(gser)};
 }

cluster::StandardClusterParamsAlg::Measure_t cluster::StandardClusterParamsAlg::StartOpeningAngle ( )

overridevirtual

Computes the opening angle at the start or end of the cluster.

Returns: angle at the start of the cluster, in radians

This algorithm returns an opening angle after weighting the hits by their charge (as defined bu Hit::Integral());

Reimplemented from cluster::ClusterParamsAlgBase.

Definition at line 85 of file StandardClusterParamsAlg.cxx.

 {
   if (NInputHits() < 3) return {0.F};
 
   algo.RefineDirection();
   return {(float)algo.GetParams().opening_angle_charge_wgt};
 }

cluster::StandardClusterParamsAlg::Measure_t cluster::StandardClusterParamsAlg::SummedADC ( )

overridevirtual

Computes the total charge of the cluster from Hit::SummedADC()

Returns: total charge of the cluster, in ADC count units

See Also: SummedADCStdDev(), Integral()

ClusterParamsAlg computes the sum from all hits.

Reimplemented from cluster::ClusterParamsAlgBase.

Definition at line 125 of file StandardClusterParamsAlg.cxx.

 {
   if (NInputHits() == 0) return {0.F};
 
   // compute all the averages
   algo.GetAverages();
   double sumADC = algo.GetParams().sum_ADC;
   return {(float)sumADC, (float)std::sqrt(sumADC)};
 }

cluster::StandardClusterParamsAlg::Measure_t cluster::StandardClusterParamsAlg::SummedADCStdDev ( )

overridevirtual

Computes the standard deviation on the charge of the cluster hits.

Returns: the standard deviation of charge of hits, in ADC count units

See Also: SummedADC()

ClusterParamsAlg computes the standard deviation of the sample of charges from all hits. Hit charge is obtained by recob::Hit::SummedADC().

Reimplemented from cluster::ClusterParamsAlgBase.

Definition at line 137 of file StandardClusterParamsAlg.cxx.

 {
   if (NInputHits() < 2) return {0.F};
 
   algo.GetAverages();
   return {(float)algo.GetParams().rms_ADC};
 }

float cluster::StandardClusterParamsAlg::Width ( util::GeometryUtilities const & gser )

overridevirtual

Computes the width of the cluster.

Returns: width of the cluster

Todo:: provide a description of the algorithm by words

Reimplemented from cluster::ClusterParamsAlgBase.

Definition at line 168 of file StandardClusterParamsAlg.cxx.

 {
   if (NInputHits() < 3) return 0.0F;
 
   algo.GetProfileInfo(gser);
   return algo.GetParams().width;
 }

Member Data Documentation

ClusterParamsAlg cluster::StandardClusterParamsAlg::algo

protected

the actual algorithm class

Definition at line 191 of file StandardClusterParamsAlg.h.

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

Public Types

Public Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation