ems::DirOfGamma Class Reference

#include "DirOfGamma.h"

Public Member Functions
	DirOfGamma (const detinfo::DetectorPropertiesData &detProp, const std::vector< art::Ptr< recob::Hit >> &src, unsigned int nbins, unsigned int idcl)
	~DirOfGamma ()
TVector2 const &	GetBaryCenterCm () const
std::vector< Hit2D * > const &	GetHits2D () const
std::vector< EndPoint > const &	GetCandidates () const
void	SetIdCandidate (size_t id)
size_t	GetIdCandidate ()
art::Ptr< recob::Hit > const &	GetFirstHit () const
std::vector< art::Ptr< recob::Hit > > const &	GetHits ()
TVector2 const &	GetFirstPoint () const
std::vector< art::Ptr< recob::Hit > > const &	GetIniHits () const
size_t	GetIdCl () const

Private Member Functions
void	FindInitialPart ()
void	FillBins ()
bool	FindCandidates ()
void	ComputeBaryCenter ()
void	ComputeMaxDist ()
void	ComputeMaxCharge ()
void	ComputeFinalValues ()

Private Attributes
size_t	fNbins
size_t	fIdCl
size_t	fCandidateID
std::vector< Hit2D * >	fPoints2D
std::vector< Bin2D >	fBins
std::vector< EndPoint >	fCandidates
art::Ptr< recob::Hit >	fStartHit
TVector2	fStartPoint
std::vector< art::Ptr< recob::Hit > >	fIniHits
std::vector< art::Ptr< recob::Hit > >	fHits
TVector2	fBaryCenter
float	fNormDist
float	fNormCharge
bool	fIsCandidateIDset

Detailed Description

Definition at line 119 of file DirOfGamma.h.

Constructor & Destructor Documentation

ems::DirOfGamma::DirOfGamma	(	const detinfo::DetectorPropertiesData &	detProp,
		const std::vector< art::Ptr< recob::Hit >> &	src,
		unsigned int	nbins,
		unsigned int	idcl
	)

Definition at line 160 of file DirOfGamma.cxx.

References ComputeBaryCenter(), ComputeMaxCharge(), ComputeMaxDist(), fBaryCenter, fBins, fHits, FillBins(), FindCandidates(), FindInitialPart(), fNbins, and fPoints2D.

  : fNbins(nbins), fIdCl(idcl), fCandidateID(0), fIsCandidateIDset(false)
{
  fHits = src;

  for (unsigned int i = 0; i < src.size(); i++) {
    Hit2D* hit = new Hit2D(detProp, src[i]);
    fPoints2D.push_back(hit);
  }

  ComputeBaryCenter();

  for (unsigned int i = 0; i < fNbins; i++)
    fBins.push_back(Bin2D(fBaryCenter));

  FillBins();
  ComputeMaxDist();
  if (FindCandidates()) {
    ComputeMaxCharge();
    FindInitialPart();
  }
}

ems::DirOfGamma::~DirOfGamma ( )

inline

Definition at line 125 of file DirOfGamma.h.

  {
    for (unsigned int i = 0; i < fPoints2D.size(); ++i)
      delete fPoints2D[i];
  }

Member Function Documentation

void ems::DirOfGamma::ComputeBaryCenter ( )

private

Definition at line 186 of file DirOfGamma.cxx.

References fBaryCenter, and fPoints2D.

Referenced by DirOfGamma().

{
  double nomx = 0.0;
  double nomy = 0.0;
  double denom = 0.0;
  for (unsigned int i = 0; i < fPoints2D.size(); i++) {
    nomx += fPoints2D[i]->GetPointCm().X() * fPoints2D[i]->GetCharge();
    nomy += fPoints2D[i]->GetPointCm().Y() * fPoints2D[i]->GetCharge();
    denom += fPoints2D[i]->GetCharge();
  }

  double bx = nomx / denom;
  double by = nomy / denom;
  fBaryCenter.Set(bx, by);
}

void ems::DirOfGamma::ComputeFinalValues ( )

private

void ems::DirOfGamma::ComputeMaxCharge ( )

private

Definition at line 276 of file DirOfGamma.cxx.

References fCandidates, and fNormCharge.

Referenced by DirOfGamma().

{
  fNormCharge = 0.0;
  for (unsigned int i = 0; i < fCandidates.size(); i++) {
    if (fCandidates[i].GetMaxCharge() > fNormCharge) {
      fNormCharge = fCandidates[i].GetMaxCharge();
    }
  }
}

void ems::DirOfGamma::ComputeMaxDist ( )

private

Definition at line 229 of file DirOfGamma.cxx.

References pma::Dist2(), fBaryCenter, fBins, fNormDist, and ems::Hit2D::GetPointCm().

Referenced by DirOfGamma().

{
  double maxdist2 = 0.0;

  for (unsigned int id = 0; id < fBins.size(); id++) {

    if (!fBins[id].Size()) continue;

    Hit2D* candidate = fBins[id].GetHits2D().front();
    if (candidate) {
      double dist2 = pma::Dist2(candidate->GetPointCm(), fBaryCenter);
      if (dist2 > maxdist2) { maxdist2 = dist2; }
    }
  }

  fNormDist = std::sqrt(maxdist2);
}

void ems::DirOfGamma::FillBins ( )

private

Definition at line 202 of file DirOfGamma.cxx.

References bin, fBaryCenter, fBins, fNbins, fPoints2D, X, and Y.

Referenced by DirOfGamma().

{
  TVector2 vstart(0, 1);

  for (unsigned int i = 0; i < fPoints2D.size(); i++) {
    TVector2 pos(fPoints2D[i]->GetPointCm().X(), fPoints2D[i]->GetPointCm().Y());
    TVector2 vecp = pos - fBaryCenter;
    float sinsign = (vstart.X() * vecp.Y()) - (vstart.Y() * vecp.X());
    float cosine = (vstart * vecp) / (vstart.Mod() * vecp.Mod());
    float theta = 180.0F * (std::acos(cosine)) / TMath::Pi();

    unsigned int id = 0;
    double bin = double(360.0) / double(fNbins);

    if (sinsign >= 0)
      id = int(theta / bin);
    else if (sinsign < 0)
      id = int(theta / bin) + (fNbins / 2);
    if (id > (fNbins - 1)) id = (fNbins - 1);

    fBins[id].Add(fPoints2D[i]);
  }

  for (unsigned int id = 0; id < fBins.size(); id++)
    fBins[id].Sort();
}

bool ems::DirOfGamma::FindCandidates ( )

private

Definition at line 247 of file DirOfGamma.cxx.

References pma::Dist2(), fBaryCenter, fBins, fCandidates, fNbins, fNormDist, fPoints2D, and ems::Hit2D::GetPointCm().

Referenced by DirOfGamma().

{
  float rad = 0.5F * fNormDist;
  unsigned int nbins = fNbins * 4;
  for (unsigned int id = 0; id < fNbins; id++) {

    if (!fBins[id].Size()) continue;

    std::vector<Hit2D*> points;
    Hit2D* candidate2D = fBins[id].GetHits2D().front();

    for (unsigned int i = 0; i < fPoints2D.size(); i++) {
      double distnorm = std::sqrt(pma::Dist2(candidate2D->GetPointCm(), fBaryCenter)) / fNormDist;
      double dist2 = pma::Dist2(candidate2D->GetPointCm(), fPoints2D[i]->GetPointCm());

      if ((distnorm > 0.5) && (dist2 < rad * rad)) points.push_back(fPoints2D[i]);
    }

    if (fBins[id].Size() > 1) {
      EndPoint ep(*candidate2D, points, nbins);
      fCandidates.push_back(ep);
    }
  }
  if (fCandidates.size())
    return true;
  else
    return false;
}

void ems::DirOfGamma::FindInitialPart ( )

private

Definition at line 286 of file DirOfGamma.cxx.

References pma::Dist2(), fBaryCenter, fCandidateID, fCandidates, fIniHits, fStartHit, and fStartPoint.

Referenced by DirOfGamma().

{
  double max_asymmetry = 0.0;
  unsigned int saveid = 0;
  bool found = false;

  double maxdist2 = 0.0;
  double maxcharge = 0.0;
  unsigned int idmaxdist = 0;
  unsigned int idmaxcharge = 0;

  for (unsigned int i = 0; i < fCandidates.size(); i++) {
    double dist2 = pma::Dist2(fCandidates[i].GetPosition(), fBaryCenter);
    double charge = fCandidates[i].GetMaxCharge();
    if (dist2 > maxdist2) {
      maxdist2 = dist2;
      idmaxdist = i;
    }
    if (charge > maxcharge) {
      maxcharge = charge;
      idmaxcharge = i;
    }
  }

  maxdist2 = 0.0;
  unsigned int idmaxdistb = 0;
  for (size_t i = 0; i < fCandidates.size(); i++) {
    if ((i == idmaxdist) || (i == idmaxcharge)) continue;

    double dist2 = pma::Dist2(fCandidates[i].GetPosition(), fCandidates[idmaxdist].GetPosition());
    if (dist2 > maxdist2) {
      maxdist2 = dist2;
      idmaxdistb = i;
    }
  }

  if (fCandidates.size() > 2) {
    for (unsigned int i = 0; i < fCandidates.size(); i++) {
      double asymmetry = fCandidates[i].GetAsymmetry();

      if ((i == idmaxdist) || (i == idmaxcharge) || (i == idmaxdistb)) {
        if (asymmetry > max_asymmetry) {
          max_asymmetry = asymmetry;
          saveid = i;
          found = true;
        }
      }
    }
  }
  else {
    for (unsigned int i = 0; i < fCandidates.size(); i++) {
      double asymmetry = fCandidates[i].GetAsymmetry();

      if ((i == idmaxdist) || (i == idmaxdistb)) {
        if (asymmetry > max_asymmetry) {
          max_asymmetry = asymmetry;
          saveid = i;
          found = true;
        }
      }
    }
  }

  if (!found)
    mf::LogError("DirOfGamma") << fCandidates.size() << "DirOfGamma - Find Initial Part problem.";

  fStartHit = fCandidates[saveid].GetHit();
  fStartPoint = fCandidates[saveid].GetPosition();
  fIniHits = fCandidates[saveid].MaxChargeBin().GetIniHits();
  fCandidateID = saveid;
}

TVector2 const& ems::DirOfGamma::GetBaryCenterCm ( ) const

inline

Definition at line 131 of file DirOfGamma.h.

{ return fBaryCenter; }

std::vector<EndPoint> const& ems::DirOfGamma::GetCandidates ( ) const

inline

Definition at line 135 of file DirOfGamma.h.

{ return fCandidates; }

art::Ptr<recob::Hit> const& ems::DirOfGamma::GetFirstHit ( ) const

inline

Definition at line 149 of file DirOfGamma.h.

{ return fStartHit; }

TVector2 const& ems::DirOfGamma::GetFirstPoint ( ) const

inline

Definition at line 153 of file DirOfGamma.h.

{ return fStartPoint; }

std::vector<art::Ptr<recob::Hit> > const& ems::DirOfGamma::GetHits ( )

inline

Definition at line 151 of file DirOfGamma.h.

{ return fHits; }

std::vector<Hit2D*> const& ems::DirOfGamma::GetHits2D ( ) const

inline

Definition at line 133 of file DirOfGamma.h.

Referenced by ems::EMShower3D::CollectShower2D().

{ return fPoints2D; }

size_t ems::DirOfGamma::GetIdCandidate ( )

inline

Definition at line 147 of file DirOfGamma.h.

{ return fCandidateID; }

size_t ems::DirOfGamma::GetIdCl ( ) const

inline

Definition at line 157 of file DirOfGamma.h.

{ return fIdCl; }

std::vector<art::Ptr<recob::Hit> > const& ems::DirOfGamma::GetIniHits ( ) const

inline

Definition at line 155 of file DirOfGamma.h.

{ return fIniHits; }

void ems::DirOfGamma::SetIdCandidate ( size_t  id )

inline

Definition at line 137 of file DirOfGamma.h.

  {

    fIsCandidateIDset = true;
    fCandidateID = id;
    fStartHit = fCandidates[id].GetHit();
    fStartPoint = fCandidates[id].GetPosition();
    fIniHits = fCandidates[id].MaxChargeBin().GetIniHits();
  }

Member Data Documentation

TVector2 ems::DirOfGamma::fBaryCenter

private

Definition at line 183 of file DirOfGamma.h.

Referenced by ComputeBaryCenter(), ComputeMaxDist(), DirOfGamma(), FillBins(), FindCandidates(), and FindInitialPart().

std::vector<Bin2D> ems::DirOfGamma::fBins

private

Definition at line 165 of file DirOfGamma.h.

Referenced by ComputeMaxDist(), DirOfGamma(), FillBins(), and FindCandidates().

size_t ems::DirOfGamma::fCandidateID

private

Definition at line 162 of file DirOfGamma.h.

Referenced by FindInitialPart().

std::vector<EndPoint> ems::DirOfGamma::fCandidates

private

Definition at line 166 of file DirOfGamma.h.

Referenced by ComputeMaxCharge(), FindCandidates(), and FindInitialPart().

std::vector<art::Ptr<recob::Hit> > ems::DirOfGamma::fHits

private

Definition at line 171 of file DirOfGamma.h.

Referenced by DirOfGamma().

size_t ems::DirOfGamma::fIdCl

private

Definition at line 161 of file DirOfGamma.h.

std::vector<art::Ptr<recob::Hit> > ems::DirOfGamma::fIniHits

private

Definition at line 170 of file DirOfGamma.h.

Referenced by FindInitialPart().

bool ems::DirOfGamma::fIsCandidateIDset

private

Definition at line 188 of file DirOfGamma.h.

size_t ems::DirOfGamma::fNbins

private

Definition at line 160 of file DirOfGamma.h.

Referenced by DirOfGamma(), FillBins(), and FindCandidates().

float ems::DirOfGamma::fNormCharge

private

Definition at line 186 of file DirOfGamma.h.

Referenced by ComputeMaxCharge().

float ems::DirOfGamma::fNormDist

private

Definition at line 185 of file DirOfGamma.h.

Referenced by ComputeMaxDist(), and FindCandidates().

std::vector<Hit2D*> ems::DirOfGamma::fPoints2D

private

Definition at line 164 of file DirOfGamma.h.

Referenced by ComputeBaryCenter(), DirOfGamma(), FillBins(), and FindCandidates().

art::Ptr<recob::Hit> ems::DirOfGamma::fStartHit

private

Definition at line 168 of file DirOfGamma.h.

Referenced by FindInitialPart().

TVector2 ems::DirOfGamma::fStartPoint

private

Definition at line 169 of file DirOfGamma.h.

Referenced by FindInitialPart().

---

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

• larreco/v09_26_01/source/larreco/DirOfGamma/DirOfGamma.h
• larreco/v09_26_01/source/larreco/DirOfGamma/DirOfGamma.cxx