pma Namespace Reference

Classes
class	bDistCenterLess2D
class	bDistCenterLess3D
class	bSegmentProjLess
struct	bTrack3DLonger
struct	bTrajectory3DDistLess
struct	bTrajectory3DOrderLess
class	Element3D
class	Hit3D
class	Node3D
class	PMAlgCosmicTagger
class	PMAlgFitter
class	PMAlgStitching
class	PMAlgTracker
class	PMAlgTrackingBase
class	PMAlgVertexing
class	ProjectionMatchingAlg
class	Segment3D
class	SortedBranchBase
class	SortedObjectBase
class	Track3D
class	TrkCandidate
class	TrkCandidateColl
class	VtxCandidate

Typedefs
typedef std::map< unsigned int, std::vector< art::Ptr< recob::Hit > > >	view_hitmap
typedef std::map< unsigned int, view_hitmap >	tpc_view_hitmap
typedef std::map< unsigned int, tpc_view_hitmap >	cryo_tpc_view_hitmap
typedef ROOT::Math::DisplacementVector2D< ROOT::Math::Cartesian2D< double > >	Vector2D
typedef recob::tracking::Vector_t	Vector3D
typedef std::map< size_t, std::vector< double > >	dedx_map
typedef std::map< size_t, pma::TrkCandidateColl >	tpc_track_map

Functions
double	Dist2 (const TVector2 &v1, const TVector2 &v2)
double	Dist2 (const Vector2D &v1, const Vector2D &v2)
double	Dist2 (const TVector3 &v1, const TVector3 &v2)
double	Dist2 (const Vector3D &v1, const Vector3D &v2)
size_t	GetHitsCount (const std::vector< pma::Hit3D * > &hits, unsigned int view)
double	GetSummedADC (const std::vector< pma::Hit3D * > &hits, unsigned int view=geo::kUnknown)
double	GetSummedAmpl (const std::vector< pma::Hit3D * > &hits, unsigned int view=geo::kUnknown)
double	GetHitsRadius3D (const std::vector< pma::Hit3D * > &hits, bool exact=false)
double	GetHitsRadius2D (const std::vector< pma::Hit3D * > &hits, bool exact=false)
double	GetSegmentProjVector (const TVector2 &p, const TVector2 &p0, const TVector2 &p1)
double	GetSegmentProjVector (const Vector2D &p, const Vector2D &p0, const Vector2D &p1)
double	GetSegmentProjVector (const TVector3 &p, const TVector3 &p0, const TVector3 &p1)
double	GetSegmentProjVector (const Vector3D &p, const Vector3D &p0, const Vector3D &p1)
TVector2	GetProjectionToSegment (const TVector2 &p, const TVector2 &p0, const TVector2 &p1)
TVector3	GetProjectionToSegment (const TVector3 &p, const TVector3 &p0, const TVector3 &p1)
double	SolveLeastSquares3D (const std::vector< std::pair< TVector3, TVector3 > > &lines, TVector3 &result)
TVector2	GetProjectionToPlane (const TVector3 &p, unsigned int plane, unsigned int tpc, unsigned int cryo)
TVector2	GetVectorProjectionToPlane (const TVector3 &v, unsigned int plane, unsigned int tpc, unsigned int cryo)
TVector2	WireDriftToCm (unsigned int wire, float drift, unsigned int plane, unsigned int tpc, unsigned int cryo)
TVector2	CmToWireDrift (float xw, float yd, unsigned int plane, unsigned int tpc, unsigned int cryo)
recob::Track	convertFrom (const pma::Track3D &src, unsigned int tidx, int pdg=0)

Typedef Documentation

typedef std::map< unsigned int, tpc_view_hitmap > pma::cryo_tpc_view_hitmap

Definition at line 22 of file PmaTrkCandidate.h.

typedef std::map< size_t, std::vector<double> > pma::dedx_map

Definition at line 31 of file Utilities.h.

typedef std::map< size_t, pma::TrkCandidateColl > pma::tpc_track_map

Definition at line 46 of file PMAlgTracking.h.

typedef std::map< unsigned int, view_hitmap > pma::tpc_view_hitmap

Definition at line 21 of file PmaTrkCandidate.h.

typedef ROOT::Math::DisplacementVector2D< ROOT::Math::Cartesian2D<double> > pma::Vector2D

Definition at line 28 of file Utilities.h.

typedef recob::tracking::Vector_t pma::Vector3D

Definition at line 29 of file Utilities.h.

typedef std::map< unsigned int, std::vector< art::Ptr<recob::Hit> > > pma::view_hitmap

Definition at line 20 of file PmaTrkCandidate.h.

Function Documentation

TVector2 pma::CmToWireDrift	(	float	xw,
		float	yd,
		unsigned int	plane,
		unsigned int	tpc,
		unsigned int	cryo
	)

Definition at line 318 of file Utilities.cxx.

References detinfo::DetectorProperties::ConvertXToTicks(), geo::TPCGeo::Plane(), geo::GeometryCore::TPC(), and geo::PlaneGeo::WirePitch().

Referenced by pma::Track3D::DriftsOfWireIntersection(), and pma::Track3D::WireDriftRange().

{
        art::ServiceHandle<geo::Geometry> geom;
        const detinfo::DetectorProperties* detprop = lar::providerFrom<detinfo::DetectorPropertiesService>();

        return TVector2(
                xw / geom->TPC(tpc, cryo).Plane(plane).WirePitch(),
                detprop->ConvertXToTicks(yd, plane, tpc, cryo)
        );
}

recob::Track pma::convertFrom	(	const pma::Track3D &	src,
		unsigned int	tidx,
		int	pdg = 0
	)

Definition at line 21 of file PMAlgTracking.cxx.

References pma::Track3D::GetDirection3D(), recob::TrajectoryPointFlags::makeMask(), and pma::Track3D::size().

Referenced by trkf::PMAlgTrajFitter::produce(), and trkf::PMAlgTrackMaker::produce().

{
    std::vector<Point_t> positions;                    positions.reserve(src.size());
    std::vector<Vector_t> momenta;                     momenta.reserve(src.size());
    std::vector<recob::TrajectoryPointFlags> outFlags; outFlags.reserve(src.size());

        for (size_t i = 0, h = 0; i < src.size(); i++)
                if (src[i]->IsEnabled())
        {
        auto const & point3d = src[i]->Point3D();
        positions.emplace_back(point3d.X(), point3d.Y(), point3d.Z());
        momenta.push_back(src.GetDirection3D(i));
        outFlags.emplace_back(h++, recob::TrajectoryPointFlags::makeMask());
        }

    int ndof = 0;
    float totChi2 = 0;

    SMatrixSym55 covStart, covEnd;
    return recob::Track(
        recob::TrackTrajectory(std::move(positions), std::move(momenta), std::move(outFlags), false),
        pdg, totChi2, ndof, covStart, covEnd, tidx);
}

double pma::Dist2	(	const TVector2 &	v1,
		const TVector2 &	v2
	)

Definition at line 19 of file Utilities.cxx.

Referenced by pma::ProjectionMatchingAlg::alignTracks(), ems::MultiEMShowers::analyze(), pma::PMAlgTracker::areCoLinear(), pma::Track3D::AttachTo(), pma::ProjectionMatchingAlg::buildMultiTPCTrack(), tss::Segmentation2D::buildSegment(), pma::ProjectionMatchingAlg::buildSegment(), tss::Cluster2D::closest(), ems::MergeEMShower3D::collectshowers(), pma::PMAlgTracker::collectSingleViewEnd(), pma::PMAlgTracker::collectSingleViewFront(), ems::DirOfGamma::ComputeMaxDist(), tss::Cluster2D::dist2(), pma::ProjectionMatchingAlg::FilterOutSmallParts(), ems::DirOfGamma::FindCandidates(), ems::DirOfGamma::FindInitialPart(), ems::MergeEMShower3D::getBestPoint(), ems::EMShower3D::GetCloseHits(), pma::ProjectionMatchingAlg::GetCloseHits(), ems::MergeEMShower3D::getCos3D(), pma::Hit3D::GetDist2ToProj(), pma::Node3D::GetDistance2To(), GetHitsRadius2D(), GetHitsRadius3D(), ems::Bin2D::GetIniHits(), ems::MultiEMShowers::getMinDist(), pma::Track3D::HitDxByView(), pma::Element3D::HitsRadius3D(), pma::Track3D::InitFromHits(), pma::VtxCandidate::JoinTracks(), pma::Track3D::Length(), tss::Cluster2D::length2(), pma::Segment3D::Length2(), pma::PMAlgVertexing::makeVertices(), pma::PMAlgTracker::mergeCoLinear(), tss::Segmentation2D::mergeDenseParts(), pma::ProjectionMatchingAlg::mergeTracks(), pma::VtxCandidate::MergeWith(), tss::bDistToPointLess::operator()(), pma::bTrack3DLonger::operator()(), pma::bDistCenterLess2D::operator()(), pma::bDistCenterLess3D::operator()(), tss::Cluster2D::outermost(), ems::MergeEMShower3D::produce(), ems::EMShower3D::Reoptimize(), tss::Segmentation2D::run(), tss::Segmentation2D::selectRing(), ems::ShowerInfo::SetP0Dist(), pma::ProjectionMatchingAlg::ShortenSeg(), tss::Segmentation2D::splitHitsNaive(), pma::Node3D::SumDist2Hits(), tss::Segmentation2D::tagDenseEnds(), pma::ProjectionMatchingAlg::TestTrk(), ems::EMShower3D::Validate(), pma::ProjectionMatchingAlg::validate(), and pma::ProjectionMatchingAlg::validate_on_adc_test().

{
        double dx = v1.X() - v2.X(), dy = v1.Y() - v2.Y();
        return dx * dx + dy * dy;
}

double pma::Dist2	(	const Vector2D &	v1,
		const Vector2D &	v2
	)

Definition at line 24 of file Utilities.cxx.

{
        double dx = v1.X() - v2.X(), dy = v1.Y() - v2.Y();
        return dx * dx + dy * dy;
}

double pma::Dist2	(	const TVector3 &	v1,
		const TVector3 &	v2
	)

Definition at line 30 of file Utilities.cxx.

{
        double dx = v1.X() - v2.X(), dy = v1.Y() - v2.Y(), dz = v1.Z() - v2.Z();
        return dx * dx + dy * dy + dz * dz;
}

double pma::Dist2	(	const Vector3D &	v1,
		const Vector3D &	v2
	)

Definition at line 35 of file Utilities.cxx.

{
        double dx = v1.X() - v2.X(), dy = v1.Y() - v2.Y(), dz = v1.Z() - v2.Z();
        return dx * dx + dy * dy + dz * dz;
}

size_t pma::GetHitsCount	(	const std::vector< pma::Hit3D * > &	hits,
		unsigned int	view
	)

Definition at line 41 of file Utilities.cxx.

References geo::kUnknown, and n.

{
        size_t n = 0;
        for (auto const& hit : hits)
                if ((view == geo::kUnknown) || (view == hit->View2D())) n++;
        return n;
}

double pma::GetHitsRadius2D	(	const std::vector< pma::Hit3D * > &	hits,
		bool	exact = false
	)

Definition at line 89 of file Utilities.cxx.

References Dist2(), and pmtana::mean().

Referenced by pma::Element3D::UpdateHitParams(), and pma::Track3D::UpdateHitsRadius().

{
        if (!exact && (hits.size() < 5)) return 0.0;

        if (hits.size() == 0) return 0.0;

        TVector2 mean(0, 0);
        for (size_t i = 0; i < hits.size(); i++)
        {
                mean += hits[i]->Point2D();
        }
        mean *= (1.0 / hits.size());

        double r2, max_r2 = pma::Dist2(hits.front()->Point2D(), mean);
        for (size_t i = 1; i < hits.size(); i++)
        {
                r2 = pma::Dist2(hits[i]->Point2D(), mean);
                if (r2 > max_r2) max_r2 = r2;
        }
        return sqrt(max_r2);
}

double pma::GetHitsRadius3D	(	const std::vector< pma::Hit3D * > &	hits,
		bool	exact = false
	)

Definition at line 67 of file Utilities.cxx.

References Dist2(), and pmtana::mean().

{
        if (!exact && (hits.size() < 5)) return 0.0;

        if (hits.size() == 0) return 0.0;

        TVector3 mean(0, 0, 0);
        for (size_t i = 0; i < hits.size(); i++)
        {
                mean += hits[i]->Point3D();
        }
        mean *= (1.0 / hits.size());

        double r2, max_r2 = pma::Dist2(hits.front()->Point3D(), mean);
        for (size_t i = 1; i < hits.size(); i++)
        {
                r2 = pma::Dist2(hits[i]->Point3D(), mean);
                if (r2 > max_r2) max_r2 = r2;
        }
        return sqrt(max_r2);
}

TVector2 pma::GetProjectionToPlane	(	const TVector3 &	p,
		unsigned int	plane,
		unsigned int	tpc,
		unsigned int	cryo
	)

Definition at line 291 of file Utilities.cxx.

References geo::TPCGeo::Plane(), geo::PlaneGeo::PlaneCoordinate(), and geo::GeometryCore::TPC().

Referenced by pma::ProjectionMatchingAlg::buildShowerSeg(), pma::VtxCandidate::ComputeMse2D(), ems::MultiEMShowers::getMinDist(), GetVectorProjectionToPlane(), and ems::EMShower3D::Validate().

{
        art::ServiceHandle<geo::Geometry> geom;

    return TVector2(geom->TPC(tpc, cryo).Plane(plane).PlaneCoordinate(p), p.X());
}

TVector2 pma::GetProjectionToSegment	(	const TVector2 &	p,
		const TVector2 &	p0,
		const TVector2 &	p1
	)

Definition at line 167 of file Utilities.cxx.

References GetSegmentProjVector().

Referenced by pma::PMAlgTracker::areCoLinear(), tss::Segmentation2D::buildSegment(), ems::MergeEMShower3D::collectshowers(), pma::VtxCandidate::Compute(), pma::VtxCandidate::JoinTracks(), ems::ShowerInfo::Pointsto(), and SolveLeastSquares3D().

{
        TVector2 v1(p1); v1 -= p0;

        double b = GetSegmentProjVector(p, p0, p1);
        TVector2 r(p0);
        r += (v1 * b);
        return r;
}

TVector3 pma::GetProjectionToSegment	(	const TVector3 &	p,
		const TVector3 &	p0,
		const TVector3 &	p1
	)

Definition at line 177 of file Utilities.cxx.

References GetSegmentProjVector().

{
        TVector3 v1(p1); v1 -= p0;

        double b = GetSegmentProjVector(p, p0, p1);
        TVector3 r(p0);
        r += (v1 * b);
        return r;
}

double pma::GetSegmentProjVector	(	const TVector2 &	p,
		const TVector2 &	p0,
		const TVector2 &	p1
	)

Definition at line 111 of file Utilities.cxx.

Referenced by GetProjectionToSegment(), pma::VtxCandidate::JoinTracks(), pma::bSegmentProjLess::operator()(), pma::ProjectionMatchingAlg::validate(), pma::ProjectionMatchingAlg::validate_on_adc(), and pma::ProjectionMatchingAlg::validate_on_adc_test().

{
        TVector2 v0(p); v0 -= p0;
        TVector2 v1(p1); v1 -= p0;

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

        return v0Norm * cosine / v1Norm;
}

double pma::GetSegmentProjVector	(	const Vector2D &	p,
		const Vector2D &	p0,
		const Vector2D &	p1
	)

Definition at line 125 of file Utilities.cxx.

{
        pma::Vector2D v0(p); v0 -= p0;
        pma::Vector2D v1(p1); v1 -= p0;

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

        return v0Norm * cosine / v1Norm;
}

double pma::GetSegmentProjVector	(	const TVector3 &	p,
		const TVector3 &	p0,
		const TVector3 &	p1
	)

Definition at line 139 of file Utilities.cxx.

{
        TVector3 v0(p); v0 -= p0;
        TVector3 v1(p1); v1 -= p0;

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

        return v0Norm * cosine / v1Norm;
}

double pma::GetSegmentProjVector	(	const Vector3D &	p,
		const Vector3D &	p0,
		const Vector3D &	p1
	)

Definition at line 153 of file Utilities.cxx.

{
        pma::Vector3D v0(p); v0 -= p0;
        pma::Vector3D v1(p1); v1 -= p0;

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

        return v0Norm * cosine / v1Norm;
}

double pma::GetSummedADC	(	const std::vector< pma::Hit3D * > &	hits,
		unsigned int	view = geo::kUnknown
	)

Definition at line 49 of file Utilities.cxx.

References geo::kUnknown.

{
        double sum = 0.0;
        for (auto const& hit : hits)
                if ((view == geo::kUnknown) || (view == hit->View2D()))
                        sum += hit->SummedADC();
        return sum;
}

double pma::GetSummedAmpl	(	const std::vector< pma::Hit3D * > &	hits,
		unsigned int	view = geo::kUnknown
	)

Definition at line 58 of file Utilities.cxx.

References geo::kUnknown.

{
        double sum = 0.0;
        for (auto const& hit : hits)
                if ((view == geo::kUnknown) || (view == hit->View2D()))
                        sum += hit->GetAmplitude();
        return sum;
}

TVector2 pma::GetVectorProjectionToPlane	(	const TVector3 &	v,
		unsigned int	plane,
		unsigned int	tpc,
		unsigned int	cryo
	)

Definition at line 298 of file Utilities.cxx.

References GetProjectionToPlane().

{
        TVector3 v0_3d(0., 0., 0.);
        TVector2 v0_2d = GetProjectionToPlane(v0_3d, plane, tpc, cryo);
        TVector2 v1_2d = GetProjectionToPlane(v, plane, tpc, cryo);

        return v1_2d - v0_2d;
}

double pma::SolveLeastSquares3D	(	const std::vector< std::pair< TVector3, TVector3 > > &	lines,
		TVector3 &	result
	)

Definition at line 187 of file Utilities.cxx.

References dir, GetProjectionToSegment(), w, x, and y.

Referenced by ems::MultiEMShowers::analyze(), ems::MergeEMShower3D::collectshowers(), pma::VtxCandidate::Compute(), and ems::ShowerInfo::Pointsto().

{
        // RS: please, ask me if you need examples/explanation of formulas as they
        // are not easy to derive from the code solely; I have Mathcad sources that
        // were used to test the solving method, weighting, etc.

        result.SetXYZ(0., 0., 0.);
        if (lines.size() < 2)
        {
                mf::LogError("pma::SolveLeastSquares3D") << "Need min. two lines.";
                return -1.0;
        }

        double m;
        std::vector< TVectorT<double> > U, P;
        for (size_t v = 0; v < lines.size(); v++)
        {
                TVector3 point = lines[v].first;
                TVector3 dir = lines[v].second;
                dir -= point; m = dir.Mag();
                if (m > 0.0)
                {
                        dir *= 1.0 / m;

                        P.push_back(TVectorT<double>(3));
                        P.back()[0] = point.X(); P.back()[1] = point.Y(); P.back()[2] = point.Z();

                        U.push_back(TVectorT<double>(3));
                        U.back()[0] = dir.X(); U.back()[1] = dir.Y(); U.back()[2] = dir.Z();
                }
                else mf::LogWarning("pma::SolveLeastSquares3D") << "Line undefined.";
        }
        if (P.size() < 2)
        {
                mf::LogError("pma::SolveLeastSquares3D") << "Need min. two lines.";
                return -1.0;
        }

        TVectorT<double> x(3), y(3), w(3);
        TMatrixT<double> A(3, 3);
        double ur, uc, pc;
        double s_uc2[3], s_ur_uc[3];
        double s_p_uc2[3], s_p_ur_uc[3];

        w[0] = 1.0; w[1] = 1.0; w[2] = 1.0;
        for (size_t r = 0; r < 3; r++)
        {
                y[r] = 0.0;
                for (size_t c = 0; c < 3; c++)
                {
                        s_uc2[c] = 0.0; s_ur_uc[c] = 0.0;
                        s_p_uc2[c] = 0.0; s_p_ur_uc[c] = 0.0;

                        for (size_t v = 0; v < P.size(); v++)
                        {
                                //w[1] = fWeights[v]; // to remember that individual coordinates can be supressed...
                                //w[2] = fWeights[v];

                                ur = U[v][r];
                                uc = U[v][c];
                                pc = P[v][c];

                                s_uc2[c] += w[r] * w[c] * (1 - uc * uc);
                                s_p_uc2[c] += w[r] * w[r] * pc * (1 - uc * uc);

                                s_ur_uc[c] += w[r] * w[c] * ur * uc;
                                s_p_ur_uc[c] += w[r] * w[r] * pc * ur * uc;
                        }

                        if (r == c)
                        {
                                y[r] += s_p_uc2[c];
                                A(r, c) = s_uc2[c];
                        }
                        else
                        {
                                y[r] -= s_p_ur_uc[c];
                                A(r, c) = -s_ur_uc[c];
                        }
                }
        }
        try { x = A.InvertFast() * y; }
        catch (...)
        {
                result.SetXYZ(0., 0., 0.); return 1.0e12;
        }

        result.SetXYZ(x[0], x[1], x[2]);

        TVector3 pproj;
        double dx, dy, dz, mse = 0.0;
        for (size_t v = 0; v < lines.size(); v++)
        {
                pproj = pma::GetProjectionToSegment(result, lines[v].first, lines[v].second);

                dx = result.X() - pproj.X(); // dx, dy, dz and the result point can be weighted
                dy = result.Y() - pproj.Y(); // here (linearly) by each line uncertainty
                dz = result.Z() - pproj.Z();

                mse += dx * dx + dy * dy + dz * dz;
        }
        return mse / lines.size();
}

TVector2 pma::WireDriftToCm	(	unsigned int	wire,
		float	drift,
		unsigned int	plane,
		unsigned int	tpc,
		unsigned int	cryo
	)

Definition at line 307 of file Utilities.cxx.

References detinfo::DetectorProperties::ConvertTicksToX(), geo::TPCGeo::Plane(), geo::GeometryCore::TPC(), and geo::PlaneGeo::WirePitch().

Referenced by pma::ProjectionMatchingAlg::FilterOutSmallParts(), ems::EMShower3D::GetCloseHits(), ems::MultiEMShowers::getMinDist(), pma::Track3D::GetUnconstrainedProj3D(), tss::Hit2D::Hit2D(), pma::Hit3D::Hit3D(), pma::Track3D::TestHitsMse(), and pma::ProjectionMatchingAlg::validate().

{
        art::ServiceHandle<geo::Geometry> geom;
        const detinfo::DetectorProperties* detprop = lar::providerFrom<detinfo::DetectorPropertiesService>();

        return TVector2(
                geom->TPC(tpc, cryo).Plane(plane).WirePitch() * wire,
                detprop->ConvertTicksToX(drift, plane, tpc, cryo)
        );
}