Inheritance diagram for reco3d::PlotSpacePoints:

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

Public Member Functions inherited from art::EDAnalyzer
	EDAnalyzer (fhicl::ParameterSet const &pset)

template<typename Config >
	EDAnalyzer (Table< Config > const &config)

std::string	workerType () const

Public Member Functions inherited from art::detail::Analyzer
virtual	~Analyzer () noexcept

	Analyzer (fhicl::ParameterSet const &pset)

template<typename Config >
	Analyzer (Table< Config > const &config)

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)

Public Member Functions inherited from art::Observer
	~Observer () noexcept

	Observer (Observer const &)=delete

	Observer (Observer &&)=delete

Observer &	operator= (Observer const &)=delete

Observer &	operator= (Observer &&)=delete

void	registerProducts (ProductDescriptions &, ModuleDescription const &)

void	fillDescriptions (ModuleDescription const &)

fhicl::ParameterSetID	selectorConfig () const

Public Member Functions inherited from art::ModuleBase
virtual	~ModuleBase () noexcept

	ModuleBase ()

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)

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

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

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

void	Plot (const std::vector< recob::SpacePoint > &pts, const std::string &suffix) const

void	Plot3D (const std::vector< recob::SpacePoint > &pts, const std::string &suffix) const

std::vector< recob::SpacePoint >	TrueSpacePoints (detinfo::DetectorClocksData const &clockData, art::Handle< std::vector< recob::Hit >> hits) const

Private Attributes
art::InputTag	fSpacePointTag

std::string	fHitLabel

std::string	fSuffix

bool	fPlots

bool	fPlots3D

bool	fPlotsTrue

Additional Inherited Members
Public Types inherited from art::EDAnalyzer
using	WorkerType = WorkerT< EDAnalyzer >

using	ModuleType = EDAnalyzer

Protected Member Functions inherited from art::Observer
std::string const &	processName () const

bool	wantAllEvents () const noexcept

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

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

	Observer (fhicl::ParameterSet const &config)

	Observer (std::vector< std::string > const &select_paths, std::vector< std::string > const &reject_paths, fhicl::ParameterSet const &config)

Protected Member Functions inherited from art::ModuleBase
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 ()

Detailed Description

Definition at line 24 of file PlotSpacePoints_module.cc.

Constructor & Destructor Documentation

reco3d::PlotSpacePoints::PlotSpacePoints ( const fhicl::ParameterSet & pset )

explicit

Definition at line 52 of file PlotSpacePoints_module.cc.

     : EDAnalyzer(pset)
     , fSpacePointTag(art::InputTag(pset.get<std::string>("SpacePointLabel"),
                                    pset.get<std::string>("SpacePointInstanceLabel")))
     , fHitLabel(pset.get<std::string>("HitLabel"))
     , fSuffix(pset.get<std::string>("Suffix"))
     , fPlots(pset.get<bool>("Plots"))
     , fPlots3D(pset.get<bool>("Plots3D"))
     , fPlotsTrue(pset.get<bool>("PlotsTrue"))
   {
     if (!fSuffix.empty()) fSuffix = "_" + fSuffix;
   }

Member Function Documentation

void reco3d::PlotSpacePoints::analyze ( const art::Event & evt )

overrideprivate

Definition at line 173 of file PlotSpacePoints_module.cc.

   {
     if (fPlots) {
       art::Handle<std::vector<recob::SpacePoint>> pts;
       evt.getByLabel(fSpacePointTag, pts);
 
       const std::string suffix = TString::Format("%s_%d", fSuffix.c_str(), evt.event()).Data();
 
       if (!pts->empty()) {
         Plot(*pts, suffix);
 
         if (fPlots3D) Plot3D(*pts, suffix);
       }
     }
 
     if (fPlotsTrue) {
       art::Handle<std::vector<recob::Hit>> hits;
       evt.getByLabel(fHitLabel, hits);
       auto const clockData =
         art::ServiceHandle<detinfo::DetectorClocksService const>()->DataFor(evt);
       const std::vector<recob::SpacePoint> pts = TrueSpacePoints(clockData, hits);
 
       const std::string suffix = TString::Format("%s_true_%d", fSuffix.c_str(), evt.event()).Data();
 
       Plot(pts, suffix);
 
       if (fPlots3D) Plot3D(pts, suffix);
     }
   }

void reco3d::PlotSpacePoints::Plot	(	const std::vector< recob::SpacePoint > &	pts,
		const std::string &	suffix
	)		const

private

Definition at line 67 of file PlotSpacePoints_module.cc.

   {
     TGraph gZX;
     TGraph gYX;
     TGraph gZY;
 
     gZX.SetTitle(";z;x");
     gYX.SetTitle(";y;x");
     gZY.SetTitle(";z;y");
 
     for (const recob::SpacePoint& pt : pts) {
       const double* xyz = pt.XYZ();
       const double x = xyz[0];
       const double y = xyz[1];
       const double z = xyz[2];
       gZX.SetPoint(gZX.GetN(), z, x);
       gYX.SetPoint(gYX.GetN(), y, x);
       gZY.SetPoint(gZY.GetN(), z, y);
     }
 
     if (gZX.GetN() == 0) gZX.SetPoint(0, 0, 0);
     if (gYX.GetN() == 0) gYX.SetPoint(0, 0, 0);
     if (gZY.GetN() == 0) gZY.SetPoint(0, 0, 0);
 
     gZX.Draw("ap");
     gPad->Print(("plots/evd" + suffix + ".png").c_str());
 
     gYX.Draw("ap");
     gPad->Print(("plots/evd_ortho" + suffix + ".png").c_str());
     gZY.Draw("ap");
     gPad->Print(("plots/evd_zy" + suffix + ".png").c_str());
   }

void reco3d::PlotSpacePoints::Plot3D	(	const std::vector< recob::SpacePoint > &	pts,
		const std::string &	suffix
	)		const

private

Definition at line 126 of file PlotSpacePoints_module.cc.

   {
     int frame = 0;
     for (int phase = 0; phase < 4; ++phase) {
       const int Nang = 20;
       for (int iang = 0; iang < Nang; ++iang) {
         const double ang = M_PI / 2 * iang / double(Nang);
 
         TGraph g;
 
         for (const recob::SpacePoint& p : pts) {
           const double* xyz = p.XYZ();
 
           double x{}, y{};
           if (phase == 0) {
             x = cos(ang) * xyz[1] + sin(ang) * xyz[2];
             y = xyz[0];
           }
           if (phase == 1) {
             x = xyz[2];
             y = cos(ang) * xyz[0] + sin(ang) * xyz[1];
           }
           if (phase == 2) {
             x = cos(ang) * xyz[2] - sin(ang) * xyz[0];
             y = xyz[1];
           }
           if (phase == 3) {
             x = -cos(ang) * xyz[0] + sin(ang) * xyz[1];
             y = cos(ang) * xyz[1] + sin(ang) * xyz[0];
           }
 
           //        const double phi = phase/3.*M_PI/2 + ang/3;
           const double phi = 0;
           g.SetPoint(g.GetN(), cos(phi) * x + sin(phi) * y, cos(phi) * y - sin(phi) * x);
         }
 
         std::string fname =
           TString::Format(("anim/evd3d" + suffix + "_%03d.png").c_str(), frame++).Data();
         g.SetTitle(fname.c_str());
         if (g.GetN()) g.Draw("ap");
         gPad->Print(fname.c_str());
       }
     }
   }

std::vector< recob::SpacePoint > reco3d::PlotSpacePoints::TrueSpacePoints	(	detinfo::DetectorClocksData const &	clockData,
		art::Handle< std::vector< recob::Hit >>	hits
	)		const

private

Definition at line 102 of file PlotSpacePoints_module.cc.

   {
     std::vector<recob::SpacePoint> pts_true;
 
     const double err[6] = {
       0,
     };
 
     art::ServiceHandle<cheat::BackTrackerService const> bt_serv;
     for (unsigned int i = 0; i < hits->size(); ++i) {
       try {
         const std::vector<double> xyz = bt_serv->HitToXYZ(clockData, art::Ptr<recob::Hit>(hits, i));
         pts_true.emplace_back(&xyz[0], err, 0);
       }
       catch (...) {
       } // some hits have no electrons?
     }
 
     return pts_true;
   }

Member Data Documentation

std::string reco3d::PlotSpacePoints::fHitLabel

private

Definition at line 40 of file PlotSpacePoints_module.cc.

bool reco3d::PlotSpacePoints::fPlots

private

Definition at line 44 of file PlotSpacePoints_module.cc.

bool reco3d::PlotSpacePoints::fPlots3D

private

Definition at line 45 of file PlotSpacePoints_module.cc.

bool reco3d::PlotSpacePoints::fPlotsTrue

private

Definition at line 46 of file PlotSpacePoints_module.cc.

art::InputTag reco3d::PlotSpacePoints::fSpacePointTag

private

Definition at line 38 of file PlotSpacePoints_module.cc.

std::string reco3d::PlotSpacePoints::fSuffix

private

Definition at line 42 of file PlotSpacePoints_module.cc.

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

larreco/larreco/SpacePointSolver/PlotSpacePoints_module.cc

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation