Inheritance diagram for CRHitRemovalByPCA:

Public Types
using	ModuleType = EDProducer

using	WorkerType = WorkerT< EDProducer >

template<typename UserConfig , typename KeysToIgnore = void>
using	Table = ProducerBase::Table< UserConfig, KeysToIgnore >

Public Member Functions
	CRHitRemovalByPCA (fhicl::ParameterSet const &pset)

virtual	~CRHitRemovalByPCA ()
	Destructor. More...

virtual void	reconfigure (fhicl::ParameterSet const &pset)

virtual void	produce (art::Event &e)

virtual void	beginJob ()
	Begin job method. More...

virtual void	endJob ()
	End job method. More...

template<typename PROD , BranchType B = InEvent>
ProductID	getProductID (std::string const &instanceName={}) const

template<typename PROD , BranchType B>
ProductID	getProductID (ModuleDescription const &moduleDescription, std::string const &instanceName) const

bool	modifiesEvent () const

template<typename T , BranchType = InEvent>
ProductToken< T >	consumes (InputTag const &)

template<typename T , art::BranchType BT>
art::ProductToken< T >	consumes (InputTag const &it)

template<typename T , BranchType = InEvent>
void	consumesMany ()

template<typename Element , BranchType = InEvent>
ViewToken< Element >	consumesView (InputTag const &)

template<typename T , art::BranchType BT>
art::ViewToken< T >	consumesView (InputTag const &it)

template<typename T , BranchType = InEvent>
ProductToken< T >	mayConsume (InputTag const &)

template<typename T , art::BranchType BT>
art::ProductToken< T >	mayConsume (InputTag const &it)

template<typename T , BranchType = InEvent>
void	mayConsumeMany ()

template<typename Element , BranchType = InEvent>
ViewToken< Element >	mayConsumeView (InputTag const &)

template<typename T , art::BranchType BT>
art::ViewToken< T >	mayConsumeView (InputTag const &it)

base_engine_t &	createEngine (seed_t seed)

base_engine_t &	createEngine (seed_t seed, std::string const &kind_of_engine_to_make)

base_engine_t &	createEngine (seed_t seed, std::string const &kind_of_engine_to_make, label_t const &engine_label)

seed_t	get_seed_value (fhicl::ParameterSet const &pset, char const key[]="seed", seed_t const implicit_seed=-1)

Static Public Member Functions
static cet::exempt_ptr< Consumer >	non_module_context ()

Protected Member Functions
CurrentProcessingContext const *	currentContext () const

void	validateConsumedProduct (BranchType const bt, ProductInfo const &pi)

void	prepareForJob (fhicl::ParameterSet const &pset)

void	showMissingConsumes () const

Private Member Functions
void	removeTaggedHits (const recob::PFParticle pfParticle, const art::Handle< std::vector< recob::PFParticle > > &pfParticleHandle, const art::FindManyP< recob::Cluster > &partToClusAssns, const art::FindManyP< recob::Hit > &clusToHitAssns, std::set< const recob::PFParticle > &taggedParticles, art::PtrVector< recob::Hit > &hitVec)

Private Attributes
std::string	fCosmicProducerLabel
	Module that produced the PCA based cosmic tags. More...

std::string	fHitProducerLabel
	The full collection of hits. More...

std::string	fPFParticleProducerLabel
	PFParticle producer. More...

double	fCosmicTagThreshold
	Thresholds for tagging. More...

int	fNumEvent
	Number of events seen. More...

int	fNumCRRejects
	Number of tracks produced. More...

Detailed Description

Definition at line 70 of file CRHitRemovalByPCA_module.cc.

Member Typedef Documentation

using art::EDProducer::ModuleType = EDProducer

inherited

Definition at line 34 of file EDProducer.h.

template<typename UserConfig , typename KeysToIgnore = void>

using art::EDProducer::Table = ProducerBase::Table<UserConfig, KeysToIgnore>

inherited

Definition at line 43 of file EDProducer.h.

using art::EDProducer::WorkerType = WorkerT<EDProducer>

inherited

Definition at line 35 of file EDProducer.h.

Constructor & Destructor Documentation

CRHitRemovalByPCA::CRHitRemovalByPCA ( fhicl::ParameterSet const & pset )

explicit

Constructor.

Arguments:

pset - Fcl parameters.

Definition at line 114 of file CRHitRemovalByPCA_module.cc.

References reconfigure().

                                                                    :
   fNumEvent(0),
   fNumCRRejects(0)
 {
     reconfigure(pset);
     produces<std::vector<recob::Hit> >();
 
     // Report.
     mf::LogInfo("CRHitRemovalByPCA") << "CRHitRemovalByPCA configured\n";
 }

CRHitRemovalByPCA::~CRHitRemovalByPCA ( )

virtual

Destructor.

Definition at line 127 of file CRHitRemovalByPCA_module.cc.

128 {}

Member Function Documentation

void CRHitRemovalByPCA::beginJob ( )

virtual

Begin job method.

Reimplemented from art::EDProducer.

Definition at line 147 of file CRHitRemovalByPCA_module.cc.

148 {

149 }

template<typename T , BranchType = InEvent>

ProductToken<T> art::Consumer::consumes ( InputTag const & )

inherited

template<typename T , art::BranchType BT>

art::ProductToken<T> art::Consumer::consumes ( InputTag const & it )

inherited

Definition at line 147 of file Consumer.h.

References art::InputTag::instance(), art::InputTag::label(), and art::InputTag::process().

 {
   if (!moduleContext_)
     return ProductToken<T>::invalid();
 
   consumables_[BT].emplace_back(ConsumableType::Product,
                                 TypeID{typeid(T)},
                                 it.label(),
                                 it.instance(),
                                 it.process());
   return ProductToken<T>{it};
 }

template<typename T , art::BranchType BT>

void art::Consumer::consumesMany ( )

inherited

Definition at line 162 of file Consumer.h.

 {
   if (!moduleContext_)
     return;
 
   consumables_[BT].emplace_back(ConsumableType::Many, TypeID{typeid(T)});
 }

template<typename Element , BranchType = InEvent>

ViewToken<Element> art::Consumer::consumesView ( InputTag const & )

inherited

template<typename T , art::BranchType BT>

art::ViewToken<T> art::Consumer::consumesView ( InputTag const & it )

inherited

Definition at line 172 of file Consumer.h.

References art::InputTag::instance(), art::InputTag::label(), and art::InputTag::process().

 {
   if (!moduleContext_)
     return ViewToken<T>::invalid();
 
   consumables_[BT].emplace_back(ConsumableType::ViewElement,
                                 TypeID{typeid(T)},
                                 it.label(),
                                 it.instance(),
                                 it.process());
   return ViewToken<T>{it};
 }

EngineCreator::base_engine_t & EngineCreator::createEngine ( seed_t seed )

inherited

Definition at line 26 of file EngineCreator.cc.

References art::EngineCreator::rng().

Referenced by evgen::CosmicsGen::CosmicsGen(), rndm::NuRandomService::createEngine(), cluster::fuzzyCluster::fuzzyCluster(), cluster::HoughLineFinder::HoughLineFinder(), art::MixFilter< T >::initEngine_(), larg4::LArG4::LArG4(), evgen::LightSource::LightSource(), evgen::NeutronOsc::NeutronOsc(), evgen::NucleonDecay::NucleonDecay(), opdet::OpMCDigi::OpMCDigi(), opdet::OptDetDigitizer::OptDetDigitizer(), phot::PhotonLibraryPropagation::PhotonLibraryPropagation(), detsim::SimDriftElectrons::SimDriftElectrons(), evgen::SingleGen::SingleGen(), evgen::SNNueAr40CCGen::SNNueAr40CCGen(), ToyOneShowerGen::ToyOneShowerGen(), and trkf::Track3DKalman::Track3DKalman().

 {
   return rng()->createEngine(placeholder_schedule_id(), seed);
 }

EngineCreator::base_engine_t & EngineCreator::createEngine	(	seed_t	seed,
		std::string const &	kind_of_engine_to_make
	)

inherited

Definition at line 32 of file EngineCreator.cc.

References art::EngineCreator::rng().

 {
   return rng()->createEngine(
     placeholder_schedule_id(), seed, kind_of_engine_to_make);
 }

EngineCreator::base_engine_t & EngineCreator::createEngine	(	seed_t	seed,
		std::string const &	kind_of_engine_to_make,
		label_t const &	engine_label
	)

inherited

Definition at line 40 of file EngineCreator.cc.

References art::EngineCreator::rng().

 {
   return rng()->createEngine(
     placeholder_schedule_id(), seed, kind_of_engine_to_make, engine_label);
 }

CurrentProcessingContext const * art::EDProducer::currentContext ( ) const

protectedinherited

Definition at line 120 of file EDProducer.cc.

References art::EDProducer::current_context_.

   {
     return current_context_.get();
   }

void CRHitRemovalByPCA::endJob ( )

virtual

End job method.

Reimplemented from art::EDProducer.

Definition at line 370 of file CRHitRemovalByPCA_module.cc.

References fNumCRRejects, and fNumEvent.

 {
     double aveCRPerEvent = fNumEvent > 0 ? double(fNumCRRejects) / double(fNumEvent) : 0.;
     
     mf::LogInfo("CRHitRemovalByPCA")
         << "CRHitRemovalByPCA statistics:\n"
         << "  Number of events = " << fNumEvent << "\n"
         << "  Number of Cosmic Rays found = " << fNumCRRejects
         << ", " << aveCRPerEvent << " average/event";
 }

EngineCreator::seed_t EngineCreator::get_seed_value	(	fhicl::ParameterSet const &	pset,
		char const	key[] = `"seed"`,
		seed_t const	implicit_seed = `-1`
	)

inherited

Definition at line 49 of file EngineCreator.cc.

References fhicl::ParameterSet::get().

Referenced by art::MixFilter< T >::initEngine_().

 {
   auto const& explicit_seeds = pset.get<std::vector<int>>(key, {});
   return explicit_seeds.empty() ? implicit_seed : explicit_seeds.front();
 }

template<typename PROD , BranchType B>

ProductID art::EDProducer::getProductID ( std::string const & instanceName = {} ) const

inlineinherited

Definition at line 123 of file EDProducer.h.

References art::EDProducer::moduleDescription_.

   {
     return ProducerBase::getProductID<PROD, B>(moduleDescription_,
                                                instanceName);
   }

template<typename PROD , BranchType B>

ProductID art::ProducerBase::getProductID	(	ModuleDescription const &	moduleDescription,
		std::string const &	instanceName
	)		const

inherited

Definition at line 56 of file ProducerBase.h.

References B, and art::ModuleDescription::moduleLabel().

Referenced by art::ProducerBase::modifiesEvent().

   {
     auto const& pd =
       get_ProductDescription<PROD>(B, md.moduleLabel(), instanceName);
     return pd.productID();
   }

template<typename T , BranchType = InEvent>

ProductToken<T> art::Consumer::mayConsume ( InputTag const & )

inherited

template<typename T , art::BranchType BT>

art::ProductToken<T> art::Consumer::mayConsume ( InputTag const & it )

inherited

Definition at line 190 of file Consumer.h.

References art::InputTag::instance(), art::InputTag::label(), and art::InputTag::process().

 {
   if (!moduleContext_)
     return ProductToken<T>::invalid();
 
   consumables_[BT].emplace_back(ConsumableType::Product,
                                 TypeID{typeid(T)},
                                 it.label(),
                                 it.instance(),
                                 it.process());
   return ProductToken<T>{it};
 }

template<typename T , art::BranchType BT>

void art::Consumer::mayConsumeMany ( )

inherited

Definition at line 205 of file Consumer.h.

 {
   if (!moduleContext_)
     return;
 
   consumables_[BT].emplace_back(ConsumableType::Many, TypeID{typeid(T)});
 }

template<typename Element , BranchType = InEvent>

ViewToken<Element> art::Consumer::mayConsumeView ( InputTag const & )

inherited

template<typename T , art::BranchType BT>

art::ViewToken<T> art::Consumer::mayConsumeView ( InputTag const & it )

inherited

Definition at line 215 of file Consumer.h.

References art::InputTag::instance(), art::InputTag::label(), and art::InputTag::process().

 {
   if (!moduleContext_)
     return ViewToken<T>::invalid();
 
   consumables_[BT].emplace_back(ConsumableType::ViewElement,
                                 TypeID{typeid(T)},
                                 it.label(),
                                 it.instance(),
                                 it.process());
   return ViewToken<T>{it};
 }

bool art::ProducerBase::modifiesEvent ( ) const

inlineinherited

Definition at line 40 of file ProducerBase.h.

References art::ProducerBase::getProductID().

     {
       return true;
     }

cet::exempt_ptr< art::Consumer > art::Consumer::non_module_context ( )

staticinherited

Definition at line 76 of file Consumer.cc.

 {
   static Consumer invalid{InvalidTag{}};
   return &invalid;
 }

void art::Consumer::prepareForJob ( fhicl::ParameterSet const & pset )

protectedinherited

Definition at line 89 of file Consumer.cc.

References fhicl::ParameterSet::get_if_present().

Referenced by art::EDProducer::doBeginJob(), art::EDFilter::doBeginJob(), and art::EDAnalyzer::doBeginJob().

 {
   if (!moduleContext_)
     return;
 
   pset.get_if_present("errorOnMissingConsumes", requireConsumes_);
   for (auto& consumablesPerBranch : consumables_) {
     cet::sort_all(consumablesPerBranch);
   }
 }

void CRHitRemovalByPCA::produce ( art::Event & evt )

virtual

Produce method.

Arguments:

evt - Art event.

This is the primary method. The goal is to produce a list of recob::Hit objects which are a "clean" subset of all hits and which are believed to be due to a neutrino interaction. It does this by considering input CosmicTag objects, relating them to PFParticles/Tracks and removing the hits associated to those objects which are believed to be Cosmic Rays.

Implements art::EDProducer.

Definition at line 164 of file CRHitRemovalByPCA_module.cc.

References anab::CosmicTag::CosmicScore(), art::PtrVector< T >::empty(), art::PtrVector< T >::end(), fCosmicProducerLabel, fCosmicTagThreshold, fHitProducerLabel, fNumEvent, fPFParticleProducerLabel, art::Ptr< T >::get(), art::DataViewImpl::getByLabel(), art::PtrVector< T >::insert(), recob::PFParticle::IsPrimary(), art::Handle< T >::isValid(), art::PtrVector< T >::push_back(), art::Event::put(), and removeTaggedHits().

 {
     ++fNumEvent;
     
     // Start by looking up the original hits
     art::Handle< std::vector<recob::Hit> > hitHandle;
     evt.getByLabel(fHitProducerLabel, hitHandle);
     
     // If there are no hits then there should be no output
     if (!hitHandle.isValid()) return;
     
     // If there are hits then we are going to output something so get a new
     // output hit vector
     std::unique_ptr<std::vector<recob::Hit> > outputHits(new std::vector<recob::Hit>);
     
     // And fill it with the complete original list of hits
     *outputHits = *hitHandle;
     
     // Recover the PFParticles that are responsible for making the tracks
     art::Handle<std::vector<recob::PFParticle> > pfParticleHandle;
     evt.getByLabel(fPFParticleProducerLabel, pfParticleHandle);
     
     // Without a valid collection of PFParticles we can't do the hit removal
     if (!pfParticleHandle.isValid())
     {
         evt.put(std::move(outputHits));
         return;
     }
     
     // Recover the clusters so we can do associations to the hits
     // In theory the clusters come from the same producer as the PFParticles
     art::Handle<std::vector<recob::Cluster> > clusterHandle;
     evt.getByLabel(fPFParticleProducerLabel, clusterHandle);
     
     // If there are no clusters then something is really wrong
     if (!clusterHandle.isValid())
     {
         evt.put(std::move(outputHits));
         return;
     }
     
     // Recover the list of cosmic tags
     art::Handle< std::vector<anab::CosmicTag> > cosmicTagHandle;
     evt.getByLabel(fCosmicProducerLabel, cosmicTagHandle);
     
     // No cosmic tags then nothing to do here
     if (!cosmicTagHandle.isValid() || cosmicTagHandle->empty())
     {
         evt.put(std::move(outputHits));
         return;
     }
     
     // Start recovering the necessary associations
     // Start by finding the associations going from cosmic tags to PFParticles
     art::FindManyP<recob::PFParticle> cosmicTagToPFPartAssns(cosmicTagHandle, evt, fCosmicProducerLabel);
     
     // From PFParticles we go to clusters
     art::FindManyP<recob::Cluster> clusterAssns(pfParticleHandle, evt, fPFParticleProducerLabel);
     
     // Likewise, recover the collection of associations to hits
     art::FindManyP<recob::Hit> clusterHitAssns(clusterHandle, evt, fPFParticleProducerLabel);
     
     // Container to contain the "bad" hits...
     art::PtrVector<recob::Hit> taggedHits;
     
     // No point double counting hits
     std::set<const recob::PFParticle*> taggedSet;
 
     // Start the identification of hits to remove. The outer loop is over the various producers of
     // the CosmicTag objects we're examininig
     for(size_t crIdx = 0; crIdx != cosmicTagHandle->size(); crIdx++)
     {
         art::Ptr<anab::CosmicTag> cosmicTag(cosmicTagHandle, crIdx);
         
             // If this was tagged as a CR muon then we have work to do!
         if (cosmicTag->CosmicScore() > fCosmicTagThreshold)
         {
             // Recover the associated PFParticle
             std::vector<art::Ptr<recob::PFParticle> > pfPartVec = cosmicTagToPFPartAssns.at(crIdx);
             
             if (pfPartVec.empty()) continue;
             
             art::Ptr<recob::PFParticle> pfParticle = pfPartVec.front();
             
             // Again, most likely needless
             if (!pfParticle) continue;
             
             // A cosmic ray must be a primary (by fiat)
             if (!pfParticle->IsPrimary()) continue;
             
             // Avoid double counting if more than one tagger running
             if (taggedSet.find(pfParticle.get()) != taggedSet.end()) continue;
             
             // Remove all hits associated to this particle and its daughters
             removeTaggedHits(pfParticle.get(), pfParticleHandle, clusterAssns, clusterHitAssns, taggedSet, taggedHits);
         }
     }
     
     // Are there any tagged hits?
     if (!taggedHits.empty())
     {
         // First order of business is to attempt to restore any hits which are shared between a tagged
         // CR PFParticle and an untagged one. We can do this by going through the PFParticles and
         // "removing" hits which are in the not tagged set.
         art::PtrVector<recob::Hit> untaggedHits;
         
         for(const auto& pfParticle : *pfParticleHandle)
         {
             if (taggedSet.find(&pfParticle) != taggedSet.end()) continue;
             
             // Recover the clusters associated to the input PFParticle
             std::vector<art::Ptr<recob::Cluster> > clusterVec = clusterAssns.at(pfParticle.Self());
             
             // Loop over the clusters and grab the associated hits
             for(const auto& cluster : clusterVec)
             {
                 std::vector<art::Ptr<recob::Hit> > clusHitVec = clusterHitAssns.at(cluster->ID());
                 untaggedHits.insert(untaggedHits.end(), clusHitVec.begin(), clusHitVec.end());
             }
         }
         
         // Filter out the hits we want to save
         FilterHits(taggedHits, untaggedHits);
         
         // The below is rather ugly but there is an interplay between art::Ptr's and the
         // actual pointers to objects that I might be missing and this is what I see how to do
         // First move all the original art::Ptr hits into a local art::PtrVector
         art::PtrVector<recob::Hit> originalHits;
 
         // Fill this one hit at a time...
         for(size_t hitIdx = 0; hitIdx != hitHandle->size(); hitIdx++)
         {
             art::Ptr<recob::Hit> hit(hitHandle, hitIdx);
             
             originalHits.push_back(hit);
         }
         
         // Remove the cosmic ray tagged hits
         FilterHits(originalHits, taggedHits);
         
         // Clear the current outputHits vector since we're going to refill...
         outputHits->clear();
         
         // Now make the new list of output hits
         for (const auto& hit : originalHits)
         {
             // Kludge to remove out of time hits
             if (hit->StartTick() > 6400 || hit->EndTick() < 3200) continue;
             
             outputHits->emplace_back(*hit);
         }
     }
     
     // Add tracks and associations to event.
     evt.put(std::move(outputHits));
 }

void CRHitRemovalByPCA::reconfigure ( fhicl::ParameterSet const & pset )

virtual

Reconfigure method.

Arguments:

pset - Fcl parameter set.

Definition at line 137 of file CRHitRemovalByPCA_module.cc.

References fCosmicProducerLabel, fCosmicTagThreshold, fHitProducerLabel, fPFParticleProducerLabel, and fhicl::ParameterSet::get().

Referenced by CRHitRemovalByPCA().

 {
     fCosmicProducerLabel     = pset.get<std::string>("CosmicProducerLabel");
     fHitProducerLabel        = pset.get<std::string>("HitProducerLabel");
     fPFParticleProducerLabel = pset.get<std::string>("PFParticleProducerLabel");
     fCosmicTagThreshold      = pset.get<double>     ("CosmicTagThreshold");
 }

void CRHitRemovalByPCA::removeTaggedHits	(	const recob::PFParticle *	pfParticle,
		const art::Handle< std::vector< recob::PFParticle > > &	pfParticleHandle,
		const art::FindManyP< recob::Cluster > &	partToClusAssns,
		const art::FindManyP< recob::Hit > &	clusToHitAssns,
		std::set< const recob::PFParticle * > &	taggedParticles,
		art::PtrVector< recob::Hit > &	hitVec
	)

private

Hit removal method

Arguments:

pfParticle - the top level PFParticle to have hits removed pfParticleHandle - handle to the PFParticle objects partToClusAssns - list of PFParticle to Cluster associations clusToHitAssns - list of Cluster to Hit associations hitVec - the current list of hits

This recursively called method will remove all hits associated to an input PFParticle and, in addition, will call itself for all daughters of the input PFParticle

Definition at line 336 of file CRHitRemovalByPCA_module.cc.

References recob::PFParticle::Daughters(), art::PtrVector< T >::end(), art::Ptr< T >::get(), art::PtrVector< T >::insert(), and recob::PFParticle::Self().

Referenced by produce().

 {
     // Recover the clusters associated to the input PFParticle
     std::vector<art::Ptr<recob::Cluster> > clusterVec = partToClusAssns.at(pfParticle->Self());
     
     // Record this PFParticle as tagged
     taggedParticles.insert(pfParticle);
     
     // Loop over the clusters and grab the associated hits
     for(const auto& cluster : clusterVec)
     {
         std::vector<art::Ptr<recob::Hit> > clusHitVec = clusToHitAssns.at(cluster->ID());
         hitVec.insert(hitVec.end(), clusHitVec.begin(), clusHitVec.end());
     }
     
     // Loop over the daughters of this particle and remove their hits as well
     for(const auto& daughterId : pfParticle->Daughters())
     {
         art::Ptr<recob::PFParticle> daughter(pfParticleHandle, daughterId);
         
         removeTaggedHits(daughter.get(), pfParticleHandle, partToClusAssns, clusToHitAssns, taggedParticles, hitVec);
     }
     
     return;
 }

void art::Consumer::showMissingConsumes ( ) const

protectedinherited

Definition at line 125 of file Consumer.cc.

Referenced by art::EDProducer::doEndJob(), art::EDFilter::doEndJob(), art::EDAnalyzer::doEndJob(), and art::RootOutput::endJob().

 {
   if (!moduleContext_)
     return;
 
   // If none of the branches have missing consumes statements, exit early.
   if (std::all_of(cbegin(missingConsumes_),
                   cend(missingConsumes_),
                   [](auto const& perBranch) { return perBranch.empty(); }))
     return;
 
   constexpr cet::HorizontalRule rule{60};
   mf::LogPrint log{"MTdiagnostics"};
   log << '\n'
       << rule('=') << '\n'
       << "The following consumes (or mayConsume) statements are missing from\n"
       << module_context(moduleDescription_) << '\n'
       << rule('-') << '\n';
 
   cet::for_all_with_index(
     missingConsumes_, [&log](std::size_t const i, auto const& perBranch) {
       for (auto const& pi : perBranch) {
         log << "  "
             << assemble_consumes_statement(static_cast<BranchType>(i), pi)
             << '\n';
       }
     });
   log << rule('=');
 }

void art::Consumer::validateConsumedProduct	(	BranchType const	bt,
		ProductInfo const &	pi
	)

protectedinherited

Definition at line 101 of file Consumer.cc.

References art::errors::ProductRegistrationFailure.

 {
   // Early exits if consumes tracking has been disabled or if the
   // consumed product is an allowed consumable.
   if (!moduleContext_)
     return;
 
   if (cet::binary_search_all(consumables_[bt], pi))
     return;
 
   if (requireConsumes_) {
     throw Exception(errors::ProductRegistrationFailure,
                     "Consumer: an error occurred during validation of a "
                     "retrieved product\n\n")
       << "The following consumes (or mayConsume) statement is missing from\n"
       << module_context(moduleDescription_) << ":\n\n"
       << "  " << assemble_consumes_statement(bt, pi) << "\n\n";
   }
 
   missingConsumes_[bt].insert(pi);
 }