cmtool::CFAlgoStartTimeCompat Class Reference

#include <CFAlgoStartTimeCompat.h>

Inheritance diagram for cmtool::CFAlgoStartTimeCompat:

Public Member Functions
	CFAlgoStartTimeCompat ()
	Default constructor. More...
virtual	~CFAlgoStartTimeCompat ()
	Default destructor. More...
virtual float	Float (const std::vector< const cluster::ClusterParamsAlg * > &clusters)
virtual void	Report ()
virtual void	Reset ()
	Function to reset the algorithm instance, called together with manager's Reset() More...
void	SetVerbose (bool on)
	Function to set verbose output. More...
void	SetTimeDist (double t)
	Function to set verbose output. More...
Public Member Functions inherited from cmtool::CFloatAlgoBase
	CFloatAlgoBase ()
	Default constructor. More...
virtual	~CFloatAlgoBase ()
	Default destructor. More...
Public Member Functions inherited from cmtool::CMAlgoBase
	CMAlgoBase ()
	Default constructor. More...
virtual	~CMAlgoBase ()
	Default destructor. More...
virtual void	EventBegin (const std::vector< cluster::ClusterParamsAlg > &clusters)
virtual void	EventEnd ()
virtual void	IterationBegin (const std::vector< cluster::ClusterParamsAlg > &clusters)
virtual void	IterationEnd ()
void	SetAnaFile (TFile *fout)
	Setter function for an output plot TFile pointer. More...

Private Attributes
double	_w2cm
double	_t2cm
bool	_verbose
double	_timeDist

Additional Inherited Members
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 CFAlgoStartTimeCompat.h.

Constructor & Destructor Documentation

cmtool::CFAlgoStartTimeCompat::CFAlgoStartTimeCompat

(

)

Default constructor.

Definition at line 10 of file CFAlgoStartTimeCompat.cxx.

                                               : CFloatAlgoBase()
  //-------------------------------------------------------
  {
    ::util::GeometryUtilities geou;
    _w2cm = geou.WireToCm();
    _t2cm = geou.TimeToCm();
    SetVerbose(false);
    SetTimeDist(5.);//in cm
  }

virtual cmtool::CFAlgoStartTimeCompat::~CFAlgoStartTimeCompat ( )

inlinevirtual

Default destructor.

Definition at line 33 of file CFAlgoStartTimeCompat.h.

{};

Member Function Documentation

float cmtool::CFAlgoStartTimeCompat::Float ( const std::vector< const cluster::ClusterParamsAlg * > &  clusters )

virtual

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

Reimplemented from cmtool::CFloatAlgoBase.

Definition at line 27 of file CFAlgoStartTimeCompat.cxx.

  {

    // Code-block by Kazu starts
    // 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;
    // Code-block by Kazu ends

    //This algorithm now works for 3 planes: find 3Dstart point from first 2 planes and find
    //How well that agrees with 3rd plane's start point location.

    //So first, make sure clusters vector has only 3 elements. If not return -1
    if ( clusters.size() != 3 )
      return -1;

    art::ServiceHandle<geo::Geometry const> geo;

    //Get Start Wire [cm], Time [cm], Plane and Channel for each cluster in the 3 planes.
    double startWirecm0 = clusters.at(0)->GetParams().start_point.w;
    double startTimecm0 = clusters.at(0)->GetParams().start_point.t;
    int startWire0 = int( startWirecm0 / _w2cm );
    unsigned char Pl0 = clusters.at(0)->GetParams().start_point.plane;
    UInt_t startChan0 = geo->PlaneWireToChannel(Pl0, startWire0);

    double startWirecm1 = clusters.at(1)->GetParams().start_point.w;
    double startTimecm1 = clusters.at(1)->GetParams().start_point.t;
    int startWire1 = int( startWirecm1 / _w2cm );
    unsigned char Pl1 = clusters.at(1)->GetParams().start_point.plane;
    UInt_t startChan1 = geo->PlaneWireToChannel(Pl1, startWire1);

    double startWirecm2 = clusters.at(2)->GetParams().start_point.w;
    double startTimecm2 = clusters.at(2)->GetParams().start_point.t;
    int startWire2 = int( startWirecm2 / _w2cm );
    unsigned char Pl2 = clusters.at(2)->GetParams().start_point.plane;
    UInt_t startChan2 = geo->PlaneWireToChannel(Pl2, startWire2);

    //Get Intersections in pairs:
    //y and z indicate detector coordinate and numbers indicate planes
    //used to generate that intersection point
    double yS01, zS01, yS02, zS02, yS12, zS12;

    //Check if Channels Intersect. If so get [Y,Z] intersection points
    bool WireIntersect01 = geo->ChannelsIntersect( startChan0, startChan1, yS01, zS01);
    bool WireIntersect02 = geo->ChannelsIntersect( startChan0, startChan2, yS02, zS02);
    bool WireIntersect12 = geo->ChannelsIntersect( startChan1, startChan2, yS12, zS12);

    // Check for 2-Plane Start-Point compatibility:
    // If Start times within pre-defined range
    // AND wires intersect --> Start Point is Good
    // Measure of compatibility is time-distance
    double compat01 = 9999, compat02 = 9999, compat12 = 9999;
    if ( (abs(startTimecm0 - startTimecm1) < _timeDist)
         and WireIntersect01 )
      compat01 = 0.01*abs(100*(startTimecm0-startTimecm1));
    if ( (abs(startTimecm1 - startTimecm2) < _timeDist)
         and WireIntersect12 )
      compat12 = 0.01*abs(100*(startTimecm1-startTimecm2));
    if ( (abs(startTimecm0 - startTimecm2) < _timeDist)
         and WireIntersect02 )
      compat02 = 0.01*abs(100*(startTimecm0-startTimecm2));

    // If no compatibility return -1
    if ( (compat01 > 1000) and (compat02 > 1000) and (compat12 > 1000) ){
      if ( _verbose ) { std::cout << "No compatibility among any 2 planes..." << std::endl; }
      return -1;
    }

    // If any of the planes were compatible, return "best"
    // This is the one where time-separation is smallest
    int bestcompat;
    double bestcompatVal=0;
    if ( (compat01 < compat02) and (compat01 < compat12) and (compat01 < 9999) ){
      bestcompat = 2;
      bestcompatVal = compat01;
    }
    if ( (compat02 < compat01) and (compat02 < compat12) and (compat02 < 9999) ){
      bestcompat = 1;
      bestcompatVal = compat02;
    }
    if ( (compat12 < compat02) and (compat12 < compat01) and (compat12 < 9999) ){
      bestcompat = 0;
      bestcompatVal = compat12;
    }

    //Now we know which two planes have a good match for start point (if any)
    //Check if the cluster on the 3rd plane is compatible with these two planes

    // To do this check that reconstructed start point from best
    // combination is inside cluster 3
    double minWire3rdClus = 9999, minTime3rdClus = 9999;
    double maxWire3rdClus = 0, maxTime3rdClus = 0;
    double hitWireTMP, hitTimeTMP;
    Double_t Start[3] = {-100., 0., 0.};
    double StartT3rdPlane = 0;
    double wireStart3rdPlane = 0;
    UInt_t Pl3rd = 0;
    std::vector<util::PxHit> hits;

    if ( bestcompat == 2 ){
      StartT3rdPlane = startTimecm2;//(startTimecm0+startTimecm1)/2.;
      Pl3rd = Pl2;
      Start[1] = yS01;
      Start[2] = zS01;
      hits = clusters.at(2)->GetHitVector();
    }
    if ( bestcompat == 0 ){
      StartT3rdPlane = startTimecm0;//(startTimecm1+startTimecm2)/2.;
      Pl3rd = Pl0;
      Start[1] = yS12;
      Start[2] = zS12;
      hits = clusters.at(0)->GetHitVector();
    }
    if ( bestcompat == 1 ){
      StartT3rdPlane = startTimecm1;//(startTimecm0+startTimecm2)/2.;
      Pl3rd = Pl1;
      Start[1] = yS02;
      Start[2] = zS02;
      hits = clusters.at(1)->GetHitVector();
    }

    for (auto& h: hits){
      hitWireTMP = h.w/_w2cm;
      hitTimeTMP = h.t;
      if ( hitWireTMP < minWire3rdClus ) { minWire3rdClus = hitWireTMP; }
      if ( hitWireTMP > maxWire3rdClus ) { maxWire3rdClus = hitWireTMP; }
      if ( hitTimeTMP < minTime3rdClus ) { minTime3rdClus = hitTimeTMP; }
      if ( hitTimeTMP > maxTime3rdClus ) { maxTime3rdClus = hitTimeTMP; }
    }
    try { wireStart3rdPlane = geo->NearestWire( Start, Pl3rd); }
    //catch ( ::larutil::LArUtilException &e) {
    catch ( ::cet::exception &e) {
      std::cout << e.what() << std::endl;
      std::cout << "Exception caught!" << std::endl;
      wireStart3rdPlane = -1;
    }

    if (_verbose){
      std::cout << "Cluster Start Times:" << std::endl;
      std::cout << "Pl0: " << startTimecm0 << "\tPl1: " << startTimecm1 << "\tPl2: " << startTimecm2 << std::endl;
      std::cout << "Best Time Separation: " << bestcompatVal << std::endl;
      std::cout << "Third Plane: " << Pl3rd << std::endl;
      std::cout << "Time Range for 3rd plane: [ " << minTime3rdClus << ", " << maxTime3rdClus << " ]" << std::endl;
      std::cout << "Start Time for 3rd Plane: " << StartT3rdPlane << std::endl;
      std::cout << "Reco z,y start point: [ " << Start[2] << ", " << Start[1] << " ]" << std::endl;
      std::cout << "Nearest Wire on 3rd Plane: " << wireStart3rdPlane << std::endl;
      std::cout << "Wire range 3rd Plane: [ " << minWire3rdClus << ", " << maxWire3rdClus << " ]" << std::endl;
    }

    //Make sure start earest wire and start time in range of 3rd cluster
    if ( wireStart3rdPlane == -1 )
      return -1;

    if ( (wireStart3rdPlane <= maxWire3rdClus) and (wireStart3rdPlane >= minWire3rdClus) and
         (StartT3rdPlane <= maxTime3rdClus) and (StartT3rdPlane >= minTime3rdClus) ){
      if (_verbose) { std::cout << "Match OK. Return: " << bestcompatVal << std::endl; }
      return 1./bestcompatVal;
    }

      if (_verbose) { std::cout << "\n\n***************************\n" << std::endl; }

    return -1;
  }

void cmtool::CFAlgoStartTimeCompat::Report ( )

virtual

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 193 of file CFAlgoStartTimeCompat.cxx.

  {

  }

void cmtool::CFAlgoStartTimeCompat::Reset ( void  )

virtual

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

Reimplemented from cmtool::CMAlgoBase.

Definition at line 21 of file CFAlgoStartTimeCompat.cxx.

  {
  }

void cmtool::CFAlgoStartTimeCompat::SetTimeDist ( double  t )

inline

Function to set verbose output.

Definition at line 60 of file CFAlgoStartTimeCompat.h.

{ _timeDist = t; }

void cmtool::CFAlgoStartTimeCompat::SetVerbose ( bool  on )

inlinevirtual

Function to set verbose output.

Reimplemented from cmtool::CMAlgoBase.

Definition at line 57 of file CFAlgoStartTimeCompat.h.

{ _verbose = on; }

Member Data Documentation

double cmtool::CFAlgoStartTimeCompat::_t2cm

private

Definition at line 64 of file CFAlgoStartTimeCompat.h.

double cmtool::CFAlgoStartTimeCompat::_timeDist

private

Definition at line 66 of file CFAlgoStartTimeCompat.h.

bool cmtool::CFAlgoStartTimeCompat::_verbose

private

Definition at line 65 of file CFAlgoStartTimeCompat.h.

double cmtool::CFAlgoStartTimeCompat::_w2cm

private

Definition at line 64 of file CFAlgoStartTimeCompat.h.

---

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

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