LArSoft libraries. More...

Classes
class	PedAlgoFixed
	Pedestal "algorithm" reading the pedestals from somewhere else. More...

class	FlashAlgoBase

class	FlashAlgoFactoryBase

class	FlashAlgoFactory
	Factory class for instantiating flash algorithm instance. More...

class	FlashFinderManager

struct	LiteOpHit_t

struct	LiteOpFlash_t

class	PECalib

class	SimpleFlashAlgo

class	SimpleFlashAlgoFactory
	A concrete factory class for pmtana::SimpleFlashAlgo. More...

class	AlgoCFD

class	AlgoFixedWindow

class	AlgoSiPM

class	AlgoSlidingWindow

class	AlgoThreshold

class	OpticalRecoException

class	PedAlgoEdges

class	PedAlgoRmsSlider

class	PedAlgoRollingMean

class	PedAlgoUB

class	PMTPedestalBase

struct	pulse_param

class	PMTPulseRecoBase

class	PulseRecoManager

class	RiseTimeCalculatorBase

class	RiseTimeThreshold

class	PMTAna

Typedefs
typedef fhicl::ParameterSet	Config_t

typedef std::vector < pmtana::LiteOpHit_t >	LiteOpHitArray_t

typedef std::vector < pmtana::LiteOpFlash_t >	LiteOpFlashArray_t

typedef std::vector< short >	Waveform_t

typedef std::vector< double >	PedestalMean_t

typedef std::vector< double >	PedestalSigma_t

typedef std::vector < pmtana::pulse_param >	pulse_param_array

Functions
std::vector< size_t >	ListOpChannels (int cryostat)

std::vector< size_t >	ListOpDets (int cryostat)

size_t	NOpDets (int cryostat)

size_t	OpDetFromOpChannel (size_t opch)

void	OpDetCenterFromOpChannel (size_t opch, double *xyz)

bool	CheckIndex (const std::vector< short > &wf, const size_t &begin, size_t &end)

double	mean (const std::vector< short > &wf, size_t start, size_t nsample)

double	edge_aware_mean (const std::vector< short > &wf, int start, int end)

double	std (const std::vector< short > &wf, const double ped_mean, size_t start, size_t nsample)

double	BinnedMaxOccurrence (const PedestalMean_t &mean_v, const size_t nbins)

int	sign (double val)

double	BinnedMaxTH1D (const std::vector< double > &v, int bins)

Variables
static const double	kINVALID_TIME = std::numeric_limits<double>::max()

static const unsigned int	kINVALID_CHANNEL = std::numeric_limits<unsigned int>::max()

static SimpleFlashAlgoFactory	__SimpleFlashAlgoFactoryStaticObject__

Detailed Description

LArSoft libraries.

C/C++ standard libraries

Typedef Documentation

typedef fhicl::ParameterSet pmtana::Config_t

Definition at line 11 of file icaruscode/icaruscode/PMT/OpReco/FlashFinder/FlashFinderFMWKInterface.h.

typedef std::vector<pmtana::LiteOpFlash_t> pmtana::LiteOpFlashArray_t

Definition at line 51 of file icaruscode/icaruscode/PMT/OpReco/FlashFinder/FlashFinderTypes.h.

typedef std::vector<pmtana::LiteOpHit_t> pmtana::LiteOpHitArray_t

Definition at line 24 of file icaruscode/icaruscode/PMT/OpReco/FlashFinder/FlashFinderTypes.h.

typedef std::vector< double > pmtana::PedestalMean_t

Definition at line 9 of file OpHitFinder/OpticalRecoTypes.h.

typedef std::vector< double > pmtana::PedestalSigma_t

Definition at line 10 of file OpHitFinder/OpticalRecoTypes.h.

typedef std::vector<pmtana::pulse_param> pmtana::pulse_param_array

Definition at line 53 of file PMTPulseRecoBase.h.

typedef std::vector< short > pmtana::Waveform_t

Definition at line 8 of file OpHitFinder/OpticalRecoTypes.h.

Function Documentation

double pmtana::BinnedMaxOccurrence	(	const PedestalMean_t &	mean_v,
		const size_t	nbins
	)

Definition at line 59 of file UtilFunc.cxx.

   {
     if(nbins<1) throw OpticalRecoException("Cannot have 0 binning");
 
     auto res = std::minmax_element(std::begin(mean_v),std::end(mean_v));
 
     double bin_width = ((*res.second) - (*res.first)) / ((double)nbins);
 
     if(nbins==1 || bin_width == 0) return ((*res.first) + bin_width /2.);
 
     //std::cout<<"Min: "<<(*res.first)<<" Max: "<<(*res.second)<<" Width: "<<bin_width<<std::endl;
 
     // Construct array of nbins
     static std::vector<size_t> ctr_v(nbins,0);
     for(auto& v : ctr_v) v=0;
     for(auto const& v : mean_v) {
 
       size_t index = int((v - (*res.first))/bin_width);
       //std::cout<<"adc = "<<v<<" width = "<<bin_width<< " ... "
       //<<index<<" / "<<ctr_v.size()<<std::endl;
 
       ctr_v[index]++;
 
     }
 
     // Find max occurrence
     auto max_it = std::max_element(std::begin(ctr_v),std::end(ctr_v));
 
     // Get the mean of max-occurrence bins
     double mean_max_occurrence = 0;
     double num_occurrence = 0;
     for(size_t bin=0; bin<ctr_v.size(); ++bin) {
 
       if(ctr_v[bin] != (*max_it)) continue;
 
       mean_max_occurrence += ((*res.first) + bin_width / 2. + bin_width * bin);
 
       num_occurrence += 1.0;
     }
 
     return (mean_max_occurrence / num_occurrence);
   }

double pmtana::BinnedMaxTH1D	(	const std::vector< double > &	v,
		int	bins
	)

Definition at line 112 of file UtilFunc.cxx.

                                                             {
 
     auto max_it = std::max_element(std::begin(v), std::end(v));
     auto min_it = std::min_element(std::begin(v), std::end(v));
 
     TH1D th("th",";;",bins,*min_it,*max_it);
 
     for (const auto & m : v) th.Fill(m);
 
     return th.GetXaxis()->GetBinCenter(th.GetMaximumBin());
   }

bool pmtana::CheckIndex	(	const std::vector< short > &	wf,
		const size_t &	begin,
		size_t &	end
	)

Definition at line 41 of file PMTPulseRecoBase.cxx.

   {
     if(begin >= wf.size() || end >= wf.size() || begin > end){
 
       std::cerr <<"Invalid arguments: waveform length = " << wf.size() << " begin = " << begin << " end = " << end << std::endl;
 
       return false;
     }
 
     if(!end) end = wf.size() - 1;
 
     return true;
   }

double pmtana::edge_aware_mean	(	const std::vector< short > &	wf,
		int	start,
		int	end
	)

Definition at line 25 of file UtilFunc.cxx.

                                                                          {
 
     auto m   = double{0.0};
     auto n_t = unsigned{0};
 
     for(int k = start; k < end; ++k) {
       if (k < 0 or k > (int)(wf.size()) - 1) continue;
       m += wf.at(k);
       ++n_t;
     }
 
     if( n_t > 0 ) m /= n_t;
     n_t = 0;
 
     return m;
   }

std::vector< size_t > pmtana::ListOpChannels ( int cryostat )

Definition at line 21 of file icaruscode/icaruscode/PMT/OpReco/FlashFinder/FlashFinderFMWKInterface.cxx.

                                                  {
     std::vector<size_t> res;
     ::art::ServiceHandle<geo::Geometry> geo;
     if(cryostat<0) {
       for(size_t opch=0; opch<geo->MaxOpChannel(); ++opch) {
         if(geo->IsValidOpChannel(opch)) continue;
         res.push_back(opch);
       }
     }else{
       auto const& bbox = geo->Cryostat(cryostat).Boundaries();
       for(size_t opch=0; opch<geo->MaxOpChannel(); ++opch) {
         if(geo->IsValidOpChannel(opch)) continue;
         auto const& pt = geo->OpDetGeoFromOpChannel(opch).GetCenter();
         if(!bbox.ContainsPosition(pt)) continue;
         res.push_back(opch);
       }
     }
     return res;
   }

std::vector< size_t > pmtana::ListOpDets ( int cryostat )

Definition at line 41 of file icaruscode/icaruscode/PMT/OpReco/FlashFinder/FlashFinderFMWKInterface.cxx.

                                              {
     std::vector<size_t> res;
     ::art::ServiceHandle<geo::Geometry> geo;
     if(cryostat<0) {
       for(size_t opdet=0; opdet<geo->NOpDets(); ++opdet) {
         res.push_back(opdet);
       }
     }else{
       auto const& bbox = geo->Cryostat(cryostat).Boundaries();
       for(size_t opdet=0; opdet<geo->NOpDets(); ++opdet) {
         auto const& pt = geo->OpDetGeoFromOpDet(opdet).GetCenter();
         if(!bbox.ContainsPosition(pt)) continue;
         res.push_back(opdet);
       }
     }
     return res;
   }

double pmtana::mean	(	const std::vector< short > &	wf,
		size_t	start,
		size_t	nsample
	)

Definition at line 13 of file UtilFunc.cxx.

   {
     if(!nsample) nsample = wf.size();
     if(start > wf.size() || (start+nsample) > wf.size())
       throw OpticalRecoException("Invalid start/end index!");
 
     double sum = std::accumulate(wf.begin()+start,wf.begin()+start+nsample,0.0) / ((double)nsample);
 
     return sum;
   }

size_t pmtana::NOpDets ( int cryostat )

Definition at line 59 of file icaruscode/icaruscode/PMT/OpReco/FlashFinder/FlashFinderFMWKInterface.cxx.

                                {
     ::art::ServiceHandle<geo::Geometry> geo;
     if(cryostat<0)
       return geo->NOpDets();
     else
       return geo->Cryostat(cryostat).NOpDet();
   }

void pmtana::OpDetCenterFromOpChannel	(	size_t	opch,
		double *	xyz
	)

Definition at line 72 of file icaruscode/icaruscode/PMT/OpReco/FlashFinder/FlashFinderFMWKInterface.cxx.

                                                           {
     ::art::ServiceHandle<geo::Geometry> geo;
     geo->OpDetGeoFromOpChannel(opch).GetCenter(xyz); 
   }

size_t pmtana::OpDetFromOpChannel ( size_t opch )

Definition at line 67 of file icaruscode/icaruscode/PMT/OpReco/FlashFinder/FlashFinderFMWKInterface.cxx.

                                          {
     ::art::ServiceHandle<geo::Geometry> geo;
     return geo->OpDetFromOpChannel(opch);
   }

int pmtana::sign ( double val )

Definition at line 104 of file UtilFunc.cxx.

                        {
 
     if (val > 0) return  1;
     if (val < 0) return -1;
     return 0;
 
   }

double pmtana::std	(	const std::vector< short > &	wf,
		const double	ped_mean,
		size_t	start,
		size_t	nsample
	)

Definition at line 42 of file UtilFunc.cxx.

   {
     if(!nsample) nsample = wf.size();
     if(start > wf.size() || (start+nsample) > wf.size())
       throw OpticalRecoException("Invalid start/end index!");
 
     double sigma = 0;
 
     for(size_t index=start; index < (start+nsample); ++index)
 
       sigma += pow( (wf[index] - ped_mean), 2 );
 
     sigma = sqrt(sigma/((double)(nsample)));
 
     return sigma;
   }

Variable Documentation

SimpleFlashAlgoFactory pmtana::__SimpleFlashAlgoFactoryStaticObject__

static

Definition at line 8 of file icaruscode/icaruscode/PMT/OpReco/FlashFinder/SimpleFlashAlgo.cxx.

const unsigned int pmtana::kINVALID_CHANNEL = std::numeric_limits<unsigned int>::max()

static

Definition at line 12 of file icaruscode/icaruscode/PMT/OpReco/FlashFinder/FlashFinderTypes.h.

const double pmtana::kINVALID_TIME = std::numeric_limits<double>::max()

static

Definition at line 11 of file icaruscode/icaruscode/PMT/OpReco/FlashFinder/FlashFinderTypes.h.

Classes

Typedefs

Functions

Variables

Detailed Description

Typedef Documentation

Function Documentation

Variable Documentation