evd::HitSelector Class Reference

#include <HitSelector.h>

Public Member Functions
	HitSelector ()
void	SaveHits (const art::Event &evt, unsigned int plane, double x, double y, double x1, double y1, double distance, bool good_plane=true)
double	SaveSeedLines (const art::Event &evt, std::vector< util::PxLine > seedline, double distance)
void	ChangeHit (const art::Event &evt, unsigned int plane, double x, double y)
std::vector< const recob::Hit * >	GetSelectedHits (unsigned int plane)
std::vector< art::Ptr < recob::Hit > >	GetSelectedHitPtrs (unsigned int plane)
void	ClearHitList (unsigned int plane)
std::vector< recob::Seed > &	SeedVector ()

Private Attributes
std::vector< recob::Seed >	fSeedVector
std::vector< std::vector < double > >	starthitout
std::vector< std::vector < double > >	endhitout

Detailed Description

Definition at line 25 of file HitSelector.h.

Constructor & Destructor Documentation

HitSelector::HitSelector

(

)

Definition at line 34 of file HitSelector.cxx.

  {
    starthitout.resize(3);
    endhitout.resize(3);
  }

Member Function Documentation

void HitSelector::ChangeHit	(	const art::Event &	evt,
		unsigned int	plane,
		double	xin,
		double	yin
	)

Save HitList to infoTransfer service

Parameters

evt	: Event handle to get data objects from
view	: Pointer to view to draw on
plane	: plane number of view

Definition at line 187 of file HitSelector.cxx.

  {
    art::ServiceHandle<evd::RecoDrawingOptions const> recoOpt;
    art::ServiceHandle<geo::Geometry const> geo;
    auto const clockData = art::ServiceHandle<detinfo::DetectorClocksService const>()->DataFor(evt);
    auto const detProp =
      art::ServiceHandle<detinfo::DetectorPropertiesService const>()->DataFor(evt, clockData);
    util::GeometryUtilities const gser{*geo, clockData, detProp};

    //get hits from info transfer, see if our selected hit is in it
    std::vector<art::Ptr<recob::Hit>> hits_saved;
    std::vector<art::Ptr<recob::Hit>> hitlist;

    double x = xin * gser.WireToCm();
    double y = yin * gser.TimeToCm();

    art::Handle<std::vector<recob::Hit>> HitListHandle;

    for (size_t imod = 0; imod < recoOpt->fHitLabels.size(); ++imod) {
      art::InputTag const which = recoOpt->fHitLabels[imod];
      evt.getByLabel(which, HitListHandle);

      for (unsigned int ii = 0; ii < HitListHandle->size(); ++ii) {
        art::Ptr<recob::Hit> hit(HitListHandle, ii);
        if (hit->WireID().Plane == plane) hitlist.push_back(hit);
      }

      util::PxHitConverter PxC{gser};
      std::vector<util::PxHit> pxhitlist;
      PxC.GeneratePxHit(hitlist, pxhitlist);

      unsigned int hitindex = gser.FindClosestHitIndex(pxhitlist, x, y);
      if (hitlist[hitindex].isNull() || (hitindex > hitlist.size())) {
        WriteMsg("no luck finding hit in evd, please try again");
        break;
      }

      art::ServiceHandle<evd::InfoTransfer> infot;
      hits_saved = infot->GetSelectedHitList(plane);
      int found_it = 0;
      for (unsigned int jj = 0; jj < hits_saved.size(); ++jj) {
        if (hitlist[hitindex]->PeakTime() == hits_saved[jj]->PeakTime()) {
          if (hitlist[hitindex]->Channel() == hits_saved[jj]->Channel()) {
            found_it = 1;
            hits_saved.erase(hits_saved.begin() + jj);
          }
        }
      }

      //if didn't find it, add it
      if (found_it != 1) { hits_saved.push_back(hitlist[hitindex]); }

      //update the info transfer list
      infot->SetHitList(plane, hits_saved);
      infot->SetTestFlag(1);
      infot->SetEvtNumber(evt.id().event());
    }
  }

void HitSelector::ClearHitList

(

unsigned int

plane

)

Definition at line 278 of file HitSelector.cxx.

  {
    art::ServiceHandle<evd::InfoTransfer> infot;
    infot->ClearSelectedHitList(plane);
  }

std::vector< art::Ptr< recob::Hit > > HitSelector::GetSelectedHitPtrs

(

unsigned int

plane

)

Save HitList to infoTransfer service

Parameters

evt	: Event handle to get data objects from
view	: Pointer to view to draw on
plane	: plane number of view

Definition at line 255 of file HitSelector.cxx.

  {
    art::ServiceHandle<evd::InfoTransfer const> infot;
    return infot->GetSelectedHitList(plane);
  }

std::vector< const recob::Hit * > HitSelector::GetSelectedHits

(

unsigned int

plane

)

Definition at line 262 of file HitSelector.cxx.

  {
    std::vector<art::Ptr<recob::Hit>> hits_saved;
    std::vector<const recob::Hit*> hits_to_draw; //draw selected hits in a different color

    art::ServiceHandle<evd::InfoTransfer const> infot;
    hits_saved = infot->GetSelectedHitList(plane);

    for (unsigned int i = 0; i < hits_saved.size(); i++)
      hits_to_draw.push_back(hits_saved[i].get());

    return hits_to_draw;
  }

void HitSelector::SaveHits	(	const art::Event &	evt,
		unsigned int	plane,
		double	xin,
		double	yin,
		double	x1in,
		double	y1in,
		double	distance,
		bool	good_plane = true
	)

Save HitList to infoTransfer service

Parameters

evt	: Event handle to get data objects from
view	: Pointer to view to draw on
plane	: plane number of view

Definition at line 88 of file HitSelector.cxx.

  {
    art::ServiceHandle<evd::RecoDrawingOptions const> recoOpt;
    art::ServiceHandle<geo::Geometry const> geo;
    auto const clockData = art::ServiceHandle<detinfo::DetectorClocksService const>()->DataFor(evt);
    auto const detProp =
      art::ServiceHandle<detinfo::DetectorPropertiesService const>()->DataFor(evt, clockData);
    util::GeometryUtilities const gser{*geo, clockData, detProp};
    std::vector<art::Ptr<recob::Hit>> hits_to_save;

    art::Handle<std::vector<recob::Hit>> HitListHandle;

    starthitout[plane].clear();
    endhitout[plane].clear();

    starthitout[plane].resize(2);
    endhitout[plane].resize(2);

    //convert input variables to cm-cm space required by GeometryUtilities
    double x = xin * gser.WireToCm();
    double x1 = x1in * gser.WireToCm();
    double y = yin * gser.TimeToCm();
    double y1 = y1in * gser.TimeToCm();

    double lslope = 0;

    if ((x - x1) != 0) lslope = (y - y1) / (x - x1);

    for (size_t imod = 0; imod < recoOpt->fHitLabels.size(); ++imod) {
      art::InputTag const which = recoOpt->fHitLabels[imod];

      std::vector<art::Ptr<recob::Hit>> hitlist;
      hitlist.clear();
      evt.getByLabel(which, HitListHandle);

      for (unsigned int ii = 0; ii < HitListHandle->size(); ++ii) {
        art::Ptr<recob::Hit> hit(HitListHandle, ii);
        if (hit->WireID().Plane == plane) hitlist.push_back(hit);
      }

      // Select Local Hit List
      util::PxHitConverter PxC{gser};
      std::vector<util::PxHit> pxhitlist;
      PxC.GeneratePxHit(hitlist, pxhitlist);
      std::vector<unsigned int> pxhitlist_local_index;
      std::vector<util::PxHit> pxhitlist_local;
      pxhitlist_local.clear();

      util::PxPoint startHit;
      startHit.plane = pxhitlist.at(0).plane;
      startHit.w = (x + x1) / 2;
      startHit.t = (y + y1) / 2;

      double orttemp = std::hypot(y1 - y, x1 - x) / 2;

      gser.SelectLocalHitlistIndex(
        pxhitlist, pxhitlist_local_index, startHit, orttemp, distance, lslope);

      for (unsigned int idx = 0; idx < pxhitlist_local_index.size(); idx++) {
        hits_to_save.push_back(hitlist.at(pxhitlist_local_index.at(idx)));
        pxhitlist_local.push_back(pxhitlist.at(pxhitlist_local_index.at(idx)));
      }

      auto const hit_index = gser.FindClosestHitIndex(pxhitlist_local, x, y);
      recob::Hit const& hit = *hits_to_save[hit_index];
      starthitout[plane][1] = hit.PeakTime();
      starthitout[plane][0] = hit.WireID().Wire;

      recob::Hit const& endhit = hit; // this obviously not correct, the fact that x,y are used for both start and end point, A.S. -> to debug
      endhitout[plane][1] = endhit.PeakTime();
      endhitout[plane][0] = endhit.WireID().Wire;
    }

    art::ServiceHandle<evd::InfoTransfer> infot;
    infot->SetHitList(plane, hits_to_save);
    infot->SetTestFlag(1);
    if (good_plane) {
      infot->SetStartHitCoords(plane, starthitout[plane]);
      infot->SetEndHitCoords(plane, endhitout[plane]);
    }
    infot->SetEvtNumber(evt.id().event());
  }

double HitSelector::SaveSeedLines	(	const art::Event &	evt,
		std::vector< util::PxLine >	seedlines,
		double	distance
	)

Save SeedLineList to infoTransfer service

Parameters

evt	: Event handle to get data objects from
view	: Pointer to view to draw on

Return value is the kinetic enegry of the track

Definition at line 50 of file HitSelector.cxx.

  {
    art::ServiceHandle<evd::InfoTransfer> infot;
    art::ServiceHandle<geo::Geometry const> geom;
    art::ServiceHandle<evd::RecoDrawingOptions const> recoOpt;
    std::map<int, std::vector<art::Ptr<recob::Hit>>> hits_to_save;

    double KineticEnergy = 0;
    art::Handle<std::vector<recob::Hit>> HitListHandle;

    for (size_t imod = 0; imod < recoOpt->fHitLabels.size(); ++imod) {
      art::InputTag const which = recoOpt->fHitLabels[imod];
      evt.getByLabel(which, HitListHandle);
    }

    infot->SetSeedList(seedlines);
    for (std::map<int, std::vector<art::Ptr<recob::Hit>>>::iterator it = hits_to_save.begin();
         it != hits_to_save.end();
         ++it) {
      infot->SetHitList(it->first, it->second);
    }
    infot->SetTestFlag(1);
    infot->SetEvtNumber(evt.id().event());

    return KineticEnergy;
  }

std::vector< recob::Seed > & HitSelector::SeedVector

(

)

Definition at line 286 of file HitSelector.cxx.

  {
    return fSeedVector;
  }

Member Data Documentation

std::vector<std::vector<double> > evd::HitSelector::endhitout

private

Definition at line 55 of file HitSelector.h.

std::vector<recob::Seed> evd::HitSelector::fSeedVector

private

Definition at line 52 of file HitSelector.h.

std::vector<std::vector<double> > evd::HitSelector::starthitout

private

Definition at line 54 of file HitSelector.h.

---

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

• HitSelector.h
• HitSelector.cxx