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...
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...

Protected Attributes
TFile *	_fout
	TFile pointer to an output file. More...

Private Attributes
double	_w2cm
double	_t2cm
bool	_verbose
double	_timeDist

Detailed Description

User implementation for CFloatAlgoBase class doxygen documentation!

Definition at line 26 of file CFAlgoStartTimeCompat.h.

Constructor & Destructor Documentation

cmtool::CFAlgoStartTimeCompat::CFAlgoStartTimeCompat

(

)

Default constructor.

Definition at line 10 of file CFAlgoStartTimeCompat.cxx.

References _t2cm, _w2cm, SetTimeDist(), SetVerbose(), util::GeometryUtilities::TimeToCm(), and util::GeometryUtilities::WireToCm().

                                               : 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 34 of file CFAlgoStartTimeCompat.h.

References Float(), Report(), and Reset().

{};

Member Function Documentation

virtual void cmtool::CMAlgoBase::EventBegin ( const std::vector< cluster::ClusterParamsAlg > &  clusters )

inlinevirtualinherited

Optional function: called at the beginning of 1st iteration. This is called per event.

Reimplemented in cmtool::CFAlgoArray, cmtool::CPAlgoArray, cmtool::CBAlgoArray, and cmtool::CBAlgoPolyShortestDist.

Definition at line 45 of file CMAlgoBase.h.

Referenced by cmtool::CMergeManager::EventBegin().

    { if(clusters.size()) return; }

virtual void cmtool::CMAlgoBase::EventEnd ( )

inlinevirtualinherited

Optional function: called at the end of event ... after the last merging iteration is over.

Reimplemented in cmtool::CFAlgoArray, cmtool::CPAlgoArray, and cmtool::CBAlgoArray.

Definition at line 51 of file CMAlgoBase.h.

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

    {return;}

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.

References _timeDist, _verbose, _w2cm, geo::GeometryCore::ChannelsIntersect(), e, hits(), geo::GeometryCore::NearestWire(), and geo::GeometryCore::PlaneWireToChannel().

Referenced by ~CFAlgoStartTimeCompat().

  {

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

virtual void cmtool::CMAlgoBase::IterationBegin ( const std::vector< cluster::ClusterParamsAlg > &  clusters )

inlinevirtualinherited

Optional function: called at the beggining of each iteration over all pairs of clusters. 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.

Reimplemented in cmtool::CFAlgoArray, cmtool::CPAlgoArray, and cmtool::CBAlgoArray.

Definition at line 59 of file CMAlgoBase.h.

Referenced by cmtool::CMatchManager::EventBegin(), cmtool::CMatchManager::IterationBegin(), and cmtool::CMergeManager::IterationBegin().

    { if(clusters.size()) return;}

virtual void cmtool::CMAlgoBase::IterationEnd ( )

inlinevirtualinherited

Optional function: called at the end of each iteration over all pairs of clusters.

Reimplemented in cmtool::CFAlgoArray, cmtool::CPAlgoArray, and cmtool::CBAlgoArray.

Definition at line 65 of file CMAlgoBase.h.

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

    {return; }

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.

Referenced by ~CFAlgoStartTimeCompat().

  {

  }

void cmtool::CFAlgoStartTimeCompat::Reset ( )

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.

Referenced by ~CFAlgoStartTimeCompat().

  {
  }

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

inlineinherited

Setter function for an output plot TFile pointer.

Definition at line 77 of file CMAlgoBase.h.

References cmtool::CMAlgoBase::_fout.

Referenced by cmtool::CMergeManager::EventBegin().

{ _fout = fout; }

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

inline

Function to set verbose output.

Definition at line 61 of file CFAlgoStartTimeCompat.h.

References _timeDist.

Referenced by CFAlgoStartTimeCompat().

{ _timeDist = t; }

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

inlinevirtual

Function to set verbose output.

Reimplemented from cmtool::CMAlgoBase.

Definition at line 58 of file CFAlgoStartTimeCompat.h.

References _verbose.

Referenced by CFAlgoStartTimeCompat().

{ _verbose = on; }

Member Data Documentation

TFile* cmtool::CMAlgoBase::_fout

protectedinherited

TFile pointer to an output file.

Definition at line 85 of file CMAlgoBase.h.

Referenced by cmtool::CMAlgoBase::CMAlgoBase(), and cmtool::CMAlgoBase::SetAnaFile().

double cmtool::CFAlgoStartTimeCompat::_t2cm

private

Definition at line 65 of file CFAlgoStartTimeCompat.h.

Referenced by CFAlgoStartTimeCompat().

double cmtool::CFAlgoStartTimeCompat::_timeDist

private

Definition at line 67 of file CFAlgoStartTimeCompat.h.

Referenced by Float(), and SetTimeDist().

bool cmtool::CFAlgoStartTimeCompat::_verbose

private

Definition at line 66 of file CFAlgoStartTimeCompat.h.

Referenced by Float(), and SetVerbose().

double cmtool::CFAlgoStartTimeCompat::_w2cm

private

Definition at line 65 of file CFAlgoStartTimeCompat.h.

Referenced by CFAlgoStartTimeCompat(), and Float().

---

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

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