cmtool::CFAlgoZOverlap Class Reference

#include <CFAlgoZOverlap.h>

Inheritance diagram for cmtool::CFAlgoZOverlap:

Public Member Functions
	CFAlgoZOverlap ()
	Default constructor. More...
float	Float (util::GeometryUtilities const &, const std::vector< const cluster::ClusterParamsAlg * > &clusters) override
void	Report () override
void	Reset () override
	Function to reset the algorithm instance, called together with manager's Reset() More...
Public Member Functions inherited from cmtool::CMAlgoBase
	CMAlgoBase ()
virtual	~CMAlgoBase ()=default
virtual void	EventBegin (const std::vector< cluster::ClusterParamsAlg > &)
virtual void	EventEnd ()
virtual void	IterationBegin (const std::vector< cluster::ClusterParamsAlg > &)
virtual void	IterationEnd ()
void	SetAnaFile (TFile *fout)
	Setter function for an output plot TFile pointer. More...
virtual void	SetVerbose (bool doit=true)
	Setter function for verbosity. More...

Protected Attributes
float	_wire_ratio_cut
Protected Attributes inherited from cmtool::CMAlgoBase
TFile *	_fout
	TFile pointer to an output file. More...
bool	_verbose
	Boolean to choose verbose mode. Turned on if CMergeManager/CMatchManager's verbosity level is >= kPerMerging. More...

Detailed Description

User implementation for CFloatAlgoBase class doxygen documentation!

Definition at line 25 of file CFAlgoZOverlap.h.

Constructor & Destructor Documentation

cmtool::CFAlgoZOverlap::CFAlgoZOverlap

(

)

Default constructor.

Definition at line 6 of file CFAlgoZOverlap.cxx.

  {
    _wire_ratio_cut = 0.1; //Preliminary cuts
  }

Member Function Documentation

float cmtool::CFAlgoZOverlap::Float ( util::GeometryUtilities const &  , const std::vector< const cluster::ClusterParamsAlg * > &  clusters  )

overridevirtual

Core function: given a set of CPANs, return a float which indicates the compatibility the cluster combination.

Reimplemented from cmtool::CFloatAlgoBase.

Definition at line 20 of file CFAlgoZOverlap.cxx.

  {

    // This ensures the algorithm works only if # clusters is > 2 (and not =2)
    // You may take out this block if you want to allow matching using clusters from only 2 planes.
    if (clusters.size() == 2) return -1;

    double wire_distance = 0;
    double ratio = 1;
    double max_wire_distance = -1;

    //Record the start/end points that retunr the maximum wire distance
    double max_end_w = -1;

    double start_w = 0;
    double end_w = 0;
    _verbose = true;

    for (auto const& c : clusters) {

      //...start_point.w in planes 0 and 1 returns a distance in slanted wire space (perp to slanted wires).
      //Rotate this to properly compare to the other planes
      if (c->Plane() != 2) {
        start_w = 0.5 * c->GetParams().start_point.w;
        end_w = 0.5 * c->GetParams().end_point.w;
        wire_distance = end_w - start_w;
      }
      else {
        start_w = c->GetParams().start_point.w;
        end_w = c->GetParams().end_point.w;
        wire_distance = c->GetParams().end_point.w - c->GetParams().start_point.w;
      }

      if (wire_distance < 0) wire_distance *= -1;

      if (max_wire_distance < wire_distance) {
        max_wire_distance = wire_distance;
        //max_plane   = c->Plane();
        //max_start_w = start_w ;
        max_end_w = end_w;
      }
    }

    //Calculate maximum z range(accounting for the slant in UV). Then compare start points. Similar
    //in this sense, to time.

    for (auto const& c : clusters) {

      if (c->Plane() != 2) {
        start_w = 0.5 * c->GetParams().start_point.w;
        end_w = 0.5 * c->GetParams().end_point.w;
        wire_distance = end_w - start_w;
      }
      else {
        start_w = c->GetParams().start_point.w;
        end_w = c->GetParams().end_point.w;
        wire_distance = c->GetParams().end_point.w - c->GetParams().start_point.w;
      }

      if (wire_distance < 0) wire_distance *= -1;

      if (start_w <= max_end_w) // && end_w+25 >=max_start_w )
        ratio *= wire_distance / max_wire_distance;
      else
        ratio *= 0.1;

      if (_verbose && ratio > _wire_ratio_cut) {
        std::cout << "\nThe wire distance for cluster in plane " << c->Plane()
                  << " is: " << wire_distance << std::endl;
        std::cout << "Max wire disatance is: " << max_wire_distance << std::endl;
        std::cout << "Ratio is: " << ratio << std::endl;
        std::cout << "Start and end points: " << start_w << ",  " << end_w << std::endl;
      }
    }
    if (_verbose && ratio > _wire_ratio_cut)
      std::cout << " FOOOOUUUUNNNND ONE WOooooooooooooooooooooooooooooooooooooooooooooooooo: "
                << ratio << std::endl;

    return (ratio > _wire_ratio_cut ? ratio : -1);
  }

void cmtool::CFAlgoZOverlap::Report ( )

overridevirtual

Optional function: called after each iterative approach if a manager class is run with verbosity level <= kPerIteration. Maybe useful for debugging.

Reimplemented from cmtool::CMAlgoBase.

Definition at line 105 of file CFAlgoZOverlap.cxx.

  {}

void cmtool::CFAlgoZOverlap::Reset ( void  )

overridevirtual

Function to reset the algorithm instance, called together with manager's Reset()

Reimplemented from cmtool::CMAlgoBase.

Definition at line 14 of file CFAlgoZOverlap.cxx.

  {}

Member Data Documentation

float cmtool::CFAlgoZOverlap::_wire_ratio_cut

protected

Optional function: called at the beginning of 1st iteration. This is called per event. Optional function: called at the end of event ... after the last merging iteration is over. Optional function: called at the beggining of each iterative loop. This provides all clusters' information in case the algorithm need them. Note this is called per iteration which may be more than once per event. Optional function: called at the end of each iterative loop.

Definition at line 76 of file CFAlgoZOverlap.h.

---

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

• larreco/larreco/RecoAlg/CMTool/CMTAlgMatch/CFAlgoZOverlap.h
• larreco/larreco/RecoAlg/CMTool/CMTAlgMatch/CFAlgoZOverlap.cxx