Inheritance diagram for hit::RawHitFinder:

Public Types
using	ModuleType = EDProducer

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

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

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)

void	fillProductDescriptions ()

void	registerProducts (ProductDescriptions &productsToRegister)

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)

std::unique_ptr< Worker >	makeWorker (WorkerParams const &wp)

template<typename T , BranchType BT>
ViewToken< T >	consumesView (InputTag const &tag)

template<typename T , BranchType BT>
ViewToken< T >	mayConsumeView (InputTag const &tag)

Protected Member Functions
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 ()

Private Member Functions
void	produce (art::Event &evt) override

Private Attributes
unsigned int	fDataSize

art::InputTag	fDigitModuleLabel

std::string	fSpillName

std::string	fCalDataModuleLabel

std::string	fHitLabelName

double	fMinSigInd

double	fMinSigCol

double	fIndWidth

double	fColWidth

double	fIndMinWidth

double	fColMinWidth

int	fMaxMultiHit

int	fAreaMethod

std::vector< double >	fAreaNorms

bool	fUncompressWithPed

bool	fSkipInd

double	fIncludeMoreTail

int	fColMinWindow

int	fIndCutoff

Detailed Description

Definition at line 36 of file RawHitFinder_module.cc.

Member Typedef Documentation

using art::EDProducer::ModuleType = EDProducer

inherited

Definition at line 17 of file EDProducer.h.

template<typename UserConfig , typename KeysToIgnore = void>

using art::detail::Producer::Table = Modifier::Table<UserConfig, KeysToIgnore>

inherited

Definition at line 26 of file Producer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 70 of file RawHitFinder_module.cc.

References recob::HitAndAssociationsWriterBase::declare_products(), fAreaMethod, fAreaNorms, fCalDataModuleLabel, fColMinWidth, fColMinWindow, fColWidth, fDigitModuleLabel, fIncludeMoreTail, fIndCutoff, fIndMinWidth, fIndWidth, fMaxMultiHit, fMinSigCol, fMinSigInd, fSkipInd, fUncompressWithPed, and art::ProductRegistryHelper::producesCollector().

                                                           : EDProducer{pset}
   {
     fDigitModuleLabel = pset.get<art::InputTag>("DigitModuleLabel", "daq");
     fCalDataModuleLabel = pset.get<std::string>("CalDataModuleLabel");
     fMinSigInd = pset.get<double>("MinSigInd");
     fMinSigCol = pset.get<double>("MinSigCol");
     fIncludeMoreTail = pset.get<double>("IncludeMoreTail", 0.);
     fIndWidth = pset.get<int>("IndWidth", 20);
     fColWidth = pset.get<double>("ColWidth");
     fIndMinWidth = pset.get<double>("IndMinWidth");
     fColMinWidth = pset.get<double>("ColMinWidth", 0.);
     fMaxMultiHit = pset.get<int>("MaxMultiHit");
     fAreaMethod = pset.get<int>("AreaMethod");
     fAreaNorms = pset.get<std::vector<double>>("AreaNorms");
     fUncompressWithPed = pset.get<bool>("UncompressWithPed", true);
     fSkipInd = pset.get<bool>("SkipInd", false);
     fColMinWindow = pset.get<int>("ColMinWindow", 0);
     fIndCutoff = pset.get<int>("IndCutoff", 20);
     mf::LogInfo("RawHitFinder_module") << "fDigitModuleLabel: " << fDigitModuleLabel << std::endl;
 
     //LET HITCOLLECTIONCREATOR DECLARE THAT WE ARE GOING TO PRODUCE
     //HITS AND ASSOCIATIONS TO RAW DIGITS BUT NOT ASSOCIATIONS TO WIRES
     //(WITH NO PARTICULAR PRODUCT LABEL).
     recob::HitCollectionCreator::declare_products(producesCollector(),
                                                   /*instance_name*/ "",
                                                   /*doWireAssns*/ false,
                                                   /*doRawDigitAssns*/ true);
   }

Member Function Documentation

template<typename T , BranchType BT>

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

protectedinherited

Definition at line 61 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::consumes().

   {
     return collector_.consumes<T, BT>(tag);
   }

ConsumesCollector & art::ModuleBase::consumesCollector ( )

protectedinherited

Definition at line 57 of file ModuleBase.cc.

References art::ModuleBase::collector_.

   {
     return collector_;
   }

template<typename T , BranchType BT>

void art::ModuleBase::consumesMany ( )

protectedinherited

Definition at line 75 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::consumesMany().

   {
     collector_.consumesMany<T, BT>();
   }

template<typename Element , BranchType = InEvent>

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

protectedinherited

template<typename T , BranchType BT>

ViewToken<T> art::ModuleBase::consumesView ( InputTag const & tag )

inherited

Definition at line 68 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::consumesView().

   {
     return collector_.consumesView<T, BT>(tag);
   }

void art::detail::Producer::doBeginJob ( SharedResources const & resources )

inherited

Definition at line 22 of file Producer.cc.

References art::detail::Producer::beginJobWithFrame(), and art::detail::Producer::setupQueues().

   {
     setupQueues(resources);
     ProcessingFrame const frame{ScheduleID{}};
     beginJobWithFrame(frame);
   }

bool art::detail::Producer::doBeginRun	(	RunPrincipal &	rp,
		ModuleContext const &	mc
	)

inherited

Definition at line 65 of file Producer.cc.

References art::detail::Producer::beginRunWithFrame(), art::RangeSet::forRun(), art::RunPrincipal::makeRun(), r, art::RunPrincipal::runID(), and art::ModuleContext::scheduleID().

   {
     auto r = rp.makeRun(mc, RangeSet::forRun(rp.runID()));
     ProcessingFrame const frame{mc.scheduleID()};
     beginRunWithFrame(r, frame);
     r.commitProducts();
     return true;
   }

bool art::detail::Producer::doBeginSubRun	(	SubRunPrincipal &	srp,
		ModuleContext const &	mc
	)

inherited

Definition at line 85 of file Producer.cc.

References art::detail::Producer::beginSubRunWithFrame(), art::RangeSet::forSubRun(), art::SubRunPrincipal::makeSubRun(), art::ModuleContext::scheduleID(), and art::SubRunPrincipal::subRunID().

   {
     auto sr = srp.makeSubRun(mc, RangeSet::forSubRun(srp.subRunID()));
     ProcessingFrame const frame{mc.scheduleID()};
     beginSubRunWithFrame(sr, frame);
     sr.commitProducts();
     return true;
   }

void art::detail::Producer::doEndJob ( )

inherited

Definition at line 30 of file Producer.cc.

References art::detail::Producer::endJobWithFrame().

   {
     ProcessingFrame const frame{ScheduleID{}};
     endJobWithFrame(frame);
   }

bool art::detail::Producer::doEndRun	(	RunPrincipal &	rp,
		ModuleContext const &	mc
	)

inherited

Definition at line 75 of file Producer.cc.

References art::detail::Producer::endRunWithFrame(), art::RunPrincipal::makeRun(), r, art::ModuleContext::scheduleID(), and art::Principal::seenRanges().

   {
     auto r = rp.makeRun(mc, rp.seenRanges());
     ProcessingFrame const frame{mc.scheduleID()};
     endRunWithFrame(r, frame);
     r.commitProducts();
     return true;
   }

bool art::detail::Producer::doEndSubRun	(	SubRunPrincipal &	srp,
		ModuleContext const &	mc
	)

inherited

Definition at line 95 of file Producer.cc.

References art::detail::Producer::endSubRunWithFrame(), art::SubRunPrincipal::makeSubRun(), art::ModuleContext::scheduleID(), and art::Principal::seenRanges().

   {
     auto sr = srp.makeSubRun(mc, srp.seenRanges());
     ProcessingFrame const frame{mc.scheduleID()};
     endSubRunWithFrame(sr, frame);
     sr.commitProducts();
     return true;
   }

bool art::detail::Producer::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
	)

inherited

Definition at line 105 of file Producer.cc.

References art::detail::Producer::checkPutProducts_, e, art::EventPrincipal::makeEvent(), art::detail::Producer::produceWithFrame(), and art::ModuleContext::scheduleID().

   {
     auto e = ep.makeEvent(mc);
     ++counts_run;
     ProcessingFrame const frame{mc.scheduleID()};
     produceWithFrame(e, frame);
     e.commitProducts(checkPutProducts_, &expectedProducts<InEvent>());
     ++counts_passed;
     return true;
   }

void art::detail::Producer::doRespondToCloseInputFile ( FileBlock const & fb )

inherited

Definition at line 44 of file Producer.cc.

References art::detail::Producer::respondToCloseInputFileWithFrame().

   {
     ProcessingFrame const frame{ScheduleID{}};
     respondToCloseInputFileWithFrame(fb, frame);
   }

void art::detail::Producer::doRespondToCloseOutputFiles ( FileBlock const & fb )

inherited

Definition at line 58 of file Producer.cc.

References art::detail::Producer::respondToCloseOutputFilesWithFrame().

   {
     ProcessingFrame const frame{ScheduleID{}};
     respondToCloseOutputFilesWithFrame(fb, frame);
   }

void art::detail::Producer::doRespondToOpenInputFile ( FileBlock const & fb )

inherited

Definition at line 37 of file Producer.cc.

References art::detail::Producer::respondToOpenInputFileWithFrame().

   {
     ProcessingFrame const frame{ScheduleID{}};
     respondToOpenInputFileWithFrame(fb, frame);
   }

void art::detail::Producer::doRespondToOpenOutputFiles ( FileBlock const & fb )

inherited

Definition at line 51 of file Producer.cc.

References art::detail::Producer::respondToOpenOutputFilesWithFrame().

   {
     ProcessingFrame const frame{ScheduleID{}};
     respondToOpenOutputFilesWithFrame(fb, frame);
   }

void art::Modifier::fillProductDescriptions ( )

inherited

Definition at line 10 of file Modifier.cc.

References art::ProductRegistryHelper::fillDescriptions(), and art::ModuleBase::moduleDescription().

   {
     ProductRegistryHelper::fillDescriptions(moduleDescription());
   }

std::array< std::vector< ProductInfo >, NumBranchTypes > const & art::ModuleBase::getConsumables ( ) const

inherited

Definition at line 43 of file ModuleBase.cc.

References art::ModuleBase::collector_, and art::ConsumesCollector::getConsumables().

   {
     return collector_.getConsumables();
   }

std::unique_ptr< Worker > art::ModuleBase::makeWorker ( WorkerParams const & wp )

inherited

Definition at line 37 of file ModuleBase.cc.

References art::ModuleBase::doMakeWorker(), and art::NumBranchTypes.

   {
     return doMakeWorker(wp);
   }

template<typename T , BranchType BT>

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

protectedinherited

Definition at line 82 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::mayConsume().

   {
     return collector_.mayConsume<T, BT>(tag);
   }

template<typename T , BranchType BT>

void art::ModuleBase::mayConsumeMany ( )

protectedinherited

Definition at line 96 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::mayConsumeMany().

   {
     collector_.mayConsumeMany<T, BT>();
   }

template<typename Element , BranchType = InEvent>

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

protectedinherited

template<typename T , BranchType BT>

ViewToken<T> art::ModuleBase::mayConsumeView ( InputTag const & tag )

inherited

Definition at line 89 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::mayConsumeView().

   {
     return collector_.mayConsumeView<T, BT>(tag);
   }

ModuleDescription const & art::ModuleBase::moduleDescription ( ) const

inherited

Definition at line 13 of file ModuleBase.cc.

References art::errors::LogicError.

Referenced by art::OutputModule::doRespondToOpenInputFile(), art::OutputModule::doWriteEvent(), art::Modifier::fillProductDescriptions(), art::OutputModule::makePlugins_(), art::OutputWorker::OutputWorker(), reco::shower::LArPandoraModularShowerCreation::produce(), art::Modifier::registerProducts(), and art::OutputModule::registerProducts().

   {
     if (md_.has_value()) {
       return *md_;
     }
 
     throw Exception{errors::LogicError,
                     "There was an error while calling moduleDescription().\n"}
       << "The moduleDescription() base-class member function cannot be called\n"
          "during module construction.  To determine which module is "
          "responsible\n"
          "for calling it, find the '<module type>:<module "
          "label>@Construction'\n"
          "tag in the message prefix above.  Please contact artists@fnal.gov\n"
          "for guidance.\n";
   }

void hit::RawHitFinder::produce ( art::Event & evt )

overrideprivatevirtual

Implements art::EDProducer.

Definition at line 100 of file RawHitFinder_module.cc.

References raw::RawDigit::ADCs(), bin, raw::RawDigit::Channel(), geo::GeometryCore::ChannelToWire(), raw::RawDigit::Compression(), DEFINE_ART_MODULE, recob::HitCollectionCreator::emplace_back(), util::end(), fColMinWindow, fDataSize, fDigitModuleLabel, fIncludeMoreTail, fIndCutoff, fMinSigCol, fMinSigInd, fSkipInd, fUncompressWithPed, art::ProductRetriever::getByLabel(), raw::RawDigit::GetPedestal(), geo::kCollection, geo::kInduction, recob::HitCollectionCreator::put_into(), recob::HitCollectionCreator::reserve(), raw::RawDigit::Samples(), geo::GeometryCore::SignalType(), and raw::Uncompress().

   {
     //GET THE GEOMETRY.
     art::ServiceHandle<geo::Geometry const> geom;
 
     //READ IN THE DIGIT LIST OBJECTS(S).
     art::Handle<std::vector<raw::RawDigit>> digitVecHandle;
 
     bool retVal = evt.getByLabel(fDigitModuleLabel, digitVecHandle);
     if (retVal == true)
       mf::LogInfo("RawHitFinder_module")
         << "I got fDigitModuleLabel: " << fDigitModuleLabel << std::endl;
     else
       mf::LogWarning("RawHitFinder_module")
         << "Could not get fDigitModuleLabel: " << fDigitModuleLabel << std::endl;
 
     std::vector<float> holder; //HOLDS SIGNAL DATA.
     std::vector<short> rawadc; //UNCOMPRESSED ADC VALUES.
 
     // ###############################################
     // ### Making a ptr vector to put on the event ###
     // ###############################################
     // THIS CONTAINS THE HIT COLLECTION AND ITS ASSOCIATIONS TO WIRES AND RAW DIGITS.
     recob::HitCollectionCreator hcol(evt, false /* doWireAssns */, true /* doRawDigitAssns */);
 
     std::vector<float> startTimes; //STORES TIME OF WINDOW START.
     std::vector<float> maxTimes;   //STORES TIME OF LOCAL MAXIMUM.
     std::vector<float> endTimes;   //STORES TIME OF WINDOW END.
     std::vector<float> peakHeight; //STORES ADC COUNT AT THE MAXIMUM.
     std::vector<float> hitrms;     //STORES CHARGE WEIGHTED RMS OF TIME ACROSS THE HIT.
     std::vector<double> charge;    //STORES THE TOTAL CHARGE ASSOCIATED WITH THE HIT.
 
     uint32_t channel = 0; //CHANNEL NUMBER.
     double threshold = 0; //MINIMUM SIGNAL SIZE FOR ID'ING A HIT.
     double totSig = 0;
     double myrms = 0;
     double mynorm = 0;
     geo::SigType_t sigType = geo::kInduction;
     std::stringstream numConv;
 
     hcol.reserve(digitVecHandle->size());
     for (size_t rdIter = 0; rdIter < digitVecHandle->size(); ++rdIter) {
       holder.clear();
 
       //GET THE REFERENCE TO THE CURRENT raw::RawDigit.
       art::Ptr<raw::RawDigit> digitVec(digitVecHandle, rdIter);
       channel = digitVec->Channel();
       fDataSize = digitVec->Samples();
 
       rawadc.resize(fDataSize);
       holder.resize(fDataSize);
 
       //UNCOMPRESS THE DATA.
       if (fUncompressWithPed) {
         int pedestal = (int)digitVec->GetPedestal();
         raw::Uncompress(digitVec->ADCs(), rawadc, pedestal, digitVec->Compression());
       }
       else {
         raw::Uncompress(digitVec->ADCs(), rawadc, digitVec->Compression());
       }
 
       //GET THE LIST OF BAD CHANNELS.
       lariov::ChannelStatusProvider const& channelStatus =
         art::ServiceHandle<lariov::ChannelStatusService const>()->GetProvider();
 
       lariov::ChannelStatusProvider::ChannelSet_t const BadChannels = channelStatus.BadChannels();
 
       for (unsigned int bin = 0; bin < fDataSize; ++bin) {
         holder[bin] = (rawadc[bin] - digitVec->GetPedestal());
       }
 
       sigType = geom->SignalType(channel);
 
       peakHeight.clear();
       endTimes.clear();
       startTimes.clear();
       maxTimes.clear();
       charge.clear();
       hitrms.clear();
 
       bool channelSwitch = false;
 
       for (auto it = BadChannels.begin(); it != BadChannels.end(); it++) {
         if (channel == *it) {
           channelSwitch = true;
           break;
         }
       }
 
       if (channelSwitch == false) {
         // ###############################################
         // ###             Induction Planes            ###
         // ###############################################
 
         //THE INDUCTION PLANE METHOD HAS NOT YET BEEN MODIFIED AND TESTED FOR REAL DATA.
         // Or for detectors without a grid plane
         //
         if (sigType == geo::kInduction && !fSkipInd) {
           threshold = fMinSigInd;
           //    std::cout<< "Threshold is " << threshold << std::endl;
           // fitWidth = fIndWidth;
           // minWidth = fIndMinWidth;
           //    continue;
           float negthr = -1.0 * threshold;
           unsigned int bin = 1;
           float minadc = 0;
 
           // find the dips
           while (bin < (fDataSize - 1)) { // loop over ticks
             float thisadc = holder[bin];
             float nextadc = holder[bin + 1];
             if (thisadc < negthr &&
                 nextadc < negthr) { // new region, require two ticks above threshold
               //                    std::cout << "new region" << bin << " " << thisadc << std::endl;
               // step back to find zero crossing
               unsigned int place = bin;
               while (thisadc <= 0 && bin > 0) {
                 //              std::cout << bin << " " << thisadc << std::endl;
                 bin--;
                 thisadc = holder[bin];
               }
               float hittime = bin + thisadc / (thisadc - holder[bin + 1]);
               maxTimes.push_back(hittime);
 
               // step back more to find the hit start time
               uint32_t stop;
               if (fIndCutoff < (int)bin) { stop = bin - fIndCutoff; }
               else {
                 stop = 0;
               }
               while (thisadc < threshold && bin > stop) {
                 //              std::cout << bin << " " << thisadc << std::endl;
                 bin--;
                 thisadc = holder[bin];
               }
               if (bin >= 2) bin -= 2;
               while (thisadc > threshold && bin > stop) {
                 //              std::cout << bin << " " << thisadc << std::endl;
                 bin--;
                 thisadc = holder[bin];
               }
               startTimes.push_back(bin + 1);
               // now step forward from hit time to find end time, area of dip
               bin = place;
               thisadc = holder[bin];
               minadc = thisadc;
               bin++;
               totSig = fabs(thisadc);
               while (thisadc < negthr && bin < fDataSize) {
                 totSig += fabs(thisadc);
                 thisadc = holder[bin];
                 if (thisadc < minadc) minadc = thisadc;
                 bin++;
               }
               endTimes.push_back(bin - 1);
               peakHeight.push_back(-1.0 * minadc);
               charge.push_back(totSig);
               hitrms.push_back(5.0);
               //            std::cout << "TOTAL SIGNAL INDUCTION " << totSig << "  5.0" << std::endl;
               // std::cout << "filled end times " << bin-1 << "peak height vector size " << peakHeight.size() << std::endl;
 
               // don't look for a new hit until it returns to baseline
               while (thisadc < 0 && bin < fDataSize) {
                 //            std::cout << bin << " " << thisadc << std::endl;
                 bin++;
                 if (bin == fDataSize) break;
                 thisadc = holder[bin];
               }
             } // end region
             bin++;
           } // loop over ticks
         }
 
         // ###############################################
         // ###             Collection Plane            ###
         // ###############################################
 
         else if (sigType == geo::kCollection) {
           threshold = fMinSigCol;
 
           float madc = threshold;
           int ibin = 0;
           int start = 0;
           int end = 0;
           unsigned int bin = 0;
 
           while (bin < fDataSize) {
             float thisadc = holder[bin];
             madc = threshold;
             ibin = 0;
 
             if (thisadc > madc) {
               start = bin;
 
               if (thisadc > threshold && bin < fDataSize) {
                 while (thisadc > threshold && bin < fDataSize) {
                   if (thisadc > madc) {
                     ibin = bin;
                     madc = thisadc;
                   }
                   bin++;
                   if (bin == fDataSize) break;
                   thisadc = holder[bin];
                 }
               }
               else {
                 bin++;
               }
 
               end = bin - 1;
 
               if (start != end) {
                 maxTimes.push_back(ibin);
                 peakHeight.push_back(madc);
                 startTimes.push_back(start);
                 endTimes.push_back(end);
 
                 totSig = 0;
                 myrms = 0;
                 mynorm = 0;
 
                 int moreTail = std::ceil(fIncludeMoreTail * (end - start));
                 if (moreTail < fColMinWindow) moreTail = fColMinWindow;
 
                 for (int i = start - moreTail; i <= end + moreTail; i++) {
                   if (i < (int)(holder.size()) && i >= 0) {
                     float temp = ibin - i;
                     myrms += temp * temp * holder[i];
 
                     totSig += holder[i];
                   }
                 }
 
                 charge.push_back(totSig);
                 mynorm = totSig;
                 myrms /= mynorm;
                 hitrms.push_back(sqrt(myrms));
 
                 //PRE CHANGES MADE 04/14/16. A BOOTH, DUNE 35T.
                 /*
                    int moreTail = std::ceil(fIncludeMoreTail*(end-start));
 
                    for(int i = start-moreTail; i <= end+moreTail; i++)
                    {
                    totSig += holder[i];
                    float temp2 = holder[i]*holder[i];
                    mynorm += temp2;
                    float temp = ibin-i;
                    myrms += temp*temp*temp2;
                    }
 
                    charge.push_back(totSig);
                    myrms/=mynorm;
                    if((end-start+2*moreTail+1)!=0)
                    {
                    myrms/=(float)(end-start+2*moreTail+1);
                    hitrms.push_back(sqrt(myrms));
                    }
                    else
                    {
                    hitrms.push_back(sqrt(myrms));
                    }*/
               }
             }
             start = 0;
             end = 0;
             bin++;
           }
         }
       }
 
       int numHits(0);                   //NUMBER OF CONSECUTIVE HITS BEING FITTED.
       int hitIndex(0);                  //INDEX OF CURRENT HIT IN SEQUENCE.
       double amplitude(0), position(0); //FIT PARAMETERS.
       double start(0), end(0);
       double amplitudeErr(0), positionErr(0); //FIT ERRORS.
       double goodnessOfFit(0), chargeErr(0);  //CHI2/NDF and error on charge.
       double hrms(0);
 
       numHits = maxTimes.size();
       for (int i = 0; i < numHits; ++i) {
         amplitude = peakHeight[i];
         position = maxTimes[i];
         start = startTimes[i];
         end = endTimes[i];
         hrms = hitrms[i];
         amplitudeErr = -1;
         positionErr = 1.0;
         goodnessOfFit = -1;
         chargeErr = -1;
         totSig = charge[i];
 
         std::vector<geo::WireID> wids = geom->ChannelToWire(channel);
         geo::WireID wid = wids[0];
 
         if (start >= end) {
           mf::LogWarning("RawHitFinder_module")
             << "Hit start " << start << " is >= hit end " << end;
           continue;
         }
 
         recob::HitCreator hit(*digitVec,    //RAW DIGIT REFERENCE.
                               wid,          //WIRE ID.
                               start,        //START TICK.
                               end,          //END TICK.
                               hrms,         //RMS.
                               position,     //PEAK_TIME.
                               positionErr,  //SIGMA_PEAK_TIME.
                               amplitude,    //PEAK_AMPLITUDE.
                               amplitudeErr, //SIGMA_PEAK_AMPLITUDE.
                               totSig,       //HIT_INTEGRAL.
                               chargeErr,    //HIT_SIGMA_INTEGRAL.
                               std::accumulate(holder.begin() + (int)start,
                                               holder.begin() + (int)end,
                                               0.), //SUMMED CHARGE.
                               1,                   //MULTIPLICITY.
                               -1,                  //LOCAL_INDEX.
                               goodnessOfFit,       //WIRE ID.
                               int(end - start + 1) //DEGREES OF FREEDOM.
         );
         hcol.emplace_back(hit.move(), digitVec);
 
         ++hitIndex;
       }
     }
 
     hcol.put_into(evt);
   }

void art::Modifier::registerProducts ( ProductDescriptions & productsToRegister )

inherited

Definition at line 16 of file Modifier.cc.

References art::ModuleBase::moduleDescription(), and art::ProductRegistryHelper::registerProducts().

   {
     ProductRegistryHelper::registerProducts(productsToRegister,
                                             moduleDescription());
   }

void art::ModuleBase::setModuleDescription ( ModuleDescription const & md )

inherited

Definition at line 31 of file ModuleBase.cc.

References art::ModuleBase::md_.

   {
     md_ = md;
   }

void art::ModuleBase::sortConsumables ( std::string const & current_process_name )

inherited

Definition at line 49 of file ModuleBase.cc.

References art::ModuleBase::collector_, and art::ConsumesCollector::sortConsumables().

   {
     // Now that we know we have seen all the consumes declarations,
     // sort the results for fast lookup later.
     collector_.sortConsumables(current_process_name);
   }

Member Data Documentation

int hit::RawHitFinder::fAreaMethod

private

Definition at line 58 of file RawHitFinder_module.cc.

Referenced by RawHitFinder().

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

private

Definition at line 59 of file RawHitFinder_module.cc.

Referenced by RawHitFinder().

std::string hit::RawHitFinder::fCalDataModuleLabel

private

Definition at line 49 of file RawHitFinder_module.cc.

Referenced by RawHitFinder().

double hit::RawHitFinder::fColMinWidth

private

Definition at line 56 of file RawHitFinder_module.cc.

Referenced by RawHitFinder().

int hit::RawHitFinder::fColMinWindow

private

Definition at line 64 of file RawHitFinder_module.cc.

Referenced by produce(), and RawHitFinder().

double hit::RawHitFinder::fColWidth

private

Definition at line 54 of file RawHitFinder_module.cc.

Referenced by RawHitFinder().

unsigned int hit::RawHitFinder::fDataSize

private

Definition at line 44 of file RawHitFinder_module.cc.

Referenced by produce().

art::InputTag hit::RawHitFinder::fDigitModuleLabel

private

Definition at line 45 of file RawHitFinder_module.cc.

Referenced by produce(), and RawHitFinder().

std::string hit::RawHitFinder::fHitLabelName

private

Definition at line 50 of file RawHitFinder_module.cc.

double hit::RawHitFinder::fIncludeMoreTail

private

Definition at line 62 of file RawHitFinder_module.cc.

Referenced by produce(), and RawHitFinder().

int hit::RawHitFinder::fIndCutoff

private

Definition at line 65 of file RawHitFinder_module.cc.

Referenced by produce(), and RawHitFinder().

double hit::RawHitFinder::fIndMinWidth

private

Definition at line 55 of file RawHitFinder_module.cc.

Referenced by RawHitFinder().

double hit::RawHitFinder::fIndWidth

private

Definition at line 53 of file RawHitFinder_module.cc.

Referenced by RawHitFinder().

int hit::RawHitFinder::fMaxMultiHit

private

Definition at line 57 of file RawHitFinder_module.cc.

Referenced by RawHitFinder().

double hit::RawHitFinder::fMinSigCol

private

Definition at line 52 of file RawHitFinder_module.cc.

Referenced by produce(), and RawHitFinder().

double hit::RawHitFinder::fMinSigInd

private

Definition at line 51 of file RawHitFinder_module.cc.

Referenced by produce(), and RawHitFinder().

bool hit::RawHitFinder::fSkipInd

private

Definition at line 61 of file RawHitFinder_module.cc.

Referenced by produce(), and RawHitFinder().

std::string hit::RawHitFinder::fSpillName

private

Definition at line 46 of file RawHitFinder_module.cc.

bool hit::RawHitFinder::fUncompressWithPed

private

Definition at line 60 of file RawHitFinder_module.cc.

Referenced by produce(), and RawHitFinder().

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

larreco/v09_25_00/source/larreco/HitFinder/RawHitFinder_module.cc

Public Types

Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation