genie::geometry::PathSegmentList Class Reference

Object to be filled with the neutrino path-segments representing geometry volume steps (generally boundary-to-boundary) along with geometry materials. Good for a single starting position and travelling along the direction of the neutrino 4-momentum. More...

#include <PathSegmentList.h>

Public Types
typedef std::list< PathSegment >	PathSegmentV_t
typedef PathSegmentV_t::const_iterator	PathSegVCItr_t
typedef std::map< const TGeoMaterial *, Double_t >	MaterialMap_t
typedef MaterialMap_t::const_iterator	MaterialMapCItr_t

Public Member Functions
	PathSegmentList ()
	PathSegmentList (const PathSegmentList &plist)
	~PathSegmentList ()
void	SetDoCrossCheck (bool doit=true)
void	SetPrintVerbose (bool doit=true)
void	SetAllToZero (void)
void	SetStartInfo (const TVector3 &pos=TVector3(0, 0, 1e37), const TVector3 &dir=TVector3(0, 0, 0))
bool	IsSameStart (const TVector3 &pos, const TVector3 &dir) const
void	AddSegment (const PathSegment &ps)
const TVector3 &	GetDirection () const
const TVector3 &	GetStartPos () const
const PathSegmentV_t &	GetPathSegmentV (void) const
size_t	size (void) const
void	FillMatStepSum (void)
const MaterialMap_t &	GetMatStepSumMap (void) const
void	CrossCheck (double &mxddist, double &mxdstep) const
void	Copy (const PathSegmentList &plist)
PathSegmentList &	operator= (const PathSegmentList &list)
void	Print (ostream &stream) const

Protected Attributes
TVector3	fStartPos
	Record, for future comparison, the path taken. More...
TVector3	fDirection
	direction (in top vol coords) More...
PathSegmentV_t	fSegmentList
	Actual list of segments. More...
MaterialMap_t	fMatStepSum
	Segment list re-evaluated by material for fast lookup of path lengths. More...
bool	fDoCrossCheck
bool	fPrintVerbose

Friends
ostream &	operator<< (ostream &stream, const PathSegmentList &list)

Detailed Description

Object to be filled with the neutrino path-segments representing geometry volume steps (generally boundary-to-boundary) along with geometry materials. Good for a single starting position and travelling along the direction of the neutrino 4-momentum.

Author

Robert Hatcher rhatc.nosp@m.her@.nosp@m.fnal..nosp@m.gov FNAL

May 26, 2009

Copyright (c) 2003-2020, The GENIE Collaboration For the full text of the license visit http://copyright.genie-mc.org

Definition at line 123 of file PathSegmentList.h.

Member Typedef Documentation

typedef std::map<const TGeoMaterial*,Double_t> genie::geometry::PathSegmentList::MaterialMap_t

Definition at line 147 of file PathSegmentList.h.

typedef MaterialMap_t::const_iterator genie::geometry::PathSegmentList::MaterialMapCItr_t

Definition at line 148 of file PathSegmentList.h.

typedef std::list<PathSegment> genie::geometry::PathSegmentList::PathSegmentV_t

Definition at line 141 of file PathSegmentList.h.

typedef PathSegmentV_t::const_iterator genie::geometry::PathSegmentList::PathSegVCItr_t

Definition at line 142 of file PathSegmentList.h.

Constructor & Destructor Documentation

PathSegmentList::PathSegmentList

(

void

)

Definition at line 188 of file PathSegmentList.cxx.

  : fDoCrossCheck(false), fPrintVerbose(false)
{

}

PathSegmentList::PathSegmentList

(

const PathSegmentList &

plist

)

Definition at line 194 of file PathSegmentList.cxx.

{
  this->Copy(plist);
}

PathSegmentList::~PathSegmentList

(

)

Definition at line 199 of file PathSegmentList.cxx.

{

}

Member Function Documentation

void genie::geometry::PathSegmentList::AddSegment ( const PathSegment &  ps )

inline

Definition at line 136 of file PathSegmentList.h.

{ fSegmentList.push_back(ps); }

void PathSegmentList::Copy

(

const PathSegmentList &

plist

)

Definition at line 245 of file PathSegmentList.cxx.

{
  fSegmentList.clear();
  fMatStepSum.clear();

  // copy the segments
  //vector<PathSegment>::const_iterator pl_iter;
  //for (pl_iter = plist.fSegmentList.begin(); pl_iter != plist.fSegmentList.end(); ++pl_iter) {
  //  this->fSegmentList.push_back( *pl_iter );
  //}

  // other elements
  fStartPos     = plist.fStartPos;
  fDirection    = plist.fDirection;
  fSegmentList  = plist.fSegmentList;
  fMatStepSum   = plist.fMatStepSum;
  fDoCrossCheck = plist.fDoCrossCheck;
  fPrintVerbose = plist.fPrintVerbose;
}

void PathSegmentList::CrossCheck	(	double &	mxddist,
		double &	mxdstep
	)		const

Definition at line 266 of file PathSegmentList.cxx.

{

  double dstep, ddist;
  mxdstep = 0;
  mxddist = 0;
  PathSegmentList::PathSegVCItr_t sitr = fSegmentList.begin();
  PathSegmentList::PathSegVCItr_t sitr_end = fSegmentList.end();
  for ( ; sitr != sitr_end ; ++sitr ) {
    const PathSegment& ps = *sitr;
    ps.DoCrossCheck(fStartPos,ddist,dstep);
    double addist = TMath::Abs(ddist);
    double adstep = TMath::Abs(dstep);
    if ( addist > mxddist ) mxddist = addist;
    if ( adstep > mxdstep ) mxdstep = adstep;
  }

}

void PathSegmentList::FillMatStepSum

(

void

)

Definition at line 229 of file PathSegmentList.cxx.

{
  fMatStepSum.clear();

  PathSegmentList::PathSegVCItr_t sitr = fSegmentList.begin();
  PathSegmentList::PathSegVCItr_t sitr_end = fSegmentList.end();
  for ( ; sitr != sitr_end ; ++sitr ) {
    const PathSegment& ps = *sitr;
    const TGeoMaterial* mat  = ps.fMaterial;
    // use the post-trim limits on how much material is stepped through
    fMatStepSum[mat] += ps.GetSummedStepRange();
  }

}

const TVector3& genie::geometry::PathSegmentList::GetDirection ( ) const

inline

Definition at line 138 of file PathSegmentList.h.

{ return fDirection; }

const MaterialMap_t& genie::geometry::PathSegmentList::GetMatStepSumMap ( void  ) const

inline

Definition at line 151 of file PathSegmentList.h.

{ return fMatStepSum; };

const PathSegmentV_t& genie::geometry::PathSegmentList::GetPathSegmentV ( void  ) const

inline

Definition at line 144 of file PathSegmentList.h.

{ return fSegmentList; }

const TVector3& genie::geometry::PathSegmentList::GetStartPos ( ) const

inline

Definition at line 139 of file PathSegmentList.h.

{ return fStartPos; }

bool PathSegmentList::IsSameStart	(	const TVector3 &	pos,
		const TVector3 &	dir
	)		const

Definition at line 223 of file PathSegmentList.cxx.

{
  return ( this->fStartPos == pos && this->fDirection == dir );
}

PathSegmentList & PathSegmentList::operator=

(

const PathSegmentList &

list

)

Definition at line 442 of file PathSegmentList.cxx.

{
  this->Copy(list);
  return (*this);
}

void PathSegmentList::Print

(

ostream &

stream

)

const

Definition at line 286 of file PathSegmentList.cxx.

{
  stream << "\nPathSegmentList [-]" << endl;
  stream << " start " << pathsegutils::Vec3AsString(&fStartPos)
         << " dir " << pathsegutils::Vec3AsString(&fDirection) << endl;

  double dstep, ddist, mxdstep = 0, mxddist = 0;
  int k = 0, nseg = 0;
  PathSegmentList::PathSegVCItr_t sitr = fSegmentList.begin();
  PathSegmentList::PathSegVCItr_t sitr_end = fSegmentList.end();
  for ( ; sitr != sitr_end ; ++sitr, ++k ) {
    const PathSegment& ps = *sitr;
    ++nseg;
    stream << " [" << setw(4) << k << "] " << ps;
    if ( fDoCrossCheck ) {
      ps.DoCrossCheck(fStartPos,ddist,dstep);
      double addist = TMath::Abs(ddist);
      double adstep = TMath::Abs(dstep);
      if ( addist > mxddist ) mxddist = addist;
      if ( adstep > mxdstep ) mxdstep = adstep;
      stream << " recalc diff"
             << " dist " << std::setw(12) << ddist
             << " step " << std::setw(12) << dstep;
        }
    stream << std::endl;
  }
  if ( nseg == 0 ) stream << " holds no segments." << std::endl;

  if ( fDoCrossCheck )
    stream << "PathSegmentList "
           << " mxddist " << mxddist
           << " mxdstep " << mxdstep
           << std::endl;

  if ( fPrintVerbose ) {
    PathSegmentList::MaterialMapCItr_t mitr     = GetMatStepSumMap().begin();
    PathSegmentList::MaterialMapCItr_t mitr_end = GetMatStepSumMap().end();
    // loop over map to get tgt weight for each material (once)
    // steps outside the geometry may have no assigned material
    for ( ; mitr != mitr_end; ++mitr ) {
      const TGeoMaterial* mat = mitr->first;
      double sumsteps         = mitr->second;
      stream << " fMatStepSum[" << mat->GetName() << "] = " << sumsteps << std::endl;
    }
  }

}

void PathSegmentList::SetAllToZero

(

void

)

Definition at line 205 of file PathSegmentList.cxx.

{
  LOG("PathS", pDEBUG) << "SetAllToZero called";

  this->fStartPos.SetXYZ(0,0,1e37); // clear cache of position/direction
  this->fDirection.SetXYZ(0,0,0);   //
  this->fSegmentList.clear();       // clear the vector
  this->fMatStepSum.clear();        // clear the re-factorized info
}

void genie::geometry::PathSegmentList::SetDoCrossCheck ( bool  doit = true )

inline

Definition at line 130 of file PathSegmentList.h.

{ fDoCrossCheck = doit; }

void genie::geometry::PathSegmentList::SetPrintVerbose ( bool  doit = true )

inline

Definition at line 131 of file PathSegmentList.h.

{ fPrintVerbose = doit; }

void PathSegmentList::SetStartInfo	(	const TVector3 &	pos = TVector3(0,0,1e37),
		const TVector3 &	dir = TVector3(0,0,0)
	)

Definition at line 216 of file PathSegmentList.cxx.

{
  this->fStartPos  = pos;
  this->fDirection = dir;
}

size_t genie::geometry::PathSegmentList::size ( void  ) const

inline

Definition at line 145 of file PathSegmentList.h.

{ return fSegmentList.size(); }

Friends And Related Function Documentation

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

friend

Definition at line 51 of file PathSegmentList.cxx.

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

Member Data Documentation

TVector3 genie::geometry::PathSegmentList::fDirection

protected

direction (in top vol coords)

Definition at line 170 of file PathSegmentList.h.

bool genie::geometry::PathSegmentList::fDoCrossCheck

protected

Definition at line 178 of file PathSegmentList.h.

MaterialMap_t genie::geometry::PathSegmentList::fMatStepSum

protected

Segment list re-evaluated by material for fast lookup of path lengths.

Definition at line 176 of file PathSegmentList.h.

bool genie::geometry::PathSegmentList::fPrintVerbose

protected

Definition at line 179 of file PathSegmentList.h.

PathSegmentV_t genie::geometry::PathSegmentList::fSegmentList

protected

Actual list of segments.

Definition at line 173 of file PathSegmentList.h.

TVector3 genie::geometry::PathSegmentList::fStartPos

protected

Record, for future comparison, the path taken.

starting position (in top vol coords)

Definition at line 169 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