#include <CommonAmps.h>

Inheritance diagram for flashmatch::CommonAmps:

Public Member Functions
	CommonAmps (const std::string name="CommonAmps")
	Default constructor. More...

	~CommonAmps ()
	Default destructor. More...

FlashMatch_t	Match (const QCluster_t &, const Flash_t &)

Public Member Functions inherited from flashmatch::BaseFlashMatch
	BaseFlashMatch (const std::string name="noname")
	Default constructor. More...

virtual	~BaseFlashMatch ()
	Default destructor. More...

Flash_t	GetEstimate (const QCluster_t &) const
	Method to call flash hypothesis. More...

void	FillEstimate (const QCluster_t &, Flash_t &) const
	Method to simply fill provided reference of flashmatch::Flash_t. More...

virtual void	SetTPCCryo (int tpc, int cryo)=0
	Sets the TPC and Cryo numbers. More...

Public Member Functions inherited from flashmatch::BaseAlgorithm
	BaseAlgorithm (const Algorithm_t type, const std::string name)
	Default constructor. More...

	~BaseAlgorithm ()
	Default destructor. More...

void	Configure (const Config_t &pset)
	Function to accept configuration. More...

Algorithm_t	AlgorithmType () const
	Algorithm type. More...

const std::string &	AlgorithmName () const
	Algorithm name. More...

Public Member Functions inherited from flashmatch::LoggerFeature
	LoggerFeature (const std::string logger_name="LoggerFeature")
	Default constructor. More...

	LoggerFeature (const LoggerFeature &original)
	Default copy constructor. More...

virtual	~LoggerFeature ()
	Default destructor. More...

const flashmatch::logger &	logger () const
	Logger getter. More...

void	set_verbosity (::flashmatch::msg::Level_t level)
	Verbosity level. More...

const std::string &	name () const
	Name getter, defined in a logger instance attribute. More...

Protected Member Functions
void	_Configure_ (const Config_t &pset)

Private Attributes
float	_percent

float	_score

float	_x_step_size

flashmatch::QCluster_t	_tpc_qcluster

flashmatch::Flash_t	_vis_array

Additional Inherited Members
Protected Attributes inherited from flashmatch::BaseFlashMatch
int	_tpc = 0
	The TPC number to use. More...

int	_cryo = 0
	The Cryostat number to use. More...

Detailed Description

Implementation of flashmatch::BaseFlashHypothesis algorithm class.
The goal of this algorithm is to compare the most prominent pieces of opflash spectra, and spectra produced by photon library given QCluster points. Tpc_point calculated as in QWeightPoint and stored based on the best match amplitudes.

Only works with Photon Library currently

Definition at line 32 of file CommonAmps.h.

Constructor & Destructor Documentation

flashmatch::CommonAmps::CommonAmps ( const std::string name = "CommonAmps" )

Default constructor.

Definition at line 13 of file CommonAmps.cxx.

     : BaseFlashMatch(name)
   {
     _percent = 0.5;
     _score   = 0.8;    
   }

flashmatch::CommonAmps::~CommonAmps ( )

inline

Default destructor.

Definition at line 40 of file CommonAmps.h.

40 {}

Member Function Documentation

void flashmatch::CommonAmps::_Configure_ ( const Config_t & pset )

protectedvirtual

Implements flashmatch::BaseAlgorithm.

Definition at line 20 of file CommonAmps.cxx.

   {
     _percent = pset.get<double>("QFracThreshold");
     _score   = pset.get<double>("ScoreThreshold");
     _x_step_size = pset.get<double>("XStepSize");
   }

FlashMatch_t flashmatch::CommonAmps::Match	(	const QCluster_t &	,
		const Flash_t &
	)

virtual

CORE FUNCTION: takes in flashmatch::QCluster_t (TPC object) and flashmatch::Flash_t (flash)
and inspect if two are consistent (i.t. matched) or not. Returns flashmatch::FlashMatch_t
which represents the compatibility of two inputs. In particular the algorithm needs to
set the "score" and "QPoint_t" of the return object. The former represents the goodness
of a match with a value larger than 0. Negative value is interpreted as no match.
The latter represents the matched TPC 3D point. The "tpc_id and "flash_id" of the return
object is set by flashmatch::FlashMatchManager, the caller of the algorithm, as it manages
the overall collection of user input flash and TPC objects.

Note it is flashmatch::FlashMatchManager configuration option to allow an assignment of the
same flash to multiple TPC object or not. If not allowed, a match with a higher "score"
in the return object is chosen.

Implements flashmatch::BaseFlashMatch.

Definition at line 27 of file CommonAmps.cxx.

   {
     
     double integral_op  = std::accumulate(std::begin(flash.pe_v),
                                           std::end(flash.pe_v),
                                           0.0);
     double maxRatio     = -1;
     //double maxX         = 0;
     
     if(_vis_array.pe_v.empty())
       _vis_array.pe_v.resize(OpDetXArray().size());
 
     // Create multimap to hold the largest amplitudes in the first slots of map
     // Normalize each PE bin to the total # of PEs-- 1/x to put highest amps in front
     std::multimap<double,int> ampToOpDet ;
 
     for(size_t k=0; k<32; k++)
       ampToOpDet.emplace(1./(flash.pe_v[k]/integral_op),k);
     
 
     // Now that opdet hits are normalized and ordered, store indices for PMTs which 
     // had the greatest hit amplitudes up to _percent of total PMT hits
     double opAmpTotal = 0;
     std::vector<int> ids ;
     ids.resize(0);
 
     for( auto const & e : ampToOpDet ){
       opAmpTotal += 1/e.first ;
       ids.push_back(e.second) ;
       if( opAmpTotal > _percent )
         break;
     }
 
     // Prepare the return values (Mostly QWeightPoint)
     FlashMatch_t f;
     if(pt_v.empty()){
       std::cout<<"Not enough points!"<<std::endl;
       return f;
     }
     
     _tpc_qcluster.resize(pt_v.size());
 
     // Get min & max x value
     double x_max = 0;
     double x_min = 1e12;
 
     for(auto const& pt : pt_v) {
       if(pt.x > x_max) x_max = pt.x;
       if(pt.x < x_min) x_min = pt.x;
     }
 
     for(double x_offset=0;
         x_offset<(250.-(x_max-x_min));
         x_offset+=_x_step_size) {
       
       // Create QCluster_t with this offset
 
       for(size_t i=0; i<_tpc_qcluster.size(); ++i) {
         _tpc_qcluster[i].x = pt_v[i].x + x_offset - x_min;
         _tpc_qcluster[i].y = pt_v[i].y;
         _tpc_qcluster[i].z = pt_v[i].z;
         _tpc_qcluster[i].q = pt_v[i].q;
       }
             
       FillEstimate(_tpc_qcluster,_vis_array);
 
       // Calculate amplitudes corresponding to max opdet amplitudes
       double visAmpTotal = 0;
       double vis_pe_sum = _vis_array.TotalPE();
 
       for(size_t i=0; i<ids.size(); i++) {
         if(_vis_array.pe_v[ids[i]]<0) continue;
         visAmpTotal += _vis_array.pe_v[ids[i]] / vis_pe_sum ;
       }
       
       double ratio = 0;
       if(opAmpTotal > visAmpTotal )
         ratio = visAmpTotal/opAmpTotal ;
       else
         ratio = opAmpTotal/visAmpTotal ;
       
       if(ratio > maxRatio) {
         maxRatio = ratio;
         //maxX     = x_offset;
 
         f.score = ratio;
         f.tpc_point.x = f.tpc_point.y = f.tpc_point.z = 0;
         f.tpc_point.q = vis_pe_sum;
 
         for(size_t pmt_index=0; pmt_index<NOpDets(); ++pmt_index) {
           double pe = _vis_array.pe_v[pmt_index];
           if(pe<0) continue;
           f.tpc_point.x += OpDetX(pmt_index) * pe / vis_pe_sum;
           f.tpc_point.y += OpDetY(pmt_index) * pe / vis_pe_sum;
           f.tpc_point.z += OpDetZ(pmt_index) * pe / vis_pe_sum;
         }
       }
     }
   
     // If min-diff is bigger than assigned max, return default match (score<0)
     if( maxRatio < _score ) {
 
       f.tpc_point.x = f.tpc_point.y = f.tpc_point.z = -1;
       f.tpc_point.q = -1;
       f.score = -1;
       return f;
     }
 
     f.hypothesis = _vis_array.pe_v;
     return f;
   }

Member Data Documentation

float flashmatch::CommonAmps::_percent

private

Definition at line 50 of file CommonAmps.h.

float flashmatch::CommonAmps::_score

private

Definition at line 51 of file CommonAmps.h.

flashmatch::QCluster_t flashmatch::CommonAmps::_tpc_qcluster

private

Definition at line 53 of file CommonAmps.h.

flashmatch::Flash_t flashmatch::CommonAmps::_vis_array

private

Definition at line 54 of file CommonAmps.h.

float flashmatch::CommonAmps::_x_step_size

private

Definition at line 52 of file CommonAmps.h.

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

Public Member Functions

Protected Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation