recob::HitRefinerAssociator Class Reference

A class handling a collection of hits and its associations. More...

#include <HitCreator.h>

Inheritance diagram for recob::HitRefinerAssociator:

Public Member Functions
void	use_hits (std::unique_ptr< std::vector< recob::Hit >> &&srchits)
	Uses the specified collection as data product. More...
void	put_into (art::Event &)
	Moves the data into the event. More...
void	put_into ()
	Moves the data into the event. More...
Constructors
	HitRefinerAssociator (art::Event &event, art::InputTag const &HitModuleLabel, std::string instance_name="", bool doWireAssns=true, bool doRawDigitAssns=true)
	Constructor: sets instance name and whether to build associations. More...
	HitRefinerAssociator (art::Event &event, art::InputTag const &HitModuleLabel, bool doWireAssns, bool doRawDigitAssns=true)
	Constructor: sets instance name and whether to build associations. More...
Public Member Functions inherited from recob::HitAndAssociationsWriterBase
size_t	size () const
	Returns the number of hits currently in the collection. More...
void	put_into (art::Event &)
	Moves the data into the event. More...
void	put_into ()
	Moves the data into the event. More...
std::vector< recob::Hit > const &	peek () const
	Returns a read-only reference to the current list of hits. More...

Protected Member Functions
void	prepare_associations (std::vector< recob::Hit > const &srchits)
	Finds out the associations for the specified hits. More...
void	prepare_associations ()
	Finds out the associations for the current hits. More...
Protected Member Functions inherited from recob::HitAndAssociationsWriterBase
	HitAndAssociationsWriterBase (art::Event &event, std::string instance_name, bool doWireAssns, bool doRawDigitAssns)
	Constructor: sets instance name and whether to build associations. More...
HitPtr_t	CreatePtr (size_t index) const
	Creates an art pointer to the hit with the specified index. More...

Protected Attributes
art::InputTag	hits_label
	Label of the collection of hits. More...
Protected Attributes inherited from recob::HitAndAssociationsWriterBase
std::string	prod_instance
	Tame of the instance for data products. More...
std::unique_ptr< std::vector < recob::Hit > >	hits
	Collection of hits. More...
std::unique_ptr< art::Assns < recob::Wire, recob::Hit > >	WireAssns
	Associations with wires. More...
std::unique_ptr< art::Assns < raw::RawDigit, recob::Hit > >	RawDigitAssns
	Associations with raw digits. More...
art::Event *	event = nullptr
	Pointer to the event we are using. More...
art::PtrMaker< recob::Hit >	hitPtrMaker
	Tool to create hit pointers,. More...

Additional Inherited Members
Static Public Member Functions inherited from recob::HitAndAssociationsWriterBase
static void	declare_products (art::ProducesCollector &collector, std::string instance_name="", bool doWireAssns=true, bool doRawDigitAssns=true)
	Declares the hit products we are going to fill. More...
Protected Types inherited from recob::HitAndAssociationsWriterBase
using	HitPtr_t = art::Ptr< recob::Hit >
	Type of art pointer to Hit. More...

Detailed Description

A class handling a collection of hits and its associations.

---

Use this object if you already have a recob::Hit data product and another collection that is going to become a data product, and you simply want the new hits associated to the wire and digit with the same channel. No hit-to-hit association is attempted (that would be, incidentally, not supported by art): the data product is used to get all the associated wires and digits, then they are associated to the new hits by channel ID. If a channel is not available, a warning is produced. If different hits on the same channel are associated to different wires or raw digits, an exception is thrown.

Definition at line 842 of file HitCreator.h.

Constructor & Destructor Documentation

recob::HitRefinerAssociator::HitRefinerAssociator	(	art::Event &	event,
		art::InputTag const &	HitModuleLabel,
		std::string	instance_name = "",
		bool	doWireAssns = true,
		bool	doRawDigitAssns = true
	)

Constructor: sets instance name and whether to build associations.

Parameters

event	the event the products are going to be put into
HitModuleLabel	label of the module used to create hits
instance_name	name of the instance for all data products
doWireAssns	whether to enable associations to wires
doRawDigitAssns	whether to enable associations to raw digits

All the data products (hit collection and associations) will have the specified product instance name.

Definition at line 513 of file HitCreator.cxx.

    : HitAndAssociationsWriterBase
        (event, instance_name, doWireAssns, doRawDigitAssns)
    , hits_label(HitModuleLabel)
  {
    hits.reset(new std::vector<recob::Hit>);
  } // HitRefinerAssociator::HitRefinerAssociator()

recob::HitRefinerAssociator::HitRefinerAssociator ( art::Event &  event, art::InputTag const &  HitModuleLabel, bool  doWireAssns, bool  doRawDigitAssns = true  )

inline

Constructor: sets instance name and whether to build associations.

Parameters

event	the event the products are going to be put into
HitModuleLabel	label of the module used to create hits
doWireAssns	whether to enable associations to wires
doRawDigitAssns	whether to enable associations to raw digits

All the data products (hit collection and associations) will have an empty product instance name.

Definition at line 874 of file HitCreator.h.

       :
      HitRefinerAssociator
        (event, HitModuleLabel, "", doWireAssns, doRawDigitAssns)
      {}

Member Function Documentation

void recob::HitRefinerAssociator::prepare_associations ( std::vector< recob::Hit > const &  srchits )

protected

Finds out the associations for the specified hits.

Definition at line 543 of file HitCreator.cxx.

  {
    if (!RawDigitAssns && !WireAssns) return; // no associations needed
    assert(event);

    // we make the associations anew
    if (RawDigitAssns) ClearAssociations(*RawDigitAssns);

    // read the hits; this is going to hurt performances...
    // no solution to that until there is a way to have a lazy read
    art::ValidHandle<std::vector<recob::Hit>> hHits
      = event->getValidHandle<std::vector<recob::Hit>>(hits_label);

    // now get the associations
    if (WireAssns) {
      // we make the associations anew
      ClearAssociations(*WireAssns);

      // find the associations between the hits and the wires
      art::FindOneP<recob::Wire> HitToWire(hHits, *event, hits_label);
      if (!HitToWire.isValid()) {
        throw art::Exception(art::errors::ProductNotFound)
          << "Can't find the associations between hits and wires produced by '"
          << hits_label << "'!\n";
      } // if no association

      // fill a map of wire vs. channel number
      std::vector<art::Ptr<recob::Wire>> WireMap;
      for (size_t iAssn = 0; iAssn < HitToWire.size(); ++iAssn) {
        art::Ptr<recob::Wire> wire = HitToWire.at(iAssn);
        if (wire.isNull()) continue;
        size_t channelID = (size_t) wire->Channel();
        if (WireMap.size() <= channelID) // expand the map of necessary
          WireMap.resize(std::max(channelID + 1, 2 * WireMap.size()), {});
        WireMap[channelID] = std::move(wire);
      } // for

      // now go through all the hits...
      for (size_t iHit = 0; iHit < srchits.size(); ++iHit) {
        recob::Hit const& hit = srchits[iHit];
        size_t channelID = (size_t) hit.Channel();

        // no association if there is no wire to associate with
        if ((channelID >= WireMap.size()) || !WireMap[channelID]) continue;

        // create an association using the same wire pointer
        WireAssns->addSingle(WireMap[channelID], CreatePtr(iHit));
      } // for hits
    } // if wire associations

    // now get the associations
    if (RawDigitAssns) {
      // we make the associations anew
      ClearAssociations(*RawDigitAssns);

      // find the associations between the hits and the raw digits
      art::FindOneP<raw::RawDigit> HitToDigits(hHits, *event, hits_label);
      if (!HitToDigits.isValid()) {
        throw art::Exception(art::errors::ProductNotFound)
          << "Can't find the associations between hits and raw digits"
          << " produced by '" << hits_label << "'!\n";
      } // if no association

      // fill a map of digits vs. channel number
      std::vector<art::Ptr<raw::RawDigit>> DigitMap;
      for (size_t iAssn = 0; iAssn < HitToDigits.size(); ++iAssn) {
        art::Ptr<raw::RawDigit> digits = HitToDigits.at(iAssn);
        if (digits.isNull()) continue;
        size_t channelID = (size_t) digits->Channel();
        if (DigitMap.size() <= channelID) // expand the map of necessary
          DigitMap.resize(std::max(channelID + 1, 2 * DigitMap.size()), {});
        DigitMap[channelID] = std::move(digits);
      } // for

      // now go through all the hits...
      for (size_t iHit = 0; iHit < srchits.size(); ++iHit) {
        recob::Hit const& hit = srchits[iHit];
        size_t channelID = (size_t) hit.Channel();

        // no association if there is no digits to associate with
        if ((channelID >= DigitMap.size()) || !DigitMap[channelID]) continue;

        // create an association using the same digits pointer
        RawDigitAssns->addSingle(DigitMap[channelID], CreatePtr(iHit));
      } // for hits
    } // if digit associations

  } // HitRefinerAssociator::put_into()

void recob::HitRefinerAssociator::prepare_associations ( )

inlineprotected

Finds out the associations for the current hits.

Definition at line 927 of file HitCreator.h.

{ prepare_associations(*hits); }

void recob::HitRefinerAssociator::put_into ( art::Event &  )

inline

Moves the data into the event.

The calling module must have already declared the production of these products with the proper instance name. After the move, the collections in this object are empty.

Deprecated:

Use the version with no arguments instead.

Definition at line 908 of file HitCreator.h.

{ put_into(); }

void recob::HitRefinerAssociator::put_into

(

)

Moves the data into the event.

The calling module must have already declared the production of these products with the proper instance name. After the move, the collections in this object are empty.

Definition at line 535 of file HitCreator.cxx.

                                      {
    prepare_associations();
    HitAndAssociationsWriterBase::put_into();
  } // HitRefinerAssociator::put_into()

void recob::HitRefinerAssociator::use_hits

(

std::unique_ptr< std::vector< recob::Hit >> &&

srchits

)

Uses the specified collection as data product.

Parameters

srchits

the collection to be used as data product

The very same collection is put into the event. This object will temporary own the collection until the hits are put into the event. If there were previous hits in the object, they are lost.

Definition at line 528 of file HitCreator.cxx.

  {
    hits = std::move(srchits);
  } // HitRefinerAssociator::use_hits()

Member Data Documentation

art::InputTag recob::HitRefinerAssociator::hits_label

protected

Label of the collection of hits.

Definition at line 921 of file HitCreator.h.

---

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

• HitCreator.h
• HitCreator.cxx