genie::geometry::ROOTGeomAnalyzer Class Reference

A ROOT/GEANT4 geometry driver. More...

#include <ROOTGeomAnalyzer.h>

Inheritance diagram for genie::geometry::ROOTGeomAnalyzer:

Collaboration diagram for genie::geometry::ROOTGeomAnalyzer:

Public Member Functions
	ROOTGeomAnalyzer (string geometry_filename)
	ROOTGeomAnalyzer (TGeoManager *gm)
	ROOTGeomAnalyzer ()
	~ROOTGeomAnalyzer ()
virtual const PDGCodeList &	ListOfTargetNuclei (void)
	implement the GeomAnalyzerI interface
virtual const PathLengthList &	ComputeMaxPathLengths (void)
virtual const PathLengthList &	ComputePathLengths (const TLorentzVector &x, const TLorentzVector &p)
virtual std::vector< std::pair< double, const TGeoMaterial * > >	ComputeMatLengths (const TLorentzVector &x, const TLorentzVector &p)
virtual const TVector3 &	GenerateVertex (const TLorentzVector &x, const TLorentzVector &p, int tgtpdg)
virtual int	GetTargetPdgCode (const TGeoMaterial *const m) const
virtual int	GetTargetPdgCode (const TGeoMixture *const m, int ielement) const
virtual void	SetScannerNPoints (int np)
	set geometry driver's configuration options
virtual void	SetScannerNRays (int nr)
virtual void	SetScannerNParticles (int np)
virtual void	SetScannerFlux (GFluxI *f)
virtual void	SetWeightWithDensity (bool wt)
virtual void	SetMixtureWeightsSum (double sum)
virtual void	SetLengthUnits (double lu)
virtual void	SetDensityUnits (double du)
virtual void	SetMaxPlSafetyFactor (double sf)
virtual void	SetTopVolName (string nm)
virtual void	SetKeepSegPath (bool keep)
virtual void	SetDebugFlags (int flgs)
virtual int	ScannerNPoints (void) const
	retrieve geometry driver's configuration options
virtual int	ScannerNRays (void) const
virtual int	ScannerNParticles (void) const
virtual bool	WeightWithDensity (void) const
virtual double	LengthUnits (void) const
virtual double	DensityUnits (void) const
virtual double	MixtureWeightsSum (void) const
virtual double	MaxPlSafetyFactor (void) const
virtual string	TopVolName (void) const
virtual TGeoManager *	GetGeometry (void) const
virtual bool	GetKeepSegPath (void) const
virtual const PathLengthList &	GetMaxPathLengths (void) const
virtual void	Local2SI (PathLengthList &pl) const
	access to geometry coordinate/unit transforms for validation/test purposes
virtual void	Local2SI (TVector3 &v) const
virtual void	SI2Local (TVector3 &v) const
virtual void	Master2Top (TVector3 &v) const
virtual void	Master2TopDir (TVector3 &v) const
virtual void	Top2Master (TVector3 &v) const
virtual void	Top2MasterDir (TVector3 &v) const
virtual GeomVolSelectorI *	AdoptGeomVolSelector (GeomVolSelectorI *selector)
	configure processing to perform path segment trimming
Public Member Functions inherited from genie::GeomAnalyzerI
virtual	~GeomAnalyzerI ()

Protected Member Functions
virtual void	Initialize (void)
virtual void	CleanUp (void)
virtual void	Load (string geometry_filename)
virtual void	Load (TGeoManager *gm)
virtual void	BuildListOfTargetNuclei (void)
virtual double	GetWeight (const TGeoMaterial *mat, int pdgc)
virtual double	GetWeight (const TGeoMixture *mixt, int pdgc)
virtual double	GetWeight (const TGeoMixture *mixt, int ielement, int pdgc)
virtual void	MaxPathLengthsFluxMethod (void)
virtual void	MaxPathLengthsBoxMethod (void)
virtual bool	GenBoxRay (int indx, TLorentzVector &x4, TLorentzVector &p4)
virtual double	ComputePathLengthPDG (const TVector3 &r, const TVector3 &udir, int pdgc)
virtual void	SwimOnce (const TVector3 &r, const TVector3 &udir)
virtual bool	FindMaterialInCurrentVol (int pdgc)
virtual bool	WillNeverEnter (double step)
virtual double	StepToNextBoundary (void)
virtual double	Step (void)
virtual double	StepUntilEntering (void)
Protected Member Functions inherited from genie::GeomAnalyzerI
	GeomAnalyzerI ()

Protected Attributes
int	fMaterial
	input selected material for vertex generation
TGeoManager *	fGeometry
	input detector geometry
string	fTopVolumeName
	input top vol [other than TGeoManager::GetTopVolume()]
int	fNPoints
	max path length scanner (box method): points/surface [def:200]
int	fNRays
	max path length scanner (box method): rays/point [def:200]
int	fNParticles
	max path length scanner (flux method): particles in [def:10000]
GFluxI *	fFlux
	a flux objects that can be used to scan the max path lengths
bool	fDensWeight
	if true pathlengths are weighted with density [def:true]
double	fLengthScale
	conversion factor: input geometry length units -> meters
double	fDensityScale
	conversion factor: input geometry density units -> kgr/meters^3
double	fMaxPlSafetyFactor
	factor that can multiply the computed max path lengths
double	fMixtWghtSum
	norm of relative weights (<0 if explicit summing required)
TVector3 *	fCurrVertex
	current generated vertex
PathLengthList *	fCurrPathLengthList
	current list of path-lengths
PathLengthList *	fCurrMaxPathLengthList
	current list of max path-lengths
PDGCodeList *	fCurrPDGCodeList
	current list of target nuclei
TGeoVolume *	fTopVolume
	top volume
TGeoHMatrix *	fMasterToTop
	matrix connecting master coordinates to top volume coordinates
bool	fMasterToTopIsIdentity
	is fMasterToTop matrix the identity matrix?
bool	fKeepSegPath
	need to fill path segment "path"
PathSegmentList *	fCurrPathSegmentList
	current list of path-segments
GeomVolSelectorI *	fGeomVolSelector
	optional path seg trimmer (owned)
TVector3	fGenBoxRayPos
TVector3	fGenBoxRayDir
int	fiface
int	fipoint
int	firay
bool	fnewpnt
double	fdx
double	fdy
double	fdz
double	fox
double	foy
double	foz
	top vol size/origin (top vol units)
double	fmxddist
double	fmxdstep
	max errors in pathsegmentlist
int	fDebugFlags

Detailed Description

A ROOT/GEANT4 geometry driver.

Author

Anselmo Meregaglia <anselmo.meregaglia \at cern.ch> ETH Zurich

Costas Andreopoulos <c.andreopoulos \at cern.ch> University of Liverpool

Robert Hatcher <rhatcher \at fnal.gov> Fermilab

Created:\n May 24, 2005

License:\n Copyright (c) 2003-2025, The GENIE Collaboration

For the full text of the license visit http://copyright.genie-mc.org

Definition at line 54 of file ROOTGeomAnalyzer.h.

Constructor & Destructor Documentation

ROOTGeomAnalyzer() [1/3]

ROOTGeomAnalyzer::ROOTGeomAnalyzer

(

string

geometry_filename

)

Constructor from a geometry file

Definition at line 71 of file ROOTGeomAnalyzer.cxx.

  : GeomAnalyzerI()
{
///
/// Constructor from a geometry file
///
  LOG("GROOTGeom", pDEBUG)
    << "ROOTGeomAnalyzer ctor \"" << geometry_filename << "\"";
  this->Initialize();
  this->Load(geometry_filename);
}

References genie::GeomAnalyzerI::GeomAnalyzerI(), Initialize(), Load(), LOG, and pDEBUG.

ROOTGeomAnalyzer() [2/3]

ROOTGeomAnalyzer::ROOTGeomAnalyzer

(

TGeoManager *

gm

)

Constructor from a TGeoManager

Definition at line 84 of file ROOTGeomAnalyzer.cxx.

  : GeomAnalyzerI()
{
///
/// Constructor from a TGeoManager
///
  LOG("GROOTGeom", pDEBUG)
    << "ROOTGeomAnalyzer ctor passed TGeoManager*";
  this->Initialize();
  this->Load(gm);
}

References genie::GeomAnalyzerI::GeomAnalyzerI(), Initialize(), Load(), LOG, and pDEBUG.

ROOTGeomAnalyzer() [3/3]

genie::geometry::ROOTGeomAnalyzer::ROOTGeomAnalyzer ( )

inline

Definition at line 59 of file ROOTGeomAnalyzer.h.

: GeomAnalyzerI() { ; } // used ONLY for derived class overloading

References genie::GeomAnalyzerI::GeomAnalyzerI().

~ROOTGeomAnalyzer()

ROOTGeomAnalyzer::~ROOTGeomAnalyzer

(

)

Definition at line 97 of file ROOTGeomAnalyzer.cxx.

{
  this->CleanUp();
 
  if ( fmxddist > 0 || fmxdstep > 0 )
    LOG("GROOTGeom",pNOTICE)
      << "ROOTGeomAnalyzer "
      << " mxddist " << fmxddist
      << " mxdstep " << fmxdstep;
}

References CleanUp(), fmxddist, fmxdstep, LOG, and pNOTICE.

Member Function Documentation

AdoptGeomVolSelector()

virtual GeomVolSelectorI * genie::geometry::ROOTGeomAnalyzer::AdoptGeomVolSelector ( GeomVolSelectorI * selector )

inlinevirtual

configure processing to perform path segment trimming

take ownership, return old

Definition at line 121 of file ROOTGeomAnalyzer.h.

  { std::swap(selector,fGeomVolSelector); return selector; }

References fGeomVolSelector.

Referenced by CreateFidSelection(), and CreateRockBoxSelection().

BuildListOfTargetNuclei()

void ROOTGeomAnalyzer::BuildListOfTargetNuclei ( void )

protectedvirtual

Determine possible target PDG codes. Note: If one is using a top volume other than the master level then the final list might include PDG codes that will never be seen during swimming through the volumes if those code are only found in materials outside the top volume.

Definition at line 806 of file ROOTGeomAnalyzer.cxx.

{
/// Determine possible target PDG codes.
/// Note:  If one is using a top volume other than the master level
/// then the final list might include PDG codes that will never
/// be seen during swimming through the volumes if those code are only
/// found in materials outside the top volume.
 
  fCurrPDGCodeList = new PDGCodeList;
 
  if (!fGeometry) {
    LOG("GROOTGeom", pFATAL) << "No ROOT geometry is loaded!!";
    exit(1);
  }
 
  TObjArray * volume_list = fGeometry->GetListOfVolumes();
  if (!volume_list) {
     LOG("GROOTGeom", pERROR)
        << "Null list of geometry volumes. Can not find build target list!";
     return;
  }
 
  std::set<Int_t> seen_mat; // list of materials we've already handled
  std::vector<TGeoVolume*> volvec; //RWH
 
  int numVol = volume_list->GetEntries();
#ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
  LOG("GROOTGeom", pDEBUG) << "Number of volumes found: " << numVol;
#endif
 
  for (int ivol = 0; ivol < numVol; ivol++) {
      TGeoVolume * volume = dynamic_cast <TGeoVolume *>(volume_list->At(ivol));
      if (!volume) {
         LOG("GROOTGeom", pWARN)
           << "Got a null geometry volume!! Skiping current list element";
         continue;
      }
      TGeoMaterial * mat = volume->GetMedium()->GetMaterial();
 
      // shortcut if we've already seen this material
      Int_t mat_indx = mat->GetIndex();
      if ( seen_mat.find(mat_indx) != seen_mat.end() ) continue;
      seen_mat.insert(mat_indx);
      volvec.push_back(volume); //RWH
 
      if (mat->IsMixture()) {
         TGeoMixture * mixt = dynamic_cast <TGeoMixture*> (mat);
         int Nelements = mixt->GetNelements();
         for (int i=0; i<Nelements; i++) {
            int ion_pdgc = this->GetTargetPdgCode(mixt,i);
            fCurrPDGCodeList->push_back(ion_pdgc);
         }
      } else {
          int ion_pdgc = this->GetTargetPdgCode(mat);
          fCurrPDGCodeList->push_back(ion_pdgc);
      }
  }
  // sort the list
  // we don't calculate this list but once per geometry and a sorted
  // list is easier to read so this doesn't cost much
  std::sort(fCurrPDGCodeList->begin(),fCurrPDGCodeList->end());
 
}

References fCurrPDGCodeList, fGeometry, GetTargetPdgCode(), LOG, pDEBUG, pERROR, pFATAL, and pWARN.

Referenced by Load().

CleanUp()

void ROOTGeomAnalyzer::CleanUp ( void )

protectedvirtual

Definition at line 718 of file ROOTGeomAnalyzer.cxx.

{
  LOG("GROOTGeom", pNOTICE) << "Cleaning up...";
 
  if ( fCurrPathSegmentList   ) delete fCurrPathSegmentList;
  if ( fCurrPathLengthList    ) delete fCurrPathLengthList;
  if ( fCurrMaxPathLengthList ) delete fCurrMaxPathLengthList;
  if ( fCurrPDGCodeList       ) delete fCurrPDGCodeList;
  if ( fMasterToTop           ) delete fMasterToTop;
}

References fCurrMaxPathLengthList, fCurrPathLengthList, fCurrPathSegmentList, fCurrPDGCodeList, fMasterToTop, LOG, and pNOTICE.

Referenced by ~ROOTGeomAnalyzer().

ComputeMatLengths()

std::vector< std::pair< double, const TGeoMaterial * > > ROOTGeomAnalyzer::ComputeMatLengths ( const TLorentzVector & x, const TLorentzVector & p )

virtual

Definition at line 207 of file ROOTGeomAnalyzer.cxx.

{
 
  // if trimming configure with neutrino ray's info
  if ( fGeomVolSelector ) {
    fGeomVolSelector->SetCurrentRay(x,p);
    fGeomVolSelector->SetSI2Local(1/this->LengthUnits());
  }
 
  TVector3 udir = p.Vect().Unit(); // unit vector along direction
  TVector3 pos = x.Vect();         // initial position
  this->SI2Local(pos);             // SI -> curr geom units
 
  if (!fMasterToTopIsIdentity) {
    this->Master2Top(pos);         // transform position (master -> top)
    this->Master2TopDir(udir);     // transform direction (master -> top)
  }
 
  this->SwimOnce(pos,udir);
 
  std::vector<std::pair<double, const TGeoMaterial*>> MatLengthList;
 
  const PathSegmentList::PathSegmentV_t& segments = fCurrPathSegmentList->GetPathSegmentV();
 
  PathSegmentList::PathSegVCItr_t sitr;
  for ( sitr = segments.begin(); sitr != segments.end(); ++sitr) {
    const PathSegment& seg = *sitr;
    double pl = seg.GetSummedStepRange();
    MatLengthList.push_back(std::make_pair(pl,seg.fMaterial));
  }
 
  return MatLengthList;
}

References fCurrPathSegmentList, fGeomVolSelector, fMasterToTopIsIdentity, genie::geometry::PathSegment::fMaterial, genie::geometry::PathSegment::GetSummedStepRange(), LengthUnits(), Master2Top(), Master2TopDir(), SI2Local(), and SwimOnce().

ComputeMaxPathLengths()

const PathLengthList & ROOTGeomAnalyzer::ComputeMaxPathLengths ( void )

virtual

Computes the maximum path lengths for all materials in the input geometry. The computed path lengths are in SI units (kgr/m^2, if density weighting is enabled)

Implements genie::GeomAnalyzerI.

Definition at line 118 of file ROOTGeomAnalyzer.cxx.

{
/// Computes the maximum path lengths for all materials in the input
/// geometry. The computed path lengths are in SI units (kgr/m^2, if
/// density weighting is enabled)
 
  LOG("GROOTGeom", pNOTICE)
     << "Computing the maximum path lengths for all materials";
 
  if (!fGeometry) {
      LOG("GROOTGeom", pFATAL) << "No ROOT geometry is loaded!!";
      exit(1);
  }
 
  //-- initialize max path lengths
  fCurrMaxPathLengthList->SetAllToZero();
 
  //-- select maximum path length calculation method
  if ( fFlux ) {
    this->MaxPathLengthsFluxMethod();
    // clear any accumulated exposure accounted generated
    // while exploring the geometry
    fFlux->Clear("CycleHistory");
  } else {
    this->MaxPathLengthsBoxMethod();
  }
 
  return *fCurrMaxPathLengthList;
}

References fCurrMaxPathLengthList, fFlux, fGeometry, LOG, MaxPathLengthsBoxMethod(), MaxPathLengthsFluxMethod(), pFATAL, and pNOTICE.

Referenced by main().

ComputePathLengthPDG()

double ROOTGeomAnalyzer::ComputePathLengthPDG ( const TVector3 & r, const TVector3 & udir, int pdgc )

protectedvirtual

Compute the path length for the material with pdg-code = pdc, staring from the input position r (top vol coord & units) and moving along the direction of the unit vector udir (top vol coord).

Definition at line 1353 of file ROOTGeomAnalyzer.cxx.

{
/// Compute the path length for the material with pdg-code = pdc, staring
/// from the input position r (top vol coord & units) and moving along the
/// direction of the unit vector udir (top vol coord).
 
  double pl = 0; // path-length (x density, if density-weighting is ON)
 
  this->SwimOnce(r0,udir);
 
  double step   = 0;
  double weight = 0;
 
  //  const TGeoVolume   * vol = 0;
  //  const TGeoMedium   * med = 0;
  const TGeoMaterial * mat = 0;
 
  // loop over independent materials, which is shorter or equal to # of volumes
  PathSegmentList::MaterialMapCItr_t itr     =
    fCurrPathSegmentList->GetMatStepSumMap().begin();
  PathSegmentList::MaterialMapCItr_t itr_end =
    fCurrPathSegmentList->GetMatStepSumMap().end();
  for ( ; itr != itr_end; ++itr ) {
    mat  = itr->first;
    if ( ! mat ) continue;  // segment outside geometry has no material
    step = itr->second;
    weight = this->GetWeight(mat,pdgc);
    pl += (step*weight);
  }
 
#ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
  LOG("GROOTGeom", pDEBUG)
       << "PathLength[" << pdgc << "] = " << pl << " in curr geom units";
#endif
 
  return pl;
}

References fCurrPathSegmentList, GetWeight(), LOG, pDEBUG, r0, and SwimOnce().

Referenced by ComputePathLengths(), and GenerateVertex().

ComputePathLengths()

const PathLengthList & ROOTGeomAnalyzer::ComputePathLengths ( const TLorentzVector & x, const TLorentzVector & p )

virtual

Computes the path-length within each detector material for a neutrino starting from point x (master coord) and travelling along the direction of p (master coord). The computed path lengths are in SI units (kgr/m^2, if density weighting is enabled)

Implements genie::GeomAnalyzerI.

Definition at line 149 of file ROOTGeomAnalyzer.cxx.

{
/// Computes the path-length within each detector material for a
/// neutrino starting from point x (master coord) and travelling along
/// the direction of p (master coord).
/// The computed path lengths are in SI units (kgr/m^2, if density
/// weighting is enabled)
 
  //LOG("GROOTGeom", pDEBUG)
  //     << "Computing path-lengths for the input neutrino";
 
#ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
  LOG("GROOTGeom", pDEBUG)
       << "\nInput nu: 4p (GeV) = " << utils::print::P4AsShortString(&p)
       << ", 4x (m,s) = " << utils::print::X4AsString(&x);
#endif
 
  // if trimming configure with neutrino ray's info
  if ( fGeomVolSelector ) {
    fGeomVolSelector->SetCurrentRay(x,p);
    fGeomVolSelector->SetSI2Local(1/this->LengthUnits());
  }
 
  TVector3 udir = p.Vect().Unit(); // unit vector along direction
  TVector3 pos = x.Vect();         // initial position
  this->SI2Local(pos);             // SI -> curr geom units
 
  if (!fMasterToTopIsIdentity) {
    this->Master2Top(pos);         // transform position (master -> top)
    this->Master2TopDir(udir);     // transform direction (master -> top)
  }
 
  // reset current list of path-lengths
  fCurrPathLengthList->SetAllToZero();
 
  //loop over materials & compute the path-length
  vector<int>::iterator itr;
  for (itr=fCurrPDGCodeList->begin();itr!=fCurrPDGCodeList->end();itr++) {
 
    int pdgc = *itr;
 
    Double_t pl = this->ComputePathLengthPDG(pos,udir,pdgc);
    fCurrPathLengthList->AddPathLength(pdgc,pl);
 
#ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
    LOG("GROOTGeom", pINFO)
      <<"Calculated path length for material: " << pdgc << " = " << pl;
#endif
 
  } // loop over materials
 
  this->Local2SI(*fCurrPathLengthList); // curr geom units -> SI
 
  return *fCurrPathLengthList;
}

References ComputePathLengthPDG(), fCurrPathLengthList, fCurrPDGCodeList, fGeomVolSelector, fMasterToTopIsIdentity, LengthUnits(), Local2SI(), LOG, Master2Top(), Master2TopDir(), genie::utils::print::P4AsShortString(), pDEBUG, pINFO, SI2Local(), and genie::utils::print::X4AsString().

Referenced by main(), MaxPathLengthsBoxMethod(), and MaxPathLengthsFluxMethod().

DensityUnits()

virtual double genie::geometry::ROOTGeomAnalyzer::DensityUnits ( void ) const

inlinevirtual

Definition at line 101 of file ROOTGeomAnalyzer.h.

{ return fDensityScale;      }

References fDensityScale.

Referenced by Local2SI().

FindMaterialInCurrentVol()

bool ROOTGeomAnalyzer::FindMaterialInCurrentVol ( int pdgc )

protectedvirtual

Definition at line 1603 of file ROOTGeomAnalyzer.cxx.

{
  TGeoVolume * vol = fGeometry -> GetCurrentVolume();
  if(vol) {
    TGeoMaterial * mat = vol->GetMedium()->GetMaterial();
    if(mat->IsMixture()) {
      TGeoMixture * mixt = dynamic_cast <TGeoMixture*> (mat);
      for(int i = 0; i < mixt->GetNelements(); i++) {
         int pdg = this->GetTargetPdgCode(mixt, i);
         if(tgtpdg == pdg) return true;
      }
    } else {
       int pdg = this->GetTargetPdgCode(mat);
       if(tgtpdg == pdg) return true;
    }
  } else {
     LOG("GROOTGeom", pWARN) << "Current volume is null!";
     return false;
  }
  return false;
}

References fGeometry, GetTargetPdgCode(), LOG, and pWARN.

Referenced by GenerateVertex().

GenBoxRay()

bool ROOTGeomAnalyzer::GenBoxRay ( int indx, TLorentzVector & x4, TLorentzVector & p4 )

protectedvirtual

Generate points in the geometry's bounding box and for each point generate random rays – a pseudo-flux – in master coordinates and SI units

Definition at line 1173 of file ROOTGeomAnalyzer.cxx.

{
/// Generate points in the geometry's bounding box and for each point generate
/// random rays -- a pseudo-flux -- in master coordinates and SI units
 
  firay++;
  fnewpnt = false;
 
  // first time through ... special case
  if ( indx == 0 )           { fiface = 0; fipoint = 0; firay = 0; fnewpnt = true; }
 
  if ( firay   >= fNRays   ) {             fipoint++;   firay = 0; fnewpnt = true; }
  if ( fipoint >= fNPoints ) { fiface++;   fipoint = 0; firay = 0; fnewpnt = true; }
  if ( fiface  >= 3 ) {
    LOG("GROOTGeom",pINFO) << "Box surface scanned: " << indx << " points/rays";
    return false; // signal end
  }
 
  if ( indx == 0 ) {
    // get geometry's bounding box
    //LOG("GROOTGeom", pNOTICE) << "Getting a TGeoBBox enclosing the detector";
    TGeoShape * TS  = fTopVolume->GetShape();
    TGeoBBox *  box = (TGeoBBox *)TS;
    //get box origin and dimensions (in the same units as the geometry)
    fdx = box->GetDX(); // half-length
    fdy = box->GetDY(); // half-length
    fdz = box->GetDZ(); // half-length
    fox = (box->GetOrigin())[0];
    foy = (box->GetOrigin())[1];
    foz = (box->GetOrigin())[2];
 
    LOG("GROOTGeom",pNOTICE)
      << "Box size (GU)   :"
      << " x = " << 2*fdx << ", y = " << 2*fdy << ", z = " << 2*fdz;
    LOG("GROOTGeom",pNOTICE)
      << "Box origin (GU) :"
      << " x = " << fox   << ", y = " << foy   << ", z = " <<   foz;
    LOG("GROOTGeom",pNOTICE)
      << "Will generate [" << fNPoints << "] random points / box surface";
    LOG("GROOTGeom",pNOTICE)
      << "Will generate [" << fNRays   << "] rays / point";
 
#ifdef VALIDATE_CORNERS
    // RWH: test that we know the coordinate transforms for the box corners
    for (int sz = -1; sz <= +1; ++sz) {
      for (int sy = -1; sy <= +1; ++sy) {
        for (int sx = -1; sx <= +1; ++sx) {
          if (sx == 0 || sy == 0 || sz == 0 ) continue;
          TVector3& pos = fGenBoxRayPos;
          pos.SetXYZ(fox+(sx*fdx),foy+(sy*fdy),foz+(sz*fdz));
          TVector3 master(fGenBoxRayPos);
          this->Top2Master(master);   // transform position (top -> master)
          this->Local2SI(master);
          TVector3 pos2(master);
          this->Master2Top(pos2);
          this->SI2Local(pos2);
          LOG("GROOTGeom", pNOTICE)
            << "TopVol corner "
            << " [" << pos[0] << "," << pos[1] << "," << pos[2] << "] "
            << "Master corner "
            << " [" << master[0] << "," << master[1] << "," << master[2] << "] "
            << " top again"
            << " [" << pos2[0] << "," << pos2[1] << "," << pos2[2] << "] ";
        }
      }
    }
#endif
 
  }
 
  RandomGen * rnd = RandomGen::Instance();
 
  double eps = 0.0; //1.0e-12; // tweak to make sure we're inside the box
 
  switch ( fiface ) {
  case 0:
    {
      //top:
      if ( firay == fNRays-1 && fipoint == fNPoints-1 )
        LOG("GROOTGeom",pINFO) << "Box surface scanned: [TOP]";
      fGenBoxRayDir.SetXYZ(-0.5+rnd->RndGeom().Rndm(),
                           -rnd->RndGeom().Rndm(),
                           -0.5+rnd->RndGeom().Rndm());
      if ( fnewpnt )
        fGenBoxRayPos.SetXYZ(fox-fdx+eps+2*(fdx-eps)*rnd->RndGeom().Rndm(),
                             foy+fdy-eps,
                             foz-fdz+eps+2*(fdz-eps)*rnd->RndGeom().Rndm());
      break;
    }
  case 1:
    {
      //left:
      if ( firay == fNRays-1 && fipoint == fNPoints-1 )
        LOG("GROOTGeom",pINFO) << "Box surface scanned: [LEFT]";
      fGenBoxRayDir.SetXYZ(rnd->RndGeom().Rndm(),
                           -0.5+rnd->RndGeom().Rndm(),
                           -0.5+rnd->RndGeom().Rndm());
      if ( fnewpnt )
        fGenBoxRayPos.SetXYZ(fox-fdx+eps,
                             foy-fdy+eps+2*(fdy-eps)*rnd->RndGeom().Rndm(),
                             foz-fdz+eps+2*(fdz-eps)*rnd->RndGeom().Rndm());
      break;
    }
  case 2:
    {
      //front: (really what I, RWH, would call the back)
      if ( firay == fNRays-1 && fipoint == fNPoints-1 )
        LOG("GROOTGeom",pINFO) << "Box surface scanned: [FRONT]";
      fGenBoxRayDir.SetXYZ(-0.5+rnd->RndGeom().Rndm(),
                           -0.5+rnd->RndGeom().Rndm(),
                           -rnd->RndGeom().Rndm());
      if ( fnewpnt )
        fGenBoxRayPos.SetXYZ(fox-fdx+eps+2*(fdx-eps)*rnd->RndGeom().Rndm(),
                             foy-fdy+eps+2*(fdy-eps)*rnd->RndGeom().Rndm(),
                             foz+fdz+eps);
      break;
    }
  }
/*
    //bottom:
      pos.SetXYZ(ox-dx+2*dx*rnd->RndGeom().Rndm(),
                 oy-dy,
                 oz-dz+2*dz*rnd->RndGeom().Rndm());
      dir.SetXYZ(-0.5+rnd->RndGeom().Rndm(),
                 rnd->RndGeom().Rndm(),
                 -0.5+rnd->RndGeom().Rndm());
    //right:
      pos.SetXYZ(ox+dx,
                 oy-dy+2*dy*rnd->RndGeom().Rndm(),
                 oz-dz+2*dz*rnd->RndGeom().Rndm());
      dir.SetXYZ(-rnd->RndGeom().Rndm(),
                 -0.5+rnd->RndGeom().Rndm(),
                 -0.5+rnd->RndGeom().Rndm());
    //back:
      pos.SetXYZ(ox-dx+2*dx*rnd->RndGeom().Rndm(),
                 oy-dy+2*dy*rnd->RndGeom().Rndm(),
                 oz-dz);
      dir.SetXYZ(-0.5+rnd->RndGeom().Rndm(),
                 -0.5+rnd->RndGeom().Rndm(),
                  rnd->RndGeom().Rndm());
*/
 
#ifdef RWH_COUNTVOLS
  curface = fiface;
#endif
 
#ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
  if ( fnewpnt )
    LOG("GROOTGeom", pNOTICE)
      << "GenBoxRay(topvol) "
      << " iface " << fiface << " ipoint " << fipoint << " iray " << firay
      << " newpnt " << (fnewpnt?"true":"false")
      << " x " << utils::print::Vec3AsString(&fGenBoxRayPos)
      << " p " << utils::print::P3AsString(&fGenBoxRayDir);
#endif
 
  if ( fnewpnt ) {
    if ( ! fMasterToTopIsIdentity) {
      this->Top2Master(fGenBoxRayPos); // transform position (top -> master)
    }
    this->Local2SI(fGenBoxRayPos);
  }
  this->Top2MasterDir(fGenBoxRayDir);  // transform direction (top -> master)
 
  x4.SetVect(fGenBoxRayPos);
  p4.SetVect(fGenBoxRayDir.Unit());
 
#ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
  LOG("GROOTGeom", pNOTICE)
    << "GenBoxRay(master) "
    << " iface " << fiface << " ipoint " << fipoint << " iray " << firay
    << " newpnt " << (fnewpnt?"true":"false")
    << " x " << utils::print::Vec3AsString(&fGenBoxRayPos)
    << " p " << utils::print::P3AsString(&fGenBoxRayDir);
#endif
 
  return true;
}

References fdx, fdy, fdz, fGenBoxRayDir, fGenBoxRayPos, fiface, fipoint, firay, fMasterToTopIsIdentity, fnewpnt, fNPoints, fNRays, fox, foy, foz, fTopVolume, genie::RandomGen::Instance(), Local2SI(), LOG, Master2Top(), genie::utils::print::P3AsString(), pINFO, pNOTICE, genie::RandomGen::RndGeom(), SI2Local(), Top2Master(), Top2MasterDir(), and genie::utils::print::Vec3AsString().

Referenced by MaxPathLengthsBoxMethod().

GenerateVertex()

const TVector3 & ROOTGeomAnalyzer::GenerateVertex ( const TLorentzVector & x, const TLorentzVector & p, int tgtpdg )

virtual

Generates a random vertex, within the detector material with the input PDG code, for a neutrino starting from point x (master coord) and travelling along the direction of p (master coord).

Implements genie::GeomAnalyzerI.

Definition at line 243 of file ROOTGeomAnalyzer.cxx.

{
/// Generates a random vertex, within the detector material with the input
/// PDG code, for a neutrino starting from point x (master coord) and
/// travelling along the direction of p (master coord).
 
  LOG("GROOTGeom", pNOTICE)
       << "Generating vtx in material: " << tgtpdg
       << " along the input neutrino direction";
 
  int nretry = 0;
  retry:  // goto label in case of abject failure
  nretry++;
 
  // reset current interaction vertex
  fCurrVertex->SetXYZ(0.,0.,0.);
 
#ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
  LOG("GROOTGeom", pDEBUG)
       << "\nInput nu: 4p (GeV) = " << utils::print::P4AsShortString(&p)
       << ", 4x (m,s) = " << utils::print::X4AsString(&x);
#endif
 
  if (!fGeometry) {
      LOG("GROOTGeom", pFATAL) << "No ROOT geometry is loaded!!";
      exit(1);
  }
 
  // calculate the max path length for the selected material starting from
  // x and looking along the direction of p
  TVector3 udir = p.Vect().Unit();
  TVector3 pos = x.Vect();
  this->SI2Local(pos);           // SI -> curr geom units
 
  if (!fMasterToTopIsIdentity) {
    this->Master2Top(pos);       // transform position (master -> top)
    this->Master2TopDir(udir);   // transform direction (master -> top)
  }
 
  double maxwgt_dist = this->ComputePathLengthPDG(pos,udir,tgtpdg);
  if ( maxwgt_dist <= 0 ) {
    LOG("GROOTGeom", pERROR)
     << "The current trajectory does not cross the selected material!!";
    return *fCurrVertex;
  }
 
  // generate random number between 0 and max_dist
  RandomGen * rnd = RandomGen::Instance();
  double genwgt_dist(maxwgt_dist * rnd->RndGeom().Rndm());
 
  LOG("GROOTGeom", pINFO)
    << "Swim mass: Top Vol dir = " << utils::print::P3AsString(&udir)
    << ", pos = " << utils::print::Vec3AsString(&pos);
  LOG("GROOTGeom", pINFO)
     << "Max {L x Density x Weight} given (init,dir) = " << maxwgt_dist;
  LOG("GROOTGeom", pINFO)
       << "Generated 'distance' in selected material = " << genwgt_dist;
#ifdef RWH_DEBUG
  if ( ( fDebugFlags & 0x01 ) ) {
    fCurrPathSegmentList->SetDoCrossCheck(true);       //RWH
    LOG("GROOTGeom", pINFO) << *fCurrPathSegmentList;  //RWH
    double mxddist = 0, mxdstep = 0;
    fCurrPathSegmentList->CrossCheck(mxddist,mxdstep);
    fmxddist = TMath::Max(fmxddist,mxddist);
    fmxdstep = TMath::Max(fmxdstep,mxdstep);
  }
#endif
 
  // compute the pdg weight for each material just once, then use a stl map
  PathSegmentList::MaterialMap_t wgtmap;
  PathSegmentList::MaterialMapCItr_t mitr     =
    fCurrPathSegmentList->GetMatStepSumMap().begin();
  PathSegmentList::MaterialMapCItr_t mitr_end =
    fCurrPathSegmentList->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 wgt = ( mat ) ? this->GetWeight(mat,tgtpdg) : 0;
    wgtmap[mat] = wgt;
#ifdef RWH_DEBUG
    if ( ( fDebugFlags & 0x02 ) ) {
      LOG("GROOTGeom", pINFO)
        << " wgtmap[" << mat->GetName() << "] pdg " << tgtpdg << " wgt " << Form("%.6f",wgt);
    }
#endif
  }
 
  // walk down the path to pick the vertex
  const genie::geometry::PathSegmentList::PathSegmentV_t& segments =
    fCurrPathSegmentList->GetPathSegmentV();
  genie::geometry::PathSegmentList::PathSegVCItr_t sitr;
  double walked = 0;
  for ( sitr = segments.begin(); sitr != segments.end(); ++sitr) {
    const genie::geometry::PathSegment& seg = *sitr;
    const TGeoMaterial* mat = seg.fMaterial;
    double trimmed_step = seg.GetSummedStepRange();
    double wgtstep = trimmed_step * wgtmap[mat];
    double beyond = walked + wgtstep;
#ifdef RWH_DEBUG
    if ( ( fDebugFlags & 0x04 ) ) {
      LOG("GROOTGeom", pINFO)
        << " beyond " << beyond << " genwgt_dist " << genwgt_dist
        << " trimmed_step " << trimmed_step << " wgtstep " << wgtstep;
    }
#endif
    if ( beyond > genwgt_dist ) {
      // the end of this segment is beyond our generation point
#ifdef RWH_DEBUG
      if ( ( fDebugFlags & 0x08 ) ) {
        LOG("GROOTGeom", pINFO)
          << "Choose vertex pos walked=" << walked
          << " beyond=" << beyond
          << " wgtstep " << wgtstep
          << " ( " << trimmed_step << "*" << wgtmap[mat] << ")"
          << " look for " << genwgt_dist
          << " in " << seg.fVolume->GetName() << " "
          << mat->GetName();
      }
#endif
      // choose a vertex in this segment (possibly multiple steps)
      double frac = ( genwgt_dist - walked ) / wgtstep;
      if ( frac > 1.0 ) {
        LOG("GROOTGeom", pWARN)
          << "Hey, frac = " << frac << " ( > 1.0 ) "
          << genwgt_dist << " " << walked << " " << wgtstep;
      }
      pos = seg.GetPosition(frac);
      fGeometry -> SetCurrentPoint (pos[0],pos[1],pos[2]);
      fGeometry -> FindNode();
      LOG("GROOTGeom", pINFO)
        << "Choose vertex position in " << seg.fVolume->GetName() << " "
         << utils::print::Vec3AsString(&pos);
      break;
    }
    walked = beyond;
  }
 
  LOG("GROOTGeom", pNOTICE)
     << "The vertex was placed in volume: "
     << fGeometry->GetCurrentVolume()->GetName()
     << ", path: " << fGeometry->GetPath();
 
  // warn for any volume overshoots
  bool ok = this->FindMaterialInCurrentVol(tgtpdg);
  if (!ok) {
    LOG("GROOTGeom", pWARN)
       << "Geometry volume was probably overshot";
    LOG("GROOTGeom", pWARN)
       << "No material with code = " << tgtpdg << " could be found at genwgt_dist="
       << genwgt_dist << " (maxwgt_dist=" << maxwgt_dist << ")";
    if ( nretry < 10 ) {
      LOG("GROOTGeom", pWARN)
        << "retry placing vertex";
      goto retry;  // yeah, I know! MyBad.
    }
  }
 
  if (!fMasterToTopIsIdentity) {
     this->Top2Master(pos); // transform position (top -> master)
  }
 
  this->Local2SI(pos);   // curr geom units -> SI
 
  fCurrVertex->SetXYZ(pos[0],pos[1],pos[2]);
 
#ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
  LOG("GROOTGeom", pDEBUG)
      << "Vtx (m) = " << utils::print::Vec3AsString(&pos);
#endif
 
  return *fCurrVertex;
}

References ComputePathLengthPDG(), fCurrPathSegmentList, fCurrVertex, fDebugFlags, fGeometry, FindMaterialInCurrentVol(), fMasterToTopIsIdentity, genie::geometry::PathSegment::fMaterial, fmxddist, fmxdstep, genie::geometry::PathSegment::fVolume, genie::geometry::PathSegment::GetPosition(), genie::geometry::PathSegment::GetSummedStepRange(), GetWeight(), genie::RandomGen::Instance(), Local2SI(), LOG, Master2Top(), Master2TopDir(), genie::utils::print::P3AsString(), genie::utils::print::P4AsShortString(), pDEBUG, pERROR, pFATAL, pINFO, pNOTICE, pWARN, genie::RandomGen::RndGeom(), SI2Local(), Top2Master(), genie::utils::print::Vec3AsString(), and genie::utils::print::X4AsString().

Referenced by main().

GetGeometry()

virtual TGeoManager * genie::geometry::ROOTGeomAnalyzer::GetGeometry ( void ) const

inlinevirtual

Definition at line 105 of file ROOTGeomAnalyzer.h.

{ return fGeometry;          }

References fGeometry.

Referenced by GetGeometry(), and main().

GetKeepSegPath()

virtual bool genie::geometry::ROOTGeomAnalyzer::GetKeepSegPath ( void ) const

inlinevirtual

Definition at line 106 of file ROOTGeomAnalyzer.h.

{ return fKeepSegPath;       }

References fKeepSegPath.

GetMaxPathLengths()

virtual const PathLengthList & genie::geometry::ROOTGeomAnalyzer::GetMaxPathLengths ( void ) const

inlinevirtual

Definition at line 107 of file ROOTGeomAnalyzer.h.

{ return *fCurrMaxPathLengthList; } // call only after ComputeMaxPathLengths() has been called

References fCurrMaxPathLengthList.

Referenced by main().

GetTargetPdgCode() [1/2]

int ROOTGeomAnalyzer::GetTargetPdgCode ( const TGeoMaterial *const m ) const

virtual

Definition at line 871 of file ROOTGeomAnalyzer.cxx.

{
  int A = TMath::Nint(m->GetA());
  int Z = TMath::Nint(m->GetZ());
 
  int pdgc = pdg::IonPdgCode(A,Z);
 
  return pdgc;
}

References genie::pdg::IonPdgCode().

Referenced by BuildListOfTargetNuclei(), FindMaterialInCurrentVol(), GetWeight(), GetWeight(), and GetWeight().

GetTargetPdgCode() [2/2]

int ROOTGeomAnalyzer::GetTargetPdgCode ( const TGeoMixture *const m, int ielement ) const

virtual

Definition at line 882 of file ROOTGeomAnalyzer.cxx.

{
  int A = TMath::Nint(m->GetAmixt()[ielement]);
  int Z = TMath::Nint(m->GetZmixt()[ielement]);
 
  int pdgc = pdg::IonPdgCode(A,Z);
 
  return pdgc;
}

References genie::pdg::IonPdgCode().

GetWeight() [1/3]

double ROOTGeomAnalyzer::GetWeight ( const TGeoMaterial * mat, int pdgc )

protectedvirtual

Get the weight of the input material. Return the weight only if the material's pdg code matches the input code. If the material is found to be a mixture, call the corresponding method for mixtures. Weight is in the curr geom density units.

Definition at line 894 of file ROOTGeomAnalyzer.cxx.

{
/// Get the weight of the input material.
/// Return the weight only if the material's pdg code matches the input code.
/// If the material is found to be a mixture, call the corresponding method
/// for mixtures.
/// Weight is in the curr geom density units.
 
  if (!mat) {
    LOG("GROOTGeom", pERROR) << "Null input material. Return weight = 0.";
    return 0;
  }
 
  bool exists = fCurrPDGCodeList->ExistsInPDGCodeList(pdgc);
  if (!exists) {
    LOG("GROOTGeom", pERROR) << "Target doesn't exist. Return weight = 0.";
    return 0;
  }
 
#ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
  LOG("GROOTGeom",pDEBUG)
       << "Curr. material: A/Z = " << mat->GetA() << " / " << mat->GetZ();
#endif
 
  // if the input material is a mixture, get a the sum of weights for
  // all matching elements
  double weight = 0.;
  if (mat->IsMixture()) {
    const TGeoMixture * mixt = dynamic_cast <const TGeoMixture*> (mat);
    if (!mixt) {
     LOG("GROOTGeom", pERROR) << "Null input mixture. Return weight = 0.";
     return 0;
    }
#ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
    LOG("GROOTGeom", pDEBUG)
      << "Material : " << mat->GetName()
          << " is a mixture with " << mixt->GetNelements() << " elements";
#endif
    // loop over elements & sum weights of matching elements
    weight = this->GetWeight(mixt,pdgc);
    return weight;
  } // is mixture?
 
  // pure material
  int ion_pdgc = this->GetTargetPdgCode(mat);
  if (ion_pdgc != pdgc) return 0.;
 
  if (this->WeightWithDensity())
    weight = mat->GetDensity(); // material density (curr geom units)
  else
    weight = 1.0;
 
#ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
  LOG("GROOTGeom", pDEBUG)
       << "Weight[mat:" << mat->GetName() << "] = " << weight;
#endif
 
  return weight;
}

References fCurrPDGCodeList, GetTargetPdgCode(), GetWeight(), LOG, pDEBUG, pERROR, and WeightWithDensity().

Referenced by ComputePathLengthPDG(), GenerateVertex(), GetWeight(), and GetWeight().

GetWeight() [2/3]

double ROOTGeomAnalyzer::GetWeight ( const TGeoMixture * mixt, int ielement, int pdgc )

protectedvirtual

Get the weight of the input ith element of the input material. Return the weight only if the element's pdg code matches the input code. Weight is in the curr geom density units.

Definition at line 1006 of file ROOTGeomAnalyzer.cxx.

{
/// Get the weight of the input ith element of the input material.
/// Return the weight only if the element's pdg code matches the input code.
/// Weight is in the curr geom density units.
 
//  int ion_pdgc = this->GetTargetPdgCode(mixt->GetElement(ielement));
  int ion_pdgc = this->GetTargetPdgCode(mixt, ielement);
  if (ion_pdgc != pdgc) return 0.;
 
  double d = mixt->GetDensity();         // mixture density (curr geom units)
  double w = mixt->GetWmixt()[ielement]; // relative proportion by mass
 
  double wtot = this->MixtureWeightsSum();
 
  // <0 forces explicit calculation of relative proportion normalization
  if (wtot < 0) {
    wtot = 0;
    for (int i = 0; i < mixt->GetNelements(); i++) {
      wtot += (mixt->GetWmixt()[i]);
    }
  }
  assert(wtot>0);
 
  w /= wtot;
  double weight = d*w;
 
  return weight;
}

References GetTargetPdgCode(), and MixtureWeightsSum().

GetWeight() [3/3]

double ROOTGeomAnalyzer::GetWeight ( const TGeoMixture * mixt, int pdgc )

protectedvirtual

Loop over the mixture elements, find the one matching the input pdgc and return its weight. Weight is in the curr geom density units.

Definition at line 955 of file ROOTGeomAnalyzer.cxx.

{
/// Loop over the mixture elements, find the one matching the input pdgc
/// and  return its weight.
/// Weight is in the curr geom density units.
 
  double weight = 0;
 
  int nm = 0;
  for (int i = 0; i < mixt->GetNelements(); i++) {
     double dw = (this->GetWeight(mixt,i,pdgc));
     if (dw>0) nm++;
     weight += dw;
  }
 
  if (nm>1) {
     for (int j = 0; j < mixt->GetNelements(); j++) {
           LOG("GROOTGeom", pWARN)
              << "[" << j << "] Z = " << mixt->GetZmixt()[j]
              << ", A = " << mixt->GetAmixt()[j]
              << " (pdgc = " << this->GetTargetPdgCode(mixt,j)
              << "), w = " << mixt->GetWmixt()[j];
     }
     LOG("GROOTGeom", pERROR)
        << "Material pdgc = " << pdgc << " appears " << nm
        << " times (>1) in mixture = " << mixt->GetName();
 
     // As of ROOT v6.25.02, isomers with the same (A,Z) are not automatically
     // combined. Thus, having two instances of the same nuclear PDG code in a
     // mixture is no longer an error condition. The code above correctly sums
     // over contributions from matching isomers (which GENIE will not care
     // about). We still print the warning so that the user can double-check
     // the material definition, but the fatal error is now removed.
     // See https://github.com/root-project/root/pull/8556 for details.
     // A new production campaign for SBND has been modestly delayed due to
     // encountering this fatal error after updating ROOT.
     // -- S. Gardiner, 19 June 2023
     //
     //LOG("GROOTGeom", pFATAL)
     //   << "Your geometry must be incorrect - Aborting";
     //exit(1);
  }
 
  // if we are not weighting with the density then the weight=1 if the pdg
  // code was matched for any element of this mixture
  if ( !this->WeightWithDensity() && weight>0. ) weight=1.0;
 
  return weight;
}

References GetTargetPdgCode(), GetWeight(), LOG, pERROR, pWARN, and WeightWithDensity().

Initialize()

void ROOTGeomAnalyzer::Initialize ( void )

protectedvirtual

Definition at line 685 of file ROOTGeomAnalyzer.cxx.

{
  LOG("GROOTGeom", pNOTICE)
                << "Initializing ROOT geometry driver & setting defaults";
 
  fCurrMaxPathLengthList = 0;
  fCurrPathLengthList    = 0;
  fCurrPathSegmentList   = 0;
  fGeomVolSelector       = 0;
  fCurrPDGCodeList       = 0;
  fTopVolume             = 0;
  fTopVolumeName         = "";
  fKeepSegPath           = false;
 
  // some defaults:
  this -> SetScannerNPoints    (200);
  this -> SetScannerNRays      (200);
  this -> SetScannerNParticles (10000);
  this -> SetScannerFlux       (0);
  this -> SetMaxPlSafetyFactor (1.1);
  this -> SetLengthUnits       (genie::units::meter);
  this -> SetDensityUnits      (genie::units::kilogram/genie::units::meter3);
  this -> SetWeightWithDensity (true);
  this -> SetMixtureWeightsSum (-1.);
 
  fMasterToTopIsIdentity = true;
 
  fmxddist = 0;
  fmxdstep = 0;
  fDebugFlags = 0;
}

References fCurrMaxPathLengthList, fCurrPathLengthList, fCurrPathSegmentList, fCurrPDGCodeList, fDebugFlags, fGeomVolSelector, fKeepSegPath, fMasterToTopIsIdentity, fmxddist, fmxdstep, fTopVolume, fTopVolumeName, genie::units::kilogram, LOG, genie::units::meter, genie::units::meter3, pNOTICE, SetDensityUnits(), SetLengthUnits(), SetMaxPlSafetyFactor(), SetMixtureWeightsSum(), SetScannerFlux(), SetScannerNParticles(), SetScannerNPoints(), SetScannerNRays(), and SetWeightWithDensity().

Referenced by ROOTGeomAnalyzer(), and ROOTGeomAnalyzer().

LengthUnits()

virtual double genie::geometry::ROOTGeomAnalyzer::LengthUnits ( void ) const

inlinevirtual

Definition at line 100 of file ROOTGeomAnalyzer.h.

{ return fLengthScale;       }

References fLengthScale.

Referenced by ComputeMatLengths(), ComputePathLengths(), GetGeometry(), Local2SI(), Local2SI(), main(), and SI2Local().

ListOfTargetNuclei()

const PDGCodeList & ROOTGeomAnalyzer::ListOfTargetNuclei ( void )

virtual

implement the GeomAnalyzerI interface

Implements genie::GeomAnalyzerI.

Definition at line 112 of file ROOTGeomAnalyzer.cxx.

{
  return *fCurrPDGCodeList;
}

References fCurrPDGCodeList.

Referenced by Load().

Load() [1/2]

void ROOTGeomAnalyzer::Load ( string geometry_filename )

protectedvirtual

Load the detector geometry from the input ROOT file

Definition at line 730 of file ROOTGeomAnalyzer.cxx.

{
/// Load the detector geometry from the input ROOT file
///
  LOG("GROOTGeom", pNOTICE) << "Loading geometry from: " << filename;
 
  bool is_accessible = ! (gSystem->AccessPathName( filename.c_str() ));
  if (!is_accessible) {
     LOG("GROOTGeom", pFATAL)
       << "The ROOT geometry doesn't exist! Initialization failed!";
     exit(1);
  }
 
// ROOT versions 6.16 - 6.24 defaulted to kG4Units [ugh]
// worse yet the interface for setting it also changed at 6.22
#if ROOT_VERSION_CODE >= ROOT_VERSION(6,16,0)
  if (TGeoManager::GetDefaultUnits() != TGeoManager::kRootUnits) {
  #if ROOT_VERSION_CODE >= ROOT_VERSION(6,22,0)
    TGeoManager::LockDefaultUnits(false);
    TGeoManager::SetDefaultUnits(TGeoManager::kRootUnits);
    TGeoManager::LockDefaultUnits(true);
  #else
    TGeoManager::SetDefaultRootUnits();
  #endif
  }
#endif
 
  TGeoManager * gm = TGeoManager::Import(filename.c_str());
 
  this->Load(gm);
}

References Load(), LOG, pFATAL, and pNOTICE.

Referenced by Load(), ROOTGeomAnalyzer(), and ROOTGeomAnalyzer().

Load() [2/2]

void ROOTGeomAnalyzer::Load ( TGeoManager * gm )

protectedvirtual

Load the detector geometry from the input TGeoManager

Definition at line 763 of file ROOTGeomAnalyzer.cxx.

{
/// Load the detector geometry from the input TGeoManager
 
  LOG("GROOTGeom", pNOTICE)
         << "A TGeoManager is being loaded to the geometry driver";
  fGeometry = gm;
 
  if (!fGeometry) {
    LOG("GROOTGeom", pFATAL) << "Null TGeoManager! Aborting";
  }
  assert(fGeometry);
 
  this->BuildListOfTargetNuclei();
 
  const PDGCodeList & pdglist = this->ListOfTargetNuclei();
 
  fTopVolume             = 0;
  fCurrPathSegmentList   = new PathSegmentList();
  fCurrPathLengthList    = new PathLengthList(pdglist);
  fCurrMaxPathLengthList = new PathLengthList(pdglist);
  fCurrVertex            = new TVector3(0.,0.,0.);
 
  // ask geometry manager for its top volume
  fTopVolume = fGeometry->GetTopVolume();
  if (!fTopVolume) {
      LOG("GROOTGeom", pFATAL) << "Could not get top volume!!!";
  }
  assert(fTopVolume);
 
  // load matrix (identity) of top volume
  fMasterToTop = new TGeoHMatrix(*fGeometry->GetCurrentMatrix());
  fMasterToTopIsIdentity = true;
 
//#define PRINT_MATERIALS
#ifdef PRINT_MATERIALS
  fGeometry->GetListOfMaterials()->Print();
  fGeometry->GetListOfMedia()->Print();
#endif
 
}

References BuildListOfTargetNuclei(), fCurrMaxPathLengthList, fCurrPathLengthList, fCurrPathSegmentList, fCurrVertex, fGeometry, fMasterToTop, fMasterToTopIsIdentity, fTopVolume, ListOfTargetNuclei(), LOG, pFATAL, and pNOTICE.

Local2SI() [1/2]

void ROOTGeomAnalyzer::Local2SI ( PathLengthList & pl ) const

virtual

access to geometry coordinate/unit transforms for validation/test purposes

convert path lengths from current geometry units to SI units

Definition at line 530 of file ROOTGeomAnalyzer.cxx.

{
/// convert path lengths from current geometry units to SI units
///
 
  double scaling_factor = this->LengthUnits();
  if (this->WeightWithDensity()) { scaling_factor *= this->DensityUnits(); }
 
#ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
  LOG("GROOTGeom", pDEBUG)
    << "Scaling path-lengths from local units -> meters "
    << ((this->WeightWithDensity()) ? "* kgr/m^3" : "")
    << " - scale = " << scaling_factor;
#endif
 
  PathLengthList::iterator pliter;
  for(pliter = pl.begin(); pliter != pl.end(); ++pliter)
  {
    int pdgc = pliter->first;
    pl.ScalePathLength(pdgc, scaling_factor);
  }
}

References DensityUnits(), LengthUnits(), LOG, pDEBUG, genie::PathLengthList::ScalePathLength(), and WeightWithDensity().

Referenced by ComputePathLengths(), GenBoxRay(), and GenerateVertex().

Local2SI() [2/2]

void ROOTGeomAnalyzer::Local2SI ( TVector3 & v ) const

virtual

convert position vector from current geometry units to SI units

Definition at line 554 of file ROOTGeomAnalyzer.cxx.

{
/// convert position vector from current geometry units to SI units
///
 
#ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
    LOG("GROOTGeom", pDEBUG)
      << "Position (loc): " << utils::print::Vec3AsString(&vec);
#endif
 
    vec *= (this->LengthUnits());
 
#ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
    LOG("GROOTGeom", pDEBUG)
      << "Position (SI): " << utils::print::Vec3AsString(&vec);
#endif
}

References LengthUnits(), LOG, pDEBUG, and genie::utils::print::Vec3AsString().

Master2Top()

void ROOTGeomAnalyzer::Master2Top ( TVector3 & v ) const

virtual

transform the input position vector from the master volume coordinate system to the specified top volume coordinate system, but not change the units.

Definition at line 591 of file ROOTGeomAnalyzer.cxx.

{
/// transform the input position vector from the master volume coordinate
/// system to the specified top volume coordinate system, but not
/// change the units.
 
#ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
    LOG("GROOTGeom", pDEBUG)
      << "Position (coord:master): " << utils::print::Vec3AsString(&vec);
#endif
 
    Double_t mast[3], top[3];
    vec.GetXYZ(mast);
    fMasterToTop->MasterToLocal(mast, top);
    vec.SetXYZ(top[0], top[1], top[2]);
 
#ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
    LOG("GROOTGeom", pDEBUG)
      << "Position (coord:top): " << utils::print::Vec3AsString(&vec);
#endif
}

References fMasterToTop, LOG, pDEBUG, and genie::utils::print::Vec3AsString().

Referenced by ComputeMatLengths(), ComputePathLengths(), genie::geometry::FidCylinder::ConvertMaster2Top(), genie::geometry::FidSphere::ConvertMaster2Top(), genie::geometry::PlaneParam::ConvertMaster2Top(), GenBoxRay(), and GenerateVertex().

Master2TopDir()

void ROOTGeomAnalyzer::Master2TopDir ( TVector3 & v ) const

virtual

transform the input direction vector from the master volume coordinate system to the specified top volume coordinate system, but not change the units.

Definition at line 614 of file ROOTGeomAnalyzer.cxx.

{
/// transform the input direction vector from the master volume coordinate
/// system to the specified top volume coordinate system, but not
/// change the units.
 
#ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
    LOG("GROOTGeom", pDEBUG)
      << "Direction (coord:master): " << utils::print::Vec3AsString(&vec);
#endif
 
    Double_t mast[3], top[3];
    vec.GetXYZ(mast);
    fMasterToTop->MasterToLocalVect(mast, top);
    vec.SetXYZ(top[0], top[1], top[2]);
 
#ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
    LOG("GROOTGeom", pDEBUG)
      << "Direction (coord:top): " << utils::print::Vec3AsString(&vec);
#endif
}

References fMasterToTop, LOG, pDEBUG, and genie::utils::print::Vec3AsString().

Referenced by ComputeMatLengths(), ComputePathLengths(), genie::geometry::FidCylinder::ConvertMaster2Top(), genie::geometry::PlaneParam::ConvertMaster2Top(), and GenerateVertex().

MaxPathLengthsBoxMethod()

void ROOTGeomAnalyzer::MaxPathLengthsBoxMethod ( void )

protectedvirtual

Generate points in the geometry's bounding box and for each point generate random rays, follow them through the detector and figure out the maximum path length for each material

Definition at line 1101 of file ROOTGeomAnalyzer.cxx.

{
/// Generate points in the geometry's bounding box and for each point
/// generate random rays, follow them through the detector and figure out
/// the maximum path length for each material
 
  LOG("GROOTGeom", pNOTICE)
    << "Computing the maximum path lengths using the BOX method";
#ifdef RWH_COUNTVOLS
  accum_vol_stat = true;
#endif
 
  int  iparticle = 0;
  bool ok = true;
  TLorentzVector nux4;
  TLorentzVector nup4;
 
  PathLengthList::const_iterator pl_iter;
 
  while ( (ok = this->GenBoxRay(iparticle++,nux4,nup4)) ) {
 
    //LOG("GMCJDriver", pNOTICE)
    //  << "\n [-] Generated flux neutrino: "
    //  << "\n  |----o 4-momentum : " << utils::print::P4AsString(&nup4)
    //  << "\n  |----o 4-position : " << utils::print::X4AsString(&nux4);
 
    const PathLengthList & pllst = this->ComputePathLengths(nux4, nup4);
 
    for (pl_iter = pllst.begin(); pl_iter != pllst.end(); ++pl_iter) {
       int    pdgc = pl_iter->first;
       double pl   = pl_iter->second;
 
       if (pl>0) {
          pl *= (this->MaxPlSafetyFactor());
 
          pl = TMath::Max(pl, fCurrMaxPathLengthList->PathLength(pdgc));
          fCurrMaxPathLengthList->SetPathLength(pdgc,pl);
       }
    }
  }
 
  // print out the results
  LOG("GROOTGeom", pDEBUG)
    << "DensWeight \"" << (fDensWeight?"true":"false")
    << "\" MixtWghtSum " << fMixtWghtSum;
  LOG("GROOTGeom", pDEBUG) << "CurrMaxPathLengthList: "
    << *fCurrMaxPathLengthList;
 
#ifdef RWH_COUNTVOLS
  // rwh
  // print statistics for average,rms of number of volumes seen for
  // various rays for each face
  for (int j = 0; j < 6; ++j ) {
    long int ns = nswims[j];
    double   x  = dnvols[j];
    double   x2 = dnvols2[j];
    if ( ns == 0 ) ns = 1;
    double avg = x / (double)ns;
    double rms = TMath::Sqrt((x2/(double)ns) - avg*avg);
    LOG("GROOTGeom", pNOTICE)
      << "RWH: nswim after BOX face " << j << " is " << ns
      << " avg " << avg << " rms " << rms
      << " never " << nnever[j];
  }
  LOG("GROOTGeom", pNOTICE)
    << "RWH: Max PathSegmentList size " << mxsegments;
  accum_vol_stat = false;
#endif
 
}

References ComputePathLengths(), fCurrMaxPathLengthList, fDensWeight, fMixtWghtSum, GenBoxRay(), LOG, MaxPlSafetyFactor(), pDEBUG, and pNOTICE.

Referenced by ComputeMaxPathLengths().

MaxPathLengthsFluxMethod()

void ROOTGeomAnalyzer::MaxPathLengthsFluxMethod ( void )

protectedvirtual

Use the input flux driver to generate "rays", and then follow them through the detector and figure out the maximum path length for each material

Definition at line 1037 of file ROOTGeomAnalyzer.cxx.

{
/// Use the input flux driver to generate "rays", and then follow them through
/// the detector and figure out the maximum path length for each material
 
  LOG("GROOTGeom", pNOTICE)
               << "Computing the maximum path lengths using the FLUX method";
 
  int iparticle = 0;
  PathLengthList::const_iterator pl_iter;
 
  const int nparticles = abs(this->ScannerNParticles());
 
  // if # scanner particles is negative, this signals that the user
  // desires to force rays to have the maximum energy (useful if the
  // volume considered changes size with neutrino energy)
  bool rescale_e = (this->ScannerNParticles() < 0 );
  double emax = fFlux->MaxEnergy();
  if ( rescale_e ) {
    LOG("GROOTGeom", pNOTICE)
      << "max path lengths with FLUX method forcing Enu=" << emax;
  }
 
  while (iparticle < nparticles ) {
 
    bool ok = fFlux->GenerateNext();
    if (!ok) {
       LOG("GROOTGeom", pWARN) << "Couldn't generate a flux neutrino";
       continue;
    }
 
    TLorentzVector   nup4  = fFlux->Momentum();
    if ( rescale_e ) {
      double ecurr = nup4.E();
      if ( ecurr > 0 ) nup4 *= (emax/ecurr);
    }
    const TLorentzVector & nux4  = fFlux->Position();
 
    //LOG("GMCJDriver", pNOTICE)
    //   << "\n [-] Generated flux neutrino: "
    //   << "\n  |----o 4-momentum : " << utils::print::P4AsString(&nup4)
    //   << "\n  |----o 4-position : " << utils::print::X4AsString(&nux4);
 
    const PathLengthList & pl = this->ComputePathLengths(nux4, nup4);
 
    bool enters = false;
 
    for (pl_iter = pl.begin(); pl_iter != pl.end(); ++pl_iter) {
       int    pdgc       = pl_iter->first;
       double pathlength = pl_iter->second;
 
       if ( pathlength > 0 ) {
          pathlength *= (this->MaxPlSafetyFactor());
 
          pathlength = TMath::Max(pathlength, fCurrMaxPathLengthList->PathLength(pdgc));
          fCurrMaxPathLengthList->SetPathLength(pdgc,pathlength);
          enters = true;
       }
    }
    if (enters) iparticle++;
  }
}

References ComputePathLengths(), fCurrMaxPathLengthList, fFlux, LOG, MaxPlSafetyFactor(), pNOTICE, pWARN, and ScannerNParticles().

Referenced by ComputeMaxPathLengths().

MaxPlSafetyFactor()

virtual double genie::geometry::ROOTGeomAnalyzer::MaxPlSafetyFactor ( void ) const

inlinevirtual

Definition at line 103 of file ROOTGeomAnalyzer.h.

{ return fMaxPlSafetyFactor; }

References fMaxPlSafetyFactor.

Referenced by MaxPathLengthsBoxMethod(), and MaxPathLengthsFluxMethod().

MixtureWeightsSum()

virtual double genie::geometry::ROOTGeomAnalyzer::MixtureWeightsSum ( void ) const

inlinevirtual

Definition at line 102 of file ROOTGeomAnalyzer.h.

{ return fMixtWghtSum;       }

References fMixtWghtSum.

Referenced by GetWeight().

ScannerNParticles()

virtual int genie::geometry::ROOTGeomAnalyzer::ScannerNParticles ( void ) const

inlinevirtual

Definition at line 98 of file ROOTGeomAnalyzer.h.

{ return fNParticles;        }

References fNParticles.

Referenced by MaxPathLengthsFluxMethod().

ScannerNPoints()

virtual int genie::geometry::ROOTGeomAnalyzer::ScannerNPoints ( void ) const

inlinevirtual

retrieve geometry driver's configuration options

Definition at line 96 of file ROOTGeomAnalyzer.h.

{ return fNPoints;           }

References fNPoints.

ScannerNRays()

virtual int genie::geometry::ROOTGeomAnalyzer::ScannerNRays ( void ) const

inlinevirtual

Definition at line 97 of file ROOTGeomAnalyzer.h.

{ return fNRays;             }

References fNRays.

SetDebugFlags()

virtual void genie::geometry::ROOTGeomAnalyzer::SetDebugFlags ( int flgs )

inlinevirtual

Definition at line 92 of file ROOTGeomAnalyzer.h.

{ fDebugFlags  = flgs; }

References fDebugFlags.

SetDensityUnits()

void ROOTGeomAnalyzer::SetDensityUnits ( double du )

virtual

Like SetLengthUnits, but for density (default units = kgr/m3)

Definition at line 439 of file ROOTGeomAnalyzer.cxx.

{
/// Like SetLengthUnits, but for density (default units = kgr/m3)
 
  fDensityScale = u / (units::kilogram / units::meter3);
  LOG("GROOTGeom", pNOTICE)
    << "Geometry density units scale factor (geom units -> kgr/m3): "
    << fDensityScale;
}

References fDensityScale, genie::units::kilogram, LOG, genie::units::meter3, and pNOTICE.

Referenced by Initialize().

SetKeepSegPath()

virtual void genie::geometry::ROOTGeomAnalyzer::SetKeepSegPath ( bool keep )

inlinevirtual

Definition at line 91 of file ROOTGeomAnalyzer.h.

{ fKeepSegPath = keep; }

References fKeepSegPath.

SetLengthUnits()

void ROOTGeomAnalyzer::SetLengthUnits ( double lu )

virtual

Use the units of the input geometry, e.g. SetLengthUnits(genie::units::centimeter) GENIE uses the physical system of units (hbar=c=1) almost throughtout so everything is expressed in GeV but when analyzing detector geometries use meters. Setting input geometry units will allow the code to compute the conversion factor. As input, use one of the constants in $GENIE/src/Conventions/Units.h

Definition at line 422 of file ROOTGeomAnalyzer.cxx.

{
/// Use the units of the input geometry,
///    e.g. SetLengthUnits(genie::units::centimeter)
/// GENIE uses the physical system of units (hbar=c=1) almost throughtout
/// so everything is expressed in GeV but when analyzing detector geometries
/// use meters. Setting input geometry units will allow the code to compute
/// the conversion factor.
/// As input, use one of the constants in $GENIE/src/Conventions/Units.h
 
  fLengthScale = u/units::meter;
  LOG("GROOTGeom", pNOTICE)
     << "Geometry length units scale factor (geom units -> m): "
     << fLengthScale;
}

References fLengthScale, LOG, genie::units::meter, and pNOTICE.

Referenced by Initialize().

SetMaxPlSafetyFactor()

void ROOTGeomAnalyzer::SetMaxPlSafetyFactor ( double sf )

virtual

Set a factor that can multiply the computed max path lengths. The maximum path lengths are computed by performing an MC scanning of the input geometry. If you configure the scanner with a low number of points or rays you might understimate the path lengths, so you might want to 'inflate' them a little bit using this method. Do not set this number too high, because the max interaction probability will be grossly overestimated and you would need lots of attempts before getting a flux neutrino to interact...

Definition at line 450 of file ROOTGeomAnalyzer.cxx.

{
/// Set a factor that can multiply the computed max path lengths.
/// The maximum path lengths are computed by performing an MC scanning of
/// the input geometry. If you configure the scanner with a low number of
/// points or rays you might understimate the path lengths, so you might
/// want to 'inflate' them a little bit using this method.
/// Do not set this number too high, because the max interaction probability
/// will be grossly overestimated and you would need lots of attempts before
/// getting a flux neutrino to interact...
 
  fMaxPlSafetyFactor = sf;
 
  LOG("GROOTGeom", pNOTICE)
    << "Max path length safety factor: " << fMaxPlSafetyFactor;
}

References fMaxPlSafetyFactor, LOG, and pNOTICE.

Referenced by Initialize().

SetMixtureWeightsSum()

void ROOTGeomAnalyzer::SetMixtureWeightsSum ( double sum )

virtual

Set it to x, if the relative weight proportions of elements in a mixture add up to x (eg x=1, 100, etc). Set it to a negative value to explicitly compute the correct weight normalization.

Definition at line 468 of file ROOTGeomAnalyzer.cxx.

{
/// Set it to x, if the relative weight proportions of elements in a mixture
/// add up to x (eg x=1, 100, etc). Set it to a negative value to explicitly
/// compute the correct weight normalization.
 
  fMixtWghtSum = sum;
}

References fMixtWghtSum.

Referenced by Initialize().

SetScannerFlux()

virtual void genie::geometry::ROOTGeomAnalyzer::SetScannerFlux ( GFluxI * f )

inlinevirtual

Definition at line 84 of file ROOTGeomAnalyzer.h.

{ fFlux       = f;  } /* flux scanner */

References fFlux.

Referenced by Initialize(), and main().

SetScannerNParticles()

virtual void genie::geometry::ROOTGeomAnalyzer::SetScannerNParticles ( int np )

inlinevirtual

Definition at line 83 of file ROOTGeomAnalyzer.h.

{ fNParticles = np; } /* flux scanner */

References fNParticles.

Referenced by Initialize(), and main().

SetScannerNPoints()

virtual void genie::geometry::ROOTGeomAnalyzer::SetScannerNPoints ( int np )

inlinevirtual

set geometry driver's configuration options

Definition at line 81 of file ROOTGeomAnalyzer.h.

{ fNPoints    = np; } /* box  scanner */

References fNPoints.

Referenced by Initialize(), and main().

SetScannerNRays()

virtual void genie::geometry::ROOTGeomAnalyzer::SetScannerNRays ( int nr )

inlinevirtual

Definition at line 82 of file ROOTGeomAnalyzer.h.

{ fNRays      = nr; } /* box  scanner */

References fNRays.

Referenced by Initialize(), and main().

SetTopVolName()

void ROOTGeomAnalyzer::SetTopVolName ( string nm )

virtual

Set the name of the top volume. This driver would ask the TGeoManager::GetTopVolume() for the top volume. Use this method for changing this if for example you want to set a smaller volume as the top one so as to generate events only in a specific part of your detector.

Definition at line 478 of file ROOTGeomAnalyzer.cxx.

{
/// Set the name of the top volume.
/// This driver would ask the TGeoManager::GetTopVolume() for the top volume.
/// Use this method for changing this if for example you want to set a smaller
/// volume as the top one so as to generate events only in a specific part of
/// your detector.
 
  if (name.size() == 0) return;
 
  fTopVolumeName = name;
  LOG("GROOTGeom",pNOTICE) << "Geometry Top Volume name: " << fTopVolumeName;
 
  TGeoVolume * gvol = fGeometry->GetVolume(fTopVolumeName.c_str());
  if (!gvol) {
     LOG("GROOTGeom",pWARN) << "Could not find volume: " << name.c_str();
     LOG("GROOTGeom",pWARN) << "Will not change the current top volume";
     fTopVolumeName = "";
     return;
  }
 
  // Get a matrix connecting coordinates of master and top volumes.
  // The matrix will be used for transforming the coordinates of incoming
  // flux neutrinos & generated interaction vertices.
  // This is needed (in case that the input top volume != master volume)
  // because ROOT always sets the coordinate system origin at the centre of
  // the specified top volume (whereas GENIE assumes that the global reference
  // frame is that of the master volume)
 
  TGeoIterator next(fGeometry->GetMasterVolume());
  TGeoNode *node;
  TString nodeName, volNameStr;
  const char* volName = fTopVolumeName.c_str();
  while ((node = next())) {
    nodeName = node->GetVolume()->GetName();
    if (nodeName == volName) {
      if (fMasterToTop) delete fMasterToTop;
      fMasterToTop = new TGeoHMatrix(*next.GetCurrentMatrix());
      fMasterToTopIsIdentity = fMasterToTop->IsIdentity();
      break;
    }
  }
 
  // set volume name
  fTopVolume = gvol;
  fGeometry->SetTopVolume(fTopVolume);
}

References fGeometry, fMasterToTop, fMasterToTopIsIdentity, fTopVolume, fTopVolumeName, LOG, pNOTICE, and pWARN.

Referenced by CreateRockBoxSelection(), and main().

SetWeightWithDensity()

virtual void genie::geometry::ROOTGeomAnalyzer::SetWeightWithDensity ( bool wt )

inlinevirtual

Definition at line 85 of file ROOTGeomAnalyzer.h.

{ fDensWeight = wt; }

References fDensWeight.

Referenced by Initialize().

SI2Local()

void ROOTGeomAnalyzer::SI2Local ( TVector3 & v ) const

virtual

convert position vector from SI units to current geometry units

Definition at line 573 of file ROOTGeomAnalyzer.cxx.

{
/// convert position vector from SI units to current geometry units
///
#ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
    LOG("GROOTGeom", pDEBUG)
      << "Position (SI): " << utils::print::Vec3AsString(&vec);
#endif
 
    vec *= (1./this->LengthUnits());
 
#ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
    LOG("GROOTGeom", pDEBUG)
      << "Position (loc): " << utils::print::Vec3AsString(&vec);
#endif
}

References LengthUnits(), LOG, pDEBUG, and genie::utils::print::Vec3AsString().

Referenced by ComputeMatLengths(), ComputePathLengths(), GenBoxRay(), and GenerateVertex().

Step()

double ROOTGeomAnalyzer::Step ( void )

protectedvirtual

Definition at line 1632 of file ROOTGeomAnalyzer.cxx.

{
  fGeometry->Step();
  double step=fGeometry->GetStep();
  return step;
}

References fGeometry.

Referenced by StepUntilEntering(), and SwimOnce().

StepToNextBoundary()

double ROOTGeomAnalyzer::StepToNextBoundary ( void )

protectedvirtual

Definition at line 1625 of file ROOTGeomAnalyzer.cxx.

{
  fGeometry->FindNextBoundary();
  double step=fGeometry->GetStep();
  return step;
}

References fGeometry.

Referenced by StepUntilEntering(), and SwimOnce().

StepUntilEntering()

double ROOTGeomAnalyzer::StepUntilEntering ( void )

protectedvirtual

Definition at line 1639 of file ROOTGeomAnalyzer.cxx.

{
  this->StepToNextBoundary();  // doesn't actually step, so don't include in sum
  double step = 0; //
 
  while(!fGeometry->IsEntering()) {
    step += this->Step();
  }
 
#ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
 
  bool isen = fGeometry->IsEntering();
  bool isob = fGeometry->IsOnBoundary();
 
  LOG("GROOTGeom",pDEBUG)
      << "IsEntering = "     << utils::print::BoolAsYNString(isen)
      << ", IsOnBoundary = " << utils::print::BoolAsYNString(isob)
      << ", Step = " << step;
#endif
 
  return step;
}

References genie::utils::print::BoolAsYNString(), fGeometry, LOG, pDEBUG, Step(), and StepToNextBoundary().

Referenced by SwimOnce().

SwimOnce()

void ROOTGeomAnalyzer::SwimOnce ( const TVector3 & r, const TVector3 & udir )

protectedvirtual

Swim through the geometry from the from the input position r0 (top vol coord & units) and moving along the direction of the unit vector udir (topvol coord) to create a filled PathSegmentList

Definition at line 1393 of file ROOTGeomAnalyzer.cxx.

{
/// Swim through the geometry from the from the input position
/// r0 (top vol coord & units) and moving along the direction of the
/// unit vector udir (topvol coord) to create a filled PathSegmentList
 
  int nvolswim = 0; //rwh
 
  if ( ! fCurrPathSegmentList ) fCurrPathSegmentList = new PathSegmentList();
 
  // don't swim if the current PathSegmentList is up-to-date
  if ( fCurrPathSegmentList->IsSameStart(r0,udir) ) return;
 
  // start fresh
  fCurrPathSegmentList->SetAllToZero();
 
  // set start info so next time we don't swim for the same ray
  fCurrPathSegmentList->SetStartInfo(r0,udir);
 
  PathSegment ps_curr;
 
  bool found_vol (false);
  bool keep_on   (true);
 
  double step    = 0;
  double raydist = 0;
 
  const TGeoVolume *   vol = 0;
  const TGeoMedium *   med = 0;
  const TGeoMaterial * mat = 0;
 
  // determining the geometry path is expensive, do it only if necessary
  bool selneedspath = ( fGeomVolSelector && fGeomVolSelector->GetNeedPath() );
  const bool fill_path = fKeepSegPath || selneedspath;
 
#ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
  LOG("GROOTGeom", pNOTICE)
    << "SwimOnce x [" << r0[0] << "," << r0[1] << "," << r0[2]
    << "] udir [" << udir[0] << "," << udir[1] << "," << udir[2];
#endif
 
  fGeometry -> SetCurrentDirection (udir[0],udir[1],udir[2]);
  fGeometry -> SetCurrentPoint     (r0[0],  r0[1],  r0[2]  );
 
  while (!found_vol || keep_on) {
     keep_on = true;
 
     fGeometry->FindNode();
 
     ps_curr.SetEnter( fGeometry->GetCurrentPoint() , raydist );
     vol = fGeometry->GetCurrentVolume();
     med = vol->GetMedium();
     mat = med->GetMaterial();
     ps_curr.SetGeo(vol,med,mat);
#ifdef PATHSEG_KEEP_PATH
     if (fill_path) ps_curr.SetPath(fGeometry->GetPath());
#endif
 
#ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
#ifdef DUMP_SWIM
       LOG("GROOTGeom", pDEBUG) << "Current volume: " << vol->GetName()
                             << " pos " << fGeometry->GetCurrentPoint()[0]
                             << " "     << fGeometry->GetCurrentPoint()[1]
                             << " "     << fGeometry->GetCurrentPoint()[2]
                             << " dir " << fGeometry->GetCurrentDirection()[0]
                             << " "     << fGeometry->GetCurrentDirection()[1]
                             << " "     << fGeometry->GetCurrentDirection()[2]
                             << "[path: " << fGeometry->GetPath() << "]";
#endif
#endif
 
     // find the start of top
     if (fGeometry->IsOutside() || !vol) {
        keep_on = false;
        if (found_vol) break;
        step = 0;
          this->StepToNextBoundary();
        //rwh//raydist += step;  // STNB doesn't actually "step"
 
#ifdef RWH_DEBUG
#ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
        LOG("GROOTGeom", pDEBUG) << "Outside ToNextBoundary step: " << step
                                 << " raydist: " << raydist;
#endif
#endif
 
        while (!fGeometry->IsEntering()) {
          step = this->Step();
          raydist += step;
#ifdef RWH_DEBUG
#ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
          LOG("GROOTGeom", pDEBUG)
              << "Stepping... [step size = " << step << "]";
          LOG("GROOTGeom", pDEBUG) << "Outside step: " << step
                                   << " raydist: " << raydist;
#endif
#endif
          if (this->WillNeverEnter(step)) {
#ifdef RWH_COUNTVOLS
            if ( accum_vol_stat ) {
              // this really shouldn't happen for the box exploration...
              // if coord transforms done right
              // could happen for neutrinos on a flux window
              nnever[curface]++;  //rwh
              if ( nnever[curface]%21 == 0 )
                LOG("GROOTGeom", pNOTICE)
                  << "curface " << curface << " " << nswims[curface]
                  << " never " << nnever[curface]
                  << " x [" << r0[0] << "," << r0[1] << "," << r0[2] << "] "
                  << " p [" << udir[0] << "," << udir[1] << "," << udir[2] << "]";
            }
#endif
            fCurrPathSegmentList->SetAllToZero();
            return;
          }
        } // finished while
 
        ps_curr.SetExit(fGeometry->GetCurrentPoint());
        ps_curr.SetStep(step);
        if ( ( fDebugFlags & 0x10 ) ) {
          // In general don't add the path segments from the start point to
          // the top volume (here for debug purposes)
          // Clear out the step range even if we keep it
          ps_curr.fStepRangeSet.clear();
          LOG("GROOTGeom", pNOTICE)
            << "debug: step towards top volume: " << ps_curr;
          fCurrPathSegmentList->AddSegment(ps_curr);
        }
 
     }  // outside or !vol
 
     if (keep_on) {
       if (!found_vol) found_vol = true;
 
       step   = this->StepUntilEntering();
       raydist += step;
 
       ps_curr.SetExit(fGeometry->GetCurrentPoint());
       ps_curr.SetStep(step);
       fCurrPathSegmentList->AddSegment(ps_curr);
 
       nvolswim++; //rwh
 
#ifdef DUMP_SWIM
       LOG("GROOTGeom", pDEBUG) << "Current volume: " << vol->GetName()
                                << " step " << step << " in " << mat->GetName();
#endif
 
#ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
       LOG("GROOTGeom", pDEBUG)
         << "Cur med.: " << med->GetName() << ", mat.: " << mat->GetName();
       LOG("GROOTGeom", pDEBUG)
         << "Step = " << step; // << ", weight = " << weight;
#endif
     } //keep on
  }
 
#ifdef RWH_COUNTVOLS
  if ( accum_vol_stat ) {
    nswims[curface]++;   //rwh
    dnvols[curface]  += (double)nvolswim;
    dnvols2[curface] += (double)nvolswim * (double)nvolswim;
    long int ns = fCurrPathSegmentList->size();
    if ( ns > mxsegments ) mxsegments = ns;
  }
#endif
 
//rwh:debug
#ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
  LOG("GROOTGeom", pDEBUG)
    << "PathSegmentList size " << fCurrPathSegmentList->size();
#endif
 
#ifdef RWH_DEBUG_2
  if ( ( fDebugFlags & 0x20 ) ) {
    fCurrPathSegmentList->SetDoCrossCheck(true);       //RWH
    LOG("GROOTGeom", pNOTICE) << "Before trimming" << *fCurrPathSegmentList;
    double mxddist = 0, mxdstep = 0;
    fCurrPathSegmentList->CrossCheck(mxddist,mxdstep);
    fmxddist = TMath::Max(fmxddist,mxddist);
    fmxdstep = TMath::Max(fmxdstep,mxdstep);
  }
#endif
 
  // PathSegmentList trimming occurs here!
  if ( fGeomVolSelector ) {
    PathSegmentList* altlist =
      fGeomVolSelector->GenerateTrimmedList(fCurrPathSegmentList);
    std::swap(altlist,fCurrPathSegmentList);
    delete altlist;  // after swap delete original
  }
 
  fCurrPathSegmentList->FillMatStepSum();
 
#ifdef RWH_DEBUG_2
  if ( fGeomVolSelector) {
    // after FillMatStepSum() so one can see the summed mass
    if ( ( fDebugFlags & 0x40 ) ) {
      fCurrPathSegmentList->SetPrintVerbose(true);
      LOG("GROOTGeom", pNOTICE) << "After  trimming" << *fCurrPathSegmentList;
      fCurrPathSegmentList->SetPrintVerbose(false);
    }
  }
#endif
 
 
  return;
}

References fCurrPathSegmentList, fDebugFlags, fGeometry, fGeomVolSelector, fKeepSegPath, fmxddist, fmxdstep, genie::geometry::PathSegment::fStepRangeSet, LOG, pDEBUG, pNOTICE, r0, genie::geometry::PathSegment::SetEnter(), genie::geometry::PathSegment::SetExit(), genie::geometry::PathSegment::SetGeo(), genie::geometry::PathSegment::SetPath(), genie::geometry::PathSegment::SetStep(), Step(), StepToNextBoundary(), StepUntilEntering(), and WillNeverEnter().

Referenced by ComputeMatLengths(), and ComputePathLengthPDG().

Top2Master()

void ROOTGeomAnalyzer::Top2Master ( TVector3 & v ) const

virtual

transform the input position vector from the specified top volume coordinate system to the master volume coordinate system, but not change the units.

Definition at line 637 of file ROOTGeomAnalyzer.cxx.

{
/// transform the input position vector from the specified top volume
/// coordinate system to the master volume coordinate system, but not
/// change the units.
 
#ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
    LOG("GROOTGeom", pDEBUG)
      << "Position (coord:top): " << utils::print::Vec3AsString(&vec);
#endif
 
    Double_t mast[3], top[3];
    vec.GetXYZ(top);
    fMasterToTop->LocalToMaster(top, mast);
    vec.SetXYZ(mast[0], mast[1], mast[2]);
 
#ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
    LOG("GROOTGeom", pDEBUG)
      << "Position (coord:master): " << utils::print::Vec3AsString(&vec);
#endif
}

References fMasterToTop, LOG, pDEBUG, and genie::utils::print::Vec3AsString().

Referenced by GenBoxRay(), and GenerateVertex().

Top2MasterDir()

void ROOTGeomAnalyzer::Top2MasterDir ( TVector3 & v ) const

virtual

transform the input direction vector from the specified top volume coordinate system to the master volume coordinate system.

Definition at line 660 of file ROOTGeomAnalyzer.cxx.

{
/// transform the input direction vector from the specified top volume
/// coordinate system to the master volume coordinate system.
 
#ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
    LOG("GROOTGeom", pDEBUG)
      << "Direction (coord:top): " << utils::print::Vec3AsString(&vec);
#endif
 
    Double_t mast[3], top[3];
    vec.GetXYZ(top);
    fMasterToTop->LocalToMasterVect(top, mast);
    vec.SetXYZ(mast[0], mast[1], mast[2]);
 
#ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
    LOG("GROOTGeom", pDEBUG)
      << "Direction (coord:master): " << utils::print::Vec3AsString(&vec);
#endif
}

References fMasterToTop, LOG, pDEBUG, and genie::utils::print::Vec3AsString().

Referenced by GenBoxRay().

TopVolName()

virtual string genie::geometry::ROOTGeomAnalyzer::TopVolName ( void ) const

inlinevirtual

Definition at line 104 of file ROOTGeomAnalyzer.h.

{ return fTopVolumeName;     }

References fTopVolumeName.

WeightWithDensity()

virtual bool genie::geometry::ROOTGeomAnalyzer::WeightWithDensity ( void ) const

inlinevirtual

Definition at line 99 of file ROOTGeomAnalyzer.h.

{ return fDensWeight;        }

References fDensWeight.

Referenced by GetWeight(), GetWeight(), and Local2SI().

WillNeverEnter()

bool ROOTGeomAnalyzer::WillNeverEnter ( double step )

protectedvirtual

Definition at line 1662 of file ROOTGeomAnalyzer.cxx.

{
// If the neutrino trajectory would never enter the detector, then the
// TGeoManager::GetStep returns the maximum step (1E30).
// Compare surrent step with max step and figure out whether the particle
// would never enter the detector
 
  if (step > 9.99E29) {
 
#ifdef __GENIE_LOW_LEVEL_MESG_ENABLED__
     LOG("GROOTGeom", pINFO) << "Wow! Current step is dr = " << step;
     LOG("GROOTGeom", pINFO) << "This trajectory isn't entering the detector";
#endif
     return true;
 
  } else
    return false;
}

References LOG, and pINFO.

Referenced by SwimOnce().

Member Data Documentation

fCurrMaxPathLengthList

PathLengthList* genie::geometry::ROOTGeomAnalyzer::fCurrMaxPathLengthList

protected

current list of max path-lengths

Definition at line 166 of file ROOTGeomAnalyzer.h.

Referenced by CleanUp(), ComputeMaxPathLengths(), GetMaxPathLengths(), Initialize(), Load(), MaxPathLengthsBoxMethod(), and MaxPathLengthsFluxMethod().

fCurrPathLengthList

PathLengthList* genie::geometry::ROOTGeomAnalyzer::fCurrPathLengthList

protected

current list of path-lengths

Definition at line 165 of file ROOTGeomAnalyzer.h.

Referenced by CleanUp(), ComputePathLengths(), Initialize(), and Load().

fCurrPathSegmentList

PathSegmentList* genie::geometry::ROOTGeomAnalyzer::fCurrPathSegmentList

protected

current list of path-segments

Definition at line 173 of file ROOTGeomAnalyzer.h.

Referenced by CleanUp(), ComputeMatLengths(), ComputePathLengthPDG(), GenerateVertex(), Initialize(), Load(), and SwimOnce().

fCurrPDGCodeList

PDGCodeList* genie::geometry::ROOTGeomAnalyzer::fCurrPDGCodeList

protected

current list of target nuclei

Definition at line 167 of file ROOTGeomAnalyzer.h.

Referenced by BuildListOfTargetNuclei(), CleanUp(), ComputePathLengths(), GetWeight(), Initialize(), and ListOfTargetNuclei().

fCurrVertex

TVector3* genie::geometry::ROOTGeomAnalyzer::fCurrVertex

protected

current generated vertex

Definition at line 164 of file ROOTGeomAnalyzer.h.

Referenced by GenerateVertex(), and Load().

fDebugFlags

int genie::geometry::ROOTGeomAnalyzer::fDebugFlags

protected

Definition at line 185 of file ROOTGeomAnalyzer.h.

Referenced by GenerateVertex(), Initialize(), SetDebugFlags(), and SwimOnce().

fDensityScale

double genie::geometry::ROOTGeomAnalyzer::fDensityScale

protected

conversion factor: input geometry density units -> kgr/meters^3

Definition at line 161 of file ROOTGeomAnalyzer.h.

Referenced by DensityUnits(), and SetDensityUnits().

fDensWeight

bool genie::geometry::ROOTGeomAnalyzer::fDensWeight

protected

if true pathlengths are weighted with density [def:true]

Definition at line 159 of file ROOTGeomAnalyzer.h.

Referenced by MaxPathLengthsBoxMethod(), SetWeightWithDensity(), and WeightWithDensity().

fdx

double genie::geometry::ROOTGeomAnalyzer::fdx

protected

Definition at line 181 of file ROOTGeomAnalyzer.h.

Referenced by GenBoxRay().

fdy

double genie::geometry::ROOTGeomAnalyzer::fdy

protected

Definition at line 181 of file ROOTGeomAnalyzer.h.

Referenced by GenBoxRay().

fdz

double genie::geometry::ROOTGeomAnalyzer::fdz

protected

Definition at line 181 of file ROOTGeomAnalyzer.h.

Referenced by GenBoxRay().

fFlux

GFluxI* genie::geometry::ROOTGeomAnalyzer::fFlux

protected

a flux objects that can be used to scan the max path lengths

Definition at line 158 of file ROOTGeomAnalyzer.h.

Referenced by ComputeMaxPathLengths(), MaxPathLengthsFluxMethod(), and SetScannerFlux().

fGenBoxRayDir

TVector3 genie::geometry::ROOTGeomAnalyzer::fGenBoxRayDir

protected

Definition at line 178 of file ROOTGeomAnalyzer.h.

Referenced by GenBoxRay().

fGenBoxRayPos

TVector3 genie::geometry::ROOTGeomAnalyzer::fGenBoxRayPos

protected

Definition at line 177 of file ROOTGeomAnalyzer.h.

Referenced by GenBoxRay().

fGeometry

TGeoManager* genie::geometry::ROOTGeomAnalyzer::fGeometry

protected

input detector geometry

Definition at line 153 of file ROOTGeomAnalyzer.h.

Referenced by BuildListOfTargetNuclei(), ComputeMaxPathLengths(), FindMaterialInCurrentVol(), GenerateVertex(), GetGeometry(), Load(), SetTopVolName(), Step(), StepToNextBoundary(), StepUntilEntering(), and SwimOnce().

fGeomVolSelector

GeomVolSelectorI* genie::geometry::ROOTGeomAnalyzer::fGeomVolSelector

protected

optional path seg trimmer (owned)

Definition at line 174 of file ROOTGeomAnalyzer.h.

Referenced by AdoptGeomVolSelector(), ComputeMatLengths(), ComputePathLengths(), Initialize(), and SwimOnce().

fiface

int genie::geometry::ROOTGeomAnalyzer::fiface

protected

Definition at line 179 of file ROOTGeomAnalyzer.h.

Referenced by GenBoxRay().

fipoint

int genie::geometry::ROOTGeomAnalyzer::fipoint

protected

Definition at line 179 of file ROOTGeomAnalyzer.h.

Referenced by GenBoxRay().

firay

int genie::geometry::ROOTGeomAnalyzer::firay

protected

Definition at line 179 of file ROOTGeomAnalyzer.h.

Referenced by GenBoxRay().

fKeepSegPath

bool genie::geometry::ROOTGeomAnalyzer::fKeepSegPath

protected

need to fill path segment "path"

Definition at line 172 of file ROOTGeomAnalyzer.h.

Referenced by GetKeepSegPath(), Initialize(), SetKeepSegPath(), and SwimOnce().

fLengthScale

double genie::geometry::ROOTGeomAnalyzer::fLengthScale

protected

conversion factor: input geometry length units -> meters

Definition at line 160 of file ROOTGeomAnalyzer.h.

Referenced by LengthUnits(), and SetLengthUnits().

fMasterToTop

TGeoHMatrix* genie::geometry::ROOTGeomAnalyzer::fMasterToTop

protected

matrix connecting master coordinates to top volume coordinates

Definition at line 169 of file ROOTGeomAnalyzer.h.

Referenced by CleanUp(), Load(), Master2Top(), Master2TopDir(), SetTopVolName(), Top2Master(), and Top2MasterDir().

fMasterToTopIsIdentity

bool genie::geometry::ROOTGeomAnalyzer::fMasterToTopIsIdentity

protected

is fMasterToTop matrix the identity matrix?

Definition at line 170 of file ROOTGeomAnalyzer.h.

Referenced by ComputeMatLengths(), ComputePathLengths(), GenBoxRay(), GenerateVertex(), Initialize(), Load(), and SetTopVolName().

fMaterial

int genie::geometry::ROOTGeomAnalyzer::fMaterial

protected

input selected material for vertex generation

Definition at line 152 of file ROOTGeomAnalyzer.h.

fMaxPlSafetyFactor

double genie::geometry::ROOTGeomAnalyzer::fMaxPlSafetyFactor

protected

factor that can multiply the computed max path lengths

Definition at line 162 of file ROOTGeomAnalyzer.h.

Referenced by MaxPlSafetyFactor(), and SetMaxPlSafetyFactor().

fMixtWghtSum

double genie::geometry::ROOTGeomAnalyzer::fMixtWghtSum

protected

norm of relative weights (<0 if explicit summing required)

Definition at line 163 of file ROOTGeomAnalyzer.h.

Referenced by MaxPathLengthsBoxMethod(), MixtureWeightsSum(), and SetMixtureWeightsSum().

fmxddist

double genie::geometry::ROOTGeomAnalyzer::fmxddist

protected

Definition at line 184 of file ROOTGeomAnalyzer.h.

Referenced by GenerateVertex(), Initialize(), SwimOnce(), and ~ROOTGeomAnalyzer().

fmxdstep

double genie::geometry::ROOTGeomAnalyzer::fmxdstep

protected

max errors in pathsegmentlist

Definition at line 184 of file ROOTGeomAnalyzer.h.

Referenced by GenerateVertex(), Initialize(), SwimOnce(), and ~ROOTGeomAnalyzer().

fnewpnt

bool genie::geometry::ROOTGeomAnalyzer::fnewpnt

protected

Definition at line 180 of file ROOTGeomAnalyzer.h.

Referenced by GenBoxRay().

fNParticles

int genie::geometry::ROOTGeomAnalyzer::fNParticles

protected

max path length scanner (flux method): particles in [def:10000]

Definition at line 157 of file ROOTGeomAnalyzer.h.

Referenced by ScannerNParticles(), and SetScannerNParticles().

fNPoints

int genie::geometry::ROOTGeomAnalyzer::fNPoints

protected

max path length scanner (box method): points/surface [def:200]

Definition at line 155 of file ROOTGeomAnalyzer.h.

Referenced by GenBoxRay(), ScannerNPoints(), and SetScannerNPoints().

fNRays

int genie::geometry::ROOTGeomAnalyzer::fNRays

protected

max path length scanner (box method): rays/point [def:200]

Definition at line 156 of file ROOTGeomAnalyzer.h.

Referenced by GenBoxRay(), ScannerNRays(), and SetScannerNRays().

fox

double genie::geometry::ROOTGeomAnalyzer::fox

protected

Definition at line 181 of file ROOTGeomAnalyzer.h.

Referenced by GenBoxRay().

foy

double genie::geometry::ROOTGeomAnalyzer::foy

protected

Definition at line 181 of file ROOTGeomAnalyzer.h.

Referenced by GenBoxRay().

foz

double genie::geometry::ROOTGeomAnalyzer::foz

protected

top vol size/origin (top vol units)

Definition at line 181 of file ROOTGeomAnalyzer.h.

Referenced by GenBoxRay().

fTopVolume

TGeoVolume* genie::geometry::ROOTGeomAnalyzer::fTopVolume

protected

top volume

Definition at line 168 of file ROOTGeomAnalyzer.h.

Referenced by GenBoxRay(), Initialize(), Load(), and SetTopVolName().

fTopVolumeName

string genie::geometry::ROOTGeomAnalyzer::fTopVolumeName

protected

input top vol [other than TGeoManager::GetTopVolume()]

Definition at line 154 of file ROOTGeomAnalyzer.h.

Referenced by Initialize(), SetTopVolName(), and TopVolName().

---

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

• ROOTGeomAnalyzer.h
• ROOTGeomAnalyzer.cxx