genie::geometry::FidSphere Class Reference

#include <FidShape.h>

Inheritance diagram for genie::geometry::FidSphere:

Collaboration diagram for genie::geometry::FidSphere:

Public Member Functions
	FidSphere (const TVector3 &center, Double_t radius)
RayIntercept	Intercept (const TVector3 &start, const TVector3 &dir) const
void	ConvertMaster2Top (const ROOTGeomAnalyzer *rgeom)
void	Print (std::ostream &stream) const
Public Member Functions inherited from genie::geometry::FidShape
	FidShape ()
virtual	~FidShape ()

Protected Attributes
TVector3	fCenter
Double_t	fSRadius
	center of the sphere More...

Detailed Description

Definition at line 107 of file FidShape.h.

Constructor & Destructor Documentation

genie::geometry::FidSphere::FidSphere ( const TVector3 &  center, Double_t  radius  )

inline

Definition at line 109 of file FidShape.h.

: fCenter(center), fSRadius(radius) { ; }

Member Function Documentation

void FidSphere::ConvertMaster2Top ( const ROOTGeomAnalyzer *  rgeom )

virtual

derived classes must implement the ConvertMaster2Top() method which transforms the shape specification from master coordinates to "top vol"

Implements genie::geometry::FidShape.

Definition at line 98 of file FidShape.cxx.

References genie::geometry::ROOTGeomAnalyzer::Master2Top().

{
  rgeom->Master2Top(fCenter);
}

RayIntercept FidSphere::Intercept ( const TVector3 &  start, const TVector3 &  dir  ) const

virtual

derived classes must implement the Intercept() method which calculates the entry/exit point of a ray w/ the shape

Implements genie::geometry::FidShape.

Definition at line 69 of file FidShape.cxx.

References dir, genie::geometry::RayIntercept::fDistIn, genie::geometry::RayIntercept::fDistOut, genie::geometry::RayIntercept::fIsHit, genie::geometry::RayIntercept::fSurfIn, and genie::geometry::RayIntercept::fSurfOut.

{
  // A new neutrino ray has been set, calculate the entrance/exit distances.
  // This sets fDistIn/fDistOut for an sphere
  RayIntercept intercept;

  TVector3 oc = fCenter - start;
  Double_t loc2 = oc.Mag2();
  Double_t r2 = fSRadius*fSRadius;
  //LOG("GeomVolSel", pNOTICE) << " loc2 = " << loc2 << " r2 " << r2;
  // if ( loc2 > r2 ) ray originates outside the sphere
  const TVector3& d = dir;
  Double_t d2 = d.Mag2();
  Double_t tca = oc.Dot(d)/d2;
  //if ( tca < 0.0 ) sphere _center_ behind the ray orgin
  Double_t lhc2 = ( r2 -loc2 )/d2 + tca*tca;
  if ( lhc2 < 0.0 ) return intercept; // ray misses the sphere
  intercept.fIsHit = true;
  Double_t lhc = TMath::Sqrt(lhc2);

  intercept.fDistIn  = tca - lhc;
  intercept.fSurfIn  = 1;
  intercept.fDistOut = tca + lhc;
  intercept.fSurfOut = 1;

  return intercept;
}

void FidSphere::Print ( std::ostream &  stream ) const

virtual

Implements genie::geometry::FidShape.

Definition at line 104 of file FidShape.cxx.

{
  stream << "FidSphere @ ["
         << fCenter.X() << ","
         << fCenter.Y() << ","
         << fCenter.Z() << "]  r = " << fSRadius;
}

Member Data Documentation

TVector3 genie::geometry::FidSphere::fCenter

protected

Definition at line 114 of file FidShape.h.

Double_t genie::geometry::FidSphere::fSRadius

protected

center of the sphere

Definition at line 115 of file FidShape.h.

---

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

• FidShape.h
• FidShape.cxx