pma::Node3D Class Reference

#include "PmaNode3D.h"

Inheritance diagram for pma::Node3D:

Public Member Functions
	Node3D ()
	Node3D (detinfo::DetectorPropertiesData const &detProp, const TVector3 &p3d, unsigned int tpc, unsigned int cryo, bool vtx=false, double xshift=0)
TVector3 const &	Point3D () const
bool	SetPoint3D (const TVector3 &p3d)
TVector2 const &	Projection2D (unsigned int view) const
double	GetDistToWall () const
bool	SameTPC (const TVector3 &p3d, float margin=0.0F) const
	Check if p3d is in the same TPC as the node. More...
bool	SameTPC (const pma::Vector3D &p3d, float margin=0.0F) const
bool	IsBranching () const
	Belongs to more than one track? More...
bool	IsTPCEdge () const
	Is the first/last in this TPC? More...
bool	IsVertex () const
	Check fIsVertex flag. More...
void	SetVertex (bool state)
void	SetVertexToBranching (bool setAllNodes)
std::vector< pma::Track3D * >	GetBranches () const
double	GetDistance2To (const TVector3 &p3d) const override
	Distance [cm] from the 3D point to the point 3D. More...
double	GetDistance2To (const TVector2 &p2d, unsigned int view) const override
pma::Vector3D	GetDirection3D () const override
TVector3	GetUnconstrainedProj3D (const TVector2 &, unsigned int) const override
	In case of a node it is simply 3D position of the node. More...
void	SetProjection (pma::Hit3D &h) const override
	Set hit 3D position and its 2D projection to the vertex. More...
double	Length2 () const override
double	SegmentCos () const
	Cosine of 3D angle between connected segments. More...
double	SegmentCosWirePlane () const
double	SegmentCosTransverse () const
double	GetObjFunction (float penaltyValue, float endSegWeight) const
	Objective function minimized during oprimization. More...
void	Optimize (float penaltyValue, float endSegWeight)
void	ClearAssigned (pma::Track3D *trk=0) override
void	ApplyDriftShift (double dx)
double	GetDriftShift () const
int	TPC (void) const
	TPC index or -1 if out of any TPC. More...
int	Cryo (void) const
	Cryostat index or -1 if out of any cryostat. More...
double	Length (void) const
const std::vector< pma::Hit3D * > &	Hits (void) const
bool	HasHit (const pma::Hit3D *h) const
pma::Hit3D &	Hit (size_t index)
void	RemoveHitAt (size_t index)
void	AddHit (pma::Hit3D *h)
size_t	NHits (void) const
unsigned int	NHits (unsigned int view) const
size_t	NEnabledHits (unsigned int view=geo::kUnknown) const
size_t	NPrecalcEnabledHits (void) const
TVector3 const &	ReferencePoint (size_t index) const
size_t	NPoints (void) const
void	AddPoint (TVector3 *p)
void	UpdateHitParams (void)
void	UpdateProjection (void)
void	SortHits (void)
double	SumDist2 (void) const
double	SumDist2 (unsigned int view) const
double	SumHitsQ (unsigned int view) const
unsigned int	NThisHits (unsigned int view) const
double	HitsRadius3D (unsigned int view) const
bool	IsFrozen (void) const
	Check if the vertex 3D position is fixed. More...
void	SetFrozen (bool state)
	Fix / relese vertex 3D position. More...
bool	SelectRndHits (size_t nmax_per_view)
bool	SelectAllHits (void)
virtual void	Disconnect (void)
virtual bool	AddNext (pma::SortedObjectBase *nextElement)
virtual int	RemoveNext (pma::SortedObjectBase *nextElement)
virtual pma::SortedObjectBase *	Next (unsigned int index=0) const
virtual unsigned int	NextCount (void) const
virtual bool	IsLast (void) const
virtual bool	IsFirst (void) const
virtual pma::SortedObjectBase *	Prev (void) const

Static Public Member Functions
static void	SetMargin (double m)
	Set allowed node position margin around TPC. More...
static float	OptFactor (unsigned int view)
static void	SetOptFactor (unsigned int view, float value)

Protected Attributes
int	fTPC
int	fCryo
bool	fFrozen
std::vector< pma::Hit3D * >	fAssignedHits
std::vector< TVector3 * >	fAssignedPoints
size_t	fNThisHits [3]
size_t	fNThisHitsEnabledAll
size_t	fNHits [3]
double	fSumHitsQ [3]
double	fHitsRadius
std::vector< pma::SortedObjectBase * >	next_vector
pma::SortedObjectBase *	next
pma::SortedObjectBase *	prev

Static Protected Attributes
static float	fOptFactors [3] = {0.2F, 0.8F, 1.0F}

Private Member Functions
bool	LimitPoint3D ()
	Returns true if node position was trimmed to its TPC volume + fMargin. More...
void	UpdateProj2D ()
double	EndPtCos2Transverse () const
double	PiInWirePlane () const
double	PenaltyInWirePlane () const
double	Pi (float endSegWeight, bool doAsymm) const
double	Penalty (float endSegWeight) const
double	Mse () const
double	MakeGradient (float penaltyValue, float endSegWeight)
double	StepWithGradient (float alfa, float tol, float penalty, float weight)
double	SumDist2Hits () const override

Private Attributes
geo::TPCGeo const &	fTpcGeo
geo::WireReadoutGeom const &	fChannelMap
double	fMinX
double	fMaxX
double	fMinY
double	fMaxY
double	fMinZ
double	fMaxZ
TVector3	fPoint3D
TVector2	fProj2D [3]
double	fDriftOffset
TVector3	fGradient
bool	fIsVertex

Static Private Attributes
static bool	fGradFixed [3] = {false, false, false}
static double	fMargin = 3.0

Detailed Description

Definition at line 31 of file PmaNode3D.h.

Constructor & Destructor Documentation

pma::Node3D::Node3D

(

)

Definition at line 29 of file PmaNode3D.cxx.

References pma::Element3D::fCryo, fProj2D, and pma::Element3D::fTPC.

  : fTpcGeo(art::ServiceHandle<geo::Geometry const>()->TPC())
  , fChannelMap(art::ServiceHandle<geo::WireReadout>()->Get())
  , fMinX(0)
  , fMaxX(0)
  , fMinY(0)
  , fMaxY(0)
  , fMinZ(0)
  , fMaxZ(0)
  , fPoint3D(0, 0, 0)
  , fDriftOffset(0)
  , fIsVertex(false)
{
  fTPC = 0;
  fCryo = 0;

  fProj2D[0].Set(0);
  fProj2D[1].Set(0);
  fProj2D[2].Set(0);
}

pma::Node3D::Node3D	(	detinfo::DetectorPropertiesData const &	detProp,
		const TVector3 &	p3d,
		unsigned int	tpc,
		unsigned int	cryo,
		bool	vtx = false,
		double	xshift = 0
	)

Definition at line 50 of file PmaNode3D.cxx.

References detinfo::DetectorPropertiesData::ConvertTicksToX(), fChannelMap, pma::Element3D::fCryo, fMaxX, fMaxY, fMaxZ, fMinX, fMinY, fMinZ, pma::Element3D::fTPC, fTpcGeo, geo::WireReadoutGeom::HasPlane(), geo::kZ, geo::BoxBoundedGeo::MaxY(), geo::BoxBoundedGeo::MaxZ(), geo::BoxBoundedGeo::MinY(), geo::BoxBoundedGeo::MinZ(), detinfo::DetectorPropertiesData::NumberTimeSamples(), SetPoint3D(), and tmp.

  : fTpcGeo(art::ServiceHandle<geo::Geometry const>()->TPC(geo::TPCID(cryo, tpc)))
  , fChannelMap(art::ServiceHandle<geo::WireReadout>()->Get())
  , fDriftOffset(xshift)
  , fIsVertex(vtx)
{
  fTPC = tpc;
  fCryo = cryo;

  unsigned int lastPlane = geo::kZ;
  while ((lastPlane > 0) && !fChannelMap.HasPlane(geo::PlaneID{fTpcGeo.ID(), lastPlane}))
    lastPlane--;

  fMinX = detProp.ConvertTicksToX(0, lastPlane, tpc, cryo);
  fMaxX = detProp.ConvertTicksToX(detProp.NumberTimeSamples() - 1, lastPlane, tpc, cryo);
  if (fMaxX < fMinX) {
    double tmp = fMaxX;
    fMaxX = fMinX;
    fMinX = tmp;
  }

  fMinY = fTpcGeo.MinY();
  fMaxY = fTpcGeo.MaxY();
  fMinZ = fTpcGeo.MinZ();
  fMaxZ = fTpcGeo.MaxZ();

  SetPoint3D(p3d);
}

Member Function Documentation

void pma::Element3D::AddHit ( pma::Hit3D *  h )

inlineinherited

Definition at line 70 of file PmaElement3D.h.

References pma::Element3D::fAssignedHits, and pma::Element3D::SetProjection().

Referenced by pma::Track3D::MakeProjection().

  {
    fAssignedHits.push_back(h);
    SetProjection(*h);
  }

bool pma::SortedBranchBase::AddNext ( pma::SortedObjectBase *  nextElement )

virtualinherited

Reimplemented from pma::SortedObjectBase.

Definition at line 69 of file SortedObjects.cxx.

References pma::SortedObjectBase::next, art::productstatus::present(), pma::SortedObjectBase::prev, and pma::SortedObjectBase::RemoveNext().

Referenced by pma::Track3D::AttachBackToSameTPC(), and pma::Track3D::AttachToSameTPC().

{
  if (!nextElement) {
    mf::LogError("pma::SortedBranchBase") << "Next == 0.";
    return false;
  }

  if (nextElement == this) {
    mf::LogWarning("pma::SortedBranchBase") << "Next == This.";
    return false;
  }

  bool present = false;
  for (size_t i = 0; i < next_vector.size(); i++) {
    if (next_vector[i] == nextElement) {
      mf::LogWarning("pma::SortedBranchBase") << "Contained.";
      present = true;
      break;
    }
  }
  if (!present) {
    if (nextElement->prev) // && (nextElement->prev != this)
      nextElement->prev->RemoveNext(nextElement);

    next = nextElement;
    next->prev = this;
    next_vector.push_back(next);
    return true;
  }
  else
    return false;
}

void pma::Element3D::AddPoint ( TVector3 *  p )

inlineinherited

Definition at line 82 of file PmaElement3D.h.

References pma::Element3D::ClearAssigned(), pma::Element3D::fAssignedPoints, and pma::Element3D::UpdateHitParams().

Referenced by pma::Track3D::MakeProjection().

{ fAssignedPoints.push_back(p); }

void pma::Node3D::ApplyDriftShift ( double  dx )

inline

Definition at line 113 of file PmaNode3D.h.

Referenced by pma::Track3D::ApplyDriftShiftInTree().

  {
    fPoint3D[0] += dx;
    fDriftOffset += dx;
  }

void pma::Node3D::ClearAssigned ( pma::Track3D *  trk = 0 )

overridevirtual

Clear hits/points vectors of this element, optionally only those which are owned by given track.

Reimplemented from pma::Element3D.

Definition at line 820 of file PmaNode3D.cxx.

References pma::Element3D::fAssignedHits, pma::Element3D::fAssignedPoints, pma::Element3D::fHitsRadius, pma::SortedBranchBase::Next(), pma::SortedBranchBase::NextCount(), pma::Segment3D::Parent(), and pma::SortedObjectBase::Prev().

Referenced by pma::Track3D::Split().

{
  if (!trk) {
    // like in the base class:
    fAssignedPoints.clear();
    fAssignedHits.clear();
  }
  else {
    std::vector<pma::Track3D*> to_check;
    pma::Segment3D* seg;
    if (Prev()) {
      seg = static_cast<pma::Segment3D*>(Prev());
      if (seg->Parent() != trk) to_check.push_back(seg->Parent());
    }
    for (unsigned int i = 0; i < NextCount(); i++) {
      seg = static_cast<pma::Segment3D*>(Next(i));
      if (seg->Parent() != trk) to_check.push_back(seg->Parent());
    }

    unsigned int p = 0;
    while (p < fAssignedPoints.size()) {
      bool found = false;
      for (size_t t = 0; t < to_check.size(); t++)
        if (to_check[t]->HasRefPoint(fAssignedPoints[p])) {
          found = true;
          break;
        }

      if (!found)
        fAssignedPoints.erase(fAssignedPoints.begin() + p);
      else
        p++;
    }

    unsigned int h = 0;
    while (h < fAssignedHits.size()) {
      bool found = false;
      pma::Hit3D* hit = fAssignedHits[h];

      for (size_t t = 0; (t < to_check.size()) && !found; t++)
        for (size_t i = 0; i < to_check[t]->size(); i++) {
          pma::Hit3D* pmaHit = static_cast<pma::Hit3D*>((*(to_check[t]))[i]);
          if (hit == pmaHit) {
            found = true;
            break;
          }
        }

      if (!found)
        fAssignedHits.erase(fAssignedHits.begin() + h);
      else
        h++;
    }
  }

  fHitsRadius = 0.0F;
}

int pma::Element3D::Cryo ( void  ) const

inlineinherited

Cryostat index or -1 if out of any cryostat.

Definition at line 37 of file PmaElement3D.h.

References pma::Element3D::fCryo, pma::Element3D::GetDirection3D(), pma::Element3D::GetDistance2To(), pma::Element3D::GetUnconstrainedProj3D(), pma::Element3D::Length2(), and pma::Element3D::SetProjection().

Referenced by pma::ProjectionMatchingAlg::mergeTracks(), pma::Segment3D::Segment3D(), and pma::Track3D::Split().

{ return fCryo; }

void pma::SortedBranchBase::Disconnect ( void  )

virtualinherited

Reimplemented from pma::SortedObjectBase.

Definition at line 62 of file SortedObjects.cxx.

References pma::SortedObjectBase::prev, and pma::SortedObjectBase::RemoveNext().

{
  while (next_vector.size())
    RemoveNext(next_vector.front());
  if (prev) prev->RemoveNext(this);
}

double pma::Node3D::EndPtCos2Transverse ( ) const

private

Definition at line 361 of file PmaNode3D.cxx.

References fPoint3D, pma::SortedObjectBase::Next(), pma::SortedObjectBase::next, pma::SortedObjectBase::Prev(), and pma::SortedObjectBase::prev.

Referenced by PenaltyInWirePlane().

{
  if (prev && next) {
    auto const& vStart = static_cast<pma::Node3D*>(prev->Prev())->fPoint3D;
    auto const& vStop = static_cast<pma::Node3D*>(next->Next())->fPoint3D;

    double dy = vStop.X() - vStart.X();
    double dz = vStop.Z() - vStart.Z();
    double len2 = dy * dy + dz * dz;
    double cosine2 = 0.0;
    if (len2 > 0.0) cosine2 = dz * dz / len2;
    return cosine2;
  }
  else
    return 0.0;
}

std::vector< pma::Track3D * > pma::Node3D::GetBranches

(

)

const

Definition at line 447 of file PmaNode3D.cxx.

References pma::SortedBranchBase::Next(), pma::SortedBranchBase::NextCount(), and pma::Segment3D::Parent().

Referenced by pma::VtxCandidate::JoinTracks().

{
  std::vector<pma::Track3D*> branches;
  if (NextCount()) {
    branches.reserve(NextCount());
    for (size_t i = 0; i < NextCount(); ++i) {
      pma::Segment3D* seg = static_cast<pma::Segment3D*>(Next(i));
      branches.push_back(seg->Parent());
    }
  }
  return branches;
}

pma::Vector3D pma::Node3D::GetDirection3D ( void  ) const

overridevirtual

Get 3D direction cosines of the next segment, or previous segment if this is the last node.

Implements pma::Element3D.

Definition at line 201 of file PmaNode3D.cxx.

References pma::Element3D::GetDirection3D(), pma::SortedObjectBase::next, and pma::SortedObjectBase::prev.

{
  pma::Element3D* seg = 0;
  if (next) { seg = dynamic_cast<pma::Element3D*>(next); }
  else if (prev) {
    seg = dynamic_cast<pma::Element3D*>(prev);
  }
  else {
    throw cet::exception("Node3D") << "Isolated vertex." << std::endl;
  }

  if (seg) { return seg->GetDirection3D(); }
  else {
    throw cet::exception("Node3D") << "Corrupted vertex." << std::endl;
  }
}

double pma::Node3D::GetDistance2To ( const TVector3 &  p3d ) const

overridevirtual

Distance [cm] from the 3D point to the point 3D.

Implements pma::Element3D.

Definition at line 176 of file PmaNode3D.cxx.

References pma::Dist2(), and fPoint3D.

Referenced by pma::Track3D::MakeFastProjection().

{
  return pma::Dist2(fPoint3D, p3d);
}

double pma::Node3D::GetDistance2To ( const TVector2 &  p2d, unsigned int  view  ) const

overridevirtual

Distance [cm] from the 2D point to the object's 2D projection in one of wire views.

Implements pma::Element3D.

Definition at line 181 of file PmaNode3D.cxx.

References pma::Dist2(), and fProj2D.

{
  return pma::Dist2(fProj2D[view], p2d);
}

double pma::Node3D::GetDistToWall

(

)

const

Definition at line 84 of file PmaNode3D.cxx.

References d, fMaxX, fMaxY, fMaxZ, fMinX, fMinY, fMinZ, and fPoint3D.

Referenced by pma::PMAlgTracker::validate().

{
  double d, dmin = fPoint3D.X() - fMinX;
  d = fMaxX - fPoint3D.X();
  if (d < dmin) dmin = d;

  d = fPoint3D.Y() - fMinY;
  if (d < dmin) dmin = d;
  d = fMaxY - fPoint3D.Y();
  if (d < dmin) dmin = d;

  d = fPoint3D.Z() - fMinZ;
  if (d < dmin) dmin = d;
  d = fMaxZ - fPoint3D.Z();
  if (d < dmin) dmin = d;

  return dmin;
}

double pma::Node3D::GetDriftShift ( ) const

inline

Definition at line 118 of file PmaNode3D.h.

Referenced by pma::Track3D::Split().

{ return fDriftOffset; }

double pma::Node3D::GetObjFunction	(	float	penaltyValue,
		float	endSegWeight
	)		const

Objective function minimized during oprimization.

Definition at line 557 of file PmaNode3D.cxx.

References Mse(), Penalty(), and PenaltyInWirePlane().

Referenced by MakeGradient(), and StepWithGradient().

{
  return Mse() + penaltyValue * (Penalty(endSegWeight) + PenaltyInWirePlane());
}

TVector3 pma::Node3D::GetUnconstrainedProj3D ( const TVector2 &  , unsigned int    ) const

inlineoverridevirtual

In case of a node it is simply 3D position of the node.

Implements pma::Element3D.

Definition at line 83 of file PmaNode3D.h.

{ return fPoint3D; }

bool pma::Element3D::HasHit ( const pma::Hit3D *  h ) const

inlineinherited

Definition at line 57 of file PmaElement3D.h.

References pma::Element3D::fAssignedHits.

  {
    for (const auto a : fAssignedHits) {
      if (h == a) return true;
    }
    return false;
  }

pma::Hit3D& pma::Element3D::Hit ( size_t  index )

inlineinherited

Definition at line 65 of file PmaElement3D.h.

References pma::Element3D::fAssignedHits.

Referenced by pma::Track3D::AddNode(), pma::Track3D::DisableSingleViewEnds(), and pma::Track3D::SortHits().

{ return *(fAssignedHits[index]); }

const std::vector<pma::Hit3D*>& pma::Element3D::Hits ( void  ) const

inlineinherited

Definition at line 55 of file PmaElement3D.h.

References pma::Element3D::fAssignedHits.

{ return fAssignedHits; }

double pma::Element3D::HitsRadius3D ( unsigned int  view ) const

inherited

Definition at line 209 of file PmaElement3D.cxx.

References pma::Dist2(), pma::Element3D::fAssignedHits, and pma::Element3D::fTPC.

Referenced by pma::Element3D::NThisHits().

{
  if (fTPC < 0) {
    if (!fAssignedHits.empty())
      mf::LogWarning("pma::Element3D") << "Hits assigned to TPC-crossing element.";
    return 0.0F;
  }

  TVector3 mean3D(0, 0, 0);
  size_t nHits = 0;
  for (auto h : fAssignedHits)
    if (h->View2D() == view) {
      mean3D += h->Point3D();
      nHits++;
    }
  if (!nHits) return 0.0;
  mean3D *= (1.0 / nHits);

  double r2, maxR2 = 0.0;
  for (auto h : fAssignedHits)
    if (h->View2D() == view) {
      r2 = pma::Dist2(h->Point3D(), mean3D);
      if (r2 > maxR2) maxR2 = r2;
    }
  return sqrt(maxR2);
}

bool pma::Node3D::IsBranching

(

)

const

Belongs to more than one track?

Definition at line 422 of file PmaNode3D.cxx.

References pma::SortedBranchBase::Next(), pma::SortedBranchBase::NextCount(), pma::Segment3D::Parent(), and pma::SortedObjectBase::prev.

Referenced by pma::PMAlgVertexing::getKinks(), and pma::PMAlgVertexing::getVertices().

{
  size_t nnext = NextCount();
  if (nnext > 1) return true; // 1 trk -> vtx -> n*trk

  if (prev && nnext) {
    pma::Segment3D* segPrev = static_cast<pma::Segment3D*>(prev);
    pma::Segment3D* segNext = static_cast<pma::Segment3D*>(Next(0));
    if (segNext->Parent() != segPrev->Parent()) // 1 trk -> vtx -> 1 trk
      return true;
  }
  return false;
}

virtual bool pma::SortedObjectBase::IsFirst ( void  ) const

inlinevirtualinherited

Definition at line 39 of file SortedObjects.h.

References pma::SortedObjectBase::prev.

{ return !prev; }

bool pma::Element3D::IsFrozen ( void  ) const

inlineinherited

Check if the vertex 3D position is fixed.

Definition at line 105 of file PmaElement3D.h.

References pma::Element3D::fFrozen.

Referenced by pma::Track3D::AddNode(), and pma::ProjectionMatchingAlg::mergeTracks().

{ return fFrozen; }

virtual bool pma::SortedBranchBase::IsLast ( void  ) const

inlinevirtualinherited

Reimplemented from pma::SortedObjectBase.

Definition at line 85 of file SortedObjects.h.

{ return !(next_vector.size()); }

bool pma::Node3D::IsTPCEdge

(

)

const

Is the first/last in this TPC?

Definition at line 436 of file PmaNode3D.cxx.

References pma::SortedBranchBase::Next(), pma::SortedBranchBase::NextCount(), pma::SortedObjectBase::prev, and pma::Element3D::TPC().

{
  if (prev && (NextCount() == 1)) {
    pma::Segment3D* segPrev = static_cast<pma::Segment3D*>(prev);
    pma::Segment3D* segNext = static_cast<pma::Segment3D*>(Next(0));

    if ((segPrev->TPC() < 0) || (segNext->TPC() < 0)) return true;
  }
  return false;
}

bool pma::Node3D::IsVertex ( ) const

inline

Check fIsVertex flag.

Definition at line 62 of file PmaNode3D.h.

Referenced by pma::PMAlgVertexing::getKinks().

{ return fIsVertex; }

double pma::Element3D::Length ( void  ) const

inlineinherited

Definition at line 53 of file PmaElement3D.h.

References pma::Element3D::Length2().

Referenced by pma::VtxCandidate::Add(), pma::Track3D::AddNode(), pma::VtxCandidate::Compute(), pma::ProjectionMatchingAlg::guideEndpoints(), Pi(), and pma::Track3D::ShiftEndsToHits().

{ return sqrt(Length2()); }

double pma::Node3D::Length2 ( void  ) const

overridevirtual

Squared sum of half-lengths of connected 3D segments (used in the vertex position optimization).

Implements pma::Element3D.

Definition at line 292 of file PmaNode3D.cxx.

References pma::SortedObjectBase::next, and pma::SortedObjectBase::prev.

Referenced by Pi(), and PiInWirePlane().

{
  double l = 0.0;
  if (next) l += (static_cast<pma::Segment3D*>(next))->Length();
  if (prev) l += (static_cast<pma::Segment3D*>(prev))->Length();

  if (next && prev)
    return 0.25 * l * l;
  else
    return l * l;
}

bool pma::Node3D::LimitPoint3D ( )

private

Returns true if node position was trimmed to its TPC volume + fMargin.

Definition at line 123 of file PmaNode3D.cxx.

References fMargin, fMaxX, fMaxY, fMaxZ, fMinX, fMinY, fMinZ, and fPoint3D.

Referenced by SetPoint3D().

{
  bool trimmed = false;

  if (fPoint3D.X() < fMinX - fMargin) {
    fPoint3D.SetX(fMinX - fMargin);
    trimmed = true;
  }
  if (fPoint3D.X() > fMaxX + fMargin) {
    fPoint3D.SetX(fMaxX + fMargin);
    trimmed = true;
  }

  if (fPoint3D.Y() < fMinY - fMargin) {
    fPoint3D.SetY(fMinY - fMargin);
    trimmed = true;
  }
  if (fPoint3D.Y() > fMaxY + fMargin) {
    fPoint3D.SetY(fMaxY + fMargin);
    trimmed = true;
  }

  if (fPoint3D.Z() < fMinZ - fMargin) {
    fPoint3D.SetZ(fMinZ - fMargin);
    trimmed = true;
  }
  if (fPoint3D.Z() > fMaxZ + fMargin) {
    fPoint3D.SetZ(fMaxZ + fMargin);
    trimmed = true;
  }

  return trimmed;
}

double pma::Node3D::MakeGradient ( float  penaltyValue, float  endSegWeight  )

private

Definition at line 562 of file PmaNode3D.cxx.

References E, fGradFixed, fGradient, fPoint3D, GetObjFunction(), pma::Segment3D::Length2(), pma::SortedObjectBase::next, pma::SortedObjectBase::prev, SetPoint3D(), and tmp.

Referenced by StepWithGradient().

{
  double l1 = 0.0, l2 = 0.0, minLength2 = 0.0;
  TVector3 tmp(fPoint3D), gpoint(fPoint3D);

  pma::Segment3D* seg;
  if (prev) {
    seg = static_cast<pma::Segment3D*>(prev);
    l1 = seg->Length2();
  }
  if (next) {
    seg = static_cast<pma::Segment3D*>(next);
    l2 = seg->Length2();
  }
  if ((l1 > 0.01) && (l1 < l2))
    minLength2 = l1;
  else if ((l2 > 0.01) && (l2 < l1))
    minLength2 = l2;
  else
    minLength2 = 0.01;

  double dxi = 0.001 * sqrt(minLength2);

  if (dxi < 6.0E-37) return 0.0;

  double gi, g0, gz;
  gz = g0 = GetObjFunction(penaltyValue, endSegWeight);

  if (!fGradFixed[0]) // gradX
  {
    gpoint[0] = tmp[0] + dxi;
    SetPoint3D(gpoint);
    gi = GetObjFunction(penaltyValue, endSegWeight);
    fGradient[0] = (g0 - gi) / dxi;

    gpoint[0] = tmp[0] - dxi;
    SetPoint3D(gpoint);
    gi = GetObjFunction(penaltyValue, endSegWeight);
    fGradient[0] = 0.5 * (fGradient[0] + (gi - g0) / dxi);

    gpoint[0] = tmp[0];
  }

  if (!fGradFixed[1]) // gradY
  {
    gpoint[1] = tmp[1] + dxi;
    SetPoint3D(gpoint);
    gi = GetObjFunction(penaltyValue, endSegWeight);
    fGradient[1] = (g0 - gi) / dxi;

    gpoint[1] = tmp[1] - dxi;
    SetPoint3D(gpoint);
    gi = GetObjFunction(penaltyValue, endSegWeight);
    fGradient[1] = 0.5 * (fGradient[1] + (gi - g0) / dxi);

    gpoint[1] = tmp[1];
  }

  if (!fGradFixed[2]) // gradZ
  {
    gpoint[2] = tmp[2] + dxi;
    SetPoint3D(gpoint);
    gi = GetObjFunction(penaltyValue, endSegWeight);
    fGradient[2] = (gz - gi) / dxi;

    gpoint[2] = tmp[2] - dxi;
    SetPoint3D(gpoint);
    gi = GetObjFunction(penaltyValue, endSegWeight);
    fGradient[2] = 0.5 * (fGradient[2] + (gi - gz) / dxi);

    gpoint[2] = tmp[2];
  }

  SetPoint3D(tmp);
  if (fGradient.Mag2() < 6.0E-37) return 0.0;

  return g0;
}

double pma::Node3D::Mse ( ) const

private

Definition at line 535 of file PmaNode3D.cxx.

References pma::SortedBranchBase::Next(), pma::SortedBranchBase::NextCount(), pma::Element3D::NPrecalcEnabledHits(), pma::SortedObjectBase::prev, and pma::Element3D::SumDist2().

Referenced by GetObjFunction().

{
  unsigned int nhits = NPrecalcEnabledHits(); //NEnabledHits();
  double mse = SumDist2();

  pma::Segment3D* seg;
  for (unsigned int i = 0; i < NextCount(); i++) {
    seg = static_cast<pma::Segment3D*>(Next(i));
    nhits += seg->NPrecalcEnabledHits(); //NEnabledHits();
    mse += seg->SumDist2();
  }
  if (prev) {
    seg = static_cast<pma::Segment3D*>(prev);
    nhits += seg->NPrecalcEnabledHits(); //NEnabledHits();
    mse += seg->SumDist2();
  }
  if (!nhits)
    return 0.0;
  else
    return mse / nhits;
}

size_t pma::Element3D::NEnabledHits ( unsigned int  view = geo::kUnknown ) const

inherited

Definition at line 34 of file PmaElement3D.cxx.

References pma::Element3D::fAssignedHits, geo::kUnknown, and n.

Referenced by pma::Track3D::AddNode(), and pma::Element3D::NHits().

{
  size_t n = 0;
  for (size_t i = 0; i < fAssignedHits.size(); i++)
    if (fAssignedHits[i]->IsEnabled() &&
        ((view == geo::kUnknown) || (view == fAssignedHits[i]->View2D())))
      n++;
  return n;
}

virtual pma::SortedObjectBase* pma::SortedBranchBase::Next ( unsigned int  index = 0 ) const

inlinevirtualinherited

Reimplemented from pma::SortedObjectBase.

Definition at line 77 of file SortedObjects.h.

Referenced by pma::Track3D::ApplyDriftShiftInTree(), pma::Track3D::AttachBackTo(), pma::Track3D::AttachBackToSameTPC(), pma::Track3D::AttachTo(), pma::Track3D::AttachToSameTPC(), ClearAssigned(), GetBranches(), pma::Track3D::GetNearestTrkInTree(), pma::Track3D::GetObjFnInTree(), pma::TrkCandidateColl::getTreeCopy(), IsBranching(), IsTPCEdge(), pma::VtxCandidate::JoinTracks(), pma::Track3D::MakeProjectionInTree(), pma::PMAlgVertexing::mergeBrokenTracks(), Mse(), pma::Track3D::NextSegment(), Penalty(), Pi(), PiInWirePlane(), pma::Track3D::ReassignHitsInTree(), pma::TrkCandidateColl::setParentDaughterConnections(), pma::TrkCandidateColl::setTreeId(), pma::Track3D::SortHitsInTree(), pma::Track3D::Split(), pma::Track3D::TuneSinglePass(), and pma::Track3D::UpdateParamsInTree().

  {
    if (next_vector.size())
      return next_vector[index];
    else
      return 0;
  }

virtual unsigned int pma::SortedBranchBase::NextCount ( void  ) const

inlinevirtualinherited

Reimplemented from pma::SortedObjectBase.

Definition at line 84 of file SortedObjects.h.

Referenced by pma::Track3D::ApplyDriftShiftInTree(), pma::Track3D::AttachBackTo(), pma::Track3D::AttachBackToSameTPC(), pma::Track3D::AttachTo(), pma::Track3D::AttachToSameTPC(), pma::Track3D::CleanupTails(), ClearAssigned(), GetBranches(), pma::Track3D::GetNearestTrkInTree(), pma::Track3D::GetObjFnInTree(), pma::TrkCandidateColl::getTreeCopy(), IsBranching(), IsTPCEdge(), pma::VtxCandidate::JoinTracks(), pma::Track3D::MakeProjectionInTree(), pma::PMAlgVertexing::mergeBrokenTracks(), Mse(), pma::Track3D::NextSegment(), Penalty(), Pi(), PiInWirePlane(), pma::Track3D::ReassignHitsInTree(), pma::Track3D::RemoveNode(), pma::TrkCandidateColl::setParentDaughterConnections(), pma::TrkCandidateColl::setTreeId(), pma::Track3D::ShiftEndsToHits(), pma::Track3D::SortHitsInTree(), pma::Track3D::TuneSinglePass(), and pma::Track3D::UpdateParamsInTree().

{ return next_vector.size(); }

size_t pma::Element3D::NHits ( void  ) const

inlineinherited

Definition at line 76 of file PmaElement3D.h.

References pma::Element3D::fAssignedHits, geo::kUnknown, and pma::Element3D::NEnabledHits().

Referenced by pma::Track3D::CleanupTails(), pma::Track3D::DisableSingleViewEnds(), PiInWirePlane(), pma::Element3D::SelectRndHits(), and pma::Track3D::SortHits().

{ return fAssignedHits.size(); }

unsigned int pma::Element3D::NHits ( unsigned int  view ) const

inlineinherited

Definition at line 99 of file PmaElement3D.h.

References pma::Element3D::fNHits.

{ return fNHits[view]; }

size_t pma::Element3D::NPoints ( void  ) const

inlineinherited

Definition at line 81 of file PmaElement3D.h.

References pma::Element3D::fAssignedPoints.

Referenced by pma::Track3D::DisableSingleViewEnds().

{ return fAssignedPoints.size(); }

size_t pma::Element3D::NPrecalcEnabledHits ( void  ) const

inlineinherited

Definition at line 78 of file PmaElement3D.h.

References pma::Element3D::fNThisHitsEnabledAll.

Referenced by Mse().

{ return fNThisHitsEnabledAll; }

unsigned int pma::Element3D::NThisHits ( unsigned int  view ) const

inlineinherited

Definition at line 100 of file PmaElement3D.h.

References pma::Element3D::fNThisHits, and pma::Element3D::HitsRadius3D().

{ return fNThisHits[view]; }

static float pma::Element3D::OptFactor ( unsigned int  view )

inlinestaticinherited

Definition at line 112 of file PmaElement3D.h.

References pma::Element3D::fOptFactors.

Referenced by pma::Element3D::SumDist2(), pma::Segment3D::SumDist2Hits(), and SumDist2Hits().

{ return fOptFactors[view]; }

void pma::Node3D::Optimize	(	float	penaltyValue,
		float	endSegWeight
	)

Optimize vertex 3D position with given penalty on connected segments angle and weight assigned to the outermost segments. Only MSE is used in case of branching nodes.

Definition at line 811 of file PmaNode3D.cxx.

References pma::Element3D::fFrozen, and StepWithGradient().

Referenced by pma::Track3D::TuneSinglePass().

{
  if (!fFrozen) {
    double dg = StepWithGradient(0.1F, 0.002F, penaltyValue, endSegWeight);
    if (dg > 0.01) dg = StepWithGradient(0.03F, 0.0001F, penaltyValue, endSegWeight);
    if (dg > 0.0) dg = StepWithGradient(0.03F, 0.0001F, penaltyValue, endSegWeight);
  }
}

double pma::Node3D::Penalty ( float  endSegWeight ) const

private

Definition at line 516 of file PmaNode3D.cxx.

References pma::SortedObjectBase::Next(), pma::SortedBranchBase::Next(), pma::SortedBranchBase::NextCount(), Pi(), pma::SortedObjectBase::Prev(), and pma::SortedObjectBase::prev.

Referenced by GetObjFunction().

{
  unsigned int nseg = 1;
  double penalty = Pi(endSegWeight, true);

  pma::Node3D* v;
  for (unsigned int i = 0; i < NextCount(); i++) {
    v = static_cast<pma::Node3D*>(Next(i)->Next());
    penalty += v->Pi(endSegWeight, false);
    nseg++;
  }
  if (prev) {
    v = static_cast<pma::Node3D*>(prev->Prev());
    penalty += v->Pi(endSegWeight, false);
    nseg++;
  }
  return penalty / nseg;
}

double pma::Node3D::PenaltyInWirePlane ( ) const

private

Definition at line 399 of file PmaNode3D.cxx.

References EndPtCos2Transverse(), fIsVertex, pma::SortedObjectBase::Next(), pma::SortedObjectBase::next, PiInWirePlane(), pma::SortedObjectBase::Prev(), and pma::SortedObjectBase::prev.

Referenced by GetObjFunction().

{
  if (fIsVertex) return 0.0;

  unsigned int nseg = 1;
  double penalty = PiInWirePlane();
  pma::Node3D* v;
  if (next) {
    v = static_cast<pma::Node3D*>(next->Next());
    penalty += v->PiInWirePlane();
    nseg++;
  }
  if (prev) {
    v = static_cast<pma::Node3D*>(prev->Prev());
    penalty += v->PiInWirePlane();
    nseg++;
  }
  if (penalty > 0.0)
    return pow(EndPtCos2Transverse(), 10) * penalty / nseg;
  else
    return 0.0;
}

double pma::Node3D::Pi ( float  endSegWeight, bool  doAsymm  ) const

private

Definition at line 460 of file PmaNode3D.cxx.

References pma::Element3D::fHitsRadius, fIsVertex, pma::Element3D::Length(), pma::Segment3D::Length2(), Length2(), pma::SortedObjectBase::Next(), pma::SortedObjectBase::next, pma::SortedBranchBase::Next(), pma::SortedObjectBase::NextCount(), pma::SortedBranchBase::NextCount(), pma::SortedObjectBase::Prev(), pma::SortedObjectBase::prev, scale, SegmentCos(), and pma::Element3D::TPC().

Referenced by Penalty().

{
  if (fIsVertex) return 0.0;

  if (prev && NextCount()) {
    pma::Segment3D* segPrev = static_cast<pma::Segment3D*>(prev);
    pma::Segment3D* segNext = static_cast<pma::Segment3D*>(Next(0));

    double scale = 1.0;
    if ((segPrev->TPC() < 0) || (segNext->TPC() < 0))
      scale = 0.5; // lower penalty on segments between tpc's

    double segCos = SegmentCos();

    double lAsymmFactor = 0.0;
    if (doAsymm) {
      double lPrev = segPrev->Length();
      double lNext = segNext->Length();
      double lSum = lPrev + lNext;
      if (lSum > 0.1) {
        double lAsymm = (1.0 - segCos) * (lPrev - lNext) / lSum;
        lAsymmFactor = 0.05 * lAsymm * lAsymm;
      }
    }

    if (fHitsRadius > 0.0F)
      return scale * (1.0 + segCos + lAsymmFactor) * fHitsRadius * fHitsRadius;
    else
      return scale * (1.0 + segCos + lAsymmFactor) * Length2();
  }
  else {
    double pi_result = 0.0;
    unsigned int nSeg = 0;
    pma::Segment3D* seg = 0;
    if (prev) {
      seg = static_cast<pma::Segment3D*>(prev);

      SortedObjectBase* prevVtx = seg->Prev();
      if (prevVtx->Prev()) nSeg++;
      nSeg += prevVtx->NextCount();
    }
    else if (next) {
      seg = static_cast<pma::Segment3D*>(next);

      SortedObjectBase* nextVtx = seg->Next();
      nSeg += nextVtx->NextCount() + 1;
    }
    else {
      mf::LogWarning("pma::Node3D") << "pma::Node3D::Pi(): an isolated vertex?";
      return 0.0;
    }
    if (nSeg == 1) pi_result = endSegWeight * seg->Length2();
    return pi_result;
  }
}

double pma::Node3D::PiInWirePlane ( ) const

private

Definition at line 378 of file PmaNode3D.cxx.

References pma::Element3D::fHitsRadius, geo::kU, Length2(), pma::SortedBranchBase::Next(), pma::SortedBranchBase::NextCount(), pma::Element3D::NHits(), pma::SortedObjectBase::prev, and SegmentCosWirePlane().

Referenced by PenaltyInWirePlane().

{
  if (prev && NextCount()) {
    pma::Segment3D* seg0 = dynamic_cast<pma::Segment3D*>(prev);
    pma::Segment3D* seg1 = dynamic_cast<pma::Segment3D*>(Next(0));
    unsigned int nInd1 = NHits(geo::kU) + seg0->NHits(geo::kU) + seg1->NHits(geo::kU);

    if (fHitsRadius > 0.0F)
      return (1.0 + SegmentCosWirePlane()) * fHitsRadius * fHitsRadius / (4 * nInd1 + 1.0);
    else
      return (1.0 + SegmentCosWirePlane()) * Length2() / (4 * nInd1 + 1.0);
  }
  else
    return 0.0;
}

TVector3 const& pma::Node3D::Point3D ( ) const

inline

Definition at line 41 of file PmaNode3D.h.

Referenced by pma::Track3D::AttachTo(), pma::VtxCandidate::Compute(), pma::Segment3D::GetDirection3D(), pma::Segment3D::GetDistance2To(), pma::PMAlgVertexing::getKinks(), pma::Segment3D::GetProjection(), pma::Segment3D::GetUnconstrainedProj3D(), pma::VtxCandidate::JoinTracks(), pma::VtxCandidate::MaxAngle(), pma::Segment3D::SetProjection(), SetProjection(), and pma::Track3D::Split().

{ return fPoint3D; }

virtual pma::SortedObjectBase* pma::SortedObjectBase::Prev ( void  ) const

inlinevirtualinherited

Definition at line 42 of file SortedObjects.h.

References pma::SortedObjectBase::prev.

Referenced by pma::Track3D::AttachBackTo(), pma::Track3D::AttachBackToOtherTPC(), pma::Track3D::AttachBackToSameTPC(), pma::Track3D::AttachTo(), pma::Track3D::AttachToSameTPC(), pma::Track3D::CanFlip(), ClearAssigned(), EndPtCos2Transverse(), pma::Track3D::Flip(), pma::VtxCandidate::JoinTracks(), pma::Track3D::MakeFastProjection(), pma::PMAlgVertexing::mergeBrokenTracks(), Penalty(), PenaltyInWirePlane(), Pi(), pma::Track3D::PrevSegment(), pma::Track3D::RemoveNode(), SegmentCos(), SegmentCosTransverse(), SegmentCosWirePlane(), pma::TrkCandidateColl::setParentDaughterConnections(), SetProjection(), pma::TrkCandidateColl::setTreeOriginAtFront(), pma::Track3D::ShiftEndsToHits(), pma::Track3D::Split(), pma::Segment3D::Start(), and pma::Element3D::UpdateHitParams().

{ return prev; }

TVector2 const& pma::Node3D::Projection2D ( unsigned int  view ) const

inline

Definition at line 47 of file PmaNode3D.h.

Referenced by pma::Segment3D::GetDistance2To(), pma::Segment3D::GetProjection(), pma::Segment3D::GetUnconstrainedProj3D(), pma::Segment3D::SetProjection(), SetProjection(), and pma::Segment3D::SumDist2Hits().

{ return fProj2D[view]; }

TVector3 const& pma::Element3D::ReferencePoint ( size_t  index ) const

inlineinherited

Definition at line 80 of file PmaElement3D.h.

References pma::Element3D::fAssignedPoints.

{ return *(fAssignedPoints[index]); }

void pma::Element3D::RemoveHitAt ( size_t  index )

inlineinherited

Definition at line 66 of file PmaElement3D.h.

References pma::Element3D::fAssignedHits.

Referenced by pma::Track3D::MakeFastProjection().

  {
    if (index < fAssignedHits.size()) fAssignedHits.erase(fAssignedHits.begin() + index);
  }

int pma::SortedBranchBase::RemoveNext ( pma::SortedObjectBase *  nextElement )

virtualinherited

Reimplemented from pma::SortedObjectBase.

Definition at line 102 of file SortedObjects.cxx.

References pma::SortedObjectBase::next, and pma::SortedObjectBase::prev.

Referenced by pma::Track3D::AttachBackToSameTPC(), pma::Track3D::AttachToSameTPC(), and pma::Track3D::Split().

{
  if (!nextElement || (nextElement == this)) return -1;

  int index = -1;
  for (unsigned int i = 0; i < next_vector.size(); i++) {
    if (next_vector[i] == nextElement) {
      index = i;
      break;
    }
  }
  if (index >= 0) {
    if (next_vector[index]->prev != this) {
      mf::LogError("pma::SortedBranchBase") << "Object structure is broken.";
      return -1;
    }
    next_vector[index]->prev = 0;
    next_vector.erase(next_vector.begin() + index);

    if (next_vector.size() > 0)
      next = next_vector.back();
    else
      next = 0;
  }
  return index;
}

bool pma::Node3D::SameTPC	(	const TVector3 &	p3d,
		float	margin = 0.0F
	)		const

Check if p3d is in the same TPC as the node.

Definition at line 103 of file PmaNode3D.cxx.

References fMaxX, fMaxY, fMaxZ, fMinX, fMinY, and fMinZ.

Referenced by pma::Track3D::GetUnconstrainedProj3D(), pma::ProjectionMatchingAlg::isContained(), pma::ProjectionMatchingAlg::validate(), pma::ProjectionMatchingAlg::validate_on_adc(), and pma::ProjectionMatchingAlg::validate_on_adc_test().

{
  if (((fMinX - margin) <= p3d.X()) && (p3d.X() <= (fMaxX + margin)) &&
      ((fMinY - margin) <= p3d.Y()) && (p3d.Y() <= (fMaxY + margin)) &&
      ((fMinZ - margin) <= p3d.Z()) && (p3d.Z() <= (fMaxZ + margin)))
    return true;
  else
    return false;
}

bool pma::Node3D::SameTPC	(	const pma::Vector3D &	p3d,
		float	margin = 0.0F
	)		const

Definition at line 113 of file PmaNode3D.cxx.

References fMaxX, fMaxY, fMaxZ, fMinX, fMinY, and fMinZ.

{
  if (((fMinX - margin) <= p3d.X()) && (p3d.X() <= (fMaxX + margin)) &&
      ((fMinY - margin) <= p3d.Y()) && (p3d.Y() <= (fMaxY + margin)) &&
      ((fMinZ - margin) <= p3d.Z()) && (p3d.Z() <= (fMaxZ + margin)))
    return true;
  else
    return false;
}

double pma::Node3D::SegmentCos

(

)

const

Cosine of 3D angle between connected segments.

Definition at line 304 of file PmaNode3D.cxx.

References fPoint3D, pma::SortedObjectBase::Next(), pma::SortedObjectBase::next, pma::SortedObjectBase::Prev(), and pma::SortedObjectBase::prev.

Referenced by Pi().

{
  if (prev && next) {
    auto const& vStop1 = static_cast<pma::Node3D*>(prev->Prev())->fPoint3D;
    auto const& vStop2 = static_cast<pma::Node3D*>(next->Next())->fPoint3D;
    pma::Vector3D v1(
      vStop1.X() - fPoint3D.X(), vStop1.Y() - fPoint3D.Y(), vStop1.Z() - fPoint3D.Z());
    pma::Vector3D v2(
      vStop2.X() - fPoint3D.X(), vStop2.Y() - fPoint3D.Y(), vStop2.Z() - fPoint3D.Z());
    double mag = sqrt(v1.Mag2() * v2.Mag2());
    double cosine = 0.0;
    if (mag != 0.0) cosine = v1.Dot(v2) / mag;
    return cosine;
  }
  else {
    mf::LogError("pma::Node3D") << "pma::Node3D::SegmentCos(): neighbours not initialized.";
    return -1.0;
  }
}

double pma::Node3D::SegmentCosTransverse

(

)

const

Cosine of 2D angle (in horizontal plane, parallel to drift) between connected segments. Should be changed / generalized for horizontal wire planes (e.g. 2-phase LAr).

Definition at line 342 of file PmaNode3D.cxx.

References fPoint3D, pma::SortedObjectBase::Next(), pma::SortedObjectBase::next, pma::SortedObjectBase::Prev(), and pma::SortedObjectBase::prev.

{
  if (prev && next) {
    auto const& vStop1 = static_cast<pma::Node3D*>(prev->Prev())->fPoint3D;
    auto const& vStop2 = static_cast<pma::Node3D*>(next->Next())->fPoint3D;
    pma::Vector2D v1(vStop1.X() - fPoint3D.X(), vStop1.Z() - fPoint3D.Z());
    pma::Vector2D v2(vStop2.X() - fPoint3D.X(), vStop2.Z() - fPoint3D.Z());
    double mag = sqrt(v1.Mag2() * v2.Mag2());
    double cosine = 0.0;
    if (mag != 0.0) cosine = v1.Dot(v2) / mag;
    return cosine;
  }
  else {
    mf::LogError("pma::Node3D") << "pma::Node3D::SegmentCosZY(): neighbours not initialized.";
    return -1.0;
  }
}

double pma::Node3D::SegmentCosWirePlane

(

)

const

Cosine of 2D angle (in plane parallel to wire planes) between connected segments. Should be changed / generalized for horizontal wire planes (e.g. 2-phase LAr).

Definition at line 324 of file PmaNode3D.cxx.

References fPoint3D, pma::SortedObjectBase::Next(), pma::SortedObjectBase::next, pma::SortedObjectBase::Prev(), and pma::SortedObjectBase::prev.

Referenced by PiInWirePlane().

{
  if (prev && next) {
    auto const& vStop1 = static_cast<pma::Node3D*>(prev->Prev())->fPoint3D;
    auto const& vStop2 = static_cast<pma::Node3D*>(next->Next())->fPoint3D;
    pma::Vector2D v1(vStop1.Y() - fPoint3D.Y(), vStop1.Z() - fPoint3D.Z());
    pma::Vector2D v2(vStop2.Y() - fPoint3D.Y(), vStop2.Z() - fPoint3D.Z());
    double mag = sqrt(v1.Mag2() * v2.Mag2());
    double cosine = 0.0;
    if (mag != 0.0) cosine = v1.Dot(v2) / mag;
    return cosine;
  }
  else {
    mf::LogError("pma::Node3D") << "pma::Node3D::SegmentCosZX(): neighbours not initialized.";
    return -1.0;
  }
}

bool pma::Element3D::SelectAllHits ( void  )

inherited

Definition at line 284 of file PmaElement3D.cxx.

References pma::Element3D::fAssignedHits.

Referenced by pma::Track3D::SelectHits(), pma::Element3D::SelectRndHits(), and pma::Element3D::SetFrozen().

{
  bool changed = false;
  for (auto h : fAssignedHits) {
    changed |= !(h->IsEnabled());
    h->SetEnabled(true);
  }
  return changed;
}

bool pma::Element3D::SelectRndHits ( size_t  nmax_per_view )

inherited

Definition at line 236 of file PmaElement3D.cxx.

References pma::Element3D::fAssignedHits, pma::Element3D::NHits(), and pma::Element3D::SelectAllHits().

Referenced by pma::Element3D::SetFrozen().

{
  if (!nmax_per_view) { return SelectAllHits(); }

  size_t nhits[3];
  for (size_t i = 0; i < 3; ++i)
    nhits[i] = NHits(i);

  int m[3], count[3];
  bool state[3];
  for (size_t i = 0; i < 3; ++i) {
    if (nhits[i] >= 2 * nmax_per_view) {
      m[i] = nhits[i] / nmax_per_view;
      state[i] = true;
    }
    else if (nhits[i] > nmax_per_view) {
      m[i] = nhits[i] / (nhits[i] - nmax_per_view);
      state[i] = false;
    }
    else {
      m[i] = 0;
      state[i] = false;
    }

    count[i] = 0;
  }

  bool b, changed = false;
  for (auto h : fAssignedHits) {
    b = h->IsEnabled();

    size_t view = h->View2D();
    if (m[view]) {
      if (count[view] % m[view] == 0)
        h->SetEnabled(state[view]);
      else
        h->SetEnabled(!(state[view]));

      ++count[view];
    }
    else
      h->SetEnabled(true);

    changed |= (b != h->IsEnabled());
  }
  return changed;
}

void pma::Element3D::SetFrozen ( bool  state )

inlineinherited

Fix / relese vertex 3D position.

Definition at line 107 of file PmaElement3D.h.

References pma::Element3D::fFrozen, pma::Element3D::SelectAllHits(), and pma::Element3D::SelectRndHits().

{ fFrozen = state; }

static void pma::Node3D::SetMargin ( double  m )

inlinestatic

Set allowed node position margin around TPC.

Definition at line 121 of file PmaNode3D.h.

References weight.

Referenced by pma::ProjectionMatchingAlg::ProjectionMatchingAlg().

  {
    if (m >= 0.0) fMargin = m;
  }

static void pma::Element3D::SetOptFactor ( unsigned int  view, float  value  )

inlinestaticinherited

Definition at line 113 of file PmaElement3D.h.

References pma::Element3D::Element3D(), pma::Element3D::fOptFactors, and value.

Referenced by pma::ProjectionMatchingAlg::ProjectionMatchingAlg().

{ fOptFactors[view] = value; }

bool pma::Node3D::SetPoint3D

(

const TVector3 &

p3d

)

Returns true if the new position was accepted; returns false if the new position was trimmed to fit insite TPC volume + fMargin.

Definition at line 166 of file PmaNode3D.cxx.

References fPoint3D, LimitPoint3D(), and UpdateProj2D().

Referenced by pma::VtxCandidate::JoinTracks(), MakeGradient(), Node3D(), and StepWithGradient().

{
  fPoint3D = p3d;

  bool accepted = !LimitPoint3D();
  UpdateProj2D();

  return accepted;
}

void pma::Node3D::SetProjection ( pma::Hit3D &  h ) const

overridevirtual

Set hit 3D position and its 2D projection to the vertex.

Implements pma::Element3D.

Definition at line 218 of file PmaNode3D.cxx.

References pma::Element3D::fFrozen, fPoint3D, pma::SortedObjectBase::Next(), pma::SortedObjectBase::next, pma::Hit3D::Point2D(), Point3D(), pma::SortedObjectBase::Prev(), pma::SortedObjectBase::prev, Projection2D(), pma::Hit3D::SetPoint3D(), pma::Hit3D::SetProjection(), and pma::Hit3D::View2D().

{
  TVector2 gstart;
  TVector3 g3d;
  if (prev) {
    pma::Node3D* vtx = static_cast<pma::Node3D*>(prev->Prev());
    gstart = vtx->Projection2D(h.View2D());
    if (!next) g3d = vtx->Point3D();
  }
  else if (next) {
    pma::Node3D* vtx = static_cast<pma::Node3D*>(next->Next());
    gstart = Projection2D(h.View2D());
    gstart -= vtx->Projection2D(h.View2D()) - Projection2D(h.View2D());
    if (!prev) {
      g3d = fPoint3D;
      g3d -= vtx->Point3D() - fPoint3D;
    }
  }
  else {
    mf::LogError("pma::Node3D") << "Isolated vertex.";
    TVector2 p(Projection2D(h.View2D()));
    h.SetProjection(p, 0.0F);
    h.SetPoint3D(fPoint3D);
    return;
  }

  TVector2 v0(h.Point2D());
  v0 -= Projection2D(h.View2D());

  TVector2 v1(gstart);
  v1 -= Projection2D(h.View2D());

  double v0Norm = v0.Mod();
  double v1Norm = v1.Mod();
  double mag = v0Norm * v1Norm;
  double cosine = 0.0;
  if (mag != 0.0) cosine = v0 * v1 / mag;

  TVector2 p(Projection2D(h.View2D()));

  if (prev && next) {
    pma::Node3D* vNext = static_cast<pma::Node3D*>(next->Next());
    TVector2 vN(vNext->Projection2D(h.View2D()));
    vN -= Projection2D(h.View2D());

    mag = v0Norm * vN.Mod();
    double cosineN = 0.0;
    if (mag != 0.0) cosineN = v0 * vN / mag;

    // hit on the previous segment side, sorting on the -cosine(prev_seg, point)  /max.val. = 1/
    if (cosineN <= cosine) h.SetProjection(p, -(float)cosine);
    // hit on the next segment side, sorting on the 1+cosine(next_seg, point)  /min.val. = 1/
    else
      h.SetProjection(p, 2.0F + (float)cosineN);

    h.SetPoint3D(fPoint3D);
  }
  else {
    float b = (float)(v0Norm * cosine / v1Norm);
    if (fFrozen) // limit 3D positions to outermose node if frozen
    {
      h.SetPoint3D(fPoint3D);
    }
    else // or set 3D positions along the line of outermost segment
    {
      g3d -= fPoint3D;
      h.SetPoint3D(fPoint3D + (g3d * b));

      p += (v1 * b);
    }
    h.SetProjection(p, -b);
  }
}

void pma::Node3D::SetVertex ( bool  state )

inline

Definition at line 63 of file PmaNode3D.h.

{ fIsVertex = state; }

void pma::Node3D::SetVertexToBranching ( bool  setAllNodes )

inline

Definition at line 64 of file PmaNode3D.h.

  {
    if (setAllNodes || !fIsVertex) fIsVertex = IsBranching();
  }

void pma::Element3D::SortHits ( void  )

inherited

Definition at line 44 of file PmaElement3D.cxx.

References pma::Element3D::fAssignedHits.

Referenced by pma::Track3D::AddNode(), pma::Track3D::SortHits(), and pma::Element3D::UpdateProjection().

{
  std::sort(fAssignedHits.begin(), fAssignedHits.end(), pma::bTrajectory3DOrderLess());
}

double pma::Node3D::StepWithGradient ( float  alfa, float  tol, float  penalty, float  weight  )

private

Definition at line 641 of file PmaNode3D.cxx.

References fGradient, fPoint3D, GetObjFunction(), MakeGradient(), SetPoint3D(), t1, t2, and tmp.

Referenced by Optimize().

{
  unsigned int steps = 0;
  double t, t1, t2, t3, g, g0, g1, g2, g3, p1, p2;
  double eps = 6.0E-37, zero_tol = 1.0E-15;
  TVector3 tmp(fPoint3D), gpoint(fPoint3D);

  g = MakeGradient(penalty, weight);
  if (g < zero_tol) return 0.0;
  g0 = g;

  //**** first three points ****//
  alfa *= 0.8F;
  t2 = 0.0;
  g2 = g;
  t3 = 0.0;
  g3 = g;
  do {
    t1 = t2;
    g1 = g2;
    t2 = t3;
    g2 = g3;

    alfa *= 1.25F;
    t3 += alfa;
    gpoint = tmp;
    gpoint += (fGradient * t3);
    if (!SetPoint3D(gpoint)) // stepped out of allowed volume
    {
      //std::cout << "****  SetPoint trimmed 1 ****" << std::endl;
      g3 = GetObjFunction(penalty, weight);
      if (g3 < g2)
        return (g0 - g3) / g3; // exit with the node at the border
      else {
        SetPoint3D(tmp);
        return 0.0;
      } // exit with restored original position
    }

    g3 = GetObjFunction(penalty, weight);

    if (g3 < zero_tol) return 0.0;

    if (++steps > 1000) {
      SetPoint3D(tmp);
      return 0.0;
    }

  } while (g3 < g2);
  //****************************//

  //*** first step overshoot ***//
  if (steps == 1) {
    t2 = t3;
    g2 = g3;
    do {
      t3 = t2;
      g3 = g2;
      t2 = (t1 * g3 + t3 * g1) / (g1 + g3);

      // small shift...
      t2 = 0.05 * t3 + 0.95 * t2;

      // break: starting point is very close to the minimum
      if (fabs(t2 - t1) < tol) {
        SetPoint3D(tmp);
        return 0.0;
      }

      gpoint = tmp;
      gpoint += (fGradient * t2);
      if (!SetPoint3D(gpoint)) // select the best point to exit
      {
        g2 = GetObjFunction(penalty, weight);
        if (g2 < g0)
          return (g0 - g2) / g2; // exit with the node at the border
        else if (g1 < g0) {
          gpoint = tmp;
          gpoint += (fGradient * t1);
          return (g0 - g1) / g1;
        }
        else if (g3 < g0) {
          gpoint = tmp;
          gpoint += (fGradient * t3);
          return (g0 - g3) / g3;
        }
        else {
          SetPoint3D(tmp);
          return 0.0;
        }
      }
      g2 = GetObjFunction(penalty, weight);

      if (g2 < zero_tol) return 0.0;
      steps++;

    } while (g2 >= g1);
  }
  //****************************//

  while (fabs(t1 - t3) > tol) {
    //*** 3-point minimization ***//
    if ((fabs(t2 - t1) < eps) || (fabs(t2 - t3) < eps)) break; // minimum on the edge
    if ((fabs(g2 - g1) < eps) && (fabs(g2 - g3) < eps)) break; // ~singularity

    p1 = (t2 - t1) * (g2 - g3);
    p2 = (t2 - t3) * (g2 - g1);
    if (fabs(p1 - p2) < eps) break; // ~linearity

    t = t2 + ((t2 - t1) * p1 - (t2 - t3) * p2) / (2 * (p2 - p1));
    if ((t <= t1) || (t >= t3)) t = (t1 * g3 + t3 * g1) / (g1 + g3);

    gpoint = tmp;
    gpoint += (fGradient * t);
    if (!SetPoint3D(gpoint)) // select the best point to exit
    {
      g = GetObjFunction(penalty, weight);
      if ((g < g0) && (g < g1) && (g < g3))
        return (g0 - g) / g; // exit with the node at the border
      else if ((g1 < g0) && (g1 < g3)) {
        gpoint = tmp;
        gpoint += (fGradient * t1);
        return (g0 - g1) / g1;
      }
      else if (g3 < g0) {
        gpoint = tmp;
        gpoint += (fGradient * t3);
        return (g0 - g3) / g3;
      }
      else {
        SetPoint3D(tmp);
        return 0.0;
      }
    }

    g = GetObjFunction(penalty, weight);
    if (g < zero_tol) return 0.0;
    steps++;
    //****************************//

    //*** select next 3 points ***//
    if (fabs(t - t2) < 0.2 * tol) break; // start in a new direction
    if (g < g2) {
      if (t < t2) {
        t3 = t2;
        g3 = g2;
      }
      else {
        t1 = t2;
        g1 = g2;
      }
      t2 = t;
      g2 = g;
    }
    else {
      if (t < t2) {
        t1 = t;
        g1 = g;
      }
      else {
        t3 = t;
        g3 = g;
      }
    }
    //****************************//
  }

  return (g0 - g) / g;
}

double pma::Element3D::SumDist2 ( void  ) const

inherited

Definition at line 164 of file PmaElement3D.cxx.

References d, pma::Element3D::fAssignedHits, pma::Element3D::fAssignedPoints, pma::Element3D::fTPC, pma::Element3D::GetDistance2To(), and pma::Element3D::SumDist2Hits().

Referenced by Mse(), and pma::Element3D::UpdateProjection().

{
  if (fTPC < 0) {
    if (!fAssignedHits.empty())
      mf::LogWarning("pma::Element3D") << "Hits assigned to TPC-crossing element.";
    return 0.0F;
  }

  double hit_sum = SumDist2Hits();

  if (fAssignedPoints.size()) {
    double d, ref_sum = 0.0F;
    for (auto p : fAssignedPoints) {
      d = sqrt(GetDistance2To(*p)) - 0.5; // guide by ref points up to ~ 3D resolution
      if (d > 0.0) ref_sum += d * d;
    }
    if (fAssignedHits.size()) { ref_sum *= 0.2 * fAssignedHits.size() / fAssignedPoints.size(); }
    hit_sum += ref_sum;
  }

  return hit_sum;
}

double pma::Element3D::SumDist2 ( unsigned int  view ) const

inherited

Definition at line 187 of file PmaElement3D.cxx.

References pma::Element3D::fAssignedHits, pma::Element3D::fTPC, pma::Element3D::GetDistance2To(), geo::kUnknown, and pma::Element3D::OptFactor().

{
  if (fTPC < 0) {
    if (!fAssignedHits.empty())
      mf::LogWarning("pma::Element3D") << "Hits assigned to TPC-crossing element.";
    return 0.0F;
  }

  double hit_sum = 0.0F;
  for (auto h : fAssignedHits) {
    if (h->IsEnabled()) {
      unsigned int hitView = h->View2D();
      if ((view == geo::kUnknown) || (view == hitView)) {
        hit_sum += OptFactor(hitView) *                   // alpha_i
                   h->GetSigmaFactor() *                  // hit_amp / hit_max_amp
                   GetDistance2To(h->Point2D(), hitView); // hit_to_fit_dist^2
      }
    }
  }
  return hit_sum;
}

double pma::Node3D::SumDist2Hits ( void  ) const

overrideprivatevirtual

Implements pma::Element3D.

Definition at line 186 of file PmaNode3D.cxx.

References pma::Dist2(), pma::Element3D::fAssignedHits, fProj2D, pma::Element3D::OptFactor(), and sum.

{
  double sum = 0.0F;
  for (auto h : fAssignedHits) {
    if (h->IsEnabled()) {
      unsigned int view = h->View2D();

      sum += OptFactor(view) *     // alpha_i
             h->GetSigmaFactor() * // hit_amp / hit_max_amp
             pma::Dist2(h->Point2D(), fProj2D[view]);
    }
  }
  return sum;
}

double pma::Element3D::SumHitsQ ( unsigned int  view ) const

inlineinherited

Definition at line 98 of file PmaElement3D.h.

References pma::Element3D::fSumHitsQ.

{ return fSumHitsQ[view]; }

int pma::Element3D::TPC ( void  ) const

inlineinherited

TPC index or -1 if out of any TPC.

Definition at line 35 of file PmaElement3D.h.

References pma::Element3D::fTPC.

Referenced by pma::Track3D::AttachBackTo(), pma::Track3D::AttachTo(), pma::Track3D::HasTPC(), IsTPCEdge(), pma::Track3D::MakeFastProjection(), pma::ProjectionMatchingAlg::mergeTracks(), Pi(), pma::Segment3D::Segment3D(), pma::Track3D::Split(), pma::ProjectionMatchingAlg::validate(), pma::ProjectionMatchingAlg::validate_on_adc(), and pma::ProjectionMatchingAlg::validate_on_adc_test().

{ return fTPC; }

void pma::Element3D::UpdateHitParams ( void  )

inherited

Definition at line 56 of file PmaElement3D.cxx.

References pma::Element3D::fAssignedHits, pma::Element3D::fHitsRadius, pma::Element3D::fNHits, pma::Element3D::fNThisHits, pma::Element3D::fNThisHitsEnabledAll, pma::Element3D::fSumHitsQ, pma::GetHitsRadius2D(), geo::kU, geo::kV, geo::kZ, pma::SortedObjectBase::Next(), pma::SortedObjectBase::NextCount(), pma::SortedObjectBase::Prev(), and r.

Referenced by pma::Element3D::AddPoint().

{
  std::vector<pma::Hit3D*> hitsColl, hitsInd1, hitsInd2;
  for (size_t i = 0; i < 3; ++i)
    fNThisHitsEnabledAll = 0;
  for (auto h : fAssignedHits) {
    if (h->IsEnabled()) fNThisHitsEnabledAll++;
    switch (h->View2D()) {
    case geo::kZ: hitsColl.push_back(h); break;
    case geo::kV: hitsInd2.push_back(h); break;
    case geo::kU: hitsInd1.push_back(h); break;
    }
  }
  fNThisHits[0] = hitsInd1.size();
  fNThisHits[1] = hitsInd2.size();
  fNThisHits[2] = hitsColl.size();

  pma::SortedObjectBase const* chain = dynamic_cast<pma::SortedObjectBase*>(this);
  pma::Element3D* el = 0;
  for (size_t b = 0; b < chain->NextCount(); b++) {
    el = dynamic_cast<pma::Element3D*>(chain->Next());
    if (el)
      for (auto h : el->fAssignedHits) {
        switch (h->View2D()) {
        case geo::kZ: hitsColl.push_back(h); break;
        case geo::kV: hitsInd2.push_back(h); break;
        case geo::kU: hitsInd1.push_back(h); break;
        }
      }
  }
  el = dynamic_cast<pma::Element3D*>(chain->Prev());
  if (el) {
    for (auto h : el->fAssignedHits) {
      switch (h->View2D()) {
      case geo::kZ: hitsColl.push_back(h); break;
      case geo::kV: hitsInd2.push_back(h); break;
      case geo::kU: hitsInd1.push_back(h); break;
      }
    }
  }

  fHitsRadius = GetHitsRadius2D(hitsColl);
  double r = GetHitsRadius2D(hitsInd2);
  if (r > fHitsRadius) fHitsRadius = r;
  r = GetHitsRadius2D(hitsInd1);
  if (r > fHitsRadius) fHitsRadius = r;

  float amp, sigmaMax = 0.0F;
  fSumHitsQ[0] = 0.0;
  fNHits[0] = hitsInd1.size();
  for (size_t i = 0; i < hitsInd1.size(); i++) {
    amp = hitsInd1[i]->GetAmplitude();
    if (amp > sigmaMax) sigmaMax = amp;
    fSumHitsQ[0] += amp;
  }
  for (size_t i = 0; i < hitsInd1.size(); i++) {
    if (sigmaMax > 0.0F) {
      amp = hitsInd1[i]->GetAmplitude();
      if (amp > 0.0F)
        hitsInd1[i]->SetSigmaFactor((float)sqrt(amp / sigmaMax));
      else
        hitsInd1[i]->SetSigmaFactor(0.01F);
    }
    else
      hitsInd1[i]->SetSigmaFactor(1.0F);
  }

  sigmaMax = 0.0F;
  fSumHitsQ[1] = 0.0;
  fNHits[1] = hitsInd2.size();
  for (size_t i = 0; i < hitsInd2.size(); i++) {
    amp = hitsInd2[i]->GetAmplitude();
    if (amp > sigmaMax) sigmaMax = amp;
    fSumHitsQ[1] += amp;
  }
  for (size_t i = 0; i < hitsInd2.size(); i++) {
    if (sigmaMax > 0.0F) {
      amp = hitsInd2[i]->GetAmplitude();
      if (amp > 0.0F)
        hitsInd2[i]->SetSigmaFactor((float)sqrt(amp / sigmaMax));
      else
        hitsInd2[i]->SetSigmaFactor(0.01F);
    }
    else
      hitsInd2[i]->SetSigmaFactor(1.0F);
  }

  sigmaMax = 0.0F;
  fSumHitsQ[2] = 0.0;
  fNHits[2] = hitsColl.size();
  for (size_t i = 0; i < hitsColl.size(); i++) {
    amp = hitsColl[i]->SummedADC();
    if (amp > sigmaMax) sigmaMax = amp;
    fSumHitsQ[2] += amp;
  }
  for (size_t i = 0; i < hitsColl.size(); i++) {
    if (sigmaMax > 0.0F) {
      amp = hitsColl[i]->SummedADC();
      if (amp > 0.0F)
        hitsColl[i]->SetSigmaFactor((float)sqrt(amp / sigmaMax));
      else
        hitsColl[i]->SetSigmaFactor(0.01F);
    }
    else
      hitsColl[i]->SetSigmaFactor(1.0F);
  }
}

void pma::Node3D::UpdateProj2D ( )

private

Definition at line 157 of file PmaNode3D.cxx.

References fChannelMap, fDriftOffset, fPoint3D, fProj2D, fTpcGeo, geo::TPCGeo::ID(), geo::Iterable< IterationPolicy, Transform >::Iterate(), and geo::vect::toPoint().

Referenced by SetPoint3D().

{
  unsigned int i = 0;
  for (auto const& plane : fChannelMap.Iterate<geo::PlaneGeo>(fTpcGeo.ID())) {
    fProj2D[i++].Set(plane.PlaneCoordinate(geo::vect::toPoint(fPoint3D)),
                     fPoint3D.X() - fDriftOffset);
  }
}

void pma::Element3D::UpdateProjection ( void  )

inlineinherited

Definition at line 89 of file PmaElement3D.h.

References pma::Element3D::fAssignedHits, pma::Element3D::SetProjection(), pma::Element3D::SortHits(), and pma::Element3D::SumDist2().

  {
    for (auto h : fAssignedHits)
      SetProjection(*h);
  }

Member Data Documentation

std::vector<pma::Hit3D*> pma::Element3D::fAssignedHits

protectedinherited

Definition at line 122 of file PmaElement3D.h.

Referenced by pma::Element3D::AddHit(), pma::Element3D::ClearAssigned(), ClearAssigned(), pma::Element3D::HasHit(), pma::Element3D::Hit(), pma::Element3D::Hits(), pma::Element3D::HitsRadius3D(), pma::Element3D::NEnabledHits(), pma::Element3D::NHits(), pma::Element3D::RemoveHitAt(), pma::Element3D::SelectAllHits(), pma::Element3D::SelectRndHits(), pma::Element3D::SortHits(), pma::Element3D::SumDist2(), pma::Segment3D::SumDist2Hits(), SumDist2Hits(), pma::Element3D::UpdateHitParams(), and pma::Element3D::UpdateProjection().

std::vector<TVector3*> pma::Element3D::fAssignedPoints

protectedinherited

Definition at line 123 of file PmaElement3D.h.

Referenced by pma::Element3D::AddPoint(), pma::Element3D::ClearAssigned(), ClearAssigned(), pma::Element3D::NPoints(), pma::Element3D::ReferencePoint(), and pma::Element3D::SumDist2().

geo::WireReadoutGeom const& pma::Node3D::fChannelMap

private

Definition at line 145 of file PmaNode3D.h.

Referenced by Node3D(), and UpdateProj2D().

int pma::Element3D::fCryo

protectedinherited

Definition at line 117 of file PmaElement3D.h.

Referenced by pma::Element3D::Cryo(), Node3D(), and pma::Segment3D::Segment3D().

double pma::Node3D::fDriftOffset

private

Definition at line 153 of file PmaNode3D.h.

Referenced by UpdateProj2D().

bool pma::Element3D::fFrozen

protectedinherited

Definition at line 121 of file PmaElement3D.h.

Referenced by pma::Element3D::IsFrozen(), Optimize(), pma::Element3D::SetFrozen(), and SetProjection().

bool pma::Node3D::fGradFixed = {false, false, false}

staticprivate

Definition at line 158 of file PmaNode3D.h.

Referenced by MakeGradient().

TVector3 pma::Node3D::fGradient

private

Definition at line 155 of file PmaNode3D.h.

Referenced by MakeGradient(), and StepWithGradient().

double pma::Element3D::fHitsRadius

protectedinherited

Definition at line 128 of file PmaElement3D.h.

Referenced by pma::Element3D::ClearAssigned(), ClearAssigned(), Pi(), PiInWirePlane(), and pma::Element3D::UpdateHitParams().

bool pma::Node3D::fIsVertex

private

Definition at line 156 of file PmaNode3D.h.

Referenced by PenaltyInWirePlane(), and Pi().

double pma::Node3D::fMargin = 3.0

staticprivate

Definition at line 159 of file PmaNode3D.h.

Referenced by LimitPoint3D().

double pma::Node3D::fMaxX

private

Definition at line 147 of file PmaNode3D.h.

Referenced by GetDistToWall(), LimitPoint3D(), Node3D(), and SameTPC().

double pma::Node3D::fMaxY

private

Definition at line 147 of file PmaNode3D.h.

Referenced by GetDistToWall(), LimitPoint3D(), Node3D(), and SameTPC().

double pma::Node3D::fMaxZ

private

Definition at line 147 of file PmaNode3D.h.

Referenced by GetDistToWall(), LimitPoint3D(), Node3D(), and SameTPC().

double pma::Node3D::fMinX

private

Definition at line 147 of file PmaNode3D.h.

Referenced by GetDistToWall(), LimitPoint3D(), Node3D(), and SameTPC().

double pma::Node3D::fMinY

private

Definition at line 147 of file PmaNode3D.h.

Referenced by GetDistToWall(), LimitPoint3D(), Node3D(), and SameTPC().

double pma::Node3D::fMinZ

private

Definition at line 147 of file PmaNode3D.h.

Referenced by GetDistToWall(), LimitPoint3D(), Node3D(), and SameTPC().

size_t pma::Element3D::fNHits[3]

protectedinherited

Definition at line 126 of file PmaElement3D.h.

Referenced by pma::Element3D::Element3D(), pma::Element3D::NHits(), and pma::Element3D::UpdateHitParams().

size_t pma::Element3D::fNThisHits[3]

protectedinherited

Definition at line 124 of file PmaElement3D.h.

Referenced by pma::Element3D::Element3D(), pma::Element3D::NThisHits(), and pma::Element3D::UpdateHitParams().

size_t pma::Element3D::fNThisHitsEnabledAll

protectedinherited

Definition at line 125 of file PmaElement3D.h.

Referenced by pma::Element3D::Element3D(), pma::Element3D::NPrecalcEnabledHits(), and pma::Element3D::UpdateHitParams().

float pma::Element3D::fOptFactors = {0.2F, 0.8F, 1.0F}

staticprotectedinherited

Definition at line 130 of file PmaElement3D.h.

Referenced by pma::Element3D::OptFactor(), and pma::Element3D::SetOptFactor().

TVector3 pma::Node3D::fPoint3D

private

Definition at line 150 of file PmaNode3D.h.

Referenced by pma::Track3D::ApplyDriftShiftInTree(), EndPtCos2Transverse(), GetDistance2To(), GetDistToWall(), LimitPoint3D(), MakeGradient(), SegmentCos(), SegmentCosTransverse(), SegmentCosWirePlane(), SetPoint3D(), SetProjection(), StepWithGradient(), and UpdateProj2D().

TVector2 pma::Node3D::fProj2D[3]

private

Definition at line 151 of file PmaNode3D.h.

Referenced by GetDistance2To(), Node3D(), SumDist2Hits(), and UpdateProj2D().

double pma::Element3D::fSumHitsQ[3]

protectedinherited

Definition at line 127 of file PmaElement3D.h.

Referenced by pma::Element3D::Element3D(), pma::Element3D::SumHitsQ(), and pma::Element3D::UpdateHitParams().

int pma::Element3D::fTPC

protectedinherited

Definition at line 117 of file PmaElement3D.h.

Referenced by pma::Element3D::HitsRadius3D(), Node3D(), pma::Segment3D::Segment3D(), pma::Element3D::SumDist2(), and pma::Element3D::TPC().

geo::TPCGeo const& pma::Node3D::fTpcGeo

private

Definition at line 144 of file PmaNode3D.h.

Referenced by Node3D(), and UpdateProj2D().

pma::SortedObjectBase* pma::SortedObjectBase::next

protectedinherited

Definition at line 53 of file SortedObjects.h.

Referenced by pma::SortedObjectBase::AddNext(), pma::SortedBranchBase::AddNext(), pma::SortedObjectBase::Disconnect(), EndPtCos2Transverse(), pma::Segment3D::GetDirection3D(), GetDirection3D(), pma::Segment3D::GetDistance2To(), pma::Segment3D::GetProjection(), pma::Segment3D::GetUnconstrainedProj3D(), pma::SortedObjectBase::IsLast(), pma::Segment3D::Length2(), Length2(), MakeGradient(), pma::SortedObjectBase::Next(), pma::SortedObjectBase::NextCount(), PenaltyInWirePlane(), Pi(), pma::SortedObjectBase::RemoveNext(), pma::SortedBranchBase::RemoveNext(), SegmentCos(), SegmentCosTransverse(), SegmentCosWirePlane(), pma::Segment3D::SetProjection(), SetProjection(), pma::SortedBranchBase::SortedBranchBase(), and pma::Segment3D::SumDist2Hits().

std::vector<pma::SortedObjectBase*> pma::SortedBranchBase::next_vector

protectedinherited

Definition at line 88 of file SortedObjects.h.

pma::SortedObjectBase* pma::SortedObjectBase::prev

protectedinherited

Definition at line 54 of file SortedObjects.h.

Referenced by pma::SortedObjectBase::AddNext(), pma::SortedBranchBase::AddNext(), pma::SortedObjectBase::Disconnect(), pma::SortedBranchBase::Disconnect(), EndPtCos2Transverse(), pma::Segment3D::GetDirection3D(), GetDirection3D(), pma::Segment3D::GetDistance2To(), pma::Segment3D::GetProjection(), pma::Segment3D::GetUnconstrainedProj3D(), IsBranching(), pma::SortedObjectBase::IsFirst(), IsTPCEdge(), pma::Segment3D::Length2(), Length2(), MakeGradient(), Mse(), Penalty(), PenaltyInWirePlane(), Pi(), PiInWirePlane(), pma::SortedObjectBase::Prev(), pma::SortedObjectBase::RemoveNext(), pma::SortedBranchBase::RemoveNext(), SegmentCos(), SegmentCosTransverse(), SegmentCosWirePlane(), pma::Segment3D::SetProjection(), SetProjection(), and pma::Segment3D::SumDist2Hits().

---

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

• larreco/v10_01_07/source/larreco/RecoAlg/PMAlg/PmaNode3D.h
• larreco/v10_01_07/source/larreco/RecoAlg/PMAlg/PmaNode3D.cxx