hit::FFTHitFinder Class Reference

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.

  {
  }

Member Function Documentation

void hit::FFTHitFinder::beginJob ( )

virtual

Reimplemented from art::EDProducer.

Definition at line 108 of file FFTHitFinder_module.cc.

  { 
  }

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.

  {
  }

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.

Referenced by art::RootOutput::beginSubRun(), art::OutputModule::doBeginRun(), art::OutputModule::doBeginSubRun(), art::OutputModule::doEndRun(), art::OutputModule::doEndSubRun(), art::ProducingService::doPostReadEvent(), art::ProducingService::doPostReadRun(), art::ProducingService::doPostReadSubRun(), art::OutputModule::doWriteEvent(), art::ProcessPackage< L >::postScheduleSignal(), art::BeginEndPackage< Level::Run >::Begin::postScheduleSignal(), art::BeginEndPackage< Level::Run >::End::postScheduleSignal(), art::BeginEndPackage< Level::SubRun >::Begin::postScheduleSignal(), art::BeginEndPackage< Level::SubRun >::End::postScheduleSignal(), art::ProcessPackage< L >::preScheduleSignal(), art::BeginEndPackage< Level::Run >::Begin::preScheduleSignal(), art::BeginEndPackage< Level::SubRun >::Begin::preScheduleSignal(), art::EventProcessor::readEvent(), art::EventProcessor::readRun(), art::EmptyEvent::readRun_(), art::EventProcessor::readSubRun(), and art::EmptyEvent::readSubRun_().

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

Member Data Documentation

int hit::FFTHitFinder::fAreaMethod

private

Type of area calculation.

Definition at line 68 of file FFTHitFinder_module.cc.

Referenced by produce(), and reconfigure().

std::vector<double> hit::FFTHitFinder::fAreaNorms

private

factors for converting area to same units as peak height

Definition at line 69 of file FFTHitFinder_module.cc.

Referenced by produce(), and reconfigure().

std::string hit::FFTHitFinder::fCalDataModuleLabel

private

Definition at line 60 of file FFTHitFinder_module.cc.

Referenced by produce(), and reconfigure().

double hit::FFTHitFinder::fColMinWidth

private

Minimum collection hit width.

Definition at line 66 of file FFTHitFinder_module.cc.

Referenced by produce(), and reconfigure().

double hit::FFTHitFinder::fColWidth

private

Initial width for collection fit.

Definition at line 64 of file FFTHitFinder_module.cc.

Referenced by produce(), and reconfigure().

double hit::FFTHitFinder::fIndMinWidth

private

Minimum induction hit width.

Definition at line 65 of file FFTHitFinder_module.cc.

Referenced by produce(), and reconfigure().

double hit::FFTHitFinder::fIndWidth

private

Initial width for induction fit.

Definition at line 63 of file FFTHitFinder_module.cc.

Referenced by produce(), and reconfigure().

int hit::FFTHitFinder::fMaxMultiHit

private

maximum hits for multi fit

Definition at line 67 of file FFTHitFinder_module.cc.

Referenced by produce(), and reconfigure().

double hit::FFTHitFinder::fMinSigCol

private

Collection signal height threshold.

Definition at line 62 of file FFTHitFinder_module.cc.

Referenced by produce(), and reconfigure().

double hit::FFTHitFinder::fMinSigInd

private

Induction signal height threshold.

Definition at line 61 of file FFTHitFinder_module.cc.

Referenced by produce(), and reconfigure().

---

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

• /cvmfs/larsoft.opensciencegrid.org/products/larreco/v07_10_02/source/larreco/HitFinder/FFTHitFinder_module.cc