dune::HitFinder35t Class Reference

Inheritance diagram for dune::HitFinder35t:

Public Member Functions
	HitFinder35t (fhicl::ParameterSet const &pset)
void	produce (art::Event &evt)
void	beginJob ()
void	endJob ()
void	reconfigure (fhicl::ParameterSet const &p)
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 Attributes
DisambigAlg35t	fDisambigAlg
TimeBasedDisambig	fTimeBasedDisambigAlg
art::ServiceHandle< geo::Geometry >	fGeom
std::string	fChanHitLabel
std::string	fAlg

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 49 of file HitFinder35t_module.cc.

Constructor & Destructor Documentation

dune::HitFinder35t::HitFinder35t ( fhicl::ParameterSet const &  pset )

explicit

Definition at line 79 of file HitFinder35t_module.cc.

                                                          :
    EDProducer(pset),
    fDisambigAlg(pset.get< fhicl::ParameterSet >("DisambigAlg")),
    fTimeBasedDisambigAlg(pset.get< fhicl::ParameterSet >("TimeBasedDisambigAlg"))
  {
    this->reconfigure(pset);
    
    // let HitCollectionCreator declare that we are going to produce
    // hits and associations with wires and raw digits
    // (with no particular product label)
    recob::HitCollectionCreator::declare_products(producesCollector(), "", true, false);
  }

Member Function Documentation

void dune::HitFinder35t::beginJob ( )

virtual

Reimplemented from art::EDProducer.

Definition at line 105 of file HitFinder35t_module.cc.

  {
    
  }

void dune::HitFinder35t::endJob ( )

virtual

Reimplemented from art::EDProducer.

Definition at line 112 of file HitFinder35t_module.cc.

  {
    
  }

void dune::HitFinder35t::produce ( art::Event &  evt )

virtual

Implements art::EDProducer.

Definition at line 119 of file HitFinder35t_module.cc.

  {
    
    auto ChannelHits = evt.getHandle< std::vector<recob::Hit> >(fChanHitLabel);
    
    // also get the associated wires and raw digits;
    // we assume they have been created by the same module as the hits
//    art::FindOneP<raw::RawDigit> ChannelHitRawDigits
//      (ChannelHits, evt, fChanHitLabel);
    art::FindOneP<recob::Wire> ChannelHitWires
      (ChannelHits, evt, fChanHitLabel);

    // this object contains the hit collection
    // and its associations to wires and raw digits:
    recob::HitCollectionCreator hcol(evt,
      /* doWireAssns */ ChannelHitWires.isValid(),
      /* doRawDigitAssns */ false
      );
    
    // Make unambiguous collection hits
    std::vector< art::Ptr<recob::Hit> >  ChHits;
    art::fill_ptr_vector(ChHits, ChannelHits);
    for( size_t h = 0; h < ChHits.size(); h++ ){
      if( ChHits[h]->View() != geo::kZ ) continue;
      
      art::Ptr<recob::Wire> wire = ChannelHitWires.at(h);
      //art::Ptr<raw::RawDigit> rawdigits = ChannelHitRawDigits.at(h);
    
      // just copy it
      hcol.emplace_back(*ChHits[h], wire);
      
    } // for
    
    // Run alg on all APAs
    if (fAlg == "TripletMatch")
      {
        auto const clockData = art::ServiceHandle<detinfo::DetectorClocksService const>()->DataFor(evt);
        auto const detProp = art::ServiceHandle<detinfo::DetectorPropertiesService const>()->DataFor(evt, clockData);
        fDisambigAlg.RunDisambig(clockData, detProp, ChHits);

        for( size_t t=0; t < fDisambigAlg.fDisambigHits.size(); t++ ){
          art::Ptr<recob::Hit>  hit = fDisambigAlg.fDisambigHits[t].first;
          geo::WireID           wid = fDisambigAlg.fDisambigHits[t].second;
      
          // create a new hit copy of the original one, but with new wire ID
          recob::HitCreator disambiguous_hit(*hit, wid);
      
          // get the objects associated with the original hit;
          // since hit comes from ChannelHits, its key is the index in that collection
          // and also the index for the query of associated objects
          art::Ptr<recob::Hit>::key_type hit_index = hit.key();
          art::Ptr<recob::Wire> wire = ChannelHitWires.at(hit_index);
          //art::Ptr<raw::RawDigit> rawdigits = ChannelHitRawDigits.at(hit_index);
      
          hcol.emplace_back(disambiguous_hit.move(), wire);
        } // for

      }
    else if (fAlg == "TimeBased")
      {
        fTimeBasedDisambigAlg.RunDisambig(ChHits);

        for( size_t t=0; t < fTimeBasedDisambigAlg.fDisambigHits.size(); t++ ){
          art::Ptr<recob::Hit>  hit = fTimeBasedDisambigAlg.fDisambigHits[t].first;
          geo::WireID           wid = fTimeBasedDisambigAlg.fDisambigHits[t].second;
      
          // create a new hit copy of the original one, but with new wire ID
          recob::HitCreator disambiguous_hit(*hit, wid);
      
          // get the objects associated with the original hit;
          // since hit comes from ChannelHits, its key is the index in that collection
          // and also the index for the query of associated objects
          art::Ptr<recob::Hit>::key_type hit_index = hit.key();
          art::Ptr<recob::Wire> wire = ChannelHitWires.at(hit_index);
          //art::Ptr<raw::RawDigit> rawdigits = ChannelHitRawDigits.at(hit_index);
      
          hcol.emplace_back(disambiguous_hit.move(), wire);
        } // for

      }
    else
      {
        throw cet::exception("HitFinder35t") << "Disambiguation algorithm name: " << fAlg << " is not supported.\n" ;
      }
    
    
    // put the hit collection and associations into the event
    hcol.put_into(evt);
    
  }

void dune::HitFinder35t::reconfigure

(

fhicl::ParameterSet const &

p

)

Definition at line 95 of file HitFinder35t_module.cc.

  {
    
    fChanHitLabel =  p.get< std::string >("ChanHitLabel");
    fAlg = p.get < std::string >("Algorithm");
    
  }  

Member Data Documentation

std::string dune::HitFinder35t::fAlg

private

Definition at line 69 of file HitFinder35t_module.cc.

std::string dune::HitFinder35t::fChanHitLabel

private

Definition at line 68 of file HitFinder35t_module.cc.

DisambigAlg35t dune::HitFinder35t::fDisambigAlg

private

Definition at line 63 of file HitFinder35t_module.cc.

art::ServiceHandle<geo::Geometry> dune::HitFinder35t::fGeom

private

Definition at line 66 of file HitFinder35t_module.cc.

TimeBasedDisambig dune::HitFinder35t::fTimeBasedDisambigAlg

private

Definition at line 64 of file HitFinder35t_module.cc.

---

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

• dunereco/dunereco/HitFinderDUNE/HitFinder35t_module.cc