ems::EndPoint Class Reference

#include <DirOfGamma.h>

Public Member Functions
	EndPoint (const Hit2D &center, const std::vector< Hit2D * > &hits, unsigned int nbins)
TVector2 const &	GetPosition () const
double	GetAsymmetry () const
double	GetMaxCharge () const
Bin2D const &	MaxChargeBin () const
std::vector< Bin2D > const &	GetBins2D () const
art::Ptr< recob::Hit > const &	GetHit () const
const std::vector< art::Ptr < recob::Hit > >	GetIniHits () const
size_t const &	GetPlane () const
size_t const &	GetTPC () const
size_t const &	GetCryo () const

Private Member Functions
void	FillBins ()
void	ComputeMaxCharge ()
void	ComputeMeanCharge ()

Private Attributes
Hit2D	fCenter2D
std::vector< Hit2D * >	fPoints2D
size_t	fNbins
double	fMaxCharge
double	fMeanCharge
std::vector< Bin2D >	fBins
size_t	fMaxChargeIdBin
size_t	fPlane
size_t	fTpc
size_t	fCryo

Detailed Description

Definition at line 104 of file DirOfGamma.h.

Constructor & Destructor Documentation

ems::EndPoint::EndPoint	(	const Hit2D &	center,
		const std::vector< Hit2D * > &	hits,
		unsigned int	nbins
	)

Definition at line 77 of file DirOfGamma.cxx.

  : fCenter2D(center), fPoints2D(hits), fNbins(nbins)
{

  for (unsigned int i = 0; i < fNbins; i++) {
    fBins.push_back(Bin2D(center.GetPointCm()));
  }

  FillBins();
  ComputeMaxCharge();
  ComputeMeanCharge();

  fPlane = center.GetHitPtr()->WireID().Plane;
  fTpc = center.GetHitPtr()->WireID().TPC;
  fCryo = center.GetHitPtr()->WireID().Cryostat;
}

Member Function Documentation

void ems::EndPoint::ComputeMaxCharge ( )

private

Definition at line 134 of file DirOfGamma.cxx.

{
  fMaxCharge = 0.0;
  unsigned int saveid = 0;
  for (unsigned int i = 0; i < fNbins; i++)
    if (fBins[i].Size() && (fMaxCharge < fBins[i].GetTotCharge())) {
      fMaxCharge = fBins[i].GetTotCharge();
      saveid = i;
    }

  fMaxChargeIdBin = saveid;
}

void ems::EndPoint::ComputeMeanCharge ( )

private

Definition at line 148 of file DirOfGamma.cxx.

{
  fMeanCharge = 0.0;
  if (fNbins == 0) return;

  unsigned int iprev, inext;

  if (fMaxChargeIdBin > 0)
    iprev = fMaxChargeIdBin - 1;
  else
    iprev = fNbins - 1;

  inext = (fMaxChargeIdBin + 1) % fNbins;

  double sumcharge = 0.0;
  for (unsigned int i = 0; i < fNbins; i++)
    if ((i != fMaxChargeIdBin) && (i != iprev) && (i != inext))
      sumcharge += fBins[i].GetTotCharge();

  fMeanCharge = sumcharge / double(fNbins);
}

void ems::EndPoint::FillBins ( )

private

Definition at line 95 of file DirOfGamma.cxx.

{
  TVector2 vstart(0, 1);

  unsigned int saveid = 0;
  bool exist = false;
  for (unsigned int i = 0; i < fPoints2D.size(); i++) {
    if (fPoints2D[i]->GetHitPtr().key() != fCenter2D.GetHitPtr().key()) {
      TVector2 pos(fPoints2D[i]->GetPointCm());
      TVector2 centre(fCenter2D.GetPointCm());
      TVector2 vecp = pos - centre;
      float sinsign = (vstart.X() * vecp.Y()) - (vstart.Y() * vecp.X());
      float cosine = (vstart * vecp) / 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]);
      fBins[(id + 1) % fNbins].Add(fPoints2D[i]);
    }
    else {
      saveid = i;
      exist = true;
    }
  }

  if (exist)
    for (unsigned int id = 0; id < fNbins; id++)
      fBins[id].Add(fPoints2D[saveid]);
}

double ems::EndPoint::GetAsymmetry

(

)

const

Definition at line 171 of file DirOfGamma.cxx.

{
  if ((fMaxCharge + fMeanCharge) == 0) return 0.0;
  return ((fMaxCharge - fMeanCharge) / (fMaxCharge + fMeanCharge));
}

std::vector<Bin2D> const& ems::EndPoint::GetBins2D ( ) const

inline

Definition at line 129 of file DirOfGamma.h.

  {
    return fBins;
  }

size_t const& ems::EndPoint::GetCryo ( ) const

inline

Definition at line 157 of file DirOfGamma.h.

  {
    return fCryo;
  }

art::Ptr<recob::Hit> const& ems::EndPoint::GetHit ( ) const

inline

Definition at line 135 of file DirOfGamma.h.

  {
    return fCenter2D.GetHitPtr();
  }

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

inline

Definition at line 141 of file DirOfGamma.h.

  {
    return MaxChargeBin().GetIniHits();
  }

double ems::EndPoint::GetMaxCharge ( ) const

inline

Definition at line 117 of file DirOfGamma.h.

  {
    return fMaxCharge;
  }

size_t const& ems::EndPoint::GetPlane ( ) const

inline

Definition at line 147 of file DirOfGamma.h.

  {
    return fPlane;
  }

TVector2 const& ems::EndPoint::GetPosition ( ) const

inline

Definition at line 109 of file DirOfGamma.h.

  {
    return fCenter2D.GetPointCm();
  }

size_t const& ems::EndPoint::GetTPC ( ) const

inline

Definition at line 152 of file DirOfGamma.h.

  {
    return fTpc;
  }

Bin2D const& ems::EndPoint::MaxChargeBin ( ) const

inline

Definition at line 123 of file DirOfGamma.h.

  {
    return fBins[fMaxChargeIdBin];
  }

Member Data Documentation

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

private

Definition at line 170 of file DirOfGamma.h.

Hit2D ems::EndPoint::fCenter2D

private

Definition at line 163 of file DirOfGamma.h.

size_t ems::EndPoint::fCryo

private

Definition at line 180 of file DirOfGamma.h.

double ems::EndPoint::fMaxCharge

private

Definition at line 167 of file DirOfGamma.h.

size_t ems::EndPoint::fMaxChargeIdBin

private

Definition at line 172 of file DirOfGamma.h.

double ems::EndPoint::fMeanCharge

private

Definition at line 168 of file DirOfGamma.h.

size_t ems::EndPoint::fNbins

private

Definition at line 165 of file DirOfGamma.h.

size_t ems::EndPoint::fPlane

private

Definition at line 178 of file DirOfGamma.h.

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

private

Definition at line 164 of file DirOfGamma.h.

size_t ems::EndPoint::fTpc

private

Definition at line 179 of file DirOfGamma.h.

---

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

• DirOfGamma.h
• DirOfGamma.cxx