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

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

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::PathSegment

(

void

)

Definition at line 76 of file PathSegmentList.cxx.

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

References fEnter, fExit, fMaterial, fMedium, fRayDist, fStepLength, and fVolume.

Referenced by operator<.

~PathSegment()

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

inline

Definition at line 60 of file PathSegmentList.h.

{ ; }

Member Function Documentation

DoCrossCheck()

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;
}

References fEnter, fExit, fRayDist, and fStepLength.

Referenced by genie::geometry::PathSegmentList::CrossCheck(), and genie::geometry::PathSegmentList::Print().

GetPosition()

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);
}

References fEnter, fExit, fStepLength, fStepRangeSet, GetSummedStepRange(), LOG, pFATAL, and pINFO.

Referenced by genie::geometry::ROOTGeomAnalyzer::GenerateVertex().

GetSummedStepRange()

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;
}

References fStepRangeSet.

Referenced by genie::geometry::ROOTGeomAnalyzer::ComputeMatLengths(), genie::geometry::PathSegmentList::FillMatStepSum(), genie::geometry::GeomVolSelectorI::GenerateTrimmedList(), genie::geometry::ROOTGeomAnalyzer::GenerateVertex(), and GetPosition().

IsTrimmedEmpty()

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

inline

Definition at line 84 of file PathSegmentList.h.

{ return fStepRangeSet.empty(); }

References fStepRangeSet.

Print()

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
 
}

References fEnter, fMaterial, fPathString, fRayDist, fStepLength, fStepRangeSet, fVolume, and genie::pathsegutils::Vec3AsString().

Referenced by operator<<.

SetEnter() [1/2]

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; }

References fEnter, and fRayDist.

SetEnter() [2/2]

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; }

References fEnter, and fRayDist.

Referenced by genie::geometry::ROOTGeomAnalyzer::SwimOnce().

SetExit() [1/2]

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]); }

References fExit.

SetExit() [2/2]

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;  }

References fExit.

Referenced by genie::geometry::ROOTGeomAnalyzer::SwimOnce().

SetGeo()

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; }

References fMaterial, fMedium, and fVolume.

Referenced by genie::geometry::ROOTGeomAnalyzer::SwimOnce().

SetPath()

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

inline

Definition at line 78 of file PathSegmentList.h.

{ fPathString = path; }

References fPathString.

Referenced by genie::geometry::ROOTGeomAnalyzer::SwimOnce().

SetStep()

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));
  }
}

References fStepLength, and fStepRangeSet.

Referenced by genie::geometry::ROOTGeomAnalyzer::SwimOnce().

Friends And Related Symbol Documentation

operator<

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

friend

Definition at line 116 of file PathSegmentList.h.

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

References fRayDist, and PathSegment().

operator<<

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

friend

Definition at line 45 of file PathSegmentList.cxx.

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

References Print().

Member Data Documentation

fEnter

TVector3 genie::geometry::PathSegment::fEnter

top vol coordinates and units

Definition at line 108 of file PathSegmentList.h.

Referenced by DoCrossCheck(), GetPosition(), PathSegment(), Print(), SetEnter(), and SetEnter().

fExit

TVector3 genie::geometry::PathSegment::fExit

top vol coordinates and units

Definition at line 109 of file PathSegmentList.h.

Referenced by DoCrossCheck(), GetPosition(), PathSegment(), SetExit(), and SetExit().

fMaterial

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

ref only ptr to TGeoMaterial

Definition at line 107 of file PathSegmentList.h.

Referenced by genie::geometry::ROOTGeomAnalyzer::ComputeMatLengths(), genie::geometry::PathSegmentList::FillMatStepSum(), genie::geometry::ROOTGeomAnalyzer::GenerateVertex(), PathSegment(), Print(), SetGeo(), and genie::geometry::GeomVolSelectorBasic::TrimSegment().

fMedium

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

ref only ptr to TGeoMedium

Definition at line 106 of file PathSegmentList.h.

Referenced by PathSegment(), SetGeo(), and genie::geometry::GeomVolSelectorBasic::TrimSegment().

fPathString

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

full path names

Definition at line 111 of file PathSegmentList.h.

Referenced by Print(), SetPath(), and genie::geometry::GeomVolSelectorBasic::TrimSegment().

fRayDist

Double_t genie::geometry::PathSegment::fRayDist

distance from start of ray

Definition at line 103 of file PathSegmentList.h.

Referenced by DoCrossCheck(), operator<, PathSegment(), Print(), SetEnter(), SetEnter(), genie::geometry::GeomVolSelectorFiducial::TrimSegment(), and genie::geometry::GeomVolSelectorRockBox::TrimSegment().

fStepLength

Double_t genie::geometry::PathSegment::fStepLength

total step size in volume

Definition at line 104 of file PathSegmentList.h.

Referenced by DoCrossCheck(), GetPosition(), PathSegment(), Print(), and SetStep().

fStepRangeSet

StepRangeSet genie::geometry::PathSegment::fStepRangeSet

collection of {steplo,stephi} pairs

Definition at line 113 of file PathSegmentList.h.

Referenced by GetPosition(), GetSummedStepRange(), IsTrimmedEmpty(), Print(), SetStep(), genie::geometry::ROOTGeomAnalyzer::SwimOnce(), genie::geometry::GeomVolSelectorBasic::TrimSegment(), genie::geometry::GeomVolSelectorFiducial::TrimSegment(), and genie::geometry::GeomVolSelectorRockBox::TrimSegment().

fVolume

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

ref only ptr to TGeoVolume

Definition at line 105 of file PathSegmentList.h.

Referenced by genie::geometry::ROOTGeomAnalyzer::GenerateVertex(), PathSegment(), Print(), SetGeo(), and genie::geometry::GeomVolSelectorBasic::TrimSegment().

---

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

• PathSegmentList.h
• PathSegmentList.cxx