dunefd::IniSegAlg Class Reference

#include <IniSegAlg.h>

Public Member Functions
	IniSegAlg (std::map< size_t, std::vector< dunefd::Hit2D > > clusters)
	IniSegAlg (std::vector< art::Ptr< recob::Track > > const &tracks, TVector3 const &mcvtx)
void	FeedwithMc (TVector2 const &vtx, TVector2 const &dir, TVector3 const &dir3d)
void	FeedwithMc (TVector3 const &dir3d)
std::map< size_t, std::vector< dunefd::Hit2D > > const &	GetSelectedCl () const
std::vector< dunefd::Hit2D > const &	GetCl () const
const bool	IsFound ()
art::Ptr< recob::Track > const &	GetTrk () const
float const &	GetDist () const
float const &	GetCos () const

Private Member Functions
void	FindClustersInRad ()
void	FindCluster ()
void	SortLess ()
void	Find3dTrack ()
TVector2	ClusterDir (std::vector< Hit2D > const &hits)

Private Attributes
std::map< size_t, std::vector< dunefd::Hit2D > >	fClusters
std::map< size_t, std::vector< dunefd::Hit2D > >	fSelCls
std::vector< art::Ptr< recob::Track > >	fSelTrks
art::Ptr< recob::Track >	fTrk
bool	fFound
TVector3	fDir3d
TVector3	fFront
TVector3	fBack
std::vector< dunefd::Hit2D >	fCl
TVector2	fMcVtx
TVector3	fMcVtx3d
TVector2	fDir
float	fRadius
float	fDistVtxCl
float const	fThrcos
float	fCos

Detailed Description

Definition at line 43 of file IniSegAlg.h.

Constructor & Destructor Documentation

dunefd::IniSegAlg::IniSegAlg

(

std::map< size_t, std::vector< dunefd::Hit2D > >

clusters

)

Definition at line 27 of file IniSegAlg.cxx.

                                                                           :
fClusters(clusters),
fRadius(10.0),
fDistVtxCl(0.0F),
fThrcos(0.9),
fCos(0.0F)
{
}

dunefd::IniSegAlg::IniSegAlg	(	std::vector< art::Ptr< recob::Track > > const &	tracks,
		TVector3 const &	mcvtx
	)

Definition at line 36 of file IniSegAlg.cxx.

                                                                                                 :
fRadius(10.0),
fDistVtxCl(0.0F),
fThrcos(0.9),
fCos(-9999.0F)
{
        fSelTrks = tracks;
        fMcVtx3d = mcvtx;
        fFound = false;
}

Member Function Documentation

TVector2 dunefd::IniSegAlg::ClusterDir ( std::vector< Hit2D > const &  hits )

private

Definition at line 80 of file IniSegAlg.cxx.

{
        TVector2 p1 = hits[0].GetPointCm();
        size_t id = 0;
        for (size_t i = 0; i < hits.size(); ++i)
                if (pma::Dist2(hits[0].GetPointCm(), hits[i].GetPointCm()) < 1.0)
                {
                        id = i; 
                }
                else break;

        TVector2 p2 = hits[id].GetPointCm(); 
        TVector2 dir = (p2 - p1) * (1 / (p2 - p1).Mod());

        return dir;     
}

void dunefd::IniSegAlg::FeedwithMc	(	TVector2 const &	vtx,
		TVector2 const &	dir,
		TVector3 const &	dir3d
	)

Definition at line 47 of file IniSegAlg.cxx.

{
        fMcVtx = vtx;
        fDir = dir;
        fDir3d = dir3d;
        FindClustersInRad();
        SortLess();
        FindCluster();
}

void dunefd::IniSegAlg::FeedwithMc

(

TVector3 const &

dir3d

)

Definition at line 57 of file IniSegAlg.cxx.

{
        fDir3d = dir3d;
        Find3dTrack();
}

void dunefd::IniSegAlg::Find3dTrack ( )

private

Definition at line 120 of file IniSegAlg.cxx.

{
        TVector3 larStart;
        TVector3 larEnd;

        double maxcos = 0.0; 
        const double thrdist = 10; // cm

        for (auto& trk: fSelTrks)
        {
                std::tie(larStart,larEnd) = trk->Direction<TVector3>(); // correct, check both directions
                TVector3 dir(larStart[0], larStart[1], larStart[2]);
        
                if (!trk->NumberTrajectoryPoints()) continue;
                TVector3 pos_p = trk->LocationAtPoint<TVector3>(0);
                TVector3 pos_end = trk->LocationAtPoint<TVector3>(trk->NumberTrajectoryPoints()-1);

                double dist = std::sqrt(pma::Dist2(pos_p, fMcVtx3d));           

                if (dist > std::sqrt(pma::Dist2(pos_end, fMcVtx3d))) continue;

                double cos = fDir3d * dir;
                
                if ((cos > maxcos) && (dist < thrdist))
                {               
                        maxcos = cos;
                        fTrk = trk;
                        fCos = cos;
                        fFound = true;
                }
        }
}

void dunefd::IniSegAlg::FindCluster ( )

private

Definition at line 97 of file IniSegAlg.cxx.

{
        double maxcos = 0;
        std::map<size_t, std::vector<dunefd::Hit2D> > chosen;
        for (auto& cl: fSelCls)
                if (cl.second.size() > 2)
                {
                        TVector2 dir = ClusterDir(cl.second);
                        double cos = fDir * dir;
                        if (cos > maxcos) 
                        {
                                maxcos = cos;
                                chosen[cl.first] = cl.second;
                                fCl.clear();
                                fCl = cl.second;
                                fDistVtxCl = std::sqrt(pma::Dist2(fCl[0].GetPointCm(), fMcVtx));
                        }
                }

        fSelCls.clear();
        fSelCls = chosen;
}

void dunefd::IniSegAlg::FindClustersInRad ( )

private

Definition at line 63 of file IniSegAlg.cxx.

{
        for (auto& cl: fClusters)
                for (size_t h = 0; h < cl.second.size(); ++h)
                        if (pma::Dist2(fMcVtx, cl.second[h].GetPointCm()) < fRadius*fRadius)
                        {
                                fSelCls[cl.first] = cl.second;
                                break;
                        }
}

std::vector< dunefd::Hit2D > const& dunefd::IniSegAlg::GetCl ( ) const

inline

Definition at line 52 of file IniSegAlg.h.

{ return fCl; }

float const& dunefd::IniSegAlg::GetCos ( ) const

inline

Definition at line 58 of file IniSegAlg.h.

{ return fCos; }

float const& dunefd::IniSegAlg::GetDist ( ) const

inline

Definition at line 57 of file IniSegAlg.h.

{ return fDistVtxCl; } 

std::map<size_t, std::vector< dunefd::Hit2D > > const& dunefd::IniSegAlg::GetSelectedCl ( ) const

inline

Definition at line 51 of file IniSegAlg.h.

{ return fSelCls; } 

art::Ptr<recob::Track> const& dunefd::IniSegAlg::GetTrk ( ) const

inline

Definition at line 55 of file IniSegAlg.h.

{ return fTrk; }

const bool dunefd::IniSegAlg::IsFound ( )

inline

Definition at line 54 of file IniSegAlg.h.

{ return fFound; }

void dunefd::IniSegAlg::SortLess ( )

private

Definition at line 74 of file IniSegAlg.cxx.

{
        for (auto& cl: fSelCls)
                std::sort(cl.second.begin(), cl.second.end(), bDistCentLess2D(fMcVtx));
}

Member Data Documentation

TVector3 dunefd::IniSegAlg::fBack

private

Definition at line 79 of file IniSegAlg.h.

std::vector< dunefd::Hit2D > dunefd::IniSegAlg::fCl

private

Definition at line 82 of file IniSegAlg.h.

std::map<size_t, std::vector<dunefd::Hit2D> > dunefd::IniSegAlg::fClusters

private

Definition at line 69 of file IniSegAlg.h.

float dunefd::IniSegAlg::fCos

private

Definition at line 90 of file IniSegAlg.h.

TVector2 dunefd::IniSegAlg::fDir

private

Definition at line 85 of file IniSegAlg.h.

TVector3 dunefd::IniSegAlg::fDir3d

private

Definition at line 76 of file IniSegAlg.h.

float dunefd::IniSegAlg::fDistVtxCl

private

Definition at line 88 of file IniSegAlg.h.

bool dunefd::IniSegAlg::fFound

private

Definition at line 75 of file IniSegAlg.h.

TVector3 dunefd::IniSegAlg::fFront

private

Definition at line 78 of file IniSegAlg.h.

TVector2 dunefd::IniSegAlg::fMcVtx

private

Definition at line 83 of file IniSegAlg.h.

TVector3 dunefd::IniSegAlg::fMcVtx3d

private

Definition at line 84 of file IniSegAlg.h.

float dunefd::IniSegAlg::fRadius

private

Definition at line 87 of file IniSegAlg.h.

std::map<size_t, std::vector<dunefd::Hit2D> > dunefd::IniSegAlg::fSelCls

private

Definition at line 70 of file IniSegAlg.h.

std::vector< art::Ptr<recob::Track> > dunefd::IniSegAlg::fSelTrks

private

Definition at line 73 of file IniSegAlg.h.

float const dunefd::IniSegAlg::fThrcos

private

Definition at line 89 of file IniSegAlg.h.

art::Ptr<recob::Track> dunefd::IniSegAlg::fTrk

private

Definition at line 74 of file IniSegAlg.h.

---

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

• dunereco/dunereco/FDSensOpt/IniSegAlg/IniSegAlg.h
• dunereco/dunereco/FDSensOpt/IniSegAlg/IniSegAlg.cxx