hough::IslandsTH2 Class Reference

#include "/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-03-06/HoughTransform/IslandsTH2.h"

Inheritance diagram for hough::IslandsTH2:

Public Member Functions
	IslandsTH2 (const char *n, const char *t, const TH2 *h, double q, bool wrapx, bool wrapy)
unsigned int	Nislands () const
void	LabelIslands (const TH2 *h)

Private Member Functions
void	Mark (const TH2 *h, int i, int j, int isl_indx, double th)
int	Neighbor (int i, int d, int imx, bool wrap)

Private Attributes
double	fPco
	Number of sigma above average peak height to use for peak threshold. More...
bool	fWrapX
	Wrap around the x boundaries? More...
bool	fWrapY
	Wrap around the y boundaries? More...
int	fNx
	Number of bins along x. More...
int	fNy
	Number of bins along y. More...
unsigned int	fNislands
	Number of islands found. More...

Detailed Description

Definition at line 12 of file IslandsTH2.h.

Constructor & Destructor Documentation

hough::IslandsTH2::IslandsTH2	(	const char *	n,
		const char *	t,
		const TH2 *	h,
		double	q,
		bool	wrapx,
		bool	wrapy
	)

Definition at line 13 of file IslandsTH2.cxx.

                                     :
    TH2S(n,
         t,
         h->GetNbinsX(),
         h->GetXaxis()->GetXmin(),
         h->GetXaxis()->GetXmax(),
         h->GetNbinsY(),
         h->GetYaxis()->GetXmin(),
         h->GetYaxis()->GetXmax()),
    fPco(q),
    fWrapX(wrapx),
    fWrapY(wrapy),
    fNx(h->GetNbinsX()),
    fNy(h->GetNbinsY()),
    fNislands(0)
  { }

Member Function Documentation

void hough::IslandsTH2::LabelIslands

(

const TH2 *

h

)

Definition at line 37 of file IslandsTH2.cxx.

References a, b, fNislands, fNx, fNy, fPco, MECModelEnuComparisons::i, calib::j, Mark(), genie::units::nb, and std::sqrt().

Referenced by hough::Hough2P::FindPeaks(), and Nislands().

  {
    
    int a;
    int b;
    double Pave = 0.0;     // average peak height
    double Ps = 0.0;       // standard deviation of average peak height
    double th = 0.0;       // threshold to use for peak cutoff
    int nb = 0;            // number of bins counted for avergae peak height
    
    // calculate average peak height to determine threshold
    // empty bins are NOT counted
    for (a = 1; a <= fNx; ++a) {
      for (b = 1; b <= fNy; ++b) {
        if (h->GetBinContent(a,b) != 0) {
          Pave += h->GetBinContent(a,b);
          Ps += (h->GetBinContent(a,b))*(h->GetBinContent(a,b));
          nb++;
        } // end if
      } // end for b
    } // end for a
    
    Pave = Pave/nb;
    Ps = sqrt(Ps/nb - Pave*Pave);
    
    th = Pave + fPco*Ps;
    
    
    int i;
    int j;
    for (i=1; i<=fNx; ++i) {
      for (j=1; j<=fNy; ++j) {
        if (this->GetBinContent(i,j) == 0 && h->GetBinContent(i,j) >= th)       {
          ++fNislands;
          this->Mark(h,i,j,fNislands,th);
        } // end if
      } // end for j
    } // end for i
    
  } // end LabelIslands

void hough::IslandsTH2::Mark ( const TH2 *  h, int  i, int  j, int  isl_indx, double  th  )

private

Definition at line 96 of file IslandsTH2.cxx.

References fNx, fNy, fWrapX, fWrapY, m, getGoodRuns4SAM::n, Neighbor(), and SetBinContent().

Referenced by LabelIslands(), and Nislands().

  {
    int m,  n;
    int ix, iy;
    
    this->SetBinContent(i,j,isl_indx);
    for (m=-1; m<=1; ++m) {
      ix = this->Neighbor(i,m,fNx,fWrapX);
      if (ix==-1) continue;
      for (n=-1; n<=1; ++n) {
        if (m==0 && n==0) continue;
        iy = this->Neighbor(j,n,fNy,fWrapY);
        if (iy==-1) continue;
        if (this->GetBinContent(ix,iy) == 0 && h->GetBinContent(ix,iy) >= th) {
          this->Mark(h,ix,iy,isl_indx,th);
        }
      }
    }
  }

int hough::IslandsTH2::Neighbor ( int  i, int  d, int  imx, bool  wrap  )

private

Definition at line 80 of file IslandsTH2.cxx.

References m, and getGoodRuns4SAM::n.

Referenced by Mark(), and Nislands().

  {
    int ix = i+m;
    if (ix==0) {
      if (wrap) return n;
      else      return -1;
    }
    if (ix>n) {
      if (wrap) return 1;
      else      return -1;
    }
    return ix;
  }

unsigned int hough::IslandsTH2::Nislands ( ) const

inline

Definition at line 20 of file IslandsTH2.h.

References d, fNislands, MECModelEnuComparisons::i, calib::j, LabelIslands(), Mark(), and Neighbor().

Referenced by hough::Hough2P::FindPeaks().

{ return fNislands; }

Member Data Documentation

unsigned int hough::IslandsTH2::fNislands

private

Number of islands found.

Definition at line 31 of file IslandsTH2.h.

Referenced by LabelIslands(), and Nislands().

int hough::IslandsTH2::fNx

private

Number of bins along x.

Definition at line 29 of file IslandsTH2.h.

Referenced by LabelIslands(), and Mark().

int hough::IslandsTH2::fNy

private

Number of bins along y.

Definition at line 30 of file IslandsTH2.h.

Referenced by LabelIslands(), and Mark().

double hough::IslandsTH2::fPco

private

Number of sigma above average peak height to use for peak threshold.

Definition at line 26 of file IslandsTH2.h.

Referenced by LabelIslands().

bool hough::IslandsTH2::fWrapX

private

Wrap around the x boundaries?

Definition at line 27 of file IslandsTH2.h.

Referenced by Mark().

bool hough::IslandsTH2::fWrapY

private

Wrap around the y boundaries?

Definition at line 28 of file IslandsTH2.h.

Referenced by Mark().

---

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

• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-03-06/HoughTransform/IslandsTH2.h
• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-03-06/HoughTransform/IslandsTH2.cxx