evd::HitSelector Class Reference

#include "HitSelector.h"

Public Member Functions
	HitSelector ()
	~HitSelector ()
void	SaveHits (const art::Event &evt, evdb::View2D *view, unsigned int plane, double x, double y, double x1, double y1, double distance, bool good_plane=true)
double	SaveSeedLines (const art::Event &evt, evdb::View2D *view, std::vector< util::PxLine > seedline, double distance)
void	ChangeHit (const art::Event &evt, evdb::View2D *view, unsigned int plane, double x, double y)
std::vector< const recob::Hit * >	GetSelectedHits (unsigned int plane)
trkf::HitPtrVec	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 47 of file HitSelector.h.

Constructor & Destructor Documentation

evd::HitSelector::HitSelector

(

)

Definition at line 34 of file HitSelector.cxx.

References endhitout, and starthitout.

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

evd::HitSelector::~HitSelector

(

)

Definition at line 41 of file HitSelector.cxx.

  {
 
  }

Member Function Documentation

void evd::HitSelector::ChangeHit	(	const art::Event &	evt,
		evdb::View2D *	view,
		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 303 of file HitSelector.cxx.

References art::EventID::event(), evd::RecoDrawingOptions::fHitLabels, util::GeometryUtilities::FindClosestHitIndex(), util::PxHitConverter::GeneratePxHit(), art::DataViewImpl::getByLabel(), evd::InfoTransfer::GetSelectedHitList(), art::Event::id(), geo::PlaneID::Plane, evd::InfoTransfer::SetEvtNumber(), evd::InfoTransfer::SetHitList(), evd::InfoTransfer::SetTestFlag(), util::GeometryUtilities::TimeToCm(), recob::Hit::WireID(), util::GeometryUtilities::WireToCm(), x, and y.

Referenced by evd::TWireProjPad::SelectOneHit().

  {
    art::ServiceHandle<evd::RecoDrawingOptions>  recoOpt;
    art::ServiceHandle<geo::Geometry>  geo;
    util::GeometryUtilities  gser;
    
    
    //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;
   // std::vector < art::Ptr < recob::Hit> > hits_to_save; 
     hitlist.clear();

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

    art::Handle< std::vector<recob::Hit> > HitListHandle;
   
    //preprocess event - load up all the hits with std::vector, as in BackTracker
    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;
     std::vector <util::PxHit> pxhitlist;
     PxC.GeneratePxHit(hitlist,pxhitlist);
      
//       art::Ptr < recob::Hit > selected_hit= hitlist[gser.FindClosestHitIndex(pxhitlist,x,y)];
      unsigned int hitindex=gser.FindClosestHitIndex(pxhitlist,x,y);
      //FindClosestHitPtr(hitlist, x, y); 
      
      // find selected hit in evD
     //art::Ptr<recob::Hit> selected_hit;
//       for(unsigned int ii = 0; ii < HitListHandle->size(); ++ii){
//      art::Ptr<recob::Hit> hit(HitListHandle, ii);
//      double time=hit->PeakTime();
//      if(hit->WireID().Plane != plane)
//        continue;
//      if(hit->WireID().Wire < (x+1) &&
//         hit->WireID().Wire > (x-1) ){
//        if(y>(time-5) && y<(time+5)){
//          selected_hit=hit;
//        }
//      }
                  
    //}
      // for c2: unsigned int hitindex cannot be < 0
      //if(hitlist[hitindex].isNull() || (hitindex<0 || hitindex > hitlist.size())){
      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 evd::HitSelector::ClearHitList

(

unsigned int

plane

)

Definition at line 432 of file HitSelector.cxx.

References evd::InfoTransfer::ClearSelectedHitList().

Referenced by evd::TWireProjPad::ClearHitList().

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

trkf::HitPtrVec evd::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 402 of file HitSelector.cxx.

References evd::InfoTransfer::GetSelectedHitList(), and trkf::HitPtrVec::Hits.

Referenced by evd::CalorPad::Draw().

  {
    art::ServiceHandle<evd::InfoTransfer>   infot;
    trkf::HitPtrVec ToReturn;
    ToReturn.Hits = infot->GetSelectedHitList(plane);
    return ToReturn;
  }

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

(

unsigned int

plane

)

Definition at line 410 of file HitSelector.cxx.

References evd::InfoTransfer::GetSelectedHitList().

  {
    std::vector < art::Ptr<recob::Hit> > hits_saved;
    //std::vector < art::Ptr < recob::Hit> > hits_to_save; 

    //std::vector < recob::Hit> hits_to_point; //needed to convert art::Ptr<recob:: Hit> to Hit* to pass to Hit2D
    std::vector< const recob::Hit*> hits_to_draw; //draw selected hits in a different color
    
    art::ServiceHandle<evd::InfoTransfer>   infot;
    hits_saved=infot->GetSelectedHitList(plane);
        
    // if(hits_saved.size()!=0){
    for(unsigned int i=0;i<hits_saved.size();i++)
      hits_to_draw.push_back(hits_saved[i].get());
    //}
    
    return hits_to_draw;     
    
  }

void evd::HitSelector::SaveHits	(	const art::Event &	evt,
		evdb::View2D *	view,
		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 170 of file HitSelector.cxx.

References endhitout, art::EventID::event(), evd::RecoDrawingOptions::fHitLabels, util::GeometryUtilities::FindClosestHitIndex(), util::PxHitConverter::GeneratePxHit(), art::DataViewImpl::getByLabel(), art::Event::id(), recob::Hit::PeakTime(), util::PxPoint::plane, geo::PlaneID::Plane, util::GeometryUtilities::SelectLocalHitlistIndex(), evd::InfoTransfer::SetEndHitCoords(), evd::InfoTransfer::SetEvtNumber(), evd::InfoTransfer::SetHitList(), evd::InfoTransfer::SetStartHitCoords(), evd::InfoTransfer::SetTestFlag(), starthitout, util::PxPoint::t, util::GeometryUtilities::TimeToCm(), util::PxPoint::w, geo::WireID::Wire, recob::Hit::WireID(), util::GeometryUtilities::WireToCm(), x, x1, y, and y1.

Referenced by evd::TWireProjPad::SaveHitList().

  {
    // unsigned int w  = 0;
    art::ServiceHandle<evd::RecoDrawingOptions>  recoOpt;
    art::ServiceHandle<geo::Geometry>  geo;
    util::GeometryUtilities  gser;
    std::vector < art::Ptr < recob::Hit> > hits_to_save; //std::vector < azart::Ptr >


    //std::vector < recob::Hit> hits_to_point; //needed to convert art::Ptr<recob:: Hit> to Hit* to pass to Hit2D
   // std::vector< const recob::Hit*> hits_to_draw; //draw selected hits in a different color

    art::Handle< std::vector<recob::Hit> > HitListHandle;
    //  recob::Hit * startHit;
    //  recob::Hit * endHit;
    //preprocess event - load up all the hits with std::vector, as in BackTracker
    
          
                  
    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;
     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=TMath::Sqrt((y1-y)*(y1-y)*gser.TimeToCm()*gser.TimeToCm() + (x1-x)*(x1-x)*gser.WireToCm()*gser.WireToCm())/2;
      double orttemp=TMath::Sqrt((y1-y)*(y1-y) + (x1-x)*(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)));
       }
    
    
   /*  gser.SelectLocalHitlist( hitlist, 
                         hits_to_save,
                           (x+x1)/2,
                           (y+y1)/2, 
                           TMath::Sqrt((y1-y)*(y1-y)*gser.TimeToCm()*gser.TimeToCm() + (x1-x)*(x1-x)*gser.WireToCm()*gser.WireToCm())/2,   
                           distance, 
                           lslope);*/           
    //const_cast<recob::Hit *> (nearHit.get())
    
     recob::Hit * hit=  const_cast<recob::Hit *> ((hits_to_save[gser.FindClosestHitIndex(pxhitlist_local,x,y)]).get());
     //gser.FindClosestHit(hits_to_save,
        //                    x, y);
      
       starthitout[plane][1] = hit->PeakTime() ;  
       starthitout[plane][0] = hit->WireID().Wire;
 
//        // this obviously not correct, the fact that x,y are used for both start and end point, A.S. -> to debug
       recob::Hit * endhit= const_cast<recob::Hit *> ((hits_to_save[gser.FindClosestHitIndex(pxhitlist_local,x,y)]).get());
       //FindClosestHit(hits_to_save,
        //                    x, y);
      
        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 evd::HitSelector::SaveSeedLines	(	const art::Event &	evt,
		evdb::View2D *	view,
		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 56 of file HitSelector.cxx.

References art::PtrVector< T >::at(), d, art::EventID::event(), evd::RecoDrawingOptions::fCaloPSet, evd::RecoDrawingOptions::fHitLabels, evd::EvdLayoutOptions::fMakeSeeds, art::DataViewImpl::getByLabel(), trkf::BezierTrack::GetCalorimetryObject(), trkf::BezierTrack::GetClosestApproaches(), art::Event::id(), geo::kCollection, anab::Calorimetry::KineticEnergy(), art::PtrVector< T >::push_back(), s, SeedVector(), evd::InfoTransfer::SetEvtNumber(), evd::InfoTransfer::SetHitList(), evd::InfoTransfer::SetSeedList(), evd::InfoTransfer::SetTestFlag(), art::PtrVector< T >::size(), and lar::dump::vector().

Referenced by evd::TWireProjPad::SaveSeedList().

  {
    art::ServiceHandle<evd::InfoTransfer>   infot;
    art::ServiceHandle<geo::Geometry>       geom;
    art::ServiceHandle<evd::RecoDrawingOptions>  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);


      art::ServiceHandle<evd::EvdLayoutOptions> evdlayoutoptions;
      if(evdlayoutoptions->fMakeSeeds)
        {
          // Code to collect hits around curve

          // First, put the hits into a PtrVector
          art::PtrVector<recob::Hit> HitVec;
          for(unsigned int i=0; i < HitListHandle->size(); ++i)
            {
              art::Ptr<recob::Hit> hit(HitListHandle,i);
              HitVec.push_back(hit);
            }

          // Get the bezier curve which was created
          if(SeedVector().size()>1)
            {
              mf::LogVerbatim("HitSelector") <<"Hit selector passing bezier track " 
                                             << SeedVector().size() << " seeds";
            }
          else if(SeedVector().size()==1)
            {
              mf::LogVerbatim("HitSelector") <<"One seed defined - splitting in half to make track";
              double Pt[3], Dir[3], Err[3];
              SeedVector().at(0).GetPoint(     Pt,  Err);
              SeedVector().at(0).GetDirection( Dir, Err);
             
              double NewPt1[3], NewPt2[3], NewDir[3];
              for(int i=0; i!=3; ++i)
                {
                  NewPt1[i] = Pt[i] - Dir[i]/2;
                  NewPt2[i] = Pt[i] + Dir[i]/2;
                  NewDir[i] = Dir[i] / 2;
                }
              SeedVector().clear();
              
              recob::Seed TheSeed1(NewPt1, NewDir, Err, Err);
              SeedVector().push_back(TheSeed1);
              recob::Seed TheSeed2(NewPt2, NewDir, Err, Err);
              SeedVector().push_back(TheSeed2);
            }
          else
            {
              mf::LogWarning("HitSelector") << "Cannot select hits, no seeds selected!";
              return 0;
            }
           
          trkf::BezierTrack BTrack(SeedVector());

          // These hold output,
          //  s = fractional distance along curve of close approach
          //  d = min distance from line
          std::vector<double> s;
          std::vector<double> d;

          // Run through the hits, getting their distances
          BTrack.GetClosestApproaches(HitVec, s, d);
         
          // Put them into the vector for output
          for(size_t ih=0; ih!=HitVec.size(); ++ih){
            if(d.at(ih)<distance){
              hits_to_save[HitVec.at(ih)->WireID().Plane].push_back(HitVec.at(ih));
            } 
          }
          art::ServiceHandle<evd::RecoDrawingOptions> recoOpt;
          calo::CalorimetryAlg calalg(recoOpt->fCaloPSet);
          for(std::map<int, std::vector<art::Ptr<recob::Hit> > >::iterator it=hits_to_save.begin(); it!=hits_to_save.end(); ++it){
            KineticEnergy += BTrack.GetCalorimetryObject(it->second, geo::kCollection,calalg).KineticEnergy();
            
          }
          mf::LogVerbatim("HitSelector") <<"Track kinetic energy : " << KineticEnergy;
        }
    }
   
   
    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 > & evd::HitSelector::SeedVector

(

)

Definition at line 439 of file HitSelector.cxx.

References fSeedVector.

Referenced by SaveSeedLines(), and evd::Display3DPad::UpdateSeedCurve().

  {
    return fSeedVector;
  }

Member Data Documentation

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

private

Definition at line 86 of file HitSelector.h.

Referenced by HitSelector(), and SaveHits().

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

private

Definition at line 83 of file HitSelector.h.

Referenced by SeedVector().

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

private

Definition at line 85 of file HitSelector.h.

Referenced by HitSelector(), and SaveHits().

---

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

• lareventdisplay/v06_22_00/source/lareventdisplay/EventDisplay/HitSelector.h
• lareventdisplay/v06_22_00/source/lareventdisplay/EventDisplay/HitSelector.cxx