Inheritance diagram for nnet::PointIdTrainingData:

Classes
struct	Config

Public Types
using	Parameters = art::EDAnalyzer::Table< Config >

using	ModuleType = EDAnalyzer

Public Member Functions
	PointIdTrainingData (Parameters const &config)

void	beginJob () override

void	endJob () override

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)

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
std::string const &	processName () const

bool	wantAllEvents () const noexcept

bool	wantEvent (ScheduleID id, Event const &e) const

Handle< TriggerResults >	getTriggerResults (Event const &e) const

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	analyze (const art::Event &event) override

int	WeightedFit (int n, std::vector< double > const &x, std::vector< double > const &y, std::vector< double > const &w, double *params) const

Private Attributes
nnet::TrainingDataAlg	fTrainingDataAlg

std::string	fOutTextFilePath

std::string	fOutNumpyFileName

bool	fDumpToRoot

bool	fDumpToNumpy

std::vector< int >	fSelectedTPC

std::vector< int >	fSelectedPlane

int	fEvent

int	fRun
	number of the event being processed More...

int	fSubRun
	number of the run being processed More...

bool	fCrop
	number of the sub-run being processed More...

int	fPatch_size_w
	crop data to event (set to false when dumping noise!) More...

int	fPatch_size_d
	patch size in wire dimension More...

double	fEm
	patch size in drift dimension More...

double	fTrk

double	fMichel

double	fNone

double	fStopTrk

double	fCleanTrk

int	nEm

int	nTrk

int	nMichel

int	nNone

int	nEm_sel

int	nTrk_sel

int	nMichel_sel

int	nNone_sel

int	nStopTrk_sel

int	nCleanTrk_sel

geo::GeometryCore const *	fGeometry

c2numpy_writer	npywriter

CLHEP::HepRandomEngine &	fEngine
	art-managed random-number engine More...

Detailed Description

Definition at line 62 of file PointIdTrainingData_module.cc.

Member Typedef Documentation

using art::EDAnalyzer::ModuleType = EDAnalyzer

inherited

Definition at line 22 of file EDAnalyzer.h.

using nnet::PointIdTrainingData::Parameters = art::EDAnalyzer::Table<Config>

Definition at line 112 of file PointIdTrainingData_module.cc.

Constructor & Destructor Documentation

nnet::PointIdTrainingData::PointIdTrainingData ( PointIdTrainingData::Parameters const & config )

explicit

Definition at line 162 of file PointIdTrainingData_module.cc.

References fGeometry, fSelectedPlane, fSelectedTPC, geo::GeometryCore::MaxPlanes(), and geo::GeometryCore::NTPC().

     : art::EDAnalyzer(config)
     , fTrainingDataAlg(config().TrainingDataAlg())
     , fOutTextFilePath(config().OutTextFilePath())
     , fOutNumpyFileName(config().OutNumpyFileName())
     , fDumpToRoot(config().DumpToRoot())
     , fDumpToNumpy(config().DumpToNumpy())
     , fSelectedTPC(config().SelectedTPC())
     , fSelectedPlane(config().SelectedView())
     , fCrop(config().Crop())
     , fPatch_size_w(config().Patch_size_w())
     , fPatch_size_d(config().Patch_size_d())
     , fEm(config().Em())
     , fTrk(config().Trk())
     , fMichel(config().Michel())
     , fNone(config().None())
     , fStopTrk(config().StopTrk())
     , fCleanTrk(config().CleanTrk())
     , fEngine(art::ServiceHandle<rndm::NuRandomService>()->registerAndSeedEngine(createEngine(0)))
   {
     fGeometry = &*(art::ServiceHandle<geo::Geometry const>());
 
     if (fSelectedTPC.empty()) {
       for (size_t tpc = 0; tpc < fGeometry->NTPC(); ++tpc)
         fSelectedTPC.push_back(tpc);
     }
 
     if (fSelectedPlane.empty()) {
       for (size_t p = 0; p < fGeometry->MaxPlanes(); ++p)
         fSelectedPlane.push_back(p);
     }
   }

Member Function Documentation

void nnet::PointIdTrainingData::analyze ( const art::Event & event )

overrideprivate

Definition at line 244 of file PointIdTrainingData_module.cc.

References util::abs(), c2numpy_float32(), c2numpy_uint16(), c2numpy_uint32(), c2numpy_uint8(), d, e, edep, art::EventID::event(), fCleanTrk, fCrop, fDumpToNumpy, fDumpToRoot, fEm, fEngine, fEvent, nnet::TrainingDataAlg::findCrop(), fMichel, fNone, fOutTextFilePath, fPatch_size_d, fPatch_size_w, fRun, fSelectedPlane, fSelectedTPC, fStopTrk, fSubRun, fTrainingDataAlg, fTrk, art::Event::id(), nCleanTrk_sel, nEm, nEm_sel, nMichel, nMichel_sel, nNone, nNone_sel, npywriter, img::DataProviderAlg::NScaledDrifts(), nStopTrk_sel, nTrk, nTrk_sel, img::DataProviderAlg::NWires(), art::EventID::run(), nnet::TrainingDataAlg::saveSimInfo(), nnet::TrainingDataAlg::setEventData(), art::EventID::subRun(), w, WeightedFit(), img::DataProviderAlg::wireData(), nnet::TrainingDataAlg::wireEdep(), nnet::TrainingDataAlg::wirePdg(), y1, y2, and y3.

   {
     fEvent = event.id().event();
     fRun = event.run();
     fSubRun = event.subRun();
 
     bool saveSim = fTrainingDataAlg.saveSimInfo() && !event.isRealData();
 
     std::ostringstream os;
     os << "event_" << fEvent << "_run_" << fRun << "_subrun_" << fSubRun;
 
     std::cout << "analyze " << os.str() << std::endl;
 
     auto const clockData =
       art::ServiceHandle<detinfo::DetectorClocksService const>()->DataFor(event);
     auto const detProp =
       art::ServiceHandle<detinfo::DetectorPropertiesService const>()->DataFor(event, clockData);
 
     CLHEP::RandFlat flat(fEngine);
 
     for (size_t i = 0; i < fSelectedTPC.size(); ++i)
       for (size_t v = 0; v < fSelectedPlane.size(); ++v) {
         fTrainingDataAlg.setEventData(
           event, clockData, detProp, fSelectedPlane[v], fSelectedTPC[i], 0);
 
         unsigned int w0, w1, d0, d1;
         if (fCrop && saveSim) {
           if (fTrainingDataAlg.findCrop(0.004F, w0, w1, d0, d1)) {
             std::cout << "   crop: " << w0 << " " << w1 << " " << d0 << " " << d1 << std::endl;
           }
           else {
             std::cout << "   skip empty tpc:" << fSelectedTPC[i] << " / view:" << fSelectedPlane[v]
                       << std::endl;
             continue;
           }
         }
         else {
           w0 = 0;
           w1 = fTrainingDataAlg.NWires();
           d0 = 0;
           d1 = fTrainingDataAlg.NScaledDrifts();
         }
 
         if (fDumpToRoot) {
           std::ostringstream ss1;
           ss1 << "raw_" << os.str() << "_tpc_" << fSelectedTPC[i] << "_view_"
               << fSelectedPlane[v]; // TH2's name
 
           art::ServiceHandle<art::TFileService const> tfs;
           TH2F* rawHist =
             tfs->make<TH2F>((ss1.str() + "_raw").c_str(), "ADC", w1 - w0, w0, w1, d1 - d0, d0, d1);
           TH2F* depHist = 0;
           TH2I* pdgHist = 0;
           if (saveSim) {
             depHist = tfs->make<TH2F>(
               (ss1.str() + "_deposit").c_str(), "Deposit", w1 - w0, w0, w1, d1 - d0, d0, d1);
             pdgHist = tfs->make<TH2I>(
               (ss1.str() + "_pdg").c_str(), "PDG", w1 - w0, w0, w1, d1 - d0, d0, d1);
           }
 
           for (size_t w = w0; w < w1; ++w) {
             auto const& raw = fTrainingDataAlg.wireData(w);
             for (size_t d = d0; d < d1; ++d) {
               rawHist->Fill(w, d, raw[d]);
             }
 
             if (saveSim) {
               auto const& edep = fTrainingDataAlg.wireEdep(w);
               for (size_t d = d0; d < d1; ++d) {
                 depHist->Fill(w, d, edep[d]);
               }
 
               auto const& pdg = fTrainingDataAlg.wirePdg(w);
               for (size_t d = d0; d < d1; ++d) {
                 pdgHist->Fill(w, d, pdg[d]);
               }
             }
           }
 
           writeAndDelete(rawHist);
           writeAndDelete(depHist);
           writeAndDelete(pdgHist);
         }
         else if (fDumpToNumpy) {
           for (size_t w = w0; w < w1; ++w) {
             int w_start = w - fPatch_size_w / 2;
             int w_stop = w_start + fPatch_size_w;
             if (w_start < int(w0) || w_start > int(w1)) continue;
             if (w_stop < int(w0) || w_stop > int(w1)) continue;
             auto const& pdg = fTrainingDataAlg.wirePdg(w);
             auto const& deposit = fTrainingDataAlg.wireEdep(w);
             auto const& raw = fTrainingDataAlg.wireData(w);
             for (size_t d = d0; d < d1; ++d) {
               int d_start = d - fPatch_size_d / 2;
               int d_stop = d_start + fPatch_size_d;
               if (d_start < int(d0) || d_start > int(d1)) continue;
               if (d_stop < int(d0) || d_stop > int(d1)) continue;
               int y0 = 0, y1 = 0, y2 = 0, y3 = 0;
               if (deposit[d] < 2e-5 || raw[d] < 0.05) { //empty pixel
                 y3 = 1;
                 ++nNone;
                 if (flat.fire() > fNone) continue;
                 ++nNone_sel;
               }
               else if ((pdg[d] & 0x0FFF) == 11) { //shower
                 y1 = 1;
                 ++nEm;
                 if ((pdg[d] & 0xF000) == 0x2000) { //Michel
                   y2 = 1;
                   ++nMichel;
                   if (flat.fire() > fMichel) continue;
                   ++nMichel_sel;
                 }
                 else {
                   if (flat.fire() > fEm) continue;
                   ++nEm_sel;
                 }
               }
               else { //track
                 y0 = 1;
                 ++nTrk;
                 //Check if the track appears to be stopping
                 int nPxlw0 = 0, nPxlw1 = 0, nPxld0 = 0, nPxld1 = 0;
                 //Try to fit a track
                 std::vector<double> wfit, dfit, chgfit;
                 double total_trkchg = 0;
                 for (int ww = w_start; ww < w_stop; ++ww) {
                   auto const& pdg1 = fTrainingDataAlg.wirePdg(ww);
                   auto const& deposit1 = fTrainingDataAlg.wireEdep(ww);
                   auto const& raw1 = fTrainingDataAlg.wireData(ww);
                   for (int dd = d_start; dd < d_stop; ++dd) {
                     if (deposit1[dd] < 2e-5 || raw1[dd] < 0.05) continue; //empty pixel
                     if ((pdg1[dd] & 0x0FFF) != 11) {                      //track pixel
                       wfit.push_back(ww);
                       dfit.push_back(dd);
                       chgfit.push_back(raw1[dd]);
                       total_trkchg += raw1[dd];
                       if (ww - w_start < 3) ++nPxlw0;
                       if (w_stop - ww - 1 < 3) ++nPxlw1;
                       if (dd - d_start < 3) ++nPxld0;
                       if (d_stop - dd - 1 < 3) ++nPxld1;
                     }
                   }
                 }
                 //Fit track pixels
                 double fit_trkchg = 0;
                 double parm[2] = {};
                 if (!wfit.empty()) {
                   if (!WeightedFit(wfit.size(), wfit, dfit, chgfit, &parm[0])) {
                     for (size_t j = 0; j < wfit.size(); ++j) {
                       if (std::abs((dfit[j] - (parm[0] + wfit[j] * parm[1])) * cos(atan(parm[1]))) <
                           3) {
                         fit_trkchg += chgfit[j];
                       }
                     }
                   }
                   else if (!WeightedFit(dfit.size(), dfit, wfit, chgfit, &parm[0])) {
                     for (size_t j = 0; j < dfit.size(); ++j) {
                       if (std::abs((wfit[j] - (parm[0] + dfit[j] * parm[1])) * cos(atan(parm[1]))) <
                           3) {
                         fit_trkchg += chgfit[j];
                       }
                     }
                   }
                 }
                 if (((nPxlw0) && (!nPxlw1) && (!nPxld0) && (!nPxld1)) ||
                     ((!nPxlw0) && (nPxlw1) && (!nPxld0) && (!nPxld1)) ||
                     ((!nPxlw0) && (!nPxlw1) && (nPxld0) && (!nPxld1)) ||
                     ((!nPxlw0) && (!nPxlw1) && (!nPxld0) && (nPxld1))) { //stopping track
                   if (flat.fire() > fStopTrk) continue;
                   ++nStopTrk_sel;
                 }
                 else if (total_trkchg && fit_trkchg / total_trkchg > 0.9) { //clean track
                   if (flat.fire() > fCleanTrk) continue;
                   ++nCleanTrk_sel;
                 }
                 else {
                   if (flat.fire() > fTrk) continue;
                   ++nTrk_sel;
                 }
               }
               c2numpy_uint32(&npywriter, event.id().run());
               c2numpy_uint32(&npywriter, event.id().subRun());
               c2numpy_uint32(&npywriter, event.id().event());
               c2numpy_uint8(&npywriter, fSelectedTPC[i]);
               c2numpy_uint8(&npywriter, fSelectedPlane[v]);
               c2numpy_uint16(&npywriter, w);
               c2numpy_uint16(&npywriter, d);
               c2numpy_uint8(&npywriter, y0);
               c2numpy_uint8(&npywriter, y1);
               c2numpy_uint8(&npywriter, y2);
               c2numpy_uint8(&npywriter, y3);
               for (int ww = w_start; ww < w_stop; ++ww) {
                 auto const& raw1 = fTrainingDataAlg.wireData(ww);
                 for (int dd = d_start; dd < d_stop; ++dd) {
                   c2numpy_float32(&npywriter, raw1[dd]);
                 }
               }
             }
           }
         }
         else {
           std::ostringstream ss1;
           ss1 << fOutTextFilePath << "/raw_" << os.str() << "_tpc_" << fSelectedTPC[i] << "_view_"
               << fSelectedPlane[v];
 
           std::ofstream fout_raw, fout_deposit, fout_pdg;
 
           fout_raw.open(ss1.str() + ".raw");
           if (saveSim) {
             fout_deposit.open(ss1.str() + ".deposit");
             fout_pdg.open(ss1.str() + ".pdg");
           }
 
           for (size_t w = w0; w < w1; ++w) {
             auto const& raw = fTrainingDataAlg.wireData(w);
             for (size_t d = d0; d < d1; ++d) {
               fout_raw << raw[d] << " ";
             }
             fout_raw << std::endl;
 
             if (saveSim) {
               auto const& edep = fTrainingDataAlg.wireEdep(w);
               for (size_t d = d0; d < d1; ++d) {
                 fout_deposit << edep[d] << " ";
               }
               fout_deposit << std::endl;
 
               auto const& pdg = fTrainingDataAlg.wirePdg(w);
               for (size_t d = d0; d < d1; ++d) {
                 fout_pdg << pdg[d] << " ";
               }
               fout_pdg << std::endl;
             }
           }
 
           fout_raw.close();
           if (saveSim) {
             fout_deposit.close();
             fout_pdg.close();
           }
         }
       }
 
   } // PointIdTrainingData::analyze()

void nnet::PointIdTrainingData::beginJob ( )

overridevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 196 of file PointIdTrainingData_module.cc.

References c2numpy_addcolumn(), C2NUMPY_FLOAT32, c2numpy_init(), C2NUMPY_UINT16, C2NUMPY_UINT32, C2NUMPY_UINT8, fOutNumpyFileName, fPatch_size_d, fPatch_size_w, nCleanTrk_sel, nEm, nEm_sel, nMichel, nMichel_sel, nNone, nNone_sel, npywriter, nStopTrk_sel, nTrk, and nTrk_sel.

   {
 
     c2numpy_init(&npywriter, fOutNumpyFileName, 50000);
     c2numpy_addcolumn(&npywriter, "run", C2NUMPY_UINT32);
     c2numpy_addcolumn(&npywriter, "subrun", C2NUMPY_UINT32);
     c2numpy_addcolumn(&npywriter, "evt", C2NUMPY_UINT32);
     c2numpy_addcolumn(&npywriter, "tpc", C2NUMPY_UINT8);
     c2numpy_addcolumn(&npywriter, "plane", C2NUMPY_UINT8);
     c2numpy_addcolumn(&npywriter, "wire", C2NUMPY_UINT16);
     c2numpy_addcolumn(&npywriter, "tck", C2NUMPY_UINT16);
     c2numpy_addcolumn(&npywriter, "y0", C2NUMPY_UINT8);
     c2numpy_addcolumn(&npywriter, "y1", C2NUMPY_UINT8);
     c2numpy_addcolumn(&npywriter, "y2", C2NUMPY_UINT8);
     c2numpy_addcolumn(&npywriter, "y3", C2NUMPY_UINT8);
     for (int i = 0; i < fPatch_size_w * fPatch_size_d; ++i) {
       c2numpy_addcolumn(&npywriter, Form("x%d", i), C2NUMPY_FLOAT32);
     }
     nEm = 0;
     nTrk = 0;
     nMichel = 0;
     nNone = 0;
     nEm_sel = 0;
     nTrk_sel = 0;
     nMichel_sel = 0;
     nNone_sel = 0;
     nStopTrk_sel = 0;
     nCleanTrk_sel = 0;
   }

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::Analyzer::doBeginJob ( SharedResources const & resources )

inherited

Definition at line 25 of file Analyzer.cc.

Referenced by art::detail::Analyzer::Analyzer().

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

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

inherited

Definition at line 68 of file Analyzer.cc.

References art::ModuleContext::scheduleID().

Referenced by art::detail::Analyzer::Analyzer().

   {
     ProcessingFrame const frame{mc.scheduleID()};
     beginRunWithFrame(std::as_const(rp).makeRun(mc), frame);
     return true;
   }

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

inherited

Definition at line 84 of file Analyzer.cc.

References art::ModuleContext::scheduleID().

Referenced by art::detail::Analyzer::Analyzer().

   {
     ProcessingFrame const frame{mc.scheduleID()};
     beginSubRunWithFrame(std::as_const(srp).makeSubRun(mc), frame);
     return true;
   }

void art::detail::Analyzer::doEndJob ( )

inherited

Definition at line 33 of file Analyzer.cc.

Referenced by art::detail::Analyzer::Analyzer().

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

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

inherited

Definition at line 76 of file Analyzer.cc.

References art::ModuleContext::scheduleID().

Referenced by art::detail::Analyzer::Analyzer().

   {
     ProcessingFrame const frame{mc.scheduleID()};
     endRunWithFrame(std::as_const(rp).makeRun(mc), frame);
     return true;
   }

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

inherited

Definition at line 92 of file Analyzer.cc.

References art::ModuleContext::scheduleID().

Referenced by art::detail::Analyzer::Analyzer().

   {
     ProcessingFrame const frame{mc.scheduleID()};
     endSubRunWithFrame(std::as_const(srp).makeSubRun(mc), frame);
     return true;
   }

bool art::detail::Analyzer::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 100 of file Analyzer.cc.

References e, and art::ModuleContext::scheduleID().

Referenced by art::detail::Analyzer::Analyzer().

   {
     auto const e = std::as_const(ep).makeEvent(mc);
     if (wantEvent(mc.scheduleID(), e)) {
       ++counts_run;
       ProcessingFrame const frame{mc.scheduleID()};
       analyzeWithFrame(e, frame);
       ++counts_passed;
     }
     return true;
   }

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

inherited

Definition at line 47 of file Analyzer.cc.

Referenced by art::detail::Analyzer::Analyzer().

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

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

inherited

Definition at line 61 of file Analyzer.cc.

Referenced by art::detail::Analyzer::Analyzer().

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

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

inherited

Definition at line 40 of file Analyzer.cc.

Referenced by art::detail::Analyzer::Analyzer().

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

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

inherited

Definition at line 54 of file Analyzer.cc.

Referenced by art::detail::Analyzer::Analyzer().

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

void nnet::PointIdTrainingData::endJob ( )

overridevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 227 of file PointIdTrainingData_module.cc.

References c2numpy_close(), nCleanTrk_sel, nEm, nEm_sel, nMichel, nMichel_sel, nNone, nNone_sel, npywriter, nStopTrk_sel, nTrk, and nTrk_sel.

   {
     std::cout << "nEm = " << nEm << std::endl;
     std::cout << "nTrk = " << nTrk << std::endl;
     std::cout << "nMichel = " << nMichel << std::endl;
     std::cout << "nNone = " << nNone << std::endl;
     std::cout << std::endl;
     std::cout << "nEm_sel = " << nEm_sel << std::endl;
     std::cout << "nTrk_sel = " << nTrk_sel << std::endl;
     std::cout << "nMichel_sel = " << nMichel_sel << std::endl;
     std::cout << "nNone_sel = " << nNone_sel << std::endl;
     std::cout << "nStopTrk_sel = " << nStopTrk_sel << std::endl;
     std::cout << "nCleanTrk_sel = " << nCleanTrk_sel << std::endl;
     c2numpy_close(&npywriter);
   }

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

Handle< TriggerResults > art::Observer::getTriggerResults ( Event const & e ) const

protectedinherited

Definition at line 75 of file Observer.cc.

References art::ProductRetriever::get(), and art::Observer::selectors_.

Referenced by art::OutputModule::doWriteEvent(), and art::Observer::wantAllEvents().

   {
     if (selectors_) {
       return selectors_->getOneTriggerResults(e);
     }
 
     // The following applies for cases where no SelectEvents entries
     // exist.
     Handle<TriggerResults> h;
     if (e.get(empty_process_name, h)) {
       return h;
     }
     return Handle<TriggerResults>{};
   }

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

string const & art::Observer::processName ( ) const

protectedinherited

Definition at line 57 of file Observer.cc.

References art::Observer::process_name_.

Referenced by art::FileDumperOutput::printPrincipal().

   {
     return process_name_;
   }

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

bool art::Observer::wantAllEvents ( ) const

inlineprotectednoexceptinherited

Definition at line 31 of file Observer.h.

References e, art::Observer::getTriggerResults(), art::Observer::wantAllEvents_, and art::Observer::wantEvent().

     {
       return wantAllEvents_;
     }

bool art::Observer::wantEvent	(	ScheduleID	id,
		Event const &	e
	)		const

protectedinherited

Definition at line 63 of file Observer.cc.

References art::Observer::rejectors_, art::Observer::selectors_, and art::Observer::wantAllEvents_.

Referenced by art::OutputModule::doEvent(), art::OutputModule::doWriteEvent(), and art::Observer::wantAllEvents().

   {
     if (wantAllEvents_) {
       return true;
     }
     bool const select_event = selectors_ ? selectors_->matchEvent(id, e) : true;
     bool const reject_event =
       rejectors_ ? rejectors_->matchEvent(id, e) : false;
     return select_event and not reject_event;
   }

int nnet::PointIdTrainingData::WeightedFit	(	int	n,
		std::vector< double > const &	x,
		std::vector< double > const &	y,
		std::vector< double > const &	w,
		double *	params
	)		const

private

Definition at line 491 of file PointIdTrainingData_module.cc.

References DEFINE_ART_MODULE, and n.

Referenced by analyze().

   {
     Double_t sumx = 0.;
     Double_t sumx2 = 0.;
     Double_t sumy = 0.;
     Double_t sumy2 = 0.;
     Double_t sumxy = 0.;
     Double_t sumw = 0.;
     Double_t eparams[2] = {};
 
     for (Int_t i = 0; i < n; i++) {
       sumx += x[i] * w[i];
       sumx2 += x[i] * x[i] * w[i];
       sumy += y[i] * w[i];
       sumy2 += y[i] * y[i] * w[i];
       sumxy += x[i] * y[i] * w[i];
       sumw += w[i];
     }
 
     if (sumx2 * sumw - sumx * sumx == 0.) return 1;
     if (sumx2 - sumx * sumx / sumw == 0.) return 1;
 
     params[0] = (sumy * sumx2 - sumx * sumxy) / (sumx2 * sumw - sumx * sumx);
     params[1] = (sumxy - sumx * sumy / sumw) / (sumx2 - sumx * sumx / sumw);
 
     eparams[0] = sumx2 * (sumx2 * sumw - sumx * sumx);
     eparams[1] = (sumx2 - sumx * sumx / sumw);
 
     if (eparams[0] < 0. || eparams[1] < 0.) return 1;
 
     eparams[0] = sqrt(eparams[0]) / (sumx2 * sumw - sumx * sumx);
     eparams[1] = sqrt(eparams[1]) / (sumx2 - sumx * sumx / sumw);
 
     return 0;
   }