vertex::AggregateVertex Class Reference

Inheritance diagram for vertex::AggregateVertex:

Public Member Functions
	AggregateVertex (fhicl::ParameterSet const &pset)
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 ()
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	produce (art::Event &evt) override
std::unique_ptr< std::vector< recob::Vertex > >	MatchV2T (art::Event &evt, art::Assns< recob::Vertex, recob::Track > &vtassn, art::Assns< recob::Vertex, recob::Shower > &vsassn, art::Assns< recob::Vertex, recob::Hit > &vhassn)

Private Attributes
std::string	fDBScanModuleLabel
std::string	fHoughModuleLabel
std::string	fTrack3DModuleLabel
std::string	fEndPointModuleLabel
art::PtrVector< recob::EndPoint2D >	feplist
art::PtrVector< recob::Track >	ftracklist
art::PtrVector< recob::EndPoint2D >	feplistStrong

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 43 of file AggregateVertex_module.cc.

Constructor & Destructor Documentation

vertex::AggregateVertex::AggregateVertex ( fhicl::ParameterSet const &  pset )

explicit

Definition at line 74 of file AggregateVertex_module.cc.

    : EDProducer{pset}
    , fDBScanModuleLabel  (pset.get< std::string >("DBScanModuleLabel"  ))
    , fHoughModuleLabel   (pset.get< std::string >("HoughModuleLabel"   ))
    , fTrack3DModuleLabel (pset.get< std::string >("Track3DModuleLabel" ))
    , fEndPointModuleLabel(pset.get< std::string >("EndPointModuleLabel"))
  {
    produces< std::vector<recob::Vertex> >();
    produces< art::Assns<recob::Vertex, recob::Hit> >();
    produces< art::Assns<recob::Vertex, recob::Track> >();
    produces< art::Assns<recob::Vertex, recob::Shower> >();
  }

Member Function Documentation

std::unique_ptr< std::vector< recob::Vertex > > vertex::AggregateVertex::MatchV2T ( art::Event &  evt, art::Assns< recob::Vertex, recob::Track > &  vtassn, art::Assns< recob::Vertex, recob::Shower > &  vsassn, art::Assns< recob::Vertex, recob::Hit > &  vhassn  )

private

Todo:

Really need to also determine the xyz position of the found vertex

Todo:

Also should determine the ID for this vertex

Todo:

Really need to also determine the xyz position of the found vertex

Todo:

Also should determine the ID for this vertex

Todo:

get a shower module label in line 169 and get a list of showers

Todo:

from the event

Definition at line 132 of file AggregateVertex_module.cc.

  {
    mf::LogInfo("AggregateVertex") << "AggregateEvent::MatchV2T(): (strong) vertexlistStrong"
                                   << " and tracklist lengths are "
                                   << feplistStrong.size() << " and " <<  ftracklist.size() << ".";

    art::FindManyP<recob::Hit> fmhst(feplistStrong, evt, fEndPointModuleLabel);

    // Bail if there are no tracks or vertices.
    //  if (!((int)vertexlistStrong.size()) || !((int)tracklist.size())) return NULL;
    if (feplistStrong.isNull() || ftracklist.isNull()) {
      return std::unique_ptr< std::vector<recob::Vertex> > (new std::vector<recob::Vertex>);
    }

    // Loop on the vertices, and all the hits in each
    std::unique_ptr< std::vector<recob::Vertex> > verts(new std::vector<recob::Vertex>);

    art::FindManyP<recob::Hit> fmht(ftracklist, evt, fTrack3DModuleLabel);
    art::FindManyP<recob::Hit> fmhs(ftracklist, evt, fTrack3DModuleLabel);

    for(size_t epctr = 0; epctr < feplistStrong.size(); ++epctr){
      art::PtrVector<recob::Track>  tlistAssoc; // Will fill with matching tracks.
      art::PtrVector<recob::Shower> slistAssoc; // Will fill with matching showers.
      std::vector<size_t> trkIdx;
      std::vector<size_t> shwIdx;

      std::vector< art::Ptr<recob::Hit> > hitvertexlistStrong = fmhst.at(epctr);

      // Should be just one hit per vtx, as per Josh, but we loop anyway.
      art::PtrVector<recob::Hit>::const_iterator hvIter = hitvertexlistStrong.begin();
      while (hvIter != hitvertexlistStrong.end())  {
        art::Ptr<recob::Hit> hitv = (*hvIter);

        // Now loop on the track hits
        for(size_t t = 0; t < ftracklist.size(); ++t){

          std::vector< art::Ptr<recob::Hit> > hittlist = fmht.at(t);

          std::vector< art::Ptr<recob::Hit> >::const_iterator htIter = hittlist.begin();
          while (htIter != hittlist.end())  {

            art::Ptr<recob::Hit> hitt = (*htIter);

            if (hitt==hitv){
              //std::cout << "AggregateEvent::MatchV2T(): BooYiggity! Vtx/Trk hit match! hitt, hitv= "
              //<< hitt << " " <<hitv << std::endl;
              //std::cout << "AggregateEvent::MatchV2T(): vtx, trk " << vtx<< " " <<trk << std::endl;
              tlistAssoc.push_back(ftracklist[t]);
              trkIdx.push_back(t);
              htIter = hittlist.end()-1; // jump to end of track hitlist, since we've satisfied match
            }
            htIter++;
          } // end hits associated to track
        } // end track loop
        hvIter++;
       }

      // Now if matching tracks were found for this vertex then create the recob::Vertex object.
      if (tlistAssoc.size() > 0){
        /// \todo Really need to also determine the xyz position of the found vertex
        /// \todo Also should determine the ID for this vertex
        double xyz[3] = {-999., -999., -999.};
        verts->push_back(recob::Vertex(xyz, verts->size()));

        // associate the tracks to the vertex
        util::CreateAssn(*this, evt, *verts, tlistAssoc, vtassn);

        //associate the track hits to the vertex
        for(auto const& t : trkIdx){
          std::vector< art::Ptr<recob::Hit> > hits = fmht.at(t);
          util::CreateAssn(*this, evt, *verts, hits, vhassn);
        }
      }// end if there are tracks to be associated

      if (slistAssoc.size() > 0){
        /// \todo Really need to also determine the xyz position of the found vertex
        /// \todo Also should determine the ID for this vertex
        double xyz[3] = {-999., -999., -999.};
        verts->push_back(recob::Vertex(xyz, verts->size()));

        // associate the showers to the vertex
        util::CreateAssn(*this, evt, *verts, slistAssoc, vsassn);

        //associate the shower hits to the vertex
        ///\todo get a shower module label in line 169 and get a list of showers
        ///\todo from the event
        for(auto const& s : shwIdx){
          std::vector< art::Ptr<recob::Hit> > hits = fmhs.at(s);
          util::CreateAssn(*this, evt, *verts, hits, vhassn);
        }
      }// end if there are showers to be associated

      //HnhinVtxes->Fill(hitvertexlistStrong.size(),1);

    } // end loop over end points

    return verts;

  } // end AggregateVertex::MatchV2T

void vertex::AggregateVertex::produce ( art::Event &  evt )

overrideprivatevirtual

Implements art::EDProducer.

Definition at line 88 of file AggregateVertex_module.cc.

  {

    art::Handle< std::vector<recob::EndPoint2D> > epListHandle;
    evt.getByLabel(fEndPointModuleLabel,epListHandle);
    for(size_t ii = 0; ii < epListHandle->size(); ++ii){
      art::Ptr<recob::EndPoint2D> ep(epListHandle, ii);
      feplist.push_back(ep); // class member
    }

    art::Handle< std::vector<recob::Track> > trackListHandle;
    evt.getByLabel(fTrack3DModuleLabel,trackListHandle);
    for(size_t ii = 0; ii < trackListHandle->size(); ++ii){
      art::Ptr<recob::Track> track(trackListHandle, ii);
      ftracklist.push_back(track); // class member
    }


    // Only match strong vertices to tracks.
    art::PtrVector<recob::EndPoint2D>::const_iterator epIter = feplist.begin();

    while (epIter != feplist.end()){
      art::Ptr <recob::EndPoint2D> ep = (*epIter);
      if (ep->ID() < 3) feplistStrong.push_back(ep); // -- cuz there's one less
      epIter++;
    }

    // We will suck up the hits out of each vertex and out of the tracks
    // and see if there's an overlap of a sufficient (>0) number of hits. If so,
    // call the track a match to the vtx. Then, .... stick the track pointer(s)
    // into the AggVertex object.  EC, 23-July-2010.
    std::unique_ptr< art::Assns<recob::Vertex, recob::Track>  > vtassn(new art::Assns<recob::Vertex, recob::Track>);
    std::unique_ptr< art::Assns<recob::Vertex, recob::Shower> > vsassn(new art::Assns<recob::Vertex, recob::Shower>);
    std::unique_ptr< art::Assns<recob::Vertex, recob::Hit>    > vhassn(new art::Assns<recob::Vertex, recob::Hit>);

    std::unique_ptr< std::vector<recob::Vertex> > vcol (MatchV2T(evt, *vtassn, *vsassn, *vhassn));

    evt.put(std::move(vcol));
    evt.put(std::move(vtassn));
    evt.put(std::move(vsassn));
    evt.put(std::move(vhassn));
  }

Member Data Documentation

std::string vertex::AggregateVertex::fDBScanModuleLabel

private

Definition at line 58 of file AggregateVertex_module.cc.

std::string vertex::AggregateVertex::fEndPointModuleLabel

private

Definition at line 61 of file AggregateVertex_module.cc.

art::PtrVector<recob::EndPoint2D> vertex::AggregateVertex::feplist

private

Definition at line 63 of file AggregateVertex_module.cc.

art::PtrVector<recob::EndPoint2D> vertex::AggregateVertex::feplistStrong

private

Definition at line 65 of file AggregateVertex_module.cc.

std::string vertex::AggregateVertex::fHoughModuleLabel

private

Definition at line 59 of file AggregateVertex_module.cc.

std::string vertex::AggregateVertex::fTrack3DModuleLabel

private

Definition at line 60 of file AggregateVertex_module.cc.

art::PtrVector<recob::Track> vertex::AggregateVertex::ftracklist

private

Definition at line 64 of file AggregateVertex_module.cc.

---

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

• larreco/larreco/VertexFinder/AggregateVertex_module.cc