dune::EmLikeHits Class Reference

Inheritance diagram for dune::EmLikeHits:

Public Member Functions
	EmLikeHits (fhicl::ParameterSet const &p)
	EmLikeHits (EmLikeHits const &)=delete
	EmLikeHits (EmLikeHits &&)=delete
EmLikeHits &	operator= (EmLikeHits const &)=delete
EmLikeHits &	operator= (EmLikeHits &&)=delete
void	reconfigure (fhicl::ParameterSet const &p)
void	produce (art::Event &e) override
Public Member Functions inherited from art::EDProducer
	EDProducer (fhicl::ParameterSet const &pset)
template<typename Config >
	EDProducer (Table< Config > const &config)
std::string	workerType () const
Public Member Functions inherited from art::detail::Producer
virtual	~Producer () noexcept
	Producer (fhicl::ParameterSet const &)
	Producer (Producer const &)=delete
	Producer (Producer &&)=delete
Producer &	operator= (Producer const &)=delete
Producer &	operator= (Producer &&)=delete
void	doBeginJob (SharedResources const &resources)
void	doEndJob ()
void	doRespondToOpenInputFile (FileBlock const &fb)
void	doRespondToCloseInputFile (FileBlock const &fb)
void	doRespondToOpenOutputFiles (FileBlock const &fb)
void	doRespondToCloseOutputFiles (FileBlock const &fb)
bool	doBeginRun (RunPrincipal &rp, ModuleContext const &mc)
bool	doEndRun (RunPrincipal &rp, ModuleContext const &mc)
bool	doBeginSubRun (SubRunPrincipal &srp, ModuleContext const &mc)
bool	doEndSubRun (SubRunPrincipal &srp, ModuleContext const &mc)
bool	doEvent (EventPrincipal &ep, ModuleContext const &mc, std::atomic< std::size_t > &counts_run, std::atomic< std::size_t > &counts_passed, std::atomic< std::size_t > &counts_failed)
Public Member Functions inherited from art::Modifier
	~Modifier () noexcept
	Modifier ()
	Modifier (Modifier const &)=delete
	Modifier (Modifier &&)=delete
Modifier &	operator= (Modifier const &)=delete
Modifier &	operator= (Modifier &&)=delete
Public Member Functions inherited from art::ModuleBase
virtual	~ModuleBase () noexcept
	ModuleBase ()
ModuleDescription const &	moduleDescription () const
void	setModuleDescription (ModuleDescription const &)
std::array< std::vector< ProductInfo >, NumBranchTypes > const &	getConsumables () const
void	sortConsumables (std::string const &current_process_name)
template<typename T , BranchType BT>
ViewToken< T >	consumesView (InputTag const &tag)
template<typename T , BranchType BT>
ViewToken< T >	mayConsumeView (InputTag const &tag)

Private Member Functions
void	removeHitsAssignedToTracks (std::vector< art::Ptr< recob::Hit > > &hitlist, const std::vector< recob::Track > &tracks, const art::FindManyP< recob::Hit > &fbp)
void	removeUnmatchedHitsCloseToTracks (detinfo::DetectorPropertiesData const &detProp, std::vector< art::Ptr< recob::Hit > > &hitlist, const std::vector< recob::Track > &tracks, const art::FindManyP< recob::Hit > &fbp)
bool	isCloseToTrack (TVector2 p, const recob::Track &trk, unsigned int view, unsigned int tpc, unsigned int cryo)

Static Private Member Functions
static double	getDist2 (const TVector2 &psrc, const TVector2 &p0, const TVector2 &p1)

Private Attributes
std::string	fHitModuleLabel
std::string	fTrk3DModuleLabel

Additional Inherited Members
Public Types inherited from art::EDProducer
using	ModuleType = EDProducer
using	WorkerType = WorkerT< EDProducer >
Public Types inherited from art::detail::Producer
template<typename UserConfig , typename KeysToIgnore = void>
using	Table = Modifier::Table< UserConfig, KeysToIgnore >
Public Types inherited from art::Modifier
template<typename UserConfig , typename UserKeysToIgnore = void>
using	Table = ProducerTable< UserConfig, detail::ModuleConfig, UserKeysToIgnore >
Static Public Member Functions inherited from art::EDProducer
static void	commitEvent (EventPrincipal &ep, Event &e)
Protected Member Functions inherited from art::ModuleBase
ConsumesCollector &	consumesCollector ()
template<typename T , BranchType = InEvent>
ProductToken< T >	consumes (InputTag const &)
template<typename Element , BranchType = InEvent>
ViewToken< Element >	consumesView (InputTag const &)
template<typename T , BranchType = InEvent>
void	consumesMany ()
template<typename T , BranchType = InEvent>
ProductToken< T >	mayConsume (InputTag const &)
template<typename Element , BranchType = InEvent>
ViewToken< Element >	mayConsumeView (InputTag const &)
template<typename T , BranchType = InEvent>
void	mayConsumeMany ()

Detailed Description

Definition at line 47 of file EmLikeHits_module.cc.

Constructor & Destructor Documentation

dune::EmLikeHits::EmLikeHits ( fhicl::ParameterSet const &  p )

explicit

Definition at line 97 of file EmLikeHits_module.cc.

                                                  : EDProducer{p}
{
        this->reconfigure(p);
        produces< std::vector<recob::Hit> >();
}

dune::EmLikeHits::EmLikeHits ( EmLikeHits const &  )

delete

dune::EmLikeHits::EmLikeHits ( EmLikeHits &&  )

delete

Member Function Documentation

double dune::EmLikeHits::getDist2 ( const TVector2 &  psrc, const TVector2 &  p0, const TVector2 &  p1  )

staticprivate

Definition at line 236 of file EmLikeHits_module.cc.

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

        TVector2 v2(psrc); v2 -= p1;
        TVector2 v3(v1);   v3 *= -1.0;

        double v0Norm2 = v0.Mod2();
        double v1Norm2 = v1.Mod2();

        double eps = 1.0E-7; // 0.001mm
        if (v1Norm2 > eps)
        {
                double mag01 = sqrt(v0Norm2 * v1Norm2);
                double cosine01 = 0.0;
                if (mag01 != 0.0) cosine01 = v0 * v1 / mag01;

                double v2Norm2 = v2.Mod2();
                double mag23 = sqrt(v2Norm2 * v3.Mod2());
                double cosine23 = 0.0;
                if (mag23 != 0.0) cosine23 = v2 * v3 / mag23;

                double result = 0.0;
                if ((cosine01 > 0.0) && (cosine23 > 0.0))
                {
                        result = (1.0 - cosine01 * cosine01) * v0Norm2;
                }
                else
                {
                        if (cosine01 <= 0.0) result = v0Norm2;
                        else result = v2Norm2;
                }

                if (result >= 0.0) return result;
                else return 0.0;
        }
        else
        {
                v1 = p0; v1 += p1; v1 *= 0.5;
                return pma::Dist2(v1, psrc);
        }
}

bool dune::EmLikeHits::isCloseToTrack ( TVector2  p, const recob::Track &  trk, unsigned int  view, unsigned int  tpc, unsigned int  cryo  )

private

Definition at line 137 of file EmLikeHits_module.cc.

{
        art::ServiceHandle<geo::Geometry> geom;
        double wirePitch = geom->TPC(tpc, cryo).Plane(view).WirePitch();

        //double driftPitch = detProp.GetXTicksCoefficient(tpc, cryo);

        double max_d2_d = 0.3 * 0.3;
        double max_d2_w = (wirePitch + 0.1) * (wirePitch + 0.1);

        bool isClose = false;
        for (size_t i = 0; i < trk.NumberTrajectoryPoints() - 1; ++i)
        {
                TVector2 p0 = pma::GetVectorProjectionToPlane(trk.LocationAtPoint<TVector3>(i), view, tpc, cryo);
                TVector2 p1 = pma::GetVectorProjectionToPlane(trk.LocationAtPoint<TVector3>(i + 1), view, tpc, cryo);
                double d2 = getDist2(p, p0, p1);

                double dpx = fabs(p0.X() - p1.X());
                double dpy = fabs(p0.Y() - p1.Y());

                if (((dpx > 0.5 * dpy) && (d2 < max_d2_d)) || (d2 < max_d2_w))
                {
                        isClose = true; break;
                }
        }
        return isClose;
}

EmLikeHits& dune::EmLikeHits::operator= ( EmLikeHits const &  )

delete

EmLikeHits& dune::EmLikeHits::operator= ( EmLikeHits &&  )

delete

void dune::EmLikeHits::produce ( art::Event &  e )

overridevirtual

Implements art::EDProducer.

Definition at line 209 of file EmLikeHits_module.cc.

{
        std::unique_ptr< std::vector< recob::Hit > > not_track_hits(new std::vector< recob::Hit >);

        std::vector< art::Ptr<recob::Hit> > hitlist;
        auto hitListHandle = evt.getHandle< std::vector<recob::Hit> >(fHitModuleLabel);
        auto trkListHandle = evt.getHandle< std::vector<recob::Track> >(fTrk3DModuleLabel);

        if (hitListHandle && trkListHandle)
        {
                art::FindManyP< recob::Hit > fbp(trkListHandle, evt, fTrk3DModuleLabel);

                art::fill_ptr_vector(hitlist, hitListHandle);
                mf::LogVerbatim("EmLikeHits") << "all hits: " << hitlist.size() << std::endl;

                removeHitsAssignedToTracks(hitlist, *trkListHandle, fbp);
                auto const detProp =
                  art::ServiceHandle<detinfo::DetectorPropertiesService const>()->DataFor(evt);
                removeUnmatchedHitsCloseToTracks(detProp, hitlist, *trkListHandle, fbp);

                for (auto const& hit : hitlist) not_track_hits->push_back(recob::Hit(*hit));
                mf::LogVerbatim("EmLikeHits") << "remaining not track-like hits: " << not_track_hits->size() << std::endl;
        }
        evt.put(std::move(not_track_hits));
}

void dune::EmLikeHits::reconfigure

(

fhicl::ParameterSet const &

p

)

Definition at line 104 of file EmLikeHits_module.cc.

{
        fHitModuleLabel = pset.get< std::string >("HitModuleLabel");
        fTrk3DModuleLabel = pset.get< std::string >("Trk3DModuleLabel");
}

void dune::EmLikeHits::removeHitsAssignedToTracks ( std::vector< art::Ptr< recob::Hit > > &  hitlist, const std::vector< recob::Track > &  tracks, const art::FindManyP< recob::Hit > &  fbp  )

private

Definition at line 111 of file EmLikeHits_module.cc.

{
        for (size_t t = 0; t < tracks.size(); t++)
                if (!(tracks[t].ID() & 0x10000))
        {
                std::vector< art::Ptr<recob::Hit> > v = fbp.at(t);
                mf::LogVerbatim("EmLikeHits") << "   track-like trajectory: " << v.size() << std::endl;
                size_t ih = 0;
                while (ih < hitlist.size())
                {
                        bool found = false;
                        for (size_t it = 0; it < v.size(); it++)
                                if (v[it].key() == hitlist[ih].key())
                        {
                                found = true; break;
                        }
                        if (found) hitlist.erase(hitlist.begin() + ih);
                        else ih++;
                }
        }
}

void dune::EmLikeHits::removeUnmatchedHitsCloseToTracks ( detinfo::DetectorPropertiesData const &  detProp, std::vector< art::Ptr< recob::Hit > > &  hitlist, const std::vector< recob::Track > &  tracks, const art::FindManyP< recob::Hit > &  fbp  )

private

Definition at line 167 of file EmLikeHits_module.cc.

{
        size_t ih = 0;
        while (ih < hitlist.size())
        {
                bool unmatched = true, close = false;;
                for (size_t t = 0; t < tracks.size(); t++)
                {
                        std::vector< art::Ptr<recob::Hit> > v = fbp.at(t);
                        for (size_t it = 0; it < v.size(); it++)
                                if (v[it].key() == hitlist[ih].key())
                                {
                                        unmatched = false; break;
                                }
                        if (!unmatched) break;
                }
                if (unmatched)
                {
                        unsigned int plane = hitlist[ih]->WireID().Plane;
                        unsigned int tpc = hitlist[ih]->WireID().TPC;
                        unsigned int cryo = hitlist[ih]->WireID().Cryostat;

                        TVector2 hcm = pma::WireDriftToCm(detProp,
                                hitlist[ih]->WireID().Wire, hitlist[ih]->PeakTime(), plane, tpc, cryo);

                        for (size_t t = 0; t < tracks.size(); t++)
                                if (!(tracks[t].ID() & 0x10000) &&
                                    isCloseToTrack(hcm, tracks[t], plane, tpc, cryo))
                        {
                                close = true; break;
                        }
                }
                if (close) hitlist.erase(hitlist.begin() + ih);
                else ih++;
        }
}

Member Data Documentation

std::string dune::EmLikeHits::fHitModuleLabel

private

Definition at line 91 of file EmLikeHits_module.cc.

std::string dune::EmLikeHits::fTrk3DModuleLabel

private

Definition at line 92 of file EmLikeHits_module.cc.

---

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

• dunereco/dunereco/HitFinderDUNE/EmLikeHits_module.cc