cmtool::CFAlgoTimeProf Class Reference

#include "CFAlgoTimeProf.h"

Inheritance diagram for cmtool::CFAlgoTimeProf:

Public Member Functions
	CFAlgoTimeProf ()
	Default constructor. More...
virtual	~CFAlgoTimeProf ()
	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...
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...
virtual void	SetVerbose (bool doit=true)
	Setter function for verbosity. More...

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

Private Member Functions
float	TProfCompare (std::vector< util::PxHit > hita, std::vector< util::PxHit > hitb)

Detailed Description

User implementation for CFloatAlgoBase class doxygen documentation!

Definition at line 25 of file CFAlgoTimeProf.h.

Constructor & Destructor Documentation

cmtool::CFAlgoTimeProf::CFAlgoTimeProf

(

)

Default constructor.

Definition at line 20 of file CFAlgoTimeProf.cxx.

                                 : CFloatAlgoBase()
  //-------------------------------------------------------
  {

  }

cmtool::CFAlgoTimeProf::~CFAlgoTimeProf ( )

virtual

Default destructor.

Definition at line 27 of file CFAlgoTimeProf.cxx.

  {
  }

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::CFAlgoTimeProf::Float ( const std::vector< const cluster::ClusterParamsAlg * > &  clusters )

virtual

Core function: given the CPAN input, return whether a cluster should be merged or not.

Reimplemented from cmtool::CFloatAlgoBase.

Definition at line 40 of file CFAlgoTimeProf.cxx.

References TProfCompare().

  {
  // We now have a vector a clusters.

  //### need a pointer to the cluster just return -1  
    for(auto const& ptr : clusters) if(!ptr) return -1;

        std::vector<util::PxHit> hits0;
        std::vector<util::PxHit> hits1;
        std::vector<util::PxHit> hits2;

        // loop over the clusters
        for(auto const& c : clusters)
        {
//      std::cout<<"Size of the xluster vector"<<clusters.size()<<std::endl;
//      auto Plane = c->Plane();
//      auto StartPoint = c->GetParams().start_point.t;
//      auto EndPoint = c->GetParams().end_point.t;
//      std::cout<<"\t RG Cluster info:\n \t plane: "<<Plane<<std::endl;
//      std::cout<<"\tStart point: "<<StartPoint<<std::endl;
//        std::cout<<"\tEnd Point: "<<EndPoint <<std::endl;

        // Get assosiations for this cluster
        if(c->Plane() ==0) hits0 = c->GetHitVector();
        if(c->Plane() ==1) hits1 = c->GetHitVector();
        if(c->Plane() ==2) hits2 = c->GetHitVector();

        }// for over the clusters

//      std::cout<<"Looking for the hits vector size"<<hits0.size()<<","<<hits1.size()<<","<<hits2.size()<<std::endl;
//        std::cout<<"############# End of loop############# "<<std::endl;
        // make an integrale over the cluster
        //bool pl0 = false;
        //bool pl1 = false;
        //bool pl2 = false;
        float tprof01 = -1;
        float tprof02 = -1;
        float tprof12 = -1;
        if(hits0.size()>0)
        {
          //pl0 = true;
          if(hits1.size()>0)
            {
              //pl1 = true;
              //we need to do a profile
              tprof01 =TProfCompare(hits0,hits1);
            }// hits1size>0
          if(hits2.size()>0)
            {
              //pl2 = true;
              //we need to do a profile
              tprof02 =TProfCompare(hits0,hits2);
                }// hits2size>0
        }// hits0size >0
        
        if(hits1.size()>0)
          {
            //pl1 = true;
            if(hits2.size()>0)
              {
                //pl2 = true;
                //we need to do a profile
                tprof12 =TProfCompare(hits1,hits2);
              }//hits2.size>0
          }// hits1size>0
        
        // This is going to be slow is we are having to re-calutlate this ever time. For now it will have to do
        
        //      std::cout<< " \t summary of timeprof Planes that are accepted :" <<pl0<<" | " <<pl1<<" | " <<pl2<<" |"<<std::endl;
        std::vector<float> tprofmatches;        
//      std::cout<< " \t                  Value of timeprof(01,02,12) :" <<tprof01<<" | " <<tprof02<<" | " <<tprof12<<" |"<<std::endl;
        tprofmatches.push_back(tprof01);
        tprofmatches.push_back(tprof02);
        tprofmatches.push_back(tprof12);
        
        float matchscore=0;
        float avgcounter=0;
//      std::cout<<"SIZE of trpfmoatch"<< tprofmatches.size()<<std::endl;
        for( unsigned int a=0;a<tprofmatches.size();a++)
        {
        if(tprofmatches[a]==-1) continue;       
        else {
                matchscore +=tprofmatches[a];
                avgcounter +=1;
                }// end of else
        }//for over the tprofmatchs
                if(avgcounter!=0){
//              std::cout << " Match Score pree "<< matchscore<<std::endl;              
//              std::cout<< " Counter "<< avgcounter;
                matchscore /= avgcounter;
//              std::cout << " Match Score "<< matchscore<<std::endl;           
                }
                else{
//              std::cout << " Match Score "<< -1<<std::endl;           
                 return -1;
                }
                return matchscore;
                
        
                
  
    //if(clusters.size()) return 1.;
    //else 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::CFAlgoTimeProf::Report ( )

virtual

Optional function: called after Bool() function is called for all possible cluster pairs by CMergeManager/CMatchManager IFF run with verbosity level kPerIteration. Maybe useful for debugging.

Reimplemented from cmtool::CMAlgoBase.

Definition at line 221 of file CFAlgoTimeProf.cxx.

  {

  }

void cmtool::CFAlgoTimeProf::Reset ( )

virtual

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

Reimplemented from cmtool::CMAlgoBase.

Definition at line 33 of file CFAlgoTimeProf.cxx.

  {

  }

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

virtual void cmtool::CMAlgoBase::SetVerbose ( bool  doit = true )

inlinevirtualinherited

Setter function for verbosity.

Reimplemented in cmtool::CBAlgoArray, cmtool::CFAlgoChargeDistrib, cmtool::CFAlgoVolumeOverlap, cmtool::CFAlgoWireOverlap, cmtool::CFAlgoStartPointCompat, cmtool::CFAlgoStartTimeCompat, cmtool::CFAlgoTimeOverlap, cmtool::CBAlgoStartInCone, cmtool::CBAlgoShortestDistNonEndpoint, and cmtool::CBAlgoTrackSeparate.

Definition at line 80 of file CMAlgoBase.h.

References cmtool::CMAlgoBase::_verbose.

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

{ _verbose = doit; }

float cmtool::CFAlgoTimeProf::TProfCompare ( std::vector< util::PxHit >  hita, std::vector< util::PxHit >  hitb  )

private

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 147 of file CFAlgoTimeProf.cxx.

References bin, detinfo::DetectorProperties::GetXTicksOffset(), detinfo::DetectorProperties::NumberTimeSamples(), and util::GeometryUtilities::TimeToCm().

Referenced by Float().

 {
   ::util::GeometryUtilities geou;
   const detinfo::DetectorProperties* detp = lar::providerFrom<detinfo::DetectorPropertiesService>();
//   int nts = larutil::DetectorProperties::GetME()->NumberTimeSamples()*larutil::GeometryUtilities::GetME()->TimeToCm();
   // Where is this?
   //int nplanes = geom->Nplanes();
   int nplanes = 3;
   double ks = 0.0;
   std::vector< std::vector<TH1D*> > signals(nplanes);
   std::vector< std::vector<TH1D*> > pulses(nplanes);
   
   double time_diff = ( detp->GetXTicksOffset((int)(hita.at(0).plane),0,0) - 
                        detp->GetXTicksOffset((int)(hitb.at(0).plane),0,0) ) * geou.TimeToCm();
   
   // First go look for the min & max of hits 
   double min_time = detp->NumberTimeSamples() * geou.TimeToCm();
   double max_time = 0;
   for(auto const& h : hita) {
     if(h.t > max_time) max_time = h.t;
     if(h.t < min_time) min_time = h.t;
   }
   for(auto const& h : hitb) {
     if( (h.t + time_diff) > max_time ) max_time = h.t + time_diff;
     if( (h.t + time_diff) < min_time ) min_time = h.t + time_diff;
   }

   TH1D histo_a("ha", "", 200, min_time-1, max_time+1);
   TH1D histo_b("hb", "", 200, min_time-1, max_time+1);
   TH1D histo_inta("hinta", "", 200, min_time-1, max_time+1);
   TH1D histo_intb("hintb", "", 200, min_time-1, max_time+1);

   // First loop over hits in A and make the hist
   // in this case let's just use plane 0,1
   for( auto const& ha : hita){
     double time = ha.t;
     //time -= larutil::DetectorProperties::GetME()->GetXTicksOffset(ha.plane)*larutil::GeometryUtilities::GetME()->TimeToCm();
     double charge = ha.charge;

     int bin = histo_a.FindBin(time);
     histo_a.SetBinContent(bin,histo_a.GetBinContent(bin)+charge);
     for (int j = bin; j<=histo_a.GetNbinsX(); ++j){
       histo_inta.SetBinContent(j,histo_inta.GetBinContent(j)+charge);
     }
   }
   if (histo_inta.Integral()) histo_inta.Scale(1./histo_inta.GetBinContent(histo_inta.GetNbinsX()));
   for( auto const& hb : hitb){
     double time = hb.t + time_diff;
     //time -= larutil::DetectorProperties::GetME()->GetXTicksOffset(hb.plane)*larutil::GeometryUtilities::GetME()->TimeToCm();
     double charge = hb.charge;
     int bin = histo_b.FindBin(time);
     histo_b.SetBinContent(bin,histo_b.GetBinContent(bin)+charge);
     for (int j = bin; j<=histo_b.GetNbinsX(); ++j){
       histo_intb.SetBinContent(j,histo_intb.GetBinContent(j)+charge);
     }
   }

/*
   std::cout 
     << "siga: "<< histo_a.GetEntries() << std::endl
     << "sigb: "<< histo_b.GetEntries() << std::endl
     << "siginta: "<<histo_inta.GetEntries() << std::endl
     << "sigintb: "<<histo_intb.GetEntries() << std::endl
     << std::endl;
*/

   if (histo_intb.Integral()) histo_intb.Scale(1./histo_intb.GetBinContent(histo_intb.GetNbinsX()));
   ks = histo_inta.KolmogorovTest(&histo_intb);

   std::cout<<ks<<std::endl;
   return ks;
 }

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

bool cmtool::CMAlgoBase::_verbose

protectedinherited

Boolean to choose verbose mode. Turned on if CMergeManager/CMatchManager's verbosity level is >= kPerMerging.

Definition at line 88 of file CMAlgoBase.h.

Referenced by cmtool::CBAlgoAngleAlign::Bool(), cmtool::CBAlgoAngleCompat::Bool(), cmtool::CBAlgoShortestDistSmallCluster::Bool(), cmtool::CBAlgoShortestDist::Bool(), cmtool::CBAlgoOutOfConeSeparate::Bool(), cmtool::CBAlgoPolyOverlap::Bool(), cmtool::CBAlgoAngleSeparate::Bool(), cmtool::CBAlgoCenterOfMassSmall::Bool(), cmtool::CBAlgoCenterOfMass::Bool(), cmtool::CBAlgoStartInPoly::Bool(), cmtool::CBAlgoStartNearEnd::Bool(), cmtool::CBAlgoProhibitAllTracks::Bool(), cmtool::CBAlgoShortestDist::CBAlgoShortestDist(), cmtool::CBAlgoShortestDistSmallCluster::CBAlgoShortestDistSmallCluster(), cmtool::CMAlgoBase::CMAlgoBase(), cmtool::CFAlgoZOverlap::Float(), cmtool::CFAlgoQRatio::Float(), cmtool::CFAlgoStartPointMatch::Float(), cmtool::CMAlgoBase::SetVerbose(), cmtool::CBAlgoShortestDist::ShortestDistanceSquared(), cmtool::CBAlgoShortestDistSmallCluster::ShortestDistanceSquared(), cmtool::CBAlgoCenterOfMass::ShortestDistanceSquared(), and cmtool::CBAlgoCenterOfMassSmall::ShortestDistanceSquared().

---

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

• /cvmfs/larsoft.opensciencegrid.org/products/larreco/v07_10_02/source/larreco/RecoAlg/CMTool/CMTAlgMatch/CFAlgoTimeProf.h
• /cvmfs/larsoft.opensciencegrid.org/products/larreco/v07_10_02/source/larreco/RecoAlg/CMTool/CMTAlgMatch/CFAlgoTimeProf.cxx