ems::ShowerInfo Class Reference

Public Member Functions
	ShowerInfo (int key, int gid, bool hasvtx, std::vector< double > vensum, recob::Track const &trk)
double	Pointsto (ems::ShowerInfo const &s1) const
double	Angleto (TVector3 const &pmapoint) const
bool	HasConPoint (void) const
int	GetKey (void) const
int	GetGid (void) const
int	GetPlaneId (void) const
double	GetAdcSum (void) const
std::vector< double >	GetEnTot (void) const
TVector3	GetFront (void) const
TVector3	GetEnd (void) const
TVector3	GetDir (void) const
std::vector< double >	GetDedx (void) const
double	GetP0Dist (void) const
void	SetP0Dist (const TVector3 p)

Private Attributes
int	fKey
int	fGId
int	fPlaneId
bool	fHasVtx
double	fAdcSum
std::vector< double >	fVentot
double	fP0Dist
TVector3	fFront
TVector3	fEnd
TVector3	fDir
std::vector< double >	fVdqdx

Detailed Description

Definition at line 53 of file MergeEMShower3D_module.cc.

Constructor & Destructor Documentation

ems::ShowerInfo::ShowerInfo	(	int	key,
		int	gid,
		bool	hasvtx,
		std::vector< double >	vensum,
		recob::Track const &	trk
	)

Definition at line 96 of file MergeEMShower3D_module.cc.

References recob::Track::DQdxAtPoint(), recob::Track::End(), fDir, fEnd, fFront, fHasVtx, fPlaneId, fVdqdx, geo::GeometryCore::Nplanes(), recob::Track::Vertex(), recob::Track::VertexDirection(), and geo::GeometryCore::View().

                                                                                                         :
fKey(key),
fGId(gid),
fVentot(vensum), 
fP0Dist(0)
{
        fDir = trk.VertexDirection();
        fFront = trk.Vertex();
        fEnd = trk.End();                       
        fHasVtx = hasvtx;

        art::ServiceHandle<geo::Geometry> geom;
        /*fVdqdx.push_back(trk.DQdxAtPoint(0, geo::kU)); 
        fVdqdx.push_back(trk.DQdxAtPoint(0, geo::kV));
        fVdqdx.push_back(trk.DQdxAtPoint(0, geo::kZ));*/ 

        
        bool isplane = false;
        unsigned int nplanes = geom->Nplanes();
        for (unsigned int p = 0; p < nplanes; ++p)
        {
                geo::View_t view = geom->View(p);
                fVdqdx.push_back(trk.DQdxAtPoint(0, view));
                if (trk.DQdxAtPoint(1, view) > 0) 
                {fPlaneId = int(p); isplane = true;}
        }
        if (!isplane) fPlaneId = -1; 

        /*if (trk.DQdxAtPoint(1, geo::kU) > 0) fPlaneId = geo::kU;
        else if (trk.DQdxAtPoint(1, geo::kV) > 0) fPlaneId = geo::kV;
        else if (trk.DQdxAtPoint(1, geo::kZ) > 0) fPlaneId = geo::kZ; 
        else fPlaneId = -1;*/
}

Member Function Documentation

double ems::ShowerInfo::Angleto

(

TVector3 const &

pmapoint

)

const

Definition at line 151 of file MergeEMShower3D_module.cc.

References fDir, and fFront.

{
        double cosine = (pmapoint - this->fFront) * (1.0 / (pmapoint - this->fFront).Mag()) * this->fDir;
        double th = 180.0F * (std::acos(cosine)) / TMath::Pi();
        
        return th;      
}

double ems::ShowerInfo::GetAdcSum ( void  ) const

inline

Definition at line 64 of file MergeEMShower3D_module.cc.

References fAdcSum.

{return fAdcSum;} 

std::vector<double> ems::ShowerInfo::GetDedx ( void  ) const

inline

Definition at line 70 of file MergeEMShower3D_module.cc.

References fVdqdx.

{ return fVdqdx; } 

TVector3 ems::ShowerInfo::GetDir ( void  ) const

inline

Definition at line 69 of file MergeEMShower3D_module.cc.

References fDir.

{return fDir;}

TVector3 ems::ShowerInfo::GetEnd ( void  ) const

inline

Definition at line 68 of file MergeEMShower3D_module.cc.

References fEnd.

{return fEnd;}

std::vector<double> ems::ShowerInfo::GetEnTot ( void  ) const

inline

Definition at line 65 of file MergeEMShower3D_module.cc.

References fVentot.

{return fVentot; } 

TVector3 ems::ShowerInfo::GetFront ( void  ) const

inline

Definition at line 67 of file MergeEMShower3D_module.cc.

References fFront.

{return fFront;}

int ems::ShowerInfo::GetGid ( void  ) const

inline

Definition at line 62 of file MergeEMShower3D_module.cc.

References fGId.

Referenced by ems::ShowersCollection::Merge(), and ems::ShowersCollection::ShowersCollection().

{return fGId;}

int ems::ShowerInfo::GetKey ( void  ) const

inline

Definition at line 61 of file MergeEMShower3D_module.cc.

References fKey.

{return fKey;}

double ems::ShowerInfo::GetP0Dist ( void  ) const

inline

Definition at line 72 of file MergeEMShower3D_module.cc.

References fP0Dist.

Referenced by ems::bDistLess::operator()().

{ return fP0Dist; }

int ems::ShowerInfo::GetPlaneId ( void  ) const

inline

Definition at line 63 of file MergeEMShower3D_module.cc.

References fPlaneId.

{return fPlaneId;}

bool ems::ShowerInfo::HasConPoint ( void  ) const

inline

Definition at line 60 of file MergeEMShower3D_module.cc.

References fHasVtx.

{return fHasVtx;}

double ems::ShowerInfo::Pointsto

(

ems::ShowerInfo const &

s1

)

const

Definition at line 130 of file MergeEMShower3D_module.cc.

References fDir, fFront, pma::GetProjectionToSegment(), and pma::SolveLeastSquares3D().

Referenced by ems::ShowersCollection::MinAngle().

{
        double cosine0 = (this->fFront - s1.fFront) * (1.0 / (this->fFront - s1.fFront).Mag()) * s1.fDir;
        double th0 = 180.0F * (std::acos(cosine0)) / TMath::Pi();

        std::vector< std::pair<TVector3, TVector3> > lines;
        lines.push_back(std::pair<TVector3, TVector3>(this->fFront, this->fFront + this->fDir));
        lines.push_back(std::pair<TVector3, TVector3>(s1.fFront, s1.fFront + s1.fDir));

        TVector3 isect, pm;
        pma::SolveLeastSquares3D(lines, isect);
        pm = pma::GetProjectionToSegment(isect, this->fFront, this->fFront +  this->fDir); 
        double cosine1 = (pm - s1.fFront) * (1.0/(pm - s1.fFront).Mag()) * s1.fDir;
        double th1 = 180.0F * (std::acos(cosine1)) / TMath::Pi();

        double thmin = th0;
        if (th1 < th0) thmin = th1;     

        return thmin;
}

void ems::ShowerInfo::SetP0Dist ( const TVector3  p )

inline

Definition at line 73 of file MergeEMShower3D_module.cc.

References pma::Dist2(), fEnd, fFront, fHasVtx, and fP0Dist.

        {
                if (fHasVtx) fP0Dist = std::sqrt( pma::Dist2(p, fFront) );
                else fP0Dist = std::sqrt( pma::Dist2(p, 0.5 * (fFront + fEnd)) );
        }

Member Data Documentation

double ems::ShowerInfo::fAdcSum

private

Definition at line 84 of file MergeEMShower3D_module.cc.

Referenced by GetAdcSum().

TVector3 ems::ShowerInfo::fDir

private

Definition at line 90 of file MergeEMShower3D_module.cc.

Referenced by Angleto(), GetDir(), Pointsto(), and ShowerInfo().

TVector3 ems::ShowerInfo::fEnd

private

Definition at line 89 of file MergeEMShower3D_module.cc.

Referenced by GetEnd(), SetP0Dist(), and ShowerInfo().

TVector3 ems::ShowerInfo::fFront

private

Definition at line 88 of file MergeEMShower3D_module.cc.

Referenced by Angleto(), GetFront(), Pointsto(), SetP0Dist(), and ShowerInfo().

int ems::ShowerInfo::fGId

private

Definition at line 81 of file MergeEMShower3D_module.cc.

Referenced by GetGid().

bool ems::ShowerInfo::fHasVtx

private

Definition at line 83 of file MergeEMShower3D_module.cc.

Referenced by HasConPoint(), SetP0Dist(), and ShowerInfo().

int ems::ShowerInfo::fKey

private

Definition at line 80 of file MergeEMShower3D_module.cc.

Referenced by GetKey().

double ems::ShowerInfo::fP0Dist

private

Definition at line 86 of file MergeEMShower3D_module.cc.

Referenced by GetP0Dist(), and SetP0Dist().

int ems::ShowerInfo::fPlaneId

private

Definition at line 82 of file MergeEMShower3D_module.cc.

Referenced by GetPlaneId(), and ShowerInfo().

std::vector<double> ems::ShowerInfo::fVdqdx

private

Definition at line 92 of file MergeEMShower3D_module.cc.

Referenced by GetDedx(), and ShowerInfo().

std::vector<double> ems::ShowerInfo::fVentot

private

Definition at line 85 of file MergeEMShower3D_module.cc.

Referenced by GetEnTot().

---

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

• larreco/v06_64_02/source/larreco/DirOfGamma/MergeEMShower3D_module.cc