genie::geometry::PathSegment Class Reference

#include <PathSegmentList.h>

Public Member Functions
	PathSegment ()
	~PathSegment ()
void	SetEnter (const TVector3 &p3enter, double raydist)
	point of entry to geometry element More...
void	SetEnter (const Double_t *p3enter, double raydist)
void	SetExit (const TVector3 &p3exit)
	point of exit from geometry element More...
void	SetExit (const Double_t *p3exit)
void	SetGeo (const TGeoVolume *gvol, const TGeoMedium *gmed, const TGeoMaterial *gmat)
	info about the geometry element More...
void	SetPath (const char *path)
void	SetStep (Double_t step, bool setlimits=true)
	step taken in the geometry element More...
bool	IsTrimmedEmpty () const
Double_t	GetSummedStepRange () const
	get the sum of all the step range (in case step has been trimmed or split) More...
TVector3	GetPosition (Double_t frac) const
	calculate position within allowed ranges passed on fraction of total More...
void	DoCrossCheck (const TVector3 &startpos, double &ddist, double &dstep) const
	perform cross check on segment, return differences More...
void	Print (ostream &stream) const

Public Attributes
Double_t	fRayDist
	distance from start of ray More...
Double_t	fStepLength
	total step size in volume More...
const TGeoVolume *	fVolume
	ref only ptr to TGeoVolume More...
const TGeoMedium *	fMedium
	ref only ptr to TGeoMedium More...
const TGeoMaterial *	fMaterial
	ref only ptr to TGeoMaterial More...
TVector3	fEnter
	top vol coordinates and units More...
TVector3	fExit
	top vol coordinates and units More...
std::string	fPathString
	full path names More...
StepRangeSet	fStepRangeSet
	collection of {steplo,stephi} pairs More...

Friends
ostream &	operator<< (ostream &stream, const PathSegment &list)
bool	operator< (const PathSegment &lhs, const PathSegment &rhs)

Detailed Description

Definition at line 56 of file PathSegmentList.h.

Constructor & Destructor Documentation

PathSegment::PathSegment

(

void

)

Definition at line 76 of file PathSegmentList.cxx.

                             :
  fRayDist(0), fStepLength(0),
  fVolume(0), fMedium(0), fMaterial(0),
  fEnter(), fExit()
{
}

genie::geometry::PathSegment::~PathSegment ( )

inline

Definition at line 60 of file PathSegmentList.h.

{ ; }

Member Function Documentation

void PathSegment::DoCrossCheck	(	const TVector3 &	startpos,
		double &	ddist,
		double &	dstep
	)		const

perform cross check on segment, return differences

Definition at line 84 of file PathSegmentList.cxx.

{
  double dist_recalc = (fEnter-startpos).Mag();
  ddist = dist_recalc - fRayDist;

  double step_recalc = (fExit-fEnter).Mag();
  dstep = step_recalc - fStepLength;
}

TVector3 PathSegment::GetPosition

(

Double_t

frac

)

const

calculate position within allowed ranges passed on fraction of total

calculate position within allowed ranges passed as fraction of trimmed segment seg.fEnter + fractotal * ( seg.fExit - seg.fEnter );

Definition at line 148 of file PathSegmentList.cxx.

{
  /// calculate position within allowed ranges passed as
  /// fraction of trimmed segment
  ///      seg.fEnter + fractotal * ( seg.fExit - seg.fEnter );
  Double_t sumrange = GetSummedStepRange();
  if ( sumrange < 0.0 ) {
    LOG("PathS", pFATAL) << "GetPosition failed fractrim=" << fractrim
                         << " because sumrange = " << sumrange;
    return TVector3(0,0,0);
  }
  Double_t target   = fractrim * sumrange;
  Double_t sum      = 0;
  for ( size_t i = 0; i < fStepRangeSet.size(); ++i ) {
    const StepRange& sr = fStepRangeSet[i];
    Double_t ds = ( sr.second - sr.first );
    sum += ds;
#ifdef RWH_DEBUG
    LOG("PathS", pINFO) << "GetPosition fractrim=" << fractrim
                        << " target " << target << " [" << i << "] "
                        << " ds " << ds << " sum " << sum;
#endif
    if ( sum >= target ) {
      Double_t overstep  = sum - target;
      Double_t fractotal = (sr.second - overstep)/fStepLength;
#ifdef RWH_DEBUG
    LOG("PathS", pINFO) << "GetPosition fractrim=" << fractrim
                        << " overstep " << overstep
                        << " fractotal " << fractotal;
#endif
      return  fEnter + fractotal * ( fExit - fEnter );
    }
  }
  LOG("PathS", pFATAL) << "GetPosition failed fractrim=" << fractrim;
  return TVector3(0,0,0);
}

Double_t PathSegment::GetSummedStepRange

(

)

const

get the sum of all the step range (in case step has been trimmed or split)

Definition at line 138 of file PathSegmentList.cxx.

{
  Double_t sum = 0;
  for ( size_t i = 0; i < fStepRangeSet.size(); ++i ) {
    const StepRange& sr = fStepRangeSet[i];
    sum += ( sr.second - sr.first );
  }
  return sum;
}

bool genie::geometry::PathSegment::IsTrimmedEmpty ( ) const

inline

Definition at line 84 of file PathSegmentList.h.

{ return fStepRangeSet.empty(); }

void PathSegment::Print

(

ostream &

stream

)

const

Definition at line 95 of file PathSegmentList.cxx.

{
  const char* vname = (fVolume)   ? fVolume->GetName()   : "no volume";
  const char* mname = (fMaterial) ? fMaterial->GetName() : "no material";
  stream << genie::pathsegutils::Vec3AsString(&fEnter) << " "
    //<< genie::pathsegutils::Vec3AsString(&fExit)
         << " " // "raydist "
         << std::setw(12) << fRayDist
         << " " // "step "
         << std::setw(12) << fStepLength << " "
         << std::left
         << std::setw(16) << vname << " '"
         << std::setw(18) << mname << "' ";
  size_t n = fStepRangeSet.size();
  const int rngw = 24;
  if ( n == 0 ) {
    stream << std::setw(rngw) << "[ ]";
  } else {
    std::ostringstream rngset;
    for ( size_t i = 0 ; i < n; ++i ) {
      const StepRange& sr = fStepRangeSet[i];
      rngset << "[" << sr.first << ":" << sr.second << "]";
    }
    stream << std::setw(rngw) << rngset.str();
  }
  stream << std::right;
#ifdef PATHSEG_KEEP_PATH
  stream << " " << fPathString;
#endif

}

void genie::geometry::PathSegment::SetEnter ( const TVector3 &  p3enter, double  raydist  )

inline

point of entry to geometry element

Definition at line 63 of file PathSegmentList.h.

     { fEnter = p3enter; fRayDist = raydist; }

void genie::geometry::PathSegment::SetEnter ( const Double_t *  p3enter, double  raydist  )

inline

Definition at line 65 of file PathSegmentList.h.

     { fEnter.SetXYZ(p3enter[0],p3enter[1],p3enter[2]); fRayDist = raydist; }

void genie::geometry::PathSegment::SetExit ( const TVector3 &  p3exit )

inline

point of exit from geometry element

Definition at line 69 of file PathSegmentList.h.

{ fExit = p3exit;  }

void genie::geometry::PathSegment::SetExit ( const Double_t *  p3exit )

inline

Definition at line 70 of file PathSegmentList.h.

     { fExit.SetXYZ(p3exit[0],p3exit[1],p3exit[2]); }

void genie::geometry::PathSegment::SetGeo ( const TGeoVolume *  gvol, const TGeoMedium *  gmed, const TGeoMaterial *  gmat  )

inline

info about the geometry element

Definition at line 74 of file PathSegmentList.h.

  { fVolume = gvol; fMedium = gmed; fMaterial = gmat; }

void genie::geometry::PathSegment::SetPath ( const char *  path )

inline

Definition at line 78 of file PathSegmentList.h.

{ fPathString = path; }

void PathSegment::SetStep	(	Double_t	step,
		bool	setlimits = true
	)

step taken in the geometry element

Definition at line 128 of file PathSegmentList.cxx.

{
  fStepLength = step;
  if (setlimits) {
    fStepRangeSet.clear();
    fStepRangeSet.push_back(StepRange(0,step));
  }
}

Friends And Related Function Documentation

bool operator< ( const PathSegment &  lhs, const PathSegment &  rhs  )

friend

Definition at line 116 of file PathSegmentList.h.

  { return ( lhs.fRayDist < rhs.fRayDist ); }

ostream& operator<< ( ostream &  stream, const PathSegment &  list  )

friend

Definition at line 45 of file PathSegmentList.cxx.

 {
   ps.Print(stream);
   return stream;
 }

Member Data Documentation

TVector3 genie::geometry::PathSegment::fEnter

top vol coordinates and units

Definition at line 108 of file PathSegmentList.h.

TVector3 genie::geometry::PathSegment::fExit

top vol coordinates and units

Definition at line 109 of file PathSegmentList.h.

const TGeoMaterial* genie::geometry::PathSegment::fMaterial

ref only ptr to TGeoMaterial

Definition at line 107 of file PathSegmentList.h.

const TGeoMedium* genie::geometry::PathSegment::fMedium

ref only ptr to TGeoMedium

Definition at line 106 of file PathSegmentList.h.

std::string genie::geometry::PathSegment::fPathString

full path names

Definition at line 111 of file PathSegmentList.h.

Double_t genie::geometry::PathSegment::fRayDist

distance from start of ray

Definition at line 103 of file PathSegmentList.h.

Double_t genie::geometry::PathSegment::fStepLength

total step size in volume

Definition at line 104 of file PathSegmentList.h.

StepRangeSet genie::geometry::PathSegment::fStepRangeSet

collection of {steplo,stephi} pairs

Definition at line 113 of file PathSegmentList.h.

const TGeoVolume* genie::geometry::PathSegment::fVolume

ref only ptr to TGeoVolume

Definition at line 105 of file PathSegmentList.h.

---

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

• Generator/src/Tools/Geometry/PathSegmentList.h
• Generator/src/Tools/Geometry/PathSegmentList.cxx