cmtool::CMManagerBase Class Reference

#include "CMManagerBase.h"

Inheritance diagram for cmtool::CMManagerBase:

Public Types
enum	CMMSGLevel_t { kPerMerging, kPerIteration, kPerEvent, kNone }
	Enum to specify message output level. More...

Public Member Functions
	CMManagerBase ()
	Default constructor. More...
virtual	~CMManagerBase ()=default
	Default destructor. More...
void	DebugMode (CMMSGLevel_t level)
	Method to enable debug mode (lots of couts) More...
void	ReportTimings (bool time_report=true)
	Method to enable timing profile cout. More...
void	Reset ()
	Method to reset itself. More...
void	AddPriorityAlgo (CPriorityAlgoBase *algo)
	Setter to add an algorithm for priority determination. More...
void	MergeTillConverge (bool doit=true)
	Switch to continue merging till converges. More...
void	SetClusters (util::GeometryUtilities const &gser, const std::vector< std::vector< util::PxHit >> &clusters)
	A simple method to add a cluster. More...
void	SetClusters (const std::vector< cluster::ClusterParamsAlg > &clusters)
	A simple method to add a cluster. More...
const std::vector< cluster::ClusterParamsAlg > &	GetInputClusters () const
	A getter for input clusters. More...
void	SetMinNHits (unsigned int n)
	A setter for minimum # of hits ... passed onto ClusterParamsAlg. More...
void	Process (util::GeometryUtilities const &gser)
	A method to execute the main action, to be called per event. More...
void	SetAnaFile (TFile *fout)
	A setter for an analysis output file. More...

Protected Member Functions
void	ComputePriority (const std::vector< cluster::ClusterParamsAlg > &clusters)
	Function to compute priority. More...
virtual void	EventBegin ()
	FMWK function called @ beginning of Process() More...
virtual void	IterationBegin ()
	FMWK function called @ beginning of iterative loop inside Process() More...
virtual bool	IterationProcess (util::GeometryUtilities const &gser)=0
	FMWK function called @ iterative loop inside Process() More...
virtual void	IterationEnd ()
	FMWK function called @ end of iterative loop inside Process() More...
virtual void	EventEnd ()
	FMWK function called @ end of Process() More...

Protected Attributes
bool	_time_report
	Timing verbosity flag. More...
unsigned int	_min_nhits
	Minimum number of hits: the limit set for ClusterParamsAlg. More...
CMMSGLevel_t	_debug_mode
	Debug mode switch. More...
std::vector< cluster::ClusterParamsAlg >	_in_clusters
	Input clusters. More...
::cmtool::CPriorityAlgoBase *	_priority_algo
	Priority algorithm. More...
TFile *	_fout
	Output analysis plot TFile. More...
std::multimap< float, size_t >	_priority
	Priority record. More...
bool	_merge_till_converge
	Iteration loop switch. More...
std::set< UChar_t >	_planes
	A holder for # of unique planes in the clusters, computed in ComputePriority() function. More...

Detailed Description

A class that instantiates merging algorithm(s) and run. The book-keeping of merged cluster sets are done by CMergeBookKeeper.

Definition at line 41 of file CMManagerBase.h.

Member Enumeration Documentation

enum cmtool::CMManagerBase::CMMSGLevel_t

Enum to specify message output level.

Enumerator
kPerMerging	Extremely verbose (cout per individual algorithm execution)
kPerIteration	Somewhat verbose (cout per merging iteration)
kPerEvent	Bit verbose (cout per event)
kNone	Normal.

Definition at line 44 of file CMManagerBase.h.

                      {
      kPerMerging,
      kPerIteration,
      kPerEvent,
      kNone
    };

Constructor & Destructor Documentation

cmtool::CMManagerBase::CMManagerBase

(

)

Default constructor.

Definition at line 16 of file CMManagerBase.cxx.

References _debug_mode, _fout, _merge_till_converge, _min_nhits, _priority_algo, _time_report, kNone, and Reset().

  {
    _fout = 0;
    _debug_mode = kNone;
    _priority_algo = nullptr;
    _min_nhits = 0;
    _merge_till_converge = false;
    Reset();
    _time_report = false;
  }

virtual cmtool::CMManagerBase::~CMManagerBase ( )

virtualdefault

Default destructor.

Member Function Documentation

void cmtool::CMManagerBase::AddPriorityAlgo ( CPriorityAlgoBase *  algo )

inline

Setter to add an algorithm for priority determination.

Definition at line 71 of file CMManagerBase.h.

Referenced by ShowerReco3D::ShowerReco3D().

{ _priority_algo = algo; }

void cmtool::CMManagerBase::ComputePriority ( const std::vector< cluster::ClusterParamsAlg > &  clusters )

protected

Function to compute priority.

Definition at line 155 of file CMManagerBase.cxx.

References _planes, _priority, _priority_algo, _time_report, and cmtool::CPriorityAlgoBase::Priority().

Referenced by cmtool::CMatchManager::IterationProcess(), and cmtool::CMergeManager::IterationProcess().

  {

    TStopwatch localWatch;
    localWatch.Start();

    _priority.clear();
    _planes.clear();

    if (!clusters.size()) return;

    // Priority is computed cluster-by-cluster. In case of two clusters having the same priority
    // value the one with lower cluster index gets the priority. Also, clusters with priority < 0
    // are not logged (assumed not to be used)

    for (size_t i = 0; i < clusters.size(); ++i) {

      size_t c_index = clusters.size() - i - 1;

      float priority = clusters.at(c_index).GetNHits();

      if (_priority_algo) { priority = _priority_algo->Priority(clusters.at(c_index)); }

      if (priority > 0) {

        _priority.insert(std::make_pair(priority, c_index));

        if (_planes.find(clusters.at(c_index).Plane()) == _planes.end())

          _planes.insert(clusters.at(c_index).Plane());
      }
    }

    if (_time_report)
      std::cout << Form("  CMManagerBase Time Report: ComputePriority = %g [s]",
                        localWatch.RealTime())
                << std::endl;
  }

void cmtool::CMManagerBase::DebugMode ( CMMSGLevel_t  level )

inline

Method to enable debug mode (lots of couts)

Definition at line 62 of file CMManagerBase.h.

{ _debug_mode = level; }

virtual void cmtool::CMManagerBase::EventBegin ( )

inlineprotectedvirtual

FMWK function called @ beginning of Process()

Reimplemented in cmtool::CMergeManager, and cmtool::CMatchManager.

Definition at line 100 of file CMManagerBase.h.

Referenced by Process().

{}

virtual void cmtool::CMManagerBase::EventEnd ( )

inlineprotectedvirtual

FMWK function called @ end of Process()

Reimplemented in cmtool::CMergeManager, and cmtool::CMatchManager.

Definition at line 112 of file CMManagerBase.h.

Referenced by Process().

{}

const std::vector<cluster::ClusterParamsAlg>& cmtool::CMManagerBase::GetInputClusters ( ) const

inline

A getter for input clusters.

Definition at line 84 of file CMManagerBase.h.

Referenced by showerreco::ShowerRecoManager::Process().

{ return _in_clusters; }

virtual void cmtool::CMManagerBase::IterationBegin ( )

inlineprotectedvirtual

FMWK function called @ beginning of iterative loop inside Process()

Reimplemented in cmtool::CMergeManager, and cmtool::CMatchManager.

Definition at line 103 of file CMManagerBase.h.

Referenced by Process().

{}

virtual void cmtool::CMManagerBase::IterationEnd ( )

inlineprotectedvirtual

FMWK function called @ end of iterative loop inside Process()

Reimplemented in cmtool::CMergeManager, and cmtool::CMatchManager.

Definition at line 109 of file CMManagerBase.h.

Referenced by Process().

{}

virtual bool cmtool::CMManagerBase::IterationProcess ( util::GeometryUtilities const &  gser )

protectedpure virtual

FMWK function called @ iterative loop inside Process()

Implemented in cmtool::CMergeManager, and cmtool::CMatchManager.

Referenced by Process().

void cmtool::CMManagerBase::MergeTillConverge ( bool  doit = true )

inline

Switch to continue merging till converges.

Definition at line 74 of file CMManagerBase.h.

References lar::dump::vector().

{ _merge_till_converge = doit; }

void cmtool::CMManagerBase::Process

(

util::GeometryUtilities const &

gser

)

A method to execute the main action, to be called per event.

Definition at line 89 of file CMManagerBase.cxx.

References _debug_mode, _in_clusters, _merge_till_converge, _time_report, EventBegin(), EventEnd(), IterationBegin(), IterationEnd(), IterationProcess(), and kPerIteration.

Referenced by cluster::ClusterMergeHelper::Process(), and showerreco::ShowerRecoManager::Reconstruct().

  {

    if (!(_in_clusters.size())) return;

    TStopwatch localWatch;

    localWatch.Start();

    EventBegin();

    if (_time_report)
      std::cout << Form("  CMManagerBase Time Report: EventBegin = %g [s]", localWatch.RealTime())
                << std::endl;

    bool keep_going = true;

    while (keep_going) {

      localWatch.Start();

      IterationBegin();

      if (_time_report)
        std::cout << Form("  CMManagerBase Time Report: IterationBegin = %g [s]",
                          localWatch.RealTime())
                  << std::endl;

      localWatch.Start();

      keep_going = IterationProcess(gser);

      if (_time_report)
        std::cout << Form("  CMManagerBase Time Report: IterationProcess = %g [s]",
                          localWatch.RealTime())
                  << std::endl;
      localWatch.Start();

      IterationEnd();

      if (_time_report)
        std::cout << Form("  CMManagerBase Time Report: IterationEnd = %g [s]",
                          localWatch.RealTime())
                  << std::endl;

      if (!_merge_till_converge) {

        if (_debug_mode <= kPerIteration)

          std::cout
            << "\033[93m  Iterative approach = OFF ... exiting from iteration loop. \033[00m"
            << std::endl;

        break;
      }
    }

    localWatch.Start();

    EventEnd();

    if (_time_report)
      std::cout << Form("  CMManagerBase Time Report: EventEnd = %g [s]", localWatch.RealTime())
                << std::endl;
  }

void cmtool::CMManagerBase::ReportTimings ( bool  time_report = true )

inline

Method to enable timing profile cout.

Definition at line 65 of file CMManagerBase.h.

References Reset().

{ _time_report = time_report; }

void cmtool::CMManagerBase::Reset

(

)

Method to reset itself.

Definition at line 27 of file CMManagerBase.cxx.

References _in_clusters, _planes, _priority_algo, and cmtool::CMAlgoBase::Reset().

Referenced by CMManagerBase(), cmtool::CMergeManager::Reset(), cmtool::CMatchManager::Reset(), and SetClusters().

  {
    _planes.clear();
    _in_clusters.clear();
    if (_priority_algo) _priority_algo->Reset();
  }

void cmtool::CMManagerBase::SetAnaFile ( TFile *  fout )

inline

A setter for an analysis output file.

Definition at line 93 of file CMManagerBase.h.

{ _fout = fout; }

void cmtool::CMManagerBase::SetClusters	(	util::GeometryUtilities const &	gser,
		const std::vector< std::vector< util::PxHit >> &	clusters
	)

A simple method to add a cluster.

Definition at line 34 of file CMManagerBase.cxx.

References _in_clusters, _min_nhits, _time_report, Initialize(), Reset(), cluster::ClusterParamsAlg::SetMinNHits(), and cluster::ClusterParamsAlg::SetVerbose().

Referenced by showerreco::ShowerRecoManager::Reconstruct().

  {

    TStopwatch localWatch;

    // Reset
    this->Reset();

    // Clear & fill cluster info

    _in_clusters.clear();

    _in_clusters.reserve(clusters.size());

    ::cluster::ClusterParamsAlg tmp_alg;
    tmp_alg.SetMinNHits(_min_nhits);
    tmp_alg.SetVerbose(false);

    for (auto const& c : clusters) {

      _in_clusters.push_back(tmp_alg);
      (*_in_clusters.rbegin()).Initialize();

      if ((*_in_clusters.rbegin()).SetHits(c) < 3) continue;
      (*_in_clusters.rbegin()).DisableFANN();
      (*_in_clusters.rbegin()).FillParams(gser, false, false, false, false, false, false);
      (*_in_clusters.rbegin()).FillPolygon(gser);
    }

    if (_time_report) {
      std::cout << Form("  CMManagerBase Time Report: SetClusters (CPAN computation) = %g [s]",
                        localWatch.RealTime())
                << " ... details below." << std::endl;

      for (auto const& c : _in_clusters)

        c.TimeReport(std::cout);
    }
  }

void cmtool::CMManagerBase::SetClusters

(

const std::vector< cluster::ClusterParamsAlg > &

clusters

)

A simple method to add a cluster.

Definition at line 75 of file CMManagerBase.cxx.

References _in_clusters, and _time_report.

  {
    TStopwatch localWatch;

    localWatch.Start();

    _in_clusters = clusters;

    if (_time_report)
      std::cout << Form("  CMManagerBase Time Report: SetClusters (copy) = %g [s]",
                        localWatch.RealTime())
                << std::endl;
  }

void cmtool::CMManagerBase::SetMinNHits ( unsigned int  n )

inline

A setter for minimum # of hits ... passed onto ClusterParamsAlg.

Definition at line 87 of file CMManagerBase.h.

References n.

{ _min_nhits = n; }

Member Data Documentation

CMMSGLevel_t cmtool::CMManagerBase::_debug_mode

protected

Debug mode switch.

Definition at line 122 of file CMManagerBase.h.

Referenced by CMManagerBase(), cmtool::CMatchManager::EventBegin(), cmtool::CMergeManager::EventBegin(), cmtool::CMergeManager::IterationBegin(), cmtool::CMergeManager::IterationEnd(), cmtool::CMatchManager::IterationProcess(), cmtool::CMergeManager::IterationProcess(), Process(), cmtool::CMergeManager::RunMerge(), and cmtool::CMergeManager::RunSeparate().

TFile* cmtool::CMManagerBase::_fout

protected

Output analysis plot TFile.

Definition at line 131 of file CMManagerBase.h.

Referenced by CMManagerBase(), and cmtool::CMergeManager::EventBegin().

std::vector<cluster::ClusterParamsAlg> cmtool::CMManagerBase::_in_clusters

protected

Input clusters.

Definition at line 125 of file CMManagerBase.h.

Referenced by cmtool::CMatchManager::EventBegin(), cmtool::CMergeManager::EventBegin(), cmtool::CMatchManager::IterationBegin(), cmtool::CMergeManager::IterationBegin(), cmtool::CMatchManager::IterationProcess(), cmtool::CMergeManager::IterationProcess(), Process(), Reset(), and SetClusters().

bool cmtool::CMManagerBase::_merge_till_converge

protected

Iteration loop switch.

Definition at line 137 of file CMManagerBase.h.

Referenced by CMManagerBase(), cmtool::CMergeManager::IterationEnd(), and Process().

unsigned int cmtool::CMManagerBase::_min_nhits

protected

Minimum number of hits: the limit set for ClusterParamsAlg.

Definition at line 119 of file CMManagerBase.h.

Referenced by CMManagerBase(), and SetClusters().

std::set<UChar_t> cmtool::CMManagerBase::_planes

protected

A holder for # of unique planes in the clusters, computed in ComputePriority() function.

Definition at line 140 of file CMManagerBase.h.

Referenced by ComputePriority(), cmtool::CMatchManager::IterationProcess(), cmtool::CMergeManager::IterationProcess(), and Reset().

std::multimap<float, size_t> cmtool::CMManagerBase::_priority

protected

Priority record.

Definition at line 134 of file CMManagerBase.h.

Referenced by ComputePriority(), cmtool::CMatchManager::IterationProcess(), and cmtool::CMergeManager::RunMerge().

::cmtool::CPriorityAlgoBase* cmtool::CMManagerBase::_priority_algo

protected

Priority algorithm.

Definition at line 128 of file CMManagerBase.h.

Referenced by CMManagerBase(), ComputePriority(), cmtool::CMatchManager::EventBegin(), cmtool::CMergeManager::EventBegin(), cmtool::CMatchManager::EventEnd(), cmtool::CMergeManager::EventEnd(), cmtool::CMatchManager::IterationBegin(), cmtool::CMergeManager::IterationBegin(), cmtool::CMatchManager::IterationEnd(), cmtool::CMergeManager::IterationEnd(), cmtool::CMatchManager::IterationProcess(), cmtool::CMatchManager::Reset(), and Reset().

bool cmtool::CMManagerBase::_time_report

protected

Timing verbosity flag.

Definition at line 116 of file CMManagerBase.h.

Referenced by CMManagerBase(), ComputePriority(), Process(), and SetClusters().

---

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

• larreco/v09_25_00/source/larreco/RecoAlg/CMTool/CMToolBase/CMManagerBase.h
• larreco/v09_25_00/source/larreco/RecoAlg/CMTool/CMToolBase/CMManagerBase.cxx