Inheritance diagram for hit::FFTHitFinder:

Public Types
using	ModuleType = EDProducer

using	WorkerType = WorkerT< EDProducer >

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

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

virtual	~FFTHitFinder ()

void	produce (art::Event &evt)

void	beginJob ()

void	endJob ()

void	reconfigure (fhicl::ParameterSet const &p)

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 Attributes
std::string	fCalDataModuleLabel

double	fMinSigInd
	Induction signal height threshold. More...

double	fMinSigCol
	Collection signal height threshold. More...

double	fIndWidth
	Initial width for induction fit. More...

double	fColWidth
	Initial width for collection fit. More...

double	fIndMinWidth
	Minimum induction hit width. More...

double	fColMinWidth
	Minimum collection hit width. More...

int	fMaxMultiHit
	maximum hits for multi fit More...

int	fAreaMethod
	Type of area calculation. More...

std::vector< double >	fAreaNorms
	factors for converting area to same units as peak height More...

Detailed Description

Definition at line 46 of file FFTHitFinder_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

hit::FFTHitFinder::FFTHitFinder ( fhicl::ParameterSet const & pset )

explicit

Definition at line 76 of file FFTHitFinder_module.cc.

References recob::HitAndAssociationsWriterBase::declare_products(), and reconfigure().

   {
     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(*this);
   }

hit::FFTHitFinder::~FFTHitFinder ( )

virtual

Definition at line 88 of file FFTHitFinder_module.cc.

89 {

90 }

Member Function Documentation

void hit::FFTHitFinder::beginJob ( )

virtual

Reimplemented from art::EDProducer.

Definition at line 108 of file FFTHitFinder_module.cc.

109 {

110 }

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 hit::FFTHitFinder::endJob ( )

virtual

Reimplemented from art::EDProducer.

Definition at line 113 of file FFTHitFinder_module.cc.

114 {

115 }

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 hit::FFTHitFinder::produce ( art::Event & evt )

virtual

Todo:

just get the integral from the fit for totSig

Todo:: need to have a disambiguation algorithm somewhere in here

Todo:

multiplicity and local_index have to be determined

Implements art::EDProducer.

Definition at line 121 of file FFTHitFinder_module.cc.

References recob::Wire::Channel(), geo::GeometryCore::ChannelToWire(), DEFINE_ART_MODULE, recob::HitCollectionCreator::emplace_back(), fAreaMethod, fAreaNorms, fCalDataModuleLabel, fColMinWidth, fColWidth, fIndMinWidth, fIndWidth, fMaxMultiHit, fMinSigCol, fMinSigInd, art::DataViewImpl::getByLabel(), raw::InvalidChannelID, geo::kCollection, geo::kInduction, recob::HitCollectionCreator::put_into(), recob::Wire::Signal(), and geo::GeometryCore::SignalType().

   { 
     
     // this object contains the hit collection
     // and its associations to wires and raw digits:
     recob::HitCollectionCreator hcol(*this, evt);
   
     // Read in the wire List object(s).
     art::Handle< std::vector<recob::Wire> > wireVecHandle;
     evt.getByLabel(fCalDataModuleLabel,wireVecHandle);
     art::ServiceHandle<geo::Geometry> geom;
    
     // also get the raw digits associated with wires
     art::FindOneP<raw::RawDigit> WireToRawDigits
       (wireVecHandle, evt, fCalDataModuleLabel);
     
     std::vector<int> startTimes;             // stores time of 1st local minimum
     std::vector<int> maxTimes;               // stores time of local maximum    
     std::vector<int> endTimes;               // stores time of 2nd local minimum
     int time               = 0;              // current time bin
     int minTimeHolder      = 0;              // current start time
     raw::ChannelID_t channel = raw::InvalidChannelID; // channel number
     bool maxFound          = false;          // Flag for whether a peak > threshold has been found
     double threshold       = 0.;             // minimum signal size for id'ing a hit
     double fitWidth        = 0.;             // hit fit width initial value
     double minWidth        = 0.;             // minimum hit width
     std::string eqn        = "gaus(0)";      // string for equation to fit
     geo::SigType_t sigType = geo::kInduction;// type of plane we are looking at
     std::stringstream numConv;
 
 
     //loop over wires
     for(unsigned int wireIter = 0; wireIter < wireVecHandle->size(); wireIter++) {
 
       art::Ptr<recob::Wire> wire(wireVecHandle, wireIter);
       startTimes.clear();
       maxTimes.clear();
       endTimes.clear();
       std::vector<float> signal(wire->Signal());
       std::vector<float>::iterator timeIter;   // iterator for time bins
       time          = 0;
       minTimeHolder = 0;
       maxFound      = false;
       channel       = wire->Channel();
       sigType       = geom->SignalType(channel);
 
       //Set the appropriate signal widths and thresholds
       if(sigType == geo::kInduction){
         threshold     = fMinSigInd;
         fitWidth      = fIndWidth;
         minWidth      = fIndMinWidth;
       }
       else if(sigType == geo::kCollection){
         threshold = fMinSigCol;
         fitWidth  = fColWidth;
         minWidth  = fColMinWidth;
       }
       // loop over signal
       for(timeIter = signal.begin(); timeIter+2 < signal.end(); timeIter++){    
         //test if timeIter+1 is a local minimum
         if(*timeIter > *(timeIter+1) && *(timeIter+1) < *(timeIter+2)){
           //only add points if already found a local max above threshold.
           if(maxFound) {
             endTimes.push_back(time+1);
             maxFound = false;
             //keep these in case new hit starts right away
             minTimeHolder = time+2; 
           }
           else minTimeHolder = time+1; 
         }
         //if not a minimum, test if we are at a local maximum 
         //if so, and the max value is above threshold, add it and proceed.
         else if(*timeIter < *(timeIter+1) && 
                 *(timeIter+1) > *(timeIter+2) && 
                 *(timeIter+1) > threshold){ 
           maxFound = true;
           maxTimes.push_back(time+1);
           startTimes.push_back(minTimeHolder);          
         }
         time++;
       }//end loop over signal vec
 
       
       //if no inflection found before end, but peak found add end point
       while(maxTimes.size()>endTimes.size()) 
         endTimes.push_back(signal.size()-1); 
       if(startTimes.size() == 0) continue;
       
       //All code below does the fitting, adding of hits
       //to the hit vector and when all wires are complete 
       //saving them 
       double totSig(0); // stores the total hit signal
       double startT(0); // stores the start time
       double endT(0);   // stores the end time
       int numHits(0);   // number of consecutive hits being fitted
       int size(0);      // size of data vector for fit
       int hitIndex(0);  // index of current hit in sequence
       double amplitude(0), position(0), width(0);  //fit parameters
       double amplitudeErr(0), positionErr(0), widthErr(0);  //fit errors
       double goodnessOfFit(0), chargeErr(0);  //Chi2/NDF and error on charge
       double minPeakHeight(0);  //lowest peak height in multi-hit fit
      
       //stores gaussian paramters first index is the hit number
       //the second refers to height, position, and width respectively
       std::vector<double>  hitSig;
 
       //add found hits to hit vector      
       while(hitIndex < (signed)startTimes.size()) {
         
         startT = endT = 0;
         numHits = 1;
         minPeakHeight = signal[maxTimes[hitIndex]];
 
         //consider adding pulse to group of consecutive hits if:
         //1 less than max consecutive hits
         //2 we are not at the last point in the signal vector
         //3 the height of the dip between the two is greater than threshold/2
         //4 and there is no gap between them
         while(numHits < fMaxMultiHit &&
               numHits+hitIndex < (signed)endTimes.size() && 
               signal[endTimes[hitIndex+numHits-1]] >threshold/2.0 &&  
               startTimes[hitIndex+numHits] - endTimes[hitIndex+numHits-1] < 2){
 
           if(signal[maxTimes[hitIndex+numHits]] < minPeakHeight) 
             minPeakHeight = signal[maxTimes[hitIndex+numHits]];
 
           ++numHits;
         }
 
         //finds the first point > 1/2 the smallest peak
         startT = startTimes[hitIndex];
 
         while(signal[(int)startT] < minPeakHeight/2.0) ++startT;
 
         //finds the first point from the end > 1/2 the smallest peak
         endT = endTimes[hitIndex+numHits-1];
 
         while(signal[(int)endT] <minPeakHeight/2.0) --endT;
         size = (int)(endT-startT);
         TH1D hitSignal("hitSignal","",size,startT,endT);
         for(int i = (int)startT; i < (int)endT; ++i)
           hitSignal.Fill(i,signal[i]);
 
         //build the TFormula
         eqn = "gaus(0)";        
 
         for(int i = 3; i < numHits*3; i+=3) {
           eqn.append("+gaus(");
           numConv.str("");
           numConv << i;
           eqn.append(numConv.str());
           eqn.append(")");
         }
 
         TF1 gSum("gSum",eqn.c_str(),0,size);
 
         if(numHits > 1) {
           TArrayD data(numHits*numHits);
           TVectorD amps(numHits); 
           for(int i = 0; i < numHits; ++i) {
             amps[i] = signal[maxTimes[hitIndex+i]];
             for(int j = 0; j < numHits;j++) 
               data[i+numHits*j] = TMath::Gaus(maxTimes[hitIndex+j],
                                               maxTimes[hitIndex+i],
                                               fitWidth);
           }//end loop over hits
           
           //This section uses a linear approximation in order to get an
           //initial value of the individual hit amplitudes 
           try{
             TMatrixD h(numHits,numHits);
             h.Use(numHits,numHits,data.GetArray());
             TDecompSVD a(h);
             a.Solve(amps);
           }
           catch(...){
             mf::LogInfo("FFTHitFinder")<<"TDcompSVD failed";
             hitIndex += numHits;
             continue;
           }
       
           for(int i = 0; i < numHits; ++i) {
             //if the approximation makes a peak vanish
             //set initial height as average of threshold and
             //raw peak height
             if(amps[i] > 0 ) amplitude = amps[i];
             else amplitude = 0.5*(threshold+signal[maxTimes[hitIndex+i]]);
             gSum.SetParameter(3*i,amplitude);
             gSum.SetParameter(1+3*i, maxTimes[hitIndex+i]);
             gSum.SetParameter(2+3*i, fitWidth);
             gSum.SetParLimits(3*i, 0.0, 3.0*amplitude);
             gSum.SetParLimits(1+3*i, startT , endT);
             gSum.SetParLimits(2+3*i, 0.0, 10.0*fitWidth);
           }//end loop over hits
         }//end if numHits > 1
         else {
           gSum.SetParameters(signal[maxTimes[hitIndex]],maxTimes[hitIndex],fitWidth);
           gSum.SetParLimits(0,0.0,1.5*signal[maxTimes[hitIndex]]);
           gSum.SetParLimits(1, startT , endT);
           gSum.SetParLimits(2,0.0,10.0*fitWidth);
         }
 
         hitSignal.Fit(&gSum,"QNRW","", startT, endT);
         for(int hitNumber = 0; hitNumber < numHits; ++hitNumber) {
           totSig = 0;
           if(gSum.GetParameter(3*hitNumber)   > threshold/2.0 && 
              gSum.GetParameter(3*hitNumber+2) > minWidth) { 
             amplitude     = gSum.GetParameter(3*hitNumber);
             position      = gSum.GetParameter(3*hitNumber+1);
             width         = gSum.GetParameter(3*hitNumber+2);
             amplitudeErr  = gSum.GetParError(3*hitNumber);
             positionErr   = gSum.GetParError(3*hitNumber+1);
             widthErr      = gSum.GetParError(3*hitNumber+2);
             goodnessOfFit = gSum.GetChisquare()/(double)gSum.GetNDF();
        int DoF = gSum.GetNDF();
 
             //estimate error from area of Gaussian
             chargeErr = std::sqrt(TMath::Pi())*(amplitudeErr*width+widthErr*amplitude);   
 
             hitSig.resize(size);
 
             for(int sigPos = 0; sigPos < size; ++sigPos){
               hitSig[sigPos] = amplitude*TMath::Gaus(sigPos+startT,position, width);
               totSig += hitSig[(int)sigPos];              
             }                       
 
             if(fAreaMethod) 
               totSig = std::sqrt(2*TMath::Pi())*amplitude*width/fAreaNorms[(size_t)sigType];              
 
             // get the WireID for this hit
             std::vector<geo::WireID> wids = geom->ChannelToWire(channel);
             // for now, just take the first option returned from ChannelToWire
             geo::WireID wid = wids[0];
 
             // make the hit
             recob::HitCreator hit(
               *wire,          // wire
               wid,            // wireID
               (int) startT,   // start_tick
               (int) endT,     // end_tick
               width,          // rms
               position,       // peak_time
               positionErr,    // sigma_peak_time
               amplitude,      // peak_amplitude
               amplitudeErr,   // sigma_peak_amplitude
               totSig,         // hit_integral
               chargeErr,      // hit_sigma_integral
               std::accumulate // summedADC
                 (signal.begin() + (int) startT, signal.begin() + (int) endT, 0.),
               1,              // multiplicity
               -1,             // local_index
               goodnessOfFit,  // goodness_of_fit
               DoF             // dof
               );
             
             // get the object associated with the original hit
             art::Ptr<raw::RawDigit> rawdigits = WireToRawDigits.at(wireIter);
             
             hcol.emplace_back(hit.move(), wire, rawdigits);
           }//end if over threshold
         }//end loop over hits
         hitIndex += numHits;    
       } // end while on hitIndex<(signed)startTimes.size()
 
     } // while on Wires
     
     // put the hit collection and associations into the event
     hcol.put_into(evt);
 
   } // End of produce()  

void hit::FFTHitFinder::reconfigure ( fhicl::ParameterSet const & p )

Definition at line 93 of file FFTHitFinder_module.cc.

References fAreaMethod, fAreaNorms, fCalDataModuleLabel, fColMinWidth, fColWidth, fIndMinWidth, fIndWidth, fMaxMultiHit, fMinSigCol, fMinSigInd, and fhicl::ParameterSet::get().

Referenced by FFTHitFinder().

   {
     fCalDataModuleLabel = p.get< std::string  >("CalDataModuleLabel");
     fMinSigInd          = p.get< double       >("MinSigInd");
     fMinSigCol          = p.get< double       >("MinSigCol"); 
     fIndWidth           = p.get< double       >("IndWidth");  
     fColWidth           = p.get< double       >("ColWidth");
     fIndMinWidth        = p.get< double       >("IndMinWidth");
     fColMinWidth        = p.get< double       >("ColMinWidth");                 
     fMaxMultiHit        = p.get< int          >("MaxMultiHit");
     fAreaMethod         = p.get< int          >("AreaMethod");
     fAreaNorms          = p.get< std::vector< double > >("AreaNorms");
   }  

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);
 }