genie::geometry::ROOTGeomAnalyzer Class Reference

A ROOT/GEANT4 geometry driver. More...

#include "/cvmfs/nova.opensciencegrid.org/externals/genie/v3_00_06_p01/Linux64bit+3.10-2.17-e17-debug/GENIE-Generator/src/Tools/Geometry/ROOTGeomAnalyzer.h"

Inheritance 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 More...
virtual const PathLengthList &	ComputeMaxPathLengths (void)
virtual const PathLengthList &	ComputePathLengths (const TLorentzVector &x, const TLorentzVector &p)
virtual const TVector3 &	GenerateVertex (const TLorentzVector &x, const TLorentzVector &p, int tgtpdg)
virtual void	SetScannerNPoints (int np)
	set geometry driver's configuration options More...
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 More...
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 More...
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 More...

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

Detailed Description

A ROOT/GEANT4 geometry driver.

Author

Anselmo Meregaglia <anselmo.meregaglia cern.ch>, ETH Zurich Costas Andreopoulos <costas.andreopoulos stfc.ac.uk>, STFC, Rutherford Lab Robert Hatcher <rhatcher fnal.gov>, Fermilab

May 24, 2005

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

Definition at line 49 of file ROOTGeomAnalyzer.h.

Constructor & Destructor Documentation

ROOTGeomAnalyzer::ROOTGeomAnalyzer

(

string

geometry_filename

)

Constructor from a geometry file

Definition at line 143 of file ROOTGeomAnalyzer.cxx.

References Initialize(), Load(), LOG, and pDEBUG.

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

ROOTGeomAnalyzer::ROOTGeomAnalyzer

(

TGeoManager *

gm

)

Constructor from a TGeoManager

Definition at line 156 of file ROOTGeomAnalyzer.cxx.

References Initialize(), Load(), LOG, and pDEBUG.

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

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

inline

Definition at line 54 of file ROOTGeomAnalyzer.h.

References ComputeMaxPathLengths(), ComputePathLengths(), GenerateVertex(), ListOfTargetNuclei(), submit_syst::x, and ~ROOTGeomAnalyzer().

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

ROOTGeomAnalyzer::~ROOTGeomAnalyzer

(

)

Definition at line 169 of file ROOTGeomAnalyzer.cxx.

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

Referenced by ROOTGeomAnalyzer().

{
  this->CleanUp();

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

Member Function Documentation

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

inlinevirtual

configure processing to perform path segment trimming

take ownership, return old

Definition at line 109 of file ROOTGeomAnalyzer.h.

References BuildListOfTargetNuclei(), CleanUp(), ComputePathLengthPDG(), fGeomVolSelector, FindMaterialInCurrentVol(), GenBoxRay(), GetTargetPdgCode(), GetWeight(), Initialize(), Load(), m, MaxPathLengthsBoxMethod(), MaxPathLengthsFluxMethod(), make_associated_cosmic_defs::p4, r(), fillBadChanDBTables::step, Step(), StepToNextBoundary(), StepUntilEntering(), std::swap(), SwimOnce(), and WillNeverEnter().

Referenced by CreateFidSelection().

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

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 827 of file ROOTGeomAnalyzer.cxx.

References exit(), fCurrPDGCodeList, fGeometry, GetTargetPdgCode(), MECModelEnuComparisons::i, LOG, pDEBUG, pERROR, pFATAL, genie::PDGCodeList::push_back(), and pWARN.

Referenced by AdoptGeomVolSelector(), and Load().

{
/// 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());

}

void ROOTGeomAnalyzer::CleanUp ( void  )

protectedvirtual

Definition at line 754 of file ROOTGeomAnalyzer.cxx.

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

Referenced by AdoptGeomVolSelector(), and ~ROOTGeomAnalyzer().

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

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 190 of file ROOTGeomAnalyzer.cxx.

References genie::GFluxI::Clear(), exit(), fCurrMaxPathLengthList, fFlux, fGeometry, LOG, MaxPathLengthsBoxMethod(), MaxPathLengthsFluxMethod(), pFATAL, pNOTICE, and genie::PathLengthList::SetAllToZero().

Referenced by main(), and ROOTGeomAnalyzer().

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

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 1362 of file ROOTGeomAnalyzer.cxx.

References fCurrPathSegmentList, genie::geometry::PathSegmentList::GetMatStepSumMap(), GetWeight(), LOG, pDEBUG, std_candles::pl, fillBadChanDBTables::step, SwimOnce(), and ana::weight.

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

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

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 221 of file ROOTGeomAnalyzer.cxx.

References genie::PathLengthList::AddPathLength(), ComputePathLengthPDG(), fCurrPathLengthList, fCurrPDGCodeList, fGeomVolSelector, fMasterToTopIsIdentity, LengthUnits(), Local2SI(), LOG, Master2Top(), Master2TopDir(), genie::utils::print::P4AsShortString(), pDEBUG, pINFO, std_candles::pl, elec2geo::pos, genie::PathLengthList::SetAllToZero(), genie::geometry::GeomVolSelectorI::SetCurrentRay(), genie::geometry::GeomVolSelectorI::SetSI2Local(), SI2Local(), and genie::utils::print::X4AsString().

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

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

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

inlinevirtual

Definition at line 89 of file ROOTGeomAnalyzer.h.

References fDensityScale.

Referenced by Local2SI().

{ return fDensityScale;      }

bool ROOTGeomAnalyzer::FindMaterialInCurrentVol ( int  pdgc )

protectedvirtual

Definition at line 1612 of file ROOTGeomAnalyzer.cxx.

References fGeometry, GetTargetPdgCode(), MECModelEnuComparisons::i, LOG, make_root_from_grid_output::pdg, and pWARN.

Referenced by AdoptGeomVolSelector(), and GenerateVertex().

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

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 1182 of file ROOTGeomAnalyzer.cxx.

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, elec2geo::pos, generate_hists::rnd, genie::RandomGen::RndGeom(), SI2Local(), Top2Master(), Top2MasterDir(), and genie::utils::print::Vec3AsString().

Referenced by AdoptGeomVolSelector(), and MaxPathLengthsBoxMethod().

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

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 279 of file ROOTGeomAnalyzer.cxx.

References ComputePathLengthPDG(), genie::geometry::PathSegmentList::CrossCheck(), exit(), fCurrPathSegmentList, fCurrVertex, fDebugFlags, fGeometry, FindMaterialInCurrentVol(), genie::utils::xml::FindNode(), fMasterToTopIsIdentity, genie::geometry::PathSegment::fMaterial, fmxddist, fmxdstep, util::frac(), genie::geometry::PathSegment::fVolume, genie::geometry::PathSegmentList::GetMatStepSumMap(), genie::geometry::PathSegmentList::GetPathSegmentV(), 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, elec2geo::pos, pWARN, generate_hists::rnd, genie::RandomGen::RndGeom(), genie::geometry::PathSegmentList::SetDoCrossCheck(), SI2Local(), Top2Master(), genie::utils::print::Vec3AsString(), wgt, and genie::utils::print::X4AsString().

Referenced by main(), and ROOTGeomAnalyzer().

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

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

inlinevirtual

Definition at line 93 of file ROOTGeomAnalyzer.h.

References fGeometry.

Referenced by GetGeometry(), and main().

{ return fGeometry;          }

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

inlinevirtual

Definition at line 94 of file ROOTGeomAnalyzer.h.

References fKeepSegPath.

{ return fKeepSegPath;       }

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

inlinevirtual

Definition at line 95 of file ROOTGeomAnalyzer.h.

References fCurrMaxPathLengthList, Local2SI(), Master2Top(), Master2TopDir(), std_candles::pl, SI2Local(), Top2Master(), Top2MasterDir(), and registry_explorer::v.

Referenced by main().

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

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

protectedvirtual

Definition at line 892 of file ROOTGeomAnalyzer.cxx.

References genie::units::A, genie::pdg::IonPdgCode(), and Z.

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

{
  int A = TMath::Nint(m->GetA());
  int Z = TMath::Nint(m->GetZ());

  int pdgc = pdg::IonPdgCode(A,Z);

  return pdgc;
}

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

protectedvirtual

Definition at line 903 of file ROOTGeomAnalyzer.cxx.

References genie::units::A, genie::pdg::IonPdgCode(), and Z.

{
  int A = TMath::Nint(m->GetAmixt()[ielement]);
  int Z = TMath::Nint(m->GetZmixt()[ielement]);

  int pdgc = pdg::IonPdgCode(A,Z);

  return pdgc;
}

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 915 of file ROOTGeomAnalyzer.cxx.

References genie::PDGCodeList::ExistsInPDGCodeList(), fCurrPDGCodeList, GetTargetPdgCode(), LOG, pDEBUG, pERROR, ana::weight, and WeightWithDensity().

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

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

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 976 of file ROOTGeomAnalyzer.cxx.

References exit(), GetTargetPdgCode(), GetWeight(), MECModelEnuComparisons::i, calib::j, LOG, Munits::nm, pERROR, pFATAL, pWARN, ana::weight, and WeightWithDensity().

{
/// 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();
     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;
}

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 1015 of file ROOTGeomAnalyzer.cxx.

References ana::assert(), d, GetTargetPdgCode(), MECModelEnuComparisons::i, MixtureWeightsSum(), w, and ana::weight.

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

void ROOTGeomAnalyzer::Initialize ( void  )

protectedvirtual

Definition at line 721 of file ROOTGeomAnalyzer.cxx.

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 AdoptGeomVolSelector(), and ROOTGeomAnalyzer().

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

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

inlinevirtual

Definition at line 88 of file ROOTGeomAnalyzer.h.

References fLengthScale.

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

{ return fLengthScale;       }

const PDGCodeList & ROOTGeomAnalyzer::ListOfTargetNuclei ( void  )

virtual

implement the GeomAnalyzerI interface

Implements genie::GeomAnalyzerI.

Definition at line 184 of file ROOTGeomAnalyzer.cxx.

References fCurrPDGCodeList.

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

{
  return *fCurrPDGCodeList;
}

void ROOTGeomAnalyzer::Load ( string  geometry_filename )

protectedvirtual

Load the detector geometry from the input ROOT file

Definition at line 766 of file ROOTGeomAnalyzer.cxx.

References exit(), shutoffs::filename, gm, LOG, pFATAL, and pNOTICE.

Referenced by AdoptGeomVolSelector(), and ROOTGeomAnalyzer().

{
/// 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);
  }
  TGeoManager * gm = TGeoManager::Import(filename.c_str());

  this->Load(gm);
}

void ROOTGeomAnalyzer::Load ( TGeoManager *  gm )

protectedvirtual

Load the detector geometry from the input TGeoManager

Definition at line 784 of file ROOTGeomAnalyzer.cxx.

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

{
/// 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

}

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 566 of file ROOTGeomAnalyzer.cxx.

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

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

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

void ROOTGeomAnalyzer::Local2SI ( TVector3 &  v ) const

virtual

convert position vector from current geometry units to SI units

Definition at line 590 of file ROOTGeomAnalyzer.cxx.

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

{
/// 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
}

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 627 of file ROOTGeomAnalyzer.cxx.

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

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

{
/// 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
}

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 650 of file ROOTGeomAnalyzer.cxx.

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

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

{
/// 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
}

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 1110 of file ROOTGeomAnalyzer.cxx.

References ComputePathLengths(), fCurrMaxPathLengthList, fDensWeight, fMixtWghtSum, GenBoxRay(), calib::j, LOG, MaxPlSafetyFactor(), ns, genie::PathLengthList::PathLength(), pDEBUG, std_candles::pl, pNOTICE, extractScale::rms, genie::PathLengthList::SetPathLength(), ana::Sqrt(), submit_syst::x, and submit_syst::x2.

Referenced by AdoptGeomVolSelector(), and ComputeMaxPathLengths().

{
/// 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

}

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 1046 of file ROOTGeomAnalyzer.cxx.

References abs(), ComputePathLengths(), emax, fCurrMaxPathLengthList, fFlux, genie::GFluxI::GenerateNext(), LOG, genie::GFluxI::MaxEnergy(), MaxPlSafetyFactor(), genie::GFluxI::Momentum(), genie::PathLengthList::PathLength(), std_candles::pl, pNOTICE, genie::GFluxI::Position(), pWARN, ScannerNParticles(), and genie::PathLengthList::SetPathLength().

Referenced by AdoptGeomVolSelector(), and ComputeMaxPathLengths().

{
/// 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++;  
  }
}

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

inlinevirtual

Definition at line 91 of file ROOTGeomAnalyzer.h.

References fMaxPlSafetyFactor.

Referenced by MaxPathLengthsBoxMethod(), and MaxPathLengthsFluxMethod().

{ return fMaxPlSafetyFactor; }

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

inlinevirtual

Definition at line 90 of file ROOTGeomAnalyzer.h.

References fMixtWghtSum.

Referenced by GetWeight().

{ return fMixtWghtSum;       }

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

inlinevirtual

Definition at line 86 of file ROOTGeomAnalyzer.h.

References fNParticles.

Referenced by MaxPathLengthsFluxMethod().

{ return fNParticles;        }

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

inlinevirtual

retrieve geometry driver's configuration options

Definition at line 84 of file ROOTGeomAnalyzer.h.

References fNPoints.

{ return fNPoints;           }

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

inlinevirtual

Definition at line 85 of file ROOTGeomAnalyzer.h.

References fNRays.

{ return fNRays;             }

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

inlinevirtual

Definition at line 80 of file ROOTGeomAnalyzer.h.

References fDebugFlags.

{ fDebugFlags  = flgs; }

void ROOTGeomAnalyzer::SetDensityUnits ( double  du )

virtual

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

Definition at line 475 of file ROOTGeomAnalyzer.cxx.

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

Referenced by Initialize(), and SetWeightWithDensity().

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

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

inlinevirtual

Definition at line 79 of file ROOTGeomAnalyzer.h.

References fKeepSegPath, and sanity_check_grl::keep.

{ fKeepSegPath = keep; }

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 458 of file ROOTGeomAnalyzer.cxx.

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

Referenced by Initialize(), and SetWeightWithDensity().

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

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 486 of file ROOTGeomAnalyzer.cxx.

References fMaxPlSafetyFactor, LOG, and pNOTICE.

Referenced by Initialize(), and SetWeightWithDensity().

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

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 504 of file ROOTGeomAnalyzer.cxx.

References fMixtWghtSum, and sum.

Referenced by Initialize(), and SetWeightWithDensity().

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

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

inlinevirtual

Definition at line 72 of file ROOTGeomAnalyzer.h.

References MakeMiniprodValidationCuts::f, and fFlux.

Referenced by Initialize(), and main().

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

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

inlinevirtual

Definition at line 71 of file ROOTGeomAnalyzer.h.

References fNParticles.

Referenced by Initialize(), and main().

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

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

inlinevirtual

set geometry driver's configuration options

Definition at line 69 of file ROOTGeomAnalyzer.h.

References fNPoints.

Referenced by Initialize(), and main().

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

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

inlinevirtual

Definition at line 70 of file ROOTGeomAnalyzer.h.

References fNRays.

Referenced by Initialize(), and main().

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

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 514 of file ROOTGeomAnalyzer.cxx.

References fGeometry, fMasterToTop, fMasterToTopIsIdentity, fTopVolume, fTopVolumeName, LOG, next(), file_size_ana::node, pNOTICE, and pWARN.

Referenced by main(), and SetWeightWithDensity().

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

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

inlinevirtual

Definition at line 73 of file ROOTGeomAnalyzer.h.

References fDensWeight, Munits::nm, SetDensityUnits(), SetLengthUnits(), SetMaxPlSafetyFactor(), SetMixtureWeightsSum(), SetTopVolName(), and sum.

Referenced by Initialize().

{ fDensWeight = wt; }

void ROOTGeomAnalyzer::SI2Local ( TVector3 &  v ) const

virtual

convert position vector from SI units to current geometry units

Definition at line 609 of file ROOTGeomAnalyzer.cxx.

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

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

{
/// 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
}

double ROOTGeomAnalyzer::Step ( void  )

protectedvirtual

Definition at line 1641 of file ROOTGeomAnalyzer.cxx.

References fGeometry, and fillBadChanDBTables::step.

Referenced by AdoptGeomVolSelector(), StepUntilEntering(), and SwimOnce().

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

double ROOTGeomAnalyzer::StepToNextBoundary ( void  )

protectedvirtual

Definition at line 1634 of file ROOTGeomAnalyzer.cxx.

References fGeometry, and fillBadChanDBTables::step.

Referenced by AdoptGeomVolSelector(), StepUntilEntering(), and SwimOnce().

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

double ROOTGeomAnalyzer::StepUntilEntering ( void  )

protectedvirtual

Definition at line 1648 of file ROOTGeomAnalyzer.cxx.

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

Referenced by AdoptGeomVolSelector(), and SwimOnce().

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

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 1402 of file ROOTGeomAnalyzer.cxx.

References genie::geometry::PathSegmentList::AddSegment(), genie::geometry::PathSegmentList::CrossCheck(), fCurrPathSegmentList, fDebugFlags, fGeometry, fGeomVolSelector, genie::geometry::PathSegmentList::FillMatStepSum(), fKeepSegPath, fmxddist, fmxdstep, genie::geometry::PathSegment::fStepRangeSet, genie::geometry::GeomVolSelectorI::GenerateTrimmedList(), genie::geometry::GeomVolSelectorI::GetNeedPath(), genie::geometry::PathSegmentList::IsSameStart(), LOG, ns, pDEBUG, pNOTICE, genie::geometry::PathSegmentList::SetAllToZero(), genie::geometry::PathSegmentList::SetDoCrossCheck(), genie::geometry::PathSegment::SetEnter(), genie::geometry::PathSegment::SetExit(), genie::geometry::PathSegment::SetGeo(), genie::geometry::PathSegment::SetPath(), genie::geometry::PathSegmentList::SetPrintVerbose(), genie::geometry::PathSegmentList::SetStartInfo(), genie::geometry::PathSegment::SetStep(), genie::geometry::PathSegmentList::size(), fillBadChanDBTables::step, Step(), StepToNextBoundary(), StepUntilEntering(), std::swap(), and WillNeverEnter().

Referenced by AdoptGeomVolSelector(), and ComputePathLengthPDG().

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

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 673 of file ROOTGeomAnalyzer.cxx.

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

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

{
/// 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
}

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 696 of file ROOTGeomAnalyzer.cxx.

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

Referenced by GenBoxRay(), and GetMaxPathLengths().

{
/// 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
}

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

inlinevirtual

Definition at line 92 of file ROOTGeomAnalyzer.h.

References fTopVolumeName.

{ return fTopVolumeName;     }

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

inlinevirtual

Definition at line 87 of file ROOTGeomAnalyzer.h.

References fDensWeight.

Referenced by GetWeight(), and Local2SI().

{ return fDensWeight;        }

bool ROOTGeomAnalyzer::WillNeverEnter ( double  step )

protectedvirtual

Definition at line 1671 of file ROOTGeomAnalyzer.cxx.

References LOG, pINFO, and fillBadChanDBTables::step.

Referenced by AdoptGeomVolSelector(), and SwimOnce().

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

Member Data Documentation

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

protected

current list of max path-lengths

Definition at line 156 of file ROOTGeomAnalyzer.h.

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

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

protected

current list of path-lengths

Definition at line 155 of file ROOTGeomAnalyzer.h.

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

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

protected

current list of path-segments

Definition at line 163 of file ROOTGeomAnalyzer.h.

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

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

protected

current list of target nuclei

Definition at line 157 of file ROOTGeomAnalyzer.h.

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

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

protected

current generated vertex

Definition at line 154 of file ROOTGeomAnalyzer.h.

Referenced by GenerateVertex(), and Load().

int genie::geometry::ROOTGeomAnalyzer::fDebugFlags

protected

Definition at line 175 of file ROOTGeomAnalyzer.h.

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

double genie::geometry::ROOTGeomAnalyzer::fDensityScale

protected

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

Definition at line 151 of file ROOTGeomAnalyzer.h.

Referenced by DensityUnits(), and SetDensityUnits().

bool genie::geometry::ROOTGeomAnalyzer::fDensWeight

protected

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

Definition at line 149 of file ROOTGeomAnalyzer.h.

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

double genie::geometry::ROOTGeomAnalyzer::fdx

protected

Definition at line 171 of file ROOTGeomAnalyzer.h.

Referenced by GenBoxRay().

double genie::geometry::ROOTGeomAnalyzer::fdy

protected

Definition at line 171 of file ROOTGeomAnalyzer.h.

Referenced by GenBoxRay().

double genie::geometry::ROOTGeomAnalyzer::fdz

protected

Definition at line 171 of file ROOTGeomAnalyzer.h.

Referenced by GenBoxRay().

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

protected

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

Definition at line 148 of file ROOTGeomAnalyzer.h.

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

TVector3 genie::geometry::ROOTGeomAnalyzer::fGenBoxRayDir

protected

Definition at line 168 of file ROOTGeomAnalyzer.h.

Referenced by GenBoxRay().

TVector3 genie::geometry::ROOTGeomAnalyzer::fGenBoxRayPos

protected

Definition at line 167 of file ROOTGeomAnalyzer.h.

Referenced by GenBoxRay().

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

protected

input detector geometry

Definition at line 143 of file ROOTGeomAnalyzer.h.

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

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

protected

optional path seg trimmer (owned)

Definition at line 164 of file ROOTGeomAnalyzer.h.

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

int genie::geometry::ROOTGeomAnalyzer::fiface

protected

Definition at line 169 of file ROOTGeomAnalyzer.h.

Referenced by GenBoxRay().

int genie::geometry::ROOTGeomAnalyzer::fipoint

protected

Definition at line 169 of file ROOTGeomAnalyzer.h.

Referenced by GenBoxRay().

int genie::geometry::ROOTGeomAnalyzer::firay

protected

Definition at line 169 of file ROOTGeomAnalyzer.h.

Referenced by GenBoxRay().

bool genie::geometry::ROOTGeomAnalyzer::fKeepSegPath

protected

need to fill path segment "path"

Definition at line 162 of file ROOTGeomAnalyzer.h.

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

double genie::geometry::ROOTGeomAnalyzer::fLengthScale

protected

conversion factor: input geometry length units -> meters

Definition at line 150 of file ROOTGeomAnalyzer.h.

Referenced by LengthUnits(), and SetLengthUnits().

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

protected

matrix connecting master coordinates to top volume coordinates

Definition at line 159 of file ROOTGeomAnalyzer.h.

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

bool genie::geometry::ROOTGeomAnalyzer::fMasterToTopIsIdentity

protected

is fMasterToTop matrix the identity matrix?

Definition at line 160 of file ROOTGeomAnalyzer.h.

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

int genie::geometry::ROOTGeomAnalyzer::fMaterial

protected

input selected material for vertex generation

Definition at line 142 of file ROOTGeomAnalyzer.h.

double genie::geometry::ROOTGeomAnalyzer::fMaxPlSafetyFactor

protected

factor that can multiply the computed max path lengths

Definition at line 152 of file ROOTGeomAnalyzer.h.

Referenced by MaxPlSafetyFactor(), and SetMaxPlSafetyFactor().

double genie::geometry::ROOTGeomAnalyzer::fMixtWghtSum

protected

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

Definition at line 153 of file ROOTGeomAnalyzer.h.

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

double genie::geometry::ROOTGeomAnalyzer::fmxddist

protected

Definition at line 174 of file ROOTGeomAnalyzer.h.

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

double genie::geometry::ROOTGeomAnalyzer::fmxdstep

protected

max errors in pathsegmentlist

Definition at line 174 of file ROOTGeomAnalyzer.h.

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

bool genie::geometry::ROOTGeomAnalyzer::fnewpnt

protected

Definition at line 170 of file ROOTGeomAnalyzer.h.

Referenced by GenBoxRay().

int genie::geometry::ROOTGeomAnalyzer::fNParticles

protected

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

Definition at line 147 of file ROOTGeomAnalyzer.h.

Referenced by ScannerNParticles(), and SetScannerNParticles().

int genie::geometry::ROOTGeomAnalyzer::fNPoints

protected

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

Definition at line 145 of file ROOTGeomAnalyzer.h.

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

int genie::geometry::ROOTGeomAnalyzer::fNRays

protected

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

Definition at line 146 of file ROOTGeomAnalyzer.h.

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

double genie::geometry::ROOTGeomAnalyzer::fox

protected

Definition at line 171 of file ROOTGeomAnalyzer.h.

Referenced by GenBoxRay().

double genie::geometry::ROOTGeomAnalyzer::foy

protected

Definition at line 171 of file ROOTGeomAnalyzer.h.

Referenced by GenBoxRay().

double genie::geometry::ROOTGeomAnalyzer::foz

protected

top vol size/origin (top vol units)

Definition at line 171 of file ROOTGeomAnalyzer.h.

Referenced by GenBoxRay().

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

protected

top volume

Definition at line 158 of file ROOTGeomAnalyzer.h.

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

string genie::geometry::ROOTGeomAnalyzer::fTopVolumeName

protected

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

Definition at line 144 of file ROOTGeomAnalyzer.h.

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

---

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

• Tools/Geometry/ROOTGeomAnalyzer.h
• Tools/Geometry/ROOTGeomAnalyzer.cxx