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)
const 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 const	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 181 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 177 of file DirOfGamma.cxx.

  : 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 187 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 204 of file DirOfGamma.cxx.

{
  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 298 of file DirOfGamma.cxx.

{
  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 249 of file DirOfGamma.cxx.

{
  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 221 of file DirOfGamma.cxx.

{
  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 268 of file DirOfGamma.cxx.

{
  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 309 of file DirOfGamma.cxx.

{
  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 194 of file DirOfGamma.h.

  {
    return fBaryCenter;
  }

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

inline

Definition at line 206 of file DirOfGamma.h.

  {
    return fCandidates;
  }

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

inline

Definition at line 229 of file DirOfGamma.h.

  {
    return fStartHit;
  }

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

inline

Definition at line 241 of file DirOfGamma.h.

  {
    return fStartPoint;
  }

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

inline

Definition at line 235 of file DirOfGamma.h.

  {
    return fHits;
  }

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

inline

Definition at line 200 of file DirOfGamma.h.

  {
    return fPoints2D;
  }

const size_t ems::DirOfGamma::GetIdCandidate ( )

inline

Definition at line 223 of file DirOfGamma.h.

  {
    return fCandidateID;
  }

size_t const ems::DirOfGamma::GetIdCl ( ) const

inline

Definition at line 253 of file DirOfGamma.h.

  {
    return fIdCl;
  }

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

inline

Definition at line 247 of file DirOfGamma.h.

  {
    return fIniHits;
  }

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

inline

Definition at line 212 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 282 of file DirOfGamma.h.

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

private

Definition at line 264 of file DirOfGamma.h.

size_t ems::DirOfGamma::fCandidateID

private

Definition at line 261 of file DirOfGamma.h.

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

private

Definition at line 265 of file DirOfGamma.h.

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

private

Definition at line 270 of file DirOfGamma.h.

size_t ems::DirOfGamma::fIdCl

private

Definition at line 260 of file DirOfGamma.h.

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

private

Definition at line 269 of file DirOfGamma.h.

bool ems::DirOfGamma::fIsCandidateIDset

private

Definition at line 287 of file DirOfGamma.h.

size_t ems::DirOfGamma::fNbins

private

Definition at line 259 of file DirOfGamma.h.

float ems::DirOfGamma::fNormCharge

private

Definition at line 285 of file DirOfGamma.h.

float ems::DirOfGamma::fNormDist

private

Definition at line 284 of file DirOfGamma.h.

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

private

Definition at line 263 of file DirOfGamma.h.

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

private

Definition at line 267 of file DirOfGamma.h.

TVector2 ems::DirOfGamma::fStartPoint

private

Definition at line 268 of file DirOfGamma.h.

---

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

• larreco/larreco/DirOfGamma/DirOfGamma.h
• larreco/larreco/DirOfGamma/DirOfGamma.cxx