ems::DirOfGamma Class Reference

#include "DirOfGamma.h"

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

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

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

Constructor & Destructor Documentation

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

Definition at line 180 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(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 122 of file DirOfGamma.h.

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

Member Function Documentation

void ems::DirOfGamma::ComputeBaryCenter ( void  )

private

Definition at line 208 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 ( void  )

private

void ems::DirOfGamma::ComputeMaxCharge ( void  )

private

Definition at line 301 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 ( void  )

private

Definition at line 248 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 ( void  )

private

Definition at line 223 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 ( void  )

private

Definition at line 271 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 ( void  )

private

Definition at line 313 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 ( void  ) const

inline

Definition at line 124 of file DirOfGamma.h.

{ return fBaryCenter; }

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

inline

Definition at line 128 of file DirOfGamma.h.

{ return fCandidates; }

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

inline

Definition at line 140 of file DirOfGamma.h.

{ return fStartHit; }

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

inline

Definition at line 144 of file DirOfGamma.h.

{ return fStartPoint; }

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

inline

Definition at line 142 of file DirOfGamma.h.

{ return fHits; }

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

inline

Definition at line 126 of file DirOfGamma.h.

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

{ return fPoints2D; }  

const size_t ems::DirOfGamma::GetIdCandidate ( void  )

inline

Definition at line 138 of file DirOfGamma.h.

{ return fCandidateID; }

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

inline

Definition at line 148 of file DirOfGamma.h.

{return fIdCl;}

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

inline

Definition at line 146 of file DirOfGamma.h.

{ return fIniHits;}

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

inline

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

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

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

private

Definition at line 156 of file DirOfGamma.h.

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

size_t ems::DirOfGamma::fCandidateID

private

Definition at line 153 of file DirOfGamma.h.

Referenced by FindInitialPart().

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

private

Definition at line 157 of file DirOfGamma.h.

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

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

private

Definition at line 162 of file DirOfGamma.h.

Referenced by DirOfGamma().

size_t ems::DirOfGamma::fIdCl

private

Definition at line 152 of file DirOfGamma.h.

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

private

Definition at line 161 of file DirOfGamma.h.

Referenced by FindInitialPart().

bool ems::DirOfGamma::fIsCandidateIDset

private

Definition at line 179 of file DirOfGamma.h.

size_t ems::DirOfGamma::fNbins

private

Definition at line 151 of file DirOfGamma.h.

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

float ems::DirOfGamma::fNormCharge

private

Definition at line 177 of file DirOfGamma.h.

Referenced by ComputeMaxCharge().

float ems::DirOfGamma::fNormDist

private

Definition at line 176 of file DirOfGamma.h.

Referenced by ComputeMaxDist(), and FindCandidates().

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

private

Definition at line 155 of file DirOfGamma.h.

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

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

private

Definition at line 159 of file DirOfGamma.h.

Referenced by FindInitialPart().

TVector2 ems::DirOfGamma::fStartPoint

private

Definition at line 160 of file DirOfGamma.h.

Referenced by FindInitialPart().

---

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

• /cvmfs/larsoft.opensciencegrid.org/products/larreco/v07_10_02/source/larreco/DirOfGamma/DirOfGamma.h
• /cvmfs/larsoft.opensciencegrid.org/products/larreco/v07_10_02/source/larreco/DirOfGamma/DirOfGamma.cxx